• 论文与研究报告 •

基于空间自相关的天然蒙古栎阔叶混交林林木胸径-树高模型

1. 1. 中国林业科学研究院资源信息研究所 北京 100091;
2. 宁波市农业科学研究院 宁波 315040;
3. 江西农业大学林学院 南昌 330045
• 收稿日期:2015-12-07 修回日期:2016-04-13 出版日期:2017-06-25 发布日期:2017-07-14
• 通讯作者: 张会儒
• 基金资助:

Individual Diameter-Height Models for Mixed Quercus mongolica Broadleaved Natural Stands Based on Spatial Autocorrelation

Lou Minghua1,2, Zhang Huiru1, Lei Xiangdong1, Li Chunming1, Zang Hao3

1. 1. Research Institute of Forest Resource Information Techniques, CAF Beijing 100091;
2. Ningbo Academy of Agriculture Sciences, Zhejiang Province Ningbo 315040;
3. College of Forestry, Jiangxi Agricultural University Nanchang 330045
• Received:2015-12-07 Revised:2016-04-13 Online:2017-06-25 Published:2017-07-14

Abstract: [Objective] Considering spatial autocorrelation among individuals, individual diameter-height models based on spatial autocorrelation were constructed. It may provide a theoretical basis for sustainable management of natural mixed forests. [Method] Three simultaneous autoregressive (SAR) models, including spatial lag model (SLM), spatial error model (SEM) and spatial Durbin model (or called spatial mixed model) (SDM) within seven spatial weight matrices, including Delaunay triangulation (DT), inverse distance raised to one power (ID1), inverse distance raised to two powers (ID2), inverse distance raised to five powers (ID5), spherical variogram (SV), gaussian variogram (GV) and exponential variogram (EV), was used to construct individual diameter at breast height and height models of mixed Quercus mongolica broadleaved natural stands in Northeast China, and treating linearization base model (BM) as a benchmark model. Model parameters of BM were estimated by ordinary least squares (OLS), model parameters of three SAR models were estimated by maximum likelihood. Model coefficients β0 and β1 of four models were tested by T-test, the autoregressive parameters ρ, γ and λ were all tested by likelihood ratio test. Moran's I (MI) was selected to compared autocorrelation of four model residuals. Three statistics, i.e. coefficient of determination (R2), root mean square error (RMSE) and Akaike information criterion (AIC), were regarded as the appropriate criteria to identify the model fitting among BM, SLM, SDM and SEM. [Result] MI values of BM residuals were larger than 1, when applying SV into BM. Therefore, SV was the unreasonable spatial weight matrix and did not regard as a spatial weight matrix in the following result analysis. MI values of BM and SLM residuals were significantly larger than the expected value I0 of MI in the all spatial weight matrices (except SV). MI values of SLM residuals were smaller than those of BM using the same spatial weight matrix. The difference between MI values of SDM residuals and I0 was not significant in other four spatial weight matrices, except GV and ID1. Similarly, the difference between MI values of SEM residuals and I0 was not significant in other five spatial weight matrices, except ID1. Three criteria of three SAR models were all better than those of BM. Using the same spatial weight matrix, MI values of SDM were very similar to those of SEM, meanwhile, MI values of SDM and SEM were both larger than those of SLM. Different spatial weight matrices (except GV) in SDM and SEM were sorted from best to worst according three criteria and the ranking was: ID2 > DT > ID > ID5 > EV. Model coefficients β1 of three SAR were very similar to those of BM, regardless of which spatial weight matrix was used. Compared with β1, model coefficients β0 of SEM were similar to those of BM, while model coefficients β0 of SDM and SLM were different to those of BM, and were changed along with the different spatial weight matrix. Among all spatial weight matrices within three SAR models, the autoregressive parameters ρ, γ and λ using ID1 were larger higher than any other spatial weight matrix. GV only applied to SEM, rather than SDM, could make the autoregressive parameter λ significant not equal to zero. The autoregressive parameters ρ, γ and λ were all not equal to zero using five spatial weight matrices (except GV).[Conclusion] Among all spatial weight matrices applied in three SAR models, SV and ID1 are the unreasonable spatial weight matrices. SLM do not remove, but reduce the spatial autocorrelation of model residuals, and slightly improve the model fitting. Model fitting of SLM was worse than those of SDM and SEM. Selecting appropriate spatial weight matrices, SDM and SEM can remove the spatial dependence of model residuals and improve the model fitting. ID2 is the best one among these selected appropriate spatial weight matrices. The diameter-height models of Quercus mongolica, Populus-Betula (Populus davidiana and Betula platyphylla) and Pinus koraiensis were constructed by species dummy variables SAR models based on ID2 and SEM.