尼罗河流域增强的复合干旱–热浪（CDHW）事件削弱生态系统水分利用效率（WUE）：GPP与ET响应、干旱记忆及气候梯度调控

《Journal of Environmental Management》：Intensifying compound drought–heatwave events weaken ecosystem water-use efficiency in the Nile River Basin

　　

复合干旱–热浪（compound drought–heatwave, CDHW）事件在尼罗河流域（Nile River Basin, NRB）加剧。热浪领先型（heatwave-leading）事件较干旱领先型（drought-leading）事件增加更快。复合事件通常较单一极端事件引起更强的总初级生产力（gross primary productivity, GPP）和蒸散发（evapotranspiration, ET）降低。多月干旱记忆（multi-month drought memory）强烈调控CDHW期间生态系统水分利用效率（ecosystem water-use efficiency, WUE）。生态系统WUE的气候控制沿干旱度梯度（aridity gradients）发生转变。

论文解读：尼罗河流域复合干旱–热浪事件对生态系统水分利用效率的影响

研究背景与问题：气候变暖显著增加了极端水文气候事件的频率、强度和持续时间，干旱和热浪已成为威胁陆地生态系统的重要极端事件。越来越多证据表明，二者可同时或相继发生，形成复合干旱–热浪（CDHW）事件，其影响可超过单独干旱或热浪影响之和。CDHW常与持续性大气阻塞和环流异常有关，抑制降水并增强地表加热；土壤水分下降减少潜热通量、增加感热通量，强化陆–气耦合，从而放大近地表增温并延长干旱。事件序列也影响生态响应：热浪领先型事件加速土壤水分耗竭并可能触发干旱，干旱领先型事件通过抑制蒸发冷却强化热浪。陆地生物圈吸收约30%人为CO₂排放，但该碳汇对气候极端事件高度脆弱。生态系统水分利用效率（WUE）定义为总初级生产力（GPP）与蒸散发（ET）之比，连接生态系统碳吸收与水消耗，是环境胁迫下生态系统韧性的重要指标。干旱通常通过土壤水分限制和气孔关闭降低WUE；热浪提高大气蒸发需求和蒸汽压亏缺（VPD），常使水损失相对碳获取增加。二者共存可产生非线性响应，显著改变生态系统碳水耦合。前期水文气候条件也构成关键不确定性，生态系统响应不仅取决于同期气候胁迫，还受持续土壤水分异常、植被驯化和延迟生理响应形成的干旱记忆影响。NRB跨越全球最强干旱–湿润气候梯度之一，从东非湿润热带源头到东北非超干旱沙漠，是研究生态系统对CDHW响应的理想自然实验室，但流域尺度上这些环境梯度如何调控WUE对CDHW的响应仍缺乏评估。因此，研究人员将复合事件刻画、生态系统碳水分析、干旱记忆评估、气候归因和生态系统轨迹分析整合到统一框架中，研究CDHW如何跨植被类型、气候梯度和事件序列重塑NRB的生态系统WUE。

研究人员开展的研究及总体结论：研究目标包括：刻画CDHW事件的时间演变和空间分布；利用回归、复合异常、百分比变化和事件分析量化WUE敏感性空间异质性；评估多时间尺度水文气候控制和干旱记忆对WUE的影响；利用相空间轨迹分析识别潜在生态系统状态转变和脆弱性变化。总体而言，CDHW事件在NRB加剧，热浪领先型事件增加快于干旱领先型事件；CDHW通常较单一干旱或热浪造成更强的GPP和ET降低；WUE响应具有明显空间异质性，受多月干旱记忆和气候梯度调控；后期CDHW暴露增加且WUE降低，表明CDHW–WUE耦合状态发生时间转变。该研究为NRB流域尺度CDHW生态影响提供综合评估，并强调在气候变化下需考虑复合极端、前期水文气候条件和气候梯度对生态系统碳水耦合及脆弱性的影响。

主要关键技术方法：研究区NRB约3.25 million km²，范围约4°S–31°N、25°–40°E，跨越11个国家。数据来源包括MODIS/Terra Version 6.1的GPP（MOD17A2HGF.061）和ET（MOD16A2GF.061），8天、500 m，2000–2025聚合为月值和0.25°；ERA5-Land提供三层土壤水分（0–7、7–28、28–100 cm，厚度加权）、2 m气温、露点、太阳辐射、风速和日最高温T_max，并计算VPD；CHIRPS v2.0提供日降水，hPET代表大气蒸发需求并用于计算标准化降水蒸散指数（SPEI）；ESA CCI土地覆盖300 m稳定像元保留农田、林地、草地和灌木。方法包括热浪识别（暖季T_max第90百分位且至少持续3天）、干旱识别（SPEI≤?1.0）、CDHW重叠识别及热浪领先、干旱领先和同步分类；WUE=GPP/ET；LOWESS去趋势、OLS回归、FDR校正、Moran's I、Welch t检验、Pearson相关与1–12月SPEI累积尺度、Pettitt变点检验、相空间轨迹分析，以及随机森林（RF）和SHAP。

研究结果：
3.1 CDHW特征的时间趋势与变异（Temporal trends and variability of CDHW characteristics）。通过年际趋势分析，2000–2025年NRB的CDHW频率、持续时间、严重性（CDHWs）和异常（CDHA）总体增加。频率增加0.92 events/decade（p=0.030），持续时间增加1.08 days/decade（p=0.026），严重性增加2.11 units/decade（p=0.046），CDHA增加0.06 units/decade（p=0.035）。频率范围为1.8–8.6 events/year，2024年达峰；流域平均空间范围为32.6%（26.0%–39.3%）；持续时间平均4.04–12.24 days；严重性3.9–20.1；CDHA 0.82–1.32，均在2024年达到最大值。

3.2 CDHW事件的时空增强与扩展（Spatiotemporal intensification and expansion of CDHW events）。空间分析显示CDHW特征高度异质，正Sen's斜率占优势。在非零Sen's斜率网格中，正趋势占CDHA 88%、持续时间87%、频率79%、严重性74%；统计显著趋势（P<0.05）分别占有效网格的39.5%、32%、26%和17.5%。2014–2025年较2000–2013年，热浪领先型CDHW空间范围增加29.5%，从23.36%升至30.26%；干旱领先型增加12.9%，从32.54%升至36.72%；总CDHW空间范围增加12.9%。空间扩展趋势从0.67±0.86% yr^?1升至1.33±1.51% yr^?1；Pettitt检验识别2014年为CDHW频率显著变点（p=0.037）。

3.3 GPP、ET和WUE对干旱、热浪和CDHW事件的响应（GPP, ET, and WUE responses to drought, heatwaves, and CDHW events）。空间上，GPP在约93.1%流域增加，ET在约92.5%流域增加，而WUE在约58.6%流域下降。流域平均GPP为8.209 g C m^?2 yr^?1、ET为4.722 mm yr^?1，均显著增加；WUE呈弱下降，为?0.009 g C m^?2 mm^?1 yr^?1。极端事件下，GPP中位异常分别为干旱?83.72%、热浪?70.27%、CDHW?84.45%；ET中位异常分别为?89.98%、?86.53%、?93.49%；WUE中位异常在干旱、热浪和CDHW下分别为17.35%、12.22%和23.07%，且差异显著。植被类型中，农田GPP在干旱、热浪和CDHW下分别下降?62.13%、?48.88%和?71.15%；林地、草地和灌木下降更强，草地CDHW下降低达?89.15%；ET在灌木和草地CDHW下降最大，分别为?96.17%和?92.22%；WUE在各植被类型中多为正响应。气候区分析显示，干旱和半干旱区GPP下降最大，CDHW下分别达?94.30%和?91.49%；湿润区下降较小，干旱下?46.52%、CDHW下?55.81%；ET呈相似格局；WUE在干旱区干旱条件下中位异常为?14.92%，在半干旱、干性半湿润和湿润区多为正响应。

3.4 生态系统WUE对CDHW的空间敏感性、复合响应和事件影响（Spatial sensitivity, composite response, and event-based impact of ecosystem WUE to CDHW）。空间回归显示WUE敏感性空间异质：45%流域为负敏感，55%为正敏感；FDR校正后47.38%有效网格的CDHW频率与WUE关系显著；残差Moran's I为0.605–0.816，均值0.706（p<0.001）；流域平均回归系数弱负为?0.04。高CDHW年与低CDHW年复合分析表明，约49%流域出现负WUE异常，流域平均异常?0.12，最小?14.88。百分比变化分析显示49%区域WUE负变化、50.81%正变化，流域平均WUE下降约2.36%，局部下降达?288.67%。事件影响评估显示67%流域为正影响、33%为负影响，流域平均弱正响应0.16，局部范围为?3.45至+1.94；负响应集中在水分限制区，正响应见于其他水文气候区。

3.5 CDHW下与生态系统WUE相关的前期水文气候时间尺度（Antecedent hydroclimatic timescale associated with ecosystem WUE under CDHWs）。通过WUE与SPEI最大累积相关分析，正相关占84.2%，负相关占15.8%；中等至强干旱–WUE关系（Rmax_cum>0.3）覆盖53.5% NRB，较强耦合主要在埃塞俄比亚高原和中东部NRB。主导时间尺度为1–12个月，长记忆8–12个月占约50.0%，短记忆1–4个月占26.2%，中等5–7个月占23.7%。植被类型中，林地干旱–WUE耦合最强（Rmax_cum=0.41），农田0.33，草地0.29，灌木0.29；林地记忆约8个月，灌木约7个月。气候区中，湿润区耦合最强（Rmax_cum=0.37），记忆约8个月；干旱区记忆约6个月。

3.6 CDHW极端下生态系统WUE的时间转变（Temporal shift in ecosystem WUE under CDHW extremes）。相空间轨迹分析显示，热浪领先型CDHW下WUE与CDHW空间范围呈强负相关（R=?0.70，P<0.001）；干旱领先型也呈显著负相关（R=?0.60，P=0.001）。轨迹和质心位移表明，2000–2013年向2014–2025年转变，从较高WUE和较低CDHW范围转向较低WUE和较高CDHW范围；Pettitt检验识别2014年为CDHW频率显著变点。结果表明CDHW–WUE耦合状态发生时间转变，后期CDHW暴露更高、WUE更低。

3.7 与生态系统WUE相关的气候预测因子（Climatic predictors associated with ecosystem WUE）。RF和SHAP分析显示，RF测试集R²为0.70–0.81。干旱区中，降水重要性0.23，其次VPD 0.18、Rs 0.17、Ta 0.15、SMC 0.14、WS 0.13；半干旱区降水0.30、VPD 0.20、WS 0.15、Rs 0.15；干性半湿润区降水0.23、Ta 0.22、Rs 0.18、VPD 0.15；湿润区降水0.24、Ta 0.21、Rs 0.18、VPD 0.18、WS 0.07。SHAP显示降水增加通常与更高WUE相关，VPD和温度升高与更低WUE相关；这些为模型预测贡献和关联，并非独立因果效应。

讨论部分总结：讨论指出，2000–2025年NRB的CDHW事件明显加剧，频率、持续时间、严重性、异常幅度和空间范围均增加，这与增暖趋势、大气蒸发需求增加、VPD升高和陆–气相互作用增强一致。热浪领先型CDHW相对干旱领先型快速扩展，表明复合极端事件由降水控制向温度驱动转变。空间上，CDHW增强在中央和南部NRB尤为明显，过渡水分条件促进强陆–气耦合。CDHW通常比单独干旱或热浪造成更大GPP和ET下降；WUE因ET下降比例更大而稳定或增加，但CDHW下WUE增强弱于单独干旱，说明极端热胁迫限制植被维持碳同化的能力。事件尺度WUE正响应与长期WUE下降反映不同时间维度：急性胁迫可暂时提高表观WUE，反复或持续胁迫则可能限制碳同化和恢复。草地和灌木GPP、ET下降最大，林地较小，农田中等；干旱和半干旱区比湿润区更敏感。前期水文气候时间尺度分析表明，林地干旱–WUE耦合最强且记忆较长，干旱区时间尺度较短；但SPEI累积尺度重叠，识别的最强相关应解释为相对水文气候条件指标，而非确定生理记忆。RF和SHAP显示，干旱和半干旱区降水对WUE预测贡献最大，湿润和亚湿润区温度、辐射和VPD贡献更均衡；降水增加与更高WUE相关，VPD和温度升高与更低WUE相关，但需视为模型关联。管理上，讨论建议将干旱、热浪、土壤水分、SPEI和WUE纳入综合预警，关注WUE持续下降和CDHW敏感区，采取灌溉效率、耐旱作物、土壤保水、植被保护与恢复等措施，并加强跨境监测、数据共享和区域合作。不确定性包括MODIS GPP和ET参数及气象输入误差向WUE传播，CHIRPS、ERA5-Land、hPET和SPEI阈值不确定性，灌溉、水库、地下水开采未定量纳入，空间依赖影响显著性估计，以及未开展正式不确定性传播和多产品集合分析。

研究结论翻译：本研究对NRB的CDHW趋势及生态系统响应进行了综合评估。研究人员记录了CDHW事件明显加剧以及向热浪主导序列的转变，这与更广泛变暖趋势和大气需求上升一致。生态系统对复合极端事件表现出显著响应：CDHW期间GPP和ET受到明显抑制，而WUE响应具有空间异质性。特别是半干旱和干性半湿润区对复合水文气候胁迫表现出明显生态系统敏感性。CDHW与WUE的耦合在很大程度上取决于前期干旱条件，多月累积SPEI在整个流域大部分地区与生态系统响应具有最强关联。相空间轨迹分析进一步揭示，后期向更高CDHW暴露和更低WUE的时间转变，表明在复合事件暴露增加下CDHW–WUE耦合状态发生变化。机器学习分析进一步揭示WUE预测贡献的气候梯度：干旱区降水贡献更大，湿润区温度、VPD和辐射贡献更均衡。总之，研究人员的发现表明，加剧的复合干旱和热浪与生态系统碳水动态的显著变化相关，强调未来生态水文评估需要考虑复合极端事件、前期水文气候条件和气候梯度。