基于非线性混合效应模型的长白落叶松人工林冠长率高效预估模型

刘帅东, 苗铮, 郝元朔, 董利虎

南京林业大学学报(自然科学版) ›› 2026, Vol. 50 ›› Issue (5) : 152-160.

PDF(1986 KB)
PDF(1986 KB)
南京林业大学学报(自然科学版) ›› 2026, Vol. 50 ›› Issue (5) : 152-160. DOI: 10.12302/j.issn.1000-2006.202509032
研究论文

基于非线性混合效应模型的长白落叶松人工林冠长率高效预估模型

作者信息 +

Efficient prediction model for crown ratio of Larix olgensis plantations based on nonlinear mixed-effects model

Author information +
文章历史 +

摘要

【目的】冠长率是反映林木活力、支撑林分生长预测的关键指标。然而,适用于长白落叶松(Larix olgensis )的高精度冠长率预测模型仍较为缺乏,其主要影响因子尚未明晰;且现有模型多依赖测量成本较高的树高与枝下高等变量。因此,本研究旨在采用易于获取的变量,构建长白落叶松人工林的非线性混合效应(NLME)冠长率模型,并确定高效的预测校正方案,以服务于精准森林经营。【方法】基于黑龙江省东京城林业局、林口林业局、孟家岗林场共463块长白落叶松标准地观测数据,以理查德(Richards)模型为基础,采用逐步回归法引入林分优势木平均高、林分公顷断面积以及与距离相关的竞争指数(Ic5)作为协变量,构建非线性混合效应模型。采用留一交叉检验和随机抽样校正方案,分析不同校正样本量(1~20株)对模型预测精度的影响。【结果】非线性混合效应模型调整决定系数(${\mathit{R}}_{\mathrm{a}}^{2}$)为0.705 4,均方根误差(RMSE)为0.086 9,赤池信息准则(AIC)为-17 147.76,表明模型拟合优度良好。参数估计显示,冠长率与林分优势木平均高、林分公顷断面积呈负相关,而与竞争指数(Ic5)呈正相关。抽样校正可显著提升预测精度,但当样本量超过5株时,精度提升效果显著减弱。【结论】本研究构建的非线性混合效应模型,在不依赖树高相关测量的前提下,有效量化了立地质量、林分密度与竞争对长白落叶松冠长率的多因素影响。推荐每个样地抽取5株样木作为最优校正方案,为黑龙江省长白落叶松人工林冠长率的低成本、高精度预测提供可靠方法与实践依据。

Abstract

【Objective】Crown ratio (CR) is a vital indicator of tree vigor and a fundamental parameter for stand growth prediction. However, high-precision CR prediction models specifically for Larix olgensis plantations remain scarce, and the key driving factors are not yet fully understood. Furthermore, existing models often rely on variables such as total tree height and height to crown base, which are costly and difficult to measure. Therefore, this study aimed to develop a nonlinear mixed-effects (NLME) crown ratio model for L. olgensis plantations using easily obtainable variables and to determine an efficient calibration strategy for prediction to support precision forest management. 【Method】The study utilized observation data from 463 sample plots of L. olgensis in the Dongjingcheng Forestry Bureau, Linkou Forestry Bureau, and Mengjiagang Forest Farm in Heilongjiang Province. A Richards function was selected as the basic model. Stepwise regression was employed to incorporate stand mean dominant height, stand basal area per hectare, and a distance-dependent competition index (Ic5) as covariates to construct the NLME model. Leave-one-out cross-validation and random sampling calibration schemes were employed to analyze the effects of different calibration sample sizes (ranging from 1 to 20 trees) on model prediction accuracy. 【Result】The developed NLME model demonstrated good goodness-of-fit, with an adjusted coefficient of determination (${\mathit{R}}_{\mathrm{a}}^{2}$) of 0.705 4, a root mean square error (RMSE) of 0.086 9, and an Akaike information criterion (AIC) of -17 147.76. Parameter estimates indicated that CR was negatively correlated with stand mean dominant height and stand basal area per hectare, but positively correlated with competition index (Ic5). Sampling calibration significantly improved prediction accuracy; however, the marginal gain in accuracy diminished markedly when the sample size exceeded five trees. 【Conclusion】The NLME model developed in this study effectively quantifies the multifactorial effects of site quality, stand density, and competition on the CR of Larix olgensis without relying on height-related measurements. A calibration sample size of five trees per plot is recommended as the optimal strategy. These findings provide a reliable method and practical tool for achieving low-cost, high-accuracy CR predictions in L. olgensis plantations in Heilongjiang Province.

关键词

冠长率 / 长白落叶松 / 非线性混合效应模型 / 抽样校正

Key words

crown ratio / Larix olgensis / nonlinear mixed-effects model / sampling calibration

引用本文

导出引用
刘帅东, 苗铮, 郝元朔, 等. 基于非线性混合效应模型的长白落叶松人工林冠长率高效预估模型[J]. 南京林业大学学报(自然科学版). 2026, 50(5): 152-160 https://doi.org/10.12302/j.issn.1000-2006.202509032
Liu Shuaidong, Miao Zheng, Hao Yuanshuo, et al. Efficient prediction model for crown ratio of Larix olgensis plantations based on nonlinear mixed-effects model[J]. Journal of Nanjing Forestry University (Natural Sciences Edition). 2026, 50(5): 152-160 https://doi.org/10.12302/j.issn.1000-2006.202509032
中图分类号: S758.5   

参考文献

[1]
郭孝玉. 长白落叶松人工林树冠结构及生长模型研究[D]. 北京: 北京林业大学, 2013.
Guo X Y. Crown structure and growth model for Larix olgensis plantation[D]. Beijing: Beijing Forestry University, 2013.
[2]
Silva B, Álava-Núñez P, Strobl S, et al. Area-wide evapotranspiration monitoring at the crown level of a tropical mountain rain forest[J]. Remote Sensing of Environment, 2017, 194:219-229. DOI:10.1016/j.rse.2017.03.023.
[3]
Hernandez-Moreno J M, Bayeur N M, Coley H D, et al. Clouds homogenize shoot temperatures,transpiration,and photosynthesis within crowns of Abies fraseri (Pursh.) Poiret.[J]. Oecologia, 2017, 183(3):667-676. DOI:10.1007/s00442-016-3799-7.
[4]
Zarnoch S J, Bechtold W A, Stolte K W. Using crown condition variables as indicators of forest health[J]. Canadian Journal of Forest Research, 2004, 34(5):1057-1070. DOI:10.1139/x03-277.
[5]
Fu L Y, Sharma R P, Hao K J, et al. A generalized interregional nonlinear mixed-effects crown width model for Prince Rupprecht larch in northern China[J]. Forest Ecology and Management, 2017, 389:364-373. DOI:10.1016/j.foreco.2016.12.034.
[6]
Fu L Y, Zhang H R, Lu J, et al. Multilevel nonlinear mixed-effect crown ratio models for individual trees of Mongolian oak (Quercus mongolica) in Northeast China[J]. PLoS One, 2015, 10(8):e0133294. DOI:10.1371/journal.pone.0133294.
[7]
Zhao D H, Kane M, Borders B E. Crown ratio and relative spacing relationships for loblolly pine plantations[J]. Open Journal of Forestry, 2012, 2(3):101-115. DOI:10.4236/ojf.2012.23014.
[8]
Lessard V C, McRoberts R E, Holdaway M R. Diameter growth models using Minnesota forest inventory and analysis data[J]. Forest Science, 2001, 47(3):301-310. DOI:10.1093/forestscience/47.3.301.
[9]
Jiang L C, Liu R L. Segmented taper equations with crown ratio and stand density for Dahurian larch (Larix gmelinii) in Northeastern China[J]. Journal of Forestry Research, 2011, 22(3):347-352. DOI:10.1007/s11676-011-0178-4.
[10]
Leites L P, Robinson A P, Crookston N L. Accuracy and equivalence testing of crown ratio models and assessment of their impact on diameter growth and basal area increment predictions of two variants of the Forest Vegetation Simulator[J]. Canadian Journal of Forest Research, 2009, 39(3):655-665. DOI:10.1139/x08-205.
[11]
Dickinson Y L, Battaglia M A, Asherin L A. Evaluation of the FVS-CR diameter growth model in structurally-heterogeneous ponderosa pine (Pinus ponderosa Douglas ex C.Lawson) stands in the Southern Rockies,and potential modifications[J]. Forest Ecology and Management, 2019, 448:1-10. DOI:10.1016/j.foreco.2019.05.031.
[12]
Temesgen H, LeMay V, Mitchell S J. Tree crown ratio models for multi-species and multi-layered stands of southeastern British Columbia[J]. The Forestry Chronicle, 2005, 81(1):133-141. DOI:10.5558/tfc81133-1.
[13]
Wykoff W, Crookston N L, Stage A R. User’s guide to the stand prognosis model[M]. Ogden,Utah: U.S.Dept.of Agriculture,Forest Service,Intermountain Forest and Range Experiment Station, 1982. DOI:10.5962/bhl.title.109367.
[14]
Hasenauer H, Monserud R A. A crown ratio model for Austrian forests[J]. Forest Ecology and Management, 1996, 84(1/2/3):49-60. DOI:10.1016/0378-1127(96)03768-1.
[15]
卢军, 李凤日, 张会儒, 等. 帽儿山天然次生林主要树种冠长率模型[J]. 林业科学, 2011, 47(6):70-76.
Lu J, Li F R, Zhang H R, et al. A crown ratio model for dominant species in secondary forests in Mao’er mountain[J]. Scientia Silvae Sinicae, 2011, 47(6):70-76.
[16]
Soares P, Tomé M. A tree crown ratio prediction equation for eucalypt plantations[J]. Annals of Forest Science, 2001, 58(2):193-202. DOI:10.1051/forest:2001118.
[17]
刘奇峰, 陈东升, 冯健, 等. 辽东山区人工红松冠长率影响因子研究[J]. 林业科学研究, 2022, 35(6):127-134.
Liu Q F, Chen D S, Feng J, et al. Influencing factors on the crown ratio of Korean pine plantation in eastern Liaoning Province[J]. Forest Research, 2022, 35(6):127-134. DOI:10.13275/j.cnki.lykxyj.2022.006.014.
[18]
Popoola F S, Adesoye P O. Crown ratio models for Tectona grandis (Linn.f) stands in osho forest reserve,oyo state,Nigeria[J]. Journal of Forest and Environmental Science, 2012, 28(2):63-67. DOI:10.7747/jfs.2012.28.2.063.
[19]
Miao Z, Widagdo F R A, Dong L H, et al. Prediction of branch growth using quantile regression and mixed-effects models:an example with planted Larix olgensis Henry trees in Northeast China[J]. Forest Ecology and Management, 2021, 496:119407. DOI:10.1016/j.foreco.2021.119407.
[20]
易达, 李凤日, 马爱云, 等. 基于混合效应模型和分位数回归的长白落叶松枝下高模型构建[J]. 应用生态学报, 2023, 34(4):1035-1042.
Yi D, Li F R, Ma A Y, et al. Model construction for height to crown base of Larix olgensis based on mixed-effects model and quantile regression[J]. Chinese Journal of Applied Ecology, 2023, 34(4):1035-1042. DOI:10.13287/j.1001-9332.202304.007.
[21]
马爱云, 李凤日, 董利虎. 应用混合效应法建立的杂种落叶松人工林单木冠幅预测模型[J]. 东北林业大学学报, 2021, 49(3):9-14,38.
Ma A Y, Li F R, Dong L H. Individual crown diameter prediction for hybrid larch based on mixed effects models[J]. Journal of Northeast Forestry University, 2021, 49(3):9-14. DOI:10.13759/j.cnki.dlxb.2021.03.002.
[22]
Monserud R A, Sterba H. Modeling individual tree mortality for Austrian forest species[J]. Forest Ecology and Management, 1999, 113(2/3):109-123. DOI:10.1016/S0378-1127(98)00419-8.
[23]
宋争, 黄朗, 胡松, 等. 基于立地混合效应的湖南丘陵平原区杉木多形立地指数模型研究[J]. 林业科学研究, 2022, 35(5):172-179.
Song Z, Huang L, Hu S, et al. Polymorphic site-index model with site mixed effects for Chinese fir plantations in Hunan hilly Plains[J]. Forest Research, 2022, 35(5):172-179. DOI:10.13275/j.cnki.lykxyj.2022.005.019.
[24]
袁慧, 杜超群, 赵小涛, 等. 基于混合效应的湖北杉木冠长率模型[J]. 中国农学通报, 2021, 37(4):31-37.
Yuan H, Du C Q, Zhao X T, et al. A crown ratio model for Cunninghamia lanceolata plantations based on mixed effect in Hubei[J]. Chinese Agricultural Science Bulletin, 2021, 37(4):31-37.
[25]
Sharma R P, Vacek Z, Vacek S. Generalized nonlinear mixed-effects individual tree crown ratio models for Norway spruce and European beech[J]. Forests, 2018, 9(9):555. DOI:10.3390/f9090555.
[26]
Dyer M E, Burkhart H E. Compatible crown ratio and crown height models[J]. Canadian Journal of Forest Research, 1987, 17(6):572-574. DOI:10.1139/x87-096.
[27]
Poudel K P, Avery S C, Granger J J. Live crown ratio models for loblolly pine (Pinus taeda) with beta regression[J]. Forests, 2021, 12(10):1409. DOI:10.3390/f12101409.
[28]
Dong L H, Pukkala T, Li F R, et al. Developing distance-dependent growth models from irregularly measured sample plot data-a case for Larix olgensis in Northeast China[J]. Forest Ecology and Management, 2021, 486:118965. DOI:10.1016/j.foreco.2021.118965.
[29]
刘兆刚, 李凤日, 于金成. 落叶松人工林单木模型的研究[J]. 植物研究, 2003, 23(2):237-244.
Liu Z G, Li F R, Yu J C. The study of individual tree model on Larix olgensis plantation[J]. Bulletin of Botanical Research, 2003, 23(2):237-244. DOI:10.3969/j.issn.1673-5102.2003.02.028.
[30]
Hegyi F. A simulation model for managing jack pine stands[C]// Fries J. Growth Models for Tree and Stand Simulation. Stockholm: Royal College of Forestry, 1974:74-90.
[31]
江希钿, 陈学文, 林纪建. 以Von Bertalanffy生长理论为基础的单木生长模型[J]. 中南林业调查规划, 1994, 13(4):5-7.
Jiang X D, Chen X W, Lin J J, et al. Single tree growth model based on Von Bertalanffy growth theory[J]. Central South Forest Inventory and Planning, 1994, 13(4):5-7.
[32]
Yan Y F, Wang J J, Liu S M, et al. Modeling the influence of competition,climate,soil,and their interaction on height to crown base for Korean pine plantations in Northeast China[J]. European Journal of Forest Research, 2024, 143:1627-1640. DOI:10.1007/s10342-024-01710-z.
[33]
Fu L Y, Sun H, Sharma R P, et al. Nonlinear mixed-effects crown width models for individual trees of Chinese fir (Cunninghamia lanceolata) in south-central China[J]. Forest Ecology and Management, 2013, 302:210-220. DOI:10.1016/j.foreco.2013.03.036.
[34]
孙宇, 李凤日, 谢龙飞, 等. 基于林分及地形因子的落叶松人工林林分生物量模型构建[J]. 南京林业大学学报(自然科学版), 2023, 47(3):129-136.
Sun Y, Li F R, Xie L F, et al. Construction of the stand-level biomass model of Larix olgensis plantations based on stand and topographic factors[J]. Journal of Nanjing Forestry University (Natural Sciences Edition), 2023, 47(3):129-136.
[35]
Xie L F, Widagdo F R A, Dong L H, et al. Modeling height-diameter relationships for mixed-species plantations of Fraxinus mandshurica rupr.and Larix olgensis Henry in Northeastern China[J]. Forests, 2020, 11(6):610. DOI:10.3390/f11060610.
[36]
Yu Y, Zhou Z Y, Sharma R P, et al. Comparing crown ratio models for spruce-fir broadleaved mixed forests using beta regression and random forest algorithm[J]. Computers and Electronics in Agriculture, 2024, 225:109302. DOI:10.1016/j.compag.2024.109302.
[37]
Zhang X Q, Wang H C, Chhin S, et al. Effects of competition,age and climate on tree slenderness of Chinese fir plantations in southern China[J]. Forest Ecology and Management, 2020, 458:117815. DOI:10.1016/j.foreco.2019.117815.
[38]
Yang Y Q, Titus S J, Huang S. Modeling individual tree mortality for white spruce in Alberta[J]. Ecological Modelling, 2003, 163(3):209-222. DOI:10.1016/S0304-3800(03)00008-5.
[39]
Kramer R D, Sillett S C, Van Pelt R, et al. Neighborhood competition mediates crown development of Picea sitchensis in Olympic rainforests:implications for restoration management[J]. Forest Ecology and Management, 2019, 441:127-143. DOI:10.1016/j.foreco.2019.03.027.

基金

国家重点研发计划青年科学家项目(2022YFD2201800)
国家自然科学基金青年项目(32301581)
第五批林草科技创新青年拔尖人才项目(2025132010)

责任编辑: 李燕文
PDF(1986 KB)

Accesses

Citation

Detail

段落导航
相关文章

/

〈 〉