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

应用生态学报 ›› 2026, Vol. 37 ›› Issue (4): 1153-1164.doi: 10.13287/j.1001-9332.202604.009

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

基于SGAM模型的小兴安岭地区林下植被碳密度空间分布格局

刘宇蒙, 贾炜玮*, 赵子鹏, 李泽霖   

  1. 东北林业大学林学院, 森林生态系统可持续经营教育部重点实验室, 哈尔滨 150040
  • 收稿日期:2026-01-06 修回日期:2026-03-03 出版日期:2026-04-18 发布日期:2026-05-29
  • 通讯作者: *E-mail: jiaww2002@163.com
  • 作者简介:刘宇蒙, 女, 2000年生, 硕士研究生。主要从事森林可持续经营研究。E-mail: 71735622@qq.com
  • 基金资助:
    十四五重点研发计划项目子课题(2022YFD2201003-02)、哈尔滨市林业和草原局直属林场森林修复技术咨询项目(HFW240100013)、伊春国土绿化示范项目和碳汇成效监测评估和森林碳汇研究项目[CS]20250026

Assessing spatial pattern of carbon density of understory vegetation in Xiaoxing’anling area based on SGAM model

LIU Yumeng, JIA Weiwei*, ZHAO Zipeng, LI Zelin   

  1. Ministry of Education Key Laboratory of Sustainable Management of Forest Ecosystem, College of Forestry, Northeast Forestry University, Harbin 150040, China
  • Received:2026-01-06 Revised:2026-03-03 Online:2026-04-18 Published:2026-05-29

摘要: 本研究以黑龙江省小兴安岭伊春地区为研究区,基于2021年1194块天然林和人工林样地实测数据,估算幼树层、灌木层和草本层的碳密度,结合地形、气候及林分因子,采用多元逐步回归模型(SMR)、广义加性模型(GAM)和空间广义加性模型(SGAM)对不同层次林下植被碳密度进行拟合,结合空间插值分析伊春地区林下植被的区域尺度分异特征。结果表明:研究区林下植被平均碳密度为0.418 t·hm-2,其中幼树层、灌木层和草本层分别为0.223、0.172和0.023 t·hm-2,幼树层和灌木层是林下植被碳密度的主要贡献层次,合计贡献94.3%,草本层贡献占比为5.7%;林分、盖度、郁闭度、年均降水量和海拔是影响林下植被碳密度的主要因子。模型比较结果显示,SGAM在3个层次上的拟合效果均优于SMR和GAM,3个模型的决定系数分别为0.79、0.71和0.62。空间插值验证结果显示,林下植被碳密度整体呈现北高南低、东多西少的空间分布格局,经验贝叶斯克里金插值在幼树层、灌木层和草本层中的均方根误差分别为5.47、3.57和1.90 t·hm-2,预测精度优于反距离权重法、径向基函数法及普通克里金法。综合考虑非线性关系与空间效应的SGAM模型及经验贝叶斯克里金插值方法,本研究揭示了林下植被碳密度驱动因素及其空间分布格局,可为森林碳汇评估和林下植被经营管理提供科学依据。

关键词: 林下植被, 碳密度, 空间广义加性模型, 空间插值

Abstract: We estimated carbon densities of the sapling, shrub, and herb layers in the Yichun region of the Xiao-xing’anling area in Heilongjiang Province, based on field measurements from 1194 natural and plantation forest plots in 2021. Combined with terrain, climate, and forest factors, we used multiple stepwise regression (SMR), generalized additive model (GAM), and spatial generalized additive model (SGAM) to fit carbon density for understory vegetation across different layers. We then analyzed the regional differentiations of understory vegetation by spatial interpolation analysis. Results showed that the average carbon density of understory vegetation was 0.418 t·hm-2, comprising 0.223, 0.172, and 0.023 t·hm-2 for the sapling, shrub, and herbaceous layers, respectively. The sapling and shrub layers were the primary contributors to understory carbon density, accounting for 94.3% of the total, while herbaceous layer contributed 5.7%. Stand density, canopy cover, closure, annual precipitation, and altitude were the primary factors influencing understory carbon density. Model comparison revealed that the SGAM model outperformed both SMR and GAM across all three strata, with determination coefficients of 0.79, 0.71, and 0.62, respectively. Spatial interpolation validation revealed an overall spatial distribution pattern of understory carbon density being higher in the north and lower in the south, with greater abundance in the east and less in the west. The root mean square error values for empirical Bayesian kriging interpolation were 5.47, 3.57, and 1.90 t·hm-2 for the sapling, shrub, and herb layers, respectively, demonstrating superior prediction accuracy compared to inverse distance weighting, radial basis function, and ordinary kriging methods. By integrating the SGAM model, which accounted for both nonlinear relationships and spatial effects, with empirical Bayesian kriging interpolation, we elucidated the drivers and spatial distribution patterns of carbon density in understory. These fin-dings would provide a scientific basis for forest carbon sink assessments and understory management practices.

Key words: understory vegetation, carbon density, spatial generalized additive model, spatial interpolation