基于稳定碳同位素的波斯栎干旱敏感性分析——以伊朗西部扎格罗斯山脉为例

《European Journal of Forest Research》：Drought sensitivity analysis of Persian oak based on stable carbon isotopes from the Zagros Mountains, western Iran

　　

摘要：耐旱的本土树种正日益受到气候变化加剧影响的威胁。在这项研究中，研究人员首次对波斯栎（Quercus brantii Lindl.）逐年分离的树木年轮开展了树木年轮稳定碳同位素和内在水分利用效率（iWUE）研究，时间跨度为1953—2022年。研究地点位于扎格罗斯山脉的三个不同海拔高度。所有样地的δ¹³C年表在1964年、1984年和2007—2008年的主要干旱事件期间均表现出较高（负值较小）的δ¹³C值，随后其各自的iWUE值也相应增加。高海拔样地的iWUE值最高。中海拔和高海拔样地的iWUE在2008年至2022年期间显著增加，表明树木对日益加剧的水分限制产生了共同的生理响应。与气候变量的相关性分析显示，在中海拔和高海拔样地，冬春季降水量与δ¹³C呈显著负相关。所有三个样地均与帕尔默干旱严重程度指数（PDSI）和标准化降水蒸散指数（SPEI）在1至12个月的累积时段内呈负相关，证实了波斯栎对上个生长季和当前生长季的水文气候条件具有高度敏感性，其中中海拔和高海拔样地的相关性尤为强烈。这些模式凸显了干旱在调节气孔导度和碳吸收中的主导作用，为扎格罗斯山脉观测到的栎树枯梢现象提供了生理学证据。总体而言，波斯栎的δ¹³C变化主要受水文气候因子（降水、相对湿度和饱和水汽压差）驱动，且这些因子的影响时间和强度在研究样地之间随空间和时间而变化。研究结果揭示，波斯栎沿海拔梯度表现出不同的适应性生理策略，从较为干旱环境（低海拔和高海拔）中的保守型水分利用行为，转变为中等湿度条件下（中海拔）更为灵活的气体交换调节。

论文解读：扎格罗斯山脉波斯栎干旱敏感性的稳定碳同位素证据

研究背景与科学问题：耐旱的本土树种在半干旱生态系统中对水资源调节和全球碳循环具有重要作用，但对气候变化的影响尤为脆弱。伊朗西部的扎格罗斯山脉拥有西亚面积最大的栎树森林，被视为气候变化热点区域。波斯栎（Quercus brantii Lindl.）是该区域的建群种，对生态系统稳定和水循环至关重要，同时支撑着约150万人的生计，提供放牧、农业和非木质林产品资源，且这一山区贡献了伊朗近40%的水资源。然而，快速变化的降水和温度已导致栎树严重枯梢和毁灭性昆虫爆发，威胁森林健康并造成社会经济影响。已有基于树轮宽度和木材解剖学的研究证实波斯栎对气候变率高度敏感，但关于其树轮纤维素稳定碳同位素（δ¹³C_TRC）的研究仍属空白，尤其缺乏沿海拔和湿度梯度对极端气候事件生理响应的认识。稳定碳同位素比值可记录植物气孔导度和光合活性相关的环境信号，并能用于计算内在水分利用效率（iWUE），从而揭示树木应对干旱胁迫的生理策略。因此，研究人员旨在通过建立δ¹³C_TRC和iWUE时间序列，分析三个不同海拔样地波斯栎对极端气候事件的生理响应，为理解栎树衰退机制及制定适应性管理策略提供依据。该论文发表在《European Journal of Forest Research》。

研究方法概述：研究区域位于扎格罗斯山脉的三个样地：伊泽（IZ，低海拔）、达劳（DA，中海拔）和萨尔阿加赛义德（SA，高海拔）。研究人员从60棵成熟健康树木中采集了120个增量钻芯，并基于与所有树木年轮宽度序列的强相关性，每个样地选取5个钻芯进行δ¹³C分析。为避免幼龄效应，选取1953—2022年共70个年轮，并用剃须刀片在双目显微镜下分离完整年轮（因干旱导致窄轮而无法分离早材/晚材）。采用多阶段α-纤维素提取法获得纤维素样品，使用元素分析仪-连续流同位素比质谱仪测量δ¹³C值，并利用已知标准进行校正。为去除大气CO₂浓度上升趋势的影响，研究人员采用McCarroll和Loader（2009）的方法对δ¹³C值进行大气趋势校正，并计算iWUE。气候分析使用最近气象站数据、CRU scPDSI格点数据、SPEI指数（1–12个月累积尺度）以及ERA5再分析资料（降水、平均温度、饱和水汽压差VPD和相对湿度rH）。统计方法包括相关分析、共线性分析、单因素方差分析、Kruskal–Wallis检验、Welch t检验、Mann–Whitney U检验及Cohen's d效应量计算。

研究结果：

（一）波斯栎树轮δ¹³C的长期变化及个体间相关性。大气趋势校正后的δ¹³C年表在各样地间表现出相对较弱的序列间相关性（Rbar为0.21–0.36）。高海拔SA样地的平均δ¹³C值（?24.04‰）显著富集于低海拔IZ和中海拔DA样地，而DA样地的平均δ¹³C值最低。两两统计比较证实，在1953—2022年期间三个样地间的树木平均校正δ¹³C存在显著差异，SA样地与另外两个样地相比效应量极大，而IZ和DA样地之间的差异仅在最近二十年显现且具有非平稳性。原始δ¹³C值均呈长期下降趋势，校正后多数树木序列仍保持同步，但部分个体表现出独特模式。

（二）波斯栎的长期植物生理响应。各样地平均δ¹³C年表显示，SA样地树木的δ¹³C值始终高于IZ和DA样地，1953—2008年间IZ样地高于DA样地。iWUE在所有样地均呈长期增加趋势，其中DA样地年变化率最高，IZ样地最低。1953—2022年期间SA样地iWUE持续最高，DA样地自1953至2008年最低，但2008年之后IZ样地转为最低。在极端干旱年份1964年、1984年和2007—2008年，所有δ¹³C年表均表现出较高（负值较小）的δ¹³C值，同期的iWUE也相应升高，表明干旱条件下树木主要通过降低气孔导度来提高iWUE。

（三）平均校正δ¹³C年表与气候数据的关系。与气象站降水数据的相关分析显示，IZ样地仅与前一年5月降水显著负相关；DA样地与前一年3—4月、当年3月及多个季节降水窗口显著负相关；SA样地与前一年1月、当年1月及若干季节降水窗口显著负相关。与温度的相关关系因样地而异：IZ样地与最高气温（T_max）和平均气温（T_mean）呈负相关，而DA和SA样地则与T_max呈正相关。与最低气温（T_min）的相关性在DA样地最强。所有样地δ¹³C与scPDSI（1953—2022）均呈显著负相关，其中DA和SA样地更为强烈。与SPEI-12的相关分析也显示相似模式，DA样地响应最强。基于ERA5数据的分析揭示了降水的负相关、rH的负相关以及VPD的相关方向在不同样地存在差异：IZ样地与当年7—9月干旱季VPD和温度显著负相关，而DA和SA样地则与VPD呈显著正相关。共线性分析显示，IZ样地中当年3月降水的纯效应贡献最大（17.46%），7月VPD与温度的共享效应为11.72%；DA样地中rH、VPD和降水共同解释了δ¹³C变异（共享效应8.44%）；SA样地前一年1月降水的纯效应高达15.87%。

讨论与结论：研究结果表明，δ¹³C的长期下降趋势既反映了大气¹³C的贫化，也反映了植物对上升CO₂的生理调整。高海拔SA和低海拔IZ样地的δ¹³C富集指示更高的iWUE，反映出干旱环境下保守型的水分利用策略；中海拔DA样地因相对湿润，气孔导度较高，碳同位素判别较大，iWUE较低。极端干旱期iWUE的同步上升表明，气孔调节是水分亏缺下碳同位素变异的主要机制。δ¹³C与降水、rH、VPD、PDSI和SPEI的关系显示，水分可利用性是控制波斯栎碳同位素判别的主要因子。不同海拔对VPD的相反响应揭示了低海拔样地以同化限制为主，而中高海拔样地以气孔调节为主。中高海拔样地对冬春降水的高度敏感说明，土壤水分补给对半干旱栎林的正常生理功能至关重要。长期iWUE上升反映出波斯栎正通过调整水分利用策略来适应大气CO₂升高和日益增强的水分胁迫，但反复出现的极端干旱仍是年际生理变率的主要驱动因子。研究结论强调，树轮δ¹³C和iWUE是评估波斯栎干旱敏感性和水分利用策略的可靠指标，揭示了该物种沿海拔梯度的高度生理可塑性。在未来气候持续暖干化背景下，中高海拔样地可能更易受冬春降水减少的影响，管理策略应优先保护土壤水分、减少人为干扰并保护干旱易发区的更新。未来结合δ¹⁸O年表、土壤水分观测和生理测量，将有助于更深入理解水源变化、气孔调节和光合作用在控制同位素信号中的相对作用，从而为扎格罗斯栎林的保护与适应性经营提供科学支撑。