-
兴安落叶松(Larix gmelinii (Rupr.) Kuzen)是松科落叶松属的落叶乔木树种,具有耐寒耐湿的特点。兴安落叶松林可促进改善东北地区生态环境,对我国东北地区森林可持续发展有着不可替代的作用[1],且对全球碳平衡有重要影响[2]。针对大兴安岭地区的天然落叶松林,多位专家学者已有一定程度的研究积累,主要有:单木层次上,有单木的生长规律研究[3],单木直径生长模型构建[4],地上生物量模型构建[5]等;林分层次上,直径与树高分布研究[6]、结构特征研究[7]、材积源生物量模型研究[8]等。上述经验模型是森林生长和收获的主要模拟手段,然而其对环境变化、人为干扰下天然落叶松林生长和产量预测具有局限性[9]。与经验模型相比,过程模型可以充分考虑林木生长生理过程的相互作用与反馈机制,考虑不同气温、养分、降水对林木生长的影响。然而,到目前为止,在东北地区应用过程模型模拟天然林林分生长的研究仍十分有限。3-PG模型(Physiological Processes Predicting Growth)是由Landsberg和Waring[10]创建的一种基于过程的生长预测模型,通过碳吸收与固定、生物量分配与损耗以及土壤水分动态平衡3个模块,构建方程来实现动态模拟,可以充分发挥过程作为过程模型的优势。目前该模型已被多个国家和地区广泛用于人工林经营管理研究[11],国内也有学者在桉树(Eucalyptus robusta Smith)[12]人工林、马尾松(Pinus massoniana Lamb.)人工林[13]进行了参数化与生长预测,结果较为可靠,但尚未见运用于天然林的报道。本研究以天然兴安落叶松林为研究对象,运用3-PG模型对胸径、树高、整株与各组分生物量进行了拟合,探索该模型在天然林中应用的可行性,为制定相关森林经营计划提供参考。
-
本研究在我国东北地区的大兴安岭山脉开展。大兴安岭是中国唯一的寒温带明亮针叶林区和仅有的寒温带生物基因库,其森林覆盖率高达 79.83%[14],属寒温带大陆性季风气候,冬长夏短,全年霜冻期长达180 d,土壤类型主要为棕色针叶林土与暗棕壤,树种主要以兴安落叶松为主,其次为白桦(Betula platyphylla Suk.)、 樟子松( Pinus sylvestris var. mongolica Litv.)、蒙古栎(Quercus mongolica Fisch. ex Ledeb.)、红皮云杉(Picea koraiensis Nakai.)等,生态系统结构丰富[15]。
-
2017年夏季,在大兴安岭山脉构建面积为1 000 m2、半径为17.85 m的样地,确保样地内胸径大于5 cm的兴安落叶松的株数占比大于80%,以此为基础构建的样地认定为纯林样地。利用罗盘仪与手持GPS获取经纬度、海拔数据,样地的经度范围为120.44°~124.51° E,纬度范围为47.53°~53.32° N;海拔范围为283.3~1 087.7 m,坐标分布见图1。对样地内所有胸径≥5 cm的兴安落叶松进行每木检尺,并记录其胸径、树高。然后使用内径为5.15 mm的生长锥对目标树进行年轮条采样,由北向南在高度1.3 m处垂直树干钻入,且钻取的深度超过髓心2~3 mm。将取出的年轮条进行编号、干燥、打磨,使用游标卡尺测量最近5年的胸径生长量,结合每木检尺的胸径数据计算2012年样地内单株木的胸径。根据样地实测数据,选取各样地树高前10%的个体作为该样地的优势木;对该部分林木的年轮条进行测算,提取树龄信息,取平均值作为林分年龄,划分林分发育阶段[16],同时满足模型对林分模拟阶段初始年龄的输入要求。
以2016年最新发布的《中华人民共和国林业行业标准中》的“立木生物量模型及碳计量参数——落叶松”为参考[17],选择适用于大兴安岭的模型,获取落叶松单株的地上、地下及分项生物量数据,以样地为单位,分别累加得到该样地各组分及总生物量。结合样地实测的胸径与树高数据,选用二元模型进行相应指标的计算。计算公式如下:
式中:
$ {M}_{A} $ ——地上生物量(kg)。D ——林木胸径(cm)。
H ——林木树高(m)。
M1、M2、M3、M4——分别为树干、树皮、树枝、树叶的生物量估计值(kg)。
g1、g2、g3——分别为树皮、树枝、树叶生物量与树干生物量的比例。
bi0, bi1, bi2——分别为二元模型的参数,其中i = 1,2,3,分别代表树皮、树枝、树叶。
MB——地下生物量估计值(kg)。
参考该区域落叶松的生物量转换因子函数对应的二元立木材积模型,得到样地内全部单木蓄积,进而累加得到各样地2012年、2017年两期的林分蓄积。
V ——林木蓄积,单位为10−3m³。
本研究所需的气象数据来源于ClimateAP软件(2.20版本)[18],从PRISM[19]和WorldClim[20]中提取月度气候为参考的栅格化数据,并将其匹配至经纬度点位置。通过软件获取本研究所需的各项气象数据,包括月平均最高气温、月平均最低气温、月平均降水量、月平均无霜天数等。通过单位转换与计算,获得用于3-PG模型运行的以下标准指标的气象格式数据。月平均气温范围为−37.8℃~26.3℃,降水主要集中在7—8月,约为111~155 mm。按照林龄将所有样地划分为3个组(见表1)。
林龄分组
Age group林龄范围
Age range/a平均胸径
DBH/cm平均树高
Height/m林分密度
Density/(trees·hm−2)总生物量
Total mass/(t·hm−2)林分蓄积
Volume/(m3·hm−2)样地个数
Plot number第一龄组 Group 1 29~46 8.87~31.18 10.18~33.90 200~670 12.93~196.95 13.69~314.56 11 第二龄组 Group 2 47~64 14.48~28.67 15.79~27.24 290~740 37.92~185.30 49.66~267.20 11 第三龄组 Group 3 67~91 12.78~29.25 16.07~28.30 120~770 27.38~159.69 36.08~215.08 12 Table 1. Forest type division and basic data
-
3-PG模型的工作原理可以简单概括为:基于林木生理参数构建一系列方程,模拟林分在环境或人为因素的影响下的动态生长过程。模型从结构上可划分为3个模块:碳吸收与固定、生物量分配与损耗以及土壤水分动态平衡。碳吸收与固定模块模拟太阳光辐射的逐级递减,被冠层吸收从而进一步转化为总初级生产力与净初级生产力,模拟碳固定的动态过程[21]。生物量分配与损耗模块基于树种的生理特性和林分密度,计算净初级生产力在林木各部分的分配,并考虑林分自疏现象。土壤水分动态平衡模块模拟降水、冠层截流、地表蒸散以及灌溉措施对土壤水分含量的动态影响。3-PG模型保留了传统经验模型与过程模型的优点,通过输入气象数据与简单的立地初始数据,即可输出模拟结果,在经过本地化的参数调整后,即可预测林分的胸径、树高、蓄积、各组分生物量等,且精度较为可靠。
-
3-PG模型运行需要的参数主要分为以下两部分:样地初始信息数据,包括林分纬度、样地调查初期的叶生物量、枝干生物量与根生物量等;模型方程组参数,包含模型的核心子模块:生物量异速生长方程、冠层吸收效率、水分利用效率、自疏与死亡相关参数等。具体参数的详细介绍可见文献[22]。所有参数的获取方式包括实测数据直接计算,方程拟合获取、查阅文献获得、类似树种参数推导以及采用模型默认的初始值。表2列出了模型运行所需的关键参数。
参数
Parameter值
Value来源
Source胸径2 cm时树叶与树干分配比 Foliage-stem partitioning ratio with DBH = 2 cm 0.8 [23] 胸径20 cm时树叶与树干分配比 Foliage-stem partitioning ratio with DBH = 20 cm 0.0513 F 枝干生物量与胸径关系的常数值 Constant in the relation of stem mass and DBH 0.083 F 枝干生物量与胸径关系的幂值 Power in the relation of stem mass and DBH 2.448 F 生长最低生长气温 Minimum temperature for growth/℃ −5 [23] 生长最适气温 Optimum temperature for growth/℃ 16 [23] 生长最高气温 Maximum temperature for growth/℃ 35 [23] fq = 0.5时水分亏缺比 Moisture ratio deficit for fq = 0.5 0.5 F 水分亏缺比的幂值 Power of moisture ratio deficit 5 F 林龄为0时的比叶面积 Specific leaf area at stand age was 0 12.93 [24] 成熟叶的比叶面积 Specific leaf area for mature leaves/(m2·kg−1) 5 [24] 比叶面积为(SLA0 + SLA1)/2时的林龄 Stand age at which specific leaf area was (SLA0 + SLA1)/2/a 8 [24] 冠层量子效率 Canopy quantum efficiency/(mol·mol−1) 0.035 [5] 定义气孔对饱和水汽压差的响应 Defines stomatal response to VPD/(1·mBar−1) 0.05 [24] 树干材积关系中常数值 Constant in the stem volume relationship 0.001125056 F 树干材积关系中胸径的幂值 Power of DBH in the stem volume relationship 1.9675 F 树干材积关系中树高的幂值 Power of height in the stem volume relationship 1.1606 F 注:F为拟合及样地实测参数。
Notes: F means fitting and observed parameters.Table 2. Modification of 3-PG model parameters
表中涉及的主要公式有:
式中:
$ {P}_{FS} $ ,$ {p}_{2},{p}_{20} $ 为生物量分配到叶与枝干的比值,$ {p}_{2},{p}_{20} $ 分别代表胸径2 cm与20 cm时的取值;$ {W}_{S} $ 为枝干生物量,$ {a}_{s},{n}_{s} $ 分别为枝干生物量与胸径关系的常数值与幂值;$ V $ 为林分蓄积,$ {a}_{V} $ 为树干材积关系常数值,$ {N}_{VB} $ 为胸径幂值,$ H $ 为树高,$ {N}_{VH} $ 为树高幂值。模型基于参数集与样地初始信息数据运行,可生成样地各指标的模拟数据(
$ \widehat{y} $ ),将其与样地实测调查的真实数据(y)分别进行简单线性回归,分析模拟值与实测值的匹配程度,使用软件为R(版本4.04)。计算指标为斜率(slope)、相关性系数(R2)、平均误差(ME)、平均绝对误差(MAE)、均方误差(RMSE)、平均相对误差(MRE)。其中相关性系数和斜率用于衡量模拟值与实测值的趋势,其余4个指标用于衡量模拟值与实测值的偏差程度。 -
将34个样地按照“五折法”分为5份,抽取80%的数据(27个样地)作为校正样本,将模型输出值与实测值做线性回归,并按照1.3.3计算各统计指标,结果如表3。
指标
Indicators胸径
DBH树高
HHeight枝干生物量
Stem mass叶生物量
Foliage mass根生物量
Root mass地上生物量
Above ground mass整株生物量
Total mass林分蓄积
Volume斜率 Slope 0.94 1.14 0.95 0.69 1.15 1.03 1.07 0.85 决定系数 R2 0.98 0.81 0.98 0.68 0.96 0.98 0.98 0.91 平均误差 ME −0.67 −2.23 −1.26 −0.07 3.29 4.96 8.27 −25.59 平均绝对误差 MAE 0.74 3.93 4.22 0.66 3.88 5.36 9.10 27.56 均方误差 RMSE 1.08 5.02 4.48 0.76 4.38 6.84 10.76 33.77 平均相对误差 MRE/% −3.33 −8.87 −2.02 −2.35 13.19 6.92 8.54 −19.77 Table 3. Results of model calibration
通过观察结果可以得知,模型各变量的决定系数范围为0.68~0.98,其中胸径、枝干生物量、地上生物量、整株生物量的相关性系数R2均为0.98;误差指标中,胸径的平均误差为−0.67 cm,平均绝对误差为0.74 cm,均方误差1.08 cm;胸径、枝干生物量、叶生物量的平均相对误差较低,控制在−3.33%~−2.02%,树高、地上生物量、总生物量的平均相对误差也控制在10%以内。综上所述,模型的校正结果较为理想。
-
按照“五折法”将建模校正样本外不重叠的20%剩余样本作为验证数据,带入建模校正后的模型获得输出结果,将模型输出值与实测值做线性回归与统计检验,结果见表4。
指标
Indicators胸径
DBH树高
Height枝干生物量
Stem mass叶生物量
Foliage mass根生物量
Root mass地上生物量
Above ground mass整株生物量
Total mass林分蓄积
Volume斜率 Slope 0.81 1.53 0.87 0.60 1.03 0.78 0.88 0.75 决定系数 R2 0.97 0.77 0.98 0.80 0.93 0.97 0.98 0.73 平均误差 ME −1.57 2.76 −2.29 −0.49 4.37 3.05 7.42 −36.48 平均绝对误差 MAE 2.20 4.68 5.78 0.88 5.21 7.28 9.52 57.39 均方误差 RMSE 2.45 5.01 6.94 1.09 6.52 9.98 11.12 67.47 平均相对误差 MRE/% −7.54 −8.80 −3.08 −15.49 11.84 3.69 6.21 −19.80 Table 4. Results of model validation
通过验证结果可以得知,枝干生物量与总生物的拟合效果依然保持在最佳状态,R2达到0.98,而枝干生物量的平均相对误差为−3.08%,是所有变量中的最小值,胸径、树高、地上生物量与总生物量的平均相对误差也控制在±10%以内。叶生物量在验证结果中表现有所提升,相关性系数R2由校正结果中的0.68上升到0.80,同时叶生物量的平均误差、平均绝对误差和均方误差均为最低值。因此可以认为模型的验证效果较为良好。
-
将所有34块样地的林木数据、气象数据与样地初始数据输入到模型中,得到输出结果并进行分析。结果表明,3-PG模型可以较好地估算林分的胸径、树高、生物量与蓄积生长情况(表5)。
指标
Indicators胸径
DBH树高
Height枝干生物量
Stem mass叶生物量
Foliage mass根生物量
Root mass地上生物量
Above ground mass整株生物量
Total mass林分蓄积
Volume斜率 Slope 0.95 1.28 0.93 0.72 1.21 0.98 1.04 1.05 决定系数 R2 0.92 0.83 0.98 0.78 0.96 0.97 0.98 0.77 平均误差 ME −1.94 0.35 −1.47 −0.21 4.73 3.33 8.06 −14.25 平均绝对误差 MAE 1.99 2.08 3.46 0.59 5.43 5.22 8.89 33.31 均方误差 RMSE 2.58 2.55 4.35 0.71 6.12 6.53 10.65 39.09 平均相对误差 MRE% −9.57 −1.63 −2.22 −7.25 7.24 4.50 7.95 −10.13 Table 5. Comparison of fitting results of 3-PG model with measured values
从表5中可以看出,将全部样地数据应用于模型中之后,枝干生物量与整株生物量的模拟值与实测值相关性系数依然保持所有变量的最高值,为0.98;树高、根生物量与蓄积的相关性系数较验证数据有所提升,其中树高相关性系数增长幅度最大,由0.77上升到0.83,同时平均相对误差减少的最大幅度也出现在树高部分,由-8.8%缩小为-1.63%。除蓄积外,所有指标的平均相对误差都控制在±10%以内,可以认为模型模拟效果良好。
-
林龄对生物量分配有显著影响。模型对不同林龄样地整株生物量的拟合结果见表6。
林龄
Forest age截距
Intercept斜率
Slope平均误差
ME平均绝对误差
MAE均方误差
RMSE平均相对误差
MRE/%样地个数
Plot number第一龄组 Group 1 2.16 0.94 4.58 6.91 8.31 4.42 11 第二龄组 Group 2 −1.43 0.95 7.53 7.76 9.44 6.56 11 第三龄组 Group 3 −2.28 0.90 11.75 11.75 13.28 13.43 12 Table 6. Simulation of total biomass in different age sample plots
对不同林龄样地进行的整株生物量拟合效果检验表明,模型对3个龄组林分的拟合效果都取得了较为理想的效果,p值均为极显著。第一龄组与第二龄组的平均相对误差控制在了6.56%以内,其中第一龄组的拟合效果最为理想,平均误差为4.58 t·hm−2,平均绝对误差为6.91 t·hm−2,均方误差为8.31 t·hm−2,平均相对误差更是控制在4.42%。随着年龄的增长,4个评价指标都出现小幅增长的趋势,其中第三龄组的平均误差为11.75 t·hm−2,平均绝对误差为11.75 t·hm−2,均方误差为13.28 t·hm−2,平均相对误差为13.43%。前两个龄组的平均相对误差控制在4.42%~6.56%之间,效果良好。模拟值与实测值的线性回归结果见图2。
-
在天然林中,由于受到环境复杂性、物种多样性、个体生长阶段差异性的影响,林分胸径在不同林龄时变化也有较大差别。因此本研究也针对模型在不同林龄的样地中胸径拟合的表现进行了分析验证。模型对不同林龄样地胸径的拟合结果如表7。
林龄
Forest age截距
Intercept斜率
Slope平均误差
ME平均绝对误差
MAE均方误差
RMSE平均相对误差
MRE/%样地个数
Plot number第一龄组 Group 1 1.20 1.04 −2.05 2.05 2.29 −9.82 11 第二龄组 Group 2 2.31 0.97 −1.75 1.75 1.96 −8.08 11 第三龄组 Group 3 0.73 1.04 −1.51 1.65 1.98 −7.02 12 Table 7. Simulation of DBH in different age sample plots
对不同林龄的样地进行的胸径拟合效果检验表明,模型对3个龄组林分的拟合效果都较为理想,p值均为极显著,且平均相对误差控制在±10%以内,模型拟合效果良好。与整株生物量拟合的结果不同,胸径的拟合效果随林龄的增长没有下降,反而有极小幅度的提升,平均误差由第一龄组的−2.05 cm下降到第三龄组的−1.51 cm,平均绝对误差由第一龄组的2.05 cm下降到第三龄组的1.65 cm。模拟值与实测值的线性回归结果见图3。
Application of 3-PG Model in the Prediction of Growth Factors in Natural Larix gmelinii Forest
- Received Date: 2021-08-30
- Accepted Date: 2021-10-29
- Available Online: 2022-02-20
Abstract: