欢迎访问林业科学,今天是

林业科学 ›› 2022, Vol. 58 ›› Issue (8): 89-98.doi: 10.11707/j.1001-7488.20220809

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

基于最粗优势木胸径生长的湖南栎类天然林立地质量评价模型

何静1,李新建2,朱晋梅3,朱光玉1,*   

  1. 1. 中南林业科技大学林学院 长沙 410004
    2. 国家林业和草原局中南调查规划院 长沙 410004
    3. 湖南省林业局 长沙 410001
  • 收稿日期:2021-06-10 出版日期:2022-08-25 发布日期:2022-12-19
  • 通讯作者: 朱光玉
  • 基金资助:
    基于林分结构效应的亚热带栎类天然混交林全林分生长模型(32271874);亚热带栎类天然混交林立地质量评价与生长预估(31570631)

Site Quality Evaluation Model of Natural Quercus Forests in Hunan Based on the Growth of the Thickest Dominant Tree Diameter at Breast Height

Jing He1,Xinjian Li2,Jinmei Zhu3,Guangyu Zhu1,*   

  1. 1. Forestry College, Central South University of Forestry and Technology Changsha 410004
    2. Central South Inventory and Planning Institute, National Forestry and Grassland Administration Changsha 410004
    3. The Forestry Department of Hunan Province Changsha 410001
  • Received:2021-06-10 Online:2022-08-25 Published:2022-12-19
  • Contact: Guangyu Zhu

摘要:

目的: 分析立地因子对林分最粗优势木胸径生长的影响, 构建含立地类型混合效应的栎类天然林优势木胸径生长模型, 导出以最粗优势木胸径为指标的基于立地分级的立地质量评价模型, 为栎类天然林立地质量评价提供一种新方法。方法: 以湖南栎类天然林为研究对象, 基于51块样地实测数据, 采用数量化方法Ⅰ筛选对优势木胸径生长影响显著的立地因子, 将立地因子按照标准分级、组合, 构成初始立地类型; 通过R语言拟合栎类天然林优势木胸径与年龄的相关关系, 筛选最优基础模型, 将初始立地类型作为随机效应加入基础模型构建混合效应模型; 应用k-means聚类将影响效果相近或相同的初始立地类型聚类成立地类型组, 并将其作为随机效应加入最优基础模型构建含立地类型组的混合效应模型; 通过导算, 得到立地质量评价模型, 采用方差分析验证林分断面积与立地指数的显著关系。结果: 对优势木胸径生长影响显著的立地因子包括海拔、坡度、坡位和坡向, 显著性顺序为海拔>坡度>坡向>坡位; 选取4种常见的树木理论生长方程进行拟合, 模型确定系数(R2)均在0.7左右, 其中Richards模型的拟合效果最好, R2为0.731 8, 平均绝对误差(MAE)为5.442 6, 均方根误差(RMSE)为6.879 1, 表达式为D=a×[(1-exp(-c×AGE)]^b; 将筛选的4种显著性立地因子按照标准分级、组合构成初始立地类型, 加入基础模型构建含立地类型的混合效应模型, R2升至0.901 6; 应用k-means聚类将初始立地类型聚类成6个立地类型组, 作为随机效应加入基础模型, 最优模型表达式为Dj=aj×[1-exp(-c×AGE)]^b+ε, R2为0.926 9, 相比基础模型提高26.7%, 赤池信息量(AIC)和贝叶斯信息量(BIC)均有所减小。采用含立地类型组的混合效应模型导出最粗优势木胸径为指标的基于立地分级的湖南栎类天然林立地质量评价模型为SQEIM-DBH=aj×[1-exp(-0.03×AGE0])^$\left\{\frac{\ln D_j}{\ln a_j \times[1-\exp (-0.03 \times \mathrm{AGE})]}\right\}$+ε, 经林分断面积验证, 立地指数SQEIM-DBH与林分断面积显著相关。结论: 立地因子对湖南栎类天然林优势木胸径生长影响显著, 以最粗优势木胸径为指标评价栎类天然林立地质量、预估林地生产力在理论上具有可行性。

关键词: 栎类天然林, 最粗优势木, 胸径生长, 混合效应, 立地质量评价模型

Abstract:

Objective: This study was carried out to analyze the effects of site factors on the DBH growth of the thickest dominant tree in the stand, construct a growth model of the dominant tree in Hunan Quercus natural forests with the mixed effects of site types and derive a site quality evaluation model based on the DBH of the thickest dominant tree as an indicator. The site quality evaluation equation was expected to propose a new method for the site quality evaluation of Quercus natural forests. Method: This paper took Hunan Quercus natural forests as the research objects. Based on the measured data of 51 plots, the site factors those having significant impacts on the growth of the dominant woods at breast height were selected by the quantitative method Ⅰ, and the site factors were classified and combined according to the standard. The initial site type, through the R language, the DBH-age of the dominant tree in the natural Quercus forests in Hunan was fitted, the optimal basic model was screened, and the initial site type was added as a random effect to the basic model to construct a mixed effect model. The k-means was used to gather the class clusters the initial site types with similar effects to form a site type group, and it was added as a random effect to the optimal basic model to construct a mixed-effect model containing the site type group. By derivative, the site quality evaluation model was obtained, and the analysis of variance was used to verify the significant relationships between forest stand area and site index. Result: The site factors those having significant effects on the growth of the thickest dominant tree diameter at breast height were altitude, slope gradient, slope position, and slope aspect. The order of significance was altitude>slope gradient>slope aspect>slope position. Four theoretical growth equations were selected for fitting, and the coefficient of determination of the model was about 0.7. The optimal basic model was the Richards model, and the coefficient of determination R2 was 0.731 8, the mean absolute error (MAE) was 5.442 6, the root mean square error (RMSE) was 6.879 1, and the expression equation was D=a×[1-exp(-c×AGE)]^b.By combining the selected four significant site factors to form the initial site type and adding it to the basic model to construct a mixed effect model with site types, the determination coefficient (R2) was increased to 0.9016. The k-means clustering was used to cluster the initial site types into 6 site type groups, which were added to the basic model as a random effect, the optimal model expression equation was Dj=aj×[1-exp(-c×AGE)]^b+ε, the model's determination coefficient (R2) was increased to 0.924 4, which was 26.7% higher than that of the basic model, whereas both AIC and BIC were reduced. By plotting a polymorphic site index curve based on a mixed effect model containing site type groups, the site quality evaluation equation of Hunan Quercus natural forests, based on site classification with the dominant wood diameter at breast height as the evaluation index, was derived as SQEIM-DBH=aj×[1-exp(-0.03×AGE0)]^$\left\{\frac{\ln D_j}{\ln a_j \times[1-\exp (-0.03 \times \mathrm{AGE})]}\right\}$+ε, and through the stand-off area verification, the site index SQEIM-DBH obtained in this study was significantly related to the stand-off area. Conclusion: Site factors might have a significant impact on the growth of the diameter at breast height of the Quercus in Hunan Province. It could be theoretically feasible to evaluate the site quality of natural Quercus forests and to predict the productivity of the natural forests based on the diameter of the thickest dominant wood as an indicator, which was expected to provide a way for the evaluation of site quality of natural forests.

Key words: Quercus natural forest, the thickest dominant tree, DBH growth, mixed effect, site quality evaluation model

中图分类号: