生态学报  2017, Vol. 37 Issue (22): 7576-7584

文章信息

周洪华, 李卫红, 李玉朋, 王玉阳, 黄湘.
ZHOU Honghua, LI Weihong, LI Yupeng, WANG Yuyang, HUANG Xiang.
基于树木年轮技术的塔里木河下游河岸胡杨林生态需水量研究
Estimation of ecological water demand of a desert riparian forest using tree rings of Populus euphratica in the Lower Tarim River
生态学报. 2017, 37(22): 7576-7584
Acta Ecologica Sinica. 2017, 37(22): 7576-7584
http://dx.doi.org/10.5846/stxb201609211901

文章历史

收稿日期: 2016-09-21
网络出版日期: 2017-07-12
基于树木年轮技术的塔里木河下游河岸胡杨林生态需水量研究
周洪华 , 李卫红 , 李玉朋 , 王玉阳 , 黄湘     
中国科学院新疆生态与地理研究所, 荒漠与绿洲生态国家重点实验室, 乌鲁木齐 830011
摘要: 生态需水量计算是干旱脆弱生态区恢复重建的一个关键问题。分析了塔里木河下游荒漠河岸林关键种-胡杨树木年轮近90年来的变化特征及对气候水文过程的响应,并基于树木年轮技术提出了维系荒漠河岸林不同恢复状态的生态需水量。研究结果表明,塔里木河下游胡杨树木年轮主要承载的是区域水文历史变化信息,可以作为定量评估生态输水工程的生态恢复效应和定量计算植被生态需水量的新指标;胡杨标准年轮宽度指数与生长年生态输水量呈显著正相关(P < 0.05),近15年生态输水对胡杨年轮宽度指数增加的平均贡献率为42.96%;若要维持塔里木河下游英苏段垂直河道300 m范围内的胡杨达到近90年来的平均生长水平,则需在生长年下泄生态需水量0.84×108 m3,若要达到断流前(1933-1974年)的平均径向生长,则每个生长年内生态需水量应达到0.91×108 m3
关键词: 胡杨     树木年轮     生态输水工程     生态需水量    
Estimation of ecological water demand of a desert riparian forest using tree rings of Populus euphratica in the Lower Tarim River
ZHOU Honghua , LI Weihong , LI Yupeng , WANG Yuyang , HUANG Xiang     
State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
Abstract: Estimation of the ecological water demand of natural plants is a key issue for the restoration and reconstruction of a fragile ecological region in arid areas. Tree rings of Populus euphratica in the Yingsu section of the lower reaches of the Tarim River and their responses to climate change and hydrological conditions were analyzed using the nonparametric Mann-Kendall trend test, Mann-Kendall-Sneyers test, moving t-test, one-way ANOVA, and Pearson correlation. The ecological water demand of a desert riparian forest in the Yingsu section was then estimated based on the tree ring characteristics using regression and stepwise regression models. The results showed that the abrupt temporal changes of temperature and precipitation were different from that of the standard chronology of P. euphratica in the lower reaches of the Tarim River. Moreover, there was no significant relationship between temperature and standard chronology of P. euphratica or between precipitation and standard chronology of P. euphratica in the lower reaches of the Tarim River. The Pearson correlation analysis also revealed no significant relationship between temperature and groundwater depth or between precipitation and groundwater depth in the lower reaches of the Tarim River. However, the width indexes of tree rings of P. euphratica were significantly and negatively related to the groundwater depth, and they were significantly and positively related to the amount of ecologically conveyed water from August to July in the Yingsu section of the lower reaches of the Tarim River. These results proved that historical information regarding water conditions could be accurately reproduced by the tree rings of P. euphratica in the lower reaches of the Tarim River. Therefore, the tree rings of P. euphratica could be used for quantitatively assessing the restoration effects of ecological water conveyance to desert riparian forests and calculating the ecological water demand for natural vegetation in the lower reaches of the Tarim River. The average contribution of ecological water conveyance to the increased tree ring width index was 42.96% in the Yingsu section of the lower reaches of the Tarim River. Estimated from the equation of the tree rings of P. euphratica and the amount of ecologically conveyed water, an annual ecological water conveyance of 0.84×108 m3 was expected to maintain the average growth (i.e., average tree ring index of 0.997 during 1933-2015) of P. euphratica in the Yingsu section within a distance of 300 m to the river channel in the lower reaches of the Tarim River. Furthermore, an annual water conveyance of 0.91×108 m3 was estimated to keep P. euphratica at an average radial growth rate of 0.972 during the period of 1933-1974 (i.e., before the main channel of the Tarim river near Yingsu section dried up).
Key words: Populus euphratica     chronology     ecological water conveyance     ecological water demand    

生态需水量系指为了维持某一环境功能或环境目标所需的水资源量[1]。生态需水量是水资源开发利用过程中必须考虑的一个重要水资源需求组分, 尤其是在水资源极端匮乏的中国西北内陆河流域, 它是脆弱生态区恢复重建的一个关键问题[2-3], 但实际水资源配置中这部分水资源的供给往往被政府部门所忽视[4]。近年来, 随着可持续发展和生态文明建设的倡议, 天然植被生态需水研究也成为了当前生态学、环境保护和水资源管理的热点问题[5-6]。目前生态需水量模型主要基于水量平衡、水热平衡、水盐平衡或水沙平衡, 没有涉及生态系统和水文系统的互作机制, 也未考虑水分与物种生产力之间的关系, 同时这些水文模型采用的数据来源是相互独立的, 由于观测和实验数据的时间序列有限, 这些水平衡模型无法诠释水文过程和生态过程的交互作用[4], 造成很多研究结果相互悖论。由此可见, 寻找长时期、易获取、可重复、同时承载植物生理过程与水文过程信息的代用资料或材料则成为了天然植被生态需水量精准计算的关键。树木是记录区域生态环境变化的绝佳档案[7-8], 树木年轮以定年准确、连续性强、分辨率高和易于获取复本等独特优点, 已成为研究环境变化的重要代用资料[9]。塔里木河流域由于过去几十年来生产用水大量挤占生态用水, 下游河道长期断流, 荒漠河岸林极度衰败, 具有战略意义的绿色走廊岌岌可危[10]。荒漠河岸林的主体是胡杨, 它是荒漠河岸林的唯一乔木和关键建群种, 且胡杨根系具有显著的水力提升能力, 群落中草本植物或浅根灌木生长发育的水分来源也主要依赖于胡杨根系提升的水分再分配[11], 因而胡杨的存亡决定着荒漠河岸林的存亡。由此可见, 在塔里木河下游要估算维持荒漠河岸林的生态需水量, 关键在于正确计算保证胡杨正常生长发育的生态需水量。基于此, 本文通过分析塔里木河下游唯一乔木——胡杨树木年轮长期生长特征, 明确生态输水量与胡杨木质部生长发育的量化关系, 试图采用树木年轮技术来定量评估生态输水对木质部形成的贡献率并计算维系塔里木河下游荒漠河岸林不同生长状态的生态需水量, 旨在为计算干旱区内陆河流域天然植物生态需水量提供新方法。

1 研究区概况

塔里木河下游为大陆性暖温带荒漠干旱气候, 降水稀少, 蒸发强烈, 年太阳总辐射5692—6360 MJ/m2, 年日照时数2780—2980 h, ≥10℃的年积温为4040—4300 ℃, 平均日差13—17 ℃左右;平均年降水量在17.4—42.0 mm之间, 平均年蒸发量(潜势)达2500—3000 mm[12]。沿河两岸以胡杨为主的荒漠河岸林是驰名中外的“绿色走廊”, 对阻止沙漠化进程、降低土壤盐渍化、保护河岸林生态系统的健康发展、维护绿洲的生态安全和持续稳定发展具有举足轻重的作用[13]。然而, 由于水资源利用不合理, 尾闾湖泊-罗布泊和台特玛湖分别于1970和1972年干涸, 下游从大西海子以下至台特玛湖共321 km河道于1972年彻底长期断流。河道断流使得地下水位大幅下降, 由地下水维系的荒漠河岸林大面积退化[13], 天然胡杨林面积从20世纪50年代的5.4×104 hm2减至目前的不足0.67×104 hm2, 减幅高达88%, 使得林间沙地活化、风蚀沙化作用加强, 土地荒漠化过程加剧, 至2000年, 该区流沙地和盐碱地面积已分别高达52.71%和12.35%, 导致塔克拉玛干和库鲁克两大沙漠有合拢趋势, 具有战略意义的绿色走廊岌岌可危[10]。塔里木河下游已成为了中国西部生态环境最脆弱的区域, 由此相继引发的严重生态问题在中国乃至世界范围都是十分典型的, 该区生态恢复与生态安全问题已引起了世界各地和社会各界及政府的高度关注。为拯救绿色走廊, 恢复生态和保护环境, 中国政府于2000年始斥资107亿元组织实施了生态输水工程。

2 材料与方法 2.1 胡杨树木年轮样芯采集

本研究的取样区域位于英苏, 地处塔里木河下游中段(图 1), 2000年始在此断面垂直于河道布设了7眼地下水监测井, 距河道的距离分别为50、150、250、350、450、750、1050 m。由于生态输水对地下水抬升的影响效应随距河道距离的增加而减小[14-15], 地下水响应最敏感的是距离河道350 m以内[14]。同时, 野外调查发现, 距河道300 m内是胡杨分布密度最高的区域。因此, 2015年9月采用生长锥(Haglof increment borers, Haglof Inc., Sweden)钻取以垂直河道300 m为直径的圆形取样区内所有胡杨树木年轮样芯, 每一样株在离树干基部1—1.5 m处取2个树芯样品, 取样时避开树结、分叉及虫蛀点, 两次取样相互垂直。由于胡杨木质部腐朽严重, 对于胸径超过50 cm而木质部腐烂后仅余≤5 cm的过短树芯视为无效取样, 共采集有效树木年轮样株29棵, 合计58个样芯。

图 1 研究区示意图 Fig. 1 The sketch map of study area
2.2 胡杨树木年轮标准年表建立

样芯经干燥、固定和打磨, 采用骨架法(Skeleton)进行交叉定年, 经树木年轮分析仪(LintabTM 6, Germany)测量样芯年轮宽度, 精度为0.001 mm。用Cofecha交叉定年质量控制程序进行交叉定年检验和生长量的订正。年表的建立采用Arstan软件完成。采用负指数函数或线性回归进行去趋势分析, 对于去趋势后与其他序列去趋势后的趋势差别很大的个别序列, 则采用步长大于总序列长度2/3的样条函数进行去趋势, 剔除相关性差、序列过短或奇异点过多的个别序列构建塔里木河下游英苏断面胡杨标准年表。利用SSS指数(Subsample Signal Strength)大于0.85来确定年表的起始年代和年表序列的可靠性。

2.3 气候水文数据收集

气温和降水数据(1958年建站以来至2010年)来自于塔里木河下游中下段的唯一气象站-铁干里克气象站(距离英苏断面约30 km), 2000—2015年生态输水数据来自于塔里木河流域管理局, 地下水埋深数据取自2000始布设的地下水监测井(每2—3个月测量一次), 将取样区内的三眼地下水监测井测得的埋深平均值视为英苏断面地下水埋深。鉴于环境因素对树木年轮宽度指数可能有一定的滞后作用, 本文在计算气温、降水量和生态输水量时将前一年8月至当年7月的平均值/累计值记为年均值/年累计量, 将胡杨年内主要生长期4—7月的平均值/累计值记为生长季平均值/累计值。

2.4 生态输水对胡杨木质部形成的贡献率及生态需水量计算

生态输水对胡杨树木年轮宽度指数的贡献率计算采用以下公式:

(1)

其中, y为生态输水对胡杨年轮宽度指数的贡献率, x1为生态输水期间胡杨年轮宽度指数, 断流期间胡杨多年平均年轮宽度指数, (x1-)为生态输水下胡杨宽度指数的增加量。式中的x是指去趋势后的标准年表宽度指数。

基于树木年轮的植被生态需水量估算方法:将2000—2015年间每年去趋势后胡杨标准年表宽度指数与断流期间去趋势后标准年表多年平均宽度指数之差作为因变量, 表征胡杨年轮宽度指数对生态输水响应的生态效应量, 将生长年内有效输水量作为自变量进行曲线拟合, 建立如下回归方程:

(2)

其中, y表征生态输水对胡杨年轮宽度指数的影响量, x为生长年内的生态需水量(×108 m3)。

2.5 统计分析方法

采用非参数Mann-Kendall单调趋势和Mann-Kandall突变检验方法检验降水和气温序列变化趋势和突变时间;采用滑动t-检验确定胡杨年轮宽度指数序列的突变时间。采用SPSS 13.0软件的单因素方差分析、相关性分析、回归分析、逐步回归分析统计胡杨年轮宽度指数与区域气温、降水、地下水埋深和生态输水量的量化关系。

3 结果与分析 3.1 胡杨树木年轮年表特征

剔除相关性差、序列过短或奇异点过多的10个序列构建了塔里木河下游英苏断面胡杨标准年表(图 2A)。年表统计特征显示(表 1), 胡杨标准年表的相关系数、平均敏感度和EPS都比较高, 表明胡杨对外界环境的敏感度较高, 年表序列包含丰富的低频和高频环境信息。胡杨标准年表的生长趋势显示, 1940s胡杨树木年轮生长较快, 1950s树木年轮宽度呈下降趋势, 1960s呈上升趋势, 1970s—1990s呈下降趋势, 2000后胡杨年轮宽度增加(图 2A)。对胡杨标准年表进行滑动t检验, 子序列长度为10 a时, 胡杨年表有4个突变点, 突变时间分别为1957年、1967年、1975年和2000年(图 2B), 表明1957—1966的10年间, 区域生态环境不利于胡杨生长, 1967—1974年间生态环境有好转趋势, 而1975—1999年之间生态环境再度恶化, 2000年生态环境有所改善。

图 2 塔里木河下游英苏断面胡杨标准年表(A)和标准年表滑动t检验(B) Fig. 2 Standard chronology of Populus euphratica (A) and Moving t-test technique of standard chronology of Populus euphratica (B) in Yingsu of the Lower Tarim River

表 1 塔里木河下游英苏胡杨标准年表统计特征 Table1 Statistical information on the standard chronologies of Populus euphratica in the Yingsu section of Lower Tarim River
年表特征
Characteristics
胡杨标准年表
Standard chronologies
of Populus euphratica
(SCP)
样本量Sampling cores in chronologies48芯/24树
时段Time span1933—2015
相关系数Seriers intercorrelation0.643
平均敏感度Mean sensitivity0.364
标准误Standard deviation0.241
自相关系数First-order autocorrelation0.520
公共区间Common period1955—2015
信噪比Signal-to-noise ratio (SNR)1.19
总体代表性Expressed population signal (EPS)0.945
样本数为最终用于年表建立的胡杨样芯样本量;时段、公共区间为样本统计值;敏感度、相关系数、自相关系数、标准误为去趋势序列的统计值;EPS、SNR为去趋势序列在公共区间内的平均统计值。SCP表示英苏断面垂直河道0—300 m胡杨标准年表; 表中字母含义下同
3.2 胡杨年轮生长与区域环境因素关系

影响树木生长的生态环境主要为气候因素和水文过程。在气候方面, 塔里木河下游1958—2010年生长季气温有显著增加趋势, 突变时间为1991年(图 3A);生长季降水在1958—2010年间并未表现出明显的突变(图 3B), 温度和降水的突变时间与胡杨标准年表的突变时间并不吻合。Pearson相关分析表明, 生长年气温、生长季气温、生长年降水量和生长季降水量对胡杨年轮宽度指数的影响都不显著(图 4);同时, 气温和降水对地下水埋深的影响也未达到显著水平(图 4)。这些结果表明, 塔里木河下游气候因素对地下水埋深和胡杨生长均没有显著影响。

图 3 塔里木河下游温度(A)和降水(B)突变点检验 Fig. 3 Mann-Kendall-Sneyers test (P>0.05) of temperature (A) and precipitation (B) in Lower Tarim River Uf, UB为原时间序列和反序列变化趋势曲线, |Uf, UB|>1.96表示上升趋势明显, |Uf, UB| < 1.96表示下降趋势明显

图 4 胡杨年表与生态环境指标(A)、地下水埋深与生态环境指标(B)及地下水埋深滞后1年与生态环境指标(C)的相关性分析 Fig. 4 Pearson correlations between ecological-environmental factors and standard chronology of Populus euphraticain (A), and grondwater depths (B) as well as grondwater depths on 1-year time lag (C) in Yingsu Section of Lower Tarim Rive 1:生长季气温Average temperature of growth season; 2:生长年气温Average temperature of growth year; 3:生长季降水Average precipitation of growth season; 4:生长年降水Average precipitation of growth year; 5:地下水埋深Average groundwater depth; 6:地下水埋深滞后1年Average groundwater depth on a 1-year time lag; 7:自然年生态输水量Ecological water conveyance; 8:生长年生态输水量Ecological water conveyance of growth year.

1972年后塔里木河下游321 km全面断流, 河岸胡杨林生存发育主要依赖于地下水, 相关分析也证实, 地下水埋深与胡杨年轮宽度指数呈显著负相关(图 4)。然而英苏断面胡杨年轮宽度指数在1975年才开始下降, 表明胡杨生长对河道断流有3年的响应滞后期, 这证实地下水位在断流期间是逐渐下降的, 经过3年时间英苏断面的地下水位就已下降到抑制胡杨生长的水位之下。2000开始, 中国政府启动了拯救塔里木河下游荒漠河岸林的应急生态输水工程, 每年为下游人工泄水, 以抬升地下水位。2000至2015年间共实施了22次以拯救塔里木河下游胡杨林为目标的应急生态输水, 其中2000、2001、2003、2005、2010、2011和2013年年内分别实施了2次输水, 16年间累计实现生态输水51.09×108 m3, 年内最大输水量8.23×108 m3、最小输水量为0.11×108 m3, 14次将水输送到了-尾闾湖泊-台特玛湖。生长年生态输水量与英苏年均地下水埋深具有显著的负相关性(图 4), 表明生长年生态输水有效地促进了地下水埋深的抬升, 显著地改善了塔里木河下游的水环境。然而值得提出的是, 自然年生态输水量与自然年地下水埋深没有显著关系, 而与滞后1年的地下水埋深呈现极显著负相关(图 4), 表明生态输水的实施与地下水埋深的抬升之间存在1年的滞后效应。但胡杨标准年轮表仅与当年的地下水埋深呈显著负相关, 与生长年有效生态输水量呈显著正相关(图 4, P < 0.05), 表明生长年有效生态输水量通过抬升地下水埋深对胡杨木质部径向增加起着显著的正向促进效应。

由此可见, 塔里木河下游胡杨的生长与区域水环境变化过程密切相关, 树木年轮主要承载的是区域水文环境的历史变化信息, 可以作为定量评估生态输水工程的生态恢复效应和定量计算植被生态需水量的新指标。

3.3 生态输水对胡杨年轮宽度指数的贡献率及胡杨林生态需水量计算

为定量辨别生态输水量对胡杨木质部生长的贡献率, 分别统计了2000—2015年生长年生态输水量(图 5), 生态输水期间(2000—2015年)多年平均年轮宽度指数(1.158)和断流期间(1975—1999年)多年平均年轮宽度指数(0.810), 经方程(1)计算显示, 生态输水后胡杨年轮宽度指数较断流期间平均增加了0.348, 即生态输水对塔里木河下游英苏断面胡杨树木年轮宽度增加的平均贡献率为42.96%;以断流期间多年平均年轮宽度指数为未实施生态输水的基准值, 则2000—2015年间生态输水对胡杨木质部形成的贡献率范围在-5.56%—142.96%之间(表 2), 其中, 负值意味着生态输水的持续中断可能会导致胡杨加速衰败, 造成生长速率快速下降, 甚至低于1975—1999年断流期间的平均水平。结合图 5表 2可知, 如遇到特大洪水季或连续洪水年、丰水年, 可根据地下水埋深抬升的滞后效应和树木年轮对河道断流的滞后效应, 采取某些生长年内进行1次或连续特大量的生态输水, 在接下来的一个干旱或枯水生长年内不再进行生态输水, 已抬升的地下水埋深足够满足次年的植物生长发育, 如2012年的生长年内共输入生态水11.15×108 m3后, 尽管2013年生长年输水并不多(仅0.14×108 m3), 但胡杨年轮宽度指数仍保持较高;再如2015年生长年内并未进行生态输水(0×108 m3), 但因2014年生长年有较大量的水量下泄(4.78×108 m3), 胡杨在2015年仍保持着较快的径向生长。然而, 2008—2010年连续中断生态输水后, 使得2009年和2010年胡杨树木径向生长重归缓慢, 甚至低于断流期间的平均生长速率, 表明即使在某一次或连续输水使得地下水埋深得到了较大提升, 后期生态输水的中断期也不能连续超过2年, 否则植被生长又将严重受抑。

图 5 2000—2015塔河下游生长年生态输水量 Fig. 5 Ecological water conveyance of growth year from 2000—2015 in the Lower Tarim River

表 2 生态输水对塔里木河下游胡杨木质部形成的贡献率 Table2 Contribution of ecological water conveyance to ring width of Populus euphratica in the lower reaches of the Tarim River
年份
Year
生态输水后与断流期间
平均年轮宽度指数差值
Difference value of chronology
between ecological water
delivery and dried up periods
生态输水的对木质部
形成的贡献率/%
Contribution rate of
ecological water conveyance
to xylem development
20000.21927.04
20010.60774.94
20020.873107.78
20031.158142.96
20040.79197.65
20050.18522.84
20060.30938.15
20070.30938.15
20080.10112.47
2009-0.045-5.56
2010-0.010-1.23
20110.68884.94
20120.18923.33
20130.0799.75
20140.0121.48
20150.10713.21

以上分析结果表明生长年内有效输水量是保证胡杨木质部生长的关键影响因素。那么, 为维持胡杨正常的生长, 每一生长年内需要从大西海子下泄多少生态水量呢?通过方程(2)计算, 如果要保证塔里木河下游英苏断面胡杨维持多年(1933—2015年)平均径向生长(树木年轮宽度指数0.972), 则每年需下泄生态水量0.84×108 m3;如果要保证塔里木河下游英苏断面胡杨达到断流前(1933—1974年)的平均径向生长(树木年轮宽度指数0.997), 则每年需下泄生态水量0.91×108 m3

4 讨论与结论 4.1 塔里木河下游胡杨树木年轮承载的环境信息

胡杨作为塔里木河下游荒漠河岸林的唯一乔木, 见证了塔里木河下游生态环境变迁的全过程, 是记录生态环境变化的绝佳“档案”, 其生长发育更能直接反应区域生态恢复的效应[16-17]。王振锡等[16]提出塔里木河下游胡杨树木年轮承载着区域水环境历史变迁信息;徐海量等[12]、安红燕等[18]、刘海兰等[19]及Yu等[20]发现生态输水与胡杨木质部径向生长有密切关系;Schilling等[17]也指出胡杨树木年轮宽度可以作为构建地下水-地表水-植物互作模型的有效参数。本文通过突变分析和相关性分析发现, 塔里木河下游温度和降水突变时间与胡杨标准年表突变时间不一致(图 2, 图 3), 气温和降水与年轮宽度指数相关性并不显著(图 4), 且气温和降水与地下水埋深相关性也并不显著(图 4), 但地下水埋深和生态输水量与胡杨年轮年表均呈显著相关(图 4), 这进一步证实塔里木河下游胡杨树木年轮并不适宜于重建和反演长时期的气温和降水, 胡杨木质部最主要承载的是区域水文环境的历史变化信息, 因而, 胡杨树木年轮可以作为评估塔里木河下游长时期内地下水埋深变化和生态输水对植被恢复的生态效应的新指标。

4.2 塔里木河下游胡杨年轮对生态输水的量化响应

生态输水对塔里木河下游胡杨树木年轮宽度增加有积极的促进作用[20-22], 本文研究结果也证实了这一结论。滑动t检验显示, 塔里木河下游英苏断面胡杨年轮标准年表在2000年开始实施生态输水工程后发生了突变, 2000年后年轮宽度指数显著增加。生态输水主要是通过抬升地下水埋深、改善区域水环境来促进胡杨木质部生长的, 但生态输水与地下水抬升存在时滞性[14-15], 本文也发现自然年生态输水与地下水埋深的抬升有1年的滞后效应(图 4), 且胡杨当年年轮宽度与当年自然年生态输水量关系不显著, 仅与生长年有效生态输水量呈显著正相关(图 4), 这与Hao等[23]研究提出的塔里木河下游植被NDVI对地下水埋深的响应有1年滞后结论是基本一致的。因此, 研究认为在分析生态输水的植被恢复生态效应评估时不能用自然年的生态输水量, 而要用生长年有效生态输水量。生长年有效生态输水的实施使得塔里木河下游英苏断面胡杨木质部年均宽度指数较断流区间增加了0.348, 生长年有效生态输水对胡杨年轮宽度指数增加的平均贡献率为42.96%。

4.3 塔里木河下游英苏荒漠河岸林生态输水量的确定

西北干旱区荒漠植被生态需水量的确定一直是荒漠植被生态系统研究的热点[4]。Fu等[24]用地下水储量、水平衡法结合遥感解译的自然植被面积, 计算得出黑河下游自然植被的生态需水量是3.91×108—4.05×108 m3;Zhang[25]用潜水蒸散发模型和地下水平衡模型计算提出为维持黑河下游哨马营断面的自然植被恢复和稳定的地下水埋深, 在2020和2030年每年需要生态需水量6.35×108 m3;在塔里木河下游, Ye等[26-27]通过潜水蒸发模型, 认为2005年塔里木河下游荒漠河岸林自然植被的生态需水量是2.4×108 m3, 同时Ye等[28]通过模拟塔里木河下游生态输水后地下水埋深的变化, 提出维持距河道500、1000、1500、2000 m受损生态系统恢复的需水量分别为1.61×108、2.48×108、3.19×108m3、3.82×108m3。这些生态需水的计算方法主要采用水量平衡模型, 没有涉及生态系统和水文系统的互制机制, 没有考虑水分与物种生产力之间的关系, 缺乏兼顾植物生理、生态和水文物理过程的系统的定量分析方法[4]。另外, 这些生态模型、水文模型或土壤水分模型采用的数据来源是相互独立的, 由于观测和实验数据的有限, 这些水平衡模型无法诠释水文过程和生态过程的交互作用。树木年轮同时承载着植物生理过程与水文过程信息, 为正确计算植物生态需水量提供了新的指标。胡杨是荒漠河岸林的唯一乔木和关键建群种, 胡杨林生态需水量基本能代表荒漠河岸林的生态需水量。通过对胡杨树木年轮在生态输水前后的生长差异与输水量关系的建模, 计算发现如果要保证塔里木河下游英苏断面垂直河道300 m范围内胡杨平均径向生长达到近90年的(1933—2015年)生长速率, 每一生长年需从大西海子为下游下泄生态需水量0.84×108 m3;如果要保证垂直河道300 m范围内胡杨达到断流前(1933—1974年)的平均径向生长, 则每一生长年下泄生态需水量0.91×108 m3。这个数值比Ye等[28]的计算结果小, 除计算方法差异外, 其可能原因还包括前者计算的是维持整个下游荒漠河岸林的生态需水, 而我们的结果仅是维持大西海子到英苏断面荒漠河岸林的生态需水。

研究显示, 胡杨具有显著的根系水力提升能力, 能通过主根和侧根将深层土壤水和地下水提升释放在离胡杨主干4 m范围内的60—120 cm浅土层[11], 这些浅土层的土壤水分正好可以供给草本植物或浅根灌木生长发育。一般而言, 在塔里木河下游, 胡杨生长的适宜生态水位是2—4 m, 胁迫水位是6 m, 限制胡杨生长的极限地下水埋深>8 m[29-31], 而胡杨根系在垂直方向上主要集中在120—170 cm [32], 野外调查发现其主根系甚至可深达6 m以上。同位素研究还显示, 胡杨根系吸水主要来源于375 cm以下的地下水[33]。因而, 根据胡杨根系分布特征和水分吸收特点, 研究建议在地下水埋深降至4—6 m时应及时进行生态输水, 以抬升地下水埋深保证胡杨根系能正常吸收到水分。

综上所述, 本文提出了采用胡杨树木年轮作为定量评估生态输水的植被恢复效应和生态需水量计算的新指标, 为定量计算干旱区断流河道天然植被生态需水量和未来研究植被生理——生态过程与水文过程互作机制提供了一个技术方法。但本文计算结果是基于仅保护大西海子水库到英苏断面距河道300 m范围内胡杨林, 并不是涵括整个塔里木河下游区域, 因此今后还需进一步扩大研究区域范围, 为塔里木河流域水资源管理和生态输水工程的科学实施提供更为详尽和科学的数据支撑。

参考文献
[1] 王根绪, 程国栋. 干旱内陆流域生态需水量及其估算——以黑河流域为例. 中国沙漠, 2002, 22(2): 129–134.
[2] Chen Y N, Li W H, Deng H J, Fang G H, Li Z. Changes in central Asia's water tower:Past, present and future. Scientific Reports, 2016, 6: 35458. DOI:10.1038/srep35458
[3] Chen Y N, Li B F, Li Z, Li W H. Water resource formation and conversion and water security in arid region of Northwest China. Journal of Geographical Sciences, 2016, 26(7): 939–952. DOI:10.1007/s11442-016-1308-x
[4] Chen Y N, Xu C C, Chen Y P, Liu Y B, Li W H. Progress, challenges and prospects of eco-hydrological studies in the Tarim River Basin of Xinjiang, China. Environmental Management, 2013, 51(1): 138–153. DOI:10.1007/s00267-012-9823-8
[5] Ling H B, Zhang P, Xu H L, Zhang G P. Determining the ecological water allocation in a hyper-arid catchment with increasing competition for water resources. Global and Planetary Change, 2016, 145: 143–152. DOI:10.1016/j.gloplacha.2016.08.012
[6] Amr A I, Kamel S, El Gohary G, Hamhaber J. Water as an ecological factor for a sustainable campus landscape. Procedia-Social and Behavioral Sciences, 2016, 6: 181–193.
[7] Farrer E C, Ashton I W, Spasojevic M J, Fu S Y, Gonzalez D J X, Suding K N. Indirect effects of global change accumulate to alter plant diversity but not ecosystem function in alpine tundra. Journal of Ecology, 2015, 103(2): 351–360. DOI:10.1111/1365-2745.12363
[8] Farrer E C, Ashton I W, Knape J, Suding K N. Separating direct and indirect effects of global change:a population dynamic modeling approach using readily available field data. Global Change Biology, 2014, 20(4): 1238–1250. DOI:10.1111/gcb.2014.20.issue-4
[9] Linderholm H W, Niklasson M, Molin T. Summer moisture variability in east central Sweden since the mid-eighteenth century recorded in tree rings. Geografiska Annaler, Series A:Physical Geography, 2004, 86(3): 277–287. DOI:10.1111/j.0435-3676.2004.00231.x
[10] 李卫红, 郝兴明, 覃新闻, 陈亚宁, 黄湘. 干旱区内陆河流域荒漠河岸林群落生态过程与水文机制研究. 中国沙漠, 2008, 28(6): 1113–1117.
[11] Hao X M, Chen Y N, Li W H, Guo B, Zhao R F. Hydraulic lift in Populus euphratica Oliv. from the desert riparian vegetation of the Tarim River Basin. Journal of Arid Environments, 2010, 74(8): 905–911. DOI:10.1016/j.jaridenv.2010.01.005
[12] 徐海量, 邓晓雅, 赵新风. 河道断流对胡杨(Populus euphratica)径向生长量的影响. 中国沙漠, 2013, 33(3): 731–736. DOI:10.7522/j.issn.1000-694X.2013.00105
[13] Chen Y N, Li W H, Xu C C, Ye Z X, Chen Y P. Desert riparian vegetation and groundwater in the lower reaches of the Tarim River basin. Environmental Earth Sciences, 2015, 73(2): 547–558. DOI:10.1007/s12665-013-3002-y
[14] 邓正波. 塔里木河下游应急输水的地下水响应与生态需水研究[D]. 南京: 河海大学, 2005. http://cdmd.cnki.com.cn/Article/CDMD-10294-2005042717.htm
[15] 白元, 徐海量, 张青青, 叶茂. 基于地下水恢复的塔里木河下游生态需水量估算. 生态学报, 2015, 35(3): 630–640.
[16] 王振锡, 潘存德, 石鑫鑫. 胡杨年轮记录的塔里木河下游54年来区域水环境历史变迁. 生态环境学报, 2010, 19(3): 637–645.
[17] Schilling O S, Doherty J, Kinzelbach W, Wang H, Yang P N, Brunner P. Using tree ring data as a proxy for transpiration to reduce predictive uncertainty of a model simulating groundwater-surface water-vegetation interactions. Journal of Hydrology, 2014, 519: 2258–2271. DOI:10.1016/j.jhydrol.2014.08.063
[18] 安红燕, 徐海量, 叶茂, 禹朴家, 龚君君. 塔里木河下游生态输水后胡杨径向生长量的时空变化. 应用生态学报, 2011, 22(1): 29–34.
[19] 刘海兰. 塔里木河不同断面胡杨径向生长对水分响应的研究[D]. 乌鲁木齐: 新疆师范大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10762-1015951287.htm
[20] Yu P J, Xu H L, Ye M, Liu S W, Gong J J, An H Y, Fu J Y. Effects of ecological water conveyance on the ring increments of Populus euphratica in the lower reaches of Tarim River. Journal of Forest Research, 2012, 17(5): 413–420. DOI:10.1007/s10310-011-0312-3
[21] 孙卫, 侯平, 李霞. 胡杨枝直径生长对塔里木河下游应急输水的响应. 生态学报, 2006, 26(9): 2913–2917.
[22] Deng X Y, Xu H L, Ye M, Li B L, Fu J Y, Yang Z F. Impact of long-term zero-flow and ecological water conveyance on the radial increment of Populus euphratica in the lower reaches of the Tarim River, Xinjiang, China. Regional Environmental Change, 2015, 15(1): 13–23. DOI:10.1007/s10113-014-0603-2
[23] Hao X M, Li W H. Impacts of ecological water conveyance on groundwater dynamics and vegetation recovery in the lower reaches of the Tarim River in northwest China. Environmental Monitoring and Assessment, 2014, 186(11): 7605–7616. DOI:10.1007/s10661-014-3952-x
[24] Fu X F, He H M, Jiang X H, Yang S T, Wang G Q. Natural ecological water demand in the lower Heihe River. Frontiers of Environmental Science & Engineering in China, 2008, 2(1): 63–68.
[25] Zhang Z H. Tree-rings, a key ecological indicator of environment and climate change. Ecological Indicators, 2015, 51: 107–116. DOI:10.1016/j.ecolind.2014.07.042
[26] Ye Z X, Chen Y N, Li W H. Ecological water demand of natural vegetation in the lower Tarim River. Journal of Geographical Sciences, 2010, 20(2): 261–272. DOI:10.1007/s11442-010-0261-3
[27] Ye Z X, Chen Y N, Li W H, Yan Y, Wan J H. Groundwater fluctuations induced by ecological water conveyance in the lower Tarim River, Xinjiang, China. Journal of Arid Environments, 2009, 73(2): 726–732.
[28] Ye Z X, Chen Y N, Li W H, Yan Y. Effect of the ecological water conveyance project on environment in the Lower Tarim River, Xinjiang, China. Environmental Monitoring and Assessment, 2009, 149(1/4): 9–17.
[29] Chen Y N, Wang Q, Li W H, Ruan X, Chen Y P, Zhang L H. Rational groundwater table indicated by the eco-physiological parameters of the vegetation:a case study of ecological restoration in the lower reaches of the Tarim River. Chinese Science Bulletin, 2006, 51(S1): 8–15. DOI:10.1007/s11434-006-8202-3
[30] Chen Y P, Chen Y N, Xu C C, Li W H. Groundwater depth affects the daily course of gas exchange parameters of Populus euphratica in arid areas. Environmental Earth Sciences, 2012, 66(2): 433–440. DOI:10.1007/s12665-011-1250-2
[31] 郝兴明, 李卫红, 陈亚宁. 新疆塔里木河下游荒漠河岸(林)植被合理生态水位. 植物生态学报, 2008, 32(4): 838–847.
[32] 叶茂, 徐海量, 王晓峰, 申瑞新. 塔里木河下游阿拉干断面胡杨根系空间分布规律研究. 西北植物学报, 2011, 31(4): 801–807.
[33] Chen Y P, Chen Y N, Xu C C, Li W H. The effects of groundwater depth on water uptake of Populus euphratica and Tamarix ramosissima in the hyperarid region of Northwestern China. Environmental Science and Pollution Research, 2016, 23(17): 17404–17412. DOI:10.1007/s11356-016-6914-8