
Increased spread of global flash droughts threatens vegetation productivity resilience来自 Increased spread of global flash droughts threatens vegetation productivity resilience | Nature Communications 基本信息题目全球爆发性干旱扩展加剧威胁植被生产力恢复力Increased spread of global flash droughts threatens vegetation productivity resilience期刊名称Nature Communications发表时间2026年3月17日Volume 17, Article 4050作者信息Renjie Guo, Xiuchen Wu, Pei Wang, Tiexi Chen, Xin Chen, Jiangtao Cai, Xiaona Wang, Zifan Zhang, Zekai Meng, Yiran LiuDOI10.1038/s41467-026-70417-z核心观点1950至2023年间全球爆发性干旱骤旱发生频次、严重程度与爆发速度显著上升且急剧从传统的湿润热点区向非热点区干旱/半干旱区扩散。植被总初级生产力GPP对骤旱的恢复力普遍低于传统缓旱在热点区的恢复力损伤尤为严重骤旱热点区恢复力相比非热点区下降 Δ -0.07 ± 0.02t -3.25p 0.001。二氧化碳CO2施肥效应虽然是维持恢复力的最主要正面因子但在骤旱引发的突发性高温和高蒸汽压差VPD胁迫下其缓解作用被明显削弱揭示出全球生态系统在骤旱加剧扩展趋势下的高脆弱性。️ 一、 引言部分简要总结 (Introduction Summary)现实与科学背景气候变暖与人类活动加剧了爆发性干旱Flash drought骤旱的发生。与演变缓慢的传统干旱Slow drought缓旱不同骤旱具有爆发速度极快、土壤水分在亚季节尺度内急剧耗竭的特征对全球农业生产和陆地碳汇造成了毁灭性打击。前人研究局限与瓶颈既有研究大多集中于土壤水分高、骤旱频繁发生的已知“热点区域”Hotspots多为湿润区对于骤旱是否正在向“非热点区域”Non-hotspots干旱与半干旱区扩展及其生态风险尚不明确。此外植被生产力恢复力Resilience在应对骤旱与缓旱时的响应分化特征以及驱动这种响应差异的底层生物物理机制仍未得到系统定量。核心科学问题本研究旨在解决三个核心问题阐明1950至2023年全球骤旱的时空演变与全球扩展特征定量评估植被生产力对骤旱与缓旱恢复力的响应分化厘清气象水文因子、CO2施肥效应、植被性状与土壤属性等对骤旱恢复力的控制路径。 二、 研究方法详细介绍与分析 (Methodology Framework)技术路线与数据整合时空范围与尺度覆盖全球陆地范围空间分辨率 1°时间尺度为 5 日/候 pentad。干旱事件识别跨越 1950–2023 年植被生产力响应分析集中于 2001–2019 年。核心数据源整合 ERA5-Land 与 GLDAS_CLSM 每日土壤水分数据计算 1 米深度加权平均土壤水百分位数。植被生产力采用 FluxSat 每日 GPP 数据与 CSIF 太阳诱导荧光SIF数据并使用 FLUXNET2015 站点通量数据进行交叉验证。气象因子气温、辐射、降水、蒸发、VPD采自 ERA5-Land植物性状与立地因子最大根深、冠层高度、树木密度、土壤阳离子交换量采自全球高分辨率生态图谱。未来预测采用 CMIP6 中 8 个地球系统模型在 SSP245 情景下的模拟数据。骤旱与缓旱判定标准骤旱Flash drought候平均土壤水分百分位数从 40% 以上降至 20% 以下且爆发阶段每候下降速率不低于 5% 百分位数持续时间至少 4 候。缓旱Slow drought土壤水分同样从 40% 以上降至 20% 以下但爆发阶段至少有一候下降速率小于 5% 百分位数持续时间至少 4 候。多维骤旱指标MFDI与热点区划分整合骤旱发生比例、频次、严重程度和爆发速度构建标准化多维骤旱综合指标MFDI均值 100标准差 10。将全球划分为 29 个独特的生态-气候-植被综合区MFDI 区域均值大于 100 的划定为骤旱热点区热点区中 99.46% 为湿润区其余为非热点区。恢复力定量与归因模型工作流恢复力定义将干旱发生后 2 年内 GPP 距平均值与多年生长季均值之差定义为恢复力Resilience Ya - Ym恢复力低于 0 表示干旱阻碍了恢复。归因模型构建包含 15 个候选因子的随机森林模型300 棵决策树5 个叶节点。使用方差膨胀因子VIF 5剔除共线性采用袋外数据置换检验Out-of-Bag Permutation Importance计算因子相对重要性利用偏依赖图Partial Dependence Plots, PDP定量分析恢复力对各主导因子的非线性响应。采用 BEAST 模型检测时间序列突变拐点。 三、 主要研究结果与充实数据分析 (Main Results Quantified Data)全球骤旱的时空演变与扩展趋势时间加速1950–2023年间全球骤旱频次距平显著上升斜率 6.63 ± 0.66 次/年z 5.01p 0.001严重程度距平显著加深斜率 0.01% ± 0.01%/年z 4.80p 0.001爆发速度距平显著加快斜率 0.02% ± 0.01%/候/年z 5.66p 0.001。突变检测显示拐点出现在 21 世纪初2002/2003年。骤旱占总干旱事件的比例呈显著加速趋势斜率 0.14% ± 0.09%/年z 6.16p 0.001证实缓旱向骤旱的转折急剧加剧。空间扩张骤旱热点区 MFDI 均值102.29 ± 10.13显著高于非热点区92.94 ± 15.57t 38.36p 0.001。然而非热点区的 MFDI 上升趋势显著高于热点区Kendall τ 趋势差异 t -7.07p 0.001表明骤旱正在向干旱及半干旱的非热点区迅速扩展。气象异常对比骤旱爆发阶段伴随着比缓旱更剧烈的降水负距平Δ -0.296 ± 0.004 mm/dt -73.58p 0.001以及更显著的气温正距平Δ 0.52 ± 0.01 °C/dt 50.97p 0.001、太阳辐射正距平Δ 5.59 ± 0.07 J/m²/dt 82.11p 0.001与 VPD 正距平Δ 0.63 ± 0.01 hPa/dt 44.51p 0.001。植被生产力恢复力的响应分化全球整体差异全球 GPP 对骤旱的恢复力显著低于缓旱恢复力差值 Δ -0.012 ± 0.004t -2.74p 0.01。热点区与非热点区对比在骤旱热点区GPP 对骤旱的恢复力高于缓旱Δ 0.04 ± 0.02t 2.45p 0.01而在非热点区GPP 对骤旱的恢复力显著低于缓旱Δ -0.05 ± 0.02t -2.38p 0.01。区域间恢复力衰退热点区植被对骤旱的恢复力显著低于非热点区Δ -0.07 ± 0.02t -3.25p 0.001热点区对缓旱的恢复力同样显著低于非热点区Δ -0.15 ± 0.02t -8.80p 0.001。底层驱动机制与控制路径CO2 施肥效应的脱节CO2 施肥效应β是骤旱和缓旱恢复力的最重要单一因子。然而偏依赖分析表明在骤旱条件下恢复力对 CO2 施肥效应增加的敏感性明显低于缓旱。这是由于骤旱引发的极端高温和高蒸汽压差限制了植被 CO2 吸收损害了光合酶活性并破坏了叶绿体结构削弱了 CO2 施肥效应的正面缓解作用。气象胁迫主导气象因子对非热点区恢复力的贡献最高骤旱贡献率 39.07% ± 1.24%缓旱贡献率 34.72% ± 0.99%。骤旱引发的突发性蒸发剧增和水分亏缺使恢复力呈陡降趋势。前置生物量的双刃剑效应干旱发生前及干旱期间的植被生长状况前置 GPP对骤旱恢复力具有更强的控制力。热点区较高的前置生物量虽能产生短期的正向遗留效应但在骤旱加剧时会加速土壤水分的蒸腾耗竭反过来导致恢复力更快衰退。上述骤旱发生频次、严重程度与爆发速度的时空演变详见 Figure 1骤旱热点区与非热点区的空间分布与扩展趋势详见 Figure 2全球及区域尺度下 GPP 对骤旱与缓旱恢复力的差异详见 Figure 3骤旱热点区与非热点区恢复力驱动因子的相对重要性及偏依赖响应分别详见 Figure 4 与 Figure 5骤旱扩散威胁植被恢复力的综合机制框架详见 Figure 6。图片名Figure 1图例说明Figure 1 | 全球爆发性干旱骤旱的时空演变格局。a-c 分别显示骤旱频次、严重程度和爆发速度的空间分布d-f 显示三者在年均蒸发量MAE与年均降水量MAP空间中的分布g-i 显示 1950–2023 年间骤旱频次、严重程度和爆发速度距平的年际变化虚线代表突变拐点BEAST检测实线为线性回归趋势Mann-Kendall检验p 0.001。图片名Figure 2图例说明Figure 2 | 骤旱热点区与非热点区的空间特征及扩展趋势。a 骤旱热点区金黄色与非热点区灰色的空间分布b 热点区与非热点区多维骤旱指标MFDI的箱线图差异t 38.36p 0.001c 1950–2023 年间热点区与非热点区 MFDI 变化趋势Kendall τ的概率密度分布t -7.07p 0.001。图片名Figure 3图例说明Figure 3 | 植被总初级生产力GPP对骤旱与缓旱恢复力的差异。a 骤旱与缓旱 GPP 恢复力差值Δ Resilience骤旱减去缓旱的空间分布b 骤旱热点区与非热点区内 GPP 对骤旱浅红和缓旱浅蓝恢复力的均值与分布对比字母不同代表存在显著差异p 0.05。图片名Figure 4图例说明Figure 4 | 骤旱热点区 GPP 恢复力驱动因子的相对重要性与偏依赖关系。a 随机森林模型导出的各类因子重要性排名红点为骤旱蓝点为缓旱b-k 为恢复力对前 10 个主导因子包含 CO2 施肥效应 β、干旱期间 GPP、最大根深、土壤水分、前置 GPP、气温、VPD、辐射、降水、蒸发的偏依赖响应曲线阴影为 95% 置信区间。图片名Figure 5图例说明Figure 5 | 骤旱非热点区 GPP 恢复力驱动因子的相对重要性与偏依赖关系。a 随机森林模型导出的各类因子重要性排名b-k 显示非热点区恢复力对前 10 个主导因子的非线性偏依赖响应曲线红线为骤旱蓝线为缓旱阴影为 95% 置信区间。图片名Figure 6图例说明Figure 6 | 骤旱全球扩展及其威胁植被生产力恢复力的概念机制图。展示了缓旱向骤旱的加速转变、骤旱向非热点区的快速扩展、气象胁迫的加剧、CO2 施肥效应敏感性的降低以及植被生产力恢复力的整体衰退过程。 四、 讨论部分的不足与未来展望 (Discussion Limitations Perspectives)研究不足与不确定性来源空间分辨率与源数据局限土壤水分融合自 ERA5-Land 与 GLDAS_CLSM 数据集受限于全球长时序列1950–2023年每日数据的处理瓶颈采用了 1° 粗网格可能平滑掉局部复杂地形与微气候对骤旱的调制作用。时间序列跨度不匹配高质量全球每日 GPP 观测数据FluxSat2001–2019年跨度 19 年相比于 74 年的气象数据较短虽然足以捕捉近期恢复力的变异但对超长期恢复力演变趋势的约束仍受限。通量站点空间分布不均用于实测验证的 FLUXNET 站点中具备完整骤旱和缓旱恢复序列的站点数量较少且全球分布不均限制了验证的代表性。地球系统模型CMIP6预测偏保守SSP245 情景下 CMIP6 模型预测的未来骤旱增长速率低于实际土壤水分观测值主因是现有模型低估了土壤水分对气温和降水变化的敏感性。静态生态性状假设归因模型中使用的最大根深、冠层高度等植物性状和土壤属性为空间静态图谱未纳入气候变暖背景下的长期动态生理与演替过程。未来研究展望高精细多源数据同化结合高空间分辨率遥感与深层土壤水分同化技术精细化识别局地尺度上的骤旱爆发过程提升骤旱早期的生态预警能力。地球系统模型改进在陆地表层模型中引入更精准的植物水力学性状如气孔导度对高温和高 VPD 的动态响应纠正对骤旱发生率及植被恢复力损伤的低估。野外操控实验与生理机制解构开展多因子增温、干旱、CO2浓度升高野外操控实验深入解构骤旱引发的水力衰竭、光合结构破坏与碳饥饿对生态系统恢复力损伤的生物物理路径。