风驱动集体动力学主导人工苏格兰松林分的树木响应

《Agricultural and Forest Meteorology》：Wind-driven collective dynamics dominate tree response in a planted Scots pine forest stand

　　

摘要：森林中由风引起的树木运动反映了大气强迫与树木特定机械特性之间的相互作用，然而，在野外条件下，林分尺度的响应究竟由相干风强迫主导还是由个体树木特征主导仍不清楚。在本研究中，研究人员利用奇异值分解（SVD）和最大协方差分析（MCA）对人工苏格兰松林分中28棵相邻树木的风致响应进行了分析，以识别耦合的风-树分量，并将树木响应分解为集体份额和个体份额。结果表明，树木响应由一个单一的风驱动耦合分量主导，所有树木的响应在很大程度上都投影到该分量上。这使得能够一致地分解为集体响应（代表共享的林分尺度动力学）和个体贡献（捕捉树木特定的偏差）。集体参与与主导风驱动分量的对齐程度系统性地相关，而个体贡献表现出更大的变异性，并且不依赖于内在动态特性或局部结构因素（如邻域配置）。偏相关分析表明，除了参与主导集体响应之外，基本摆动频率和拥挤指数均不能解释风-树耦合的额外变异。这些发现表明，在结构均匀的林分中，树木响应主要由林分尺度的风荷载主导，表现为一个单一的风驱动耦合模式，该模式解释了对齐后树木响应场中约82%的方差。个体树木特性主要调节偏离该主导集体响应的偏差。这为人工苏格兰松林分中的风-树相互作用提供了一个降阶、基于过程的视角，并突显了自然风条件下集体动力学的核心作用。

风扰动是森林动态的重要自然驱动因素，极端风事件通过拔根和茎干断裂导致大规模树木死亡，重塑森林结构、改变物种组成，并影响碳储存和养分循环等关键生态过程。在人工林中，风害是一个持续风险因素，林分易损性取决于环境强迫和森林结构属性。现有研究通常从个体树木力学和林分集体行为两个层面解释风害，但动态联系——即风荷载与多棵树集体响应之间的耦合机制——仍缺乏深入研究。早期实验和建模表明，林分结构强烈影响个体树木的风荷载，例如间距增大导致树冠更大、风荷载更高，但抵抗倾覆能力未成比例增加，而密度较高的林分中树木受益于邻居的机械支撑和遮蔽。然而，对于林分尺度响应是否主要受相干风强迫或个体树木特性支配，在野外条件下仍不明确。本研究基于对人工苏格兰松（Pinus sylvestris L.）林分中28棵相邻树木风致运动的先前分析，从林分尺度动力学系统视角重新分析数据，旨在确定树木响应是否主要由共同的风驱动响应或个体树木行为解释。研究人员通过量化时间滞后、使用奇异值分解（SVD）和最大协方差分析（MCA）分析耦合风-树场，并将对齐后的树木响应场分解为集体和个体份额，以提供降阶、基于过程的视角。论文发表在《Agricultural and Forest Meteorology》。

本研究所用关键技术方法包括：样本来自德国弗赖堡大学运营的哈特海姆森林研究站（47°56′04″N, 7°36′02″E，海拔201 m），选取28棵相邻人工苏格兰松，测量其茎倾斜（10 Hz采样率）和风速（四个超声风速仪，安装高度z/h=1.2）。主要方法包括：1) 奇异值分解（SVD）识别风荷载和树木响应中的主导时空模式；2) 最大协方差分析（MCA）量化风-树场耦合，提取主导耦合分量；3) 通过交叉相关分析确定各树木相对于主导风荷载信号的响应时间滞后（τ_i），并利用线性回归拟合空间滞后结构，计算视速度（v_prog）；4) 将时间对齐后的树木响应矩阵分解为集体份额（由第一SVD分量重建）和个体份额（残差）；5) 傅里叶分析计算基本摆动频率（f₀）；6) 基于高斯核函数定义拥挤指数（CI）量化局部几何嵌入。

研究结果包括以下六个子部分：

**3.1 风-树相互作用的低维结构**：通过SVD分析，风荷载第一分量解释65.3%的方差，而树木响应第一分量解释82.4%的方差，表明树木响应高度集中于一个主导分量，反映加载驱动的同步性，而非共振型动力学。这为分离集体和个体贡献提供了定量基础。

**3.2 风荷载的空间滞后结构和推进**：通过交叉相关分析，各树木的残差滞后（τ_i）与沿最佳拟合推进方向的位置呈显著线性关系（r²=0.75），对应视速度5.4 m/s。推进方向与独立测量的平均风向仅差13.1°，且小时尺度上中位角度差10.5°，95.6%的估计值差异小于30°。视速度与小时平均风速强相关（r²=0.79），表明滞后结构反映了动态风场推进。

**3.3 对齐和时间一致性的验证**：时间序列显示，对齐后树木响应矩阵的第一主分量（PC1）在整个观测期间始终主导，振幅与风速同步变化，第二分量（PC2）幅值较小且交替正负，可能代表恢复运动。这证实对齐过程成功隔离了相干动力学。

**3.4 分解为集体和个体分量**：集体份额中位数为0.82，范围约0.70–0.90，个体份额相应较小，但空间分布不均匀，表明集体响应主导但存在局部变异。

**3.5 集体动力学与风对齐的耦合**：集体份额与对齐度（r_aligned）呈正相关，而个体份额无正相关甚至微弱负相关。偏相关分析显示，在控制集体份额后，拥挤指数（CI）和基本摆动频率（f₀）均不能显著解释风-树耦合的额外变异（r=0.03, p=0.86; r=0.24, p=0.23）。这表明集体响应直接反映风驱动耦合，个体特性主要调节偏离。

**3.6 内在树木动力学对响应分区的有限影响**：集体份额在整个f₀范围内保持较高，无系统性关系，个体份额也无显著依赖。树高与集体份额有弱正相关（r=0.33, p=0.082），但林分均匀性限制了尺寸差异的影响。

讨论部分总结：研究结果表明，在结构均匀的人工苏格兰松林分中，风致树木响应主要表现为风驱动的低维过程，集体响应占主导（约82%方差），个体特性仅调节偏差。空间滞后结构证实集体响应相干推进，方向与风速关系支持风场推进解释。这提示风-树动力学不能仅通过个体树木特性理解，而需考虑林分尺度组织。对森林管理而言，风害减轻策略应关注林分尺度风荷载暴露，而非仅优化个体稳定性，例如间伐会同时改变结构、风暴露和邻居相互作用，其效果取决于时间和空间模式。未来研究需探究管理引起的结构异质性变化如何影响集体与个体响应的平衡及风害风险。

结论部分翻译如下：本研究证明，在相对均匀的人工苏格兰松林分中，风致树木运动可描述为以风驱动为主的低维过程。树木响应一致地收敛于一个主导集体分量，表明相干风荷载在林分尺度上主导动力学，而个体树木特性主要调节偏离该共同响应的偏差。空间滞后模式进一步表明集体响应相干线内推进，推进方向与风向一致，速度与风速相关，证实剩余小时间偏移与动态风荷载结构通过林分相关联。风-树耦合与局部林分结构、内在动态特性或基本摆动频率无独立关系，进一步支持在调查条件下树木响应主要由宽带相干风荷载主导，而非频率选择性放大。该结论适用于调查的均匀林分和中度风条件，但共振效应在异质林分或强风下不可排除。这些发现表明，风致树木运动不能通过孤立考虑个体树木完全理解，而需反映其在林分尺度风驱动响应中的参与，同时保留树木特定偏差。从实践角度，森林管理策略应兼顾林分尺度风荷载暴露，而非仅优化个体稳定性。间伐作为关键例子，同时改变结构、风暴露和邻居相互作用，其效果取决于强度、空间模式和时间。未来工作需探究管理引起的结构异质性变化及后续适应如何影响集体与个体响应平衡及风害风险。