Spatial and temporal distribution characteristics and discoloration law of knots in the Mytilaria laosensis plantation

中国林业科学研究院热带林业实验中心(22°11′N，106°75′E)地处广西最南部，与越南谅山市接壤，试验区下设哨平、青山、伏波和白云4个分场。米老排人工林位于青山分场，海拔约440 m。该地区属于亚热带季风气候，雨量充沛，夏无酷热，冬无严寒，年平均气温为 22.5 ℃，年平均降水量约1 210 mm，雨季集中在4—10月，约占全年降水量的3/4，全年无霜期约345 d，年平均日照约1 620 h。土壤主要为砖红壤，呈酸性(pH 4.5)，土壤有机质质量分数为24.65 g·kg⁻¹，全氮质量分数为1.02 g·kg⁻¹，全磷质量分数为 0.47 g·kg⁻¹，有效氮质量分数为117.20 mg·kg⁻¹，有效磷质量分数为1.20 mg·kg⁻¹，有效钾质量分数为84.20 mg·kg⁻¹。

米老排人工林建于2008年，采用来源于同一优良母树的种子培育的实生苗进行造林，初始种植密度为1 667 株·hm⁻²(株行距2 m×3 m)，该样地在培育过程中没有进行间伐和施肥等处理。2019年2月，在该样地中设置了5块20 m×30 m的标准地，对每块样地内的米老排进行编号，并测定每棵树的胸径、树高、枝下高、冠幅等因子。各样地调查数据如表1所示。

依据每木检尺数据，在每个标准地中分别选出1株优势木、1株亚优势木和2株标准木进行节子解析试验，共计20株。用喷漆标记树干的北方向位置，并用胸径尺测定其胸径后，用油锯将样树伐倒，再次测定其树高、枝下高和冠幅等因子(表2)。以2 m为区分段，用油锯进行切割，分别统计整株节子的数量及其所在高度。随后，将各区分段有节子的部分锯切成 20~30 cm 的圆盘，随后使用油锯沿节子中心部分进行纵切，得到节子剖面(图1)，抛光至清晰，最终得到1101个完整的节子。

Figure 1.  Longitudinal section of the occluded branch (modified from HEIN^[12]) and the actuallongitudinal section of a knot

将所测定的节子指标作为自变量(x₁~x₄)，节子变色长度为因变量(y)。进行逐步回归分析，筛选出影响节子变色长度的关键因素。并构建多元线性归回模型：y=a₀+a₁x₁+a₂x₂+a₃x₃+a₄x₄，其中a₀、a₁、a₂、a₃、a₄为参数。

采用SPSS 25.0对数据进行统计分析，使用最小显著性差异法(LSD)进行多重比较，显著性水平为α=0.05，用Excel制表，利用Origin软件制图。

图2显示：米老排大多数枝条在第1~2年开始生长；大量枝条在第3~4年死亡，较少会在第5~8年死亡；第5~8年是节子愈合的高峰期，约占所有节子的87.2%，少部分节子会在第9年痊愈。在枝条死亡脱落后，枝条周围的形成层逐渐将树干内部的枯枝残留包围形成节子，76.7%的节子需要3~6 a的时间才能完全包裹，而能在1~2 a被完全包裹的节子数量较少，仅为10.4%，12.9%的节子需要7 a以上的时间才能完成包裹。

Figure 2.  Temporal distribution of branch occlusion of M. laosensis

如图3所示：在垂直水平上，米老排节子数量随树干高度增加呈先上升后下降的趋势。以2.0 m为区分段，在树干高度为2.0~4.0 m内节子数量百分比最高，约占整株节子总数的17.8%。在2.0~4.0 m之后，随着高度的增加，节子数量开始逐渐减少。从节子数量累计百分比看，61.9%节子集中在0~8.0 m。

Figure 3.  Vertical distribution of M. laosensis knots

由表3可知：节子的平均变色长度达102.5 mm。由图4可知：节子的变色长度与节子直径呈正比。除直径40.0~50.0 mm与直径50.0~60.0 mm间的节子变色长度无显著差异外，其他不同节子直径间变色长度均存在显著差异(P＜0.05)。节子直径为0~10.0 mm时节子平均变色长度最小，仅为58.4 mm，而当节子直径为50.0~60.0 mm时，节子平均变色长度最大，为136.0 mm。

Figure 4.  Effects of knot diameter on discoloration length

由图5可见：米老排节子变色长度随着分枝角度的增加，呈先增大而后减小的趋势，其中节子分枝角区间在80.0°~100.0°的节子变色长度显著小于其他分支角区间(P＜0.05)。当节子的分枝角度在80.0°~100.0°内时，节子平均变色长度最小，仅为60.0 mm；当分枝角度在20.0°~40.0°内时，节子平均变色长度最大，达110.2 mm。

Figure 5.  Effects of branch insertion angle on discoloration length

由图6可见：米老排节子变色长度与愈合时间呈正相关，且当节子愈合时间不超过6 a时，变色长度随愈合时间的增加而缓慢增大，而当愈合时间超过6 a，变色长度随愈合时间的增加增速变快。

Figure 6.  Effects of knot occlusion time on discoloration length

将枝条出现年龄、枝条死亡年龄、节子愈合时间、分枝角度、死节长度、节子着生高度以及节子直径等7个指标作为自变量，与节子变色长度进行逐步回归分析，得到描述节子变色长度(y_DL)的最优模型：y_DL=24.388+1.729x_KD+0.350x_DBL−6.635x_YB+5.451x_OT。其中：x_KD为节子直径，x_DBL为死节长度，x_YB为枝条出现年龄，x_OT为节子愈合时间。模型中各参数均具有统计学意义，模型的R²为0.649(表4)。模型的残差分布均匀，表明方程模拟效果较好(图7)。

Figure 7.  Residual distribution of discoloration length model

本研究发现：米老排枝条在第2~4年达到死亡的高峰期，树木上枯枝残留在后期被树干完全包裹形成节子。而水曲柳Fraxinus mandshurica节子形成高峰期在第6~15年^[24]，格木Erythrophleum fordii在第11~20年是控制节子的关键时期^[25]。这些差异主要来源于以下几个原因：一是树种差异，不同树种之间树干径向生长速率不同^[13]；二是生长立地条件不同，在高密度种植情况下，不同个体之间竞争有限的周围环境资源，会导致个体分配到枝条生长的能量受到影响^[26]。本研究将枝条死亡年龄作为米老排修枝的判断依据，在米老排人工林生长第2~4年可对其进行首次修枝。

树干上的枝条死亡后，死枝残桩有时会残留在树干内，阻止树干外部的愈合，延长了枝条愈合所需时间，增大树干内部与外界水分和真菌等接触的时间，从而增加了木材内部受到侵染而发生腐坏的概率^[27]。本研究发现：米老排人工林多数节子愈合需要3~6 a，甚至更长的时间。有研究表明：节子愈合时间与节子的直径呈正比^[28]，与树干径向生长速率呈反比^[12]，因此，可适当对林木进行间伐或者施肥等措施，以促进树干的径向生长速率，缩短节子愈合所需时间。在米老排垂直分布上，其半数以上节子主要集中在树高8.0 m以下部分，这一结果与楸树Catalpa bungei^[29]、水曲柳^[24]的研究结果相似，据此可以将8.0 m作为米老排修枝的参考高度。

在本研究中：节子变色长度随节子直径的增大而增大，当节子直径不超过50.0 mm时，节子变色长度会随着节子直径的增加而显著增大。关追追等^[30]研究发现：当水曲柳节子直径超过15.0 mm时，节子变色长度会随着节子直径的增加而明显增大。而DǍNESCU等^[31]认为：欧亚槭Acer pseudoplatanus节子直径超过50.0~60.0 mm时，木材变色率才会明显提高。说明米老排人工林节子在真菌等微生物侵染的情况下更容易发生变色。节子变色长度随分枝角度的增加而减小，且在本研究中当分枝角度大于60.0°时，节子变色长度随分枝角度的增加而显著降低。这与对西南桦^[32]的研究结果一致，而与水曲柳^[24]的研究结果相反。在本研究中，米老排节子愈合时间超过6 a时，变色长度明显增大。本研究运用逐步回归分析，筛选出节子直径(x_KD)、死节长度 (x_DBL)、枝条出现年龄(x_YB)和节子愈合时间(x_OT) 4个关键因子，得到描述节子变色长度(y_DL)的最优模型：y_DL=24.388+1.729x_KD+0.350x_DBL−6.635x_YB+5.451x_OT，相关性达到了极显著水平(P＜0.01)。计算结果与实际测量结果之间差异极小，因此可以利用该模型预测米老排枝条形成节子后对木材的影响状况。下一步研究可以对不同年龄、密度和立地条件的林分进行取样分析。

米老排人工林在自然整枝条件下，大多数枝条的死亡时间在第2~4年，节子愈合年龄主要集中在第5~8年。节子直径、节子愈合时间、死节长度和枝条出现年龄是决定节子变色长度的关键因子。当节子直径在不超过50.0 mm和愈合时间不超过10 a时，节子变色长度均随直径和年龄的增加显著增加。因此，为了减少节子变色长度，提高米老排木材的质量和出材率，可以在米老排生长的第2~4年，对树干8.0 m以下的米老排人工林进行首次修枝。