怎么才是闽楠（欧阳勋志教授团队论文）

原创 单凯丽等 南京林业大学学报

论文推荐

赣中天然闽楠单木冠幅预测模型的研究

单凯丽，臧 颢，潘 萍，宁金魁，欧阳勋志*

鄱阳湖流域森林生态系统保护与修复国家林业和草原局重点实验室

树冠是树木积累有机物的主要场所,它决定了树木的生产力和活力,其大小的评判指标有冠幅、树冠横断面积、树冠长度、树冠率、冠基高度等，其中冠幅是最易准确测量和最常用的指标之一,常被作为单木生长模型、生物量模型、枯损模型等中的预测变量,因此,越来越多国内外学者开始构建冠幅预测模型。近年来,有学者利用树木因子、林分因子和立地条件等建立了冠幅模型,如Thorpe等以不列颠哥伦比亚省中北部的3个树种为研究对象,用空间邻域法比较竞争对冠幅模型的影响,发现添加与相邻木大小和距离相关的竞争因子可极大提高模型精度;雷相东等以落叶松云冷杉林9个树种为研究对象,发现林分密度指数、林木竞争指数、树冠比和枝下高等因子对冠幅模型有显著影响;贺梦莹等以长白落叶松(Larix olgensis )-水曲柳(Fraxinus mandshurica)混交林为研究对象,发现与距离无关的竞争因子对冠幅具有直接影响;覃阳平等以杉木(Cunninghamia lanceolata) 为研究对象,构建与单木大小和单木竞争相关的冠幅模型,结果发现与竞争有关的冠幅模型精度显著高于与竞争无关的冠幅模型。总体看来,不少学者研究表明竞争对冠幅预测有显著影响。

构建冠幅模型的方法多从传统最小二乘法向混合效应模型转变,混合效应模型处理重复测量数据和纵向数据有一定优势,此模型的随机效应参数可以解释不同区域或样地间的冠幅曲线变化，固定效应参数可以解释冠幅的平均水平变化。如Raptis 等以黑松(Pinus thunbergii)为研究对象，构建样地水平的非线性混合效应冠幅预测模型,结果表明混合效应模型的预测能力优于传统回归模型;符利勇等和韩艳刚等研究树木冠幅模型时均得出相同结果。近年来,由于机器学习具有对自变量和因变量不需要做分布假定等优点,在生态学、林学等各领域受到广泛关注,如分析气候变量与木材密度的关系、预测林分生物量、立地质量评价等,主要形式有增强回归树(boostedregression trees， BRT)和随机森林(random forest，RF)等,然而应用机器学习建立冠幅预测模型的研究目前尚鲜见报道。

闽楠(Phoebe bournei)为樟科( Lauraceae)楠属(Phoebe)的常绿乔木树种,其材质优良,纹理美观，具有重要的经济和生态价值。受各种因素的影响导致其近于枯竭,现为国家二级濒危保护植物,我国仅在福建、江西、湖南、湖北、浙江等海拔1 000 m以下的山地常绿阔叶林中有零星分布。为了解闽楠的生长规律和提高闽楠天然次生林的林分质量，不少学者逐渐对现存的天然闽楠开展了相关模型研究:曹梦等建立了闽楠的单木生长模型和相容性生物量模型,褚欣等建立了闽楠直径分布预估模型,欧建德建立了林下闽楠更新层生境质量模型等,但对闽楠冠幅预测模型的研究还鲜见报道。本期论文推荐的作者以江西省吉安市的天然闽楠为研究对象,采用普通最小二乘法(OILS模型)、混合效应模型、增强回归树(BRT)和随机森林(RF)4种建模方法添加某一种竞争指标,构建单木冠幅预测模型,选出拟合和预测最优的模型形式来分析竞争对冠幅模型的影响,为提高冠幅的预测精度和科学经营提供参考。

下面跟学报君一探究竟！

作者简介

通讯作者

欧阳勋志，男，博士、二级教授，博士生导师，主要从事森林资源监测与经营管理等方面的研究。国家林业和草原局森林经营工作专家组委员，江西省百千万人才工程人选、江西省高等学校中青年学科带头人、江西省林业科技推广智库专家、江西农业大学森林经理学学科带头人。兼任江西省林业经济学会理事长、江西省林学会常务理事、中国林学会森林经理分会常务理事、中国林学会林业计算机应用分会常务理事、中国系统工程学会林业系统工程专业委员会常务委员等职。2011年以来主持国家自然科学基金4项以及国家林业公益性行业科研专项等20余项科研项目。获江西省科技进步二等奖2项、梁希林业科学技术奖二等奖2项等奖项。发表学术论文150余篇，出版著作3部。

第一作者

单凯丽，女，1995年11月生，江西农业大学林学院林业硕士研究生，主要从事森林资源监测与管理研究。

关键词：闽楠；冠幅模型；竞争指标；混合效应模型；刀切法

基金项目：国家林业公益性行业科研专项项目(201504301) ;国家自然科学基金项目(31760207)。

引文格式：单凯丽,臧颢,潘萍,等.赣中天然闽楠单木冠幅预测模型的研究[J].南京林业大学学报(自然科学版),2022,46(1):89-97.SHAN K L,ZANG H,PAN P,et al.Individual tree crown width prediction models for natural Phoebe bournei in central Jiangxi Province[J].Journal of Nanjing Forestry University (Natural Sciences Edition), 2022,46(1):89-97.DOI:10.12302/j.issn.1000-2006.202006045.

1目的

探讨竞争指标和建模方法对天然闽楠(Phoebe bournei)单木冠幅预测模型的影响,以期为精准预测天然闽楠单木的冠幅提供参考。

2方法

以江西省赣中25块闽楠天然次生林典型样地中的1011株闽楠为研究对象,采用普通最小二乘法(OLS模型)仅考虑样地水平的混合效应模型、增强回归树和随机森林4种建模方法建立单木冠幅模型,分别添加林分每公顷断面积、大于对象木的断面积之和、林分密度指数和简单竞争指数4种竞争指标,分析竞争指标对冠幅模型的影响,采用决定系数(R²)、均方根误差(RMSE)、平均相对误差绝对值(RMA)和平均绝对误差(MAE)确定最佳模型。

2.1 研究区概况

研究区位于江西省的吉安市(113°48’~115°56'E,25°58'~27°57'N),属亚热带季风气候，年平均气温17.1~18.6 ℃ ,年平均降水量1 360~1 577 mm,年无霜期281 d,地形以丘陵、山地为主，母岩以花岗岩、砂岩为主,土壤以红壤为主。植物资源丰富,森林覆盖率达到67.6% ,该地植被类型主要以常绿阔叶林为主,针叶林主要是以杉木、马尾松(Pinus massoniana)为主。

▲现场定样地范围

▲量胸径

2.2 数据获取

本研究对赣中吉安市的安福县、井冈山市及遂川县的闽楠天然次生林分布区域进行踏查,选取人为活动干扰程度相对较轻且林分在分布区域内具有代表性的地块设置典型样地进行调查。样地大小依据其分布地形和群落分布情况等因素而定,样地面积为400 ㎡(20 m×20 m)或600 ㎡(20 m×30 m),总共25块样地。样地调查主要内容包括胸径2 cm 以上林木的胸径(DBH)、树高(H)、位置(X和Y坐标)以及东南、西北2个方向的冠幅半径等基本因子。共调查样木1 778株,其中闽楠1011株,本研究所指单木冠幅( crown width , WC)为东西、南北两个方向的平均值。样地基本概况见表1。

▼表1 样地数据统计表

2.3 模型选取

2.3.1 OLS模型

冠幅模型的建立是为了预测单木冠幅时更加简便和准确,因此需要选用容易获取的测树因子作为自变量。在以前的研究中,最常见的冠幅基础模型仅以胸径为解释变量,本研究选用11种常用的冠幅(WC)-胸径(DBH)模型(包括线性、非线性方程)作为构建闽楠单木冠幅模型的备选基础模型(表2),通过评价指标选出精度最优的基础模型。

▼表2 备选基础模型

注:WC为单木冠幅;DBH为胸径DBH;a0 b0 c0为模型参数

2.3.2 基于样地水平的混合效应模型

用混合效应模型可以有效地处理重复测量的数据,并且克服一般线性模型中反应变量必须具有独立和等方差的缺点。本研究建立基于样地水平的单水平混合效应模型,在选取最优的OLS基础模型中添加竞争因子,模型的所有参数均视为混合效应参数,以不同的组合形式进行混合效应模型拟合和检验,比较模型的精度,其一般形式为:

式中:β为固定效应向量;uk为k样地的随机效应向量;DBH为林木胸径;εk为误差项,且与参数uk相互独立;P表示uk或εk服从正态分布;D为样地对应的随机参数的方差协方差结构矩阵;Qk为k样地内误差效应的方差协方差矩阵。

2.3.3 增强回归树

增强回归树(BRT)结合了统计和机器学习的两种算法,其核心思想是通过不断地随机选择和自我学习生成多重分类回归树,再将这些回归树集合起来,构成一个更大的分类树,以此来提高稳定性和精度。在模型拟合时,共有4个参数需要设置:①损失函数的形式为“gaussian”;②收缩参数(shrinkage) ,决定了单棵决策树在模型建立过程中的相对贡献率,本研究设置为0.05;③树的棵数，控制了模型的迭代次数,过大或过小会有过于拟合或欠拟合的风险,本研究设置为1 000;④交互深度(interaction depth)采用交叉验证。

2.3.4 随机森林

采用随机森林(RF)这一集成学习算法组合若干互不相关的决策树。在模型拟合时,共有2个参数需要设置:①树的棵数(n),一般而言,当棵树在n>500时整体误差率便趋于稳定,但仍需依据具体数据而定,为保障预估结果的可靠性且不会影响计算效率,本研究设置n=1 000;②树节点抽选的变量个数(m) ,基于n= 1 000,通过选取由不同的m值对模型进行调优,通常取值为所有变量数的1/3,本研究中该参数设置为2。

3结果

不添加竞争指标时,模型的预测能力从大到小表现为:混合效应模型>OLS模型>增强回归树>随机森林;添加竞争指标时,最优模型表现为:混合效应模型>OIS模型>随机森林>增强回归树。OLS模型中添加大于对象木的断面积之和竞争指标时预测能力最佳;增强回归树中添加固定半径为7 m 的简单竞争指数时预测能力最佳;随机森林中添加林分每公顷断面积竞争指标时预测能力最佳;混合效应模型不添加竞争指标时预测能力最佳(RMSE为0.846 0,RMA 为0.211 1,MAE为0.650 1) ,并且都优于其他模型。

3.1 基础模型的确定

利用刀切法进行建模和检验,得出模型最终评价指标值的平均值和标准差(表3)。

▼表3 基础模型的评价指标值

选择的11个基础模型在描述冠幅-胸径关系时, M8、M10和M11拟合不收敛,表3中已剔除。拟合收敛模型中,线性函数(M1)和二次项函数(M2)的拟合效果与检验效果都较其他模型略好，从检验结果看M2的RMSE值比M1的降低0.24% ,RMA值相同,MAE值降低0.66%。因此选择M2作为基础模型。

3.2 竞争指标对冠幅－胸径模型的影响

3.2.1 对OLS模型的影响

将竞争指标分别添加到基础模型中进行冠幅模型拟合和检验,最终选出每个竞争指标的最优模型进行比较。拟合结果见表4,结果表明,每个竞争指标的添加都至少能保留一个参数,说明与冠幅存在竞争关系,且模型的拟合精度有所不同。ABL明显提高了模型精度,模型检验的RMSE、RMA、MAE值较未添加竞争指标分别降低了2.03% 、1.11%和1.85% ;AB和I SD添加则略微降低了模型精度,但不明显;而添加Hegyi的模型而言,模型精度相较不添加任何竞争指标模型略低,且随着Hegyi值的逐渐增加,几乎没有影响模型精度的变化(图1)。固定半径等于1 m时,模型的精度相对较高，与未添加竞争指标模型相比,模型检验的 RMSE、RMA、MAE值分别提高了0.27%、0.32%和0.21%。总的来说,添加ABL时OLS模型的精度最高,其表达式如下:

▼表4 含竞争指标的４种模型的评价指标

注:AB为林分每公顷断面积;ABL为胸径大于对象木的树木断面积之和;ISD为林分密度指数;Hegyi为简单竞争指标。Hegyi(R=*)为固定半径法得到的最优模型,Hegyi(N=*)为固定株数法得到的最优模型。

▲图1 添加Ｈegyi的４种模型评价指标值

式中:W Cki和DBHki分别是k样地内对象木i的冠幅和胸径值; a0、b0、c0和c1为模型参数,εki为模型误差。

3.2.2 对混合效应模型的影响

以M2为基础模型构建冠幅-胸径混合效应模型。拟合结果(表4)表明,添加ABL时,模型精度略微降低,添加AB和I SD模型精度则出现明显降低，与未添加竞争指标相比RMSE值分别提高了5.77%和5.87% , RMA值分别提高了3.74%和3.69% ,MAE值分别提高了6.28%和6.41% ;而就添加Hegyi的模型而言,随着固定半径和固定株数的增加,模型精度相较不添加任何竞争指标模型略低,且模型间精度没有明显变化(图1)。总的来说,不添加任何竞争指标的混合效应模型使得冠幅-胸径模型精度最高,其表达式如下:

式中: ui1、ui2和ui3为随机参数。

3.2.3 对BRT的影响

在BRT基础上添加竞争指标,拟合结果(表4)表明,添加竞争指标AB、ABL和I SD,模型的解释能力都未得到提高,与未添加竞争指标相比,模型预测的RMSE值提高了1.60%~3.54% ,RMA值提高了3.88%~ 12.74%，MAE值提高了1.62%~6.20% ;而就添加Hegyi的模型而言(图1),随着固定半径的增加,模型精度呈波动变化,但模型精度均比不添加竞争指标高,添加固定半径为7 m时，模型精度最高,与未添加竞争指标相比RMSE值降低了2.80% , RMA值降低了2.02% , MAE值降低了1.51%。随着固定株数的增加,模型精度也呈波动变化,添加固定株数在4株以内的竞争指标,模型精度比不添加竞争指标模型低,当添加固定株数为4株时,模型精度较高,与未添加竞争指标模型相比,RMSE值降低了1.71% , RMA值降低了1.57% ,MAE值降低了0.96% ,总的来说,添加固定半径为7m的Hegyi时BRT模型的精度最高。

3.2.4 对RF的影响

在RF基础上添加竞争指标,模型的解释能力都有得到提高。拟合结果(表4)表明,就添加竞争指标AB、ABL和I SD而言,模型的精度表现为:AB>ABL>I SD,AB的添加明显提高了模型精度,与未添加指标模型相比,模型检验的RMSE、RMA、MAE值分别降低了10.65% 、5.88%和7.25%;而就添加Hegyi的模型而言(图1),模型精度随着Hegyi值的逐渐增加,有明显的波动变化,且最优的模型均为添加AB的模型精度高。

3.3 模型精度比较

通过对比各个模型在检验数据集上的表现可知,4种模型对冠幅的预测能力均很好,但结果存在一定的差异。不添加竞争指标时,4种模型的预测能力表现为:混合效应模型>OIS模型>BRT>RF。从检验数据来看,混合效应模型的RMSE、RMA和 MAE值要明显低于其他3种模型,RMSE值分别降低了32.62% 、43.68%和46.03% , RMA值分别降低了24.53% 、32.23%和37.95% ,MAE值分别降低了31.16% ,39.11%和43.02%。可见，不添加竞争指标时,混合效应模型的冠幅-胸径模型的精度是最高的;在添加竞争指标后,除混合混合效应模型外,其余3种模型方法的最优模型精度都得到提升。将添加竞争指标后4种建模方法的最优模型进行比较,模型预测能力发生变化:混合效应模型>OLS模型>RF>BRT,混合效应模型的精度依然明显高于其他模型,同时,RF在添加AB后，模型精度高过BRT的最优模型。但无论是是否添加竞争指标,OILS模型和混合效应模型的精度都明显高于BRT和 RF。将全部数据代入OIS最优模型和混合效应最优模型中,计算固定参数值、随机参数协方差组成以及模型拟合统计量,结果见表5,混合效应模型具有较高的精度,可以很好地模拟闽楠的冠幅和胸径的关系。

▼表5 模型参数、方差估计值及评价指标

4结论

研究结果对天然闽楠单木冠幅生长进行精准预测,以及为提高闽楠天然次生林的林分质量提供理论依据。

关注我们，更多精彩

,