塔里木河下游河道输水对沿岸土壤含水率的影响

冯 斐

（塔里木河流域管理局，新疆库尔勒841000）



塔里木河下游河道输水对沿岸土壤含水率的影响

冯 斐

（塔里木河流域管理局，新疆库尔勒841000）

摘要：利用塔里木河下游土壤含水量5年期的监测资料，根据变异系数和方差分析方法研究了输水后塔里木河下游土壤含水率时空变化特性。结论显示出塔里木河下游输水对河道沿岸的土壤含水率有明显影响，输水改善了土壤中的湿度进而影响了植被的生长。

关键词：河道输水；土壤体积含水率；变异系数；塔里木河下游

土壤水是干旱、半干旱地区天然植被生存、生长发育的基本条件，决定植被分布和动态变化，是生态系统结构与功能的关键因素［1～3］。土壤湿度变化是降水、冠层截留、植物蒸腾、土壤蒸发、地表径流、地下渗漏等多种因素综合作用的结果，它对土壤-植被-大气传输体（SVAT）中的物质迁移过程具有重要影响［4～6］。由于中子仪测定法具有快速、省力，不破坏土壤结构等的优点，并且其系统误差和性能及其测定土壤水分的精度已经得到普遍认可［7］，因此利用中子仪研究土壤剖面含水量的变化，掌握干旱区退化生态系统的土壤水的时空变化规律对塔里木河下游受损生态系统的恢复有极为重要的意义。

塔里木河下游自2000年实施生态输水工程以来，断流30多年的下游河段恢复了水文过程的完整性。塔里木河下游两岸分布的自然植被主要是非地带性的隐域植被，它们的生存不依赖于大气降水，而主要依赖于地下水供给其蒸腾和蒸发。但地下水是通过毛细管作用上升补给土壤水分，从而被植被吸收利用。因此，在一定意义上说，地下水是通过改变土壤体积含水率来影响植被的生长。本文利用塔里木河下游不同监测剖面连续近5年的土壤水动态变化数据，试图掌握塔里木河下游不同退化程度区土壤体积含水率的时空分布特征，深化对塔里木河下游生态输水后生态水文过程理解，为塔里木河下游退化天然植被恢复提供科学依据。

1　材料和方法

1.1研究区概况

研究区位于塔里木河下游大西海子至台特玛湖间，地理坐标介于东经87.6O～88.5O，北纬N39.5O～40.6O间（图1），属暖温带荒漠干旱气候，这里降水稀少，多年平均降水量仅为20～50 mm；而年平均蒸发量却高达2500～3000 mm。该地区的地带性植被是温性灌木和半灌木。在河漫滩及两岸的低阶地发育着大面积非地带性的草甸植被，形成由胡杨（PoPuluseuPhratica）、灌木和草本植物组成的面积广阔的乔、灌、草带［8］，塔里木河下游两岸主要土壤类型有绿洲潮土、风沙土、草甸土、盐土等。近50年来，由于人类不合理的水土资源开发，造成塔里木河下游出现了极其严重的生态退化［9～10］。

1.2研究方法

（1）样点设置和数据获取。为全面、准确地掌握塔里木河下游输水后土壤水的变化过程，在离塔里木河下游离输水水源-大西海子水库由近及远不同距离处，选取英苏和依干布及麻两个典型断面，分别垂直于输水河道，按一定间距布设2～3个长期固定土壤水监测剖面。土壤含水量测定采用中子水分仪（CNC503DR型，北京核子仪器有限公司，该仪器测定的是土壤非饱和带的容积水分含量，其主要技术参数为：测量范围是0%～70%；当含水量0%～40%时，观测精度为2%，当含水量40%～70%时，观测精度为3%）。土壤水含水量测定时间是每年每月1次，在输水期间和输水后每月监测2次。文中采用的数据是2002年8月～2006年12月。

（2）数据处理。统计学表明，土壤体积含水率的平均值代表研究尺度内土壤体积含水率的集中趋势，标准差和变异系数（Cv）表明土壤体积含水率对中心位置的偏离（变异性）程度，它们分别为土壤体积含水率数据分布位置、分散情况、分布形状的量度［11～12］。Hlllel（1980）认为当某一土壤特性参数的变异系数Cv≤10%时为弱变异性；当10%＜Cv＜100%时为中等变异性；当Cv≥100%时为强变异性［13］。

图1　塔里木河下游简图

式中：S表示一组数据的标准差；表示该组数据的平均数。

利用单因素方差分析（one-wayANOVA）和多重比较（LSD）的方法计算各断面离河道不同距离间不同土层的土壤体积含水率差异，并对其变化差异进行单一样本T检验。

2　结果与分析

2.1土壤含水率的空间变异性特征

8次输水后（2006年11月），塔里木河下游典型断面的土壤体积含水率在垂直方向上变化的统计特征如表1所示。在地表以下0～280cm土层范围内，整体上随着土壤深度增加，土壤体积含水率呈现先减少后增加趋势。而各层土壤体积含水率变异系数存在明显差异。输水后在地表60cm以下，土壤体积含水率的空间差异性极显著（p＜0.001）。

表1　塔里木河下游输水后主要断面土壤体积含水率（%）统计特征值

在水平方向上，各断面离河道不同距离土壤体积含水率差异性明显（图2左）。在0～280cm土层范围内，离水源较近的英苏断面的土壤平均含水率要明显高于离水源较远的依干布及麻断面的土壤平均含水率，在同一断面离河道较近的土壤平均含水率要高于离河道较远的土壤平均含水率。各断面离土壤体积含水率的空间变异系数规律性也明显（图2右）。

整体上英苏断面的土壤体积含水率空间变异程度低于依干布及麻断面，但土壤体积含水率变化规律不同的是，在离河道不同距离的土壤体积含水率的空间差异性表现出相反的变化规律。英苏断面的离河道较近C1号监测孔，属于中等变异性；在离河道较远的C2号监测孔的变异程度加大。但在依干布及麻断面，离河道较近的H1监测孔土壤体积含水率变异性很高，属于强变异性；离河道较远H2监测孔的土壤体积含水率空间变异性则较低，。这说明，离水源越近，近河道的土壤体积含水率变异性低于远河道的土壤体积含水率变异性，离水源越远，近河道的土壤体积含水率空间变异程度要高于远河道的变异程度。

图2　塔里木河下游输水后英苏和依干布及麻断面监测孔的土壤体积含水率和变异系数变化

通过对各断面不同距离的含水量进行单因素的方差分析（表2），不同土层含水率之间除20cm土层的差异不显著外，其他各土层之间的差异呈极显著（p＜0.001）。对于同一断面，英苏断面除了20cm土层外，其它各土层的C1和C2的土壤含水率差异性都很显著（a =0.05）；而依干布及麻断面的H1和H2的土壤含水率仅在40～60cm和240～280cm土层范围内，两者差异性显著（a =0.05），在80～220cm土层范围内，两者的差异性不显著（a =0.5）。对于不同断面，英苏和依干布及麻断面，除了C2和H2在20～40cm和280cm土层外，四个监测孔间的土壤含水率差异性显著（a =0.05）。

表2　不同监测样点之间的方差分析与多重比较

从不同土壤体积含水率和变异系数随土壤深度的变化曲线可以看出：各剖面土壤体积含水率变化在地表以下0～3m深度范围内上有3个明显不同特征的变化层：水分散失层、水分增加层和水分变化层（表3）。以英苏断面为例，C1号监测孔，0～60cm土层的土壤体积含水率较低，变化幅度较小，变异系数较低；60～100cm土层的含水率显著增加，增幅最大，变异系数较高；100～280cm土层的含水率仍在上升，但增幅有所减缓，变异程度明显下降。

表3　塔里木河下游典型断面各土层土壤体积含水率变化（%）

2.2土壤含水率的时序性特征

图3　英苏和依干布及麻断面监测孔土壤体积含水率时间变化

受间歇性输水的影响，塔里木河下游各断面土壤体积含水率随时间变化规律十分明显（图3）。以英苏断面为例，离河道较近C1监测孔土壤体积含水率随着间歇性输水变化十分明显。在0～140cm土层内，土壤体积含水率在输水期间增加幅度很快，如在第4次输水后（2002年12月）土壤平均含水率为0.03 cm3/cm3，第五次（I）输水后（2003年7月）增加到0.09 cm3/cm3，第7次（I）输水后（2005年8月）增加到0.12 cm3/cm3。在停水期间的土壤体积含水率略有下降，但整体上随着输水次数和时间的延续土壤体积含水率呈现增加趋势。在140cm土层以下，随着深度增加，土壤体积含水率随输水时间变化的升降幅度减弱。在160～280cm土层范围内，土壤体积含水率随着时间延续整体上开始呈现下降趋势，但在整个输水过程中其升降幅度明显，尤其是在240cm土层以下，土壤体积含水率的升降变化与输水时间和次数关系密切。

离水源距离较远的依干布及麻断面土壤体积含水率随时间变化明显晚于C断面。以H1监测孔为例，0～120cm土层在2002年9月～2003年12普遍较低。在2004年2月～2006年12月期间，土壤体积含水率增加幅度很大；120～200cm土层土壤体积含水率增加的时间比0～120cm土层要早。120～160cm土层在2003年12月5日左右土壤体积含水率达到高峰，而180cm土层和200cm土层分别2003年11月和10月中旬达到最高值。此后至2004年2月27日，即在第六次输水之前，其各土层土壤体积含水率降低到最低；200～280cm土层的土壤体积含水率在第5次（I）和第5次（II）输水中出现两个高峰期。此后各层土壤体积含水率随着输水时间的延续幅度也较低。

3　结论

（1）塔里木河下游8年输水后土壤体积含水率在空间上有很强的变异性。垂直方向上，土壤体积含水率随着土壤深度增加整体呈先减少后增加趋势，且各层土壤体积含水率变异性存在明显差异。

（2）根据各断面土壤体积含水率和变异系数随土壤深度的变化特征，塔里木河下游主要断面土壤体积含水率变化在地表以下0～3m深度范围内上有3个明显不同特征的变化层：水分散失层、水分增加层和水分变化层。各断面在同一土层深度下土壤体积含水率差异显著。

（3）塔里木河下游输水后土壤体积含水率变化在时间上具有不稳定性。随着土层深度增加和间歇性输水的延续，离河道较近的土壤体积含水率对输水过程的响应要明显高于离河道较远的土壤体积含水率对输水过程的响应；在离水源距离较近的土壤体积含水率变化明显早于离水源较远的土壤体积含水率变化。

参考文献

［1］赵传燕，冯兆东，南忠仁.黄土高原西部土壤水分时空变化模拟研究-以安家坡流域为例.冰川冻土，2007，29（05）：785-794.

［2］Baird AJ，WilbyRL，Ecohydrology：Plantsand waterin terrestrial and aquaticenvironments［M］.London：Routledge，1998.

［3］李王成，冯绍元，康绍忠，等.石羊河中游荒漠绿洲区土壤水分的分布特征.水土保持学报，2007，21（03）：138-157.

［4］王云强，张兴昌，从伟，等.黄土区不同土地利用方式坡面土壤含水率的空间变异性研究，农业工程学报，2006，22（12）：65-71.

［5］邱扬，傅伯杰，王军，等.土壤水分时空变异及其与环境因子的关系.生态学杂志，2007，26（01）：100-107.

［6］肖洪浪，李锦秀，赵良菊，等.土壤水异质性研究进展与热点.地球科学进展，2007，22（09）：955-959.

［7］毛飞，任三学，刘庚山，等.中子仪测量农田土壤水分精度的比较研究.中国生态农业学报，2005，13（04）：103-106.

［8］Xu HL，YeM，SongYD，etal.Thenaturalvegetation resPonsesto thegroundwaterchangeresultingfrom ecologicalwaterconveyances totheLowerTarimRiver.Environ MonitAssess，2007，131：37-48.

［9］宋郁东，樊自立，雷志栋.中国塔里木河流域水资源与生态问题研究［M］.新疆人民出版社，1999：398.

［10］王让会，樊自立.利用遥感和GIS研究塔里木河下游阿拉干地区土地沙漠化［J］.遥感学报，1998（02）：137-142.

［11］刘玉玫.统计学基础［M］.北京：中国统计出版社，2002：101.

［12］吕贻忠，胡克林，李保国.毛乌素沙地不同沙丘土壤水分的时空变异.土壤学报，2006，43（01）：152-154.

［13］胡顺军，康绍忠，宋郁东，等.渭干河灌区土壤水盐空间变异性研究.水土保持学报，2004，18（02）：10-20.

中图分类号：S152.7

文献标识码：A

文章编号：1008-1305（2016）01-0079-05

DOI：10.3969 /j.issn.1008-1305.2016.01.028

收稿日期：2015-08-24

基金项目：水利部公益项目和优秀青年科技人才培养项目（2013721032）项目资助。

作者简介：冯斐（1981年—），男，工程师。