地中海松林干旱诱导衰退的多因素诊断：树轮、碳氧稳定同位素与养分失衡的整合研究

《Journal of Ecology》：Multifactorial diagnosis of drought-induced decline in Mediterranean conifer forests

　　

干旱诱导的树木死亡是一个复杂过程，涉及树木水分、碳和养分关系之间相互依赖的变化，同时伴随生长和性状变异的调整。在干旱加剧背景下预测衰退，需要对其潜在生理过程获得机制性认识。研究人员应用树木个体水平的多指标方法，以识别位于西班牙东北部地中海 Pinus sylvestris、Pinus pinaster 和 Pinus halepensis 森林死亡率高发区中与活力下降相关的因素。整合树木年轮指标（tree-ring metrics）、叶部碳和氧稳定同位素（δ¹³C、?¹⁸O）、元素组成（elemental composition）和土壤性质（soil properties），以评估干旱对生长（DGIi）、内在水分利用效率（iWUE）、由?¹⁸O推断的气孔导度（G_?¹⁸O）和树木养分状况的影响。衰退树木的 DGIi 始终高于非衰退树木。在 P. sylvestris 和 P. pinaster 中，DGIi 升高与 G_?¹⁸O 降低和 iWUE 升高相关，表明衰退树木遭受更严重的干旱胁迫。衰退的 P. sylvestris 和 P. halepensis 还表现出显著养分失衡，尤其是 P、K、Zn、Fe、S 和 Ca，且所有物种均显示 Fe:Cu 比值紊乱。物种间活力差异还与针叶大小和形态、树高和土壤 pH 变异相关，其中衰退的 P. sylvestris 和 P. halepensis 在单位面积叶质量（LMA）上表现出相反调整。综合：这些结果突出表明，水力功能损害、养分失衡、叶部衰老和不利土壤条件之间的相互作用，放大了地中海松林研究种群中的气候胁迫。此类潜在反馈可能加速树冠脱叶，并最终在干旱加剧下促进树木死亡。

研究背景与问题：复合高温—干旱事件引发的森林衰退和树木死亡已在全球森林生物群系中被记录。大气蒸发需求上升和干旱加剧威胁森林持续，因此需要改进森林衰退和树木死亡预测，并要求理解树木碳—水平衡机制。水力失败和碳饥饿被提出为解释干旱死亡的两个主要机制。由于环境、生物和树木内在特征交互，机制认识仍不足。采用同时考虑树木水分、碳和养分状态的多指标方法可推进理解。树木在干旱下通常降低生长，衰退或死亡树常死亡前生长下降或停滞。更脆弱树木和种群常表现更强气候—生长耦合，可作为早期预警信号。树木年轮提供年分辨记录。形态性状和树大小揭示获取型与保守型资源利用策略权衡，但种内模式不清。δ¹³C 和 δ¹⁸O 测量整合气孔调节、光合速率和水分利用效率。水分利用与养分状态紧密相关：干旱诱导气孔关闭降低蒸腾驱动的土壤养分吸收；表层土干燥使根系向深层吸水，而水与养分沿土壤/基岩剖面常垂直解耦。干旱还影响土壤养分可用性、离子移动、微生物活性、矿化及细根和菌根功能。养分对光合、水力功能和氧化还原代谢重要，从而支撑抗旱性和恢复力。地中海针叶林受气候变化严重影响，养分缺乏风险高。西班牙东北部松林死亡率高发区是研究重点。研究目标为：表征各物种健康与衰退树形态、干旱对生长影响和叶同位素组成；辨别干旱诱导衰退中叶养分变化及其与水分利用策略和土壤性质关系；确定与树冠落叶相关的结构、生理、营养和土壤因素。假设：衰退树更严重干旱胁迫；衰退与P、K及微量元素失衡相关；解释活力变异的因素依赖物种特征和立地条件。

研究开展与主要结论意义：研究人员在西班牙东北部沿干旱梯度选择四个死亡率高发天然林，包括两个 Pinus sylvestris 林、一个 Pinus pinaster 林和一个 Pinus halepensis 林。共选90株优势/共优树，每物种30株。按树冠落叶分非衰退（ND，<40%）和衰退（DD，>60%）。采集叶、木芯和土壤，测定形态、树轮、叶δ¹³C、δ¹⁸O和元素、土壤性质。结论：衰退树DGIi更高；P. sylvestris和P. pinaster衰退树G_?¹⁸O降低、iWUE升高；P. sylvestris和P. halepensis养分失衡，所有种Fe:Cu紊乱；物种活力差异还与针叶大小形态、树高和土壤pH相关。综合表明水力损害、养分失衡、叶衰老和不利土壤条件相互作用放大气候胁迫，潜在反馈可加速树冠脱叶并导致死亡。意义是提供多指标机制诊断框架。

主要关键技术方法（不超过250字）：在西班牙东北部四个死亡率高发天然林（两个 P. sylvestris、一个 P. pinaster、一个 P. halepensis）开展，共90株优势/共优树，每物种30株。根据树冠落叶分ND和DD。叶和P. sylvestris木芯于2023年夏采集，P. pinaster和P. halepensis年轮样于2020年秋采集，并采集0–15 cm土样。树木年轮交叉定年计算年轮宽度指数（RWI）、断面积增量（BAI）、干旱生长影响指数（DGIi）和标准化降水蒸散指数（SPEI）。叶片测δ¹³C和δ¹⁸O，推算iWUE和G_?¹⁸O；用电感耦合等离子体发射光谱（ICP-OES）测叶元素并计算N:P、N:K、K:Ca、Fe:Cu；测土壤质地、pH、C/N和有效养分。统计用线性模型、主成分分析（PCA）、置换多元方差分析（PERMANOVA）、偏最小二乘回归（PLS）和随机森林。

研究结果：
3.1 形态性状与生长对干旱的响应：通过树冠落叶调查、树轮和SPEI关系。P. sylvestris死亡率高发区ND落叶约33%、DD约58%；P. pinaster ND约19%、DD约79%；P. halepensis ND约15%、DD约73%。年龄无显著差异。P. sylvestris ND胸径（DBH）和树高小于DD；P. pinaster和P. halepensis ND更粗更高。DD P. pinaster和P. halepensis针叶更小；DD P. sylvestris LMA更大，P. halepensis LMA更小。自1960年BAI下降，P. pinaster和P. halepensis DD下降更陡；最近25年ND的BAI更高。生长对SPEI响应：P. sylvestris时间尺度最长8–10个月；P. pinaster 5个月，ND/DD无差异；P. halepensis 2个月，DD更强。RWI-SPEI相关在P. sylvestris和P. halepensis的DD更强。DGIi在P. halepensis最高。结论：衰退树生长下降更强、干旱影响更大，但种间响应时间尺度不同。

3.2 叶部同位素组成与养分含量：通过叶同位素、元素、PCA和土壤相关。iWUE在P. sylvestris最低，P. pinaster最高；ND P. sylvestris和P. pinaster低于DD。?¹⁸O在P. sylvestris最高，其次P. pinaster、P. halepensis；P. sylvestris和P. pinaster DD的?¹⁸O高于ND，指示DD的G_?¹⁸O更低；P. halepensis无显著差异。iWUE与?¹⁸O在P. sylvestris和P. pinaster正相关。元素：P. sylvestris DD的N、Ca、S、Fe及N:P、Fe:Cu更高，Zn更低；P. pinaster仅Fe:Cu在DD更高；P. halepensis DD的N:K和Fe:Cu更高，K、Cu、Zn更低。PCA显示P. sylvestris和P. halepensis叶多元素组成在活力级间显著不同，P. pinaster不显著。土壤：P. pinaster较高pH与P、N、Zn、Ni较高相关，较高有效P与Cu、Mg、Fe较低相关；P. halepensis有效P与多数叶养分正相关，pH相反。结论：P. sylvestris和P. halepensis存在养分失衡，所有种Fe:Cu改变，同位素和养分与水分策略及土壤性质相关。

3.3 变量协变及其对树冠落叶的影响：通过PLS和随机森林。PLS解释校准方差R²=0.51–0.70，交叉验证Q²=0.25–0.34，RMSEP=14%–18%。唯一三种松树一致显著的是BAI与落叶负相关。P. sylvestris中LMA和iWUE与落叶正相关，Fe及N:P、N:K、Fe:Cu也正相关。P. pinaster和P. halepensis中树高与落叶负相关；P. pinaster中iWUE、?¹⁸O和叶S与落叶正相关，土壤有效P正相关。P. halepensis中针叶大小和LMA与落叶负相关，Ca和Fe:Cu正相关，Cu和Zn负相关，土壤pH正相关，有效P负相关。跨物种随机森林支持这些关联。结论：径向生长是跨种最一致指标，其他因素具种特异性。

讨论部分总结：4.1 生长和同位素信号与冠层状况相关。长期径向生长下降和低生长率是干旱诱导死亡指标。P. pinaster和P. halepensis ND与DD生长轨迹在二三十年前分化，可能由1994–1995和1986干旱触发；P. sylvestris缺乏分化可能反映极端事件导致快速水力失败。DD更负生长趋势、更强干旱—生长耦合和更高DGIi。P. sylvestris和P. pinaster DD的δ¹³C推断iWUE更高，iWUE与?¹⁸O正相关，说明DD气孔调节更严格。2022年极端复合干旱热浪可能放大这些差异。P. halepensis同位素与落叶无关，可能全林分慢性水分胁迫，形态和养分调节。4.2 树木营养受气孔调节对树水流影响和土壤性质驱动。P. sylvestris和P. halepensis元素差异明显，P. pinaster弱，说明后者死亡主要受水力/碳限制。所有种DD的Fe:Cu改变，P. sylvestris还积累Fe、Ca、S，P. pinaster S高，P. halepensis Cu低。Fe:Cu可能为营养胁迫新指标。机制包括蒸腾驱动养分质量流降低、韧皮部运输和形成层活动下降、低移动元素封存及土壤pH影响。N:P和N:K失衡与P. sylvestris和P. halepensis落叶相关；P、K不足损害水力光合，Zn低影响膜稳定和抗旱。4.3 多因素途径诊断衰退。径向生长最一致。P. sylvestris和P. pinaster落叶与iWUE和?¹⁸O正相关；P. sylvestris生长与二者负相关，符合气孔关闭限碳、非结构性碳水化合物（NSC）消耗和叶衰老脱落。P. sylvestris中LMA和iWUE与落叶正相关，显示硬叶和严格气孔；P. halepensis相反，高LMA与低落叶相关，高iWUE与弱干旱影响相关，健康树K、Zn、Cu更高。P. pinaster中iWUE与P、K正相关。树高在P. pinaster和P. halepensis与落叶负相关。土壤pH在P. halepensis较高与较差冠层相关；P. pinaster酸性土上较高pH与N、P、Zn较高相关，有效P与落叶正相关可能反映衰退树凋落物富集表土。方法上树轮提供长期，叶为单年快照但机制可解释。

研究结论翻译：研究结果突出表明，地中海松林干旱诱导衰退源于相互依赖的水力、碳和营养约束。在所研究的松树种群中，衰退树木一致表现出径向生长受损和更强的干旱生长影响。衰退的 P. sylvestris 和 P. pinaster 中，iWUE升高和G_?¹⁸O降低的同位素证据揭示了慢性干旱胁迫及相关的碳限制。与养分谱一起评估时，这些同位素特征与一系列相关胁迫响应一致，其中持续气孔关闭和土壤驱动限制与养分失衡（Fe、Ca、S、Zn及N:P:K化学计量）共同发生，进而促进叶片衰老加速、水力功能损害并最终导致树冠落叶。在 P. halepensis 林分中，养分失衡以及叶形态和生理变化（如LMA和iWUE）似乎主导干旱胁迫响应。所有物种衰退树木均显示Fe:Cu比值改变，提示其可能是干旱相关营养胁迫的新指标。P. pinaster 中相对较弱的营养变化表明，该种群衰退主要由碳—水动态和/或立地条件驱动。总体而言，结果揭示一种潜在级联序列：干旱驱动的树冠损失和死亡源于气孔调节、碳耗竭、营养损害、氧化胁迫和土壤限制之间的反馈。结果还突出，只有多指标方法才能捕捉这些机制如何共同发生和相互作用。该整合视角提供了机制性指标，有助于检测树木死亡背后的生态生理模式，并为干旱加剧下的适应性森林管理提供信息。未来研究应纳入更广泛的环境和土壤条件，澄清关键微量元素作用，并将养分失衡与蒸腾动态和根系深度等水分利用性状联系起来。