广东省常绿阔叶林生物量生长模型

曹 磊，刘晓彤，李海奎，雷渊才

(中国林业科学研究院资源信息研究所，北京 100091)

在全球气候变化的大背景下，森林作为陆地上最大的生态系统，其吸收、储存二氧化碳的机制和过程被广泛认为是最有效的固碳方式，其碳汇作用受到越来越多国家和地区的重视[1-3]。迄今为止，国内外已建立了各种形式的生物量模型[4-6]。已建立的生物量模型大多采用基于幂函数形式的异速生长方程[7-9]。异速模型中作为解释变量的胸径和树高等测树因子，其本身就是生长的结果，所建模型虽精度较高，但模型的建立与建模所选样本的范围和区域有关，缺乏一定的生物学意义[10-11]，且在实际应用中，胸径和树高等测树因子往往需要样木的野外复测，再代入到相应的异速模型中，费时费力。

近年来，越来越多的学者将时间变量考虑到方程中，构建了含有年龄变量的异速生长方程。如薛春泉等[12]基于樟树、木荷和枫香3 个广东省主要乡土阔叶树种的生物量实测数据，构建了含有年龄和胸径的单木生物量异速模型，同时还考虑了单木干材、干皮、树枝和树叶与地上部分总量相兼容的问题。随后，薛春泉等[13]又基于Schumacher，Chapman-Richards, Logistic 和Korf 等4 种理论生长方程构建了不同起源（天然林和人工林）的地上与地下部分单木生物量模型，并通过联立方程组控制总量的方法解决了地上部分各组分生物量的相容性问题。考虑时间因素的异速模型生物量始终保持幂指数生长的趋势，生物量随年龄变化生长的生物学意义难以体现，且模型仅局限在林龄已知的人工林和碳汇造林的生物量估计，模型的应用和推广有一定的限制。具有生物学意义的理论生长方程在森林生长与收获中应用较广[14-16]。理论生长方程多用来描述胸径和树高等测树因子的生长趋势[17-19]，很少用来描述生物量与时间的理论生长关系。由于观测数据获取的困难性和立地条件的复杂性等，使得关于生物量生长模型的研究鲜有报道，尤其是区域尺度生物量生长模型的相关研究。

广东省的森林类型主要以常绿阔叶林为主，主要包括栎类（Quercusspp.） 和木荷（Schima superbaGardn. et Champ.）等阔叶树种，蓄积量大，分布范围广，是广东省碳汇造林的主要树种。本研究基于广东省1997 至2017 年的森林资源连续清查数据，筛选有五期连续生长30 株以上保留木的固定样地30 个（共计1 412 株样木），构建了基础生物量生长模型、含有林分特征的生物量生长模型、含有林分特征和立地条件的生长模型等不同形式的林分生物量生长模型，以期为广东省区域尺度森林生物量和碳储量的估算提供模型和方法，提高区域尺度森林生物量和碳储量估算的准确性，也为生物量生长模型在其他区域的应用提供参考。

1 研究区概括和数据

1.1 研究区概括

30 个固定样地的样木数据来源于广东省森林资源连续清查数据，固定样地面积约为0.067 hm2，植被类型属于常绿阔叶林，包括栎类和木荷等常绿阔叶树种，样地林分起源全部为天然林，地类为乔木林地，样地海拔高度在300 ～ 1 300 m 之间，地形主要以低山丘陵为主，土壤质地以壤土和砂壤土为主。30 个固定样地在广东省的分布如图1 所示。

图1 30 个样地在广东省的分布Fig. 1 Distribution of 30 fixed plots in Guangdong province

1.2 数据

表1 给出了30 个固定样地五期相应的林分统计量， 包括林分平均胸径Dg（ 定义为均方直径）、林分胸高断面积（Ba）以及林分密度指数（SDI）[20]等。样木类型包括进界木和保留木等。表2 中列出了30 个样地建模样本的株数和林分平均胸径（均方直径）的统计量。

表1 30 个固定样地不同年份的样木数据统计量Table 1 Descriptive statistics of sample trees in 30 permanent plots at different years

表2 建模样本数量和林分平均胸径统计量Table 2 Descriptive statistics of number and stand diameter at breast height of sample treesfor model data

2 方法

2.1 样地平均年龄的确定

同一时间序列下，对于保留木来说，只需确定第一期样地的林分平均年龄，再加上相应时间序列的时间间隔，即可推算出每期每块样地的林分平均年龄。根据相应时间序列首期和末期样地平均胸径的平均生长速率v，每块样地第一期（1997）林分平均年龄t 按以下公式确定：

2.2 样地生物量的计算

我国已发布的13 个常见树种立木生物量模型和碳计量参数的行业标准提供了不同区域单株立木生物量的异速模型，包括一元模型和二元模型以及地上、地下生物量模型[21-23]，其他文献关于广东省区域内相应树种的立木生物量异速模型也均有报道[24]。

2.3 基础生物量生长模型

本研究选择Richards 理论生长方程作为基础方程[25]，模型结构如下：

式中，Y 为样地生物量，t 为样地年龄，a, b, c 为模型参数。

2.4 含有林分特征的生物量生长模型

为解释不同林分密度状况对林分生物量生长的影响，刻画不同林分结构条件下生物量生长的差别，基于Richards 基础方程，将方程中相关的参数描述为与林分平均胸高断面积和林分密度指数相关的参数方程，将式（3）扩展成含有林分特征的生物量生长模型，模型形式如下：

式中，Ba 为林分胸高断面积（m2·hm－2），SDI 为林分密度指数，k 为常数，其他参数意义同式（3）。

2.5 含有林分特征和立地条件的生物量生长模型

式（3）中参数a 和c 是生长曲线的最大值和相应的形状，分别对参数a 和c 进行有序样本聚类，将30 块样地依次划分为3 种类别，以此将区域内不同林分立地状况分为3 类，代表3 种不同立地等级和生产力（优、良和差），依次求得每个样地参数a 和c 的所属类别，将相应的哑变量加入到式（4）中，得到含有林分特征和立地条件的分级生物量生长模型，其形式如下[26-27]：

式中， S11, S12是 参数a 的哑变量， S21, S22是参数c的哑变量。a1, a2, a3, c1, c2, c3分别是相应的参数值。

2.6 模型检验与评价

采用决定系数（R2）、调整决定系数（Ra2），估计值的标准误（SEE），总相对误差（TRE）、平均系统误差（ASE）、平均百分标准误差（MPSE）、赤池信息准则AIC 和BIC 评价生物量生长模型的拟合效果，公式详见文献[28]。关于样本数据是否要划分为建模样本和检验样本的问题，仍具有争议[28-31]。本研究构建的林分生物量模型是以30 个样地保留木数据为基础，为了充分利用30 个样地林分生物量信息，不区分建模样本和检验样本[28,32-33]，将30 个样地生物量的建模样本作为检验样本对所建的模型进行检验和评价。

3 结果与分析

表3 给出了3 种不同形式生物量生长模型拟合结果的统计指标。基础生物量生长模型的调整决定系数Ra2不到0.48，平均相对误差（MPSE）高达37.40%，AIC 和BIC 值也均达到了1 000 以上，拟合效果最差。将相应林分特征加入基础生物量生长模型后，决定系数R2提高到了0.836，估计值的标准差（ SEE） 下降到了36.92 kg， 总相对误差（TRE）、平均系统误差（ASE）、平均百分标准误差（MPSE）和平均预估误差（MPE）分别下降了 0.93%， 4.17%， 16.36% 和 3.25%， AIC 和BIC 也均明显降低，拟合效果得到一定程度改善。含有林分特征的生长模型中通过哑变量反映不同的立地条件后，拟合效果得到进一步提高，决定系数R2达到了0.961，TRE、均ASE 接近于0，SEE、MPSE、MPE、AIC 以及BIC 值均达到最小。

图2 给出了3 种不同形式生物量生长模型的残差图，可以看到式（3）模型拟合结果的残差（图2（A））表现出明显的 “喇叭状” 规律，随着样地平均年龄的增大，生物量估计值的变动情况变大，呈现出明显的异方差性，拟合效果最差。式（4）模型拟合结果残差图（图2（B））“喇叭状”虽没有式（3）明显，但也呈现出一定程度的异方差性。式（5）模型拟合的残差图（图2（C））在x 轴两侧的分布较为均匀，残差随样地年龄的增加变动幅度较小，无明显的异方差现象，拟合效果最优。

表3 3 种不同形式生物量生长模型拟合效果的评价统计指标 Table 3 Statistics indices of goodness-of-fit for three different forms of biomass growth function

图3 给出了不同生长模型实际值与观测值之间的关系。式（5）模型拟合的估算值与实际值的坐标点均匀分布在坐标轴对角线两侧（图3（C）），模型的估算值与实际值相差较小，模型拟合效果最优。而式（3）与式（4）模型拟合的相应坐标点零散分布在坐标轴对角两侧（图3（A）、（B）），模型估算值与实际值相差较大，尤其是式（3）模型，明显高估了生物量较小的样地，而严重低估了生物量较大的样地，拟合效果最差。

基础生长模型仅仅考虑了林木自身的生长特性，仅以此建立的生物量生长方程不能很好地解释林分生物量的累积变化，拟合效果最差，且出现严重的异方差现象；林分胸高断面积和林分密度指数是广泛使用的林分密度指标，将相应林分密度指标体现在模型参数上，很好地解释了林分内林木的拥挤状况等林分特征对林木生长的影响，使得模型拟合效果得到了很大提高，一定程度上改善了模型的异方差现象；另外，由于林分所处地理环境的差异，林分本身供给林木生长的养分不同，即不同的立地条件下，林分具有不同的立地生产力。基于哑变量方法，通过划分不同立地等级，在模型参数上体现了不同立地的生产力，解释了由于立地生产力的不同而导致林分生长变化的影响，大大提高了生长模型的拟合精度（表4）。

4 讨论

单木生物量的估算一般基于异速生长方程，对样地所有单木进行汇总或者取标准木的方法，得到林分水平上的生物量[34]。而区域尺度上生物量的估算一般基于生物量与材积的比值，包括IPCC 法，转换因子连续函数法和加权生物量回归模型法[35-36]，但建模过程复杂，需要大规模重复的野外测量数据。本研究基于Richards 理论生长方程，通过构建含有林分年龄的生物量生长模型，从而实现了对特定时间下区域尺度上的生物量预估，可用于区域尺度生物量和碳储量的估算和碳汇潜力的评估。

图2 3 种不同形式生物量生长模型的残差分布（A）、（B）和（C）表示式（3）、式（4）和式（5）Fig. 2 Residuals of three different forms of biomass growth functions (A), (B), (C) represents Eq. (3), Eq. (4) and Eq. (5)

图3 3 种不同形式生物量生长模型实际值与估算值关系（A）、（B）和（C）表示式（3）、式（4）和式（5）Fig. 3 Relationofactual value and prediction for three biomass growth models. (A), (B), (C) represents Eq. (3), Eq. (4) and Eq. (5)

生物量理论生长方程通过建立区域尺度生物量与时间（年龄）之间的关系，反映了区域尺度下生物量生长的变化趋势，可用于估算出某个年份下区域的森林生物量。区域尺度上林分基准年龄的确定是关键。龙时胜等[36]通过设置首期不同年龄采集样地内各树种的直径与年龄坐标信息，并基于树高生长至1.3 m 处的年龄，采用树高生长方程准确定位，从而估算异龄林林分年龄。随后，龙时胜等[37]又根据相同直径的优势木和被压木的关系，通过年龄变化范围内的随机模拟，估计初始阶段的立木年龄推算林分年龄。本研究确定初始阶段林分年龄与龙时胜等的方法类似。

表4 最优生物量生长模型参数估计值的统计量Table 4 Statistics indices of parameter estimate for the optimal biomass growth model

Richards 基础生长模型不能很好的解释生物量生长过程，需要考虑林分特征和立地等因素解释其生物量的生长变化。林分密度是影响林分生长的重要因素之一，也是营林工作中能够有效控制的因子。慧淑荣等[38]在Richards 生长方程的基础上，将立地指数和林分密度引入到方程中，构建了日本落叶松林分胸径、树高、断面积和蓄积的理论生长方程，精度检验和估测效果均较好；Sharma 等[19]分龄级和森林类型建立了包含林分密度，林分胸高断面积以及地位指数等林分密度指标在内的北美短叶松和黑云杉胸径-树高Richards 理论生长方程。

混合效应和哑变量模型兼顾了不同区域的现实条件，是考虑到了不同情况相互兼容的两种建模方法。但有研究表明，利用哑变量模型的拟合效果优于混合效应模型[39]。本研究针对不同的立地条件，通过构建哑变量反映不同样地类别的划分，兼顾到了不同样地之间的差异。这种差异可能来源于不同的立地因子如坡度、坡向和海拔等，立地因子反映了立地生产力的不同，在模型上主要体现在相应哑变量参数值的不同。理论生长方程中参数a 是方程的最大值，一般认为与立地条件和生产力有着密切关系。划分的3 种不同立地代表了30 个样地的3 种不同立地状况，同时也体现了3 种不同类别立地潜在生产力。划分不同立地等级，同时考虑林分特征变量，兼顾了不同立地的潜在生产力，拟合精度较高，预估效果较好，可满足区域尺度森林生物量和碳储量的估算要求，另一方面，林分密度变量可以为森林经营管理者提供控制林分密度的参考，为合理评价区域碳汇潜力和提高固碳速率提供一定借鉴。

应用在区域尺度下林分生物量模型，样本单元数至关重要，大尺度范围应用的立木生物量模型，样本单元数的增加可以提高模型的准确度和通用性[6]。本研究仅选取了广东省森林资源清查30 个固定样地数据进行建模，在样本单元数上仅仅满足了区域大样本的要求[6]，为得到更准确的生物量模型，需要覆盖全省范围尽可能多的固定样地数据。

5 结论

本研究基于五期（1997, 2002, 2007, 2012, 2017）广东省森林资源连续清查资料，通过30 个固定样地五期保留木数据，采用Richards 理论生长方程，分别构建和比较了基础生长方程、含有林分特征的生长方程和含有林分特征和立地条件的生长方程等不同形式的生物量生长模型。在林分生物量生长模型中，引入林分胸高断面积和林分密度指数等指标以及反映不同立地生产力差异的哑变量，极大地改善了模型的拟合效果，提升了模型预测能力，为广东省区域尺度常绿阔叶林生物量和碳储量估算提供了借鉴和参考。