复合闪旱–热浪极端事件下植被绿度响应的滞后特征

《Agriculture》：Delayed Vegetation Greenness Response to Compound Flash Drought–Heatwave Extremes

　　

复合闪旱–热浪极端事件（FDHW）使植被暴露于迅速发展的水分与热量胁迫之下，然而区域评估常将事件检测与植被响应混为一谈，且鲜有研究解析延迟的冠层轨迹。研究人员利用ERA5-Land土壤湿度和温度阈值，量化了1995–2024年生长季中国东北黑土区（NBSR）的FDHW事件；采用时空图神经网络（ST-GNN）对阈值提取的事件掩膜进行空间正则化，并将AVHRR NDVI（Advanced Very High Resolution Radiometer Normalized Difference Vegetation Index，高级甚高分辨率辐射计归一化植被指数）保留用于独立的事后影响评估。研究发现，闪旱和FDHW频率呈现强烈的年际变率，而非显著的单调趋势。FDHW发生频率从1995–2005年的3.8 d/生长季增加至2016–2024年的4.8 d/生长季，但Theil–Sen趋势接近于零（0.05 d/10年）。地表–大气合成分析表明，FDHW发生前土壤湿度逐渐耗竭，且在最高温度异常出现前约三天，地表能量释放由潜热（latent heat）向感热（sensible heat）转换。NDVI合成分析揭示了延迟的绿度响应：异常在事件后两周内维持负值，于参考FDHW格网点日后约一周达到最小值，随后在16–30天期间部分恢复。平均NDVI抑制幅度较小（短期?0.009；长期?0.006），但南部农田主导轨迹中有12.1%、北部林–农交错带轨迹中有10.7%维持持续的负异常。这些结果表明，NBSR地区FDHW的影响更少表现为事件频率的稳步上升，而更多地表现为延迟且空间异质的植被胁迫，凸显了设立事后监测窗口和跨传感器验证以支撑农业风险评估与适应规划的必要性。

## 一、研究背景与问题提出

东北黑土区（NBSR，约112–133° E，43–52° N，涵盖黑龙江、吉林、辽宁三省及内蒙古东部地区）作为全球四大黑土带之一，其肥沃的黑土层支撑着中国玉米和大豆生产的核心地带，对国家粮食安全具有战略意义。然而，该区域地处气候敏感的中纬度地带，升温速率超过全球平均水平，以雨养农业为主的耕地、局部灌溉农田及自然植被正日益暴露于极端天气事件，尤其是干旱和高温的复合威胁之下。

闪旱（flash drought）对该区域尤为相关，因其发展过于迅速，许多农业管理措施来不及响应。与传统干旱数月乃至数季的形成过程不同，闪旱以数天至数周内土壤水分亏缺的快速 intensification 为特征，常伴随高温和高水汽压 deficit（vapor pressure deficit, VPD）。快速 onset 使植被生理调节和管理适应的时间窗口极为有限；气孔关闭和光合作用下调可能在作物发育出更深根系或进入较不敏感的物候阶段之前就已发生。观测证据表明，随着气候变暖，闪旱频率在中纬度多地呈增加态势。闪旱常与热浪（heatwave）重叠，形成复合干旱–热浪事件（FDHW），其水热复合胁迫程度往往超过单一灾害之和。全球分析显示，复合干热极端事件已经 intensified，尤其在土壤 moisture 常处于能量限制与水分限制蒸散发过渡带的气候过渡区。

尽管叶片和冠层尺度上干旱与高温复合胁迫影响植被的生理路径已得到充分记录——土壤水分亏缺驱动气孔关闭，减少潜热通量并升高冠层温度；高温则直接损害Rubisco活性和光合电子传递——但这些机制的区域化表达仍难以准确度量。基于阈值的FDHW检测方法可将快速演变的事件分割为碎片化斑块，导致受影响面积和植被指数异常的估算存在噪声。更严重的是，当植被指数既用于识别干旱又用于度量其后果时，将出现循环论证问题：绿度下降可能成为事件定义的一部分，随后又被计为事件影响。此外，植被响应的时间滞后性增添了新的复杂性：若冠层绿度在气象 onset 数天后才响应，事件期间的NDVI可能错过主要响应，而后期测量则可能已反映恢复或额外胁迫。

卫星研究已记录到中国东北地区植被对闪旱和复合极端事件的响应存在空间异质性，技术综述表明光学植被指数仍是生态系统胁迫监测的核心工具，但需注意传感器、大气及指数特异性限制。近期工作在从历史和未来气候情景中检测复合干旱–热浪事件的深度学习框架方面取得进展，但农田与半自然植被在复合事件后恢复轨迹的定量比较仍然缺乏，事件 delineation 的空间一致性如何影响后续植被影响估算也尚待研究。这些不足使得跨土地覆盖类型和事件检测方法的植被响应比较变得困难。

基于此，本研究分析了1995–2024年生长季NBSR地区的FDHW事件。研究通过快速表层土壤湿度下降识别闪旱，通过格网点特异性温度百分位数识别热浪，再利用闪旱尺度的干热指数（Compound Dry-Hot Index）提取复合事件。阈值定义的掩膜经ST-GNN空间正则化处理后，进行NDVI响应分析，但NDVI本身被排除在事件定义之外。研究利用地表能量通量合成分析和逐像素NDVI异常，回答三个核心问题：（i）1995年以来FDHW频率和空间范围是否发生变化；（ii）事件 onset 前地表能量分配如何演化；（iii）NBSR南部农田主导区与北部林–农交错带在复合事件后的植被恢复有何差异。

## 二、主要技术方法

研究采用的关键技术方法包括：（1）基于ERA5-Land日值再分析数据（0.1°空间分辨率，1995–2024年5–9月）的物理阈值法识别闪旱和热浪，其中闪旱通过表层土壤湿度（SWVL1，0–7 cm）5天滑动窗口内从高于40th百分位数快速下降至低于20th百分位数判定，热浪通过日最高温度超过90th百分位数且持续≥3天判定，复合事件通过标准化干热 severity 指数（CDHF）≥1判定；（2）ST-GNN（Spatiotemporal Graph Neural Network，时空图神经网络）作为空间正则化工具，以距离阈值0.5°构建静态空间邻接矩阵，结合图卷积网络（GCN）提取空间特征、2层长短期记忆网络（LSTM，隐藏层大小64）捕获30天序列的时间动态，采用Focal Loss处理约1:31.9的极端类别不平衡，通过优化验证集F1-score确定决策阈值；（3）合成分析法（composite analysis），以FDHW affected格网点的 onset 日为Day 0，提取前后31天标准化异常序列，分析土壤湿度、潜热通量（SLHF）、感热通量（SSHF）和最高温度（MX2T）的演化；（4）基于NOAA AVHRR NDVI CDR Version 5产品（0.05°分辨率，重采样至0.1°）的独立事后植被影响评估，构建逐像素30年气候学基线，定义短期（1–15天）和长期（16–30天）两个事后评估窗口，计算NDVI异常（ΔNDVI = NDVI_affected ? NDVI_baseline），并区分南部农田主导区（48° N以南）和北部林–农交错带（48° N以北）进行恢复曲线对比，通过非恢复概率（P?(t)）及其半衰期和渐近非恢复分数（P?_∞）量化恢复特征；（5）水热相位空间分解，以土壤湿度异常和最高温度异常（均为Z-score）构建二维应力空间，分析不同胁迫组合下的延迟绿度响应。

## 三、研究结果

### 3.1 复合事件的发生趋势、空间格局与地表–大气耦合

研究结果表明，1995–2024年间三种极端事件类型呈现截然不同的变化轨迹。区域平均闪旱发生频率在1995–2005年为2.5 d/生长季，2016–2024年为2.9 d/生长季，Theil–Sen斜率为0.07 d/10年（95%置信区间?0.13至0.28；Mann–Kendall P=0.48）。FDHW频率同期分别为3.8和4.8 d/生长季，斜率为0.05 d/10年（95% CI ?0.21至0.38；P=0.59）。独立热浪频率斜率约为0.03 d/10年（P=0.82）。均无可检测的显著单调趋势，年际变率占主导。季节分布上，FDHW日数贯穿整个5–9月，而非集中于仲夏，且存在轻微前移趋势：7–8月占比从约48%（1995–2005）降至43%（2016–2024），而5月占比从约18%升至25%。

像素级趋势分析揭示三个地理上 distinct 的区域。西南半干旱–半湿润过渡区（43–46° N，115–120° E）呈现最强的正趋势，FDHW频率以0.5–1.5 d/10年增加，约97%的像素呈增加态势，但仅约四分之一达到P<0.05显著性。该区域年降水量400–500 mm，经历了广泛的草地转农田过程，土壤 moisture 常处于水分限制蒸散发的临界阈值附近。相反，黑龙江和吉林中部部分地区显示?0.2至?1.0 d/10年的下降趋势，与东亚夏季风降水增强和灌溉设施发展空间重合。北部大兴安岭和长白山山前地带趋势微弱（?0.05至0.05 d/10年），高森林覆盖和充足基线降水缓冲了变化。

合成分析显示了一致的标准化异常演化序列。Day ?15至Day ?10期间，土壤湿度异常以约?0.03σ/day逐渐下降，潜热通量异常维持正值而感热通量异常接近零，温度异常无系统偏离。Day ?5至Day ?3为关键转折期：土壤湿度耗竭加速至?0.1σ/day，潜热和感热通量异常曲线交叉，潜热转负、感热转正。该交叉点约在Day ?3，即温度峰值前约三天。此后，最高温度异常以约0.25σ/day快速 intensified，于Day 0达到+2.1σ。Day 0至Day +5，四变量维持极端异常值，反映干土壤、减少的潜热冷却、增强的感热加热和高温之间的正反馈。Day +5后约10–15天逐渐恢复气候学状态。

### 3.2 ST-GNN空间正则化性能

在2016–2024年独立测试集上，完整ST-GNN的精确率（precision）为0.31，召回率（recall）为0.62，F1-score为0.42（优化阈值0.35）。这些指标反映了极端类别不平衡下稀有事件分类的固有难度。ST-GNN的核心贡献在于空间正则化而非逐像素分类提升：相比原始物理阈值输出的碎片化掩膜，经图邻域聚合后空间连续性显著改善，主要事件 footprint 得以保留，孤立检测被抑制。

消融实验阐明了时空组件的功能分工。去除LSTM的GNN-only变体导致判别能力近乎完全丧失（F1=0.04），确认FDHW识别本质上是时间累积问题；去除GCN的LSTM-only变体获得略高的逐像素F1（0.43），但空间碎片化更严重。完整ST-GNN以约1个F1点的代价换取了更空间连贯的事件掩膜。主要植被响应结论（10–15天延迟抑制、南北带恢复对比方向）在原始阈值、LSTM-only和ST-GNN三种事件定义下均成立，表明ST-GNN并非检测这些生态模式的必要条件，而是为合成分析和影响评估提供更清晰的空间边界。

### 3.3 植被绿度对复合事件的响应

30年记录中，物理阈值标准识别了847,000个FDHW affected格网点日。短期窗口（1–15天）和长期窗口（16–30天）均显示相对于多 年NDVI基线的系统性负异常偏移。短期平均NDVI_Anom为?0.009（95% CI ?0.010至?0.007），约?0.13σ；长期平均为?0.006（95% CI ?0.008至?0.004），约?0.09σ。逐日追踪显示，事件后两周内平均NDVI异常维持负值，约第7天达最深（约?0.012），第15天后恢复。

以2007年为典型重灾年 illustrate 高时间分辨率响应：NDVI异常在事件后约5–8天变为可检测负值，10–15天抑制最深（平均最低NDVI_Anom≈?0.025，约?0.35σ），约20天内恢复至近气候学水平。按物候阶段分层， peak 生长季（7–8月）事件的滞后NDVI响应最为显著，早期（5–6月）响应较弱且变异更大。就绝对幅度而言，平均效应仅占区域典型生长季NDVI范围的约2–3%，主要生态信号是延迟响应和异常分布的展宽，而非整体性植被崩溃。异常分布随时间展宽：短期标准差（0.049）小于长期（0.055），Levene's test P<0.001，反映恢复与不恢复格网点的生态分化。

分带恢复曲线显示，南部农田主导区的渐近非恢复比例略高于北部林–农交错带：P?_∞为12.1%（95% CI 11.8–12.4%）vs. 10.7%（95% CI 10.4–11.0%），差异约1.4个百分点。从事件终止起算的实证恢复半衰期在两区均较短（南部约1.6天，北部约1.4天）：约三分之一 depressed 格网点一天内恢复至基线，约一半两天内恢复，约84%一周内恢复，残余11–12%在45天视界内维持 depressed。该 zonal 差异幅度小且对事件定义敏感，作为方向性一致的 pattern 而非精确估计呈现。

水热相位空间分解显示：wet/cool象限平均NDVI_Anom为?0.003（n≈1.2×10⁵），作为气候学基线；dry/cool象限为+0.002（n≈4.1×10⁵），可能反映适度土壤干燥下减少的非生产性营养生长或深层土壤 moisture 的持续供应；wet/hot象限最负（?0.011，n≈3.0×10⁵），与 moist 条件下高温直接损害光合功能、且高基线植被 growth 导致较大偏离预期的现象一致；dry/hot象限为?0.006（n≈2.3×10⁵，含ST-GNN边界计数）。dry/hot象限较wet/hot不更负， likely 反映了基础NDVI和事件群体的系统性差异：复合事件发生于基线NDVI已部分抑制的 dry 条件，而最 severe 的wet/hot事件发生于有利于极端温度 buildup 和高基线植被 growth 的 moist 条件。四象限均值阶乘对比（wet/cool + dry/hot ? dry/cool ? wet/hot）约?0.012 NDVI单位，但由于各象限均值受基线NDVI和事件群体差异主导，该相位空间仅作为描述性分析，不支持多plicative synergy 的结论。

## 四、讨论与结论

讨论部分围绕三个主题展开。第一，关于闪旱与热浪在复合事件中的相对作用：FDHW频率的年际波动更多追踪快速土壤湿度 depletion 的变化而非热极端的变化，对早期预警的启示在于土壤 moisture 监测、灌溉 timing 和水分保持措施应与热浪防备同等重视。空间异质性趋势凸显了区域平均统计不足以指导空间 target 适应，西南过渡区作为地表–大气耦合热点应优先纳入干旱监测和水资源规划。

第二，关于ST-GNN在复合事件 delineation 中的角色：其核心贡献是空间正则化而非独立事件检测，应被理解为有学习的空间后处理步骤。关键方法论贡献在于将事件识别与影响评估基本分离，减少了传统干旱–植被研究中循环论证的风险。

第三，关于植被监测与生态评估启示：约一至两周的NDVI响应延迟是最直接可用的结果之一，意味着当前在极端事件期间或 immediately after 评估植被状况的监测系统可能 substantially 低估绿度下降。事件后至少两周的持续监测窗口、以及基于分布尾部的指标而非仅区域平均NDVI，对农业胁迫 screening 更为有效。南部农田区略高的 persistent 异常率可能与单一作物系统的功能冗余较低、根系 profile 较浅、以及收获管理干扰等因素有关，但48° N边界也分离了气候、物候和管理强度差异，故该对比不 isolate 土地覆盖效应。NDVI仅为冠层绿度信号而非直接生理损害测量，与SIF、GPP、LAI及县域产量统计的验证是未来优先方向。

研究结论部分指出：利用30个生长季的ERA5-Land和AVHRR NDVI数据，本研究考察了中国东北黑土区的闪旱–热浪复合事件及其后续冠层绿度响应。区域平均闪旱和复合FDHW频率在1995–2024年间无统计显著趋势，分别围绕2.5–2.9和3.8–4.8 d/生长季波动，年际变率显著。空间上，西南半干旱过渡区呈现最正局部趋势，中部平原部分下降，与区域水分 availability、灌溉和土地利用历史差异一致。FDHW affected格网点中，土壤湿度在事件前下降，标准化地表通量异常在峰值温度前约三天转向感热主导。NDVI在气象 onset 后不立即下降，最强负异常出现在约两周内（最深约一周），事件结束后中度负异常仍更频繁。南部农田主导区的 persistent 负NDVI异常率略高于北部林–农交错带，但该对比幅度小且对所用事件定义敏感。鉴于最大NDVI响应存在延迟，监测应持续至FDHW onset 后至少两周。不使用NDVI定义FDHW也避免了循环归因问题。这些NDVI异常是否对应产量下降或田间作物胁迫尚待解决，未来需借助SIF、GPP、LAI、产量统计和管理记录等独立证据加以验证。