Turn off MathJax
Article Contents

ZHANG Dongyan, WANG Dongzhi, LI Xiao, GAO Yushan, LI Tianyu, CHEN Jing. Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression[J]. Journal of Zhejiang A&F University. doi: 10.11833/j.issn.2095-0756.20190461
Citation: ZHANG Dongyan, WANG Dongzhi, LI Xiao, GAO Yushan, LI Tianyu, CHEN Jing. Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression[J]. Journal of Zhejiang A&F University. doi: 10.11833/j.issn.2095-0756.20190461

Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression

doi: 10.11833/j.issn.2095-0756.20190461
  • Received Date: 2019-08-05
  • Rev Recd Date: 2019-11-22
  • Available Online: 2019-12-17
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(3)  / Tables(3)

Article views(487) PDF downloads(12) Cited by()

Related
Proportional views

Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression

doi: 10.11833/j.issn.2095-0756.20190461

Abstract:   Objective  With the employment of nonlinear quantile regression method using dummy variables, the current study is aimed to establish a prediction model for the relationship between height and diameter at breast height (DBH) in Larix principis-rupprechtii and Betula platyphylla mixed forest so as to better predict the tree structure and site productivity of mixed forests.   Method   Taking L. principis-rupprechtii and B. platyphylla mixed forest of Saihanba Mechanised Tree farm in Hebei Province as the research object, with 83 pieces of standard land survey data used and dummy variables created, this paper adopted the least square method and nonlinear quantile regression method respectively in the construction of the relationship model of tree height and DBH of different species.   Result   The accuracy of the nonlinear quantile regression prediction model based on dummy variables was higher than that of the one constructed using the least square method. Specifically, when the least square method was used to fit the tree height and DBH relationship model of different tree species, the determination coefficient, average difference and average absolute error of different tree species models were within the range of 0.787−0.814, 1.581−1.877 and 2.447−2.654 respectively. When the nonlinear quantile regression method was used, the coefficient, average deviation, and average absolute error were within the range of 0.839−0.921, 0.213−1.469, 0.561−2.322. In accordance with the residual analysis, when the quantiles of τ is 0.7, the relationship model of tree height and DBH of different species demonstrated a higher accuracy.   Conclusion   To sum up, compared with the one constructed employing the least square method, the prediction model of tree height and DBH relationship of different tree species adopting the nonlinear quantile regression method has higher prediction accuracy. [Ch, 3 fig. 3 tab. 33 ref.]

ZHANG Dongyan, WANG Dongzhi, LI Xiao, GAO Yushan, LI Tianyu, CHEN Jing. Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression[J]. Journal of Zhejiang A&F University. doi: 10.11833/j.issn.2095-0756.20190461
Citation: ZHANG Dongyan, WANG Dongzhi, LI Xiao, GAO Yushan, LI Tianyu, CHEN Jing. Relationship between height and diameter at breast height(DBH) in mixed coniferous and broadleaved forest based on quantile regression[J]. Journal of Zhejiang A&F University. doi: 10.11833/j.issn.2095-0756.20190461
  • 树高和胸径不仅是用来预测林分蓄积量[1]、生物量[2]、立地生产力[2-3]及林分结构[4]的重要变量,而且是森林资源调查及经营效果评价的重要因子。在标准地调查过程中,树高测量难度较大且观测成本高,其观测误差也相对较大,给精准林业质量提升带来了一定困难[2,4],而胸径观测方便且精度较高。因此,根据标准地调查数据,建立树高与胸径关系预测模型,可降低调查成本,提高预测精度[5],这对于森林质量精准提升具有重要意义。当前林业研究多采用线性或非线性树高与胸径关系模型来模拟预测两者之间的关系,其参数估计方法多采用最小二乘法来模拟,然而最小二乘法是基于均值回归,利用变量均值来拟合模型参数[6],该方法要求调查数据需满足独立正态同分布等条件。在林业调查中,树高与胸径观测数据不能满足该要求,而分位数回归对调查数据没有严格要求[7],利用变量条件分位数来建模[6],对具有尖峰、厚尾、异方差显著的数据拟合效果更加稳健[8-10]。分位数回归理论框架是由KOENKER[11]提出,已在医学[12]、经济学[13]、教育与政策[14-15]及自然资源管理等领域进行了研究与应用。在林业相关研究中,分位数回归被应用于模拟林分自疏边界线[16]、直径分布规律[17]、林分密度指数[18]及森林病虫害[19]等方面研究。ÖZÇELIK等[20]基于分位数回归建立了树高与胸径关系模型,高慧淋等[21]采用此方法建立了长白落叶松Larix olgensis人工林最大林分密度线模型,提高了模型预测精度及适用性。然而在华北暖温带针阔混交林中,如何基于一个分位数回归模型,预测不同树种树高与胸径关系是亟待解决的科学问题。在混交林中为了描述树种结构对树木生长影响,部分学者[22-25]采用哑变量方法构建了不同间伐方式、不同地域树高曲线及生长量预测模型。然而基于包含哑变量的非线性分位数回归方法来构建不同树种树高与胸径关系模型的研究较少。因此,本研究以河北省塞罕坝华北落叶松Larix principis-rupprechtii与白桦Betula platyphylla针阔混交林为研究对象,基于哑变量方法和分位数回归相结合方法,构建混交林不同树种分位数回归模型,为精确描述树高与胸径的关系提供依据。

  • 河北省塞罕坝机械林场(41°22′~42°58′N,116°53′~118°31′E)位于河北省最北部,地势北高南低,属华北暖温带立地类型区,林场总面积约9.2×104 hm2,总蓄积约8.1×106 m3。土壤类型以褐色森林土、棕色森林土及风沙土等为主;成土母质主要为坡积物、残积物及洪积物等;极端最高气温为33.4 ℃,最低气温−43.3 ℃,年均气温−1.3 ℃,年均无霜期64 d,年均降水量460.3 mm,是典型的半干旱半湿润寒温性大陆季风气候。研究区植被类型丰富,主要乔木树种有华北落叶松、白桦、樟子松Pinus sylvestris、云杉Picea asperata等,主要灌木树种有山刺玫Rosa daverica、胡枝子Lespedeza bicolor、沙棘Hippophae rhamnoides等,主要草本植物有地榆Sanguisorba offcinalis、唐松草Thalictrum aquilegifolium、曼陀罗Datura stramonium等。

  • 在北曼店、大唤起、阴河、千层板和第三乡等5个林场设立了83块标准地(30 m×30 m),对标准地内各林分因子(林分密度、平均高、平均胸径、树种断面积、林分总断面积、优势高等)和立地因子(海拔、坡度、坡向、坡位、土层厚度等)进行调查,共调查立木10 104株(华北落叶松5 258株,白桦4 846株),林分年龄分布为24~45 a,不同标准地混交度分布为0.39~0.62。研究过程中,分树种将观测数据分别按3∶1分为建模数据(62块标准地)和检验数据(21块标准地),基本信息如表1表2所示。

    统计量海拔/m林分断面积/hm2密度/(株·hm–2)华北落叶松白桦
    胸径/cm树高/m胸径/cm树高/m
    平均值1 424.00208.871 070.0015.5010.6014.3010.10
    最大值1 672.00452.191 725.0034.6017.9029.2016.20
    最小值1 177.00 19.96 675.00 5.00 4.20 5.10 4.60
    标准差 149.63 78.69 241.00 4.25 2.34 3.92 1.77

    Table 1.  Data of establishment model

    统计量海拔/m林分断面积/hm2密度/(株·hm–2)华北落叶松白桦
    胸径/cm树高/m胸径/cm树高/m
    平均值1 464.00254.001 162.0017.5011.0012.50 8.90
    最大值1 540.00387.402 000.0034.7017.8027.2015.60
    最小值1 396.00157.83 850.00 5.00 3.20 5.10 4.80
    标准差 39.28 73.41 269.76 6.26 2.70 4.47 2.16

    Table 2.  Data of test model

  • 在描述树木生长及树高与胸径关系的近百种不同模型中,Richard方程不但具有可解释的生物学意义,而且具有易收敛且灵活性高等特性。部分研究基于Richard方程构建了不同林分类型树高与胸径关系的预测模型,均取得了较好的预测结果[18, 20, 26-29]。因此,本研究以Richard方程作为构建华北落叶松与白桦针阔混交林树高与胸径关系基础模型,模型表达如式(1)所示。

    式(1)中:Hij为第i个样地第j株树的树高(m);dij为第i个样地第j株树的胸径(cm);abc为基础模型的参数;εij为误差项。

  • 为了解决模型预测精度的影响,可以在模型中加入哑变量[23, 30-32]。包含哑变量的树高与胸径关系预测模型,不仅可以实现模型对不同树种相容性,而且在一定程度上可以提供模型预测精度及适用性,包含哑变量的树高与胸径关系预测模型表达如式(2)所示。

    式(2)中:Mi为哑变量,当M1=1、M2=0时为华北落叶松,当M1=0、M2=1时为白桦;aibici为模型参数;εij为误差项。

  • 由于分位数回归对模型误差不需要严格假设条件,因此本研究基于Richard方程,选取5个分位点(τ=0.1、0.3、0.5、0.7、0.9)构建不同树种的树高与胸径关系预测模型,利用加权最小一乘法可以得到不同分位点参数估计值,具体见式(3)。

    式(3)中:S为不同分位点估计值;${\hat H_\tau }$Hij分别为第i个样地第j株树在不同分位点τ树高预测值与树高值(m);diji个样地第j株树胸径(cm);τ为分位点。

  • 统计分析均基于SPSS 24.0和SAS 9.4中的PROC NLIN和PROC NLP完成,基于模型确定系数(R2)、平均差(MD)、平均绝对误差(MAD)对模型拟合精度及适用性进行评价与比较。

    式(4)~(6)中:Hij${\hat H_{i\!j}}$${\overline H_{i\!j}}$分别为树高观测值、预测值和平均值;m为标准地数量;n为标准地株数。

  • 图1为不同树种胸径与树高的关系。华北落叶松和白桦的树高分别为4~18 m、6~16 m,胸径分别为6~32 cm、6~28 cm。

    Figure 1.  Height-diameter distribution of different tree species

  • 表3可见:在华北落叶松与白桦针阔混交林中,基于分位数回归的不同树种不同分位点确定系数均大于传统回归方法,平均差及平均绝对误差均小于传统回归方法。在确定的5个分位点中,当分位点τ=0.7时,华北落叶松与白桦的树高与胸径关系预测模型精度最高。基于不同分位点预测模型建立了不同树种在各分位点残差分布图(图2),确定当分位点位为τ=0.7时,华北落叶松与白桦树高与胸径关系模型能够更好描述两者之间的关系。

    .
    模型参数评价指标
    abcMDMADR2
    OLSL20.2990.0010.1551.5182.4470.814
    OLSB16.3010.0320.9061.8772.6540.787
    τL=0.117.4450.0330.9780.9622.0320.848
    τL=0.317.5760.0380.9820.7251.7090.855
    τL=0.517.6570.0430.9920.4651.3420.872
    τL=0.717.9980.0511.0020.2130.5610.912
    τL=0.918.0050.0651.0051.4322.2060.839
    τB=0.113.9560.0430.8920.9922.0430.848
    τB=0.314.5590.0520.9810.6811.8820.861
    τB=0.514.7030.0590.9970.4281.2690.882
    τB=0.714.9960.0711.0010.2350.8520.921
    τB=0.915.0110.0921.0091.4692.3220.838
    说明:L表示华北落叶松,B表示白桦;OLS表示最小二乘法

    Table 3.  Fitting and evaluation of basic model and quantile regression model

    Figure 2.  Residual distribution of each quantile of different tree species

  • 基于不同分位点预测模型,分别对华北落叶松和白桦的树高与胸径关系进行了模拟(图3)。不同树种在不同分位点树高与胸径关系预测趋势及范围基本一致,表明包含哑变量的分位数回归模型预测效果较好。

    Figure 3.  Relationship between tree height and DBH of different quantiles

  • 在描述树高与胸径关系的多种线性和非线性预测模型中,通常采用具有生物学意义且灵活性较高的Richard方程来研究不同林分类型树高与胸径的关系[26, 33],因此,本研究将Richard方程作为构建华北落叶松与白桦混交林树高与胸径关系基础模型。对红松Pinus koraiensis人工林、土耳其松Pinus brutia和黎巴嫩雪松Cedrus libani混交林的研究[20, 26-28]表明:Richard方程是描述其树高与胸径关系的最优模型。

    在构建树木生长及生物量预测模型中,包含哑变量模型具有更高的预测精度[22, 25, 31-32]。本研究在华北落叶松与白桦针阔混交林中,基于非线性分位数回归构建的包含哑变量树高与胸径的关系模型,其精度高于传统回归预测模型,这与人工林最大密度线确定[21]及树高与胸径关系模型[20]的研究结论相一致,表明非线性分位数回归较传统回归方法更加稳定,可用于人工林和混交林立地潜在生产力的评价。

  • 本研究以塞罕坝华北落叶松与白桦针阔混交林调查数据为基础,确定Richard模型为描述不同树种树高与胸径关系的基础模型,在基础模型中构造一个表示树种的哑变量,并利用分位数回归在一个模型中同时估计不同树种及不同分位点的树高与胸径关系模型参数,经过检验不同树种分位数回归模型均能较好反映树高与胸径的关系,当分位点τ=0.7时,分位数回归模型预测精度最高,拟合效果最好。

Reference (33)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return