广西大青山西南桦人工林立地指数表编制

唐 诚，王春胜，庞圣江，贾宏炎，曾 杰，雷渊才

（1.石河子大学 农学院，新疆 石河子 832003； 2.中国林业科学研究院 热带林业研究所，广东 广州510520； 3.中国林业科学研究院 热带林业实验中心，广西 凭祥536000；4.中国林业科学研究院 资源信息研究所，北京100091）

立地质量评价一直是森林可持续经营的一个主要议题。生产实践中，评价立地质量、预测林分生产力对于树种选择、林分经营决策至关重要［1］。同龄林中，立地指数早已作为立地质量评价的一个指标得到全球公认，且作为立地指数基本产量表的关键变量得以应用［2-3］。在小范围、局域尺度开展立地质量评价时，国内外学者们常采用导向曲线法研制立地指数表或者立地指数曲线，用以评价林木生长状况以及林分生产力［4-7］。西南桦Betula alnoides是中国热带、南亚热带地区的一个珍贵乡土用材树种，生长迅速，适应性强，木材纹理细致，质地均匀，密度适中，不翘不裂，加工性能优良，已广泛应用于木地板和家具制作以及房屋装饰［8-9］，树皮提取物具有消炎、减肥及降血脂作用［10-11］。近年来，中国西南桦迅猛发展，云南、广西、贵州、广东、福建等地均有栽培，其人工林面积已逾15万hm2［9］。然而，在西南桦人工林快速发展过程中，局部地区因未能适地适树，盲目造林，出现林木生长不良、林分生产力低下等现象［12］，严重影响西南桦人工林规模发展和可持续经营。广西大青山是中国西南桦栽培历史最早的地区，自20世纪70年代末即开始西南桦驯化栽培研究，并逐步推广应用，具有从幼龄至成熟龄各个阶段的西南桦林分。西南桦立地指数研究是掌握其立地生产力的基础，对其人工林可持续经营具有重要意义。然而，至今尚未见有关西南桦人工林立地指数表的报道。因此，本研究以广西大青山西南桦人工林为对象，基于样地调查与优势木树干解析，应用导向曲线法编制立地指数表，为该地区及类似地区西南桦人工林立地质量评价、生产力预估提供科学依据。

1 研究地概况

研究地位于广西大青山林区中国林业科学研究院热带林业实验中心（21°57′～22°19′N， 106°39′～106°59′E），属南亚热带季风气候，年均气温为20～22℃，年均降水量约1 500 mm，主要集中在4-10月，年均日照1 260 h。研究地海拔为130～1 045 m，地形地貌复杂，自然条件优越，是开展热带南亚热带地区珍贵阔叶树种试验示范的理想场所。

2 研究方法

2.1 样地调查及解析木数据整理

考虑西南桦林分的立地因子、林分密度等特征，于热带林业实验中心的伏波、青山和白云实验场设置49块临时样地，每样地面积为600 m2。对各样地进行常规立地调查和生长观测，每样地选取5株干形完整、树高最大的树作为优势木，计算优势木平均树高为优势高，并选取1株树高最接近优势高的优势木进行树干解析。

对解析数据按年龄进行整理（表1），共有825对树高-年龄数据，计算平均高及标准差。因有些样木存在不规则生长时期，剔除超出平均值±3倍标准差范围的异常数据，对剩余的816对数据重新整理、统计，用于编制立地指数表［13］。

说明：*优势木在树干解析过程中1～5 a的部分数据缺失

2.2 导向曲线拟合与选择

采用常用的9个模型对树高-年龄数据进行拟合，应用决定系数（R2），平均绝对误差（EAMR），平均相对误差（ERMR）和均方根误差（ERMSE）等4个指标对各模型拟合质量进行评价［14-15］，选择最优曲线作为导向曲线。计算上述统计指标时，将对数双曲线式和对数曲线式方程进行变换，保证所有模型的因变量统一为树高（H）。

2.3 基准年龄与指数级距确定

通过优势高平均生长量和连年生长量变化趋势以及优势高龄阶变动系数分析确定西南桦人工林基准年龄［5］。根据基准年龄时的优势高变幅确定指数级距。应用标准差调整法［13］编制立地指数表。

2.4 立地指数表检验

应用落点检验［16］、拟合显著性和预测精度检验［16-17］对所编立地指数表的精确性和适用性进行检验。运用Excel和R软件（3.3.2版）进行数据整理分析。

3 结果与分析

3.1 导向曲线筛选

由表2可知：9个方程中，对数曲线式的R2最小，其EAMR，ERMR和ERMSE均最大，说明其拟合效果最差；其他8个方程中，Richards和Weibull方程表现几乎一致，以Richards方程的R2最大，EAMR，ERMR和ERMSE最小，拟合效果最好。因此，选择Richards方程作为最优导向曲线方程。其方程表达式为：H=24.25（1-e-0.118t）1.278。 其中：H为优势高，t为年龄。

3.2 基准年龄及指数级距

由图1可知：5 a左右时，优势高连年和平均生长量达到高峰，6 a时两者相交，15 a后优势高连年生长量大致趋于稳定；优势高变异系数的变动幅度在15 a后也基本稳定。因此将西南桦人工林基准年龄确定为15 a。

表2 9个方程表达式及其拟合结果Table 2 Nine equations and their fitting results

中国多以1.0～2.0 m作为指数级距［13］。以15 a作为基准年龄时，优势高为15.2～25.8 m，其绝对变幅为10.6 m，西南桦较为速生，确定其指数级距为2.0 m。可划分为6个指数级（16.0， 18.0， 20.0， 22.0， 24.0 和 26.0 m）。

3.3 立地指数表

应用标准差调整法编制西南桦人工林立地指数表，结果见表3。森林经营中，调查获取林分优势高及年龄，可查此表获得所属立地指数级，进行立地质量评价。

3.4 立地指数表检验

3.4.1 落点检验 将编表所用49株优势木树高值作散点图，绘制于立地指数曲线簇上（图2）。由图2可以看出：3个点落在立地指数曲线簇外，即所编表能解释93.8%的优势木的生长状况，落点检验合格。

3.4.2 显著性检验 利用卡方（χ2）检验进行显著性检验。本研究中立地指数定为16～26，从建模数据中随机抽取3株年龄大于基准年龄的解析木，其立地指数分别为16，20和24 m，χ2计算结果见表 4。查 χ2临界值表分别为 30.14，22.36和30.14。通过对比分析可知：其χ2值远小于其相应临界值，说明经标准差调整后立地指数曲线与优势木树高生长趋势之间无显著差异。

3.4.3 预报精度检验 不同年龄的立地指数估计误差（Si）和不同立地指数级的估计误差（Sz）计算结果见表4。由表4可知：各年龄的估计误差为0.07～1.55 m。3～7 a的估计误差在1.00 m以上，以3和4 a的估计误差为最大；而大于7 a后，各年龄的估计误差基本上小于1.00 m，表明此立地指数表对较大年龄的预报精度高。各立地指数级之间比较，立地指数级为20.0和24.0时估计误差均小于1.00 m，而立地指数级为16.0时略大于1.00 m。由此可见，所编立地指数表整体上预报精度较高。

图1 广西大青山西南桦人工林优势高平均生长量和连年生长量变化曲线Figure 1 Curve of dominant height growth in Betula alnoides plantations at Daqing Mountain，Guangxi

图2 西南桦人工林立地指数曲线落点检验Figure 2 Test of falling point for site index curves of Betula alnoides plantations

4 讨论

利用导向曲线法编制立地指数表，曲线的选取至关重要。导向曲线既要满足对数据的高质量拟合，又要符合树种的优势高生长规律。大多数树种的高生长过程为慢—快—慢，表现为 “S”形；速生树种的高生长常表现为早期生长快，此后生长逐渐减缓并趋于生长停止，即 “快—慢”的 “肩形”曲线。学

者们常用Richards，Weibull等具上渐近线的方程拟合树高生长过程［18-19］，且Richards方程具有很强的可塑性，通过引入形状参数c而使Richards方程对树木生长具有广泛的适应能力。西南桦为速生珍贵树种，本研究中Richards和Weibull方程的拟合效果均较好，以Richards方程表现略优，最终选择Richards方程作为导向曲线方程。

表3 广西大青山西南桦人工林立地指数表Table 3 Site index table of Betula alnoides at Daqing Mountain，Guangxi

表4 立地指数表显著性和预报精度检验Table 4 Significance and prediction accuracy of site index table

西南桦为喜阳树种，大量营林实践及科学研究证实，林龄5～8 a内树高生长迅速且波动性大，而林分优势高的平均和连年生长量则在此后均开始下降且波动性逐渐变小［20-21］。本研究中，6 a左右时西南桦人工林优势高生长达最大值，此后逐渐下降（图1），其波动性也开始减弱，不同年龄时方程拟合误差亦表现出随着树高生长波动性减弱而降低的趋势，因而7 a以上林龄的立地指数估计误差小于7 a内各年龄。

数据源同样会影响模型拟合质量的高低［22］，应用各龄阶平均年龄和平均优势高进行模型拟合，虽然能够显著提高模型的决定系数，但是极大地降低了模型的自由度［6］。本研究针对西南桦树高生长快的生物学特性，以1 a为一个龄阶，用解析木原始数据进行模型拟合，提高了模型自由度，其决定系数约0.90。当然，应用林分优势木的解析数据编表常会造成某些年龄阶段优势高生长的有偏估计，主要是由于林分生长是一个动态过程，临时样地中选取的优势木很难保证其真正一直处于林分优势木地位［23］，从本研究中西南桦优势高生长的波动性规律亦可见一斑。

5 结论

本研究以广西大青山林区西南桦人工林为对象，通过样地调查选取49株平均优势木进行树干解析，获得816对优势高-年龄数据，应用9个常见数学模型进行拟合，筛选出Richards方程作为导向曲线。依据优势高连年生长量和各年龄优势高变异系数的变化趋势，确定西南桦基准年龄为15 a；根据15 a时优势高变动幅度确定以2.0 m为指数级距；应用标准差调整法编制了立地指数表。拟合显著性和预报精度检验表明：立地指数表所反映的优势高生长过程与实际的优势高生长曲线一致，所编立地指数表预报精度较高，尤其是7 a以上林分的预测误差小于1.00 m，此表可用于广西大青山及类似地区西南桦人工林立地质量评价和生长潜力预估。