欢迎访问《应用生态学报》官方网站,今天是

应用生态学报 ›› 2026, Vol. 37 ›› Issue (3): 731-740.doi: 10.13287/j.1001-9332.202603.002

• 研究论文 • 上一篇    下一篇

基于广义加性混合模型的长白落叶松人工林树高-胸径模型

李志超, 梅雪松, 董灵波*   

  1. 东北林业大学林学院, 森林生态系统可持续经营教育部重点实验室, 哈尔滨 150040
  • 收稿日期:2025-11-23 接受日期:2026-01-10 出版日期:2026-03-18 发布日期:2026-09-18
  • 通讯作者: *E-mail: farrell0503@126.com
  • 作者简介:李志超, 男, 1998年生, 硕士研究生。主要从事森林培育与经营研究。E-mail: 1433766912@qq.com
  • 基金资助:
    国家自然科学基金项目(32471856)

Tree height-diameter modeling of Larix olgensis plantations based on generalized additive mixed models

LI Zhichao, MEI Xuesong, DONG Lingbo*   

  1. Ministry of Education Key Laboratory of Sustainable Management of Forest Ecosystems, College of Forestry, Northeast Forestry University, Harbin 150040, China
  • Received:2025-11-23 Accepted:2026-01-10 Online:2026-03-18 Published:2026-09-18

摘要: 本研究基于东北林业大学帽儿山实验林场和哈尔滨市胜利实验林场71块固定样地数据,构建仅包含林分密度指数(SDI)或立地指数(SI)的主效应广义加性模型(GAM)、包含二者交互作用的广义加性混合模型(GAMM),以及引入胸径对数(logD)调整斜率的变系数模型(VC-GAMM)。采用赤池信息准则(AIC)和贝叶斯信息准则(BIC)进行模型筛选,利用调整决定系数(Radj2)、均方根误差百分比(RMSE)和平均绝对误差百分比(MAE)评价模型性能,量化林分密度和立地条件及其交互作用对长白落叶松人工林单木树高-胸径曲线的影响,为该地区人工林的可持续经营提供理论依据。结果表明:VC-GAMM模型为最优树高-胸径模型,其Radj2值为0.733,较GAM模型提升9.6%~13.6%,较GAMM模型提升0.14%,RMSE和MAE分别为12.4%和9.7%,与其他模型相比处于最低水平;边际效应结果表明,在立地条件极差(标准化立地指数SIs=-1.89)时,标准化林分密度指数(SDIs)每增加1个标准差仅可提升树高0.34 m,而在极好立地条件(SIs =+2.15)下可提高0.78 m,树高增益相差1.9%;树高最大值分布于中等SDIs和较高SIs的组合区域;林分密度和立地条件的交互作用在小径级阶段(胸径25%分位数)最强,其树高相对变异率高出中、大径级阶段5%。林分密度和立地条件间存在显著的交互作用,其解释树高变异的能力优于单一因子;在良好立地条件下增加林分密度可促进树高生长,而贫瘠立地条件下,高林分密度会加剧竞争并抑制树高生长;树高对林分密度和立地条件交互作用的敏感性在林分早期阶段更强。

关键词: 林分密度, 立地条件, 交互作用, 树高-胸径模型, 广义加性模型

Abstract: Based on data from 71 permanent plots established in the Maoershan Experimental Forest Farm of Northeast Forestry University and the Shengli Experimental Forest Farm in Harbin City, we developed three types of height-diameter models for Larix olgensis plantations. A generalized additive model (GAM) included only the main effects of either the stand density index (SDI) or the site index (SI). A generalized additive mixed model (GAMM) incorporated the interaction between SDI and SI. A varying-coefficient generalized additive mixed model (VC-GAMM) adjusted the slope by the logarithm of diameter at breast height (logD). Model selection was based on the Akaike information criterion (AIC) and Bayesian information criterion (BIC). Model performance was evaluated using the adjusted coefficient of determination (Radj2), root mean square error percentage (RMSE), and mean absolute error percentage (MAE). The effects of stand density, site conditions, and their interaction on individual-tree height-diameter curves of L. olgensis plantations were quantified to provide a theoretical basis for sustainable plantation management. The results showed that the VC-GAMM was the optimal height-diameter model, with a Radj2 value of 0.733, representing an improvement of 9.6%-13.6% over the GAM models and 0.14% over the GAMM model, while its RMSE and MAE were 12.4% and 9.7%, respectively, which were the lowest among all kinds of models. Marginal effect analysis showed that under extremely poor site conditions (standardized site index SIs=-1.89), a one-standard-deviation increase in the standardized stand density index (SDIs) increased tree height by only 0.34 m, whereas under excellent site conditions (SIs=+2.15), tree height increased by 0.78 m, corresponding to 1.9% difference in relative height gain. The maximum tree height occurred in regions characterized by moderate SDIs and relatively high SIs. The interaction between stand density and site conditions was strongest at the small-diameter stage (25th percentile of DBH), where the relative variation in tree height was 5% higher than that at the medium- and large-diameter stages. There was a significant interaction between stand density and site conditions, and their combined explanatory power for tree height variation exceeded that of either factor alone. Increasing stand density promoted height growth under favorable site conditions. Under poor site conditions, high stand density intensified competition and suppressed height growth, with tree height being more sensitive to this interaction during the early developmental stages of the stand.

Key words: stand density, site quality, interaction effect, tree height-diameter model, generalized additive model