黄土高原北部水蚀风蚀交错区是黄河下游河床粗泥沙的主要来源区^[1], 也是土壤严重侵蚀与干旱并存的特殊区域, 生态环境非常脆弱。为了防止水土流失, 改善该区域脆弱的生态环境, 近几十年, 黄土高原地区开展了大量的退耕还林还草措施^[2]。基于遥感影像分析, 黄土高原植被覆盖度从1999年的31.6%增加到了2013年的59.6%, 植被恢复在该区域综合治理水土流失和改善生态环境中发挥了重要作用^[3-4]。柠条和苜蓿具有生态适应性强、根系发达以及保水固土效果显著的特点, 又是优质饲料, 是黄土高原地区人工植被建设的主要灌木和草本种类^[5]。因此, 柠条和苜蓿在研究区广泛种植, 是该区保持水土和植被恢复的典型植被。由于柠条和苜蓿的多年生、深根系和强耗水的特性, 其持续生长对土壤水分消耗十分强烈, 导致黄土高原深层土壤干燥化加剧、土壤水分生态环境恶化及人工植被逐渐退化^[6-8]。研究典型人工灌草植被建设下土壤水分的消耗和干燥化过程, 以及土壤干层的形成发育, 对黄土高原植被管理与生态环境重建具有重要的参考价值。

基于不同植被类型下土壤干层与生长年限的相互关系研究, 一些学者发现, 苜蓿、柠条、油松和沙打旺等植被下的土壤干层受植物生长年限影响显著^[9-11]。李军等^[8]研究表明, 3年生苜蓿草地出现中度干燥化, 土壤干化深度可达5.0—7.6 m, 4年生以上的苜蓿草地已达到严重干燥化, 土壤干层厚度超过10.0 m。程积民等^[12]研究发现柠条幼龄期(0—6年)在0—8.0 m土层中未形成干层, 中龄期(6—14年)形成土壤干层, 其厚度为2.4—2.6 m, 老龄期(大于14年)土壤干层厚度达到7.0 m, 土壤干燥化的程度和干层厚度均随着柠条生长年限的延长而增加。然而, 以往关于不同年限植被下土壤干层的研究大部分是借助空间代替时间的方法, 时间序列不连续甚至年限差别较大, 且空间代替时间的方法忽略了土壤属性、地形地貌等对土壤水分的影响差异。因此, 由于数据的不连续以及长期土壤水分观测数据的缺乏, 限制了对不同植被类型下土壤干燥化过程的准确认识。

本文在黄土高原北部的神木六道沟小流域, 基于长期连续的土壤水分定位观测数据, 研究柠条灌木和苜蓿草地土壤干层的形成发育过程及其动态变化规律, 探讨不同生长年限柠条和苜蓿对土壤水分的消耗特征, 确定该区柠条和苜蓿的适宜种植年限, 为黄土高原北部脆弱生态区灌草植被管理和土壤干层调控提供科学依据。

1 材料与方法 1.1 研究区概况

研究区位于黄土高原北部水蚀风蚀交错区的神木市以西六道沟小流域(110°21′—110°23′E, 38°46′—38°51′N)。流域面积6.89 km², 海拔1094.0—1273.9 m, 北依长城, 地处毛乌素沙漠的边缘, 属于黄土高原水蚀风蚀交错带的强烈侵蚀中心。该区属于中温带半干旱气候, 年平均蒸发力785.4 mm, 多年均降水量415.4 mm(1971—2014), 年际降水分布不均匀(图 1), 且年内降水分布主要集中在6—9月, 其降雨量占全年降水量的77.4%, 年均气温8.4℃。地貌为片沙覆盖的黄土丘陵, 主要土壤类型有绵砂土、硬黄土、红黄土和风砂土。流域位于森林草原向干旱草原过渡带, 其植被类型属于灌丛草原类型。随着国家退耕还林还草工程的推进和区域经济的发展, 坡耕地逐步转变为林草地, 特别是苜蓿草地、柠条灌木地和果园林地面积显著增加。

1.2 实验设计

本研究在一远离侵蚀沟、坡度为12—14°、坡向为西北方向的完整坡面, 建设人工柠条和苜蓿植被小区。坡面地形、土壤理化性质和土壤质地等较为一致, 试验小区内的土壤类型均为黄绵土, 从而消除了土壤质地在空间上的差异。小区规格为5.0 m×61.0 m, 并于2004年栽植1年生柠条幼苗, 植被密度为12000丛/hm²；苜蓿小区播种种植, 植株行距为50.0 cm×50.0 cm。为实现土壤水分的长期定位观测, 在每个试验小区从坡顶向坡底依次布设11根中子管, 相互间隔距离为5.0 m, 监测土壤深度为4.0 m。

1.3 数据观测与处理

对于柠条和苜蓿小区, 观测时间均为2004年至2014年, 柠条对应的生长年限为2—12年, 苜蓿对应的生长年限为1—11年。利用CNC503DR型中子仪对0—4.0 m土壤剖面体积含水量进行测定, 其中0—1.0 m土层测定间隔为10 cm, 1.0 m以下土层测定间隔为20 cm。测定频率为每半月测定一次。发生较大降水后加密测定土壤含水量。本研究中所用的土壤含水量数据为不同植被下11个监测样点的平均值。

柠条地上生物量用代表样枝方法测定, 在植被小区外采取与小区内相似的标准枝9枝, 先以110℃高温杀青1小时, 然后在80℃恒温烘箱里烘72小时, 最后用精度为0.01 g电子天平称重。在小区内测定柠条的株高和基径利用植被生长(基径和高)与生物量间的Allometry关系计算柠条生物量。采用全收获法或1.0 m×1.0 m样方刈割法计算苜蓿地上生物量, 苜蓿每年在7月底和10月底刈割两次。

1.4 土壤干层

王云强^[13]研究定义土壤干层的上边界土层为1.0 m, 因为0—1.0 m土层的土壤水分通常能够通过当年的降水补给得到及时的补充, 认为该土层不属于土壤干层的范畴。在本研究中, 我们主要研究的是具有较强稳定性的土壤干层的发生发育过程, 因此, 我们将土壤干层的起始深度定为1.0 m。基于黄土高原北部土壤的基本物理属性, 土壤稳定含水量通常被认为是约等效于60%的田间持水量, 并认为土壤含水量低于土壤稳定含水量时, 即形成土壤干层^[14]。本研究将土壤稳定含水量值确定为田间持水量的60%, 作为土壤干层的上限, 并以此作为判断干层是否出现的标准。

为了便于定量描述各类土地利用方式下土壤干燥化程度, 李军等^[8]定义土壤干燥化指数—SDI(soil desiccation index)为某一土层土壤实际有效含水量占该层土壤稳定有效含水量比值的百分数, 公式表达为：

式中, SDI为土壤干燥化指数, SM为土壤含水量, WM为凋萎含水量, SSM为土壤稳定含水量。

土壤干燥化指数SDI值越小, 表明土壤干燥化强度越高, 可供植物吸收利用的土壤有效含水量也越少。依据土壤干燥化指数SDI值的大小, 土壤干燥化强度可划分为6级(表 1)。

2 结果与分析 2.1 柠条和苜蓿植被对土壤水分的消耗特征

柠条和苜蓿地中的剖面土壤含水量的时间变化趋势相似, 随着植被的生长, 大量土壤水分被消耗, 剖面土壤含水量逐渐降低(图 2)。2—8年生柠条对1.0—4.0 m剖面土壤水分消耗强烈, 剖面含水量均值表现出随着植被生长年限的延长呈现迅速下降趋势, 剖面含水量均值由观测初期的21.4%(2年生)减小到9.4%(8年生), 含水量减少速率每年可达2.0%；之后9—12年生柠条地中剖面含水量均值处于相对稳定状态, 变化幅度很小(8.2%—9.0%)。不同生长年限苜蓿对土壤水分的消耗特征与柠条相似, 1—7年生苜蓿地剖面土壤含水量均值先迅速下降, 含水量均值由观测初期的21.3%(1年)减小到8.6%(7年), 含水量减少速率每年可达2.1%；之后苜蓿地剖面含水量均值相对稳定, 在8.5%—10.5%之间浮动(图 2)。

观测初期, 柠条和苜蓿地土壤水分在1.0—4.0 m剖面上分布比较均匀, 1.0—2.0 m, 2.0—3.0 m和3.0—4.0 m土层储水量接近, 随着植被生长3个土层储水量呈现随深度增加而减小的趋势, 且随着生长年限增大3个土层储水量逐年减小(图 3)。2—8年生柠条和1—7年生苜蓿, 土壤储水量减少量和消耗速率由大到小依次为1.0—2.0 m, 2.0—3.0 m和3.0—4.0 m, 之后各土层储水量趋于接近。由图 3可知, 柠条和苜蓿地土壤水分消耗可分为3个阶段, 第一阶段：快速消耗阶段(种植前5年), 2—6年生柠条地1.0—2.0 m, 2.0—3.0 m和3.0—4.0 m土层储水量的消耗速率分别为：24.3 mm/a, 21.3 mm/a, 16.0 mm/a；1—5年生苜蓿地3个土层对应的储水量消耗速率分别为23.9 mm/a, 23.45 mm/a, 20.0 mm/a。第二阶段：缓慢消耗阶段(种植6—8年)柠条地各层储水量的消耗速率分别为：1.7 mm/a, 10.2 mm/a, 12.5 mm/a；苜蓿地中3个土层对应的的储水量消耗速率分别为3.7 mm/a, 3.1 mm/a, 7.2 mm/a。第三阶段：相对稳定阶段(种植9—11年), 柠条和苜蓿地各层土壤水分都处于相对稳定状态, 且3个土层储水量差异很小, 在80.0—100.0 mm之间变化。在第三阶段, 1.0—2.0 m土层水分对年降水量有一定的响应, 尤其在2013年, 降水较多, 1.0—2.0 m土层水分得到补给, 且苜蓿地中补给量明显大于柠条地。2.0 m以下土层储水量变化幅度较小。

2.2 不同生长年限柠条和苜蓿地土壤干燥化特征

2—12年生柠条地土壤剖面平均含水量、储水量和SDI值均表现出逐年减小的趋势(表 2)。由1.0—4.0 m剖面平均SDI可知, 2—7年生柠条地在1.0—4.0 m剖面上无干燥化现象。但3年生的柠条地土壤剖面已经出现水分亏缺现象, 4年生柠条地开始出现干层, 厚度为0.4 m。6年生柠条地土壤干层深度可达2.4 m, 干层的厚度为1.4 m。8年生柠条1.0—4.0 m土壤剖面出现中度干燥化, 1.0—3.0 m土层干化严重, 3.0—4.0 m土层含水量也接近稳定土壤含水量(表 2和图 4)。9年生柠条1.0—4.0 m土壤剖面全部形成干层, 9—12年生柠条土壤干层深度均超过4.0 m。12年生柠条地1.0—4.0 m剖面处于严重干燥化, 土壤含水量分布于6.5%—9.4%范围内。柠条地土壤干化程度随着生长年限的增长而增强。

1—11年生苜蓿地土壤剖面平均含水量和储水量均表现出逐年减小的趋势(表 3)。1—4年生苜蓿剖面土壤处于无干燥化状态, 但2—4年生土壤剖面均出现水分亏缺, 且亏缺量逐年增大；5年生苜蓿生长量迅速增大, 土壤剖面开始出现轻度干燥化现象, 生长季末干层深度达3.6 m, 干层厚度为2.6 m(图 4)；除10年生苜蓿外, 6—11年生苜蓿土壤剖面干燥化程度逐年加重, 由轻度干燥化到严重干燥化状态, 且7年以后土壤干层的深度超过4.0 m, 厚度大于3.0 m。10年生苜蓿地中土壤水分亏缺量和干层的厚度均减小, 土壤干燥化程度减弱, 主要是由于2013年降水较多(图 1), 降水补给1.0—2.0 m土层。

2.3 苜蓿与柠条生物量随生长年限的变化

柠条和苜蓿地上最大干生物量随生长年限增大, 其变化趋势不同(图 5)。总体上2—12年生柠条地上最大干生物量随柠条生长年限的增长呈逐渐增大趋势。6年和9年生柠条生物量均小于前一年生物量, 主要是由于这两个年份(2008年和2011)均为枯水年, 降水较少影响了植被生物量。已有的研究表明柠条生长被分为三个阶段：幼龄期(0—6年), 中龄期(7—13年)和老龄期(≥14年), 且地上生物量在13年之前随植物生长增加^[12]。因此2—12年生的柠条经历了幼龄期和中龄期, 还未进入老龄期, 所以柠条地上最大生物量随植被年限的增加呈不断增加趋势。

1—11年生苜蓿地上最大干生物量随生长年限的增长呈先增大后减小, 即中间高两头低的趋势(图 5)。苜蓿自播种后生长迅速, 生长至第6年其地上生物量达到高峰, 之后苜蓿产量开始下降, 并逐渐开始退化。在第10年(2013年), 即使年降水量较大, 处于丰水年, 生物量依旧减少。11年生苜蓿地多数苜蓿植株趋于生理衰败, 加之深层水分的大量消耗殆尽, 生物产量极低。研究表明苜蓿地上部生物量与生长年限呈显著的负相关关系, 苜蓿草地生长高峰期为第5—6年, 苜蓿生长至6—7年之后, 产量迅速下降, 其经济生长年限一般为10年^[15-16]。研究期间的苜蓿生物量较柠条低, 但两者的耗水量差异不大, 因此, 柠条的水分利用效率高于苜蓿。白永红等^[17]研究表明, 黄土区柠条的水分利用效率显著高于苜蓿(P＜0.05), 是苜蓿的1.5倍, 与我们的研究结果一致。

3 讨论

黄土高原近50年来的人工林草建设, 人为引进很多外来高耗水植物物种, 人工植被过度耗水打破降水和原生植被之间的动态平衡, 普遍出现了以土壤干层为代表的新的生态环境问题^{[6, 9, 18-19]}。近年来已有大量研究发现, 人工植被大量消耗土壤储水导致土壤水分减少^[20], 几乎所有的人工乔、灌林和紫花苜蓿等多年生豆科牧草都能使土壤形成干燥化土层, 土壤干层现象分布相当广泛^[21-23]。王力等^[14]研究表明, 黄土高原地区人工植被用地与自然植被恢复土地相比, 其土壤水分总是较低且使土壤干化。黄土区, 柠条和苜蓿的耗水量是降雨量的1.3倍和1.2倍, 降雨无法满足柠条和苜蓿的需求, 它们只能过度消耗土壤水分来满足生长, 最终造成土壤干化^[17]。在该研究区, 也有研究表明, 三种典型植被恢复类型土壤水分由大到小依次为撂荒地>柠条地>苜蓿地。观测初期, 三种植被恢复用地含水量剖面分布相似, 但随观测年限增加, 柠条和苜蓿地中土壤水分消耗深度加深, 耗水程度增强, 尤其是在1.0 m以下土层更为明显, 并在种植2年后, 1.0 m以下土层开始出现干化^[24]。在黄土区, 多年种植柠条和苜蓿会导致土壤产生干层, 且干层深度和干燥化程度随着柠条和苜蓿生长年限的延长而增大^[25-26]。

本研究表明, 柠条和苜蓿的大量耗水使得其均在种植5年(6年生柠条和5年生苜蓿)后剖面土壤水分显著减少并形成土壤干层；且剖面土壤水分变化速率在种植5年时存在一个“拐点”, 在这之前土壤含水量和储水量大量减少, 5年之后土壤剖面水分减少速率降低, 深层土壤水分变化较小。种植8年后, 土壤中的储水被消耗殆尽, 土壤水分处于相对稳定状态, 这表明低于该年限所对应的土壤含水量时, 柠条和苜蓿不能吸收利用, 只能依靠当年的降水生存生长。

2—6年生柠条, 随着生长年限的增长, 剖面土壤含水量迅速减少；在6—12年生长期间, 柠条生物量不断增大, 柠条地剖面土壤含水量仍然随着生长年限的增长而呈减少趋势, 但深层土壤剖面水分变化较小, 表明在柠条生长6年之后, 剖面土壤水分随着生长年限降低的速率减小, 深层土壤含水量在不同年限间变化相对较小。这与李小芳等^[27]研究结果一致, 其研究表明随着柠条生长年限延长和柠条根系的下伸, 土壤干层逐年加深加厚, 5年生柠条干层深度达3.0 m, 7年生超过4.0 m, 10年生干层深度接近7.0 m, 11年以后, 1.0—9.0 m土层土壤湿度基本稳定, 接近土壤萎蔫湿度。程积民等^[28]将柠条的生长发育期分为幼龄期(0—6年)、中龄期(6—14年)和老龄期(大于14年)。柠条在幼龄期生长主要依靠天然降水和土壤水分补给；在中龄期, 土壤储水量较少, 植被主要依靠天然降水补给维持植被生长。因此, 本研究建议柠条种植的生长年限在5—6年左右, 最大不要超过6年, 其对应的地上最大生物量为5050 kg/hm²。

本研究中苜蓿在1—6年, 生物量随着生长年限的增长而增大, 为苜蓿生长的旺盛期, 对土壤水分需求增加, 导致土壤水分迅速减少；种植6年后深层土壤水分已产生严重亏缺, 水分供应受到限制, 从而使得土壤剖面水分减小速率变小。因此, 在苜蓿生长过程中, 1—6年生长前期土壤干燥化进程明显快于7—12年生长后期, 前6年耗水深度和干层厚度逐年增加, 7—12年生苜蓿生长衰退, 耗水深度和干层厚度有所减小^[29-30]。苜蓿生长至第5年, 1—3.6 m土层出现干层, 在第10年(2013年), 即使处于丰水年条件下, 仅1.0—2.0 m土层水分得到补给, 其干层得到短时间的恢复, 2.0 m以下干层难以恢复^[31]。表明在研究区苜蓿种植5年后就会产生干层, 且难以恢复。因此, 建议苜蓿在研究区种植年限不要超过5年, 其对应的地上最大生物量为1980 kg/hm²。或者在苜蓿种植的第4年采取一定的水分调控措施, 如条带状减少苜蓿种植, 苜蓿与农作物或自然恢复植被间作, 维持土壤水分平衡, 以防止深层土壤干化层的形成。

4 结论

本文研究了黄土高原北部水蚀风蚀交错区六道沟小流域内两种典型人工灌草植被——柠条和苜蓿, 在不同生长年限下土壤水分时空动态变化特征及土壤干层的形成过程, 取得如下结论：

(1) 种植多年生的柠条和苜蓿会导致深层土壤干燥化, 且土壤干化的程度随生长年限增大而加剧。柠条和苜蓿的耗水速率在种植5年时存在一个“拐点”, 在这之前土壤含水量和储水量大量减少, 5年之后剖面土壤水分下降速率减小。种植8年后, 土壤水分处于相对稳定的低水平状态, 为无效水, 不能被柠条和苜蓿根系吸收利用, 植被只能依靠当年的降水维持生存。

(2) 土壤干层的形成速率和厚度与植被类型有关。苜蓿从第5年开始, 1.0 m以下出现干层, 干层厚度达到2.6 m；柠条从第6年开始形成干层, 厚度达1.4 m。因此, 为减少深层土壤干层的发生发育, 建议苜蓿在该地区的种植年限不超过5年, 其对应的地上最大生物量为1980 kg/hm²；柠条的生长年限不超过6年, 其对应的地上最大生物量为5050 kg/hm²。