《Earth's Future》：Intensifying Flash Droughts Aggravate Vegetation Ecosystem Vulnerability in the Monsoon Region of China: Evidence From Resistance and Resilience Responses
编辑推荐：
植被的抵抗力（resistance）与恢复力（resilience）决定了其在骤发干旱（flash drought）胁迫下维持生态系统功能稳定性的能力。然而，量化植被抵抗力、恢复力与脆弱性（vulnerability）事件尺度响应特征的研究框架仍显不足。为填补
植被的抵抗力（resistance）与恢复力（resilience）决定了其在骤发干旱（flash drought）胁迫下维持生态系统功能稳定性的能力。然而，量化植被抵抗力、恢复力与脆弱性（vulnerability）事件尺度响应特征的研究框架仍显不足。为填补这一空白，研究人员开发了一种高时间分辨率框架，以评估植被对骤发干旱多维度属性的响应，量化了1982至2019年中国植被抵抗力、恢复力与脆弱性的时空响应特征。分析表明，过去四十年间，中国季风区植被对骤发干旱的脆弱性以每年0.019的速率增加。一方面，这主要由于骤发干旱胁迫下植被生态系统抵抗力的平均下降速率（每年0.074）显著高于其恢复力的平均增加速率（每年0.010）。另一方面，植被生态系统抵抗力与恢复力之间的权衡（trade-off）关系发生了转变，尤以长江中下游流域为甚。在季风与非季风区，可解释机器学习模型表明，温度是影响骤发干旱事件中植被抵抗力与恢复过程的最显著胁迫因子，SHAP值对整体模型解释的贡献介于20%至50%之间。这些发现凸显了在当前全球变暖背景下，关注日益频繁的骤发干旱对陆地生态系统碳汇功能稳定性影响的迫切需求。
研究背景与立项依据：在全球气候变化与人为活动驱动下，干旱事件的复杂性持续增强，其中骤发干旱（flash drought）因次季节尺度上的快速爆发而成为全球干旱监测与预警的关键挑战。相较于传统季节性干旱，骤发干旱对作物产量与植被生态系统总初级生产力（GPP）影响更剧烈，可通过土壤湿度快速耗竭抑制气孔导度、降低碳同化效率，甚至削弱陆地生态系统碳汇功能。已有研究多聚焦骤发干旱的演变特征、形成机制与生态影响，但针对生态系统在骤发干旱推进过程中抵抗力机制、恢复过程与脆弱性响应的系统化定量研究仍有限，尤其缺乏在一致分析结构下逐事件量化抵抗力、恢复力与脆弱性的框架。为此，研究人员以1982至2019年中国不同农业气候区为对象，利用高时间分辨率GPP数据构建多属性评估框架，揭示季风区植被脆弱性上升的驱动机制。
主要技术方法概览：研究人员采用BEPS模型模拟的日尺度GPP（1982至2019年）作为核心生产力指标，并以EC-LUE与GLASS的8天尺度GPP数据交叉验证；土壤湿度识别基于GLDAS-Noah根区（0至100 cm）数据，经三重比对与SMCI1.0地面站基准校验。骤发干旱识别采用全生命周期骤发干旱识别框架（FLCFDIF），划分发展、持续、恢复与终止四阶段。植被响应量化依托“双生命周期曲线”（土壤湿度百分位与标准化GPP异常），定义抵抗力为抵抗期GPP异常下降率之倒数、恢复力为恢复期GPP异常回升率、脆弱性为二者倒数复合指标，并以Pearson相关系数判识抵抗力与恢复力权衡或协同。驱动因子解析集成CatBoost等树基机器学习与SHAP值解释，以及地理探测器交互检测，变量涵盖气象、地表通量（感热SSHF、潜热SLHF）及陆气耦合指数（SM_LE、LE_LCL）。
研究结果：
3.1 植被生态系统抵抗力、恢复力与脆弱性的三维响应关系：基于像素与事件尺度核算显示，中国季风区（长江中下游、云贵高原东部、东北东部沿海）为骤发干旱风险热点，抵抗力与恢复力呈空间逆向分布——热点区抵抗力低（长江中下游多年均3.3±2.0）而恢复力高（1.3±0.2）。季风区脆弱性（0.7±0.6）显著高于非季风区（0.5±0.6），高脆弱区对应低抵抗高恢复格局。极旱区抵抗力最低但恢复力最高；高脆弱集中于干旱指数0.7至1.0的湿润带。混合林较纯林稳定，草甸与高寒生态系统恢复力强且脆弱性低。农业生态系统因灌溉未纳入GLDAS-Noah模拟，其内在抵抗力可能被高估。
3.2 生态系统生产力抵抗力、恢复力与脆弱性的时空动态：Sen斜率与Mann-Kendall检验表明，长江中下游、四川盆地周边、云贵高原抵抗力显著下降（年率-0.0936、-0.1042、-0.2751），恢复力同步上升（0.0192、0.0289、0.0209）；脆弱性显著增加区集中于东南季风区，超95%变化区呈显著上升，多源GPP结论一致。
3.3 抵抗力与恢复力权衡加剧脆弱性：脆弱性上升区抵抗力年均降速0.074远超恢复力增速0.010。1982至2019年网格尺度抵抗力与恢复力总体呈显著权衡（P<0.05）；长江中下游2000年前权衡减弱（0.0047/年），2000后转协同。概率密度矩阵显示低抵抗（0至25分位）配任意恢复力格点萎缩，低恢复力（0至25分位）配任意抵抗力格点扩张，反映分布结构偏移。
3.4 生产力抵抗力与恢复力关键驱动：CatBoost在季风区恢复力模拟R达0.776、抵抗力0.650。SHAP显示温度对抵抗力与恢复力解释贡献均超50%，潜热SLHF等地表因子对抵抗力贡献超34%，陆气耦合因子对季风区恢复力贡献超10%。温度上升使抵抗力SHAP由正转负。地理探测器中，干旱强度（最低土壤湿度百分位）与持续时间交互对抵抗力解释q=0.40，与土壤类型交互对恢复力q=0.55。
讨论与结论翻译：讨论指出，高时间分辨率数据对解码骤发干旱下光合调控不可或缺——日尺度BEPS捕捉到0至20天内GPP负异常始发，而8天数据仅约10个观测点，会简化机制解释。季风区高脆弱性源于海气陆耦合高频触发与抵抗力-恢复力权衡向协同偏移，使植被陷入“低抵抗-低恢复”恶性循环。结论部分表明：该研究提出抵抗-恢复-脆弱多属性框架，耦合高时频植被数据与可解释机器学习，揭示1982至2019年季风区脆弱性年增0.019，高温异常为首要胁迫，潜热通量标志能-水胁迫转换；未来需拓展至全球并耦合GCM多情景模拟，以支撑生态系统保护与农业适应。论文发表于《Earth's Future》。
