绿洲化使干旱内流盆地产生响应补偿型恢复力

《Earth's Future》：Oasisification Yields Response-Compensated Resilience in an Arid Endorheic Basin

　　

快速绿洲化被广泛用于加强干旱内流盆地的区域生计和生态安全，但其在水资源日益短缺和生态权衡的背景下进行。这种“绿化”是否真正增强了生态系统恢复力——以及通过哪些组成部分——仍不清楚。在此，研究人员将压力-状态-响应（PSR）核算与可解释机器学习（XGBoost–TreeSHAP）相结合，以刻画塔里木河流域（2000–2020年）不同景观中压力、生态系统状态、响应能力和综合恢复力（ER）的协同变化。绿洲覆盖面积从12.06%扩大至14.32%，流域尺度ER从0.50小幅上升至0.54。然而，这一增益与全流域86.9%的区域响应组分的改善同时发生，而压力在75.6%的区域上升——最显著的是干流走廊（98.8%）——并且仅有65.8%的区域呈现净生态改善。因此，综合恢复力似乎表现为响应补偿型而非压力缓解型。TreeSHAP归因揭示了不同的预测因子结构：净初级生产力（NPP）是绿洲模型中的主要贡献因子（26%–43%），而蒸发需求和水分相关变量在非绿洲模型中贡献更大。PSR象限分析进一步定位了干流沿线的高压、低状态条件。这些结果表明，综合恢复力的上升可能掩盖限水内陆河流系统中持续存在的压力，且需要协调监测压力、状态和响应能力，以支持在气候和人类强迫加剧下的水土资源综合管理。

干旱内流河流域是封闭的内陆水系，其径流主要产生于周边山地的降水和冰雪融水，流经山前冲积扇后大量消耗于灌溉引水和蒸散发。内流河流域约占全球陆地面积的五分之一，却集中了近一半的缺水地区，水量平衡对气候强迫和人类调控高度敏感。自工业革命以来，人为变暖、用水强度增加与绿洲快速扩张相互叠加，导致尾闾湖萎缩、生境退化、盐渍化加剧等问题。全球灌溉面积从1961年的1.47×10⁸公顷增至2020年的3.43×10⁸公顷，绿洲扩张往往伴随下游生态流量减少和湖泊退化。生态系统恢复力（ER）被定义为生态系统在干扰下维持关键结构与功能的能力，涵盖抵抗、吸收、恢复和长期适应。然而，在干旱内流河流域，绿洲化过程中压力、状态与响应如何协同变化，以及不同景观类型中ER的驱动机制是否不同，仍不清楚。

为此，研究人员以塔里木河流域为研究对象，该流域是中国最大的内流河流域，具有典型的山地-绿洲-荒漠镶嵌格局，并经历了快速的绿洲扩张。研究集成2000–2020年多源栅格数据，构建了基于压力-状态-响应（PSR）框架的综合恢复力指数，并耦合XGBoost回归与TreeSHAP可解释机器学习，系统分析了ER的时空格局、组分关系及驱动因子。研究发现，研究期间绿洲覆盖面积从12.06%增至14.32%，流域尺度ER从0.50上升到0.54，但全流域75.6%的区域压力上升，86.9%的区域响应改善，仅65.8%的区域呈现净生态改善。这说明综合恢复力的提升主要来自响应补偿而非压力缓解，表现为“响应补偿型”恢复力。该研究为水资源约束下干旱内流盆地的绿洲开发与生态安全平衡提供了可量化、可比较的依据，成果发表在《Earth's Future》。

在方法上，研究人员整合了归一化植被指数（NDVI）、植被覆盖度（FVC）、太阳诱导叶绿素荧光（SIF）、净初级生产力（NPP）、降水（PRE）、气温（TEM）、潜在蒸散（POT）、蒸气压亏缺（VPD）、土壤湿度（SM）、帕尔默干旱强度指数（PDSI）、人类足迹指数（HFP）等多源数据，统一重采样至1千米网格。基于InVEST模型计算产水量、土壤保持、生境质量和碳储存等生态系统服务，运用熵权法确定指标权重，分别计算压力（P）、状态（S）和响应（R）维度得分，并通过加权求和得到综合ER指数。驱动分析中，分别对绿洲和非绿洲区域建立XGBoost回归模型，采用时空分块交叉验证和TreeSHAP进行特征归因，以识别各预测因子的贡献方向和强度。

结果如下：

3.1 塔里木河流域绿洲扩张动态：研究人员分析了九个子流域2000–2020年绿洲面积比例，发现所有子流域绿洲均呈扩张趋势，但程度