基于净辐射的生态系统水分—能量限制框架揭示温室增暖下生态系统水分限制无净增加

《Earth's Future》：No Net Increase in Ecosystem Water Limitation Under Greenhouse Warming

　　

能量和水分可利用性是控制陆地生态系统功能（如蒸散发）的关键因素。准确识别蒸散发受水分或能量限制的主导机制，对于理解生态系统对气候变化的响应至关重要。近期基于地球系统模式（ESM）预估的研究表明，在全球增暖背景下，生态系统将广泛从能量限制型转变为水分限制型。研究人员证明，采用热力学上更为合理的能量指标将从根本上改变这一预估结果。地表用于蒸散发及相关生态系统功能的能量可利用性主要由净辐射（Rn）决定，而非下行短波辐射或气温。基于这一以Rn为核心的生态系统限制框架，研究人员发现，到本世纪末，蒸散发和蒸腾的水分限制均无显著的净增加趋势。理想化的ESM CO2试验表明，这一格局归因于CO2升高对植被产生的生理效应——主要是气孔调节和水分利用效率（WUE）增强——在很大程度上抵消了温室增暖下增强的辐射强迫所导致的水分限制加剧。关键生理变量如净初级生产力（NPP）和叶面积指数（LAI）在过渡地带对基于Rn的水分—能量限制格局转变也表现出协调一致且机制上可解释的响应。这些发现与其他基于生态水文研究的结论一致，并得到CO2升高条件下植被对干旱条件敏感性降低的经验证据支持。该研究弥合了生态学与物理学理论之间的鸿沟，改进了生态系统水分—能量限制分析，并为未来生态系统动态提供了清晰的机制理解。

论文解读

一、研究背景与科学问题

陆地生态系统的功能（如蒸散发、初级生产力等）同时受能量和水分可利用性的共同限制。准确判断生态系统处于水分限制还是能量限制主导的状态，是预测生态系统对气候变化响应的核心前提。然而，学术界在“全球增暖是否导致陆地生态系统水分限制加剧”这一问题上长期存在激烈争论。传统上，基于Penman、Thornthwaite和Penman-Monteith等潜在蒸散发（PET）方程构建的干旱指数（AI）普遍预测全球变暖将导致广泛的大气干燥趋势和干旱区扩张。但从生态水文学视角出发，观测到的植被绿化现象以及CO₂升高导致的气孔导度下降，表明陆地生态系统正变得更加高效地利用土壤水分，对大气干旱的抵御能力也在增强。

这一争论的关键症结在于：在生态系统限制分析中，究竟应该用什么变量来表征能量可利用性。近期多项研究使用下行短波辐射（SW_in）和气温（T_a）作为能量指标，基于生态系统限制指数（ELI）预测了全球生态系统将广泛从能量限制转向水分限制。然而，根据李比希最小因子定律，一个变量要成为生态系统功能的限制因子，必须同时满足与响应变量正相关且对其构成上界约束这两个条件。观测数据表明，SW_in或T_a与蒸散发（ET）之间存在广泛的负异常相关，违反了这一定律。因此，采用SW_in作为能量指标可能严重高估了未来生态系统的水分限制趋势。

二、研究内容与核心发现

针对上述问题，研究人员提出了以净辐射（R_n）为能量指标的生态系统限制框架，构建了基于R_n的生态系统限制指数（RELI），利用16个CMIP6模式的历史及SSP5-8.5情景模拟数据（1981–2100年），系统分析了全球暖区陆地生态系统水分与能量限制格局的时空演变趋势及其驱动机制。研究得出以下核心结论：全球尺度上，基于R_n的ET和蒸腾（T）限制指数在整个21世纪均无统计显著的净增加趋势；水分限制和能量限制区域之间呈现双向转换格局，净变化不显著；CO₂的生理效应（气孔调节和水分利用效率增强）在很大程度上抵消了辐射强迫增暖所导致的水分限制加剧。

三、主要技术方法

研究人员采用的主要技术方法包括：其一，构建基于净辐射的生态系统限制指数（RELI），以土壤湿度（SM）或降水（P）为水分可利用性指标、以R_n为能量可利用性指标，通过计算各变量月尺度异常值之间的Kendall秩相关系数之差来判定水分限制（RELI>0）或能量限制（RELI<0）状态；其二，使用CMIP6框架下16个模式的历史和SSP5-8.5情景模拟数据，数据统一重采样至2°×2°分辨率；其三，利用三组理想化CO₂试验（1pctCO2、1pctCO2-rad和1pctCO2-bgc）分离CO₂辐射强迫与生理强迫的独立贡献，该分析基于6个输出完整的模式（BCC-CSM2-MR、CMCC-ESM2、CNRM-ESM2-1、IPSL-CM6A-LR、MIROC-ES2L和UKESM1-0-LL）；其四，采用Theil–Sen斜率估计和Mann–Kendall检验分析趋势显著性，并应用Benjamini-Hochberg错误发现率（FDR）校正控制多重检验的假阳性率；其五，使用GLASS AVHRR叶面积指数（LAI）卫星数据集（1981–2018年）对模式历史模拟进行独立验证。

四、主要研究结果

3.1 全球水分与能量限制趋势

基于ET的RELI（RELI_{Rn_ET}）在1981–2100年间净增加仅为0.015±0.0453（p>0.05），无统计显著性；基于T的RELI（RELI_{Rn_T}）净增加为0.0216±0.0534（p>0.05），同样不显著。相比之下，使用SW_in作为能量指标的原始ELI（ELI_{SWin_ET}）则呈现显著上升趋势（0.0582±0.0538，p<0.05）。研究发现，ELI_{SWin_ET}的上升主要来源于cor(SW_in′, ET′)的下降，而这反映了SW_in与SM在干旱区的耦合以及SW_in与饱和水汽压差（VPD）在湿润区的正协方差所带来的混淆效应，即“双重计数”问题。

3.2 水分与能量限制的驱动因素

全球暖区陆地平均R_n在1981–2100年间增加了0.255 mm/d，远大于SW_in的增加幅度（0.0533 mm/d）。其中净长波辐射（LW_n）增加0.151 mm/d，主要源于温室气体增暖导致的下行长波辐射（LW_in）增强；净短波辐射（SW_n）增加0.104 mm/d，与地表反照率降低（冰雪消融和植被绿化）有关。土壤湿度（SM）持续下降，但ET持续上升，表明ET并未在全球尺度上日益受到SM限制。蒸腾（T）在整个世纪基本无变化，其与ET的差异反映了地面蒸发和冠层截留蒸发的增加。CO₂诱导的气孔关闭（生理节水效应）在叶片尺度上抵消了LAI扩张带来的结构效应和大气蒸发需求增强，从而稳定了T；而物理蒸发（ET-T）受大气蒸发需求驱动持续增加，推高了总ET。

3.3 水分与能量限制格局及其转换

基于RELI_{Rn_ET}，约38.4±16.2%的暖区陆地面积保持水分限制状态，多数区域保持能量限制状态。两个时期（1981–2010年与2071–2100年）的对比揭示了双向转换格局：半干旱区（如萨赫勒和美国西部）出现从水分限制向能量限制的转换（占暖区面积8.9±5%），而湿润区边缘（如巴西东北部和南非）出现从能量限制向水分限制的转换（占8±2.6%）。净效应为水分限制面积仅减少0.9±5.5%（p>0.05），无统计显著性。这一格局与基于CO₂校正的湿波文比干旱指数（AI_{Zhou_CO2}）以及Berg和McColl（2021）估算的全球旱地面积净减少趋势一致。相比之下，ELI_{SWin_ET}预测69.9±7.6%的面积处于水分限制状态，且能量向水分限制的单向转换显著增加5.5±4.2%（p<0.05），明显高估了水分限制程度。在过渡地带，净初级生产力（NPP）、LAI和水分利用效率（WUE）均与RELI呈负异常相关，表明水分限制加剧会抑制生态系统功能，而WUE的下降幅度相对较小，反映了CO₂升高下气孔调节对植被干旱敏感性的缓解作用。

3.4 RELI对CO₂生理与辐射强迫的机制归因

理想化CO₂试验表明，辐射强迫（1pctCO2-rad）持续加剧水分限制，表现为ET和T对SM的依赖性增强，这是增暖导致大气蒸发需求升高的典型特征；在无CO₂气孔调节的情况下，T逐渐与能量可利用性脱钩而更直接受SM约束。生理强迫（1pctCO2-bgc）则显著缓解水分限制，尤其是T基RELI显著下降，这归因于气孔关闭和WUE增强使T对SM的依赖性降低。在完全CO₂强迫（1pctCO2）下，生理效应与辐射效应及其交互作用相互拮抗，导致两种RELI均无显著长期水分限制增强趋势（RELI_{Rn_ET}：0.015±0.0256，p>0.05；RELI_{Rn_T}：0.0044±0.050，p>0.05）。生理效应通过气孔调节有效抑制了蒸腾对大气蒸发需求的响应，而辐射增暖通过加速物理蒸发成为推高总ET的主要物理驱动力。

五、讨论与结论

本研究的核心贡献在于揭示了以往基于SW_in或气温的生态系统限制分析存在根本性的方法缺陷——这些变量与ET之间存在负相关混淆，导致水分限制被系统性高估。采用热力学上更为合理的R_n作为能量指标后，全球暖区陆地生态系统在21世纪不会出现显著的水分限制净增加。这一结论与Berg和McColl（2021）关于无旱地扩张的生态水文学证据、Greve等人（2014，2019）关于有限大陆干燥化的研究以及Milly和Dunne（2016）的发现一致，共同挑战了传统Penman型PET和AI指数所预测的广泛陆地表层干旱化趋势。理想化CO₂试验提供的机制证据表明，CO₂升高通过气孔调节和水分利用效率增强等生理途径，在很大程度上抵消了辐射强迫增暖所预期加剧的水分限制。该研究强调了在评估温室增暖下植被陆地水分限制或干旱化趋势时，必须充分考虑生态系统对CO₂升高的生态水文响应，这一发现对生态系统管理和气候适应策略具有深远意义，尤其是在限制格局转换最可能发生的过渡地带。