Hydraulic traits govern opposing range shifts of montane trees under warming来自 Hydraulic traits govern opposing range shifts of montane trees under warming | Nature Climate Change 基本信息题目水力学性状调控变暖背景下高山树种相反方向的分布区偏移Hydraulic traits govern opposing range shifts of montane trees under warming期刊名称Nature Climate Change发表时间2026年8月6日Volume 16作者信息张先亮Xianliang Zhang、刘鸿雁Hongyan Liu、Rubén D. Manzanedo、Justin T. Maxwell、Loïc DOrangeville、杜明超Mingchao Du、Tim RademacherDOI10.1038/s41558-026-02726-6核心观点基于半球尺度树木年轮数据45个树种、3,057个样点、121,743株树木、全球山地树种海拔分布区偏移观测102个树种以及11个水力学功能性状研究揭示了水力学性状水力效率与水力安全的权衡是调控山地树种在气候变暖和干旱下分布区偏移方向与速率的核心生理机制。气候敏感、水力效率高如较高的叶/枝条比导水率的树种倾向于快速向高海拔迁移以追踪适宜气候而耐受干旱、水力安全边际大如较负的栓塞水势 P50 和较低的膨压丧失点 TLP的树种则倾向于在低海拔存留或向低海拔扩张。️ 一、 引言部分简要总结 (Introduction Summary)现实与科学背景山地覆盖了全球 25% 的陆地面积具有显著的垂直带谱并孕育了极高的生物多样性。在气候变暖驱动下高山树种正经历加速的海拔分布区偏移。前人研究局限性与学术争议传统假说基于热生态位理论预测物种将统一向高海拔温度较低处迁移。然而实测表明树种的响应存在高度异质性虽然欧洲及北半球多数树线呈上升趋势但仍有约 33% 的树种表现出向低海拔偏移的逆向现象且物种迁移速率普遍滞后于气候等温线的移动速率。仅凭气候变暖速率等外部暴露度指标无法解释物种间迁移方向与速率的分化缺乏将器官尺度的生理机制与大尺度生物地理分布相联系的预测框架。核心科学问题本研究旨在回答树木叶片与枝条尺度的水力学性状水力安全与水力效率的权衡如何调控树木径向生长对温度与干旱的海拔敏感性水力学性状如何直接或间接决定树种海拔分布区偏移的方向与速率 二、 研究方法详细介绍与分析 (Methodology Framework)数据源与样本筛选山地范围与样点依据全球山地生物多样性评估GMBA v2.0划定的 8,616 个山脉边界筛选位于山地范围内的样点。树木年轮数据整合国际树木年轮数据库ITRDB与树木生态网络DendroEcological Network提取 1901–1990 年间全球 3,057 个样点、45 个山地树种10 个裸子植物35 个被子植物共 121,743 株树木的径向生长年轮数据。分布区偏移数据从 BioShifts 全球物种分布区重塑数据库中提取 102 个山地树种的海拔偏移速率m/yr。水力学性状从 TRY 植物性状数据库及文献中提取 45 个树种的 11 个水力学性状包含最大树高 Hmax、实际树高 Hact、枝条比导水率 Ks、叶比导水率 Kl、木质部栓塞 50% 导水率损失水势 P50、预晓水势 Ppd、正午水势 Pmd、膨压丧失点 TLP、叶面积/边材面积比 AsAl、导管/管胞直径 Vdia 以及边材密度 WD。气候数据采用 CRU v4.05 格点数据集分辨率 0.5°提取月平均/最高/最低气温、降水、日较差 DTR 以及自校准帕尔默干旱严重程度指数scPDSI。数据清洗与时间动态评估使用 dpR 包中的 Friedman 去趋势方法α 5构建样点树轮宽度标准化年表。设置 40 年滑动窗口1901–1940 至 1951–1990以 5 年为步长评估生长-气候关系的海拔依赖性在 20 世纪的动态稳定性。为消除不同山脉间宏观气候如大气干燥度、土壤类型的混淆效应在单一山脉内部开展树种特异性的海拔梯度重新评估。统计模型与归因模型工作流海拔气候敏感性计算利用 Bootstrap 重采样1,000 次拟合 1951–1990 年间各树种径向生长对气候ScPDSI、降水、温度敏感性随海拔变化的斜率。Deming 回归模型考虑到水力性状与迁移速率均存在测量误差采用 Deming 回归Deming regression评估水力学性状与物种分布区迁移速率的关系避免普通最小二乘法的无偏估计偏误。主成分分析与结构方程模型采用主成分分析PCA将 11 个水力学性状降维为 PC1水力效率轴Kl, Ks, Vdia解释 47.4% 方差和 PC2水力安全轴P50, TLP, AsAl, Ppd解释 25.8% 方差。构建结构方程模型SEMSPSS AMOS解析“水力性状 - 海拔气候敏感性 - 迁移速率”的直接与间接级联路径。 三、 主要研究结果与充实数据分析 (Main Results Quantified Data)生长-气候关系随海拔的响应分化与时间解耦在 1951–1990 年间不同树种的干旱与温度敏感性随海拔呈现显著分化。18% 的树种干旱敏感性随海拔显著增加13% 显著降低13% 的树种温度敏感性随海拔显著降低。约三分之一33%的树种在 20 世纪后半叶表现出生长-干旱敏感性海拔依赖关系的时间解耦或反转例如 Quercus alba 从正相关转为负相关。水力学性状调控海拔气候敏感性水力效率性状Kl, Ks, Vdia与干旱敏感性的海拔斜率呈显著负相关但与温度敏感性的海拔斜率呈显著正相关如 Kl 与温度敏感性斜率 R² 0.33, P 0.05Ks 与温度敏感性斜率 R² 0.27, P 0.05。水力安全性状P50, TLP, Ppd, Pmd则与干旱敏感性的海拔斜率呈正相关。水力学性状、生长-气候敏感性与物种海拔偏移的概念机制框架详见 Figure 145 个树种在 1951–1990 年间生长-气候相关系数随海拔变化的斜率分布详见 Figure 2关键水力学性状对海拔气候敏感性的调控关系详见 Figure 3。图片名Figure 1图例说明Figure 1 | 水力学性状调控变暖和干旱下高山树种海拔分布区偏移的机制框架图。展示了水力学性状水力安全与效率权衡直接影响分布区偏移方向与速率以及通过调控树木径向生长的海拔气候敏感性间接驱动分布区重塑的双重路径。图片名Figure 2图例说明Figure 2 | 1951–1990年间 45 个高山树种树木径向生长与气候因子日较差 DTR、scPDSI 干旱指数、降水量、温度相关系数随海拔变化斜率的 Bootstrap 估计分布。橙色代表被子植物绿色代表裸子植物。图片名Figure 3图例说明Figure 3 | 水力学性状对生长-气候海拔敏感性的调控作用。a 各水力性状与气候敏感性斜率的交叉相关矩阵b-e 显示叶比导水率 Kl (b)、枝条比导水率 Ks (c)、栓塞水势 P50 (d) 以及膨压丧失点 TLP (e) 对温度、降水量和干旱指数海拔敏感性的影响。水力学性状直接决定分布区偏移方向与速率山地树种的海拔迁移速率存在极大差异介于向上快速迁移5 m/yr如 Acer negundo, Populus alba与向下扩张3 m/yr如 Quercus coccifera, Acacia dealbata之间。样点气候变暖速率与物种分布区偏移速率之间无显著相关性R² 接近 0, P 0.05表明单纯的气候暴露度无法解释迁移分化。Deming 回归表明叶比导水率 KlR² 0.25, P 0.003, 直线斜率 1.00、膨压丧失点 TLPR² 0.14, P 0.05, 直线斜率 1.37以及栓塞水势 P50R² 0.12, P 0.05, 直线斜率 0.57与物种海拔迁移速率呈显著正相关。水力效率高高 Kl和水力安全边际小较不负的 P50 和较高 TLP的树种倾向于快速向高海拔迁移而水力安全策略保守较负的 P50 和较低 TLP的树种倾向于向低海拔扩展或在低海拔存留。各水力学性状与物种海拔迁移速率的 Deming 回归关系详见 Figure 4。图片名Figure 4图例说明Figure 4 | 水力学功能性状对树种海拔迁移速率与方向的驱动作用。a 各水力性状与迁移速率的相关系数b-d 采用 Deming 回归拟合的叶比导水率 Kl (b)、膨压丧失点 TLP (c) 及栓塞水势 P50 (d) 与海拔迁移速率的显著正相关关系。SEM 级联路径与机制解析SEM 模型解释了物种迁移速率 37.2% 的方差。PCA 第一主成分 PC1水力效率通过强烈增强海拔温度敏感性路径系数 0.716, P 0.001并抑制海拔干旱敏感性路径系数 -0.291, P 0.05间接驱动树种向高海拔快速迁移温度敏感性对迁移速率路径系数 0.785, P 0.01。PC2水力安全则主要通过调控降水敏感性路径系数 -0.216, P 0.05和干旱敏感性路径系数 -0.283, P 0.05影响物种在低海拔地区的存留与重塑。水力学性状、PCA 主成分与生长敏感性及迁移速率的结构方程模型 SEM 路径详见 Figure 5。图片名Figure 5图例说明Figure 5 | 结构方程模型SEM解析水力学性状调控树种海拔迁移的直接与间接路径。a 关键水力性状的路径分析b 水力性状的主成分分析PC1 代表水力效率PC2 代表水力安全c 基于 PCA 主成分的 SEM 路径图模型整体 R² 37.2%蓝色为正路径红色为负路径箭头宽度对应标准化路径系数。 四、 讨论部分的不足与未来展望 (Discussion Limitations Perspectives)研究不足与不确定性来源非气候与生物间相互作用的混淆受限于全球尺度数据未完全偶联种间竞争、天敌/病虫害、土地利用历史遗留效应以及氮沉降梯度对树种迁移的调制作用。高海拔生境面积收缩限制Spatial Squeeze山地地形随着海拔升高呈非线性面积收缩2,000 m 以上区域仅占全球山地总面积的不到 5%向高海拔迁移的水力敏感型树种面临剧烈的生境压缩与高强度的种间竞争风险。数据时空跨度与末端滞后树轮与 BioShifts 数据主要覆盖 1901–1990 年1990 年以后的最新加速变暖响应未完全包含植被分布区尾端trailing edge存在迁移惯性Migration inertia与滞后效应。未来研究展望建立分海拔梯度的精细化保护策略依据树种水力学性状划分保护优先级对高水力效率、高海拔收缩风险的树种建立高海拔气候庇护所Climate refugia对水力安全型树种关注低海拔耐受机制。结合高分辨率微气候与生理表达结合微地形Microtopography微气候观测与树木水力学性状表观可塑性Phenotypic plasticity提高对复杂地形下树种分布区重塑的预测能力。整合动态植被与生物地理模型将水力学性状驱动的生长-干旱/温度敏感性机制无缝并轨至下一代地球系统模型ESM和动态全局植被模型DGVM中。