-
林分密度指标是衡量林木竞争的一个重要指标,确定合适的密度指标是开展林分密度研究的重要前提。许多学者从不同的角度出发提出了各种密度指标来衡量一个林分的密度情况,例如Reinek林分密度指数(SDI)[1]、Nilson密度(SD指数)[2]、优势高-营养面积比(Z指数)[3]、相对植距(RS指数)[4]以及每公顷株数(N)等。因此,选择合适的林分密度指标对于提高生长模型的应用以及模型精度有重要意义。
森林生长模型是实现森林质量精准提升的关键技术。在预估年生长量时,各学者提出了不同方法,包括固定生长率法、内插法[5]、迭代法[6]和可变生长率法[7]。固定生长率法是假设在森林的整个生长期内,林木的年生长率是固定不变的,显而易见,并不符合林木的生长规律,因为随着林分生长,用于预估模型中的林分因子和单木因子时刻都在发生变化,必定会对单木生长量乃至林分蓄积产生影响。Cao等[8]在迭代法的基础上提出可变生长率法,解决了以上存在的问题,且比迭代法更简单,计算用时更少。张雄清等[9]利用可变生长率法建立的单木生长模型考虑了林分因子(林分优势高、林分断面积)和单木因子(树高、胸径),比利用固定生长率法建立的模型误差更小,预测精度更高,更具代表性。之后,Zhang等[10]基于可变生长率法建立了林分断面积、林分平均直径和林分直径标准差模型。
林分蓄积量是评价生产力的一个重要指标,所以选择合适的林分密度指标来构建杉木(Cunninghamia lanceolata (Lamb.) Hook.)林分蓄积量模型就显得尤为重要。本研究基于可变生长率法构建杉木林分蓄积量生长模型,并在模型中引入不同林分密度指标,选出能够应用于杉木林分蓄积量模型中的最优林分密度指标,以期实现杉木林分蓄积量的精准预测。
-
样地设置在福建武夷山北部的邵武市,位于117°43′ E, 27°05′ N。地貌特征主要是高山和低丘陵地区,海拔为250~700 m,坡度为25°~35°,为亚热带季风气候,年平均温度17.7 ℃,1月平均温度6.8 ℃,7月平均温度28 ℃,最低极端气温为−7.9 ℃,年日照时间1 740.7 h,平均霜冻期为95 d。年降水量1 768 mm,年平均相对湿度为82%。气候条件适合杉木生长。
试验林使用1年生苗木于1982年造林,完全随机区组设计,分5种造林密度:A:2 m × 3 m(1 667 株·hm−2),B:2 m × 1.5 m(3 333 株·hm−2),C:2 m × 1 m(5 000 株·hm−2),D:1 m × 1.5 m(6 667 株·hm−2)和E:1 m × 1 m(10 000 株·hm−2)。每个样地大小为20 m × 30 m,每种造林密度均重复3次,总计15个样地。标记了一共4 800棵树,在冬季测量树高和胸径,从1984至1990年,每年进行1次测量;从1992年至2010年,每隔1或者2年进行1次测量。在每个样地中,选取最高的6株树,并计算其平均树高作为林分优势高。具体统计数据见表1。研究数据随机抽取60%用于模型构建,剩余40%用于模型验证。
变量
Attribute最小值
Min.最大值
Max.平均值
Mean标准差
Standard deviation林龄
Age/a4.0 28.0 14.3 7.6 株数
Number /(trees·hm−2)1 175 10 000 4 617 2 557 优势高
Dominant height/m2.9 26.9 12.7 6.1 胸径
Tree diameter/cm3.6 39 9.3 5.1 断面积
Basal area/(m2·hm−2)0.8 81.1 39.3 22.2 蓄积量
Stand volume/(m3·hm−2)0.9 518.0 204.3 144.2 Table 1. Summary statistics of stand and tree variables of Chinese fir plantation
-
利用常用的5种林分密度指标,分别是每公顷株数密度N、林分密度SDI指数、SD指数、优势高—营养面积比Z指数和相对植距RS指数构建杉木林分蓄积量生长模型。
-
(1) SDI指数
SDI指数是Reineke基于完满立木度提出的一种林分密度指标[1]。其优点是不仅能反映林地内的株数,而且能反映林木的大小,同时不受林龄和立地条件的影响,具有容易测算、使用简单的特点[11]。其数学表达式为:
式中:Dq为林分平方平均胸径cm;D0为林分标准直径,杉木一般取值20 cm;β为Reinek自稀疏系数。Sun等[12]认为杉木同龄林中林木的自然死亡率达到2%时才开始启动林分自稀疏,因此在剔除自稀疏死亡率小于2%的样本后,根据Reineke的自稀疏理论,对单位面积株数N和林分平均胸径Dq进行对数转化,即lnN和lnDq计算杉木自然稀疏线斜率β。据此分析,计算得到自然稀疏斜率β=−1.897 2。
(2) SD指数
根据Nilson[2]对林木之间的平均距离L和平方平均胸径Dq的定义,林分相对密度SD指数的公式如下:
式中:a,b为参数,k=−a,D0取值20 cm。
(3)优势高—营养面积比Z指数
张连金等[3]采用单位面积内株数与优势木的平均高来反映林分中林木的相对密度的大小。其表达式如下:
式中: Hd为优势木平均高。该指数与林龄、立地指数等无关,计算相对较为简单。
(4) 相对植距RS指数
相对植距包含了林分优势高和每公顷株数,其表达式如下[11]:
-
本研究选取了优势高-营养面积比Z指数、每公顷株数密度N、Reineke密度指数SDI、Nilson密度SD指数和相对植距RS构建杉木林分蓄积量模型;
式子中:M1,M2分别为第1期和第2期的林分蓄积量;
A1为第1期的林分年龄;
Hd1第1期的林分平均优势高;
K1为第1期时的5种林分密度指标,包括:优势高-营养面积比Z指数,株数N,林分密度指数SDI,林分密度SD指数,相对植距RS指数;
x1~x4,α1~α3为模型参数;
z1~z4为5种不同造林密度设置的哑变量(造林密度1 667 株·hm−2:z1=1, z2=z3=z4=0;造林密度3 333 株·hm−2:z2=1, z1=z3=z4=0;造林密度5 000 株·hm−2:z3=1, z1=z2=z4=0;造林密度6 667 株·hm−2:z4=1, z1=z2=z3=0;造林密度10 000 株·hm−2:z1=z2=z3=z4=0)。
由于调查的间隔并不都是1年1次,因此为了提高模型的预测精度,可变生长率法被引入到本研究中。可变生长率法考虑间隔期内,林分优势高的变化引起的林分蓄积量生长的变化。杉木林分蓄积量年生长模型利用递归的方式推导过程如下:
(t + 1)年时:
(t + 2)年时:
(t + q)年时:
式中,t为调查年,q为调查间隔期。
研究中构建的杉木优势高模型形式如下:
其中β1、β2为待估参数。
本研究中, 杉木林分蓄积量模型的参数估计均利用 SAS中非线性回归模块来完成。利用可变生长率法估计模型参数时, 应用了循环运算。
-
林分优势高模型和林分蓄积生长模型通过统计量平均绝对偏差(MAD)、均方根误差( RMSE)和决定系数(R2) 进行评价。
-
杉木林分优势高模型的参数估计结果以及模型的均方根误差、决定系数见表2。结果表明,模型的参数估计值均为有效值,并且模型的决定系数很高,达到了0.953 9,均方根误差RMSE较小,为1.343 1。因此模型估计得到的林分优势高是可信的,可以用于构建杉木林分蓄积量模型。
参数
Parameter估计值
Estimate标准误差
Standard dev. error决定系数
R2均方根误差
RMSEβ1 1.583 2 0.360 1 0.953 9 1.343 1 β2 1.470 1 0.325 3 Table 2. Parameter estimates and model evaluation of stand dominant height of Chinese fir
分别用5种不同林分密度指标,利用可变生长率法建立了6种杉木林分蓄积量年生长模型,其中包括一种不含密度指标的生长模型。参数估计和标准误差及模型评价见表3。从表3中可以看到,5类含有密度指标的模型决定系数R2均在0.979以上,模型精确度较高,且大于不含密度指标的对照组模型(模型决定系数为0.972 8)。林分密度对于杉木蓄积量具有不可忽视的影响,因此将林分密度指标引入杉木林分蓄积量模型可以提升模型精确度,优化模型表现。
密度指标
Density
index参数
Parameter估计值
Estimate标准误差
Deviation
error决定
系数
R2平均绝对
偏差
MAD均方根
误差
RMSE优势高
营养面积
Zα1 2.591 9 0.313 7 0.979 3 10.452 5 20.651 2 α2 −0.255 7 0.044 3 α3 0.109 1 0.038 2 x1 1.086 1 0.338 1 x2 0.860 4 0.288 0 x3 0.644 9 0.236 0 x4 0.645 0 0.205 6 每公顷
株数
Nα1 −0.066 0 1.316 1 0.979 9 10.268 7 20.363 4 α2 −0.085 8 0.031 0 α3 0.000 3 0.000 1 x1 2.892 5 1.001 7 x2 2.281 8 0.808 9 x3 1.787 2 0.627 9 x4 1.334 6 0.441 5 相对植距
RSα1 4.589 9 0.369 5 0.979 9 10.046 1 20.369 4 α2 −0.300 8 0.047 4 α3 −3.354 2 2.107 2 x1 0.650 9 0.199 8 x2 0.453 9 0.188 1 x3 0.342 7 0.182 2 x4 0.452 2 0.177 5 林分密
度指数
SDIα1 3.440 2 0.236 3 0.979 4 10.431 1 20.630 9 α2 −0.326 2 0.061 4 α3 0.000 5 0.000 2 x1 0.721 9 0.228 7 x2 0.524 8 0.204 4 x3 0.418 8 0.192 7 x4 0.472 4 0.183 1 Nilson密
度指数
SDα1 3.825 1 0.259 6 0.979 0 10.613 2 20.811 7 α2 −0.344 3 0.075 5 α3 0.000 5 0.000 2 x1 0.448 3 0.188 8 x2 0.335 0 0.189 9 x3 0.289 9 0.190 2 x4 0.404 2 0.185 2 不含密
度指标
No densityα1 3.512 5 0.038 8 0.972 8 10.925 1 21.123 1 α2 −0.180 5 0.001 9 x1 0.379 9 0.190 4 x2 0.330 8 0.194 5 x3 0.310 9 0.195 5 x4 0.417 3 0.190 8 Table 3. Parameter estimation and model evaluation of stand volume model with different density indices
从模型决定系数方面来看,所有模型的R2数值由高到低顺序为:包含每公顷株数N的林分蓄积量模型(0.979 9)、相对植距模型(0.979 9)、林分密度指数SDI模型(0.979 4)、优势高营养面积比Z模型(0.979 3)、Nilson密度指数模型(0.979 0)以及不含密度指标模型(0.972 8)。除去因为参数估计异常而被舍弃的2号和3号模型,杉木林分蓄积量模型中表现最好的是以林分密度指数SDI为密度指标的模型。同时,SDI模型也取得了最小的平均绝对偏差和均方根误差。因此,从本研究结果中可以认为,林分密度指数SDI模型是估计杉木林分蓄积量的最佳模型。而不含密度指数模型则是表现最差的林分蓄积量模型,这也说明了在拟合和估计杉木林分蓄积量时,缺少密度指数将会影响模型的准确度。SD指数模型也因为在所有密度指数模型中表现最差,在估计杉木林分蓄积方面成为最不适合的密度指标。我们还发现,根据哑变量参数的估计值可得,在低造林密度(1 667~3 333 株·hm−2),林分蓄积生长量要大于中高造林密度的林分(5 000~10 000 株·hm−2)。
基于研究得到的预估杉木林分蓄积量的最佳模型,SDI模型,得到了模型残差分布图(图1)。从图1可以看到,都倾向于低估杉木的林分蓄积量。在林分蓄积量处于0~100 m3·hm−2的时候,残差点基本分布于0~−20的范围内;当实际林分蓄积量达到100 m3·hm−2以上时,残差点分布范围则有所扩大,基本上均匀分布在−10~−50 m3·hm−2的区间中。当然在实际林分蓄积量处于1~10 m3·hm−2和较大蓄积的时候,模型也有少量高估蓄积量的情况发生。
Develop Annual Stand Volume Growth Model of Chinese fir Including Different Stand Density Indices
- Received Date: 2021-04-13
- Accepted Date: 2022-04-11
- Available Online: 2022-10-20
Abstract: