《Frontiers in Plant Science》：Dynamic stability of climate-growth relationships in Picea crassifolia in the Qilian Mountains: the modulating role of elevation gradient and time-varying characteristics
编辑推荐：
本研究聚焦于祁连山干旱半干旱区，探讨在全球气候变暖背景下，优势树种青海云杉(Picea crassifolia)径向生长对气候变化的响应。研究发现，其气候-生长关系具有显著的海拔差异性和时间不稳定性，中高海拔位点（3100-3300 米）尤为敏感。这揭示了高山森林生态系统对气候变化的响应是动态的，而非静态的，为理解山区生态系统的适应机制及制定适应性管理策略提供了新视角。
材料与方法
本研究在祁连山排露沟小流域内，沿海拔梯度（2900、3000、3100、3200、3300 米）设置了五个采样点。在每个海拔点建立一个20 米 × 20 米的永久样地，共采集了75棵健康优势或亚优势青海云杉的150个树芯样本。所有样本在实验室经过干燥、固定、打磨等标准树木年代学处理流程。使用LINTAB? 6高精度测量系统测量年轮宽度，并使用COFECHA程序进行交叉定年质量控制。最终，利用ARSTAN软件包，通过负指数曲线或负斜率的线性回归去除树木自身的年龄趋势，构建了标准化的残差年表用于后续分析。
气候数据来源于CRU TS（气候研究单位时间序列）网格化数据集，提取了1960年至2021年研究区域（38.52°-39.02°N, 100.27°-100.77°E）的月平均、最高、最低温度和月总降水量。统计分析方面，通过计算五个海拔年表间的滑动Pearson相关系数以评估生长同步性；通过静态Pearson相关和30年滑动窗口相关分析，评估了年表与从前一年9月到当年10月的月气候因子之间的相关强度及其时间稳定性。年表可靠性通过信噪比(SNR)和表达群体信号(EPS)等指标评估，所有年表的EPS值均大于0.85，满足分析要求。
结果
3.1 不同海拔青海云杉年表特征
年表统计特征显示，海拔3300米的年表平均敏感度(MS=0.155)最高，表明树线位置树木生长对年际环境波动最敏感。海拔3200米的年表一阶自相关(AC1=0.760)最高，表明前一年气候对当年生长的遗留效应最强。海拔3100米的年表表现出最高的序列间平均相关系数、信噪比(SNR)和表达群体信号(EPS)，表明其包含最强的共同气候信号，最适于检测区域气候信号。
3.2 不同海拔青海云杉树轮宽度与气候因子的关系
生长同步性分析表明，海拔2900米与3000米年表间同步性最强(r=0.547, p<0.001)，而海拔3000米与3200米年表间无显著关联，表明可能存在海拔阈值或生态过渡带使生长响应解耦。
气候-生长响应存在显著海拔分异。在低海拔（2900-3000米），树木生长主要受水分有效性限制，与当年1月和9月降水呈显著正相关(p<0.01)，而与12月温度呈显著负相关。在高海拔（3200-3300米），树木生长对冬季和早夏条件更敏感，与前一年12月和当年6月降水呈显著负相关(p<0.05)，与当年6月温度也呈显著负相关，表明冬季温暖和早夏过多水分可能抑制生长。海拔3100米处，生长与冬、春季温度呈一致负相关，表明早季变暖可能通过促进雪融或增加蒸发需求而产生负面影响。
3.3 气候-生长关系的时间稳定性
30年滑动窗口相关分析揭示了气候-生长关系普遍存在时间非平稳性，且非平稳性的程度和季节表现沿海拔梯度系统变化。在中高海拔位点（3100-3300米），气候生长响应，特别是在生长季（6-9月）及相邻月份，表现出显著的时间变异性。例如，在3100米，对关键水分补给期（9月）和生长季（6月）的降水响应，以及对前一年9-12月和当年6、8、10月温度的敏感性均表现出显著波动。这种不稳定性可能源于气候变暖驱动的日益加剧的水热失衡。相比之下，在树线位置（3300米），非平稳性相对受限，影响的变量和季节窗口较少，温度作为限制因子的作用表现出比低海拔更大的时间稳定性。
讨论
4.1 海拔梯度对青海云杉生长的影响
传统观点认为，低海拔生长受水分限制，高海拔生长受温度限制。然而，本研究发现，在整个海拔梯度上，青海云杉径向生长与前一年12月、当年1-2月及5-6月的温度均呈负相关，这与经典范式相悖。在高海拔，低温不仅直接限制光合作用，还通过干扰水分和养分吸收、土壤冻结、碳分配等生理过程，以及冻融循环造成的细胞膜损伤来制约生长。此外，高海拔强烈的太阳辐射加剧了蒸腾失水，在祁连山有限的降水和高蒸发需求共同作用下，水分胁迫成为树线附近树木生长的重要限制因素。因此，高海拔观测到的负温度-生长相关主要归因于水分胁迫和冻融动态等间接机制，而非直接的热限制。
在中海拔（约3100米），温度和水分的协同作用更有利于生长，表现出生长最优，符合生态最适性理论。树木形成层的活动直接响应这些环境信号，决定了年轮宽度。
4.2 不同海拔青海云杉径向生长对气候响应的差异
本研究揭示的气候-生长关系时间非平稳性，在中高海拔位点最为明显。这种非平稳性在全球多地高山生态系统中均有报道。其主要机制涉及温度与水分有效性的复杂相互作用。在高海拔地区，关键物候事件（如融雪时间、生长季开始）以及晚春霜冻等，会扰乱形成层活动。温度介导的融雪物候变化改变了土壤水分的时空可用性，从而放大了生长响应的变异性。
在中海拔生态过渡带，温度介导的干旱胁迫和水分有效性相互交织，导致生长响应波动剧烈。在高海拔范围（3200-3300米），生长对冬季寒冷和春季变暖响应的非平稳性，与寒冷温带及高山地区的研究结果一致，表明生长响应不仅受平均气候条件控制，也受晚春霜冻或过早融雪等极端季节性事件的频率和强度影响。
需要指出的是，本研究基于CRU TS网格数据（0.5° × 0.5°分辨率）得出的非平稳性推论存在局限性。尽管该数据在区域尺度上具有良好的一致性，但其较粗的空间分辨率可能无法解析排露沟小流域内复杂地形驱动的局地微气候动态。因此，观测到的部分非平稳性可能源于气候输入数据本身未解决的偏差或平滑效应，而不仅仅是树木生理或生态响应的真实变化。
结论
本研究证明，祁连山青海云杉的气候-生长关系是高度动态且非平稳的，沿海拔梯度变化显著。在低海拔（2900-3000米），树木生长主要受水分限制，与非生长季降水呈正相关。在高海拔（3200-3300米），生长更多受温度制约，与冬季降水和夏季温度呈负相关，反映了对热约束和水分胁迫的敏感性。生长响应的时间非平稳性在中高海拔（3100-3300米）尤为突出，在变暖趋势下，温度与水分胁迫之间的相互作用日益难以预测。这些发现强调了高山森林中生长-气候相互作用日益复杂，表明温度和水分有效性（包括其季节动态）都是树木生长的关键决定因子。本研究凸显了构建适应性森林管理框架的必要性，该框架需明确纳入短期气候极端事件和长期主导气候驱动因素的转变。
