-
生物量一般通过材积信息和生物量扩展因子进行估计,也可利用生物量模型进行估计[1]。从最新研究结果知,对生物量的估计,采用生物量模型法准确性最高,而林木水平的生物量模型其预精度要高于林分水平的模型[2]。根据国内外文献,目前所建林木水平的生物量方程已经涵盖了全球不同生态地理区域的主要树种[3-12]。从建模方法看,已经从原来拟合单个的独立的生物量方程逐渐发展成为构建相容的可加的生物量联立方程组[13-27]。但这些模型通常只是解决总量与分量的可加性问题[14-16, 20, 24-25]或生物量与材积的相容性问题[17, 19, 21],同时解决相容性和可加性问题的文献很少[23, 26]。此外,还有学者利用哑变量模型或混合模型方法,建立不同区域、不同起源的立木生物量模型[17, 23, 28-30]。
云杉(Picea spp.)是我国重要的针叶树种,广泛分布于东北、西北和西南诸省,其面积蓄积仅次于落叶松[31]。在新疆的森林蓄积中,落叶松林蓄积只占近1/4,而云杉林蓄积占了一半以上。本研究将利用采集自新疆自治区的云杉立木生物量样本数据,研建一体化立木生物量模型系统,主要研究目标包括:(1)利用误差变量联立方程组方法,同时解决地上生物量、地下生物量、立木材积与生物量转换因子和根茎比等因子之间的相容性问题,以及地上生物量与各分量之间的可加性问题;(2)利用哑变量模型方法,解决在地上生物量和地下生物量样本单元数不相等的情况下,如何建立一体化联立模型系统的问题,并分析按不同区域建立一元立木生物量模型是否能显著提高预估精度。
HTML
-
本研究所用数据来自国家森林资源连续清查生物量调查建模项目组,包括230株样木的地上生物量和干、皮、枝、叶各个分量及立木材积实测数据,其中78株样木还实测了地下生物量。按区域分,150株样木采集自阿尔泰山林区,80株样木采集自天山林区。采集的样本在新疆云杉分布区域内具有广泛的代表性。
全部样木株数按10个径阶(2、4、6、8、12、16、20、26、32 cm以及38 cm以上)均匀配置,每个径阶约23株样木,各个径阶的样木也要求尽可能均匀配置在不同的树高级。每一株样木伐前测量胸径和冠幅,伐后测量树高、冠长和年龄等因子。然后分上、中、下3段测量树干总鲜质量和3个圆盘鲜质量,并将树干从0、0.05、0.1、0.2、0.3、0.4、0.5、0.6、0.7、0.8和0.9相对树高处分成11段,根据测量的每段基部直径用Smalian公式计算树干材积,即每段材积等于大头和小头断面积平均值乘以每段长度。再分上、中、下3层测量树枝总鲜质量,并从各层挑选样枝,根据针叶和无叶枝条的鲜质量比,推算整个样木树枝和树叶鲜质量,同时选取树枝和树叶样品称其鲜质量。从全部样木中选取约1/3(共78株)同时测量地上和地下生物量。根据样木大小,采用机械或人工方法将整个根系挖出,将根上的泥土清除干净后,分别树桩(根蔸)、粗根(大于10 mm)和细根(210 mm,不包括2 mm以下须根)称其总鲜质量,并选取样品称质量。所有样品带回实验室,在85℃温度下烘干至恒质量。根据干质量和鲜质量之比,算出各个部分的生物量及样木的总生物量。样木的选取,干材、树皮、树枝、树叶和树根等各个部分的测量、取样和干质量测定,都严格执行林业行业标准[30]。表 1给出了新疆云杉分区域的立木生物量建模样本的统计特征值。
区域
Region项目
Item样本量
Sample size变量
Variable平均值
mean最小值
min.最大值
max.标准差
SD变动系数
CV/%阿尔泰山林区
Altai mountainousregion地上生物量
Above-ground biomass150 胸径 Diameter/cm 17.3 1.3 65.5 13.6 78.9 树高 Height /m 12.8 1.5 36.8 8.4 67.5 立木材积 Volume/(dm3·株-1) 356.5 0.6 4 163.4 611.1 171.4 地上生物量 Above-ground biomass /(kg·株-1) 176.9 0.4 1 540.8 260.5 147.3 生物量转换因子 Biomass conversion factor 0.686 9 0.164 6 1.585 9 0.265 1 38.6 地下生物量
Below-ground biomass50 胸径 Diameter/cm 17.1 1.3 52.0 12.8 75.0 树高 Height/m 12.9 1.5 33.6 8.5 65.8 地上生物量 Above-ground biomass/(kg·株-1) 168.5 1.1 832.0 225.2 133.6 地下生物量 Below-ground biomass/(kg·株-1) 37.3 0.3 260.9 56.4 151.2 根茎比 Root-to-shoot ratio 0.295 8 0.139 7 0.605 8 0.128 3 43.4 天山林区
Tianshan mountainous region地上生物量
Above-ground biomass80 胸径 Diameter /cm 17.2 1.0 57.0 13.1 76.0 树高 Height/m 13.5 1.5 44.0 9.0 66.5 立木材积 Volume/(dm3·株-1) 368.5 0.6 5 248.8 725.6 196.9 地上生物量 Above-ground biomass/(kg·株-1) 185.8 0.7 1 090.2 254.6 137.0 生物量转换因子 Biomass conversion factor 0.656 5 0.324 1 1.160 4 0.202 6 30.9 地下生物量
Below-ground biomass28 胸径 Diameter/cm 16.2 2.0 38.9 11.4 70.7 树高 Height/m 12.2 1.9 35.0 8.2 76.6 地上生物量 Above-ground biomass/(kg·株-1) 162.6 1.6 726.9 208.2 128.1 地下生物量 Below-ground biomass/(kg·株-1) 54.2 0.2 269.1 76.4 141.0 根茎比 Root-to-shoot ratio 0.259 1 0.052 5 0.522 4 0.109 2 42.2 Table 1. Statistics of individual tree biomass samples of Picea spp. in Xinjiang
-
式中,y为生物量(kg)或树干材积(dm3);xj为反映林木大小的变量,如胸径(D,cm)和树高(H,m);βj为待估参数;ε为误差项,假定其服从正态分布。因为生物量和材积数据都表现为异方差性,在参数估计时应当采取消除异方差的措施。本研究采用加权回归估计方法,其权函数根据最小二乘回归模型的残差拟合而来[16, 34-35]。因为基于胸径D的一元模型和基于胸径D与树高H的二元模型都得到广泛应用,本文将建立一元和二元两套模型。地上生物量、地下生物量和树干材积方程的二元模型分别表述如下:
式中,Ma和Mb分别为地上生物量和地下生物量(kg),V为树干材积(dm3),ai、bi和ci为待估参数,其它符号同前。
-
在生物量建模实践中,相容性一般是指地上生物量与立木材积方程及生物量转换因子模型之间的参数相关性和协调一致性[17, 19, 21],也包括地下生物量与地上生物量方程及根茎比模型之间的参数相关性和协调一致性[36-38]。因此,为解决相容性问题,可通过生物量转换因子(BCF)和根茎比(RSR)将地上生物量、立木材积及地下生物量三者关联起来。可加性是指地上生物量或总生物量与各个分项生物量模型之间的协调一致性,以确保各分量估计值之和等于总量[16, 25-26]。另外,由于测定地下生物量与地上生物量的样木数量不相等,前者仅为后者的1/3,通常的做法是分别建立地上生物量和地下生物量联立方程组[33, 39]。为了解决这一问题,此处引入一个哑变量x,对于测定了地下生物量的样木,令x=1;对于未测定地下生物量的样木,令x=0[26, 40]。为了保证地上生物量与干、皮、枝、叶各个分量的可加性,本研究拟按比例控制和代数控制2种方法进行建模对比[20, 25-26],其二元模型的一体化联立方程组分别为:
式中,BCF和RSR分别为生物量转换因子和根茎比,M1、M2、M3、M4分别为干、皮、枝、叶生物量,$\hat{M}_{a}$、$\hat{M}_{b}$、$\hat{V}$、$\hat{M}_{2}$、$\hat{M}_{4}$分别为地上生物量、地下生物量、立木材积、树皮生物量、树叶生物量的估计值,x为标识是否有地下生物量数据的指示变量,fi、gi和hi为待估参数,其它符号同(2)~(4)式。如果删除(5)、(6)式中所有包含H的指数项,二元模型就变成了一元模型,为省篇幅,这里不再列出。
需要补充说明一点,关于(6)式中4个分项生物量方程,把所占比例较小的树皮和树叶生物量设计为独立模型,把比例较大的干材和树枝生物量设计为代数方程之差,是考虑降低传递误差的相对影响,不至于造成代数方程之差出现负值。为了分析传递误差的影响大小,本研究对(6)式中的干材和树枝生物量方程略作调整,变为以下形式:
二者的差异是,在(6)式中f0Df1Hf2表示树冠生物量(包括树枝和树叶),而在(7)式中f0Df1Hf2表示树干生物量(包括干材和树皮)。
估计上述联立方程组参数的方法有多种,其中常用的有似乎不相关回归(SUR)和误差变量(EIV)联立方程组方法。有研究表明EIV方法比SUR方法更具优势和更有潜力[41],本研究采用EIV建模方法,通过ForStat软件来估计(5)、(6)式中的参数,其中Ma、Mb、M1、M2、M3、M4、V、BCF和RSR作为误差变量,D和H作为无误差变量[42]。
-
为了分析不同区域对生物量和树干材积变动的影响,在前述联立方程组的基础上,再引入一个哑变量z。当样木属于阿尔泰山林区时z=1,属于天山林区时z=0。与(7)式对应的含区域因子的一元哑变量模型系统如下:
式中, a1、a3、b1、b3、c1、c3、f1、f3、g1、g3、h1和h3为哑变量z的待估参数,其它符号同前。若其中部分参数在统计上与0无显著差异,可将其从模型中剔除。
-
用于评价立木生物量模型或材积模型的统计指标有很多[43-44],本研究选用以下4项,即:确定系数(R2)、估计值的标准误(SEE)、平均预估误差(MPE)和平均百分标准误差(MPSE)。另外,为分析不含区域因子模型(以下称总体平均模型)与包含区域因子模型(即哑变量模型)之间是否存在显著差异,采用F检验方法[17, 45-46]。利用计算的F值与临界值进行比较,以确定两个模型之间是否存在显著差异。一元和二元模型之间的差异也可采用F统计指标进行检验。
1.1. 数据采集
1.2. 模型建立
1.2.1. 一体化联立方程组
1.2.2. 含区域因子的哑变量模型系统
1.3. 模型评价
-
利用新疆云杉的全部生物量样本,通过ForStat软件中的“非线性误差变量联立方程组”方法,先拟合(5)、(6)、(7)式共3组模型,以下分别称之为比例法、代数法A和代数法B。表 2列出了以上3种建模方法的参数估计值,表 3列出了3组模型各目标变量的统计指标。
参数
Parameter不同建模方法的参数估计值
Parameter estimates of different methods比例法
Proportion method代数法A
Algebraic method A代数法B
Algebraic method Ba0 0.096 26 0.094 69 0.096 33 a1 2.085 86 2.067 71 2.088 19 a2 0.372 32 0.398 62 0.370 08 b0 0.022 23 0.021 69 0.022 21 b1 2.367 41 2.342 59 2.368 90 b2 0.106 51 0.142 74 0.105 45 c0 0.084 31 0.083 35 0.084 66 c1 1.744 14 1.740 86 1.753 66 c2 1.009 76 1.017 23 0.998 74 f0 0.354 66 0.060 09 0.040 91 f1 0.024 36 2.611 94 1.752 23 f2 -0.289 55 -0.420 42 0.883 12 g0 1.345 68 0.011 24 0.010 86 g1 0.999 69 1.753 16 1.763 87 g2 -1.499 88 0.652 32 0.654 19 h0 0.880 81 0.023 61 0.023 55 h1 0.931 59 2.563 80 2.557 86 h2 -1.576 40 -0.463 10 -0.456 05 Table 2. Parameter estimates of two-variables simultaneous equations of different methods
方法
Method指标
Statistics目标变量 Target variable Ma Mb V M1 M2 M3 M4 比例法
Proportion methodR2 0.832 0 0.667 8 0.978 1 0.794 9 0.881 2 0.703 5 0.609 2 SEE 106.17 37.58 96.77 69.26 9.18 40.17 18.46 MPE/% 7.62 19.52 3.47 9.08 7.30 11.14 12.86 MPSE/% 21.20 55.25 14.96 23.71 27.71 37.29 44.55 代数法A
Algebraic method AR2 0.831 6 0.663 4 0.978 5 0.792 6 0.886 2 0.705 4 0.607 6 SEE 106.32 37.83 96.03 69.17 8.80 39.50 18.12 MPE/% 7.63 19.65 3.44 9.07 7.00 10.95 12.63 MPSE/% 21.48 56.13 15.06 35.07 25.26 37.82 45.54 代数法B
Algebraic method BR2 0.830 2 0.667 2 0.977 6 0.798 5 0.883 8 0.706 9 0.611 0 SEE 106.76 37.62 97.87 67.73 8.89 39.66 18.05 MPE/% 7.67 19.53 3.51 8.88 7.07 10.99 12.57 MPSE/% 21.16 55.24 14.91 21.73 25.66 38.87 45.64 Table 3. Statistics of two-variables simultaneous equations of different methods
从表 2可以看出,3种建模方法的前9个参数都是很接近的;2种代数法的后6个参数也很接近,中间的3个参数(f0、f1和f2)差异大,是因为在代数法A中是树冠生物量模型参数,而在代数法B中是树干生物量模型参数。再从表 3中各项统计指标的对比可知,对地上、地下生物量和立木材积的估计,3种建模方法差异很小;对干、皮、枝、叶4个分项生物量的估计,代数法B与比例法不相上下,但要略优于代数法A。也就是说,代数法和比例法这2种满足可加性的方法,其预估效果基本相当,难分优劣。而代数法的2种模型设计,(7)式的效果更好,它可以明显提高干材生物量(M1)的预估精度,而对树枝生物量(M3)预估精度的影响很小。因此,后续的代数法特指模型(7)式。
-
除了与前述比例法和代数法对应的一元联立方程组以外,还进一步拟合了含区域因子的哑变量模型系统。表 4为总体平均模型的参数估计值及根据哑变量模型得出的不同林区的参数估计值,表 5为基于两套一元模型系统算出的各项统计指标。
方法
method参数
Parameter总体平均模型
Population average model哑变量模型
Dummy-variable model新疆自治区
Xinjiang autonomous region阿尔泰山林区
Altai mountainous region天山林区
Tianshan mountainous region比例法
Proportion methoda0 0.166 58 0.168 45 0.173 08 a1 2.239 11 2.221 85 2.250 24 b0 0.039 39 0.050 94 0.027 70 b1 2.288 40 2.152 65 2.492 38 c0 0.166 97 0.166 74 0.172 51 c1 2.434 74 2.419 35 2.453 05 f0 0.286 07 0.303 05 0.245 06 f1 -0.176 29 -0.174 29 -0.163 18 g0 0.677 28 0.662 07 0.730 20 g1 -0.132 17 -0.106 38 -0.192 43 h0 0.469 32 0.588 37 0.321 69 h1 -0.298 54 -0.378 16 -0.164 16 代数法
Algebraic methoda0 0.166 13 0.167 93 0.172 30 a1 2.239 96 2.222 95 2.251 50 b0 0.039 10 0.050 67 0.027 45 b1 2.291 21 2.154 44 2.495 89 c0 0.168 55 0.166 95 0.176 54 c1 2.430 46 2.418 57 2.442 99 f0 0.091 17 0.089 11 0.099 03 f1 2.291 36 2.282 74 2.292 54 g0 0.020 57 0.021 01 0.019 81 g1 2.140 90 2.135 99 2.148 55 h0 0.031 54 0.039 89 0.023 17 h1 2.047 60 1.942 53 2.193 20 Table 4. Parameter estimates of population average and dummy-variable models of different methods
方法
Method模型
Model指标
Statistics目标变量 Target variable Ma Mb V M1 M2 M3 M4 比例法
Proportion method总体平均模型 R2 0.866 7 0.715 6 0.813 0 0.714 3 0.848 6 0.785 8 0.700 4 SEE 94.38 34.54 282.41 81.01 10.27 33.83 16.02 PA
modelMPE/% 6.78 17.94 10.12 10.62 8.17 9.38 11.16 MPSE/% 20.72 44.08 30.46 31.37 27.77 39.55 47.21 哑变量模型 R2 0.877 7 0.820 4 0.829 0 0.718 4 0.848 6 0.781 0 0.735 8 SEE 90.80 27.82 271.24 81.91 10.46 34.84 15.32 DV
modelMPE/% 6.52 14.45 9.72 10.74 8.32 9.66 10.67 MPSE/% 20.83 37.44 30.22 31.24 27.64 39.68 44.93 代数法
Algebraic method总体平均模型 R2 0.866 3 0.714 6 0.814 1 0.717 7 0.848 6 0.786 3 0.696 1 SEE 94.50 34.61 281.61 79.80 10.13 33.65 15.92 PA
modelMPE/% 6.79 17.97 10.09 10.46 8.06 9.33 11.09 MPSE/% 20.71 44.16 30.52 31.50 27.78 39.56 47.07 哑变量模型 R2 0.877 2 0.819 9 0.828 7 0.721 9 0.849 2 0.789 1 0.733 2 SEE 90.98 27.86 271.46 79.92 10.15 33.88 14.98 DV
modelMPE/% 6.53 14.47 9.73 10.48 8.08 9.39 10.44 MPSE/% 20.81 37.47 30.29 31.42 27.64 39.62 44.78 Table 5. Statistics of population average and dummy-variable models of different methods
从表 5与表 3的对比看,对立木材积和干材、树皮生物量的估计,二元模型均明显优于一元模型,F检验结果也显示其差异极为显著;然而对地上、地下生物量及树枝、树叶生物量的估计,二元模型反而不如一元模型,其差异在统计上也极为显著。
从表 5中总体平均模型与哑变量模型的对比看,对地上、地下生物量及立木材积的估计,哑变量模型均要好于总体平均模型,其差异在统计上极为显著;但是对干、皮、枝、叶4个分项生物量而言,除树叶生物量外,对其它3项的估计,哑变量模型与总体平均模型在统计上没有显著差异(见表 6)。
方法 method 目标变量 Target variables Ma Mb V M1 M2 M3 M4 比例法 Proportion method 10.16** 21.57** 10.58** 0.39 0.00 -0.58 3.58* 代数法 Algebraic method 10.01** 21.65** 9.68** 0.83 0.43 0.49 15.74** 注:“**”表示在0.01水平下显著,“*”表示在0.05水平下显著,。 Table 6. F-values of population average and dummy-variable models of different methods
2.1. 二元联立方程组
2.2. 一元联立方程组
-
我国云杉立木生物量建模的研究成果很少[47-48],但云杉是全世界广为分布的树种,为了分析区域之间的差异性,这里选用国内外主要文献中的云杉一元立木地上生物量模型进行对比,表 7列出了不同模型的参数及文献来源,表 8列出了不同模型对各个径阶的地上生物量估计值。
序号
No地上生物量模型
Above-ground biomass models样本量
Sample size应用区域
Region for application文献来源
Reference0 M0=0.166 6D2.239 1 230 新疆 Xinjiang 本研究 1 M1=0.155 6D2.258 8 150 新疆(阿尔泰林区)Xinjiang(Altai) [48] 2 M2=0.105 5D2.389 3 150 云南 Yunnan [48] 3 M3=0.144 4D2.260 4 195 加拿大(新不伦瑞克)Canada(NB) [3] 4 M4=0.125 3D2.332 3 212 美国 USA [4] 5 M5=0.116 0D2.360 0 70 欧洲(北部地区)Europe(Boreal) [7] 6 M6=0.255 0D2.174 0 144 欧洲(温带地区)Europe(Temperate) [7] Table 7. The above-ground biomass models for comparison
D/cm M0 M1 M2 M3 M4 M5 M6 5 6.12 5.90 4.94 5.49 5.35 5.18 8.44 10 28.89 28.23 25.86 26.30 26.93 26.57 38.07 15 71.62 70.55 68.14 65.77 69.33 69.19 91.91 20 136.39 135.12 135.49 126.01 135.63 136.42 171.78 25 224.78 223.67 230.92 208.68 228.23 230.99 279.04 30 338.11 337.64 356.98 315.10 349.17 355.19 414.76 35 477.49 478.27 515.94 446.46 500.24 511.04 579.89 40 643.89 646.65 709.83 603.76 683.02 700.36 775.21 45 838.20 843.74 940.53 787.93 898.96 924.78 1 001.44 50 1 061.22 1 070.44 1 209.76 999.81 1 149.37 1 185.84 1 259.22 合计 Total 3 826.70 3 840.22 4 198.38 3 585.31 4 046.23 4 145.57 4 619.75 相差 RD/% / 0.35 9.71 -6.31 5.74 8.33 20.72 注:M0~M6各列为表 7中各地上生物量模型的估计值;相差行数据是指M1~M6各列合计数与M0列合计数之间的差异百分数。 Table 8. The estimates of above-ground biomass models for different diameter classes
从表 8可以看出,与本研究所建模型相比,6个模型中有5个模型的相对差异在±10%以内,只有最后一个模型的相对差异较大,超过了20%。该模型的2个参数与其它模型的参数相差明显,尤其是幂指数参数(2.17)与通用性模型的幂指数(2.33)[49]相差也很大。从我国、北美和欧洲3组对应模型的数据看,相同直径林木的地上生物量估计值,高纬度地区要小于低纬度地区,即新疆小于云南,加拿大小于美国,欧洲北部地区小于温带地区。其深层次原因有待于作进一步研究。