• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖
Volume 33 Issue 5
Oct.  2020
Article Contents
Turn off MathJax

Citation:

Study on Age-independent Tree Model: Taking the Average Height of Larix gmelinii as an Example

  • Corresponding author: LEI Yuan-cai, yclei@caf.ac.cn
  • Received Date: 2020-01-08
    Accepted Date: 2020-04-15
  • Objective To construct an age-independent theoretical growth model for simulating forest growth. Method The Richards model was used to fit the average height of Larix gmelinii stand in Jilin Province by using the difference form of implicit age. The growth model independent of age was constructed, and the free parameters were expressed as functions of stand correlation factor and site factors. Result The Jackknife method was used to validate the model. The age-independent model correction determination coefficient was between 0.92 and 0.94, while the age-dependent model correction determination coefficient was 0.85, and the age-independent model correction determination coefficient was increased by 8%~10%. The root mean square error (RMSE) of the age-independent model was between 1.30 and 1.59, which was much less than that of the age-dependent model (20.48). Comparing the age-independent models, introducing climate factor model to correct the RMSE and total relative error (TRE) of the determination coefficients was significantly lower than that without introducing, which showed that the introduction of climate factor improved the simulation effect of the model. Conclusion The age-independent theoretical models can be applied to stand where the age is unavailable or unrelated. In this paper, the age-independent site factors and climate factors are introduced into the model fitting in the way of free parameters. By comparing the two kinds of models with fixed parameters and free parameters, it is found that the model with free parameters performs better. This paper argues that the model is more effective by introducing site factors and climate factors into the average height growth model.
  • 加载中
  • [1] 王明亮, 唐守正. 标准树高曲线的研制[J]. 林业科学研究, 1997, 10(3):259-264. doi: 10.3321/j.issn:1001-1498.1997.03.006

    [2] 张连金, 惠刚盈, 孙长忠, 等. 木麻黄生长与收获模型系统的研究[J]. 林业科学研究, 2014, 27(1):59-65.

    [3] 李风日, 吴俊民, 鲁胜利. Richards函数与Schnute生长模型的比较[J]. 东北林业大学学报, 1993, 21(4):15-24.

    [4]

    Martin G L, Ek A R. A comparison of competition measures and growth models for predicting plantation red pine diameter and height growth[J]. For Sci, 1984, 30: 731-743.
    [5]

    Monserud R A, Sterba H. A basal area increment model for individual trees growing in even-and uneven-aged stands in Austria[J]. For Ecol Manage, 1996, 80: 57-80. doi: 10.1016/0378-1127(95)03638-5
    [6]

    Trasobares A, Tomé M, Miina J. Growth and yield model for Pinus halepensis Mill. in Catalonia, north-east Spain[J]. For Ecol Manage, 2004, 203: 49-62. doi: 10.1016/j.foreco.2004.07.060
    [7]

    Vanclay J K. Modelling forest growth and yield — applications to mixed tropical forests[M]. CAB International, Oxon, UK. 1994.
    [8]

    Vanclay J K. Aggregating tree species to develop diameter increment equations for tropical rainforests[J]. For Ecol Manage, 1991, 42: 143-168. doi: 10.1016/0378-1127(91)90022-N
    [9]

    Huang S, Titus S J. An age-independent individual tree height prediction model for boreal spruce-aspen stands in Alberta[J]. CAN J FOR RES, 1994, 24(7): 1295-1301. doi: 10.1139/x94-169
    [10]

    Bailey R L, Clutter J L. Base-age invariant polymorphic site curves[J]. For Sci, 1974, 20: 155-159.
    [11]

    Borders B E, Bailey R L, Ware K D. Slash pine siteindex from polymorphic model by joining (splining) nonpolynomial segments with an algebraic difference method [J]. For Sci,, 1984, 30: 411-423.
    [12]

    Amaro A, Reed D, Tomé M, et al. Modelling dominant height growth: Eucalyptus plantations in Portugal[J]. For Sci, 1998, 44: 37-46.
    [13]

    Cieszewski C J, Bailey R. Generalized algebraic difference approach: theory based derivation of dynamic equations with polymorphism and variable asymptotes[J]. For Sci, 1998, 46: 116-126.
    [14]

    Rivas J J C, González J G, González A, et al. Compatible height and site index models for five pine species in El Salto, Durango (Mexico)[J]. For Ecol Manage, 2004, 201: 145-160. doi: 10.1016/j.foreco.2004.05.060
    [15]

    Tomé J, Tomé M, Barreiro S, et al. Age-independent difference equations for modelling tree and stand growth[J]. Can J For Res, 2006, 36(7): 1621-1630. doi: 10.1139/x06-065
    [16]

    Gea-lzquierdo G, Canellas I, Montero G. Site index in agroforestry systems: age-dependent and age-independent dynamic diameter growth models for Quercus ilex in Iberian open oak woodlands[J]. Can J For Res, 2008, 38: 110-113.
    [17] 葛宏立, 项小强, 何时珍, 等. 年龄隐含的生长模型在森林资源连续清查中的应用[J]. 林业科学研究, 1997, 10(4):420-424. doi: 10.3321/j.issn:1001-1498.1997.04.014

    [18] 邹奕巧, 杜 群, 葛宏立. 有年龄生长模型应用于无年龄情况研究[J]. 浙江农林大学学报, 2012, 29(6):889-896.

    [19] 高东启, 邓华锋, 程志楚, 等. 北京市蒙古栎单木与年龄无关的生长预测模型研究[J]. 福建林学院学报, 2013, 33(4):305-309. doi: 10.3969/j.issn.1001-389X.2013.04.004

    [20] 邱思玉, 曹元帅, 孙玉军, 等. 杉木人工林与年龄无关的优势高生长模型[J]. 南京林业大学学报: 自然科学版, 2019, 43(5):121-127.

  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

Figures(2) / Tables(6)

Article views(4042) PDF downloads(18) Cited by()

Proportional views

Study on Age-independent Tree Model: Taking the Average Height of Larix gmelinii as an Example

    Corresponding author: LEI Yuan-cai, yclei@caf.ac.cn
  • Research Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China

Abstract:  Objective To construct an age-independent theoretical growth model for simulating forest growth. Method The Richards model was used to fit the average height of Larix gmelinii stand in Jilin Province by using the difference form of implicit age. The growth model independent of age was constructed, and the free parameters were expressed as functions of stand correlation factor and site factors. Result The Jackknife method was used to validate the model. The age-independent model correction determination coefficient was between 0.92 and 0.94, while the age-dependent model correction determination coefficient was 0.85, and the age-independent model correction determination coefficient was increased by 8%~10%. The root mean square error (RMSE) of the age-independent model was between 1.30 and 1.59, which was much less than that of the age-dependent model (20.48). Comparing the age-independent models, introducing climate factor model to correct the RMSE and total relative error (TRE) of the determination coefficients was significantly lower than that without introducing, which showed that the introduction of climate factor improved the simulation effect of the model. Conclusion The age-independent theoretical models can be applied to stand where the age is unavailable or unrelated. In this paper, the age-independent site factors and climate factors are introduced into the model fitting in the way of free parameters. By comparing the two kinds of models with fixed parameters and free parameters, it is found that the model with free parameters performs better. This paper argues that the model is more effective by introducing site factors and climate factors into the average height growth model.

  • 了解和掌握森林资源的动态变化,对制定森林抚育经营方案和措施,促进森林质量精准提升具有重要的意义。森林生长模型是了解和掌握森林资源动态变化的一项重要工具。这类方程的形式很多,可以分为理论生长函数、经验生长函数和Box-Lucas生长函数,而著名的理论生长函数有Richards模型,Schumacher模型,Schnute模型等[1-3]。这些模型一直被广泛应用,如模拟树高生长,优势高生长,断面积生长等。森林生长模型在构建和应用研究中,通常使用年龄作为重要因子,但林木年龄因子往往存在两方面的问题,一是测定困难,一般需要用生长锥测定具有清晰年轮树种的从髓心开始的整个年轮的数据,对于没有清晰年轮的树种,测定往往非常困难,而绝大多数阔叶树不具有清晰的年轮;二是年龄不一定有效,在树木年龄差异较大的林分中,年龄不能用于预测林木或林分生长和产量,特别是对于天然林的生长和动态监测更加依赖于年龄无关的生长模型。因此,年龄无关的森林生长模型的研究变得越来越必要。

    国内外许多学者开展了一些与年龄无关的森林生长模型的相关研究。Martin等[4],Monserud等[5],Trasobares等[6]通过使用与年龄无关的经验线性模型模拟了单木的生长,但是这类模型可能不能确保单木或林分生长有适合的生长曲线特性,即生长曲线达到最大生长量后生长逐渐下降。而应用理论生长函数这个特性进行外推预测林木生长时,可以保证单木或林分的生物学特性。

    Vanclay[7-8]考虑单木或林分生长曲线的生物学特性,构建了具有这种性质的线性经验函数,并已经被许多作者成功地使用。这类生长函数的表达式如下:

    其中Y为林木变量生长量,dYdt为林木变量生长率,k可以取为1或2,常数β0可以表达为环境变量和(或)竞争变量的函数。然而,这些参数没有特定生物学意义的含义。

    Huang和Titus[9]为了模拟阿尔伯达北部混交林中白云杉(Picea glauca (Moench) Voss)的直径增量生长变化,比较了前面生长函数的非线性形式与能描述单木或林分生长曲线形状的非线性Box-Lucas函数,结论是Box-Lucas函数模拟直径生长效果比非线性的经验生长函数好。Box-Lucas函数表达式为:

    其中:θ1θ2是估计增量曲线形状的参数,dYdt为林木变量生长率。这些函数和其他类似的函数可以用来对单木或林分变量的增量进行建模。这些参数也不具有特定生物学意义的含义。

    理论生长函数通常建立在公认的生物学关系的基础上,因此它们的模型参数可以从生物学的角度进行解释。不管树种和林分类型,树木的生长都遵循某种S型或凸型生长曲线,它显示出树木寿命的3个阶段。模拟树木和林分生长的一种非常成功的方法是使用用差分方程表示的生长函数,如Bailey等[10]和Borders等[11]提出的,已被许多学者使用[12-14]。差分方程表示一组生长函数,除了一个“自由”参数外,所有参数都是公共的。但多数这类差分方法并没有隐去年龄。Tomé等[15]利用隐去年龄的差分方法将Lundqvist模型和Richards模型分别应用于葡萄牙桉树优势高生长和软木栎胸径生长的研究;Gea-lzquierdo等[16]基于Hossfeld IV和korf模型的差分形式模拟树木直径的生长。葛宏立[17]以Johnso -Schumacher 模型为基础,导出年龄隐含的模型应用于森林资源连续清查数据的生长研究,邹奕巧等[18]用丽水市的马尾松样木胸径数据,导出了Richards模型隐去年龄的形式,并采用单点估计法和双点估计法2种方法进行研究。但是,此类研究目前为止相对较少,而且这些研究往往没有或者较少考虑到环境因子对于生长的影响。本研究以理查德模型为基础研究隐含年龄的差分形式,并建立差分方程的系数参数与环境因子和林分因子的相关函数,以便更好的解决年龄不可及或者不是有效变量的问题和更全面解释单木和林分生长的生物学特性,为混交林质量精准提升提供一种有效的经营工具。

1.   材料和方法
  • 研究采用2004和2009年吉林省两期177块落叶松固定样地数据,测量间隔为5 a。固定样地为20 m×30 m,主要调查因子有胸径、林分年龄、平均树高、郁闭度、地位级、坡位、坡向、坡度、海拔高度、土层厚度和土壤类型等(见表1)。落叶松样地的选取原则是树种组成中,落叶松的胸径断面积超过70%,平均高通过每块样地3株平均木的树高取平均得到,平均年龄通过生长锥钻取3株平均木的树芯计算的年龄平均值得到。研究区气候数据从WorldClim数据库下载得到(网址:www.worldclim. org)。每个样地可以下载19个用于构建气候敏感的生物量模型的候选气候因子(表2),分辨率为1 km×1 km。这些因子可以分为两大类,即温度有关的变量10个和降水量有关的变量9个。每个气候因子是1990—2000年的平均值。

    变量名
    Variables
    平均胸径
    Mean DBH/cm
    平均年龄
    Mean Age/a
    平均高
    Mean Height/m
    BIO12/
    mm
    BIO13/
    mm
    BIO8/℃土壤厚度
    Soil depth/cm
    断面积
    Basal area/(m2·hm−2
    平均数Mean 12.82 34.05 12.68 5.87 98.70 19.00 43.23 11.24
    标准差
    Stand deviation
    4.51 25.95 5.08 1.85 4.22 1.71 10.03 7.07
    最小值 Min 5.47 6.00 3.50 2.02 84.44 15.01 20.00 0.28
    最大值 Max 28.75 70.00 33.00 10.00 109.87 22.05 80.00 37.50
      注:平均温度BIO8、年降水雨量BIO12和最湿季度降水量BIO13
      Note:Mean temperature of wettest quarter BIO8, annual precipitation BIO12 and precipitation of wettest month BIO13

    Table 1.  Stand and climate factors

    利用双重筛选逐步回归方法,建立两期平均高和19个气候因子的线性关系,入选的因子为最湿季度的平均温度BIO8、年降水量BIO12和最湿季度的降水量BIO13,这3个气候因子在建模中最为敏感,本研究予以选用。具体含义见表2

    气候参数
    Climatic factors
    描述
    Descriptions
    BIO8/℃最湿季度的平均温度
    Mean temperature of wettest quarter
    BIO12/mm年降水量
    Annual precipitation
    BIO13/mm最湿月的降水量
    Precipitation of wettest month

    Table 2.  Implications for selected climate factors

  • 理查德方程认为一个器官的生长可以表达成有机物合成和降解之差。假设合成可以被表达成质量的异速生长关系,而降解和质量是成比率的,因此生长可以被表达为:

    式中:β1β2m是参数。生长方程是不基于年龄的,因此对于异龄林样地或者年龄未知的异龄林样地数据结构模拟单木或林分生长是非常适用的。实际应用中,上式瞬时生长量用单木或林分收获量表达即形式可以由式(3)积分为:

    式中:参数可以分别由下式表达:

    A是当t趋于无穷时的渐近线,k是和生长率相关的参数,m是形状参数。$ c={\mathrm{e}}^{-k{t}_{0}} $t0Y=0时的年龄。大多数应用中,如果t0=0时C=1。这里取C=1.

    与年龄无关的差分形式可以通过类似于通常用于将生长函数公式化为差分方程的方法(如Amaro等[12]的过程来获得)。考虑到吉林省177块天然落叶松林的连续清查数据缺乏年龄,为了建立单木或林分的树高生长模型,本研究基于理查德方程(4)推导不含年龄的差分生长方程。根据式(4),第i年和第i+a年的平均高为:

    从而得到

    式(9)和式(10)两个式子相减得到

    因此,第i+a年的差分生长方程就可以写成:

    在这里,a为预测期,或为连续清查时间间隔,Yi+a是预测期(或连续清查时间间隔)的单木或平均高,Yi是连续清查时间初期的单木或平均高,A是当t趋于无穷时的渐近线,k是和生长率相关的参数,m是形状参数。式(12)表达比较复杂,但是当拟合单木树高或平均树高生长时参数很容易收敛。

    为了模拟平均树高生长,本研究分别尝试采用k为固定参数和自由参数两种形式。也考虑了m为自由参数的情况,但获得的结果是m受林分因子、气候因子的影响不明显。当k为自由参数时,将k表达为林分因子和气候因子的线性函数关系。

    其表达形式分别有如下4个形式,并进行比较。

    在这里,$ {b}_{0} $~$ {b}_{5} $是估计参数,BA是样地断面积,S是土壤厚度,BIO8为最湿季度的平均温度,BIO12是年降水量,BIO13是最湿月份的降水量。

  • 选择平均误差($\bar e$),均方根误差(RMSE),总相对误差(TRE),确定系数(${R^2}$)4种统计量来检验模型的拟合优度和误差,并观察残差分布是否有异质性。利用修正确定系数AjustedR2对模型进行检验和评价(式22)。

    其中:${y_t}$${{\overset{\frown} y}_t}$,nk分别表示第t个观测值、第t个估计值、观测值数和回归模型参数个数。

  • 对模型间的差异采用F统计检验,F统计指标按照下式计算:

    其中:ssesdfs分别为Richards模型1的回归模型残差平方和和自由度,ssecdfc分别为Richard模型2的回归模型残差平方和和自由度。

    利用上式计算的值与自由度分别为dfs-dfcdfsF临界值进行比较,以确定两个模型之间是否存在显著差异。

2.   结果
  • 表3是利用最小二乘法拟合的模型式(13)~(15)、式(12)和式(4)的参数值。从表中可以看到,所有模型的渐近线A值在17.53~22.72之间,年龄无关的平均高生长模型式(13)~(15)和式(12)的渐近线在21.85~22.72之间,更符合落叶松最大高的实际情况。而年龄相关的模型式(4)的渐近线只有17.53,明显偏低。

    模型 ModelsAkcm
    式13 Formula 1323.44530.030~0.057 40.110 1
    式14 Formula 1421.853 20.017~0.065 80.122 3
    式15 Formula 1521.933 70.0218~0.058 00.135 6
    式16 Formula 1622.000 00.0210~0.05 90.114 6
    式12 Formula 1222.705 20.03380.157 1
    式4 Formula 417.705 20.05560.0150.985 2

    Table 3.  Model parameters

    表4中,利用刀切法对模型进行验证,每次从建模样本中随机抽取50%的样本,使用同一模型拟合,重复10次,进行模型参数的灵敏度分析来评价模型;拟合结果使用验证样本年龄无关的模型式(13)~(16)和式(12)的修正确定系数在0.92~0.94之间,而年龄相关的模型修正确定系数为0.85,时间无关的模型修正确定系数提高了8%~10%;模型式(13)~(16)和式(12)的均方根误差RMSE在1.30~1.59之间,远远小于年龄相关的模型式(4)的20.48。如图1a图1b所示,时间和平均高的关系很难用一条曲线来描述,而用年龄无关的模拟拟合的平均高和实测平均高相关系数却高达0.97。

    Figure 1.  a The relationship between mean height and Age b The relationship between fitting data and measured data of formula 16

    模型
    Models
    式13
    Formula 13
    式14
    Formula 14
    式15
    Formula 13
    式16
    Formula 16
    b0 0.0192 0.109 285 0.163 349 0.159 205
    b1 −5.3E-05 −2.3E-05
    b2 0.009 78 0.009 982 −4.3E-05
    b3 0.000 615 0.000 451 0.001 235
    b4 0.001 268 0.001 228 −0.001 39
    b5 −0.001 06 −0.001 41 0.114 556

    Table 4.  Parameters of linear function

    年龄无关的几个模型进行对比,从表5F值来看,F0.05在2.42~3.89之间,F值在0.15~81.05之间。从表中可以看出,模型式(15)和式(16)的F值是0.15,F0.05是3.89,说明这两个模型差异不显著;同样模型式(15)和式(13)的F值是2.98,F0.05是3.89,说明这两个模型差异也不显著。其它模型两两对比,都差异显著。

    模型
    Models
    式13
    Formula 13F/F0.05
    式14
    Formula 14 F/F0.05
    式15
    Formula 15F/F0.05
    式14
    Formula 14F/F0.05
    式12
    Formula 12 F/F0.05
    式13
    Formula 13
    式14
    Formula 14
    5.97/2.65
    式15
    Formula 15
    2.98/3.8912.22/3.89
    式16
    Formula 16
    5.82/3.046.21/3.040.15/3.89
    式12
    Formula 12
    15.37/3.897.94/3.0420.84/3.0410.36/3.04
    式4
    Formula 4
    161.78/3.8940.73/2.4254.03/2.6581.05/3.04
      注:−表示两个模型自由度相同,无法计算F
      Note:− means that the two models have the same degree of freedom, so the value of F cannot be calculated

    Table 5.  F parameter of models

    模型式(13)和式(14)对比,式(14)引入了气候因子,气候因子的选用参考了José Tomé等[12]的文章,从F值来看,差异显著,同时式(14)的修正决定系数0.941 9,也高于式(13)的0.935 2(表6),而式(14)的均方根误差RMSETRE较式(13)也有明显的降低,说明气候因子的引入提高了模型的模拟效果。图2是6个模型的残差分布图,从图上看出年龄无关的模型(式12~16)的残差分布均优于年龄相关的模型(式4)。

    Figure 2.  Distribution of residuals errors (a. Formula 12, b. Formula 13, c. Formula 14, d. Formula 15, e. Formula 16, f. Formula 4)

    模型 Models$\bar e$RMSETREAjuestR2
    式13 Formula 130.000 30.513 31.450 40.935 2
    式14 Formula 140.002 50.486 01.300 40.941 9
    式15 Formula 150.001 00.504 31.400 40.937 4
    式16 Formula 160.000 30.504 61.401 60.937 4
    式12 Formula 12−5.881 356E-060.537 61.590 90.928 9
    式4 Formula 40.000 31.929 020.483 70.851 5
      注: $\bar e$,平均误差; RMSE, 均方根误差; TRE, 总相对误差; AjustedR2,修正确定系数
      Note: $\bar e$, mean bias; RMSE, root mean square error; TRE, total relative error; AjustedR2, Correction coefficient of determination

    Table 6.  Model statistics

    本研究在模型式(13)~(16)中采用自由参数k的形式,将k表达成立地因子和气候因子的线性形式,而模型式(12)采用固定参数的形式,对比两类模型发现,模型式(13)~(16)的修正可绝系数均大于式(12),而误差均小于式(12)。

3.   讨论
  • 本研究基于Richards函数建立了年龄无关的差分方程,从而可用于年龄未知或者年龄不是有效因子的情况下的林分。选择年龄已知的固定样地,对比了年龄无关模型和年龄有关模型的模拟精度。年龄无关的生长模型比年龄相关的模型精度高而误差小。葛宏立等(17)采用年龄无关模型方法研究马尾松生长模型时预测精度高达99.80%,高东启等[19]用年龄无关的模型拟合蒙古栎生长,模型的R 都大于0.98,预估精度在99.40%~99.72%之间,邱思玉等[20]分别用不同的年龄无关树高模型建立优势高模型也得到了较高的预估精度,研究结果与本研究基本一致。分析原因,可能是因为年龄无关的生长模型直接以单木自身的观测值为自变量来预测后期的生长情况,这在一定程度上淡化了林相、立地、竞争等因子的影响。

    另外,从拟合得到的渐近线来看,年龄无关模型拟合的渐近线更符合落叶松高的实际情况,分析原因,作者认为一方面年龄无关的模型具有较好的预测能力。另一方面外业调查获得的林分年龄存在一定的偏差,这也和年龄因子难于获得有关。

    同时,本研究采用自由参数的方式将与年龄无关的立地因子和气候因子引入模型的拟合中,对比固定参数和自由参数的两类模型,可以发现自由参数的模型表现更好。因此认为参数k与立地因子和气候因子相关性强,这和Tomé等[15]的结论一致。引入气候因子的模型比仅引入林分因子的模型表现更好,说明考虑气候对生长的影响更有利于优势高生长的模拟。

    但是,以5 a为预测间隔期是比较短的,树木生长短期内不会有太大的生长变化,使得模型的预测误差较小,从而保证了模型有较高的预测精度。模型仍然需要三期及以上的数据进行进一步的检验,从而分析在年龄间隔较大时能否保证较高的预测精度。

Reference (20)

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return