黄艳丽，李占斌 ，苏 辉，柏兰峰，孙宝洋，刘晨光

(1.西北农林科技大学水土保持研究所，黄土高原土壤侵蚀与旱地农业国家重点实验室，杨凌 712100；2.河南理工大学测绘与国土信息工程学院，焦作 454000；3.西安理工大学西北水资源与环境生态教育部重点实验室，西安 710048；4.中国科学院水利部水土保持研究所，杨凌 712100；5.新乡学院生命科学技术学院，新乡 453003）

0 引 言

土壤水分通常是指非饱和带中的水分[1]，是土壤质量的重要标志，更是黄土高原植物生命活动的直接水源[2]。土壤水的存贮、补给、消耗、更新和平衡，对农牧业生产与生态环境都有重要意义，对干旱和半干旱地区的水资源平衡影响巨大[3]。作为退化生态系统结构与功能的关键因子[4]，土壤水分对黄土高原生态恢复与重建影响巨大。

1999年，随着“山川秀美”工程的实施，退耕还林（草）在改善黄土高原生态环境、减少水土流失[5]的同时，引起土壤水分严重亏缺[6-7]并形成干层[8-9]。当前，黄土高原在立地[10]、坡面[11]、小流域[12]、区域[13-15]尺度上均开展了土壤水分研究。立地层面土壤水分受土地利用、地形地貌影响[16-17]，但土地利用对深层土壤含水量（Soil water content，SWC）影响更显著[18-19]。坡面土壤水分分布受地形[20]与土地利用[21]的综合影响，从坡顶到坡下先减小后增大[22]，阴坡＞山脊＞阳坡，40～200 cm土层垂向分布先减小后增大；变异程度坡上＞坡中＞坡下，阳坡＞阴坡，垂向变化与 SWC明显正相关[23]。小流域土壤水分受坡向、坡位、土地利用类型影响显著[24]，但各因素影响具有时空变异性[11,25]。贾小旭等[26]认为黄土高原土壤蓄水量从南向北递减，随深度增加SWC空间变异性增强、时间变异性减小，干燥度、黏粒、归一化植被指数和坡度是土壤水分的重要影响因素。王亚飞等[27]指出植被恢复方式是土壤水分变化的主要影响因素，生长年限的影响相对较弱，多种植被恢复方式的土壤储水量大小顺序为弃耕地＞荒草地＞杏树地＞苜蓿地。土壤水分研究范畴与深度的拓展，增强了对黄土高原土壤水分分布、变化、影响因素及变化环境下动态响应的科学认识，但已有研究主要集中于不同植被[28]、立地、小流域与区域土壤水分状态[29]、储量[30]的监测、模拟及其分布格局、变化过程的探索，缺乏对不同生态治理方式影响下流域尺度特别是小流域尺度土壤水分的系统对比研究。而不同生态治理小流域土壤水分的静态分布与动态过程的对比研究既是黄土高原以小流域为基本单元的生态建设规划、布局的基础[31]，又是科学认识土壤水文变化的关键，而自然恢复与人工林治理分别作为产生土壤水分正负效应的2种典型治理方式尤其缺乏小流域尺度上的对比分析。因此，本研究以甘肃省庆阳市人工刺槐小流域杨家沟及其相邻自然恢复小流域董庄沟为研究对象，在小流域尺度上以固定时间步长对照监测两流域不同区段、坡向、坡位、深度的土壤水分，对比研究两者的空间分异与时间变化，并分析各因素对流域土壤水分空间与时间变化的作用，为小流域生态建设规划、布局与土壤水资源调控、预测提供科学依据。

1 材料与方法

1.1 研究区概况

杨家沟、董庄沟是甘肃省庆阳市南小河沟流域（107°30′～107°37′E，35°41′～35°44′N）的 2 条分支流域，处蒲河下游、董志原西侧，属黄土高塬沟壑区典型（图1），是黄河水利委员会 20世纪设立的水土综合治理试验流域。多年平均气温 9.3℃，最高气温39.6℃，最低气温-22.6℃，最大日温差23.7℃，年温差62.2 ℃，无霜期155 d。多年平均降雨量556.5 mm，最大降水量805.2 mm，最小降水量319.8 mm。两流域面积分别为0.87、1.15 km2, 沟长1500、1600 m，沟道比降10.67%、8.93%，沟道走向由北向南，谷坡一般是40～60°的陡坡和大于60°的悬崖、立壁；沟床一般是黄土或黑垆土；塬面、山坡为黄土，干容重1.4g/cm3，为粉砂壤土。作为两条相邻支沟，杨家沟与董庄沟土地利用、地形地貌、土壤条件极为相似，土壤侵蚀也同样严重。20世纪50年代杨家沟开始通过人工载植刺槐进行生态治理，而董庄沟一直作为杨家沟的对比流域实施封育、自然恢复，历经60a治理，两流域分别发育形成了稳定的人工林草（刺槐、山杏、沙棘、马牙草、冰草、艾蒿）与天然荒草地植被群落（马牙草、冰草、艾蒿），常被用作生态效用研究的对比小流域。

图1 研究流域位置图Figure.1 Location map of study area

1.2 试验方法

1.2.1 样点选择原则与方法

在杨家沟的上游（沿流域走向距沟头200～300 m）、下游（距沟口200～300 m）分别选择一组半阴坡（东坡，45～135º）与半阳坡（西坡，225～315º）作为典型坡面（4条），坡面各自然条件保持一致：海拔1150～1250 m、坡度 35～40º、坡长 25～30 m间，植被包括 50龄刺槐（60%～75%）、灌丛（30%～50%）与草被（80%～95%）。为保证两流域SWC的可比性，根据杨家沟各典型坡面的区位、坡向、坡度、坡长、海拔等特点在董庄沟一一选择、确定对照坡面（4条），坡面天然草被盖度在70%～80%间。然后在两流域各典型坡面的中央条带上顺坡向下每隔8～10 m依次确定坡肩、坡腰、坡脚的采样点。杨家沟各样点都与刺槐与灌木植株保持3 m以上距离以减小冠层截留对土壤水分的影响。

1.2.2 采样方法

为对比研究两流域土壤水分，选择林木蒸腾耗水需求最旺盛的 8月份汛期作为采样时期。在杨家沟与董庄沟的对照坡面的对照坡位、每 4 d同时采土样一次，从2016年8月2日至8月30日共采样8次。用4 cm土钻在0～120 cm土深内分9层取样（0～60 cm土深每10 cm、60～120 cm土深每20 cm作为一层），每层取3个重复样品。每次取样后尽量恢复样点原貌并做标记，下次取样应错开之前历次取样点。

1.2.3 样品测试方法

采样后立即测试土壤样品湿质量，带回实验室、放入烘箱中持续105℃高温烘烤12h至恒质量，测干质量并计算土壤样品的质量含水量，取 3个重复样品平均值作为测试值。

1.3研究方法

1.3.1 数据处理与分析

采用Excel2010进行数据整理与统计。通过SPSS18.0的单因素方差分析与独立样本 T检验进行土壤含水量的差异性检验，选择最小显著性检验方法（LSD法）进行多重比较，借助球形检验与多变量方差模型分析小流域各空间因素对 SWC空间与时间变化的影响。使用Origin9.0制图。

1.3.2 相关指标与计算公式

1）剖面含水率

式中，θin是样点i第n次采样的剖面含水量，%；m是土壤采样划分的层次；θinj是样点i第n次采样第j层样品的SWC，hj是土壤采样第 j层的土层深度，cm；H是土壤采样总深度，cm。

2）变异系数

式中，CV是SWC的变异系数，%；SD是SWC的标准偏差；MN是SWC平均值，%。

3）土壤水分流域变差系数

为了定量分析杨家沟与董庄沟对照样点间的土壤水分差，构造了土壤水分流域变差系数指标。土壤水分流域变差系数是两流域各区位、坡向、坡位、土层的SWC之差与其平均值的相对百分数，既能定量反映两对照土壤样品的土壤水分差异，通过其平均值标准化又能适当规避土壤湿度高低对土壤水分变化的影响。为保持指标一致性，始终以董庄沟SWC作为参照。计算公式如下：

其中，DRSM是杨家沟（Y）相对于董庄沟（D）的土壤水分流域变差系数，θY、θD分别是杨家沟与董庄沟某一土地类型的SWC。

2 结果与分析

2.1 杨家沟土壤水分特征

杨家沟2016年8月各样次剖面SWC依次为13.29%、11.93%、10.71%、10.49%、10.52%、10.26%、13.41%与11.46%（图2a）。随气候变化，整个监测期可以分为4个阶段：第一阶段（8月2—10日），前期降雨影响下相对较高SWC的快速下降；第二阶段（8月10—22日），持续干旱影响下相对较低SWC的震荡稳定；第三阶段（8月22—26日），8月23、24日累计超过60 mm降雨影响下土壤水分快速补充；第四阶段（8月26—30日），相对湿润状态下SWC快速下降。

图2 土壤水分及其空间变异系数的时间变化Fig.2 Mean value and coefficient of spatial variation(CVS) of SWC in Y and D

杨家沟SWC空间变异性中等；空间变异程度上游＞下游，东坡＜西坡，各坡位变异程度基本相当，坡腰稍高于坡肩与坡脚；剖面内，变异系数在25.15%～30.62%之间上下浮动，50～60 cm土层最大。时间变异系数上游＜下游，东坡＜西坡，坡肩＜坡脚＜坡腰，除东坡与坡肩属弱变异外均为中等变异，坡腰与西坡SWC时间稳定性最差；各土层时间变异系数随土深增加不断下降（表1）。

杨家沟剖面SWC除8月22、26日上游稍高于下游外，其余下游＞上游，上、下游间差异不显著。变化过程中，第一、四阶段SWC下游远高于上游；第二阶段，两者相对平衡；第三阶段，上游＞下游。总体上，入渗雨水在土壤中的再分配会扩大上游与下游SWC差距，而干旱与降雨补充都会减小差距。除8月26日外剖面SWC东坡＞西坡，降雨补给后（8月2、26日）两坡向SWC最接近，持续干旱、蒸发影响下（8月10、22日）差距最大。不同坡位SWC相比，坡脚＞坡腰＞坡肩，分别是流域均值的1.25、0.91与0.84倍，除8月30日外其余各样次坡下显著高于其他坡位，坡面土壤水分分布不均。剖面内，SWC由土表向下至土深100 cm不断下降，100～120 cm土层转而上升（表1）。

2.2 董庄沟沟土壤水分特征

董庄沟各样次剖面SWC依次为15.18%、14.99%、13.59%、13.42%、13.32%、12.67%、15.86%与15.90%（图2a）；4个阶段总体上经历了下降（8月2-22日）与上升（8月22-30日）2个过程，但2个过程中速度变化一致，都是先疾（降、升）后缓（降、升）。

SWC上游＞下游，除降雨产流后的8月26、30日各样次上、下游间差异显著。不同坡向SWC相比，东坡＞西坡，除8月10、26日外各样次两坡向差异显著（P＜0.05）。各坡位SWC相比坡脚＞坡腰＞坡肩，分别是流域均值的1.18、1.01、0.81倍，但差异不显著。（P＜0.05）。SWC沿剖面向下先短暂下降（0～20 cm）后不断上升（20～120 cm）（表 2）。

SWC空间变异系数除东坡＜20%外，其余各坡面均在20%～40%间变化，上游＜下游，东坡＜西坡，坡肩＜坡脚＜坡腰，沿剖面向下先增大（0～60 cm）再减小（60～120 cm）。下游、坡肩、40 cm以下土层SWC时间变异系数大于 10%，属中等变异，而流域上游、东（西）坡、坡腰（脚）与40 cm以上土层均属弱变异，土壤水分空间格局较稳定。

表1 杨家沟土壤含水量的平均值与变异系数Table 1 Mean value and coefficient of variation(CV) of SWC in Yangjiagou %

2.3 土壤水分的对比分析

杨家沟剖面 SWC小于董庄沟，相当于后者的72.04%～87.56%，除 8月 2日外两流域均差异显著（P＜0.05,图 2a）。杨家沟 8月 2日 SWC相当于董庄沟的87.56%，与之最接近； 8月3、6日依次0.8、5.1 mm的降雨未能改变两流域SWC下降的趋势，8月6日SWC分别下降1.36%与0.19%，前者减少量是后者的7.16倍，下降到后者的79.58%，两者差距进一步扩大；第二阶段杨家沟SWC在董庄沟的78.15～80.94%间浮动。第三阶段降雨使杨家沟、董庄沟SWC较前次分别增加了3.15%与3.18%，前者迅速提高到后者的84.58%；但在后续降雨再分配与29日0.8 mm降雨的影响下，前者8月30日SWC较前次下降了1.96%，而后者不降反有微量上升，这使前者进一步下降到后者的72.04%，差距达到最大。

除8月10日外，SWC空间变异系数杨家沟＜董庄沟；两者变化趋势基本一致，但干旱影响下董庄沟空间变异性增大趋势明显（图2b）。

两流域SWC时间变异系数分别为10.99%与8.76%，杨家沟土壤水分空间分布格局的时间稳定性稍弱于董庄沟。

在地形等微域分异影响下，流域不同土地类型的SWC差异显著[32]；在植被作用与微域因素的交互影响下人工林与自然恢复小流域的SWC差距在不同区段、坡向、坡位与土层上应该有不同表现。

表2 董庄沟土壤含水量的平均值与变异系数Table 2 Mean value and coefficient of variation(CV) of SWC in D %

2.3.1 各区段土壤水分对比

杨家沟与董庄沟不同区段SWC相比，上游差异显著，下游差异不显著（P＞0.05）。上游土壤水分流域变差系数在23.02%～44.75%间变化，8月26日最小，8月10日最大，均值33.98%；下游在-7.79%～26.60%之间，均值仅有9.80%，是上游的28.84%，8月2日杨家沟下游SWC高于董庄沟，8月30日两者下游SWC相差最多（表3）。由不同区段SWC流域变差系数的时间变化可知，杨家沟下游SWC部分时间会接近甚至超过董庄沟，相较上游，下游土壤水分亏缺不明显。SWC的空间变异性相比，上游与下游均是杨家沟＜董庄沟，但下游差异明显；时间变异性均是上游较弱，但各区段前者高于后者。

董庄沟剖面SWC上游＞下游，主要表现在20 cm以下土层，且上游坡面坡脚相比坡腰与坡肩SWC与下游坡面同坡位的差距更大。这是由于天然草被消耗土壤水分相对较少，使流域深层SWC处于较高水平，汛期暴雨影响下流域较快产流，沟道径流常在上游低洼处形成短期滞水，后期不断回馈坡面特别是坡脚深层土壤水分而保持上游土壤水分的较高水平，而下游相对开阔的沟道环境不仅增加了土壤水分的蒸散发消耗，也加快了沟道水的排泄，造成董庄沟SWC上游＞下游。而杨家沟由于林木的高耗水作用使土壤水分经常处于亏缺状态，暴雨产流几率大大减少，下游的土壤浅层（0～40 cm）含水量稍高于上游，且浅层SWC的下游优势在降雨补给后的干旱初期尤其明显，可能在于下游相对上游良好的植被条件产生的更强的蒸发调节作用更多地减少了浅层土壤水分的蒸发消耗。

2.3.2 各坡向土壤水分对比

同一坡向SWC杨家沟＜董庄沟，除8月10日外两流域东坡SWC差异显著，而西坡各样次差异均不显著（P＜0.05）。东坡土壤水分流域变差系数平均值28.62%，西坡仅有14.83%（表3），这种东、西坡向间土壤水分流域变差系数的高低悬殊的主要原因：1）人工林对不同坡向土壤蒸发的调节能力不同，2）土壤水分差异诱发的阴阳坡生境梯度上人工林群落密度、组成及根系呈现坡向性分布差异。由于林木过度耗水杨家沟SWC低于董庄沟，但林木丰富的冠层、枯落物等明显削弱了土面蒸发，尤其是削弱了辐射时间相对更长的半阳坡（西坡）的土面蒸发，即植被对半阳坡的土壤蒸发调节量相对更大，这种差异性蒸发调节使原由植被耗水差异造成的流域间土壤水分差在东、西坡向间获得了不同程度的削弱，一定程度上导致两流域东坡土壤水分差异显著而西坡不显著（P＜0.05）。同时，阳坡植物较强的抗逆性引起植被密度、组成等一系列自适应调整，降低了土壤水分消耗，这种变化的诱因是土壤水分[33]，也就是说与东坡相比，杨家沟西坡更加匮乏的土壤水分诱发的植被抗逆反应降低了流域西坡的土壤水分消耗，也在一定程度上削弱了与董庄沟西坡的SWC差。两流域SWC空间变异系数相比，东坡前者＞后者，西坡相反，后者不同坡向SWC空间变异性差距较大；时间变异系数杨家沟西坡最高，杨家沟东坡与董庄沟东、西坡相近，均属于弱变异。

2.3.3 各坡位土壤水分对比

各坡位 SWC杨家沟＜董庄沟，坡脚、坡腰、坡肩SWC流域变差系数平均值分别为 16.22%、35.23%与19.09%，坡腰明显大于坡肩与坡脚。由各坡位SWC流域变差系数的时间变化来看，8月2日坡脚、坡腰最小，8月26日坡肩最小，降雨入渗能在一定程度上减弱人工林坡面土壤水分亏缺（表3）。

对比两流域各坡位的空间变异系数，除8月10日坡腰、坡肩杨家沟＞董庄沟外，前者其余各样次各坡位空间变异性均小于后者同一样次的对照坡位。前者空间变异性坡腰＞坡肩＞坡脚，后者却是坡腰＞坡脚＞坡肩，与胡伟等[23]结论不同，原因是文献[23]研究40 m×280 m退耕裸坡地40～200 cm土层的土壤水分分布，研究坡面的特征、利用方式与剖面层次均不相同。作为林（草）坡面水分运动中转站的坡腰水分运动复杂、多变，是其SWC空间变异性各坡位最大的原因。由于坡面雨水的顺坡向下汇聚，坡肩土壤水分获得降雨入渗补给的时间与数量均小于坡腰与坡脚，尤其小于坡脚，因此坡肩SWC大多低于坡腰与坡脚，这种坡面SWC的不均衡分布格局在土壤水分亏缺的杨家沟尤其明显。坡肩SWC显著较低（P＜0.05），这造成杨家沟坡肩林木与草被呈斑块分布、个体差异巨大，植被的空间变异引起土壤蒸散发的空间变异，并最终引起SWC空间变异性增大。而坡脚汇聚坡面雨水由此获得更长时间与更多数量的降雨补给，相对较高的SWC较大程度地保障了人工林的耗水需求，一般是坡面上人工林生长最旺盛、植被分布最均匀的坡位，因此其土壤水分空间分布变异性弱于坡腰与坡肩，变异系数最低。而董庄沟坡肩与其坡腰、坡脚相比较低的SWC并不会局限耗水需求相对较低的天然草被地生长，由此造成的SWC空间变异微乎其微；但坡脚土壤水分常与频发的沟道径流互馈。因此，董庄沟土壤水分空间变异系数坡脚＞坡肩。

两流域不同的植被与土壤水文特征也决定了两者各坡位SWC的时间变化，杨家沟坡腰时间稳定性最差，董庄沟坡脚时间稳定性较差。

表3 杨家沟和董庄沟的土壤水分流域变差系数的时间变化Table 3 Time variation of difference ratio of SWC (DRSM) between Y and D %

2.3.4 各土层土壤水分对比

杨家沟与董庄沟各土层SWC相比，表层（0～10 cm）前者显著高于后者，深层（60～120 cm）后者显著高于前者（P＜0.05），这是沿剖面向下前者 SWC呈下降趋势、后者呈上升趋势的必然结果（图3），与植被措施影响下土壤浅层较深层湿润[15]、人工林过度消耗深层土壤水分的研究结论一致。

根据土壤水分流域变差系数的土层变化（表3），除董庄沟超渗产流的8月26日，20cm以上土层DRSM均是负值，即20 cm土层内SWC杨家沟＞董庄沟；且随干旱愈烈，DRSM呈负值的土层愈深，绝对值愈大，8月22日40 cm以上土层DRSM＜0，0～10 cm与10～20 cm土层甚至达到-48.00%与-43.60%，杨家沟40 cm 以上土层不存在土壤水分亏缺。

8月26日由于流域产流董庄沟土壤浅层大多处于滞水状态，造成0～20 cm土层DRSM＞0，即浅层SWC董庄沟微高于杨家沟。8月30日在土深50cm范围内DRSM＜0，即经过降雨入渗补充后0～50 cm土层的SWC杨家沟＞董庄沟，说明50 cm土深内杨家沟的持水能力强于董庄沟。监测期内50～120 cm土层内DRSM＞0，深层 SWC杨家沟始终小于董庄沟；80～100 cm 土层DRSM平均值最高，是监测范围内杨家沟缺水最严重的土层。

SWC空间变异性相比，除表层（0～10 cm）外杨家沟＜董庄沟，但两流域均是50～60 cm土层最强。同样降雨条件下（8月 23、24日）董庄沟更快产流，8月26日两流域表层（0～10 cm）SWC空间变异性均最小，董庄沟（13.32%）＜杨家沟（17.52%）。

杨家沟与董庄沟各土层SWC时间变异系数的垂向变化趋势一致，随土深增加SWC空间分布格局时间稳定性不断增强，这与贾晓旭等[26,34]土壤水分时空变异性随土层深度下降的结论一致。0～20 cm土层后者＞前者，20～120 cm土层前者＞后者，这是由于浅层土壤水分主要受气候、地形与植被影响，在相同的气候与地形条件下杨家沟的林木冠层与枯落物通过截留降雨、抑制蒸发有效减弱了浅层SWC的时间变异性；而刺槐人工林相较荒草地扎根深、耗水多，提高了剖面深层SWC空间分布的时间变异性。

图3 两个小流域不同土层的土壤含水量随时间变化Fig.3 Variations of SWC of D and Y with time in different soil layers

3 讨论

3.1 人工林小流域土壤水分亏缺的空间异质性

黄土高原林草建设过度消耗土壤水分、造成土壤干化是普遍的共识，但如同小流域土壤水分的空间异质性，人工林小流域的土壤水分亏缺与干化是否也存在空间异质性？杨家沟与董庄沟 SWC对比研究中，下游、西坡、坡脚与坡肩、0～60 cm土层的土壤水分流域变差系数相对较小，而上游、东坡、坡腰、60～120 cm土层相对较大；以土壤水分可持续利用的董庄沟为参照[28]，杨家沟下游、西坡、坡脚与坡肩、0～60 cm土层SWC亏缺相对较小，而上游、东坡、坡腰与60～120 cm土层是造成流域土壤水分亏缺的主要区域或层次。这种流域土壤水分相对亏缺程度的空间异质性是人工林影响下小流域土壤水分空间再分配而形成的土壤水分新格局，可以作为生态建设中植被选择与空间配置的直接依据。可见，不同植被与地形因素组合产生不同的土壤水分条件，以解决水土矛盾为目标的黄土高原小流域生态治理必须适地（水）适树（草）、在遵循土壤水分分布的地带性特征的同时研究并遵循非地带性特征，尤其要研究并遵循林草种类与土地类型交互影响下所形成的流域土壤水分分布特征。

3.2 主要因素对小流域土壤水分时空变化的影响

在分析各空间因素对 SWC时空变化的影响效用前，先对8次SWC监测数据的相关性进行球形检验，检验结果P＜0.05，说明8次数据具有相关性，测量数据不满足球形分布假设，不能进行单因素方差分析，因此对 SWC数据进行多变量方差分析。主体内效用分析结果表明杨家沟 SWC的各样次差异具有统计学差异，不同坡向、坡位SWC的时间变化趋势是一致的，但不同区段、土层SWC的时间变化趋势不同；主体间效用分析结果则表明坡位、土层是影响杨家沟SWC分布的重要空间因素，而坡向[35]、区段对SWC分布的影响不显著（P＜0.05）。董庄沟各样次SWC同样具有统计学差异，但不同坡向、区段、土层SWC的时间变化趋势不同、这3个空间层面的SWC分布格局随时间而变化，仅不同坡位 SWC的时间变化趋势是相同的，因此，流域总体 SWC空间分布的时间变化趋势性不强；区段、坡向、坡位、土层 4个空间因素对流域 SWC分布影响均显著（P＜0.05）。

综上，人工林的水文调节功能弱化了地形、地貌等环境分异引起的SWC的空间变异性；而8月汛期林木旺盛的蒸腾耗水需求与频繁的降雨补给使杨家沟土壤处于暴干暴湿的循环交替中，提高了SWC的时间变异性。

4 结论

通过对人工刺槐小流域杨家沟与自然恢复小流域董庄沟2016年8月不同区段、坡向、坡位、土层土壤水分的变化及其对比分析，结论如下：1）董庄沟SWC上游＞下游，东坡＞西坡，坡脚＞坡腰＞坡肩，沿剖面向下先短暂降低（0～20 cm）又持续增大（20～120 cm），不同区段、坡向间差异显著（P＜0.05）；2）在林木过度耗水及其冠层、枯落物蒸发调节的影响下，杨家沟与董庄沟相比SWC整体下降（P＜0.05）、土壤水分格局发生变化，下游与东坡土壤较湿润，SWC坡脚显著高于坡腰与坡肩（P＜0.05）、剖面垂向分布先持续下降（0～100 cm）后转而上升（100～120 cm）；3）区段、坡向、坡位与土层对董庄沟 SWC空间变化的影响均显著，仅坡位、土层对杨家沟 SWC空间变化影响显著，SWC空间变异性董庄沟（29.19%）＞杨家沟（23.51%）；杨家沟不同坡向、坡位的SWC时间变化趋势是一致的，董庄沟仅不同坡位的SWC时间变化趋势呈一致性，SWC时间变异性杨家沟（10.99%）＞董庄沟（8.76%）；4）杨家沟下游、西坡、坡脚与坡肩、0～60 cm土层SWC亏缺相对较小，而上游、东坡、坡腰与60～120 cm土层是造成流域土壤水分亏缺的主要区域或层次，人工林影响下 SWC的相对丰缺格局可以作为小流域生态建设植被配置与布局的直接依据。

以土壤水分流域变差系数定量评价人工林相对自然恢复小流域的土壤水分亏缺程度，对认识黄土高原小流域生态水文变化及指导生态治理的效果评价、规划布局均具有重要意义。然而，本研究集中对比了杨家沟与董庄沟8月份上（下）游、东（西）坡谷坡坡面120 cm剖面内的SWC，研究的时间与空间范围有限，导致研究结果的局限性。因此，进一步提高土壤水分监测的时空范围，系统对比不同生态治理小流域土壤水分的时空动态变化是今后的研究重点。