中国西北干旱区植被绿化可能加剧水资源短缺风险

《Journal of Hydrology: Regional Studies》：Vegetation greening in arid Northwest China may exacerbate water scarcity risks

　　

植被通过蒸腾作用（transpiration）在调节水循环中发挥着关键作用。因此，理解植被如何与陆地水储量（terrestrial water storage, TWS）在空间、时间和土壤深度上进行相互作用，对于制定干旱区可持续水资源策略至关重要。一个主要挑战在于，干旱区植被-水分关系具有高度异质性、非线性且常存在滞后效应，这使得识别绿化何时缓解或加剧水分胁迫变得困难。尽管如此，这些地区植被与TWS之间时空异质性且相互耦合的关系仍知之甚少。研究人员在此阐明了中国西北干旱区（Northwest China, NWC）水资源与植被动态耦合的机制及阈值行为。研究发现，在历史时期（1981–2023年），海拔2000米以上区域表现出更高的陆地水储量异常（terrestrial water storage anomalies, TWSA），伴随广泛的植被绿化。相比之下，塔里木盆地和柴达木盆地等低洼盆地呈现结构性干旱（总土壤水分接近零且植被覆盖极少）。值得注意的是，研究人员识别了植被对不同深度总土壤水分（total soil moisture, TSM）的非线性响应：（1）初期植被绿化导致浅层土壤水分（shallow soil moisture, SM）快速耗竭；（2）超过关键阈值（NDVI ≈ 0.133，LAI ≈ 4.667）后的持续绿化依赖于降水和融雪对根区土壤水分的滞后补给。在1.5–4°C（2024–2099年）的变暖情景下，植被趋势预计将在37°N附近出现逆转。具体而言，在3°C变暖情景下，NWC北部由绿化转为退化，而南部由干旱转为恢复。这种空间逆转与降水格局改变密切相关：北疆和黄土高原面临干旱加剧和植被损失，而南疆和青藏高原边缘则因湿度增加和植被再生而受益。通过明确量化深度依赖性阈值、滞后补给效应及未来植被响应的空间逆转，该研究为界定干旱区生态恢复的水资源限制提供了基于过程的科学基础。

中国西北干旱区（Northwest China, NWC）是全球典型的水分约束区域，其植被-水文耦合机制对理解气候变化下的生态系统响应具有重要科学价值。该地区地形从高山到内陆盆地呈显著梯度变化，水资源主要依赖高山冰雪融水补给。近年来，全球变暖和大气蒸发需求增加加剧了干旱区的水分胁迫，而大规模生态恢复工程（如三北防护林工程）又改变了区域水循环格局。然而，植被与陆地水储量（terrestrial water storage, TWS）、根区土壤水分（root-zone soil moisture）之间的复杂相互作用，以及这些相互作用如何随深度和时间变化，仍缺乏系统的阈值框架进行研究。为此，研究人员整合多源遥感与再分析数据，基于1981至2023年的长时序记录，采用广义可加模型（generalized additive models, GAMs）、格兰杰因果检验（Granger causality tests）及偏最小二乘回归（partial least squares regression, PLSR）等方法，建立了植被变化与深度依赖性土壤水分约束、滞后补给过程及未来变暖轨迹相联系的阈值框架，以期为干旱区生态恢复的水资源承载力上限提供定量依据。

研究使用的主要关键技术方法包括：基于季节趋势分解（Seasonal and Trend decomposition using Loess, STL）的时间序列预处理技术，用于分离长期趋势、季节周期和残差成分；滑动窗口最大相关分析方法，用于确定植被指数与水分变量之间的最优滞后时间和最大相关系数；广义可加模型，用于捕捉植被-水分关系的非线性特征；分段线性回归（piecewise linear regression），用于识别植被响应的关键阈值及其95%置信区间；双向格兰杰因果检验，用于评估植被与水分变量之间的时序因果关系和主导方向；偏最小二乘回归，用于量化多气候驱动因素对陆地水储量变化的相对贡献；基于CMIP6多模式集合模拟的全球变暖水平阈值判定方法，用于识别1.5°C、2°C、3°C和4°C变暖情景下对应的模拟时段和未来预估。样本数据来源于ERA5再分析资料（0.25°×0.25°分辨率）、AVHRR卫星遥感植被数据集（~0.05°分辨率）、Li等人重建的TWSA数据集（0.5°×0.5°分辨率），以及8个全球气候模式的SSP情景模拟输出。

研究结果部分包含以下核心发现：

**2.1. 气候变暖驱动山地生态系统显著绿化**。1981至2023年间，NWC植被指数呈现显著的空间异质性趋势。高海拔山地如祁连山、阿尔泰山和天山北坡持续绿化， fraction of absorbed photosynthetically active radiation（FAPAR）、总初级生产力（gross primary productivity, GPP）和净初级生产力（net primary productivity, NPP）稳步增长；而塔里木盆地、柴达木盆地和河西走廊等内陆热盆地则出现归一化植被指数（normalized difference vegetation index, NDVI）和叶面积指数（leaf area index, LAI）显著下降。时间序列分析表明，GPP、NPP和FAPAR自1980年代以来持续上升，2000年后增速加快，而LAI和NDVI增幅微弱。这种"生产力上升-绿化停滞"的分离现象归因于：持续升温和大气干旱化抑制叶片扩展；CO₂施肥效应、增温诱导的春季物候提前及地形抬升带来的融水输入共同提升单位叶面积光合效率；以及人类活动在东部的退耕还林、灌溉等措施与盆地地下水超采的差异化影响。土壤水分四层分析揭示，高海拔区域各层水分充足（约0.682 m³ m^-3），而塔里木和柴达木盆地自表层至根区均接近零值，形成深层结构性干旱。

**2.2. 根区水分和外部补给驱动植被-生态生产力平衡**。最大相关和最优滞后分析表明，深层土壤水分（swvl3–swvl4）和TWSA与植被的季节至年际关联强于浅层水分，植被动态更受深层和综合水储量调控而非近地表短期波动。双向滞后分析揭示：浅层土壤水分主要参与快速短期相互作用，深层土壤水分和TWSA主导更强且滞后更长的双向耦合；植被生产力指标（FAPAR、GPP、NPP）对水分可用性的敏感性强于冠层结构指标（NDVI、LAI）。GAM分析显示生产力相关指标对TWSA和深层土壤水分的解释力更高（最大R²=0.14），但仍存在大量未解释变异。分段回归识别出关键阈值行为：低至中等植被水平时，植被生长主要耗水，水储量随植被增加而下降；超过阈值后斜率减弱甚至转正，表明水-植被协同变化。该阈值行为普遍存在于TWSA和四个土壤层中，浅层土壤水分阈值前斜率最陡、响应快速且易受耗竭，深层土壤水分和TWSA则显示更平滑持久的过渡。格兰杰因果检验证实水-植被耦合的不对称性：水资源对植被的影响强于反向影响，且随土壤深度增加滞后延长；植被对水资源的因果影响主要限于浅层土壤水分（滞后小于4个月），对深层水和TWSA的影响因累积传输周期长而较弱。

**2.3. 地表水分流失与深层补给的双通道机制**。1981至2023年间，气温以每十年0.456°C的速率上升。蒸发与TWSA及浅层土壤水分呈显著负相关（r约-0.145至-0.168），降水亦呈负相关（约-0.15至-0.20），表明大气水分需求增加快速耗竭近地表土壤水。随深度增加，温度效应的最优滞后期缩短，三类气候变量对深层土壤水分的影响转为正相关（降水与swvl4的相关系数约0.15）。GAM的R²值深层高于浅层，反映深层土壤主要通过季节至年尺度的渗透等积分过程响应气候。PLSR分析表明，对TWSA和浅层土壤水分而言，降水是最强负向驱动因子（系数-0.115至-0.218），蒸发次之（-0.114至-0.183）；深层土壤水分（swvl3–swvl4）中，蒸发成为主导正向贡献因子（0.132–0.143），温度呈中等正向影响。偏相关分析显示，在控制植被效应后，蒸发的标准化系数均值为-0.54，降水为-0.38。这些发现揭示NWC水资源受两个相互关联机制调控：快速的近地表耗竭通道和延迟的深层-流域尺度补给通道，后者由降水、融雪和冰川融水渗透驱动。

**2.4. 植被生长加剧水资源短缺风险**。基于CMIP6多模式集合的未来预估显示，随全球变暖从1.5°C升至4°C，植被指数和上层土壤水分（mrsos）的时序趋势增强，相关性由显著负转正。空间上，LAI、GPP和NPP在37°N附近发生格局逆转：从北增南减转为北减南增。低变暖情景（1.5–2°C）下，祁连山、阿尔泰山、天山和青藏高原植被指数下降而新疆盆地上升，mrsos在大部分区域显著增加，这是因为升温抬升雪线、加速融雪暂时提高地表水分，随后被增强的蒸散发和饱和水汽压差（vapor pressure deficit, VPD）耗竭。3°C变暖时，下降区域急剧收缩，仅东天山保持显著下降，西藏和黄土高原出现广泛增长。4°C变暖时，北部盆地和戈壁地区VPD更高、热浪更频繁，植物维持呼吸和微生物呼吸增加，GPP和NPP由增长转为停滞或下降，植物采取保守策略导致LAI受限，下降区域扩大3.63倍。各指标对变暖的敏感性差异显著：LAI在3–4°C时增加，GPP和NPP在中低变暖下可通过生理调节维持，但强变暖下在东天山及其以南、准噶尔盆地等区域下降；mrsos虽在全变暖水平下以增加为主，但增速随变暖增强而减缓。

讨论部分，研究人员将NWC的发现置于全球干旱区变化背景中进行分析。自1980年代以来，全球干旱区普遍经历气温上升、降水减少和蒸发增加，导致陆地水储量、表层土壤水分（surface soil moisture, SSM）和根区土壤水分（RZSM）广泛持续下降。与全球其他干旱区类似，NWC呈现出绿化与水资源下降并存的现象——全球约40%–50%的绿化伴随TSM下降，主因是植被扩展增强蒸腾耗水。NWC高海拔山区的持续绿化主要由融雪和冰川融水补给维持，而非仅依赖人工造林；低海拔盆地则因缺乏外部水源而深陷结构性干旱。研究提出的阈值框架具有全球意义：低于阈值时植被绿化以消耗水资源为代价，超过阈值后需稳定的高山融水补给和节水措施支撑。这与西非的"水文悖论"、萨赫勒地区的降雨-径流耦合变化等现象相呼应。研究强调，未来NWC生态恢复的可持续性不能仅凭绿化判断，即使三北防护林工程持续至2050年，若变暖驱动的蒸散发超过有效补给，植被收益将难以维持。对比地中海、加利福尼亚和中亚等地区的经验，干旱区恢复策略应优先降低大气水分需求、保障高海拔水源稳定输入、增强根区水储存能力，并将水资源承载力阈值作为生态恢复的上限依据。

研究结论指出：1981–2023年间，NWC高海拔山区植被活动持续增加，盆地则近乎零增长，2000米以上区域绿化与水可用性共存，而塔里木和柴达木盆地持续水分受限。植被-水分耦合具有强烈非线性和尺度依赖性：初期绿化伴随浅层土壤水分快速耗竭，持续生长依赖深层土壤水分和陆地水储量的延迟补给。PLSR分析表明气候变化与表层土壤水分负相关，而深层土壤水分与降水和温度在年际尺度上正相关。未来1.5°C至4°C变暖下，37°N附近的格局逆转将NWC分为相反轨迹——北部由绿转退、南部由干转复。尤其在4°C情景下，北部盆地GPP和NPP比LAI下降更快。总体而言，NWC的植被绿化并非总是生态改善的标志，而可能加剧有限水资源的压力，这取决于海拔、土壤水分深度和降水变化。未来工作需进一步区分气候变化、生态恢复、灌溉和地下水利用对植被-水分相互作用的相对贡献，并纳入更高分辨率观测和过程模型以更好约束持续变暖下的阈值行为和水承载极限。