• 中国中文核心期刊
  • 中国科学引文数据库(CSCD)核心库来源期刊
  • 中国科技论文统计源期刊(CJCR)
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

广东省针叶树种蓄积量和生物量生长模型研究

黄金金 刘晓彤 张逸如 李海奎

引用本文:
Citation:

广东省针叶树种蓄积量和生物量生长模型研究

    作者简介: 黄金金,硕士研究生。主要研究方向:林业统计与生物数学模型。E-mail:Huang202098@163.com.
    通讯作者: 李海奎, lihk@ifrit.ac.cn
  • 中图分类号: S718.55

Growth Models of Stand Volume and Biomass of Coniferous Species in Guangdong Province

    Corresponding author: LI Hai-kui, lihk@ifrit.ac.cn
  • CLC number: S718.55

  • 摘要: 目的 定量化表征森林蓄积量和生物量与年龄的关系,为预测森林蓄积量和生物量提供依据。 方法 基于广东省3个针叶树种5期连清数据,通过保留木林分平均胸径和林龄之间的关系,分别以生长潜力和生长速度分级、是否考虑竞争、分步和联合建模共8种组合,构建林分蓄积量和生物量生长模型,并评价模型拟合优度;以独立的4期连清数据为验证样本,评价模型的适用性。 结果 (1)不同种模型系确定系数皆在0.980以上,平均预估误差在±1%以内,总体相对误差在±0.5%以内。综合模型拟合和预测效果,以生长速度分级、不考虑竞争的分步模型最优,竞争对模型的适用性影响不大。(2)最优模型蓄积量和生物量的估计误差最大分别为10.36%和10.22%,模型适用性较好,4期估计误差表现为中期高于首末两期。(3)马尾松生长潜力最大,杉木最小,对杉木的估计效果优于马尾松和湿地松。(4)立地质量等级越高,生长量的极值越大,达到最高峰所需年限也更短;同一立地质量下马尾松生长量最大,其次为湿地松,杉木最小。 结论 含参数分级和林分特征的蓄积量和生物量模型,可以反映立地质量对广东省针叶树种蓄积量和生物量预测的影响,为精确估测森林储量提供方法学支持,也为其他地区林分水平生长模型的构建提供借鉴。
  • 图 1  建模数据固定样地分布示意图

    Figure 1.  Distribution map of fixed sample plots of modeling data

    图 2  8种蓄积量生长模型的区域尺度估计误差

    Figure 2.  Regional-scale estimation errors of eight volume growth models

    图 3  8种生物量生长模型的区域尺度估计误差

    Figure 3.  Regional-scale estimation errors of eight biomass growth models

    图 4  等级2~4的连年生长量和平均生长量

    Figure 4.  Annual increment and average increment from class 2 to class 4

    表 1  建模数据统计量

    Table 1.  Modeling data statistics

    树种 Tree species  统计量 Statistics  1997年2002年2007年2012年2017年
    马尾松 Pinus massoniana 胸径 DBH/cm 10.2±3.2 13.4±3.9 15.6±4.2 17.8±4.5 20.4±5.0
    蓄积量 Volume/dm3 43.7±51.3 83.0±75.0 120.5±92.8 164.3±112.6 228.3±143.1
    生物量 Biomass/kg 42.8±46.2 79.7±69.1 113.9±85.8 154.0±105.1 213.3±134.8
    湿地松
    Pinus elliottii
    胸径 DBH/cm 8.2±2.0 12.1±2.5 14.7±3.1 17.0±3.5 19.2±4.1
    蓄积量 Volume/dm3 19.2±13.5 54.5±30.1 91.1±47.6 132.5±66.7 183.7±99.6
    生物量 Biomass/kg 21.5±12.8 53.8±27.2 85.0±40.9 120.5±56.9 164.3±85.2
    杉木
    Cunninghamia lanceolata
    胸径 DBH/cm 8.8±2.8 11.7±3.3 13.5±3.6 15.0±3.7 16.7±4.1
    蓄积量 Volume/dm3 33.1±42.4 66.4±66.5 93.8±82.6 121.2±91.2 159.0±112.7
    生物量 Biomass/kg 22.7±25.1 43.0±38.7 59.2±47.9 75.2±52.7 97.2±65.0
      注:数据为平均值±标准差。
      Note: The data is average ± standard deviation.
    下载: 导出CSV

    表 2  验证数据统计量

    Table 2.  Verification data statistics

    树种 Tree species  统计量 Statistics2002年2007年2012年2017年
    马尾松 Pinus massoniana 胸径 DBH/cm 9.8±3.5 12.4±4.0 14.9±4.6 18.0±5.4
    蓄积量 Volume/dm3 37.1±43.0 65.6±61.1 102.2±84.9 162.6±120.7
    生物量 Biomass/kg 37.3±39.4 64.0±55.8 97.2±77.9 152.5±111.0
    湿地松 Pinus elliottii 胸径 DBH/cm 8.4±1.7 12.1±2.8 15.4±3.6 18.8±4.4
    蓄积量 Volume/dm3 19.2±11.9 54.1±33.5 101.8±56.5 169.1±94.5
    生物量 Biomass/kg 21.2±11.3 54.1±30.2 94.2±48.6 151.8±80.9
    杉木 Cunninghamia lanceolata 胸径 DBH/cm 8.8±2.6 11.2±2.9 13.2±3.2 15.2±3.3
    蓄积量 Volume/dm3 29.9±26.8 54.7±42.5 81.4±57.4 116.1±73.3
    生物量 Biomass/kg 20.7±16.4 36.0±25.1 51.9±33.5 72.3±42.4
      注:数据为平均值±标准差。
      Note: The data is average ± standard deviation.
    下载: 导出CSV

    表 3  8种模型系

    Table 3.  Eight model systems

    模型 Models方程 Equation模型 Models方程 Equation
    Model 1 $ {D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*\left({age}_{i} + {ad}_{j}\right)}) $
    $ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
    $ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
    Model 5 $ {D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)}) $
    $ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
    $ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
    Model 2 $ {D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)}) $
    $ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
    $ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
    Model 6 $ {D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)}) $
    $ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
    $ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
    Model 3 $ \left\{{array}{c}{D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $ Model 7 $ \left\{{array}{c}{D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $
    Model 4 $ \left\{{array}{c}{D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $ Model 8 $ \left\{{array}{c}{D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $
    下载: 导出CSV

    表 4  模型拟合评价指标

    Table 4.  Model fitting evaluation index

    模型
    Models
    胸径 DBH 蓄积量 Volume 生物量 Biomass
    R2SEE/cmTRE/%MPE/% R2SEE/dm3TRE/%MPE/% R2SEE/kgTRE/%MPE/%
    Model 1 0.984 0.77 0.00 0.36 0.985 13.42 −0.54 0.83 0.987 10.97 −0.08 0.79
    Model 2 0.985 0.75 0.00 0.35 0.982 14.54 −0.71 0.90 0.985 11.93 −0.27 0.85
    Model 3 0.981 0.74 −0.35 0.35 0.986 12.74 −0.90 0.79 0.989 10.32 −0.48 0.74
    Model 4 0.981 0.75 −0.20 0.35 0.986 12.67 −0.64 0.78 0.989 10.28 −0.14 0.74
    Model 5 0.985 0.74 0.00 0.35 0.986 12.83 −0.11 0.79 0.988 10.84 0.11 0.78
    Model 6 0.985 0.74 0.00 0.35 0.986 13.01 −0.12 0.80 0.987 11.02 0.10 0.79
    Model 7 0.981 0.74 −0.01 0.35 0.987 12.27 0.15 0.76 0.989 10.26 0.31 0.73
    Model 8 0.982 0.73 −0.04 0.34 0.988 12.15 0.00 0.75 0.989 10.17 0.17 0.73
    下载: 导出CSV

    表 5  树种哑变量参数和竞争指数参数

    Table 5.  Dummy variable parameters of tree species and competition index parameters

    模型
    Models
    因变量
    Dependent variables
    参数a
    Parameter a
    参数c
    Parameter c
    b0
    123123
    Model 5 胸径 DBH 58.15 41.57 36.59
    蓄积量 Volume 1 827.54 853.40 761.84 2.709 1 2.857 1 2.926 5
    生物量 Biomass 1 892.76 760.24 460.70 2.596 8 2.546 1 2.593 4
    Model 6 胸径 DBH 57.97 41.75 36.44 −0.281 0
    蓄积量 Volume 1 745.73 836.66 745.45 2.738 6 2.924 5 3.030 5
    生物量 Biomass 1 784.94 739.66 447.88 2.631 5 2.615 2 2.695 1
      注:— 表示该方程中没有这个参数;1、2、3分别表示马尾松、湿地松和杉木。
      Note: — indicates that there is no such parameter in the equation. 1, 2 and 3 represent Pinus massoniana, Pinus elliottii and Cunninghamia lanceolata, respectively.
    下载: 导出CSV

    表 6  基于参数b分级的立地质量等级哑变量参数

    Table 6.  Site quality grade dummy variable parameters based on parameter b classification

    模型
    Models
    因变量
    Dependent variables
    参数b哑变量估计值
    Estimated values of dummy variable of parameter b
    123456
    Model 5 胸径 DBH 0.007 7 0.012 6 0.018 0 0.024 9 0.035 3 0.052 7
    蓄积量 Volume 0.010 3 0.016 7 0.023 9 0.032 8 0.047 0 0.072 9
    生物量 Biomass 0.009 3 0.015 1 0.021 7 0.029 8 0.042 8 0.066 5
    Model 6 胸径 DBH 0.007 8 0.012 9 0.018 3 0.025 4 0.035 6 0.054 2
    蓄积量 Volume 0.010 6 0.017 2 0.024 6 0.033 8 0.048 3 0.076 6
    生物量 Biomass 0.009 7 0.015 6 0.022 5 0.031 0 0.044 3 0.070 6
    下载: 导出CSV
  • [1] 欧光龙, 胥 辉. 森林生物量模型研究综述[J]. 西南林业大学学报:自然科学, 2020, 40(1):1-11.

    [2] 曾伟生, 孙乡楠, 王六如, 等. 东北林区10种主要森林类型的蓄积量、生物量和碳储量模型研建[J]. 北京林业大学学报, 2021, 43(3):1-8. doi: 10.12171/j.1000-1522.20200058

    [3] 张雄清, 张建国, 段爱国. 基于单木水平和林分水平的杉木兼容性林分蓄积量模型[J]. 林业科学, 2014, 50(1):82-87.

    [4]

    Chave J, Condit R, Aguilar S, et al. Error propagation and scaling for tropical forest biomass estimates[J]. Philosophical Transactions of the Royal Society of London B:Biological Sciences, 2004, 359(1443): 40-420.
    [5] 罗云建, 张小全, 王效科, 等. 森林生物量的估算方法及其研究进展[J]. 林业科学, 2009, 45(8):129-134. doi: 10.3321/j.issn:1001-7488.2009.08.023

    [6] 董利虎, 李凤日. 三种林分生物量估算方法的比较[J]. 应用生态学报, 2016, 27(12):3862-3870.

    [7] 张会儒, 唐守正, 王奉瑜. 与材积兼容的生物量模型的建立及其估计方法研究[J]. 林业科学研究, 1999, 12(1):53-59. doi: 10.3321/j.issn:1001-1498.1999.01.010

    [8] 傅 煜, 雷渊才, 曾伟生. 单木生物量模型估计区域尺度生物量的不确定性[J]. 生态学报, 2015, 35(23):7738-7747.

    [9] 秦立厚, 张茂震, 钟世红, 等. 森林生物量估算中模型不确定性分析[J]. 生态学报, 2017, 37(23):7912-7919.

    [10] 杨惠滨, 国庆喜. 地形与竞争因子对红松胸径与年龄关系的影响[J]. 生态学报, 2016, 36(20):6487-6495.

    [11]

    Zeide B. Analysis of growth equations[J]. Forest Science, 1993, 39(3): 594-616. doi: 10.1093/forestscience/39.3.594
    [12]

    Takahashi K, Homma K, Vetrova V P, et al. Stand structure and regeneration in a Kamchatka mixed boreal forest[J]. Journal of Vegetation Science, 2001, 12(5): 627-634. doi: 10.2307/3236902
    [13]

    Stewart B J, Neily P D, Quigley E J, et al. Selected Nova Scotia old-growth forests: Age, ecology, structure, scoring[J]. The Forestry Chronicle, 2003, 79(3): 632-644. doi: 10.5558/tfc79632-3
    [14] 胡云云, 亢新刚, 赵俊卉. 长白山地区天然林林木年龄与胸径的变动关系[J]. 东北林业大学学报, 2009, 37(11):38-42. doi: 10.3969/j.issn.1000-5382.2009.11.013

    [15] 曹 磊. 基于多期保留木实测胸径估计吉林省蒙古栎天然林年龄[D]. 北京: 中国林业科学研究院, 2020.

    [16] 曹 磊, 刘晓彤, 李海奎, 等. 广东省常绿阔叶林生物量生长模型[J]. 林业科学研究, 2020, 33(5):61-67.

    [17] 龙时胜, 曾思齐, 甘世书, 等. 基于林木多期直径测定数据的异龄林年龄估计方法[J]. 中南林业科技大学学报, 2018, 38(9):1-8.

    [18]

    Lhotka J M, Loewenstein E F. An individual-tree diameter growth model for managed uneven-aged oak-shortleaf pine stands in the Ozark Highlands of Missouri, USA[J]. Forest Ecology and Management, 2011, 261(3): 770-778. doi: 10.1016/j.foreco.2010.12.008
    [19] 彭 娓, 李凤日, 董利虎. 黑龙江省长白落叶松人工林单木生长模型[J]. 南京林业大学学报:自然科学版, 2018, 42(3):19-27.

    [20]

    Mailly D, Turbis S, Pothier D. Predicting basal area increment in a spatially explicit, individual tree model: a test of competition measures with black spruce[J]. Revue Canadienne De Recherche Forestière, 2003, 33(3): 435-443.
    [21] 中华人民共和国农林部. 立木材积表 (LY208−77)[S]. 北京: 中国标准出版社, 1977.

    [22] 国家林业局. 立木生物量模型及碳计量参数——马尾松(LY/T 2263-2014 )[S]. 北京: 中国标准出版社, 2014.

    [23] 国家林业局. 立木生物量模型及碳计量参数——湿地松(LY/T 2261-2014 )[S]. 北京: 中国标准出版社, 2014.

    [24] 国家林业局. 立木生物量模型及碳计量参数——杉木(LY/T 2264-2014)[S]. 北京: 中国标准出版社, 2014.

    [25] 惠刚盈, 胡艳波, 赵中华, 等. 基于交角的林木竞争指数[J]. 林业科学, 2013, 49(6):68-73. doi: 10.11707/j.1001-7488.20130610

    [26] 曾伟生, 唐守正. 立木生物量方程的优度评价和精度分析[J]. 林业科学, 2011, 47(11):106-113. doi: 10.11707/j.1001-7488.20111117

    [27] 薛春泉, 徐期瑚, 林丽平, 等. 基于异速生长和理论生长方程的广东省木荷生物量动态预测[J]. 林业科学, 2019, 55(7):86-94. doi: 10.11707/j.1001-7488.20190709

    [28] 雷相东, 李希菲. 混交林生长模型研究进展[J]. 北京林业大学学报, 2003, 25(3):105-110. doi: 10.3321/j.issn:1000-1522.2003.03.022

    [29]

    Zeng W S, Zhang H R, Tang S Z. Using the dummy variable model approach to construct compatible single-tree biomass equations at different scales — a case study for Masson pine (Pinus massoniana) in southern China[J]. NRC Research Press, 2011, 41(7): 1547-1554.
    [30] 吴宏炜, 田 意, 黄光灿, 等. 基于非线性度量误差的湿地松生长模型[J]. 林业资源管理, 2019(6):71-76.

    [31]

    Biging G S, Dobbertin M. A comparision of distance-dependent competition measures for height and basal area growth of individual conifer trees[J]. For Sci, 1992, 38(3): 695-720.
    [32] 臧 颢, 刘洪生, 黄锦程, 等. 竞争和气候及其交互作用对杉木人工林胸径生长的影响[J]. 林业科学, 2021, 57(3):39-50.

    [33]

    Daniels R F, Burkhart H E, Clason T R. A comparison of competition measures for predicting growth of loblolly pine trees[J]. Canadian Journal of Forest Research, 1986, 16(6): 1230-1237. doi: 10.1139/x86-218
    [34]

    Kuehne C, Weiskittel A R, Waskiewicz J. Comparing performance of contrasting distance-independent and distance-dependent competition metrics in predicting individual tree diameter increment and survival within structurally-heterogeneous, mixed-species forests of Northeastern United States[J]. Forest Ecology and Management, 2019, 433: 205-216. doi: 10.1016/j.foreco.2018.11.002
    [35] 孟宪宇. 测树学[M]. 3 版. 北京: 中国林业出版社, 2006, 176, 290-291.

    [36] 赵 菡, 雷渊才, 符利勇. 江西省不同立地等级的马尾松林生物量估计和不确定性度量[J]. 林业科学, 2017, 53(8):81-93. doi: 10.11707/j.1001-7488.20170810

  • [1] 曹磊刘晓彤李海奎雷渊才 . 广东省常绿阔叶林生物量生长模型. 林业科学研究, 2020, 33(5): 61-67. doi: 10.13275/j.cnki.lykxyj.2020.05.008
    [2] 徐金良毛玉明郑成忠范荣德周世水陈永辉成向荣虞木奎 . 抚育间伐对杉木人工林生长及出材量的影响. 林业科学研究, 2014, 27(1): 99-107.
    [3] 陈少雄杨民胜王理平 . 尾叶桉造林密度与蓄积量、抗风和材性关系研究. 林业科学研究, 1998, 11(4): 435-438.
    [4] 张会儒唐守正王奉瑜 . 与材积兼容的生物量模型的建立及其估计方法研究. 林业科学研究, 1999, 12(1): 53-59.
    [5] 谭德远郭泉水王春玲马超 . 寄生植物肉苁蓉对寄主梭梭生长及生物量的影响研究. 林业科学研究, 2004, 17(4): 472-478.
    [6] 黄兴召孙晓梅张守攻陈东升 . 辽东山区日本落叶松生物量相容性模型的研究. 林业科学研究, 2014, 27(2): 142-148.
    [7] 曾伟生 . 我国3种针叶林的材积源生物量模型研建. 林业科学研究, 2021, 34(4): 49-57. doi: 10.13275/j.cnki.lykxyj.2021.04.006
    [8] 朱春全王世绩王富国张启钮利民袁国文 . 六个杨树无性系苗木生长、生物量和光合作用的研究*. 林业科学研究, 1995, 8(4): 388-394.
    [9] 周再知郑海水尹光天杨曾奖陈康泰 . 橡胶树生物量估测的数学模型*. 林业科学研究, 1995, 8(6): 624-629.
    [10] 李意德曾庆波吴仲民周光益陈步峰 . 我国热带天然林植被C贮存量的估算. 林业科学研究, 1998, 11(2): 156-162.
    [11] 陆新育陈绍信李淼泉常显明 . 泡桐生物量的研究*. 林业科学研究, 1990, 3(5): 422-426.
    [12] 郑松发郑德璋廖宝文李云 . 广东省潮间带污染与红树林造林. 林业科学研究, 1997, 10(6): 639-646.
    [13] 贾小容苏志尧 . 广东省山茶科物种丰富度水平地带性格局分析. 林业科学研究, 2008, 21(4): 469-474.
    [14] 廖宝文郑德璋郑松发 . 海桑林生物量的研究*. 林业科学研究, 1990, 3(1): 47-54.
    [15] 王翰琛张雄清张建国屈彦成姜丽 . 杉木人工林不同密度间伐林分生长优势的变化规律分析. 林业科学研究, 2021, 34(5): 32-38. doi: 10.13275/j.cnki.lykxyj.2021.005.004
    [16] 何贵平陈益泰胡炳堂封剑文刘化桐蔡宏明 . 杉木与马褂木、檫树混交林及其纯林生物量和土壤肥力研究. 林业科学研究, 2001, 14(5): 540-547.
    [17] 廖宝文郑德璋郑松发陈步峰 . 木榄林生物量的灰色动态预测. 林业科学研究, 1991, 4(4): 360-367.
    [18] 周国模吴家胜应叶青姚建祥 . 喜树种源苗期生物量研究. 林业科学研究, 1999, 12(4): 386-391.
    [19] 陈建新王明怀殷祚云林军陈木发梁胜耀 . 广东省秃杉引种栽培效果及栽培区划分研究. 林业科学研究, 2002, 15(4): 399-405.
    [20] . 园林绿化树种灰色综合评价与分级选择的研究. 林业科学研究, 2009, 22(3): -.
  • 加载中
图(4) / 表(6)
计量
  • 文章访问数:  3659
  • HTML全文浏览量:  2787
  • PDF下载量:  96
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-10-21
  • 录用日期:  2021-12-31
  • 网络出版日期:  2022-03-11
  • 刊出日期:  2022-06-20

广东省针叶树种蓄积量和生物量生长模型研究

    通讯作者: 李海奎, lihk@ifrit.ac.cn
    作者简介: 黄金金,硕士研究生。主要研究方向:林业统计与生物数学模型。E-mail:Huang202098@163.com
  • 1. 中国林业科学研究院资源信息研究所,北京 100091
  • 2. 国家林业和草原局森林经营与生长模拟重点实验室,北京 100091

摘要:  目的 定量化表征森林蓄积量和生物量与年龄的关系,为预测森林蓄积量和生物量提供依据。 方法 基于广东省3个针叶树种5期连清数据,通过保留木林分平均胸径和林龄之间的关系,分别以生长潜力和生长速度分级、是否考虑竞争、分步和联合建模共8种组合,构建林分蓄积量和生物量生长模型,并评价模型拟合优度;以独立的4期连清数据为验证样本,评价模型的适用性。 结果 (1)不同种模型系确定系数皆在0.980以上,平均预估误差在±1%以内,总体相对误差在±0.5%以内。综合模型拟合和预测效果,以生长速度分级、不考虑竞争的分步模型最优,竞争对模型的适用性影响不大。(2)最优模型蓄积量和生物量的估计误差最大分别为10.36%和10.22%,模型适用性较好,4期估计误差表现为中期高于首末两期。(3)马尾松生长潜力最大,杉木最小,对杉木的估计效果优于马尾松和湿地松。(4)立地质量等级越高,生长量的极值越大,达到最高峰所需年限也更短;同一立地质量下马尾松生长量最大,其次为湿地松,杉木最小。 结论 含参数分级和林分特征的蓄积量和生物量模型,可以反映立地质量对广东省针叶树种蓄积量和生物量预测的影响,为精确估测森林储量提供方法学支持,也为其他地区林分水平生长模型的构建提供借鉴。

English Abstract

  • 森林蓄积量和生物量是反映森林资源水平的基本指标,精准估测森林蓄积量和生物量,可以为林业资源管理决策提供科学依据[1]。国内外关于蓄积量和生物量模型的研建主要集中于单木水平,而林分水平的较少,尚未形成我国林业行业标准[2]。林分水平的生长模型能直接预测林分蓄积量和生物量,减少单木的误差积累[3-4]。目前估计区域尺度蓄积量和生物量的建模方法有多种,应用较多的一是以异速生长方程为基础,根据林分总胸高断面积、林分优势高、林分平均高等估计林分以及区域的地上生物量[5-6];二是生物量转换因子法,张会儒等[7]将全国森林类型按照林分优势树种归并为9类,构建了全国主要树种组林木与材积相容的生物量回归模型;三是通过单木模型推算,傅煜等[8]通过Monte Carlo法反复模拟由单木生物量模型推算区域尺度地上生物量,估计了江西省杉木地上总生物量,秦立厚等[9]分析量化单木生物量模型中各种不确定性来源,估计区域地上生物量。

    大量研究表明,林木胸径与林木年龄之间存在紧密的生长关系,为建立有关的胸径生长模型奠定了基础[10-14]。曹磊[15-16]和龙时胜等[17]采用Richards生长方程拟合树木直径与年龄的生长关系,模型的拟合优度较高,其中龙时胜等[17]利用模型估计林分平均年龄时,拟合年龄与真实年龄的绝对误差最大值为2年,相对误差最大值为8.2%,两者间的差异性较小,说明利用林木多期直径估计异龄林林分平均年龄具有较高的可行性。大量研究[18-20]认为,林分蓄积量、生物量生长主要受树木自身大小、竞争因子和立地条件的影响。

    本研究以广东省3个针叶树种155块样地为对象,利用5期连清数据,基于林分平均胸径与年龄的生长关系估计样地林龄,进而建立林分水平的蓄积量和生物量生长模型,以体现生长潜力的参数a分级和与生长速度有关的参数b分级、考虑林分竞争与否、分步建模(一元非线性回归)和联合建模(非线性误差变量联立方程组)共8种组合,构建模型系,评价模型的拟合优度和预测效果,筛选最优模型,以期实现对广东省针叶林蓄积量和生物量的准确估计。

    • 广东省位于20°13′~25°31′ N和109°39′~117°19′ E,处于东亚季风区,包括了热带、中亚和南亚热带3种气候,年均气温22.3℃,月降水量约313 mm,年日照时数约1 746 h,太阳辐射总量约4 200~5 400 MJ·m−2,生长期水热充足,利于植被生长。地貌类型复杂多样,以低山、丘陵为主,地势北高南低。

    • 建模数据:从广东省第5—9次连清数据中,选择马尾松(Pinus massoniana Lamb.)、湿地松(Pinus elliottii Engelm.)、杉木(Cunninghamia lanceolata (Lamb.) Hook.)3个针叶树种,筛选5次清查中皆存在且保留木大于5株的样地,共155块(马尾松林91块、湿地松林16块、杉木林48块),样地分布见图1。固定样地为0.0667 hm2,样木调查因子包括样木号、胸径,单株材积按广东省地方一元材积表[21]计算获得,单株生物量按行业标准生物量模型[22-24]获得。为避免株数差异造成样地蓄积量和生物量的差异,本研究中样地蓄积量和生物量为样地内5期保留木的蓄积量和生物量算数平均值,统计量如表1所示。

      图  1  建模数据固定样地分布示意图

      Figure 1.  Distribution map of fixed sample plots of modeling data

      表 1  建模数据统计量

      Table 1.  Modeling data statistics

      树种 Tree species  统计量 Statistics  1997年2002年2007年2012年2017年
      马尾松 Pinus massoniana 胸径 DBH/cm 10.2±3.2 13.4±3.9 15.6±4.2 17.8±4.5 20.4±5.0
      蓄积量 Volume/dm3 43.7±51.3 83.0±75.0 120.5±92.8 164.3±112.6 228.3±143.1
      生物量 Biomass/kg 42.8±46.2 79.7±69.1 113.9±85.8 154.0±105.1 213.3±134.8
      湿地松
      Pinus elliottii
      胸径 DBH/cm 8.2±2.0 12.1±2.5 14.7±3.1 17.0±3.5 19.2±4.1
      蓄积量 Volume/dm3 19.2±13.5 54.5±30.1 91.1±47.6 132.5±66.7 183.7±99.6
      生物量 Biomass/kg 21.5±12.8 53.8±27.2 85.0±40.9 120.5±56.9 164.3±85.2
      杉木
      Cunninghamia lanceolata
      胸径 DBH/cm 8.8±2.8 11.7±3.3 13.5±3.6 15.0±3.7 16.7±4.1
      蓄积量 Volume/dm3 33.1±42.4 66.4±66.5 93.8±82.6 121.2±91.2 159.0±112.7
      生物量 Biomass/kg 22.7±25.1 43.0±38.7 59.2±47.9 75.2±52.7 97.2±65.0
        注:数据为平均值±标准差。
        Note: The data is average ± standard deviation.

      样地胸径采用林分平均胸径:

      $ {D}_{g}=\sqrt{\frac{1}{n}{\sum }_{i=1}^{n}{d}_{i}^{2}} $

      (1)

      式中,$ {D}_{g} $为林分平均胸径,n为样木株数,di为样木胸径。

      验证数据:选取2002年达到起测胸径、存在于4次连清中,且保留木大于5株的样地,共126块,其中马尾松59块、湿地松22块、杉木45块。各统计量如表2所示。

      表 2  验证数据统计量

      Table 2.  Verification data statistics

      树种 Tree species  统计量 Statistics2002年2007年2012年2017年
      马尾松 Pinus massoniana 胸径 DBH/cm 9.8±3.5 12.4±4.0 14.9±4.6 18.0±5.4
      蓄积量 Volume/dm3 37.1±43.0 65.6±61.1 102.2±84.9 162.6±120.7
      生物量 Biomass/kg 37.3±39.4 64.0±55.8 97.2±77.9 152.5±111.0
      湿地松 Pinus elliottii 胸径 DBH/cm 8.4±1.7 12.1±2.8 15.4±3.6 18.8±4.4
      蓄积量 Volume/dm3 19.2±11.9 54.1±33.5 101.8±56.5 169.1±94.5
      生物量 Biomass/kg 21.2±11.3 54.1±30.2 94.2±48.6 151.8±80.9
      杉木 Cunninghamia lanceolata 胸径 DBH/cm 8.8±2.6 11.2±2.9 13.2±3.2 15.2±3.3
      蓄积量 Volume/dm3 29.9±26.8 54.7±42.5 81.4±57.4 116.1±73.3
      生物量 Biomass/kg 20.7±16.4 36.0±25.1 51.9±33.5 72.3±42.4
        注:数据为平均值±标准差。
        Note: The data is average ± standard deviation.
    • 采用Mitscherlich理论生长方程,拟合5期林分平均胸径和相应时间间隔之间的关系;以Richards理论生长方程,拟合样地平均蓄积量和生物量与相应时间间隔之间的关系。

      $ \begin{split} & {D}_{gij}=a1*(1-{e}^{-b1*({age}_{i} + {ad}_{j})}) (i=\mathrm{1,2},\cdots ,155;\\ & j=\mathrm{1,2},\cdots ,5) \end{split}$

      (2)

      $ \begin{split} & {V}_{ij}=a2*(1-{e}^{-b2*({age}_{i} + {ad}_{j})}{)}^{c2} (i=\mathrm{1,2},\cdots ,155;\\ & j=\mathrm{1,2},\cdots ,5) \end{split} $

      (3)

      $ \begin{split} & {BIO}_{ij}=a3*(1-{e}^{-b3*({age}_{i} + {ad}_{j})}{)}^{c3} (i=\mathrm{1,2},\cdots ,155;\\ & j=\mathrm{1,2},\cdots ,5) \end{split} $

      (4)

      式中:$ {D}_{gij} $$ {V}_{ij} $$ {BIO}_{ij} $分别为第i个样地第j期的样地平均胸径、蓄积量和生物量,$ {age}_{i} $是第$ i $个样地的初期林龄,$ {ad}_{j} $为第$ j $期与初期(1997年)的时间间隔,$ a1、a2、a3、b1、b2、b3、c2、c3 $为参数。

      基于上述胸径、蓄积量和生物量3个生长方程,分别以体现生长潜力的参数a分级和与生长速度有关的参数b分级、考虑竞争与否、分步或联合建模等8种组合,构建模型系(表3)。

      表 3  8种模型系

      Table 3.  Eight model systems

      模型 Models方程 Equation模型 Models方程 Equation
      Model 1 $ {D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*\left({age}_{i} + {ad}_{j}\right)}) $
      $ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
      $ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
      Model 5 $ {D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)}) $
      $ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
      $ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
      Model 2 $ {D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)}) $
      $ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
      $ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
      Model 6 $ {D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)}) $
      $ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}} $
      $ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}} $
      Model 3 $ \left\{{array}{c}{D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $ Model 7 $ \left\{{array}{c}{D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $
      Model 4 $ \left\{{array}{c}{D}_{gij}={a1}_{m}*(1-{e}^{-{b1}_{k}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{m}*(1-{e}^{-{b2}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{m}*(1-{e}^{-{b3}_{k}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $ Model 8 $ \left\{{array}{c}{D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)})\\ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}{array}\right. $
    • 对于林分水平模型来说, 与树冠有关的竞争指数能较好地反映树木间的竞争关系,因此选用基于交角且对象木与竞争木距离小于10 m的竞争指数$ ju10 $,详见惠刚盈等[25]。由于林木竞争对单木胸径的生长具有较大的影响,但对林分平均胸径的影响并不算太大,为解释林木竞争对林分生长的影响,故将竞争指数引入了胸径生长方程:

      $ {D}_{gij}=a1*(1-{e}^{-b1*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}*\left({age}_{i} + {ad}_{j}\right)}) $

      (5)

      式中:$ {b}_{0} $为竞争指数参数,$ {ju10}_{ij} $为第i个样地第j期的竞争指数(样地竞争指数为样地内所有保留木竞争指数的算数平均值),其余同上。

    • 为了兼顾区域尺度上立地质量的多样性以提高模型在大尺度范围内的适用性,对参数ab分级代替立地质量的差异,以参数a分级为例:

      将式(2)作为立地质量划分等级的分级方程,已知林分平均胸径$ {D}_{gij} $和时间间隔$ {ad}_{j} $,以样地号和树种号构成的唯一编号作为参数aage的哑变量,树种作为参数b的哑变量,拟合方程(6)。将参数${a1}_{i}$从低到高排序,并按照有序样本聚类分为6类。

      $ {D}_{gij}={a}_{i}*(1-{e}^{-{b1}_{k}*\left({age}_{i} + {ad}_{j}\right)}) $

      (6)

      式中: $ {a1}_{i}(i=\mathrm{1,2},\dots ,155) $是第i个样地的参数,$ {b1}_{k}(k=\mathrm{1,2},3 $时分别代表马尾松、湿地松和杉木)是第$ k $个树种的参数,其余同上。

    • 以Model 5为例,基于式(2),以参数a区别树种,以参数b区别立地质量差异,以唯一编号作为初期林龄age(未知)的哑变量,根据5期林分平均胸径$ {D}_{gij} $和时间间隔$ {ad}_{j} $之间的关系来估计样地初期林龄$ {age}_{i} $,见式(7)。将$ {age}_{i} $作为式(3)和(4)的自变量,以参数c区别树种,其余参数处理同上,分别拟合样地平均蓄积量$ {V}_{ij} $和生物量$ {BIO}_{ij} $与林龄$ {age}_{i} $的关系,即式(8)和(9)。

      $ {D}_{gij}={a1}_{k}*(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)}) $

      (7)

      $ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*(ag{e}_{i} + {ad}_{j})}{)}^{{c2}_{k}} $

      (8)

      $ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*(ag{e}_{i} + {ad}_{j})}{)}^{{c3}_{k}} $

      (9)

      式(7)~(9)中:$ {a1}_{k}、{a2}_{k}、{a3}_{k}\left(k=\mathrm{1,2},3\right) $分别为胸径、蓄积量和生物量生长方程中的第k个树种的参数,$ {b1}_{m}、{b2}_{m}、{b3}_{m}\left(m=\mathrm{1,2},\cdots,6\right) $分别为胸径、蓄积量和生物量生长方程中的第m个立地质量等级的参数,$ {c2}_{k}、{c3}_{k}\left(k=\mathrm{1,2},3\right) $分别为蓄积量和生物量生长方程中的第k个树种的参数,其余同上。

    • 为了比较建模方法对生长模型拟合优度和估测精度的差异,除分步建模法外,本研究还采用了联合建模法,以Model 7为例。由于林龄未知,且3个方程中年龄参数$ {age}_{i} $是一致的,基于式(2)~(4),参数处理同分步建模法,根据林分平均胸径、蓄积量和生物量与时间间隔之间的关系,采用非线性误差变量联立方程组法同时拟合3个方程,即式(10)。

      $ \left\{\begin{array}{c}{D}_{gij}={a1}_{k}*\left(1-{e}^{-{b1}_{m}*\left({age}_{i} + {ad}_{j}\right)}\right)\\ {V}_{ij}={a2}_{k}*(1-{e}^{-{b2}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c2}_{k}}\\ {BIO}_{ij}={a3}_{k}*(1-{e}^{-{b3}_{m}*\left({age}_{i} + {ad}_{j}\right)}{)}^{{c3}_{k}}\end{array} \right. $

      (10)
    • 考虑各模型参数不同,为了检验8种模型系的拟合效果,采用以下4项模型评价指标:确定系数$ {R}^{2} $、估计值标准误差SEE、平均预估误差MPE、总体相对误差TRE,公式详见文献[26]

    • (1)确定立地质量等级和林龄(以参数b分级为例)

      根据胸径生长方程反推林龄计算公式,如下:

      $ \begin{split} & pag{e}_{ij}=\frac{-\mathrm{ln}(1-\frac{{D}_{gij}}{{a1}_{k}})}{{b1}_{m}*{\left({e}^{{ju10}_{ij}}\right)}^{{b}_{0}}} \\ & (i=\mathrm{1,2},\cdots ,126;\;j=\mathrm{1,2},\cdots ,4) \end{split}$

      (11)

      式中,$ pag{e}_{ij} $为第i个样地第j期的样地林龄,当不引入竞争指数时$ {b}_{0}=0 $,其余同上。

      由于验证样本所在样地的立地质量等级m和初期林龄皆未知,需先确定立地质量等级,再由式(11)求出林龄,才可估计蓄积量和生物量。当立地质量分级时,至少需要两期胸径才能确定样地的立地质量等级,薛春泉等[27]发现长周期确定的等级比短周期稳定,因此本研究采用期初(2002年)和期末(2017年)2期胸径来确定样地的立地质量等级,确定方法如下:

      $ {OSC}_{i}=\mathrm{m}\mathrm{i}\mathrm{n}|pag{e}_{i4}-pag{e}_{i1}-15| $

      (12)

      式中,$ pag{e}_{i1}\mathrm{和}pag{e}_{i4} $为第i个样地期初和期末的样地林龄,该方法的直观解释就是使推算的2002年和2017年林龄间隔与实际林龄间隔(15年)相差最小的立地质量等级,即为该样地的立地质量等级。

      根据式(12)确定样地的立地质量等级后可知参数$ {b}_{m} $,再由式(11)计算初期(2002年)林龄age,则2007年、2012年、2017年时样地的林龄分别由初期林龄相应加上5、10、15年。

      (2)模型估测效果评价指标

      由蓄积量和生物量生长方程估计样地平均蓄积量和生物量,样地总蓄积量和总生物量可由平均蓄积量和生物量乘以株数获得,按式(13)~(14)计算各树种4个时期下区域尺度蓄积量和生物量估计的总相对误差EVEBIO

      $ EV=\frac{{\sum }_{i=1}^{n}\left(P{V}_{i}-{V}_{i}\right)*{N}_{i}}{{\sum }_{i=1}^{n}{V}_{i}*{N}_{i}} $

      (13)

      $ EBIO=\frac{{\sum }_{i=1}^{n}(PBI{O}_{i}-BI{O}_{i})*{N}_{i}}{{\sum }_{i=1}^{n}BI{O}_{i}*{N}_{i}} $

      (14)

      式中:n为样地数,$ {N}_{i} $为第i个样地的保留木株数,$ P{V}_{i} $$ {V}_{i} $分别为第i个样地平均蓄积量的估计值和真实值,$ PBI{O}_{i} $$ BI{O}_{i} $分别为第i个样地平均生物量的估计值和真实值。

    • 连年生长量函数$ Z\left(t\right) $为总生长过程曲线$ {y}\left({t}\right) $的一阶导数,为树木在某年的实际生长速度,即连年生长量随年龄t变化的函数。Richards生长方程的连年生长量方程为:

      $ Z\left(t\right)=\frac{dy\left(t\right)}{t}=a\cdot b\cdot c\cdot (1-{e}^{-b\cdot t}{)}^{c-1}\cdot {e}^{-b\cdot t} $

      (15)

      平均生长量函数$ \theta \left(t\right) $用来表示树木在某一时刻的平均生长速度,总生长过程曲线$ y\left(t\right) $除以年龄t则可得平均生长量依年龄变化的函数。Richards生长方程的平均生长量方程为:

      $ \theta \left(t\right)=\frac{y\left(t\right)}{t}=\frac{a}{t}\cdot (1-{e}^{-b\cdot t}{)}^{c} $

      (16)
    • 表4列出了8种模型系的拟合评价指标。整体上看,8种模型系的3个生长方程拟合确定系数$ {R}^{2} $皆达0.980以上,总相对误差TRE在±0.5%以内,平均预估误差MPE在±1%以内,拟合效果良好,但各模型间的4项指标差异不大。当控制建模方法和立地质量分级的参数相同、仅比较竞争影响时,引入竞争指数对模型拟合优度的提升没有明显的作用,如Model 5和Model 6拟合指标间差异甚微、Model 2较Model 1的蓄积量和生物量生长方程的拟合效果甚至出现了一定程度的下降;类似地,可以发现参数b分级模型的拟合优度较参数a分级模型稍好;当控制竞争影响和分级参数相同时,相对于分步建模,采用含度量误差的联立方程组时,蓄积量和生物量生长方程的拟合效果有所提高,但是胸径生长方程的拟合效果稍差。

      表 4  模型拟合评价指标

      Table 4.  Model fitting evaluation index

      模型
      Models
      胸径 DBH 蓄积量 Volume 生物量 Biomass
      R2SEE/cmTRE/%MPE/% R2SEE/dm3TRE/%MPE/% R2SEE/kgTRE/%MPE/%
      Model 1 0.984 0.77 0.00 0.36 0.985 13.42 −0.54 0.83 0.987 10.97 −0.08 0.79
      Model 2 0.985 0.75 0.00 0.35 0.982 14.54 −0.71 0.90 0.985 11.93 −0.27 0.85
      Model 3 0.981 0.74 −0.35 0.35 0.986 12.74 −0.90 0.79 0.989 10.32 −0.48 0.74
      Model 4 0.981 0.75 −0.20 0.35 0.986 12.67 −0.64 0.78 0.989 10.28 −0.14 0.74
      Model 5 0.985 0.74 0.00 0.35 0.986 12.83 −0.11 0.79 0.988 10.84 0.11 0.78
      Model 6 0.985 0.74 0.00 0.35 0.986 13.01 −0.12 0.80 0.987 11.02 0.10 0.79
      Model 7 0.981 0.74 −0.01 0.35 0.987 12.27 0.15 0.76 0.989 10.26 0.31 0.73
      Model 8 0.982 0.73 −0.04 0.34 0.988 12.15 0.00 0.75 0.989 10.17 0.17 0.73
    • 8种蓄积量和生物量生长方程的模型验证指标分别如图2图3所示。总体而言,各模型的预测效果都较好,从图2可知,区域尺度蓄积量估计误差最大的仅有17.62%,说明蓄积量估计精度可达82.38%以上;从图3可知,区域尺度生物量估计误差最大的仅有12.51%,则生物量估计精度可达87.49%以上。8个模型中以基于参数b分级的分步模型(即Model 5和Model 6)预测效果最优,其对3个树种4期蓄积量和生物量的估计误差较其余模型稍小,且在不同树种和不同时期的表现都较为稳定,说明模型适用性较好。以Model 5为例,3个树种中对杉木的估计效果稍好,估计误差最大为−5.25%,而估计马尾松和湿地松时误差最大为10.00%和10.36%。8个模型对3个针叶树种的4期区域尺度蓄积量和生物量估计的误差,大都表现出中期误差高于首末2期,初期估计精度较高,但明显出现低估的现象,伴随着样地林龄的增大,中期出现高估且精度稍低,但到末期(2017年)时估计误差降低,精度有所提高。

      图  2  8种蓄积量生长模型的区域尺度估计误差

      Figure 2.  Regional-scale estimation errors of eight volume growth models

      图  3  8种生物量生长模型的区域尺度估计误差

      Figure 3.  Regional-scale estimation errors of eight biomass growth models

    • 综合8个模型的拟合效果(表4)和验证效果(图2图3),Model 5最优,Model 6次之。表5表6列出了基于参数b分级的分步模型的模型参数估计值。由表5可知,马尾松胸径、蓄积量和生物量的生长潜力均为最大,杉木最小,其中马尾松蓄积量和生物量的生长极限值可达湿地松和杉木的2倍以上。由表6可知,立地质量等级越高,树木的生长速度越快,但是各等级下的林分蓄积量生长速度均为最快,其次为生物量,胸径生长速度最慢。

      表 5  树种哑变量参数和竞争指数参数

      Table 5.  Dummy variable parameters of tree species and competition index parameters

      模型
      Models
      因变量
      Dependent variables
      参数a
      Parameter a
      参数c
      Parameter c
      b0
      123123
      Model 5 胸径 DBH 58.15 41.57 36.59
      蓄积量 Volume 1 827.54 853.40 761.84 2.709 1 2.857 1 2.926 5
      生物量 Biomass 1 892.76 760.24 460.70 2.596 8 2.546 1 2.593 4
      Model 6 胸径 DBH 57.97 41.75 36.44 −0.281 0
      蓄积量 Volume 1 745.73 836.66 745.45 2.738 6 2.924 5 3.030 5
      生物量 Biomass 1 784.94 739.66 447.88 2.631 5 2.615 2 2.695 1
        注:— 表示该方程中没有这个参数;1、2、3分别表示马尾松、湿地松和杉木。
        Note: — indicates that there is no such parameter in the equation. 1, 2 and 3 represent Pinus massoniana, Pinus elliottii and Cunninghamia lanceolata, respectively.

      表 6  基于参数b分级的立地质量等级哑变量参数

      Table 6.  Site quality grade dummy variable parameters based on parameter b classification

      模型
      Models
      因变量
      Dependent variables
      参数b哑变量估计值
      Estimated values of dummy variable of parameter b
      123456
      Model 5 胸径 DBH 0.007 7 0.012 6 0.018 0 0.024 9 0.035 3 0.052 7
      蓄积量 Volume 0.010 3 0.016 7 0.023 9 0.032 8 0.047 0 0.072 9
      生物量 Biomass 0.009 3 0.015 1 0.021 7 0.029 8 0.042 8 0.066 5
      Model 6 胸径 DBH 0.007 8 0.012 9 0.018 3 0.025 4 0.035 6 0.054 2
      蓄积量 Volume 0.010 6 0.017 2 0.024 6 0.033 8 0.048 3 0.076 6
      生物量 Biomass 0.009 7 0.015 6 0.022 5 0.031 0 0.044 3 0.070 6
    • 基于参数b分级不考虑竞争的分步建模(即Model 5)最优,为节约篇幅,选择样地立地质量等级占比较多的等级2、3和4,绘制Model 5下蓄积量、生物量的连年生长量和平均生长量的关系曲线图(图4)。蓄积量和生物量的连年生长量、平均生长量均显示出立地质量越好,生长量能达到的极限值越大,且所需年限更短,各相邻等级间的差异随着等级的增加而逐渐加大。对于不同的树种来说,同一立地质量等级下的生长量均为马尾松>湿地松>杉木,生长量达到最大时的林龄为杉木>湿地松>马尾松,各树种间生长量差异随着等级的增加而逐渐加大。立地质量等级越高,平均生长量和连年生长量曲线相交时的样地年龄越小。

      图  4  等级2~4的连年生长量和平均生长量

      Figure 4.  Annual increment and average increment from class 2 to class 4

    • 处理多树种问题是建立林分模型的一个难点,通常是将组成树种进行分组或合并[28]。曾伟生等[2]建立了东北林区林分水平的三储量模型,筛选出10种主要森林类型并引入哑变量体现差异后,蓄积量模型、生物量模型和碳储量模型的确定系数分别从0.945、0.805、0.839提高到0.959、0.949、0.951。本研究对存在多个组成树种的样地按优势树种分属不同的样地,并构建哑变量区别广东省的3个主要针叶树种,从表5图4可知,各树种间的生长潜力和生长量有较大的差异,表现为马尾松>湿地松>杉木,其中马尾松远超2倍左右,这是因为建模数据(表1)中马尾松保留木大多为天然林,而杉木的起源大多为人工林,Zeng等[29]研究结果表明同一直径的单木地上生物量估计值在不同起源上存在一定差异,天然林的生物量大于人工林。此外,在同一立地质量下,不同树种间生长量的差异主要是由于树种特性的差异引起的,不同树种对同一区域的生长适应程度并不一致,马尾松可能更适宜广东的立地条件。各树种间区域尺度蓄积量和生物量的预测精度也有较大的差异,对杉木的估计效果较马尾松和湿地松稍好,可能原因是不同树种的建模样本量和检验样本量不一致,另一个可能的原因是不同的树种有不同最优的基础生长模型,如吴宏炜等[30]研究发现,Schumacher模型是湿地松林分蓄积量模型的最优基础模型。对于估计误差,大都表现出中期高于首末两期的现象,一方面可能是由于采用首末两期胸径数据来确定验证样本所在样地的立地质量等级,相较于中期,推算的各样地立地质量等级更贴合期初和期末时的实际情况,因此用于林龄确定、估计蓄积量和生物量时,3个生长方程中反映立地质量差异的参数更为稳定,模型的预测效果会更好。另一方面随着林龄增大,蓄积量和生物量逐渐积累,公式(13)和(14)分母变大,这可能使后期的估计误差降低,估计效果更好。

      非线性独立回归估计法的估计效果优于含度量误差变量联立方程组法,因为分步建模时,3个方程都是独立拟合给出参数的最优估计,而采用联立方程组法联立求解时,由于方程中林龄一致的约束,为保证方程之间的联系,得到更优化的参数估计,各方程较独立拟合时都有了不同的变化。从拟合评价指标(表4)可以看出,对比分步建模法,采用联立方程组法可以提高蓄积量和生物量生长方程的拟合效果,但是却降低了胸径生长方程的拟合优度,导致林分年龄的估计精度有所下降,使用验证数据集检验模型的适用性时,由于林分年龄未知,确定首期林龄成为关键的一步,所以胸径生长方程的拟合效果直接影响了验证时的估计精度。由于林分大部分为混交林,具有年龄不一、树种结构关系复杂等特点,因此探究林木间竞争作用对林分结构和生长变化的影响至关重要。一般而言,模型中的参数个数越多,拟合效果会更优,但是引入竞争指数后对于模型的拟合效果(表4)和估计精度(图23)没有明显的提高。虽然从理论上讲包含树木空间位置的竞争指数对预测林木生长更有效, 但Biging等[31]研究发现与距离无关的竞争指数比与距离有关的竞争指数预测单木生长更好,臧灏等[32]探究人工杉木林胸径生长对竞争的响应时也发现对比于与距离有关的单木竞争指标,采用与距离无关的单木竞争指标的生长模型估计精度更高。竞争指数表现的预测能力依赖于树种、立地条件等因素,这些因素的变化会导致不同类型的竞争指数表现不一致,且不具有决定性[33-34]。本研究引入竞争指数对林分模型的适用性没有明显影响,可能与选取的竞争指数有关,也有可能和立地的限制有关,因为在较大范围内密度对林分蓄积量、生物量的影响不大[35]

      模型的应用范围和模型的预测精度相关联,当模型应用于大尺度的森林资源测算时,若忽略立地质量的差异必然导致模型估计结果的不准确。立地质量是指在某一立地上既定森林或其他植被类型的生产潜力[35],既定立地上的不同树种和相同的树种在不同的立地上的立地质量都是有差异的。本研究为体现广东省样地间的立地质量多样性,分别对与生长潜力有关的参数a和与生长速度有关的参数b进行分级,结果表明以参数b分级体现立地质量的差异优于参数a分级。薛春泉等[26]研究发现生长速度分级的胸径生长方程与地上生物量的异速生长方程构成的模型系,不仅拟合效果达到了0.95以上,区域尺度生物量的估计误差也下降到了8%以内,这与本研究参数b分级模型预测效果较好具有一致性。另外,林分蓄积量和生物量的生长量随着立地质量的提高而增加,说明立地质量等级越高,森林蓄积量和生物量密度越大。立地质量与森林生产力相关,不同的立地质量代表着不同的生产力,立地质量越高,林分的生长速率越快,达到最大生长量所需的年限也更短。赵菡等[36]用优势木树高等级代替立地等级估计江西省不同立地质量等级的马尾松林生物量,研究发现不同立地等级的地上生物量均值估计结果随着立地质量等级的升高而增大,这与本研究结果一致。

    • 本研究以广东省3种针叶优势树种为研究对象,利用1997—2017年的5期森林资源连续清查数据,基于林分平均胸径、蓄积量和生物量与相应调查时间间隔的关系建立林分生长模型,实现了对广东省蓄积量和生物量的准确估计。本研究结论如下:(1)基于与生长速度有关的参数b分级不考虑竞争的分步模型最优,引入竞争指数与否差异不大。(2)最优模型的区域尺度蓄积量和生物量估计误差最大分别为10.36%和10.22%,模型的适用性较好,4期的估计误差表现为中期高于首末两期。(3)不同树种中以马尾松的生长潜力最大,杉木最小,但是模型对杉木的估计效果较马尾松和湿地松稍好。(4)立地质量等级越高,生长量极值越大,达到最高峰所需年限也更短;同一立地质量下马尾松的生长量最大,其次为湿地松,杉木最小。

参考文献 (36)

目录

    /

    返回文章
    返回