树木生长监测中径向变化仪的机遇与局限：基于微芯样本的木质部物候学评估

本刊推荐：为评估径向变化仪在监测木材物候关键日期中的可靠性，研究人员对6种具不同树轮结构（无孔、散孔、环孔）和树皮类型（光滑、鳞状、裂隙）的树种，同步采用径向变化仪（DBH测量）与微芯取样（xylogenesis观察）法进行对比。结果表明，径向变化仪的估算存在显著偏差，其准确性受物种生活性状（如树轮结构、树皮类型）和环境因素（如降水、海拔）影响。该研究发表于《Quantitative Plant Biology》，强调需结合木质部形成监测以校正径向变化仪的局限性，对气候变化下森林生产力与碳汇评估具有重要方法论意义。

随着全球气候变化对森林生态系统的影响日益显著，树木的生长动态与碳汇功能成为生态学研究的热点。木材形成是树木适应环境的核心过程，其物候时序——即形成层细胞分裂、膨大、壁增厚和成熟的起止时间——直接关系到年轮结构与生物量积累。传统上，研究者通过微芯取样（microcoring）结合组织学观察来精确追踪木质部形成过程，但该方法耗时耗力且对树木有损伤。为便于大尺度长期监测，径向变化仪（dendrometer）被广泛用于记录树干半径或周长变化，进而推估生长物候。然而，径向变化仪所测信号实为木质部、韧皮部及树皮水分变化的综合响应，其能否真实反映木材形成的临界日期一直缺乏系统验证。

为厘清径向变化仪在木材物候监测中的准确性与局限性，Ignatius Kristia Adikurnia与Cyrille Rathgeber*在《Quantitative Plant Biology》上发表研究，通过对6种温带树种（银冷杉、欧洲云杉、欧洲落叶松、欧洲赤松、欧洲山毛榉和无梗花栎）同步开展径向变化仪测量与微芯取样，对比分析了两种方法对生长季开始（细胞膨大起始日b_E）与结束（细胞膨大终止日c_E）的估算差异。这些树种覆盖了关键的生活性状差异：系统发育类群（裸子植物与被子植物）、树轮结构（无孔均质、无孔异质、散孔、环孔）及树皮类型（光滑、鳞状、裂隙），并跨越了不同海拔、气候条件的站点（Donon山、Mélézin梯度、Hesse森林）。

研究团队通过三个独立数据集（Donon、Mélézin、Hesse）获取了265棵树年的同步监测数据。微芯样本经脱水、石蜡包埋、切片染色后，在显微镜下统计细胞分化阶段；径向变化仪数据则采用五参数逻辑斯蒂生长模型（RDendrom包）拟合，以5%（DOY.05）和95%（DOY.95）年生长量完成度作为生长季起止的代理指标。通过绝对差异（AD）、平均绝对差异（MAD）和简化主轴回归（RMA）等统计方法，评估两种方法的一致性。

关键方法概述

研究整合了来自法国和爱沙尼亚的三种数据集（Donon、Mélézin、Hesse），涵盖6树种265棵树年的同步微芯取样与手动带状径向变化仪（EMS DB20）测量。微芯样本经Trephor®工具每周采集，通过显微解剖学观察木质部细胞分化阶段；径向变化仪数据通过逻辑斯蒂模型拟合DBH变化，提取DOY.05与DOY.95作为物候代理指标。统计分析包括残差评估、RMA回归及环境因子（降水、海拔）与生活性状（树轮结构、树皮类型）的关联检验。

结果1：逻辑斯蒂模型的拟合优度

整体上，径向变化仪记录的S形生长曲线能被五参数模型较好拟合（平均RMSE=0.17±0.06 mm），但物种间差异显著。欧洲山毛榉拟合最佳（RMSE=0.10 mm），而无梗花栎最差（RMSE=0.24 mm）。残差分析显示，模型在生长季边界（尤其是5%和95%阈值附近）的拟合精度较低，且裸子植物与被子植物对阈值响应的模式不同。例如，落叶松和栎树在生长季末期常出现树干收缩或持续增粗，导致95%阈值估算偏差增大。

结果2：径向变化仪与微芯估算值的比较

径向变化仪对银冷杉和欧洲云杉的生长季起止估算与微芯法无显著差异，但对其他物种存在系统性偏差。生长季开始时间上，径向变化仪普遍滞后（落叶松、山毛榉、赤松、栎树的MAD分别为9±24、10±6、13±13、27±11天）；生长季结束时间上，径向变化仪对落叶松和山毛榉提前估算（MAD=0±62、12±10天），对赤松和栎树则滞后（MAD=28±48、8±23天）。RMA回归进一步揭示，仅山毛榉和银冷杉在生长季起止均存在显著线性关系，但存在偏移或斜率偏差；落叶松、赤松和云杉在生长季结束时的估算甚至无显著关联。

结果3：生活性状对估算准确性的影响

被子植物在生长季开始时的MAD显著高于裸子植物，而结束时相反。树轮结构异质性越高（如环孔栎树、异质无孔赤松），MAD越大；裂隙树皮（赤松、落叶松、栎树）的MAD显著高于光滑或鳞状树皮。物种聚类分析显示，生长季开始时MAD从冷杉/云杉/山毛榉→落叶松/赤松→栎树递增；结束时则从山毛榉/栎树/冷杉/云杉→赤松/落叶松递增。树皮厚度与MAD无显著关联。

结果4：环境因子的调节作用

春季降水增多会显著增加生长季开始时的MAD（尤其在冷杉、山毛榉和栎树中），但秋季降水对结束时的MAD无一致影响。海拔效应在落叶松中显著：从低海拔（1350 m）到高海拔（2300 m），MAD逐渐减小，表明径向变化仪在寒冷站点准确性更高。

结论与讨论

本研究系统揭示了径向变化仪在估算木材物候关键日期中的局限性：其信号受木质部生长、韧皮部活动及树皮水分变化的共同干扰，且准确性高度依赖物种生活性状（树轮结构、树皮类型）和环境条件（降水、温度）。尽管逻辑斯蒂模型能良好拟合生长曲线，但基于阈值的临界日期提取方法在生长季边界可靠性低。尤其对于裂隙树皮物种或湿润春季，径向变化仪易低估生长开始时间；而对平滑树皮物种，生长结束时间可能被高估。这些偏差质疑了仅靠径向变化仪评估气候变化下森林物候响应的结论。

作者建议，未来研究应结合微芯法的直接观测与径向变化仪的大尺度部署优势，通过校正模型（如广义加性模型GAM）或分析昼夜径变周期（而非累积生长曲线）来提升物候估算精度。该研究为木质部形成监测提供了方法论警示，对精准评估全球变暖背景下森林碳动态具有重要实践意义。