论文部分内容阅读
本文以黑龙江省佳木斯市孟家岗林场落叶松人工林为研究对象,利用7块落叶松人工林标准地的49株解析样木树干解析数据,通过年轮图像分析系统WinDENDORTMV6.5查数年轮所获得共1186个圆盘数据。分析了落叶松树皮厚度(BT)和树皮因子(DIB/DOB)随树高(HT)和胸径(DBH)的变化规律,并利用SAS 9.4软件中的PROCMIXED模块,在所选最优基础模型的基础上,考虑树木效应、样地效应,分别构建了落叶松人工林树皮厚度(树皮因子、任意高度处树皮厚度)的线性混合效应预测模型。模型评价指标选用赤池信息准则(AIC)、贝叶斯信息准则(BIC)、-2倍的对数似然值(-2LL)及似然比检验(LRT)。1)树皮因子的结果表明:树皮因子沿树干垂直方向随相对高度增加基本保持恒定。当相对高度超过80%后,树皮因子随相对高度增加而减小。树皮因子与胸径无显著相关性。任意参数组合的混合效应模型的拟合效果均好于传统模型。对于树皮因子模型,基于树木效应时含b1,b2,b4即相对高度(h/HT)和相对高度的平方((h/HT)2)、胸径处去皮直径与胸径的比值(DBHIB/DBH)组合的为最优混合模型。调整后确定系数R2为0.734,平均绝对误差MAE为0.022,均方根误差RMSE为0.032;基于样地效应时含b1,b2即相对高度(h/HT)和相对高度的平方((h/HT)2)随机参数组合的为最优混合模型。调整后确定系数R2为0.696..平均绝对误差MAE为0.022,均方根误差RMSE为0.034;所有最优模型在具有无结构(UN)方差-协方差矩阵时拟合效果最好。2)任意高度处树皮厚度的结果表明:树皮厚度随相对高度增加呈S型曲线变化,整体与相对高度呈负相关,相对高度为0～20%区间内,树皮厚度随高度增加减小的速率很快,在20%～80%区间内,树皮厚度随高度增加减小速率变缓,相对高度80%到梢头,树皮厚度减小速率又加快。任意高度处树皮厚度与胸径存在显著的正相关性。任意参数组合的混合效应模型的拟合效果均好于传统模型。对于任意高度处树皮厚度模型,基于树木效应时含b1,b2即相对高度(h/HT)和任意高度处带皮直径(dob)组合的为最优混合模型。调整后确定系数R2为0.878,平均绝对误差MAE为0.053,均方根误差RMSE为0.083;基于样地效应时含b0,62,b3即截距项、任意高度出带皮直径(dob)、径高比(DBIH/HT)组合的为最优混合模型。调整后确定系数R2为0.820,平均绝对误差MAE为0.061,均方根误差RMSE为0.093;所有最优模型在具有无结构(UN)方差-协方差矩阵时拟合效果最好。不论是树皮因子还是树皮厚度模型,树木效应对模型的影响是最大的。混合效应模型的预估精度相比于传统回归模型有明显提高。树皮因子和树皮厚度混合效应模型的构建可以更加准确预测立木去皮材积,树皮材积,为实际森林经营和木材生产提供帮助和指导。