海拔梯度下大黄石生态区柔松（Pinus flexilis）千年生长的水文气候约束分异

《Journal of Biogeography》：Elevation-Dependent Hydroclimatic Constraints on Millennial Growth of Pinus flexilis Across Its Ecological Amplitude in the Greater Yellowstone Ecoregion

　　

目的：检验柔松（Pinus flexilis）在其占据的区域海拔梯度上，气候—生长关系、时间一致性和对复合气候极端事件的响应如何变化。地点：大黄石生态区（Greater Yellowstone Ecoregion），美国怀俄明州西北部。时间范围：定年树轮材料跨度为514–2020 ce，详细分析聚焦于过去千年。主要研究分类单元：柔松（Pinus flexilis James）。方法：研究人员在高、中、低海拔样点建立了三条千年长度的轮宽年表，将径向生长与器测及重建气候资料进行比较，利用主成分分析（PCA）和移动窗口分析量化共享的与海拔依赖的生长格局，并评估生长对复合干热与冷湿极端事件的响应。结果：各海拔生长均对水分敏感，并在暖夏季条件下下降，表明热量在整个梯度上放大了水分限制。各年表共享宽泛的水文气候信号，但高海拔与低海拔样点在生长幅度和生长响应的时间上常有差异。高海拔Wolf Knob在冷期小冰期条件下暂时与主导生长格局解耦，而低海拔Outlaw Ridge在干热极端事件后表现出最持久的生长抑制。主要结论：在大黄石生态区所采样的生态幅度内，柔松生长既反映共享的水文气候变率，也反映海拔特异的限制。未来增温可能在某些年份减弱林线附近的低温限制，但很可能加剧低海拔的水分限制与干旱遗留效应。

研究背景：植物对环境变化的响应取决于气候强迫如何转化为生理限制。在山地森林中，这种转化很少在物种分布区内均匀发生，因为温度、生长季长度、积雪、蒸发需求和土壤水分随海拔急剧变化。持续增温预计将加剧干旱胁迫和复合干热极端事件，但高海拔环境（低温与短生长季限制生产力）与低海拔环境（水分限制和大气水分需求更强）的后果可能不同。北美西部干旱胁迫已是森林的主要压力，预计的干旱胁迫水平可能超过过去1000年的大旱，1997年以来干旱已影响约2000万公顷森林，持续增温预计将增加未来森林死亡。柔松耐旱、具可塑性且耐受气候胁迫，分布于内陆西部从低林线、干性山地、亚高山到高林线的广泛地形与气候环境，近几十年死亡率上升并受山松甲虫暴发扰动。大黄石生态区具有陡峭海拔梯度、半干旱大陆性气候和强烈的温度与降水季节对比，上下林线均对气候变化敏感，且柔松个体常超过1000年，其树轮记录可揭示环境限制、时间一致性和干旱恢复如何随反复气候胁迫而变化。以往研究表明柔松可在包括大旱和大湿在内的重大气候极端中持续存在，且几乎无上坡迁移证据，但这种持续性是否在各海拔上表现相似尚不清楚，特别是该区域采样生态幅度内的生长究竟反映共享水文气候信号、海拔特异限制，还是在温度与水分限制交互方式不同的时期发生时间重组，仍是待解决的生理与生态问题。因此研究人员建立三条千年长度柔松树轮年表，提出三个问题：（1）不同海拔种群的径向生长对温度、降水和干旱的响应如何变化；（2）样点间生长格局随时间的一致性如何；（3）复合干热与冷湿极端是否产生海拔依赖的生长遗留效应。

研究方法关键技术：研究人员在大黄石生态区选取Wolf Knob（2939 m）、Sunlight Basin（1956 m）和Outlaw Ridge（1600 m）三个柔松样点，2014—2021年间采集活树与死树样本，共265株（Outlaw 72株、Sunlight 104株、Wolf 89株）；活树用5 mm生长锥胸高处取芯，每株两芯，死树用链锯取薄横截面。室内样品经逐级砂纸打磨、3200 dpi扫描后用CooRecorder测量轮宽，视觉交叉定年并以COFECHA统计检验，用R包dplR以年龄依赖平滑样条在无信号框架下标准化去趋势。统计分析将年表与月降水、月最高温度（T_max）、Palmer干旱强度指数（PDSI）及自校准Palmer干旱强度指数（scPDSI）比较，器测数据取自PRISM与KNMI Climate Explorer，长期比较使用重建的8月T_max和北美干旱图集（NADA）重建6—8月PDSI，共同区间为968–2003 ce；采用50年滚动相关、主成分分析（PCA）、100年移动窗口PCA、同步性与平均绝对差（MAD）分析、年际成对生长差异及20年LOESS平滑，并用双自举叠加纪元分析（SEA）和双变量事件分析（BEA）评估复合极端事件响应。

研究结果：

3.1 年表建立与气候—生长关系：建立Wolf 1281年、Sunlight 1047年、Outlaw 1272年的三条千年长度年表；Sunlight的样点间相关（RBAR=0.75）和平均敏感度（0.46）最高，Outlaw次之（RBAR=0.70，敏感度0.39），Wolf最低（RBAR=0.61，敏感度0.26）。月气候相关显示显著但季节多变的温度响应：Wolf是唯一与上一年最高温度显著正相关的年表，尤其上一年6—9月（r=0.29–0.47）；当年生长季最高温度与Wolf 6月和7月（r=?0.36）、Outlaw 6月（r=?0.32）、Sunlight 7月（r=?0.31）生长负相关。降水响应更具季节性和样点特异性：Wolf和Outlaw均与当年6、7月降水显著正相关（r=0.28–0.31），Outlaw还与上一年6月降水正相关（r=0.35）；Sunlight最强降水响应在5月（r=0.37）及上一年10月—当年9月水文年（r=0.59）；Wolf与上一年7月和9月降水显著负相关（r=?0.30至?0.39）。scPDSI相关最强且最一致，水文年scPDSI相关在Wolf为0.51、Sunlight为0.74、Outlaw为0.59。

3.2 共同与海拔依赖的生长模态：PCA将共同期生长变率分离为两个主导模态，合计解释79.7%方差。PC1解释47.2%，三个年表载荷均为正，其与格点6—8月降水和scPDSI正相关、与最高温度负相关。PC2解释32.5%，分离出海拔对比：Wolf强正载荷、Outlaw负载荷、Sunlight接近零。50年滚动相关显示样点间一致性随时间变化，涉及Sunlight的两对样点最一致（Sunlight–Outlaw平均r≈0.29，Wolf–Sunlight≈0.25），Wolf–Outlaw平均较弱（r≈0.13）且更常不显著或为负。移动窗口PCA显示共同生长信号的强度和结构随时间变化，PC1平均解释49.7%（范围38.0%–61.8%），PC2平均30.4%（范围21.5%–37.6%）；多数窗口三个年表PC1载荷为正，但Wolf在1368–1557 ce的持续窗口区间与主导成分解耦，最强负载荷出现在1418–1517 ce窗口，同期Sunlight和Outlaw在PC1上仍强正。

3.3 生长同步性与成对生长差异：同步性—MAD分析显示不同样点对间生长时间与幅度关系不同。Wolf–Outlaw的50年滚动同步性与平均绝对生长差异基本无关；Wolf–Sunlight呈显著负关系，Sunlight–Outlaw负关系更强。成对生长差异记录显示海拔间相对生长反复发生数十年尺度转换，968–2003 ce间三对差异均值接近零，但20年平滑显示持续偏离零。涉及Outlaw的对比显示最清晰的长期转变：13世纪部分时段和20世纪末至21世纪初，Wolf–Outlaw与Sunlight–Outlaw差异多为正；15世纪部分时段低海拔样点生长多于Wolf，表明相对生长优势在共同区间内反复转移。

3.4 对复合气候极端的生长响应：SEA和BEA分析聚焦968–2003 ce共同区间，干热年（重建8月T_max≥75百分位且NADA PDSI≤10百分位，n=41）和冷湿年（T_max≤25百分位且PDSI≥90百分位，n=44）。Wolf对干热年响应最弱，仅在滞后+1年显著生长下降，冷湿年在事件年及滞后+1年产生显著正偏离。Sunlight对两类复合极端均有显著响应，干热事件年生长显著降低并在滞后+1和+3年持续受抑，冷湿年从滞后?1至+2年产生显著正偏离。Outlaw在复合极端后表现出最清晰的多年的响应，干热年后在滞后+1、+2和+3年显著生长下降，冷湿事件从事件年至滞后+3年产生显著正偏离。BEA显示Wolf在滞后1年超过95%零包络但更长滞后未超过；Sunlight在滞后1–10年超过95%零包络并在更长滞后接近上包络；Outlaw在短至十年尺度滞后显示高观测聚集，但曲线未持续超过95%零包络。

讨论部分总结：讨论指出三条年表表明大黄石生态区径向生长强烈受水文气候变率组织，各点scPDSI正相关和当年夏季最高温度负响应表明水分限制并不局限于最低海拔种群，水分可利用性在整个采样生态幅度上限制生长，而热量在生长季放大了这一限制，这与径向生长整合形成层活动、碳分配和水分状态的生理预期一致；柔松虽常被描述为耐旱且生理可塑，但耐旱并不意味着对水分限制不敏感，其他北美西部柔松记录也报告了类似的 moisture-sensitive 生长。PCA结果识别出两个互补维度：PC1为共享水文气候模态，PC2为Wolf与Outlaw之间的海拔依赖对比，这种共同强迫与局地过滤的组合是山地树轮网络的反复特征；相邻样点对尤其Sunlight–Outlaw同步性与生长差异呈更强负关系，Sunlight接近零的PC2载荷和强气候相关提示其处于过渡位置。移动窗口PCA显示生长一致性结构非平稳，Wolf在1368–1557 ce窗口与主导模态解耦，与附近Beartooth高原冰缘体独立古气候证据一致，该记录中约1400–1550 ce间同位素偏负，尽管定年不确定性和100年窗口限制逐年归因，但对应关系支持Wolf在早期小冰期条件下周期性转向更强的冷季、积雪或生长季长度限制的解释，说明高海拔柔松并非低海拔树木的“少受干旱胁迫版本”，而可能占据独特环境生态位，这有助于解释未来增温的梯度不对称效应。复合极端分析显示干热与冷湿年留下海拔特异的生长印记：Wolf干热抑制有限并从冷湿异常中受益，Sunlight在BEA中与干热事件时间耦合最强，Outlaw在SEA中干热后多年遗留最清晰，表明对复合极端的脆弱性并非沿单一单调梯度表达。Outlaw干热后长达三年的生长下降与干性森林中报告的干旱遗留效应一致，可能反映水力修复不完全、非结构性碳水化合物耗竭、形成层活动减弱或冠层与根系恢复延迟，树轮数据无法直接分离这些机制，但滞后生长响应识别了这类生理携带效应在野外表达最强的位置。更广泛的含义是长期树轮记录可揭示环境生理如何在空间和时间上运作：区域气候变率与一致生长格局相关，而海拔决定主导限制是水分亏缺、热放大的干旱胁迫还是冷季限制；对柔松而言未来增温可能减弱林线附近某些热限制，同时增加低海拔的蒸发需求、干旱遗留和生长抑制，长期恢复力不仅取决于物种水平的耐旱性，也取决于塑造生长、恢复和种群间一致性的海拔特异生理限制。

结论翻译：本文表明，大黄石生态区柔松生长由一个强的、千年尺度的水文气候信号所组织，并持续被海拔依赖的限制所修饰。Wolf、Sunlight和Outlaw的生长变率反映主要由水分可利用性和夏季温度驱动的广泛区域一致性，而幅度、敏感性和持续性的差异沿海拔梯度出现。这些格局表明区域气候强迫使各海拔生长同步，但局地生理限制决定树木如何响应和恢复于气候胁迫。在长期背景下，柔松已在数百年气候变率中持续存在，但本文记录的海拔差异表明其对未来变化的脆弱性并不均匀。高海拔种群相对缓冲于短期干旱胁迫，而低海拔种群在干热期表现出强烈且持久的生长抑制。随着区域增温和干旱加剧，这些限制很可能增强，增加低海拔胁迫同时改变高海拔样点的限制平衡。这些结果表明森林对气候变化的响应随海拔系统性变化，并强调在评估长期森林恢复力时需同时考虑区域气候一致性和局地生态背景。