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

林业科学 ›› 2019, Vol. 55 ›› Issue (11): 105-116.doi: 10.11707/j.1001-7488.20191112

• 论文与研究报告 • 上一篇    下一篇

基于地理加权回归模型的大兴安岭中部天然次生林更新分布

张凌宇, 刘兆刚   

  1. 东北林业大学林学院 森林生态系统可持续经营教育部重点实验室 哈尔滨 150040
  • 收稿日期:2019-03-07 修回日期:2019-09-02 出版日期:2019-11-25 发布日期:2019-12-21
  • 基金资助:
    国家重点研发计划项目(2017YFC0504103)。

Regeneration and Distribution of Natural Secondary Forests in the Central Part of Daxing'an Mountains Based on Geographically Weighted Regression Model

Zhang Lingyu, Liu Zhaogang   

  1. Key Laboratory of Sustainable Forest Ecosystem Management, Ministry of Education School of Forestry, Northeast Forestry University Harbin 150040
  • Received:2019-03-07 Revised:2019-09-02 Online:2019-11-25 Published:2019-12-21

摘要: [目的] 研究大兴安岭地区天然次生林不同位置间森林更新的空间相关性和空间分布模式,探索尺度效应对空间自相关性的影响,对森林更新的影响因子进行分析,从更深层次了解更新动态中潜在的规律性,为该地区天然次生林的经营决策提供理论依据和技术支持。[方法] 以黑龙江省大兴安岭新林林业局翠岗林场为研究区,基于2018年7-8月在研究区建立的45块固定样地数据,以林分因子、地形因子、林分空间结构、土壤厚度和物种多样性5方面的9个因子为自变量,建立全局泊松模型和以地理加权回归模型为基础的4种尺度(5、10、15和20 km)地理加权泊松模型(GWPR)对该地区天然次生林更新状况进行模拟,利用全局Moran I和局域Moran I分别对模型残差的全局空间自相关性和空间分布状况进行描述,评价全局模型和各尺度局域模型的拟合效果,对尺度效应下各局域模型之间的差异进行说明,采用5 km尺度局域模型绘制研究区森林更新的空间分布,对研究区森林更新状况进行评价和分析。[结果] 5 km尺度局域模型具有最好的模型残差局域化空间分布效果,可形成不同模型残差聚集分布的理想状态,模型变量的参数估计值产生跨度最大的变化范围,模型稳定性最好,随着空间尺度逐渐增大模型稳定性不断下降,但总体上要好于全局模型;同时,处于该尺度下的局域模型,模型残差的空间自相关性最低。局域模型的拟合效果好于全局模型,其中5 km尺度局域模型的MSE和AIC在所有模型中最小。研究区内更新株数呈南高北低、东西差异不明显的分布趋势。[结论] 在构建局域模型时,应考虑空间尺度的影响,5 km尺度局域模型可以很好模拟研究区天然次生林更新的空间分布状况,有效降低甚至去除空间自相关性。研究区的林分更新株数主要集中在1 000~2 000株·hm-2之间,更新等级处于不良水平,森林天然更新能力整体较弱,可采用人工促进天然更新等措施进行森林经营。

关键词: 森林更新, 空间自相关性, 地理加权泊松模型, 空间尺度, 大兴安岭

Abstract: [Objective] Through the analysis of the spatial correlation and spatial distribution pattern of forest regeneration in different locations of natural secondary forests in Daxing'an Mountains, this study was implemented to explore the impacts of scale effect on the spatial autocorrelation, to understand the potential regularity in regeneration dynamics from a deeper level by analyzing the influencing factors of forest regeneration, and finally to provide theoretical basis and technical support for the operation and decision of natural secondary forest in this area.[Method] We took Cuigang forestry station of Xinlin Forestry Bureau of Daxing'an Mountains in Heilongjiang Province as the research area. Based on the data of 45 permanent sample plots established in the research area from July to August 2018, we selected 9 factors in 5 aspects including stand factors, topographic factors, forest stand spatial structure, soil thickness and species diversity as the independent variables, and established the global Poisson regression model and geographically weighted Poisson regression(GWPR)model under 4 scales(5 km, 10 km, 15 km and 20 km)on the basis of geographically weighted regression model to simulate the regeneration status of natural secondary forest in this area. Global Moran I and local Moran I were used to respectively describe the global spatial autocorrelation and spatial distribution of model residuals, to evaluate the fitting effects of global model and of local models under different scales, and to explain the differences among the each local model under different scales. Finally, the local model under 5 km was adopted to draw the spatial distribution plan of forest regeneration in the research area so as to evaluate and analyze the forest regeneration in the research area.[Result] The local model under 5 km made the best local spatial distribution of model residuals, formed the ideal distribution state of aggregated distribution of different model residuals, the parameter estimates of the model variables produced the largest range of variation, and had the best stability. With the gradual increase of the spatial scale, the stability of the model declined gradually, but still generally better than that of the global model. Meanwhile, the local model under this scale showed the lowest spatial autocorrelation of model residuals. The fitting effects of local models were better than those of global models, where, the local model under 5 km had the minimum MSE value and AIC value. In the research area, the number of regeneration individuals in the south part was larger than the one in the north part, while the differences between the east and the west were not obvious.[Conclusion] The influences of spatial scale shall be taken into consideration when the local model is established. The local model under 5 km adopted in this study could well simulate the spatial distribution of natural secondary forest regeneration in the research area, and could effectively reduce or even remove the spatial autocorrelation. In the research area, the number of stand regeneration individuals is mainly between 1 000-2 000 hm-2, that is, the overall regeneration is at a bad level, and the natural forest regeneration ability is generally weak, therefore, management measures such as artificial promoting of natural regeneration shall be taken for forest management.

Key words: forest regeneration, spatial autocorrelation, geographically weighted Poisson regression model, spatial scale, Daxing'an Mountains

中图分类号: