当绿化过度时：生态恢复削弱了黄土高原的水文抵抗力

《Forest Ecosystems》：When greening goes too far: Ecological restoration undermines hydrological resistance on the Loess Plateau

　　

大规模生态恢复日益被提倡作为增强景观恢复力的策略，但其水文后果，特别是在加剧的干旱和大气干燥度下，仍然未被充分理解。聚焦于中国的黄土高原，那里退耕还林工程（GGP）推动了广泛的植被恢复，研究人员开发了一个多方法框架，该框架整合了一个干旱–水指标——水消耗阈值（WDT），与因果机器学习和路径分析，以量化植被绿化如何改变对干旱的水文抵抗力。研究结果表明，密集的植被绿化削弱了径流和土壤水的干旱抵抗力。径流的平均处理效应（ATE）为0.026，土壤水为0.055，绿化引起的影响分别发生在研究区域的53.6%和70.6%。在空间上，最强的削弱集中在半湿润地区，那里植被恢复最为显著。路径分析进一步表明，绿化主要通过直接和蒸散发介导的途径削弱水文抵抗力。随着气候极端事件加剧，恢复工作必须在植被目标与长期水可持续性之间取得平衡。

**论文解读：生态恢复对黄土高原水文抵抗力的影响机制**

**研究背景与问题**
全球范围内，大规模生态恢复被广泛视为提升景观恢复力的关键战略，如《波恩挑战》和联合国生态系统恢复十年等倡议旨在恢复数亿公顷退化土地。然而，植被恢复不可避免地增加水需求，尤其在干旱和半干旱地区，可能加剧水文变异性。随着气候极端事件（如干旱和大气干燥度增强）的加剧，大规模绿化引发的水文权衡尚不明确，亟需系统评估。黄土高原作为东亚最脆弱的生态区之一，其深厚的黄土层、半干旱气候和长期水土流失历史，为研究提供了理想场所。中国自1999年实施的退耕还林工程（GGP）显著提升了植被盖度，带来了碳汇增强和土壤保持等生态效益，但同时也引发了区域水安全担忧——土壤水分以每年3.18 mm的速度下降，地下水储量减少高达60%。现有研究多聚焦长期趋势或平均状态，忽略了干旱胁迫下植被如何调节水动态这一关键缺口，限制了评估生态恢复在气候极端事件下是增强还是削弱水文抵抗力的能力。

**研究目的与意义**
本研究旨在量化GGP前后黄土高原植被变化对水文抵抗力的影响，评估干旱条件下水文抵抗力的时空格局，并揭示植被绿化削弱水文抵抗力的机制。该研究发表在《Forest Ecosystems》上，为平衡生态恢复目标与长期水可持续性提供了科学依据，强调了过度绿化可能加速干旱期水耗竭的风险。

**主要技术方法**
研究人员采用多方法集成框架，主要关键技术包括：①基于水消耗阈值（WDT）的水文响应建模，通过分段线性回归和自举重采样识别干旱强度与径流/土壤水异常的非线性关系；②因果机器学习框架，使用双稳健学习器（DR-learner）结合倾向性评分加权和结果回归建模，以GGP后叶面积指数（LAI）显著上升的像元为处理组，其他为对照组，估计植被绿化的平均处理效应（ATE）和条件平均处理效应（CATE）；③路径分析，将蒸散发和降水作为中介变量，检验绿化对WDT的直接影响和间接效应。数据来源包括GIMMS LAI4g（1982–2020）和MODIS LAI（2021–2024）融合的植被指标、ERA5-Land再分析数据（径流、土壤水、蒸散发等）、SPEIbase v2.11干旱指数，以及NASADEM和CHILI地形数据。

**研究结果**

**3.1 植被绿化的时空动态**
GGP实施后，黄土高原65.06%的区域绿化速率超过0.001 m²·m^?2·year^?1，区域平均LAI上升趋势从0.003增至0.017 m²·m^?2·year^?1。气候带分析显示，从干旱区到半湿润区，绿化程度逐渐增强，半湿润区年均LAI变化达0.021 m²·m^?2·year^?1，显著高于干旱区（0.006）和半干旱区（0.013），Kruskal-Wallis检验表明差异显著（p < 0.01）。

**3.2 径流与土壤水及水文抵抗力的变化**
GGP后，土壤水出现显著下降趋势（斜率 ?2.298 mm·year^?1，p < 0.05），且超过76%的区域土壤水显著减少；径流趋势不显著。WDT变化（ΔWDT）在空间上呈现东南–西北对比：东南部ΔWDT升高（抵抗力下降），占研究区59.16%（径流）和46.91%（土壤水）；西北部ΔWDT降低（抵抗力增强）。半湿润区平均ΔWDT呈上升趋势，而干旱区呈下降趋势，表明绿化与气候带共同影响抵抗力重分布。

**3.3 植被绿化对水文抵抗力的因果效应**
因果推断表明，植被绿化对径流和土壤水的WDT产生了正向平均处理效应，ATE分别为0.026和0.055，置信区间均大于零，说明绿化普遍削弱了水文抵抗力。空间上，CATE正向值覆盖研究区53.6%（径流）和70.6%（土壤水），主要分布在东南和中部半湿润区；西北干旱区则呈现负向CATE（抵抗力增强）。气候带差异显著，半湿润区CATE最高，表明绿化对抵抗力的削弱最强烈。

**3.4 植被绿化影响水文抵抗力的机制路径**
路径分析模型拟合良好（RMSEA = 0.009，SRMR = 0.003）。绿化对径流WDT和土壤水WDT的直接效应显著，标准化路径系数分别为0.30和0.37。绿化还通过增强蒸散发（系数0.59）间接增大WDT，蒸散发对土壤水WDT的间接效应（0.195）强于对径流（0.035）。降水路径不显著，表明水汽压亏缺（VPD）虽负向影响蒸散发和降水，但降水变化并非主导WDT的途径。因此，水文抵抗力削弱主要归因于绿化的直接效应和蒸散发介导的间接效应。

**讨论与结论**
讨论部分指出，大规模绿化导致的水文补偿机制不足以抵消干旱期的水耗竭。尽管前期研究显示绿化可能通过陆–气反馈增加降水，但本研究通过WDT量化发现，干旱期植被高水需求会侵蚀水文缓冲能力，尤其是在半湿润区。因果模型稳健性通过十折交叉验证和SHAP分析得到确认，径流和土壤水模型的ATE在所有折中均保持正值。研究局限性包括：ERA5-Land再分析数据可能存在系统性偏差，且空间因果推断中未考虑相邻像元间的溢出效应。未来需结合原位观测、空间显式因果模型及网络化处理定义来改进评估。

**结论**
本研究表明，退耕还林工程相关的植被绿化削弱了黄土高原的水文干旱抵抗力。因果推断结果显示，植被绿化使径流和土壤水的WDT分别增加0.026和0.055，CATE正向区域分别占53.6%和70.6%，土壤水受影响更广泛。削弱效应在植被绿化最显著的半湿润区最强。路径分析表明，水文抵抗力下降主要与植被诱导的蒸散发增强有关，而降水途径不显著。这些发现表明，黄土高原的生态恢复产生了水文权衡：虽然植被恢复改善了生态条件，但过度绿化会加速干旱期的水耗竭。因此，未来的恢复工作应超越单纯追求最大植被覆盖，转而考虑区域水承载能力、气候带差异和适宜的植被密度。