Critical slowing down of semiarid vegetation resilience is amplified by intensifying heatwaves来自 Critical slowing down of semiarid vegetation resilience is amplified by intensifying heatwaves | Nature Communications基本信息中文题目热浪加剧放大全球半干旱区植被弹性临界减速效应英文题目Critical slowing down of semiarid vegetation resilience is amplified by intensifying heatwaves期刊名称Nature Communications发表时间2026年7月3日在线发表作者信息付莉、黄建平、李昶豫、魏文敏、彭明、张国龙、任裕、邓建明兰州大学草地农业生态系统国家重点实验室、生态学院、大气科学学院DOI10.1038/s41467-026-75130-5核心观点本研究聚焦于全球半干旱区植被生态系统利用多源卫星遥感指标定量评估了2001至2024年间植被弹性的动态变化趋势。研究表明尽管在二氧化碳富集等因素驱动下全球半干旱区呈现广泛的“变绿”植被绿度增加现象但超过半数的区域实际上正在经历临界减速Critical Slowing Down, CSD过程表现为系统一阶时间自相关系数lag-1 TAC持续显著上升意味着系统受扰动后的恢复速度放缓、生态弹性实质性下降。空间分析和数学模型进一步揭示近年来不断加剧的热浪尤其是最大热浪强度与气候干旱化发生复合作用显著放大了这种弹性衰退。这一发现确证了植被“变绿”掩盖了其内部稳定性受损的科学事实强调必须将弹性诊断指标纳入常规生态监测体系中。️一、 引言部分简要总结(Introduction Summary)研究背景与科学问题基于遥感观测的大量文献表明过去数十年内全球陆地表面呈现出广泛的“变绿”Greening趋势半干旱区亦不例外。然而仅依赖平均绿度或生物量指标的单调上升无法客观准确地映射生态系统内部的真实稳定性特征。随着全球气候变暖极端气候事件的频发与加剧如高温热浪、极端干旱等正对脆弱的半干旱区生态系统施加急性热应力和水分胁迫。热浪常伴随高饱和差VPD会强烈抑制植物的光合作用与结构生长。在此现实下探明半干旱区植被系统能否在扰动后有效恢复即维持生态弹性以及极端热浪事件如何在背景气候变异下调制这种弹性成为了全球变化生态学亟待解决的核心科学问题。前人研究局限性与核心瓶颈既往关于植被弹性的研究多侧重于探讨生态系统对平均气候条件或常态化温湿变率的响应且高关注度往往集中在亚马逊雨林等热带湿润复杂生态系统。半干旱区作为全球陆地干燥度扩张的主要贡献区域其生态脆弱性极高但针对其生态恢复力、内部稳态演变规律的研究较为匮乏。更为关键的是传统评估手段容易将“结构性相变”如灌木入侵导致的绿度指标虚高等同于生态系统的健康好转从而忽视了绿度趋势与系统内部弹性之间可能存在的脱钩Decoupling现象。热浪特征频率、历时、强度在空间上如何协同背景水文气候因子驱动系统稳态失衡其物理与生物化学响应路径至今仍不清晰。二、 研究方法详细介绍与分析(Methodology Framework)空间范围界定与样本过滤研究区域依据干燥度指数Aridity Index, AI P/PET进行严格界定。采用1961—1990年气候基准期的 0.5度 空间分辨率数据锁定 AI 值处于 0.20 到 0.50 之间的典型半干旱区。为剔除高纬度冰雪覆盖导致的遥感反演噪声将空间纬度截取在 60°S 至 60°N 之间。在样本清洗阶段使用 MODIS MCD12C1 土地覆盖产品仅保留 2001 至 2024 年间土地覆盖类型IGBP 1-10 类包含森林、灌丛、稀树草原及草地始终保持绝对未变的原生自然植被像元。同时基于 GFEDv5 烧失面积数据集计算 2002 至 2022 年间的累积火烧比例剔除超出 95 百分位的高频火烧干扰像元从而确保时间序列分析免受非气候性剧烈干扰如毁林转垦、火灾重创等的混淆。核心数据源与环境因子提取为了全方位捕捉植被的生理与结构特征本研究综合提取了三种互补的卫星监测指标1生理活性指标采用 GOSIF v2 月尺度 0.05度 指标基于 OCO-2、MODIS 及气象重分析数据训练而成直接表征光合生理胁迫2结构绿度指标采用 MODIS MOD13C2Collection 6.1月尺度 NDVI 数据剔除多年平均 NDVI 0.1 的荒漠像元3冠层结构指标采用 MODIS MOD15A2H 8天 500米 LAI 产品并利用双线性插值将其合成为月尺度 0.05度 数据。在气象因子方面利用 ERA5 逐小时 2米 气温数据反演日最高气温Tmax用于识别热浪。降水数据则提取自月尺度 TerraClimate 降水数据集。核心数学模型与算法工作流信号剥离去趋势STL分解为准确捕获临界减速信号对月尺度植被指数时间序列进行基于Loess局部加权回归的季节趋势分解STL算法将原始序列精确切分为长期趋势项、季节项和残差项Anomaly。随后将残差项用于弹性指数计算彻底摒弃了由于均值漂移和年际物候波动对时间自相关系数带来的虚假人为抬升。临界减速指标测度在去趋势后的残差序列上应用 60个月5年的右对齐滑动窗口通过一阶自回归模型 $X_{t1} \phi X_{t} \epsilon_{t}$ 计算一阶时间自相关系数lag-1 TAC, 即模型中的参数 $\phi$。$\phi$ 值的增大即科学定义上的“临界减速”反映了恢复速度放缓。同时利用 36个月 和 84个月 窗口作为稳健性检验。弹性变化趋势则通过对时间序列项 $\phi$ 进行普通最小二乘线性回归的斜率Slope进行度量。极端热浪事件识别热浪采用相对阈值法定义即当某像元的日最高气温连续 3天 或以上超过该日历天结合前后7天共15天滑动窗口在 2001—2024 年基准期历史分布的 90 百分位数即记为一次热浪。据此统计年尺度的热浪频次Count、平均单次历时Duration、累积超额强度Cumulative Intensity及最大超额强度Maximum Intensity。空间重分析与多元因果归因空间误差模型Spatial Error Model, SEM针对传统 OLS 回归中残差存在强空间自相关的问题Morans I 检验显著利用残差变差函数Variogram拟合出 93公里 空间自相关空间范围。以此构建 row-standardized 空间权重矩阵 W引入自回归误差项模型将 SIF lag-1 TAC 的长期趋势作为响应变量不仅将热浪各项指标与以降水趋势表征的背景气候变率作为解释变量更创新地引入了“最大热浪强度 $\times$ 降水趋势”的交互效应项。偏最小二乘结构方程模型PLS-SEM利用 R 语言 plspm 包在各变绿/弹性象限内部构建路径网络。将最大热浪强度趋势、降雨趋势、SIF平均趋势以及SIF lag-1 TAC趋势纳入统一体系通过1000次 自助抽样法Bootstrap估计路径系数用以剖析限定范围内的多变量条件协变网络。三;主要研究结果与充实数据分析(Main Results Quantified Data)半干旱区植被变绿与弹性衰退的广泛脱钩数据分析表明在 2001 至 2024 年间全球半干旱区呈现出显著的绿度上升与弹性下降并存态势。基于日光诱导叶绿素荧光SIF的测算高达78.4%的半干旱植被区表现为光合活性增长SIF 呈正趋势趋势然而与此同时有56.6%的陆地像元其 SIF lag-1 TAC 呈现出统计学上显著的上升趋势空间平均增长斜率 trend 0.0016 /年确证了临界减速现象的广泛存在。这一结论在以 NDVI 和 LAI 测算的模型中获得了高度一致的稳健性交叉验证。图片名Figure 1图例说明Figure 1 | 2001—2024年全球半干旱区植被弹性与绿度空间趋势图。图1a 显示全球半干旱区 56.6% 的区域 SIF lag-1 TAC 显著上升揭示系统处于临界减速状态。图1c 显示全球空间平均 lag-1 TAC 的时间演变其长期线性趋势持续走高。图1e 与 1f 定量给出了不同植被群落的响应非对称性落叶阔叶林、稀树草原以及密灌丛的临界减速幅度最为剧烈而草地系统相对平缓同时变绿Mean SIF 上升与弹性衰退Mean lag-1 TAC 上升发生了大范围脱钩。双维象限框架下的演变特征与群落非对称性利用 SIF 趋势与 SIF lag-1 TAC 趋势的联合分布构建二维双维象限。空间统计表明47.5%的研究区域被划归为“象限1绿度增加但弹性衰退”Greening and declining resilience该比例大幅超越了表现为良好稳态的“象限4绿度增加且弹性增强”区域占33.8%。区域异质性上南美半干旱区SAM和南亚半干旱区SAS的临界减速趋势最为严重详见 Figure 1d。群落非对称性表现为象限1弹性受损区中木本成分如开阔灌丛、稀树草原、阔叶灌丛的占比明显高于象限4而象限4则由纯草地系统绝对主导草地占比达91.4%相较于象限1的81.1%见 Figure 3c。这表明木本植物入侵灌丛化虽然推高了卫星遥感捕获的绝对绿度但其自身高昂的蒸腾需求在水文限制下反而削弱了系统的恢复力。图片名Figure 2图例说明Figure 2 | 全球半干旱植被绿度—弹性双维变化象限及条件概率分布图。图2b 核密度分析确证象限1变绿兼弹性下降占研究区的绝对多数47.5%。图2c 条件概率差值 $\Delta$ 森林图表明在澳大利亚AUS、北非NAF、欧洲EU和北美NAM等地区负值结果意味着弹性衰退在“变绿”背景下发生的概率显著高于“变黄”背景。极端热浪对弹性减速的放大机制与背景调控时序分析显示2001—2024年全球半干旱区年均热浪发生频次、持续历时以及最大、累积强度均呈现明显的非线性上升并在2011年前后达到阶段性峰值详见 Figure 3a同时年降水量整体步入震荡下降通道呈现显著的气候干旱化Drying态势。空间误差模型SEM结果指出在所有热浪因子中最大热浪强度的跃升与植被临界减速lag-1 TAC 上升展现出最为稳健和强烈的正相关关系。更为关键的是背景水文条件对该效应具备显著的调制路径详见 Figure 4c 的交互效应在降水减少的严重干旱化背景下-1 SD 降水趋势最大热浪强度对系统临界减速的边际正效应显著增强回归斜率陡峭上升而在背景水文条件充沛1 SD 降水趋势即 Wetting 趋势的区域该项正相关关联则大幅削弱甚至趋于水平平缓。图片名Figure 3图片名Figure 4图例说明Figure 4 | 全球半干旱区最大热浪强度与降水趋势的空间分布及其在空间误差模型SEM中的交互效应。图4a2 与 4b2 分别定量给出了标准化后的最大热浪强度趋势显著正相关与降水趋势显著负相关对临界减速斜率的边际效应影响。图4c 交互效应曲线上清晰展示当背景水文趋于变干-1 SD红线时最大热浪强度驱动 lag-1 TAC 上升弹性减速的斜率最强揭示了复合干热极端事件对生态稳态崩溃的叠加放大机制。内部象限多元路径网络解析在各象限内部实施的偏最小二乘结构方程模型PLS-SEM进一步展示了条件协变网络。在象限4良性稳态变绿中降水趋势与 lag-1 TAC 呈显著负相关路径系数 -0.060而最大热浪强度与 lag-1 TAC 呈显著正相关路径系数 0.042支持了水分促进弹性、高温侵蚀弹性的普适物理路径。而在已经发生脱钩的象限1中最大热浪强度不仅表现出与生理绿度SIF的正向共变路径系数 0.083同时与弹性指标的直接关联减弱这有力地表明极端热应力在驱动植物群落结构调整如促使部分耐热木本成分扩张从而维持表面绿度的同时已悄然使系统内部结构极化增加了潜在的生态稳态转换风险。图片名Figure 5图例说明Figure 5 | 四个绿度—弹性象限内部的偏最小二乘结构方程模型PLS-SEM路径结构。图展示了在剔除空间共线性与多因子干扰后降水趋势、最大热浪强度趋势与植被绿度趋势、弹性退化趋势之间的象限内部条件路径系数标有星号 * 表示95%置信区间不含零具显著性定量直观地刻画了极端高温在不同演化稳态下的复杂控制规律。四、 讨论部分的不足与未来展望(Discussion Limitations Perspectives)研究不足与不确定性来源弹性代理指标的间接性与局限基于 CSD 理论框架提取的一阶时间自相关系数lag-1 TAC属于生态系统弹性变化的间接动力学诊断指标。该指标易受到非平稳随机过程、长期趋势残留项以及背景噪声结构变异的干扰。虽然本研究通过多种窗口期和不同遥感源进行了交叉稳健性清洗但仍无法完全等同于由控制火试或大尺度天然破坏后获取的绝对物理恢复力。背景水文因子的精细度不足研究中仅采用年降水量趋势作为背景水文气候变干Drying的单一代理变量在物理机制上略显简化。单一降水指标未能全面表征土壤表层/深层水分动态、地下水埋深波动亦无法完全涵盖陆气相互作用中饱和差VPD引起的强蒸腾损耗且未能充分考量极端事件发生时存在的气候—植被间驱动响应时滞效应。群落结构解译的空间分辨率限制利用 0.05度约5公里分辨率的混合像元进行长序列分析无法精确在个体或微观群落尺度上彻底区分草本与木本植物各自独立的生理贡献率。因此对由于“灌丛入侵”导致绿度与弹性脱钩的生物学机制推论仍有赖于未来更高分辨率的群落物种谱系制图来提供直接判据。未来研究重点与展望方向多源干旱度诊断指标并轨未来的深化研究应当积极引入多维度的综合干旱与水分胁迫指标如标准化降水蒸散指数SPEI、标准化降水指数SPI以及卫星微波反演得到的浅表层土壤湿度数据。同时需要定量计算极端事件对植被稳态影响的非线性时滞窗口以期更精准地刻画生态系统底层的水分胁迫归因路径。多尺度野外控制实验重塑亟待在全球不同典型半干旱区野外长期定位站开展人工模拟热浪强度的增温控制实验配合微观层面的叶片热耐受性叶片热安全裕度和群落内部水分利用效率测量在微观机理上确证热浪强度与背景水分亏缺对植被稳态转换的联合控制阈值。弹性预警系统与Routine生态监测重构鉴于植被绿度增加往往对内部生态脆弱性具有强烈的掩盖效应文章最后前瞻性地指出未来的生态环境常规监测体系不能再盲目依赖传统绿度如单调 NDVI 趋势作为单一健康指针。必须开发并轨全时段、常规化的植被稳态弹性损失在线诊断与早期预警系统Early-Warning Systems以此科学指导 Compound 极端干热背景下全球及区域性精准生态保护策略的制定。