严重干旱后美国西南部树木生长轨迹的复杂性

《Ecosphere》：Complex trajectories of tree growth in the southwestern United States after severe drought

　　

气候变化驱动的极端干旱对森林生态系统具有重大影响，包括大规模树木死亡和初级生产力下降，并进一步反馈影响全球碳循环。然而，极端干旱事件对森林死亡、生态系统响应以及恢复/干旱后轨迹的长期影响目前尚不清楚。本研究结合美国西南部五种主要树种的逐年树轮宽度数据以及长期森林清查与监测样地数据，分析了2002年极端干旱事件对树木后续生长的影响。研究人员量化了存活树木因林分密度降低带来的潜在资源增加而生长加速的程度，以及因干旱胁迫持续影响而导致生长减缓的程度。研究发现，不同物种的干旱后生长轨迹呈现多样化模式：部分物种如颤杨（Populus tremuloides）出现生长大幅增加的明显迹象，而其他物种如美国黄松（Pinus ponderosa）则表现出清晰的生长抑制，反映了显著的干旱遗留效应。总林分断面积（total basal area）是预测干旱后生长响应的最佳指标，但其回归效应（正向或负向）因物种而异；例如，美国黄松在较高密度林分中生长低于预期，而云杉在较高密度林分中生长则高于预期。气候水分亏缺（climatic water deficit, CWD）和林分年龄也被发现是多种物种干旱后生长轨迹的重要驱动因素。本研究结果有助于阐明美国西南部不同森林类型对未来干旱事件的响应方式及其对该区域碳循环的影响。

研究背景与问题

气候变化正驱动全球范围内极端干旱事件的频率和强度持续增加，这对森林生态系统构成了严重威胁。严重干旱往往因温度升高和大气干燥度增加而加剧，可导致树木生长下降、大规模树木死亡，并影响森林碳汇功能。尽管干旱诱导的树木死亡在全球多种森林类型中已有广泛观测，但当前科学界对于这些死亡事件在长期生态和碳循环效应方面的认识仍存在不确定性。特别是极端干旱后森林恢复动态、物种特异性生长响应以及气候对恢复过程的影响机制，尚需深入研究。森林作为陆地碳汇的重要贡献者，其响应气候变化的能力直接关系到全球碳循环的走向，而改进动态全球植被模型（dynamic global vegetation models, DGVMs）以准确预测森林碳通量，是当前生态系统研究的关键任务之一。

与此同时，干旱后存活树木的生长轨迹可能呈现复杂模式：一方面，邻近树木死亡可能通过竞争释放（competitive release）效应增加存活树木的资源可获得性；另一方面，干旱可能造成的水力损伤（hydraulic damage）可能在树木存活后长期持续，导致生长抑制。这种"干旱遗留效应"（drought legacy effects）的方向和强度因物种功能特性和林分条件而异，但相关机制尚不明确。为此，研究人员以美国西南部2002年严重干旱事件为研究对象，旨在量化干旱后树木生长轨迹的复杂性，并解析气候变量、物种特异性状和林分特征如何调控干旱后生长遗留效应的方向与幅度。

关键技术方法

研究人员整合了美国林务局森林清查与分析（USFS Forest Inventory and Analysis, FIA）项目在亚利桑那州、科罗拉多州和犹他州设立的33个永久样地的数据，以及五种主要树种（美国黄松、恩氏云杉、颤杨、单籽松、犹他圆柏）的树轮宽度记录。气候数据来源于TerraClim高分辨率全球数据集（~4 km空间分辨率），提取了气候水分亏缺、最高温度、年降水量、雪水当量、水汽压亏缺、实际蒸散量、土壤含水量和帕默尔干旱 severity指数等变量。研究采用线性混合效应模型建立1960–1989年基线期树轮宽度与气候变量的关系，以此预测2000–2014年干旱后的"预期生长"，通过比较观测生长与预期生长的差异（即残差）来量化干旱遗留效应。进一步构建允许截距和斜率随物种变化的混合效应模型（Model B），分析了气候水分亏缺、总林分断面积和林分平均年龄对干旱后生长轨迹的调控作用。

研究结果

统计气候生长模型与树轮宽度残差分析

研究人员为每个物种建立了最佳气候–生长预测模型。基线期模型显示，帕默尔干旱 severity指数是所有五个物种最终模型中的共同变量，表明该指数对树木生长具有一致的重要预测作用。将基线模型应用于2002年干旱后年份时，发现了显著的物种间差异：颤杨在干旱后3年、5年和10年均表现出显著的正残差（p < 0.001），暗示竞争释放效应；犹他圆柏在干旱后3年和5年同样显示正残差（p < 0.01），但10年后效应消失；单籽松各时期残差均不显著；而美国黄松和恩氏云杉在所有三个时期均显示显著负残差（p < 0.001），表明存在持续的干旱遗留损伤。时间序列比较显示，恩氏云杉和美国黄松的生长抑制随时间略有减弱（10年残差幅度小于3年和5年），但仍未恢复至预期水平。

干旱后生长未解释变异的混合效应模型分析

在三个解释变量（气候水分亏缺、总林分断面积、平均林分年龄）的分析中，允许物种间截距和斜率均变化的Model B均为最优模型。气候水分亏缺模型的条件R²为0.73，该变量对恩氏云杉、美国黄松和颤杨的残差有显著影响。总林分断面积模型的条件R²最高（0.78），对恩氏云杉、美国黄松、单籽松和颤杨均有显著或边际显著影响。林分年龄模型的条件R²为0.70，对恩氏云杉和美国黄松有显著影响。

死亡物种身份对遗留效应的影响

研究还发现，样地中死亡物种的身份部分调控了存活树木的干旱后生长变化，尤其在松–桧混交林中表现明显。当死亡物种为异种（通常为单籽松）时，存活的犹他圆柏表现出 higher than expected 的生长；而颤杨在同种个体死亡为主时，其生长增幅更大。值得注意的是，美国黄松在未记录死亡的样地中反而表现出最大的生长下降，这可能与样地调查半径限制未能捕捉邻近死亡率有关，也可能反映了竞争效应——即全部存活的林分中树木因冠层未稀疏而生长较差。相反，恩氏云杉在仅记录有云杉死亡时，其预期生长下降幅度最大。

讨论与结论

研究人员指出，美国西南部森林的干旱后树木生长轨迹复杂且高度物种特异性，反映了物种特性、干旱胁迫和竞争环境之间的广泛相互作用。与采伐研究中的竞争释放效应相比，本研究中正向生长异常的幅度通常较小，这可能源于2002年和2006年干旱事件的特征差异、物种间干旱耐受性变异，以及干旱造成的长期生理损伤。

物种特异性响应方面，更耐旱的犹他圆柏相较于单籽松表现出更高的相对干旱后生长，这可能是因为单籽松在干旱中承受了更多的功能性木质部损伤，无法充分利用改善后的水分条件。颤杨的正向生长异常则可能归因于其快速产生新木质部的能力，从而促进干旱后恢复。然而，美国黄松和恩氏云杉的持续生长抑制表明，干旱遗留损伤在某些物种中可能长期存在。

驱动因子分析揭示了林分结构因子的关键作用。总林分断面积作为竞争强度的代理指标，被发现是预测干旱后生长响应的最佳指标，但其效应方向因物种而异：恩氏云杉和颤杨在较高密度林分中生长高于预期，可能反映了竞争缓解效应；而美国黄松则呈现相反模式。林分年龄对美国黄松和恩氏云杉的干旱后生长有显著介导作用，但方向相反——恩氏云杉呈现正相关，而美国黄松呈现负相关，这与美国黄松在人工林管理中较早间伐的实践一致，也反映了老龄未管理林中竞争动态的独特影响。气候水分亏缺与物种干旱后生长的关系同样呈现位点依赖性：恩氏云杉在较湿润（低气候水分亏缺）位点表现出更大的生长下降，可能与高海拔森林的能量限制特征有关；而美国黄松则在较干旱位点遭受更严重的生长抑制。

研究结论部分指出，干旱后树木生长轨迹的过程复杂且高度物种特异性。更好地理解复杂生态水文过程及其在干旱期间和之后对树木水分资源的影响，是未来研究的重要方向。水分循环过程的变化——如土壤蒸发、径流、地下水补给、冠层截留和蒸腾——可能导致进一步加剧干旱死亡风险的反馈效应。另一方面，由于邻近树木死亡导致的竞争释放所产生的正向生长异常，可能有助于抵消未来干旱和虫害对存活树木的影响。本研究揭示了美国西南部物种间和种内干旱响应的变异性，其结果有助于指导进一步研究，聚焦于大规模生物物理过程和反馈如何由个体树木对干旱的响应所塑造。