基于滞后模型分析树木生理阈值响应干旱的转移机制

《Plant, Cell & Environment》：Employing a Hysteresis Approach to Analyze Shifts in Tree Physiological Thresholds in Response to Drought

针对气候变化下树木干旱死亡风险加剧、现有生理模型忽视内部储水等关键环节的问题，本研究创新性地采用滞后（hysteresis）建模方法，分析了两个对比生态系统（云雾林与美国温带森林）中液流（Js）、茎干直径波动（SDF）和茎干体积含水量波动（SVWCF）与饱和水汽压差（VPD）的日循环数据。研究量化了滞后环参数（质心、面积、滞后、旋转角），揭示了树木在胁迫下表现出的“茎干水分提取”与保守的“避旱”策略，证实滞后衍生参数是诊断干旱胁迫、预测树木水力恢复力及用水策略的有效早期指标。

　　

在全球气候变暖的背景下，极端干旱事件正变得越来越频繁，这对森林生态系统构成了前所未有的威胁。树木的生存与否，往往取决于其应对水分胁迫的精细生理机制。然而，长期以来，科学家们用于预测树木响应的生理模型存在一个明显的短板——它们常常忽略了树木内部储水这一关键缓冲机制。这就好比我们只看汽车的油耗（蒸腾作用），却忽略了油箱里到底还剩多少油（内部储水）。当土壤水分枯竭时，树木并非坐以待毙，它们会动用储存在茎干木质部中的“救命水”，或者采取保守策略关闭气孔以减少损失。但这种从“吸收”到“解吸”的动态转换过程极其复杂，传统的线性模型难以捕捉这种非线性的、带有时间延迟的生理响应。

为了攻克这一难题，揭示树木在干旱下的真实应对策略，一项发表在《Plant, Cell & Environment》上的研究另辟蹊径。研究人员采用了一种名为“滞后（hysteresis）”的建模方法，试图通过分析树木生理指标的日循环轨迹，来解码它们在水分胁迫下的“行为密码”。这项研究在两个截然不同的生态系统中展开：一个是物种丰富、湿度极高的秘鲁云雾林（cloud forest），另一个是物种组成相对简单的美国温带森林。通过在这两个地点设置穿透雨减少（throughfall reduction）实验，模拟自然发生的干旱条件，研究团队收集了包括液流密度（Js）、茎干直径波动（SDF）、茎干体积含水量波动（SVWCF）以及大气饱和水汽压差（VPD）在内的高频小时数据。利用广义超越方程，他们量化了五种日滞后环，并从中提取出质心（centroids）、面积（area）、滞后（lag）和旋转角（angle of rotation）等关键参数。最终，研究不仅证实了滞后参数作为诊断工具的有效性，还成功识别出了树木在面对干旱时采取的两种截然不同的用水策略：一种是依赖内部储水的“提取策略”，另一种是维持正水平衡的“保守避旱策略”。

为了达成上述目标，作者主要运用了以下几项关键技术方法：首先，在两个对比鲜明的生态系统（秘鲁云雾林和美国温带森林）设立了穿透雨排除实验样地，分别监测了 Weinmannia bangii 和 Quercus spp. 等树种。其次，利用热比率法（heat-ratio method）构建液流传感器测量液流密度（Js_hour），使用点测微计（point dendrometer）测量茎干直径波动（SDF_hour），并利用经校准的探针测量茎干体积含水量波动（SVWCF_hour）。此外，结合气象数据计算饱和水汽压差（VPD）。最后，核心分析采用了 R 语言中的“hysteresis”包，基于 Lapshin 的广义超越方程拟合椭圆模型，从每小时聚合的数据中提取滞后环参数，并使用线性混合效应模型比较干旱与控制处理间的差异。

3.1 观测到的滞后环（Observed Hysteresis Loops）

通过对 Js_hour、SDF_hour和 SVWCF_hour的日周期分析，研究发现这些变量与 VPD 之间呈现出显著的非线性滞后环，清晰地描述了树木的吸水（absorption）与解吸（desorption）动态。具体来说，解吸阶段始于黎明后 VPD 升高导致 Js 增加，消耗茎干储水导致茎干收缩（SDF 负值）；而吸收阶段则发生在午后至夜间，随着 VPD 降低，根系的吸水作用使得茎干水分得到补充，茎干膨胀（SDF 正值）。尽管干旱处理和对照处理的滞后模式不同，但在不同物种间，滞后环的循环方向（顺时针或逆时针）保持一致。

3.2 滞后环建模与阈值定义（Modelling Hysteresis Loops and Thresholds Definition）

研究使用的超越方程能够准确拟合野外实测数据，观测值与拟合值的 R²在 0.60 到 0.83 之间（p < 0.01）。虽然经验滞后环并不总是完美的椭圆（例如 SDF 环有时呈“鱼形”），但模型依然成功捕捉到了关键的生理阈值（x₁, x₂, y₁, y₂），这些阈值定义了从吸水到解吸以及反之的转换点。即便模型在某些 Js:VPD 环中预测出了负的液流值，但其揭示的生理阈值的生物学意义依然成立。

研究结果（续）

通过对比控制组和干旱组的椭圆衍生系数（表 4 和表 5），研究量化了干旱引起的具体变化。例如，在 Js:VPD 环中，干旱处理显著改变了质心（c_x, c_y）和饱和度点（b_x, b_y）的数值。更重要的是，研究发现干旱导致了滞后环几何形状的改变，这直接反映了树木生理状态的转变。

结论与讨论

这项研究的核心贡献在于建立了一套基于滞后分析的量化框架，用于解析树木生理阈值对干旱的响应转移。研究结果表明，滞后衍生参数（如质心、面积、滞后时间和旋转角）不仅仅是数学拟合的产物，它们是具有明确生理学意义的诊断指标。

首先，研究明确了树木在干旱胁迫下表现出的两种泾渭分明的用水策略。一种是“茎干水分提取策略”（stem water-extractive strategy），这类树木在胁迫下表现出强烈的茎干收缩和水分耗竭，意味着它们高度依赖内部储存的水来维持蒸腾，这是一种以牺牲储水为代价换取短期生存的激进策略。另一种则是“保守的避旱策略”（conservative drought-avoidance strategy），这类树木努力维持正的茎干水分平衡，通过更早或更严格的气孔调节来减少水分损失。这两种策略在滞后环的形态上表现为截然不同的特征，例如不同的滞后面积和旋转角。

其次，该研究强调了不能仅凭滞后“面积”来判断树木胁迫状态。以往的研究多关注面积的变化，但本文指出，生理阈值的突变（如气孔关闭或木质部空穴化）会导致非线性的生理响应，单纯的面积变化可能掩盖了背后的不同机制。通过引入质心、旋转角以及 x₁, x₂, y₁, y₂等阈值参数，研究提供了一个多维度的视角。例如，滞后（lag）参数反映了水力电容（hydraulic capacitance）的大小，即树木利用内部储水缓冲外界大气需求（VPD）的时间延迟；而 coercivity（矫顽力）则衡量了触发生理响应变化所需的驱动力强度。

最后，这项研究将植物生理学中的水分动力学过程与人体肺力学的滞后现象进行了巧妙的类比（如图 1 所示）。正如肺部的吸气-呼气循环依赖于粘弹性特性一样，树木的吸水-解吸循环也依赖于其内部水力结构的物理属性。这种跨学科的类比不仅加深了我们对植物水力系统的理解，也为开发下一代更精准的生理模型提供了理论基础。通过整合液流（Js）、茎干动态（SDF）和茎干含水量（SVWCF）的多源数据，该研究为预测气候变化背景下树木的存活率和水资源利用策略提供了强有力的新工具，具有重要的生态学意义和全球变化应用前景。