林业科学  2015, Vol. 51 Issue (4): 78-88   PDF    
DOI: 10.11707/j.1001-7488.20150410
0

文章信息

吴恒, 党坤良, 田相林, 孙帅超, 陈书军, 赵鹏祥, 曹田健
Wu Heng, Dang Kunliang, Tian Xianglin, Sun Shuaichao, Chen Shujun, Zhao Pengxiang, Cao Tianjian
秦岭林区天然次生林与人工林立地质量评价
Evaluating Site Quality for Secondary Forests and Plantation in Qinling Mountains
林业科学, 2015, 51(4): 78-88
Scientia Silvae Sinicae, 2015, 51(4): 78-88.
DOI: 10.11707/j.1001-7488.20150410

文章历史

收稿日期:2014-06-10
修回日期:2014-10-20

作者相关文章

吴恒
党坤良
田相林
孙帅超
陈书军
赵鹏祥
曹田健

秦岭林区天然次生林与人工林立地质量评价
吴恒1, 2, 党坤良1, 田相林1, 2, 孙帅超1, 2, 陈书军1, 3, 赵鹏祥1, 曹田健1, 2     
1. 西北农林科技大学林学院 杨凌 712100;
2. 生态仿真优化实验室 杨凌 712100;
3. 陕西秦岭森林生态系统国家野外科学观测研究站 宁陕 711600
摘要:【目的】 以秦岭林区典型的松栎林带为例,采用临时样地、固定样地和解析木3套数据建立立地指数模型和立地形模型,比较天然次生林与人工林立地质量评价差异,为天然次生林立地质量评价提供科学依据和评价方法。【方法】 首先对样地直径分布进行SW正态性检验和株数积累分布曲线检验,论证天然次生林编制立地指数表的可行性; 然后拟合优势树高生长方程,根据R2和标准估计误差选择各树种导向曲线模型; 依据树高连年生长量曲线与平均生长量曲线相交且趋于平缓、树高变异系数趋于平稳确定基准年龄和基准胸径; 分别采用相对优势高法、标准差调整法和变动系数调整法形成立地指数和立地形曲线簇,编制油松、华山松、锐齿栎和落叶松立地指数表和立地形表。运用卡方检验、落点检验和相关性检验分析不同树种立地指数表和立地形表精度差异, 采用平均误差、平均绝对误差和平均相对误差比较通用地位级表与立地指数表、立地形表评价立地质量偏差。【结果】 各树种卡方检验值均小于卡方检验临界值,符合精度要求。油松和落叶松立地指数卡方检验平均值低于立地形,锐齿栎立地指数卡方检验平均值高于立地形,华山松立地指数与立地形卡方检验平均值差异较小。卡方检验、落点检验和相关性检验综合分析结果表明,立地形适用于天然次生林立地质量评价,而立地指数适用于人工林。与通用地位级表的比较表明,地位级评价立地质量平均误差、平均绝对误差、平均相对误差大于立地指数和立地形,地位级评价落叶松人工林立地质量误差显著大于立地指数。【结论】 编制的油松、落叶松立地指数表,华山松、锐齿栎立地形表满足立地质量评价精度要求,能够反映林区天然次生林和人工林立地质量差异,可用于林区立地质量评价和森林经营管理。同时如何解决次生林立地质量评价中混交和异龄的特点、建立更具普遍性的立地指数和立地形模型仍是一个值得深入研究的问题。
关键词立地质量评价    天然次生林    人工林    立地指数    立地形    秦岭    
Evaluating Site Quality for Secondary Forests and Plantation in Qinling Mountains
Wu Heng1, 2, Dang Kunliang1, Tian Xianglin1, 2, Sun Shuaichao1, 2, Chen Shujun1, 3, Zhao Pengxiang1, Cao Tianjian1, 2     
1. College of Forestry, Northwest A & F University Yangling 712100;
2. Laborary of Ecological Optimization of Simulation Yanling 712100;
3. Qinling National Forest Ecosystem Research Station, Shaanxi Province Ningshan 711600
Abstract: 【Objective】 Site quality is essential for stand development. Stands established with natural or artificial regeneration methods requiring effective approaches to evaluate the site productivity in forestry practice. The purpose of this study was to evaluate the site quality for secondary forests and plantation using site index and site form models. 【Method】The permanent and temporary inventory data, and the stem analysis data, were collected from pine-oak stands in Qinling Mountains. Feasibility of applying site index models for secondary forests was validated by cumulative frequency curve and SW normality tests. Reference age and DBH were determined by analyzing of mean height increment and current annual height increment. Methods for formulating site index and site form tables, such as ratio method, standard deviation adjustment, and coefficient variation adjustment were applied by species after fitting dominant height growth equations. Chi-Square test, falling point test and correlation test were used for accuracy analysis of site index and site form for secondary forest and plantation. MRES, AMRES, MRES% were applied for error analysis of site class, site index and site form. 【Result】Chi-square values of each species were lower than critical values. Chi-square values of site index of Pinus tabulaedormis and Larix principis-rupprechtii were lower than that of site form. Chi-square value of site index of Quercus aliena var. acuteserrata was higher than that of site form. Little difference among site index and site form of Pinus armandii was observed. In general, site form appears superior for secondary forests, while site index is suitable for plantation. Site index and site form were more accurate than site class in evaluating the site quality. 【Conclusion】 Establishing site index tables and site form tables could be used for site evaluating in forest management. Further study of mixed and uneven-aged secondary forest site quality evaluation is required.
Key words: forest site evaluation    secondary forest    plantation    site index    site form    Qinling Mountains    

森林立地是森林生产力的基础,对森林更新、树种选择、地力维持和经营管理至关重要。在造林规划设计、森林抚育经营和森林经营规划等林业生产实践中,立地质量评价是基础,科学的立地质量评价方法对提高林分生长收获预估准确性、优化抚育经营设计、实现森林可持续经营具有极其重要的意义。

林分生产力的立地质量评价方法可以分为生物因子法和地理因子法两大类,其中生物因子法的地位级和立地指数是立地质量评价最常用的方法(Batho et al.,2006)。立地指数多用于人工林立地质量评价,少见用于天然林。研究表明,天然次生林多为相对同龄林,亦可用立地指数评价其立地质量(Hann,1998孟宪宇等,2001郑聪慧等,2013);但是,由于天然林树种组成和林分结构的复杂性,天然林立地质量评价更多采用立地形法(Vanclay et al.,1988Huang et al.,1993马建路等,1995陈永富等,2000)。

采用不同方法评价立地质量各有优缺点。地理因子法易于分类(刘建军等,1996),但缺乏作为立地条件影响林分生长的生物学解释。地位级法简便易行(李悦黎等,1993),但其精度和准确性相对低于立地指数。立地指数的无偏估计要求准确的年龄测量值,人工林根据造林记录可查出准确的林分年龄,天然次生林则需要通过年轮分析法确定。立地指数评价天然次生林立地质量时易受林分结构的影响(Lynch,1958Cieszewski et al.,1993MacFarlane et al.,2000),影响的大小因树种差异而有所不同(Flewelling et al.,2001)。相较基于年龄和树高的地位级和立地指数而言,立地形评价立地质量时对基准胸径的选择较为敏感(Goelz et al.,1992),立地形取决于胸径和树高的关系。人工林由于胸径与树高的关系受造林密度、经营措施等的影响(Wang,1998),使得立地形的评价效果较差。

不同立地质量评价方法特点各异,采用单一评价方法容易忽略不同方法对林分立地质量评价造成的误差,系统比较不同方法在人工林和天然次生林立地质量评价中的差异能加深对森林立地质量的认识,提高立地评价的准确性。

秦岭林区面积占陕西省有林地总面积的35.8%,森林覆盖率达56.3%,天然林面积和蓄积分别占林区总面积和总蓄积的86.7%和98.3%,是天保工程的重要组成部分(韩福利,2006)。然而,该林区长期以来缺少对森林立地质量的准确评价,严重制约了林区森林抚育经营技术的发展与提高。因此,本研究选取秦岭林区具有代表性的松栎林带,比较不同森林立地质量评价方法,编制普适性广、应用性强的地位指数表和立地形表。通过比较不同编表方法及统计学检验,证明多元化评价森林立地质量的科学性和必要性。

1 研究区概况与研究方法 1.1 研究区概况

研究地区位于西北农林科技大学火地塘教学试验林场。该林区地处秦岭南坡中段宁陕县境内,林区地形复杂,山势东高西低,坡度在20°~25°之间,海拔1 420~2 474 m。林区年降水量1 000 mm,年积温2 200~3 100 ℃,属于北亚热带气候类型。土壤主要有山地棕壤、暗棕壤和山地草甸土。油松(Pinus tabulaeformis)林分分布在海拔1 000~2 000 m地带,在坡度较大、土层较薄的山坡、山脊占绝对优势并形成纯林,有一定的天然更新能力,是该林区最常见的针叶林。华山松(Pinus arm and ii)林分分布在海拔1 400~2 300 m地带,上限多与青杄(Picea wilsonii)、秦岭冷杉(Abies chensiensis)、红桦(Betula albo-sinensis)和铁杉(Tsuga chinensis)混交,下限多与锐齿栎(Quercus aliena var. acuteserrata)和油松等混交,常见于山坡、谷底或溪旁。锐齿栎林分在海拔1 300~1 900 m形成优势群落,多与松属、其他阔叶树混交,生产力较高。华山松和锐齿栎林天然更新均表现较好。落叶松人工林主要是人工更新的华北落叶松(Larix principis-rupprechtii),分布在海拔1 000~2 400 m之间。

1.2 数据来源和特征

本研究采用临时样地(圆形样方半径10 m)、固定样地(方形20 m×20 m)和解析木多源数据。其中油松临时样地378块,固定样地6块,解析木10株; 华山松临时样地1 434块,固定样地20块,解析木6株; 落叶松临时样地203块,固定样地34块,解析木7株; 锐齿栎临时样地944块,固定样地16块,解析木22株。解析木数据除来源于西北农林科技大学火地塘教学试验林场外,还来源于陕西省宁东林业局旬阳坝林场、新矿林场和高桥林场。龄阶按5年、径阶按4 cm整化,以各龄阶、径阶优势木平均高为准,±3倍标准差为界剔除该龄阶或径阶范围之外的异常点后进行描述性统计(表 1)。

表 1 各树种优势木特征值 Tab.1 Characteristics of dominant trees by tree species
1.3 编制天然次生林立地指数表可行性验证

在同龄纯林中,林分直径结构近似遵循正态分布,平均直径在株数积累分布曲线上在55%~64%范围内,一般近似于60%(孟宪宇,2006)。孟宪宇等(2001)研究认为,天然次生林的径阶分布如果与正态分布没有显著性差异,便可视为“相对同龄林”,该方法的可行性可通过直径频次直方图来验证(郑聪慧等,2013)。本文按照单个树种株数占样地总株数65%以上的标准选取了5块油松固定样地、8块华山松固定样地和7块锐齿栎固定样地进行直径分布Shapiro-Wilk统计量正态性检验,用三次抛物线方程拟合肩形曲线,计算相对直径为1时的株数积累百分数。当显著性水平为0.05时,有3块油松样地、4块华山松样地和2块锐齿栎样地直径结构与正态分布没有显著性差异,17块样地株数积累百分数在55%~64%范围内(表 2),表明本研究中的天然次生林符合“相对同龄林”特征。

表 2 天然次生林直径结构正态性检验和株数积累百分数检验 Tab.2 Tests of diameter distribution and accumulation for secondary forests
1.4 立地指数和立地形模型的建立 1.4.1 优势木选择

常用的优势木选择方法是选择林分中最高的林木作为优势木,并要求编表与使用表时选测方法一致(孟宪宇,2006)。本文对20 m×20 m固定样地选择4株最高的优势木作为高选样本(鄢志明等,1983); 对临时样地在半径10 m的圆形样方内选择3株最高的优势木作为高选样本(孟宪宇,2006詹昭宁等,1978)。将异常数据剔除后的样本数据80%用于导向曲线拟合,20%用于检验。

1.4.2 导向曲线的拟合

导向曲线的选择直接影响模型对立地质量评价的准确性,这就需要导向曲线的形式既符合树高生长的生物学规律,又要能对数据进行最优化拟合。良好的导向曲线应该呈平滑的“S”形,并且具有上限渐近线(Weiskittel et al.,2011)。常用的导向曲线有抛物线式、双曲线式、单分子式和理查兹式等,本文采用了模型1-5(李海奎等,2011)进行胸径-树高生长拟合、模型6-19(孟宪宇,2006Perin et al.,2013)进行年龄-树高生长拟合(表 3)。拟合时为了不减少模型的自由度,不采用各龄阶(或径阶)年龄(或胸径)与树高平均值拟合,而采用所有建模数据点进行拟合,根据决定系数R2、SEE选择导向曲线模型。

表 3 树高生长备选模型 Tab.3 Models for tree height growth
1.4.3 基准年龄、基准胸径和指数级距的确定

基准年龄和基准胸径对立地指数和立地形模型的影响十分显著,选择不当会对立地质量评价造成偏差。本研究在确定基准年龄时,利用解析木资料分析树高的生长过程,分别对不同立地条件下解析木的连年生长量和平均生长量求平均值并绘制曲线图;同时计算建模样本各龄阶树高变异系数CV和树高变异系数变化幅度ΔCV并绘制曲线图,根据曲线图上树高生长趋于平缓且能灵敏反映立地质量的原则确定基准年龄。基准胸径的确定方法与基准年龄的确定方法一致。指数级距C根据基准年龄和基准胸径时的树高绝对变幅和指数级个数概算。公式如下:

$ CV = {S_{\rm{i}}}/{H_i}; $ (1)
$ \Delta CV = C{V_{i + 1}}/C{V_i}; $ (2)
$ C = \Delta H/k。 $ (3)

式中: CV为变异系数; ΔCV为变异系数变化幅度; Si为第i龄阶树高标准差;Hi为第i龄阶主曲线树高值; C为指数级距; ΔH为基准年龄时的树高绝对变幅; k为指数级个数。

1.4.4 立地指数表和立地形表的编表方法

立地指数常用的编表方法有标准差调整法[st and ard deviation adjustment,SDA,式(4)]、变动系数调整法[coefficient variation adjustment,CVA,式(5)]和相对优势高法[ratio method/RM,式(6)](孟宪宇,2006)。标准差调整法多用于阔叶树种,而变动系数调整法普遍用于针叶树种,相对优势高法精度相对较低(陶吉兴等,1990)。立地形的调整方法与立地指数一致(马建路等,1995陈永富等,2000)。本文对每个树种分别采用3种方法建立立地指数模型,将年龄换成胸径、龄阶换成径阶同理建立立地形模型:

$ {H_{ij}} = {H_{ik}} \pm [({H_{0j}} - {H_{0k}})/{S_{A0}} \bullet {S_{Ai}}]; $ (4)
$ {H_{ij}} = {H_{ik}}[1 \pm [({H_{0j}} - {H_{0k}})/({H_{0k}} \bullet {C_{Ai}})]; $ (5)
$ {H_{ij}} = ({H_{0j}}/{H_{0k}} \times 100\% ) \bullet {H_{ik}}。 $ (6)

式中: Hij为第i龄阶第j指数级调整后的树高; Hik为第i龄阶的导向曲线树高; H0j为基准年龄时第j指数级的树高; H0k为基准年龄时导向曲线树高; SA0为基准年龄所在龄阶树高标准差理论值; SAi为第i龄阶树高标准差理论值; CA0为基准年龄所在龄阶树高变动系数理论值; CAi为第i龄阶树高变动系数理论值。

1.5 模型的评价与检验

对编制的立地指数表和立地形表,运用检验样本进行卡方检验、落点检验和相关性检验。卡方检验根据检验样本的年龄(或胸径)和优势树高值查出对应的立地指数(或立地形),再求出该龄阶(或径阶)的树高理论值Ĥ,然后与实际树高值H进行卡方检验,检验公式为:

$ {\chi ^2} = \sum [{(H - \hat H)^2}/\hat H]。 $ (7)

落点检验将检验样本做散点图,绘制到立地指数和立地形曲线簇中,算出散点落在曲线簇内的概率,即为立地指数表或立地形表能够解释优势木生长的概率值。运用SPSS19.0软件对优势树高理论值与实际值进行相关性分析。卡方检验值作为评价模型的主要指标,根据临界值表判断模型精度是否符合要求。卡方检验值均低于临界值时,检验值越小、落点检验概率值越大和相关系数越大,表明树高理论值与实际值之间的偏离程度越小,模型精度就越高。

根据各树种的小班平均年龄和平均高查秦巴林区通用地位级表,得到对应的地位级,再求出平均树高理论值,油松、华山松和落叶松采用针叶树地位级表,锐齿栎采用阔叶树地位级表。同样的方法得到小班树高理论值,用实际值H与理论值Ĥ进行误差统计量分析,检验统计量包括平均误差(MRES)、平均绝对误差(AMRES)和平均相对误差(MRES%),检验公式为:

$MRES = \sum (H - \hat H)/n; $ (8)
$ AMRES = \sum (|H - \hat H|)/n; $ (9)
$ MRES\% = 100 \times \sum (H - \hat H)/\sum H。 $ (10)

平均误差绝对值和平均绝对误差小于0.5、相对误差绝对值小于5%,说明树高理论值与真实值的偏差小于0.5 m,模型检验值符合精度要求。误差绝对值越小,模型精度越高。

2 结果与分析 2.1 立地指数和立地形模型拟合结果

导向曲线拟合结果(表 4)表明,落叶松人工林拟合决定系数R2高于其他天然次生林,这与人工林数据分布比较集中、天然次生林数据分布比较分散有关。各树种的拟合决定系数R2均低于0.9,这与模型拟合时采用所有数据点而非各龄阶(或径阶)年龄(或胸径)和树高平均值有关。

表 4 各树种导向曲线拟合结果 Tab.4 Fitting results of guide curves
2.2 基准年龄、基准胸径与指数级距

以油松为例,当年龄30年左右、胸径18 cm左右时,树高的连年生长量曲线与平均生长量曲线相交,之后生长趋于平缓,树高变异系数趋于平稳,变异系数变化幅度趋于1,从而确定油松的基准年龄为30年,基准胸径为18 cm。同样的方法确定华山松、落叶松和锐齿栎的基准年龄分别是30,20和20年,基准胸径分别是20,16和12 cm。根据基准年龄和基准胸径时的树高变化幅度确定油松、华山松、落叶松和锐齿栎的指数级距分别为2,2,1和2 m,指数级个数都为5。以油松、落叶松立地指数和华山松、锐齿栎立地形为例展开得到立地指数和立地形曲线簇,见图 1

图 1 油松(a)和落叶松(b)立地指数、华山松(c)和锐齿栎(d)立地形曲线簇 Fig. 1 Examples of site index and site form curves by tree species
2.3 立地指数与立地形模型检验评价结果

表 5比较了立地指数表和立地形表的统计检验结果。各树种卡方检验值均小于临界值,油松、落叶松卡方检验临界值为77.9(df=60-1,α=0.05),华山松卡方检验临界值为178.5(df=150-1,α=0.05),锐齿栎卡方检验临界值为1399(df=115-1,α=0.05),符合精度要求。油松、落叶松和华山松立地指数卡方检验平均值分别比立地形低2.83,4.56和0.17,锐齿栎立地指数卡方检验平均值比立地形高0.80。油松、落叶松和锐齿栎立地指数表落点检验平均值分别高立地形表 3.89%,10.55%和3.19%,华山松落点检验立地指数表与立地形表差异较小,平均值差为1.56%,相对优势高法落点检验平均值分别高标准差调整法和变异系数调整法2.45%和1.08%。油松和落叶松立地指数表相关系数平均值分别高于立地形表0.018和0.014,华山松差异较小为0.002,锐齿栎立地指数表相关系数略低于立地形表为0.004,标准差调整法和变异系数调整法的相关系数差异不大且都大于相对优势高法,平均值分别为0.976,0.977和0.969。

表 5 各树种立地指数与立地形模型卡方、落点和相关性统计检验 Tab.5 Chi-square, falling point, and correlation statistics in site index and site form models by species

以卡方检验值为主要指标,综合考虑落点检验和相关性检验结果确定林分立地质量评价方法分别是落叶松和油松采用立地指数、华山松和锐齿栎采用立地形; 编表调整方法分别是油松和锐齿栎选择CVA、华山松和落叶松选择SDA;导向曲线选取分别是油松选取Gaussian、华山松选取Korf、落叶松选取Sigmoidal、锐齿栎选取Weibull。编成油松、落叶松立地指数表和华山松、锐齿栎立地形表见附表 1~4

附表 1 油松立地指数表(基准年龄: 30年) Appendix 1 Site index table of P.tabulaeformis(reference age: 30 a)
附表 2 华山松立地形表(基准胸径: 20cm) Appendix 2 Site form table of P.armandii(reference DBH: 20 cm)
附表 3 落叶松立地指数表(基准年龄: 20年) Appendix 3 Site index table of L.principis-rupprechtii(reference age: 20 a)
附表 4 锐齿栎立地形表(基准胸径: 12 cm) Appendix 4 Site form table of Q.aliena var.acuteserrata(reference DBH: 12 cm)
2.4 不同立地质量评价方法的比较

立地指数、立地形与地位级的统计误差分析如表 6所示。从表 6可以看出,油松和落叶松地位级平均误差(MRES)大于立地指数平均误差,华山松和锐齿栎地位级与立地形平均误差差异较小; 立地指数和立地形平均绝对偏差(AMRES)均在0.54以下,地位级在0.8左右,落叶松最大为1.614; 立地指数和立地形平均相对误差(MRES%)均小于1%,落叶松地位级相对误差最大为15.5%。结果表明,地位级评价立地质量的偏差大于立地指数和立地形,落叶人工林采用地位级的立地质量评价方法误差明显大于立地指数。通用地位级表不区分树种、林分平均高也受树种组成等影响都可能降低地位级评价立地质量的准确性。

表 6 立地指数、立地形与地位级误差检验 Tab.6 Fit statistics for MRES, AMRES, MRES% of site index, site form, and site class
3 结论与讨论

本研究采用固定样地、临时样地和解析木高选样本数据,运用3种编表方法对林区4个主要树种编制立地指数表和立地形表。油松和落叶松适合采用立地指数表,华山松和锐齿栎适合采用立地形表,各树种地位指数表符合统计检验要求,能用于林分立地质量的评价。然而,立地质量评价准确性会受建模数据、评价方法等因素的影响。

3.1 数据对立地质量评价模型的影响

林分最高木与最粗木之间存在着极显著相关性,高选和径选方法之间存在着高度的统一(赵美丽等,1994)。鄢志明等(1983)在编制油松地位指数表时对高选样本和径选样本的差异性检验表明,2种方法选择得到的样本没有显著性差异,树高相对于胸径不受林分密度的影响,而胸径则与密度相关。因此,本研究未考虑径选而直接采用高选优势木。

解析木数据虽被广泛地用于立地指数模型(Curtis,1964),但也存在其局限性。使用解析木建立的优势树高生长模型是否会高估树高生长仍有争议(Raulier et al.,2003Perin et al.,2013)。固定样地数据是建模最理想的数据,因为它能比较样地组内和组间的统计值去验证模型的准确性。与解析木数据相比,固定样地数据拟合的优势树高生长曲线会更早地达到渐近线水平(Raulier et al.,2003);但实际中由于各种原因较难获取大量固定样地的数据,临时样地数据亦可用于生长建模。单独使用临时样地数据会造成模型的估计偏差,Perin等(2013)在比较挪威云杉(Picea abies)的树高生长时发现,使用临时样地建立的树高生长模型会低估树高生长。由于临时样地数据是在短期内测量完成的,使得测量数据更容易受极端气候变化的影响,Keen等(1937)发现干旱地区西黄松(Pinus ponderosa)对降雨极其敏感,从而导致优势树高的数据误差。本文采用Garcia(2005)Perin等(2013)基于多样数据的方法编制立地位指数表和立地形表,力图降低单一数据类型造成的立地质量评价误差;然而,本研究调查数据集中于秦岭南坡中段,有待于补充秦岭东西段和北坡数据,以提高本文模型的精度和普适性。

3.2 导向曲线和基准胸径对立地质量评价模型的影响

人工林常采用各龄阶年龄和优势树高的平均值进行导向曲线拟合(孟宪宇,2006段劼等,2009),平均年龄和优势木平均高拟合能显著提高决定系数R2,但也极大地减少了模型的自由度;特别是天然次生林由于优势木年龄不一致,以及数据的分布较人工林存在明显差异,采用各龄阶平均年龄和优势木平均高建模更容易造成模型估计误差。本文采用所有建模数据点进行拟合,提高了树高生长模型的自由度。

立地形作为评价混交林和异龄林的有效方法,需要确定一个树高生长趋于稳定且能灵敏反映立地质量差异的基准胸径(Weiskittel et al.,2011)。根据林木调查数据中出现频次较多的胸径值为基准胸径的方法(Vanclay et al., 1988)缺少合理的生物学解释。令胸径-树高方程二阶导数为零时得到的基准胸径(陈永富,2000)是树高连年生长量最大时的胸径,根据基准年龄对应的胸径值确定基准胸径(Huang et al.,1993)则需要准确的年龄-胸径方程。本文根据树高生长量趋于稳定且能灵敏反映立地质量差异时对应的胸径为基准胸径(马建路等,1995)满足反映立地质量差异的要求;但由于缺乏固定样地和解析木数据,特别是落叶松和锐齿栎解析木数据的短缺,使得这2个树种基准年龄和基准胸径确定的准确性受到影响。

3.3 立地质量评价模型检验方法

如果对立地指数表的落点检验采用编表数据进行检验,容易引起检验样本和编表样本的耦合(段劼等,2009)。在立地质量评价模型检验中,通用式卡方检验不适合检验连续型随机变量的实际值与理论值的拟合优度(吴再明,1988)。本文将样本数据分为编表样本和检验样本,其中检验样本占总体样本的1/5,这样可降低编表样本数据检验立地指数表造成的耦合现象。本文采用的检验样本是由固定样地数据、临时样地数据和解析木数据随机选取,可用通用式卡方检验进行有效检验。本文可在相同显著水平上,比较预测值与真实值的相关系数大小,也可比较各种立地质量评价方法的优劣。

3.4 不同方法对立地质量评价准确性的影响

油松和锐齿栎立地指数曲线簇与地位级曲线簇相比,成熟林时地位级曲线树高生长仍然处于上升状态,相对高估了成熟林时树高生长(图 2)。误差统计量分析结果表明,通用地位级表评价立地质量的误差大于立地指数或立地形表,其中落叶松人工林采用地位级法评价立地质量误差最大(表 6)。采用林分平均高或者优势高评定立地质量均会受人为干扰而产生误差,地位级法不适用于采用下层伐的林分,立地指数和立地形不适用于使用上层伐的林分(孟宪宇,2006)。秦岭林区在禁伐令前多有间伐作业,然而,历史间伐强度和间伐方式等准确数据难以获得,这在一定程度上也影响了本文立地质量评价的准确性。

图 2 地位级与立地指数曲线簇 Fig. 2 Curves of site class and site index
参考文献(References)
[1] 陈永富,杨彦臣,张怀清,等. 2000. 海南岛热带天然山地雨林立地质量评价研究. 林业科学研究, 13(2):134-140.
(Chen Y F, Yang Y C, Zhang H Q, et al. 2000. A study on site quality evaluation of natural tropical mountainous rain forest in Hainan Island. Forest Research, 13(2):134-140 [in Chinese]).(3)
[2] 段 劼,马履一,贾黎明,等. 2009.北京低山地区油松人工林立地指数表的编制及应用.林业科学, 45(3):7-12.
(Duan J, Ma L Y, Jia L M, et al. 2009. Establishment and application of site index table for Pinus tabulaeformis plantation in the low elevation area of Beijing. Scientia Silvae Sinicae, 45(3):7-12 [in Chinese]).(2)
[3] 韩福利. 2006.秦岭林区森林资源保护的意义和可持续经营利用. 陕西林业科技, (4):69-72.
(Han F L. 2006. The significance of forest resources conservation in Qinling Mountainous area as well as sustained management and utilization. Shaanxi Forest Science and Technology, (4):69-72 [in Chinese]).(1)
[4] 李海奎,法 蕾. 2011.基于分级的全国主要树种树高-胸径曲线模型.林业科学, 47(10): 83-90.
(Li H K, Fa L. 2011. Height-diameter model for major tree species in China using the classified height method. Scientia Silvae Sinicae, 47(10):83-90 [in Chinese]).(6)
[5] 李悦黎,杜纪山. 1993.火地塘教学实验林场森林资源的数据分析及经营对策.西北林学院学报, 8(3):53-58.
(Li Y L, Du J S. 1993. Analysis and management strategy of forest resources of Huoditang teaching and experimental forest farm. Journal of Northwest Forestry University, 8(3):53-58 [in Chinese]).(1)
[6] 刘建军,雷瑞德,尚廉斌,等. 1996.火地塘林场立地分类的初步研究.西北林学院学报, 11(增):31-36.
(Liu J J, Lei R D, Shang L B, et al. 1996. Site classification at Huoditang forest farm. Journal of Northwest Forestry University, 11(S1):31-36 [in Chinese]).(1)
[7] 马建路,宣立峰,刘德君. 1995.用优势树全高和胸径的关系评价红松林的立地质量. 东北林业大学学报, 23(2):20-27.
(Ma J L, Xuan L F, Liu D J. 1995. Site quality estimation for natural Korean pine forest using total height and diameter of dominant tree. Journal of Northeast Forestry University, 23(2):20-27 [in Chinese]).(3)
[8] 孟宪宇. 2006.测树学. 北京: 中国林业出版社.
(Meng X Y. 2006. Forest mensuration. Beijing: Chinese Forestry Publishing House. [in Chinese])(20)
[9] 孟宪宇,陈东来. 2001.山杨次生林地位指数表编制方法的研究. 北京林业大学学报, 23(3):47-51.
(Meng X Y, Chen D L. 2001. Study on constructive method of site index table for poplar secondary forest. Journal of Beijing Forestry University, 23(3):47-51 [in Chinese]).(2)
[10] 陶吉兴,唐明荣. 1990.常规立地指数表的误差来源与分析. 浙江林学院学报, 7(4):99-103.
(Tao J X, Tang M R. 1990. A source and analysis of errors in common site index table. Journal of Zhejiang Forestry College, 7(4):99-103 [in Chinese]).(1)
[11] 吴再明. 1988.论X2检验在地位指数曲线检验中的应用. 云南林业调查规划, (2):1-6.
(Wu Z M. 1988. Application of Chi-square test in site index curve. Yunnan Forest Inventory and Planning, (2):1-6. [in Chinese])(1)
[12] 鄢志明,彭志鹏. 1983. 油松地位指数表的编制. 陕西林业科技, (4):7-12.
(Yan Z M, Peng Z P. 1983. Site index table establishment of Pinus tabulaeformis. Shaanxi Forest Science and Technology, (4):7-12 [in Chinese]).(2)
[13] 詹昭宁,郝祖渊. 1978.林分中三株最高树树高平均值与蓄积量相关关系的初步探讨. 林业资源管理, (4):9-12.
(Zhan Z N, Hao Z Y. 1978. Relationship analysis of stand volume and average top height of 3 dominant trees. Forest Resources Management, (4):9-12 [in Chinese]).(1)
[14] 赵美丽,王才旺. 1994.林分优势高测定方法的探讨. 内蒙古林业调查设计, (2):10-16.
(Zhao M L, Wang C W. 1994. Methods of dominant height mensuration. Inner Mongolia Forestry Investigation and Design, (2):10-16 [in Chinese]).(1)
[15] 郑聪慧,贾黎明,段 劼,等. 2013.华北地区栓皮栎天然次生林地位指数表的编制. 林业科学, 49(2):79-85.
(Zheng C H, Jia L M, Duan J, et al. 2013. Establishment of site index table for Quercus variabilis natural secondary forest in north China. Scientia Silvae Sinicae, 49(2):79-85 [in Chinese]).(2)
[16] Batho A, Garcia O. 2006. De Perthuis and the origins of site index: a historical note. FBMIS, 1: 1-10.(1)
[17] Cieszewski C J, Bella I E. 1993. Modeling density-related Iodgepole pine height growth, using Czarnowski's stand dynamics theory. Canadian Journal of Forest Research, 23(12): 2499-2506.(1)
[18] Curtis R O. 1964, A stem-analysis approach to site-index curves. Forest Science, 10(2): 241-256.(1)
[19] Flewelling J W, Collier R, Gonyea B, et al. 2001. Height-age curves for planted stands of Douglas fir, with adjustments for density. College of Forest Resources, University of Washington.(1)
[20] Garcia O. 2005. Comparing and combining stem analysis and permanent sample plot data in site index models. Forest Science, 51(4): 277-283.(1)
[21] Goelz J C G, Burk T E. 1992. Development of a well-behaved site index equation: jack pine in north central Ontario. Canadian Journal of Forest Research, 22(6): 776-784.(1)
[22] Hann D W. 1998. Extending southwest Oregon's Douglas-fir dominant height growth equation to older ages.(1)
[23] Huang S, Titus S J. 1993, An index of site productivity for uneven-aged or mixed-species stands. Canadian Journal of Forest Research, 23(3): 558-562.(2)
[24] Keen F P. 1937. Climatic cycles in eastern Oregon as indicated by tree RINGS1. Monthly Weather Review, 65: 175.(1)
[25] Lynch D W. 1958. Effects of stocking on site measurement and yield of second-growth ponderosa pine in the Inland Empire. Intermountain Forest and Range Experiment Station.(1)
[26] MacFarlane D W, Green E J, Burkhart H E. 2000. Population density influences assessment and application of site index. Canadian Journal of Forest Research, 30(9): 1472-1475.(1)
[27] Perin J, Hébert J, Brostaux Y, et al. 2013. Modelling the top-height growth and site index of Norway spruce in Southern Belgium. Forest Ecology and Management, 298: 62-70.(17)
[28] Raulier F, Lambert M C, Pothier D, et al. 2003. Impact of dominant tree dynamics on site index curves. Forest ecology and management, 184(1): 65-78.(2)
[29] Vanclay J K, Henry N B. 1988. Assessing site productivity of indigenous cypress pine forest in southern Queensland. Commonwealth Forestry Review, 67(1): 53-64.(2)
[30] Wang G G. 1998. An ecologically based model for site index conversion among species. Canadian Journal of Forest Research, 28(2): 234-238.(1)
[31] Weiskittel A R, Hann D W, Kershaw Jr J A, et al. 2011. Forest Growth and Yield Modeling. Wiley.(2)