稳定同位素揭示中东欧五种原生阔叶树种对比鲜明的干旱响应策略

《Frontiers in Forests and Global Change》：Stable isotope insights into contrasting drought-response strategies of five native broadleaf tree species in East Central Europe

　　

摘要 气候变化背景下的未来森林动态，在很大程度上取决于树种对环境胁迫加剧的生理驯化能力。研究人员将碳、氧稳定同位素分析与生长测量相结合，以评估大气CO2与干旱胁迫对匈牙利中部五种种群共存树种——田野枫（Acer campestre）、花白蜡（Fraxinus ornus）、土耳其栎（Quercus cerris）、柔毛栎（Quercus pubescens）和银毛椴（Tilia tomentosa）——的内在水分利用效率（intrinsic water-use efficiency, iWUE）与次生生长的影响。广义可加模型（generalized additive model, GAM）解释了iWUE的大部分变异（校正后R2 = 0.77）与生长的大部分变异（校正后R2 = 0.75）。同位素与iWUE年表显示，这些树种经历了日益加剧的水热气候胁迫，主要由温度上升与水汽压亏缺（vapor pressure deficit, VPD）增大所驱动。由上年9月至当年8月计算的12个月水分平衡（12-month water balance, WB12）是最强的气候预测因子，与碳同位素判别值（Δ13C）（r = 0.55–0.73）和iWUE（r = ?0.60 至 ?0.80）均呈强相关。尽管各树种的生理适应存在差异，但iWUE升高、生长停滞或下降的总体格局表明，大气CO2浓度升高未能完全抵消增暖与变干对树木功能造成的不利影响。与其他树种不同，土耳其栎（Q. cerris）的iWUE在2010年前后仅小幅上升，此后增幅更为显著。iWUE与生长之间的关系揭示了树种间在干旱响应策略上的根本差异：银毛椴（T. tomentosa）与花白蜡（F. ornus）采取更为保守的用水策略，iWUE急剧上升而生长下降，这可能有利于短期存活却增加长期脆弱性；相比之下，柔毛栎（Q. pubescens）与田野枫（A. campestre）在iWUE升高的同时维持了生长，提示其具有更宽的水力安全边际（hydraulic safety margin）或更深的扎根深度；土耳其栎（Q. cerris）呈中间型格局，生长起初稳定，而在严重干旱下显著下降。除这些生理差异外，立地条件与景观结构似乎也通过局地小气候效应（包括蒸散需求）调节树木的响应。这些发现表明，在中欧较干旱地区，森林组成可能发生演替，因为共存树种在抵御日益加剧的干旱胁迫的能力上差异显著。

论文解读

研究背景与立项依据

中东欧森林正同时承受三重压力：持续增暖、变干以及大气CO₂浓度不断上升。已有的树木年轮研究记录到径向生长下降、死亡率上升与衰退事件频发，但这些指标主要反映碳分配结果，难以直接揭示光合作用与蒸腾作用等生理过程。稳定同位素则可提供机理层面的信息：在有利环境下，同位素往往比径向生长对气候变异更敏感。iWUE（intrinsic water-use efficiency，单位失水所固定的碳）可通过树轮稳定碳同位素组成（δ¹³C）估算，而碳同位素判别值Δ¹³C在扣除大气δ¹³C长期下降趋势后，能分离出树木自身的生理信号；氧同位素组成（δ¹⁸O）则主要受叶片水分富集控制，可帮助区分δ¹³C变异究竟源于气孔关闭还是光合能力改变。全球尺度研究已证实iWUE长期上升，但这一上升可能来自CO₂同化增强、气孔导度下降或二者共同作用，而长期严重干旱下气孔关闭可能完全抵消高CO₂的有利效应。

关键的知识缺口在于：中欧许多次要落叶阔叶树种的生理响应知之甚少，且在同一气候背景下对多个共存阔叶树种开展同位素对比评估的研究十分稀缺；此外，立地条件与周边景观结构通过局地小气候对树木响应的调制作用长期被忽视。该论文发表于《Frontiers in Forests and Global Change》，正是针对上述不足开展研究。

主要技术方法

研究人员在匈牙利中部的两个干旱混合阔叶林站（G?d?ll? Hills—GOD，47.58°N/19.42°E；Somogy Hills—SOM，46.86°N/17.94°E）设置样地，选取五种种群共存的优势木作为样本队列：GOD采样田野枫（A. campestre）、花白蜡（F. ornus）、柔毛栎（Q. pubescens）、土耳其栎（Q. cerris）4种，SOM采样银毛椴（T. tomentosa）与土耳其栎（Q. cerris），每个种—站点组合选取5株优势木。方法上整合了四类技术：其一，利用生长锥在胸高（1.3 m）取一个芯，经交叉定年（COFECHA验证）后测定年轮宽度并计算基础面积增量（basal area increment, BAI）；其二，按年解剖分离早材与晚材，仅取当年光合产物为主的晚材（田野枫因无法可靠分离而用整轮），混池后提取α-纤维素并同步测定δ¹³C与δ¹⁸O；其三，依据Farquhar模型由Δ¹³C计算iWUE；其四，利用Mann-Kendall趋势检验（趋势自由预白化TFPW、Sen斜率）与Pearson相关分析气候—生理关系，并以广义可加模型（GAM）量化iWUE与生长的驱动因素。

研究结果

3.1 气候、Δ¹³C、δ¹⁸O、iWUE与BAI的格局

近几十年年均温显著上升（Kendall's τ = 0.48–0.49，p < 0.05），尤其在1990年后，而VPD亦显著上升（τ = 0.39–0.50）；WB12年际波动大但无明显趋势。作为局地小气候代理的实际蒸散（actual evapotranspiration, AET）显示，2000–2003年严重干旱期GOD站点及邻域夏季AET比SOM高约10%（105 mm vs 96 mm）。各树种Δ¹³C均值相近，仅花白蜡在两地显著偏高；Mann-Kendall检验显示除GOD的土耳其栎Δ¹³C显著上升外，其余趋势多不显著。δ¹⁸O在GOD多稳定或下降、在SOM则显著上升；平均iWUE以花白蜡最低，所有树种iWUE均显著上升，但土耳其栎直到约2010年前增幅平缓（Sen斜率GOD 0.21 year^?1、SOM 0.29 year^?1），而银毛椴（0.64 year^?1）与田野枫（0.63 year^?1）增幅最大。δ¹⁸O对极端气候的响应具明显种—站点差异：GOD的银毛椴在干旱年δ¹⁸O显著升高，田野枫与柔毛栎变化克制，土耳其栎居中；SOM的银毛椴与土耳其栎则呈干旱年升高、湿润年显著偏负的格局。

3.2 Δ¹³C、δ¹⁸O和iWUE与气候变量的关系

相关分析表明，Δ¹³C与春夏季降水一致呈正相关（r = 0.25–0.50），与温度关系较弱；5–8月水分平衡（water balance, WB）与Δ¹³C中等相关（r最高0.5），而WB12为最强预测因子（r = 0.55–0.73），仅田野枫更受当年夏季WB影响。δ¹⁸O在SOM与夏季气候的关系更强、更一致，GOD的田野枫甚至无显著关系。iWUE与WB12呈强负相关（r = ?0.60 至 ?0.80），与VPD亦强正相关（r = 0.64–0.75）。Δ¹³C与δ¹⁸O总体负相关（r = ?0.63 至 ?0.57为最强），提示两同位素体系在部分树种中紧密耦合。

3.3 跨树种iWUE的驱动因素及其对生长的影响

GAM分析显示，iWUE模型校正后R² = 0.77、解释78%偏差，受气候水分供应、大气CO₂浓度与树种身份显著影响，树龄仅有 modest 的非线性贡献；WB12呈强负线性效应，VPD效应较弱且边际显著，而CO₂ × 树种（CO₂ × SPEC）交互显著，说明各树种对CO₂升高的iWUE响应幅度不同。生长模型校正后R² = 0.75、解释76.6%偏差，树龄效应显著（edf ≈ 5.4）；iWUE对各树种生长的光滑项均显著，多数接近线性（edf ≈ 1）。偏效应曲线揭示了对比鲜明的响应：银毛椴生长随iWUE升高急剧下降，花白蜡在低iWUE段响应较弱而高iWUE段陡降，反映了水力安全与碳同化之间的权衡与阈值行为。

讨论

4.1 Δ¹³C、δ¹⁸O与iWUE的时间变化及其与气候变量的关系

Δ¹³C稳定或上升意味着生理判别信号保持恒定或增强，残差变异主要捕捉树木对水分与蒸发条件的响应；其与夏季气候的强正相关（尤其在GOD）表明气孔导度是碳同位素判别的关键驱动，气孔调节在干旱期维持叶片水分状态。δ¹⁸O的站点分异（GOD稳定/下降、SOM上升）凸显了生态水文条件对叶片蒸发过程的塑造：SOM的δ¹⁸O与夏季温度、降水、WB关系更强，反映叶片水蒸发富集与大气干燥度紧密耦合，而GOD关系较弱，可能源于更深土层水分利用或叶片水周转的种间差异。这一差异可归结为景观结构——GOD处于更大、更连续的森林之中，森林覆盖缓冲局地极端小气候、降低蒸发需求；SOM则被农业用地环绕，更易受暖干气流影响。AET在严重干旱期GOD比SOM高约10%支持了这一解释，但直接测量小气候数据缺失，加之两站点均无地下水补给、海拔相近、相距约110 km，降水源δ¹⁸O差异有限，故景观缓冲作用仍属合理推断而非实验验证。

4.2 种特异性的iWUE驱动因素及其对生长的影响

花白蜡属典型"耗水型"策略，两林分iWUE均较低、气孔较开放，有利于水分供应有利时持续同化与生长；2010年后随着干旱频率与VPD上升，其策略转为保守，iWUE陡升。相反，银毛椴与田野枫的iWUE上升更渐进、更一致，暗示气孔层面对变干更早、更敏感的响应。δ¹⁸O格局同样具种—站点特异性：GOD的银毛椴同位素敏感性最高，田野枫与柔毛栎响应平衡，土耳其栎居中；SOM的银毛椴与土耳其栎则波动更大。iWUE与生长的解耦——2010年后多树种iWUE陡升而生长停滞或下降——说明高iWUE并不必然转化为高碳增益或高生产力，而更可能代表以牺牲碳同化为代价维持叶片水分的保守策略，可作为干旱加剧与未来生长下降的早期预警。Δ¹³C与δ¹⁸O的耦合程度也随树种与站点变化，银毛椴与花白蜡在GOD耦合紧密，而田野枫与柔毛栎在GOD关系较弱，提示后者Δ¹³C受非气孔（生化）限制的影响更大。iWUE模型中WB12的主导作用确认土壤水分亏缺是iWUE上升的主驱动（气孔调节机制），而VPD效应仅轻微；显著的CO₂ × 树种交互则反映不同水力结构、气孔敏感性与光合能力的差异。生长模型进一步表明，水分利用效率与生长的关系并非一致：银毛椴生长随iWUE陡降，花白蜡呈非线性阈值响应，支持"CO₂施肥效应与干旱的博弈决定iWUE是否转化为生产力收益"的判断，也凸显根深度、水力策略等功能性状在预测未来森林响应中的重要性。

4.3 局限性

δ¹⁸O信号受树龄影响较小，故识别的生理格局不太可能完全由个体发生效应驱动。景观介导的小气候缓冲作用需谨慎解释：两站点区域气候相似、仅在周边森林覆盖与干旱年AET上存在差异，但缺乏林分内气温、相对湿度、土壤水分等连续直接测量，该效应属间接证据支持的可能解释，需未来在更多站点开展连续微气候监测加以检验。此外，由于无直接气体交换测量，iWUE上升不能唯一归因于气孔导度下降，也可能反映光合能力下降或二者兼有。

结论（翻译）

长期的同位素与生长分析表明，在中东欧较干旱地区，所研究的五种原生阔叶树种已因温度上升与干旱频发而承受相当程度的气候胁迫。尽管过去几十年所有树种的iWUE均有所提高，但这些变化并未伴随持续而增大的生长，提示大气CO₂升高未能补偿增暖与变干造成的不利影响。各树种表现出多样的干旱响应策略：银毛椴与花白蜡采取更为保守的用水策略，iWUE急剧上升而生长明显放缓，这可能有利于短期存活却可能导致碳耗竭并增加脆弱性；相比之下，柔毛栎与田野枫能够在iWUE升高的同时维持生长，提示其水力特性或扎根深度缓和了干旱的负面效应；土耳其栎呈中间型格局，在中等胁迫下生长近乎恒定，而在严重干旱下显著下降。除种间差异外，研究结果提示立地条件与周边景观结构也可能通过局地小气候效应（包括蒸散需求）影响干旱响应，尽管树种身份仍是所观测格局的首要驱动因子。综上，这些发现突出了树种间在水力与生理干旱响应策略上的强烈差异。