-
森林立地质量评价是实现科学造林、合理高效利用林地的重要保证。掌握森林立地质量可为森林的可持续发展制定更科学的经营措施及提高森林的经济、生态、社会效益提供理论基础[1]。传统的森林立地质量评价采用的是地位指数或地位级作为重要指标，即根据林分优势木的平均高和林分平均高2种直接评价方法评定立地质量[1-3]。虽然这2种方法在同龄林的质量评价中精度较高，但需要进行大量的实地调查、采伐并制作解析木，导致经济成本和生态成本较高，永久性样地和代表性试验样地缺乏，因而传统森林立地质量评价所需的大尺度数据相对缺乏[3-4]，实用性较差。随着近年来遥感技术应用于立地质量的评价研究[5-6]，大尺度的立地分类逐渐克服了传统森林立地质量评价中生境制图时的耗时、耗人力的缺点。此外，统计学方法用于森林立地质量评价和分类的实践也取得了良好效果，如聚类分析、多元回归分析、数量化理论[7-12]等各种统计方法的应用使得定量研究在立地分类和立地质量评价综合系统的构建中有了更广泛的空间[9-14]。然而，尽管这些方法可以增加森林立地质量评价和分类结果的经济产出，并在一定条件下提高精度，但他们基本上仍沿用传统的评价指标，在成本与效益之间没有突出改善，难以适应大范围森林立地质量评价[15-17]。国家森林资源连续清查数据是覆盖范围广、精度可靠的数据，但在较大范围的森林立地评价研究应用却非常少，相关案例主要是结合遥感影像以样地总生物量为评价指标。而样地总生物量会受年龄、人为活动等多种干扰因素的影响，难以准确评价森林立地质量。本研究利用近期5期浙江省国家森林资源清查（national forest inventory，NFI）样地数据，提取杉木Cunninghamia lanceolata林样地优势木的最大胸径生长率，作为备选立地质量评价指标，结合数字高程模型（DEM）和土壤数据，运用数量化理论Ⅰ方法建立浙江省杉木人工林立地质量评价数量化理论模型。
-
浙江省地处中国东南沿海长江三角洲南翼，属亚热带季风气候，季风显著，四季分明，年气温适中，光照较多，雨量丰沛，空气湿润，雨热季节变化同步，气候资源类型多样，气象灾害繁多。年平均气温15.0~18.0 ℃，年平均降水量为980~2 000 mm，年平均日照时数1 710~2 100 h。
全省陆域面积为1 055×104 hm2，其中山地和丘陵占74.6%。全省林地面积667.97×104 hm2，森林面积584.42×104 hm2（其中85.95×104 hm2为国家特别规定灌木林）；杉木人工林纯林面积96.74×104 hm2，占森林面积的19.43%。全省大部分地区适宜杉木生长，南部属于杉木中心产区。
-
立地因子（包括地形地貌、土壤、气候、植被、生物等）对树木的生长有重要影响。坡度、坡向、海拔、土壤、坡位、湿度、太阳辐射等地形及环境因素会不同程度影响杉木的生长[18-19]，其中地形因子主要包括海拔、坡位、坡向、坡度等。湿度及太阳辐射由于数据较难获取[20]，本研究考虑其与地区所在的纬度关系较大，因此将纬度列为立地因子之一。地形因子可由DEM处理得到，浙江省DEM由地理空间数据云网站上下载得到。DEM提取纬度、海拔、坡位、坡向、坡度等信息在Arc GIS软件中完成。
-
浙江省1989-2009年共5期NFI样地数据，包括样地和样木数据。由于浙江省在1989年后有部分样地位置移动，1989年清查数据仅作分析参考用，实际应用为1994-2009年共4期数据。4期NFI样地数据中，有连续复位样木的杉木人工林样地529个，用于建模；另有用于模型检验的临安地区2004-2009年杉木林固定样地51个。所有样地数据包含样地编号、样地位置、行政区、地理地貌信息、土壤信息和样地内植被信息等。样木数据包含样地编号、样木编号、树种、胸径等。由于样地是固定（永久）的，样地号和样木号在不同调查期是固定的，因而可实现每株样木生长的连续监测。
-
数量化理论Ⅰ是一种类似多元回归的分析方法[15-16]，用于自变量都是定性变量、基准变量是定量变量的因素分析与预测问题，采用说明性多变量模拟线性表示式中基准变量的定量变化。假定有如下的线性模型：
$$ y_{i}=\sum\limits_{u=1}^{h} b_{u} x_{i}(u)+\sum\limits_{j=1}^{m} \sum\limits_{k=1}^{r_j} \delta_{i}(j, k) b_{j k}+\varepsilon_{i}, i=1, 2, \cdots, n。 $$ (1) 假设自变量中有h个是定量的，它们在第i个样本中的数据为xi（u）（u=1, 2, …, n），有m个是定性的，即m个项目，其中第j项目有rj个类目，它们在第i个样本中的反应式基准变量的数据为yi（i=1, 2, …, n）。
利用最小二乘法，即Bu（u=1, 2, …, h）和bjk（j=1, 2, …, m; k=1, 2, …, rj）的最小二乘估计构造正规方程：
$$ X^{\prime} X \hat{b}=X^{\prime} y。 $$ (2) 求解正规方程组（2）得预测方程, 并对方程精度进行检验，同时通过偏相关系数、方差比和范围评价每个立地因子对因变量作用的大小。
数量化理论Ⅰ相对于一般回归分析，其不同之处在于可把定性变量与定量变量同时纳于回归模型。在数量化理论模型中，因变量和自变量被分别称为基准变量和说明变量。说明变量中定性变量称为项目，每个定性变量的划分称为类目，定性变量的类目类似于定量变量的取值区间。
-
根据数据易采集、便于立地质量评价模型推广的原则，参考前人研究成果，选择立地因子。本研究选取影响杉木生长的地形因子与土壤类型作为立地因子。地形因子包括地貌、纬度、海拔、坡度、坡向、坡位等，其中纬度、海拔、坡度是定量变量，地貌、坡向、坡位与土壤类型是定性变量。根据数量化理论Ⅰ对变量的约定，结合样地立地因子数据，对定性立地因子进行类目划分（表 1），并对项目（因子）的类目利用示性函数（0, 1）数量化方法对定性因子进行定量化处理，即:
表 1 定性立地因子类目划分与量化编号
Table 1. Classification and quantitative numbering of qualitative site factors
立地因子 类目 划分标准 编号 地貌 中山 海拔为1 000 m以上的山地 1 低山 海拔500~999 m的山地 2 丘陵 海拔＜500 m，相对高差100 m以下，没有明显的脉络 3 坡位 山脊 山脉的分水线及其两侧各下降垂直高度15 m的范围 1 上坡 从脊部以下至山谷范围内的山坡3等分后的最上等分部位 2 中坡 从脊部以下至山谷范围内的山坡3等分后的中坡位 3 下坡 从脊部以下至山谷范围内的山坡3等分后的最下等分部位 4 山谷 汇水线两侧的谷地 5 平地 平台或台地的地段 6 坡向 东 方位角68°~112° 1 南 方位角158°~202° 2 西 方位角248°~292° 3 北 方位角338°~23° 4 东北 方位角23°~67° 5 东南 方位角113°~157° 6 西南 方位角203°~247° 7 西北 方位角293°~337° 8 无坡向 坡度小于5° 9 土壤类型 红壤 按实际土壤类型划分（数据来源为浙江省土壤类型调查数据） 1 黄壤 2 水稻土 3 紫色土 4 其他 5 $$ \delta_{i}(j, k)=\left\{\begin{array}{l} {1}&当i样本第j项目取第k类目时\\ {0}&当i样本第j项目不取第k类目时 \end{array}\right. 。 $$ (3) 将每个样地的定性因子代入式（3）中，与定量因子一起编制立地因子数量化反应表（表 2），用于立地质量模型建立与评价。
表 2 立地因子反应表
Table 2. Response table of site factors
样地号 坡位（X1） … 土壤类型（X3） 纬度 海拔
(协变量)坡度 山脊 上坡 中坡 下坡 山谷 平地 … 红壤 黄壤 水稻土 紫色土 其他 22570 0 0 1 0 0 0 … 1 0 0 0 0 31 500 40 22599 0 0 0 1 0 0 … 0 0 1 0 0 31 540 30 22602 0 0 1 0 0 0 … 1 0 0 0 0 31 480 35 … … … … … … … … … … … … … … … … 22549 0 0 1 0 0 0 … 0 0 1 0 0 31 850 30 说明：纬度用公里网（百公里）表示 -
反映森林立地质量最主要和直接的指标是生产力。可观测的森林立地生产力指标主要有特定年龄的树高、胸径和蓄积量，或者其生长量，应用最多的是特定年龄的优势木树高。然而，优势木树高和年龄的准确测量不仅成本很高，而且难度特大。胸径生长率是最能直观反映立地质量的森林立地质量评价因子。现有的NFI数据信息量丰富，且数据可靠性好。胸径是NFI调查中最为重要的测树因子，但在NFI中不调查单株样木的年龄。从NFI中提取的胸径生长率只能代表在某一清查期内某棵树胸径的年均生长率，并不能直接反映所在地段的立地质量好坏。为了充分利用NFI数据开展森林立地质量评价，本研究提出1个新的立地质量评价指标：最大胸径生长率指标。以最大胸径生长率指标为立地评价指标，与地形地貌、环境、土壤等易获取因子建立模型，可以用于估计包含无林地在内的所有林地的潜在生产力。一般来说，林木的生长规律是从幼龄林到成熟、过熟林，胸径和树高的生长速度都呈逐渐降低的趋势，即在树木的生长过程中生长率都会出现一个高峰。因此，在NFI清查次数足够多的情况下, 每个样地复位样木的生长率都会在某一时期出现最高值，即最大胸径生长率，它代表了该地块的最大生产潜力。若以R表示最大胸径生长率，则：
$$ R=\max \left\{r_{k}, k=1, 2, \cdots, n\right\}。 $$ (4) 式（4）中：k为清查期，rk为第k期的胸径生长率，n为森林资源清查总期数。
NFI清查间隔期为5 a，杉木的龄级期也是5 a，5期NFI清查数据即可覆盖5个龄级期的生长率。假定杉木轮伐期为25 a，近熟林进入成熟林时全部采伐，则5次清查能包含全部最大胸径生长率。因此，NFI与龄级之间的关系如图 1所示。图 1中，样地样木的龄级为Ⅰ~Ⅴ，成熟林全部采伐（由于现实森林的林龄不同，因而有A、B、C、D、E等5种情形）。
图 1 清查期次与样木龄级分布
Figure 1. Inventory period and age class distribution of sample trees
样木胸径生长率的计算方法如下：
$$ r_{i j k}=\frac{D_{i j(k+1)}-D_{i j k}}{5 D_{i j k}} 。 $$ (5) 式（5）中：rijk为第i个样地、第j株复位样木、第k期清查的胸径生长率；Dij（k+1）为k+1期样木胸径；Dijk为第k期样木胸径。
样地胸径生长率计算建立在样木生长率的基础上，设第i个样地、第k期清查的胸径生长率为rik，则：
$$ r_{i k}=\frac{1}{n} \sum\limits_{j=1}^{j=n}\left\{r_{i j k}\right\}。 $$ (6) 样地优势木生长率比样地生长率更能代表样地所在位置的立地生产潜力，样地优势木胸径生长率计算公式如下：
$$ r_{ik\mathrm{D}}=\frac{1}{3} \sum\limits_{j=1}^{j=3}\left\{r_{{i} j k}\right\}。 $$ (7) 式（7）中：rikD即第i样地第k期3株优势木生长率的均值。
根据样地的胸径生长率，从每个样地多期调查数据所生成的胸径生长率时间系列中选取最大者：
$$ r_{i\text {D }}=\max \left\{r_{ik\text {D } }, k=1, 2, \cdots, p\right\}。 $$ (8) 式（8）中：riD相当于样地i位置幼龄林优势木的平均胸径生长率。
-
复相关系数是估计回归精度的重要指标之一，复相关系数越大，则估测相关越密切，说明模型对优势木最大胸径生长率均值的估计效果越好，所得的数量化得分表就越可靠[21]。根据复相关系数检验线性关系是否显著进行F检验，可确定坡位、坡向、土壤类型、纬度位置、海拔、坡度等变量的线性关系是否显著[22]。
-
模型预测精度检验采用均方根误差（RMSE）作为模型的评价指标。均方根误差是用来衡量实际测量值同模型预测值之间的偏差，它可作为衡量测量精度的一种数值指标[23]。在模型建立后，为对模型拟合性能进行检验，运用浙江省杭州市临安区的实际测量值代入模型计算最大胸径生长率进行实地检验，采用平均绝对误差、RMSE值进行评估，反映模型估测的精度。
-
采用方差比方法进行。该检验方法具有精确的样本分布，而不依赖于大样本渐近极限分布；在通常的时间序列数据非正态分布的情况下，这种非参数方差比检验具有更高的检验功效。方差比检验对于一段时间内某一时间间隔的方差比例效果较好。方差比检验是指各因子（项目）方差占模型方差的百分比，表示立地因子（坡位、坡向、土壤类型、纬度位置、海拔、坡度）对因变量（最大胸径生长率）的作用是否显著。
-
根据NFI数据分析，杉木生长率从幼林到成熟林呈下降趋势。图 2A和图 2B显示了1994-1999年和2004-2009年调查间隔期内杉木单株生长率按胸径分布的变化。首先，胸径生长率随胸径增大而降低，2个时期趋势相同。其次，其趋势为单调下降。再次，曲线在纵轴上分布较宽，特别是当胸径大约为20~30 cm时，生长率出现较大的变化范围。图 2B（2004-2009年）最为明显。它们代表了多条单调下降的曲线，其本质是不同立地质量的生长率差异。这一特点图 2B较图 2A突出，原因是2004年这有部分位置移动后的新样地加入。2幅图同时显示，随着胸径增大，其生长率分化明显，而在较小胸径阶段生长率随胸径增大而较快下降，且在纵轴上分布区间较窄。
图 2 杉木单株生长率比较
Figure 2. Comparison of growth rate of C. lanceolata
-
单株胸径生长率连续变化趋势与年龄有关。在无年龄的情况下，这种变化在连续多期之间呈单调曲线形式上升或下降。图 3为单株木胸径生长率随胸径增长而的变化情况。为了图示清晰，将每个胸径值所对应的多株样木胸径生长率取均值，每条曲线代表 1个胸径级。图 3中D8代表第1次清查时胸径为8 cm，D12代表第1次清查时胸径为12 cm，……。余类推。从图 3可以发现：小径级样木胸径生长率大，大径级样木胸径生长率小。以后一直保持这一规律。对于同一株样木来说，清查次数增加就是年龄的增加，随着年龄的增长，样木胸径生长率下降，最后趋于平缓。结果表明：在杉木生长过程中，单株胸径生长率存在最大值，这个最大值代表所在地段的森林立地生产力。因此，最大胸径生长率可以被用作森林立地质量评价的直接指标。
图 3 杉木不同胸径径级平均生长率随时间变化趋势
Figure 3. Trend of average growth rate of different DBH grades of C. lanceolata with time
-
通过多期复位样木分析，提取每个样地的最大胸径生长率。由于样地采伐、灾害等因素影响，用样地全部样木胸径生长率的平均值代表样地生产潜力的状态不适合所有样地。表 3显示：用3株最大胸径生长率取平均与样地的立地因子有更好的相关性。因此，为了避免偶然因素的干扰，每个样地的最大胸径生长率取3棵优势木胸径生长率的均值。共提取529个样地的最大胸径生长率。
表 3 最大胸径生长率与立地因子之间皮尔森相关系数
Table 3. Pearson correlation coefficient between the maximum DBH growth rate and the site factor
纬度 海拔 地貌 坡向 坡位 坡度 土壤类型 最大胸径生长率（3株） -0.420 4 0.558 3 -0.447 8 -0.216 3 -0.114 4 0.068 6 0.145 2 最大胸径生长率（全部） -0.061 5 -0.059 6 0.025 7 -0.019 6 -0.018 2 0.071 7 -0.009 2 -
以优势木最大胸径生长率均值作为因变量，坡位、坡向、土壤类型、纬度位置、海拔、坡度作为自变量，采用数量化理论Ⅰ建立杉木林立地质量评价模型。4期529个NFI的样地数据参与数量化拟合，运用软件SPSS 19.0完成线性回归分析模型建立与评价。具体结果见表 4。
表 4 杉木林立地因子数量化理论Ⅰ得分及t检验表
Table 4. Score and t-test table of the quantification theory Ⅰ model of C. lanceolata
项目 类目 得分值 偏相关系数 t检验 截距 30.971 6.118 坡位 山脊 0.035 0.229 5.185 上坡 -0.186 中坡 0.000 下坡 0.076 山谷 2.737 平地 -2.064 坡向 东 -0.241 0.437 7.887 南 0.000 西 -0.908 北 -0.289 东北 -0.232 东南 -0.749 西南 0.171 西北 0.120 无坡向 1.781 土壤 红壤 0.000 0.147 3.388 黄壤 0.076 水稻土 -0.597 紫色土 -0.461 其他 -1.391 纬度位置 -0.797 0.137 3.155 海拔/m 0.006 -0.038 -0.801 坡度 -0.011 0.080 1.835 复相关系数R=0.606 F=14.723 根据表 4数据，本研究最终建立浙江省杉木林立地质量评价预测模型公式为：
$$ Y= 30.971-0.797δ（1）+0.006δ（2）-0.011δ（3）+\\0.035δ（4, 1）-0.186δ（4, 2）+0.076δ（4, 4）+2.737δ（4, 5）-\\2.064δ（4, 6）-0.241δ（5, 1）-0.908δ（5, 3）\\-0.289δ（5, 4）-0.232δ（5, 5）-0.749δ（5, 6）+0.171δ（5, 7）+\\ 0.120δ（5, 8）+1.781δ（5, 9）+0.076δ（6, 2）-0.597δ（6, 3）\\-0.461δ（6, 4）-1.391δ（6, 5）。 $$ (9) 式（9）中：δ（1）为样地纬度位置，δ（2）为样地海拔，δ（3）为样地坡度，δ（4, i）为样地坡位第i个项目，δ（5, i）为样地坡向第i个项目，δ（6, i）为样地土壤类型第i个项目。表 4显示模型统计检验F=14.723，达极显著水平，模型复相关系数R为0.606。方差分析（表 5）表明：各项目对模型贡献均达显著水平。
表 5 各类立地因子影响因素方差比
Table 5. Variance ratio of the influencing factors of various site factors
第j类因素 1（坡位） 2（坡向） 3（土壤类型） 4（纬度） 5（海拔） 6（坡度） nσy2 方差（nσy2） 265.498 3 341.505 3 152.350 1 64.841 0 20.612 2 84.161 7 928.968 7 方差比$\left(\frac{\sigma_{j}^{2}}{\sigma_{y}^{2}}\right)$ 0.285 8 0.367 6 0.164 0 0.069 8 0.022 2 0.090 6 -
方差比是各因子（项目）方差占模型方差的百分比，表示因子对因变量（最大胸径生长率）的作用是否显著。按公式计算各影响因素的方差，结果见表 5。结果表明：对胸径生长率的贡献程度从大到小依次是坡向、坡位、土壤类型、坡度、纬度、海拔。其中坡向对杉木林胸径生长率的贡献率最大，坡向会影响太阳光照量和太阳辐射强度，造成不同坡向的日照强度和热量不同，从而对杉木林的生长产生较大影响。坡位对杉木林的胸径生长率的贡献率较大，不同坡位具有相应的小气候环境，影响着湿度、水分和温度等条件。土壤类型对胸径生长率的贡献同样较大的原因是由于土壤类型不同，杉木生长的土壤营养环境不同，从而对杉木林的生长产生一定影响。从表 4和表 5可以看出：纬度和海拔的影响较小，主要是因为本研究所选择的研究区为浙江省，纬度变化范围较小；同时浙江省内山地平均海拔一般在200~1 000 m，而本研究杉木人工林样地分布在400~800 m较窄的范围内，两者的差距对于杉木林来说不能形成较大区分度，因此纬度和海拔的对于胸径生长率的贡献率较小。
-
在模型建立后，为对模型拟合性能进行检验，运用浙江省杭州市临安区的实际测量值代入模型计算最大胸径生长率进行实地检验，采用平均绝对误差、RMSE值进行评估，对全部样本和临安地区样本的观测值和预测值之间的误差进行评价，反映模型估测的精度。从表 6可以看出：临安地区样本和全部样地作为模型估计结果检验的平均绝对误差分别为3.407和2.125，整体平均绝对误差相对较小，而且随着样本的增加，模型的拟合精度提高，表明利用最大生长率进行数量化理论拟合是精度相对较高的方法。均方根误差能够反映预测值和测量值之间的偏差程度，数据显示临安区样本和全部样本检验的均方根误差为0.518和1.421，而计算两者的相对均方根误差的值分别为0.073和0.062，相对均方根误差越小，相对均值的离散程度越小，模型的精度也就越高。
表 6 模型拟合误差分析表
Table 6. Error analysis table of the model-fitting
拟合统计量 平均绝对误差 均方根误差 相对均方根误差 临安地区样本 3.407 0.518 0.073 全部样本 2.125 1.421 0.062 -
基于NFI的最大胸径生长率是杉木人工林立地质量评价的可靠指标，在经济上节约、高效。从建模因子的贡献率来看，对胸径生长率的贡献程度从大到小依次是坡向、坡位、土壤类型、坡度、纬度、海拔。利用最大胸径生长率，结合DEM数据建立数量化理论Ⅰ模型，模型拟合效果和检验结果均比较理想。在缺乏树龄和树高数据的情况下，可以利用基于多期NFI数据的样地最大胸径生长率作为评价指标进行森林立地质量评价。不仅可以用于现有杉木人工林立地质量评价，还可以用于无林地立地生产潜力预测。
本研究提出基于NFI的最大胸径生长率作为杉木人工林立地质量评价指标，结合地理信息建立立地质量评价模型，结果符合研究区内杉木林生长立地条件现状。但限于研究区NFI数据积累的时间跨度有限，模型精度还有进一步提高的空间。由于NFI每5 a清查1次，其数据的时间分辨率为5 a。因此，NFI清查期数越多，在同一个样地所测样木的年龄阶段越多，从中提取的最大胸径生长率对相应立地生产力的代表性越接近真实。随着NFI数据的不断积累，基于NFI的森林立地质量评价精度会不断提高。同时，基于数量化理论Ⅰ的杉木最大胸径生长率估计模型所解决的问题是得到所有森林立地上的假设最大胸径生长率，只是一个以生长率代表立地质量的初步展现，在此基础上还可以有进一步的分类、聚类及其检验等分析。
Site quality evaluation of Cunninghamia lanceolata plantations in Zhejiang Province based on maximum DBH growth rate
-
摘要:
目的 森林立地质量评价是指对立地的宜林性或潜在的生产力进行判断和预测，目的在于建立浙江省杉木Cunninghamia lanceolata人工林立地质量评价数量化理论模型。 方法 提出1个新的立地质量评价指标：最大胸径生长率指标。以最大胸径生长率指标为立地评价指标，以浙江省1994-2009年国家森林资源清查（NFI）数据为基础，从NFI复位样木时间序列数据中提取529个杉木人工林样地优势木最大胸径生长率；结合易于获取的纬度、海拔、坡位、坡向、坡度、土壤类型等立地因子，基于数量化理论Ⅰ方法，建立杉木人工林立地质量评价模型。 结果 模型统计检验值为14.723，模型复相关系数为0.606，达极显著水平，模型效果较理想。基于NFI的最大胸径生长率是立地质量评价的合适指标；坡向、坡位、土壤类型、坡度、纬度、海拔等立地因子对最大胸径生长率贡献依次减弱。 结论 基于NFI和数量化理论Ⅰ的浙江省杉木人工林立地质量评价模型总体指标合理，克服了传统立地质量评价对年龄的依赖，能有效利用NFI历史数据，适用于同龄林和异龄林；构建以最大胸径生长率为指标的立地质量评价模型，能对研究区立地质量进行有效评价，不仅可对现有杉木人工林生产力进行评价，也可以对无林地的生产潜力进行预测。 Abstract:Objective Forest site quality evaluation refers to the judgement and prediction of site suitability or potential productivity. The aim is to establish a site quality evaluation model of Cunninghamia lanceolata plantaion. Method A new site quality evaluation index, the maximum growth rate index of DBH, was proposed. Taken as the site evaluation index, and based on the National Forest Inventory(NFI) data of Zhejiang Province from 1994 to 2009, the maximum DBH growth rate for 529 permenant plots of C. lanceolata plantations was extracted from the NFI time series data as the site quality evaluation index of C. lanceolata plantations. Integrated with readily available site factors such as latitude, altitude, slope position, slope direction, slope, soil types, etc., and based on the quantitative theory Ⅰ method, a site quality evaluation model of C. lanceolata plantations was established. Result The statistical test value of the model is 14.723, and the complex correlation coefficient of the model is 0.606, which reaches a very significant level, and the effect of the model is satisfactory. The maximum DBH growth rate based on NFI is an appropriate index for site quality evaluation. The contribution of site factors such as slope direction, slope position, soil type, slope, latitude and altitude to the maximum growth rate of DBH decreases sequentially. Conclusion Based on NFI and quantitative theory Ⅰ, the overall index of site quality evaluation model of C. lanceolata plantations in Zhejiang Province is reasonable. It overcomes the dependence of traditional site quality evaluation on age, and can effectively use NFI historical data, which is suitable for forests of the same age and different age. The site quality evaluation model can effectively evaluate the site quality of the study area. It can not only evaluate the productivity of existing C. lanceolata plantations, but also predict the potential productivity of non-forest land. -
图 3 杉木不同胸径径级平均生长率随时间变化趋势
Figure 3 Trend of average growth rate of different DBH grades of C. lanceolata with time
表 1 定性立地因子类目划分与量化编号
Table 1. Classification and quantitative numbering of qualitative site factors
立地因子 类目 划分标准 编号 地貌 中山 海拔为1 000 m以上的山地 1 低山 海拔500~999 m的山地 2 丘陵 海拔＜500 m，相对高差100 m以下，没有明显的脉络 3 坡位 山脊 山脉的分水线及其两侧各下降垂直高度15 m的范围 1 上坡 从脊部以下至山谷范围内的山坡3等分后的最上等分部位 2 中坡 从脊部以下至山谷范围内的山坡3等分后的中坡位 3 下坡 从脊部以下至山谷范围内的山坡3等分后的最下等分部位 4 山谷 汇水线两侧的谷地 5 平地 平台或台地的地段 6 坡向 东 方位角68°~112° 1 南 方位角158°~202° 2 西 方位角248°~292° 3 北 方位角338°~23° 4 东北 方位角23°~67° 5 东南 方位角113°~157° 6 西南 方位角203°~247° 7 西北 方位角293°~337° 8 无坡向 坡度小于5° 9 土壤类型 红壤 按实际土壤类型划分（数据来源为浙江省土壤类型调查数据） 1 黄壤 2 水稻土 3 紫色土 4 其他 5 
下载: 导出CSV
表 2 立地因子反应表
Table 2. Response table of site factors
样地号 坡位（X1） … 土壤类型（X3） 纬度 海拔
(协变量)坡度 山脊 上坡 中坡 下坡 山谷 平地 … 红壤 黄壤 水稻土 紫色土 其他 22570 0 0 1 0 0 0 … 1 0 0 0 0 31 500 40 22599 0 0 0 1 0 0 … 0 0 1 0 0 31 540 30 22602 0 0 1 0 0 0 … 1 0 0 0 0 31 480 35 … … … … … … … … … … … … … … … … 22549 0 0 1 0 0 0 … 0 0 1 0 0 31 850 30 说明：纬度用公里网（百公里）表示
下载: 导出CSV
表 3 最大胸径生长率与立地因子之间皮尔森相关系数
Table 3. Pearson correlation coefficient between the maximum DBH growth rate and the site factor
纬度 海拔 地貌 坡向 坡位 坡度 土壤类型 最大胸径生长率（3株） -0.420 4 0.558 3 -0.447 8 -0.216 3 -0.114 4 0.068 6 0.145 2 最大胸径生长率（全部） -0.061 5 -0.059 6 0.025 7 -0.019 6 -0.018 2 0.071 7 -0.009 2
下载: 导出CSV
表 4 杉木林立地因子数量化理论Ⅰ得分及t检验表
Table 4. Score and t-test table of the quantification theory Ⅰ model of C. lanceolata
项目 类目 得分值 偏相关系数 t检验 截距 30.971 6.118 坡位 山脊 0.035 0.229 5.185 上坡 -0.186 中坡 0.000 下坡 0.076 山谷 2.737 平地 -2.064 坡向 东 -0.241 0.437 7.887 南 0.000 西 -0.908 北 -0.289 东北 -0.232 东南 -0.749 西南 0.171 西北 0.120 无坡向 1.781 土壤 红壤 0.000 0.147 3.388 黄壤 0.076 水稻土 -0.597 紫色土 -0.461 其他 -1.391 纬度位置 -0.797 0.137 3.155 海拔/m 0.006 -0.038 -0.801 坡度 -0.011 0.080 1.835 复相关系数R=0.606 F=14.723
下载: 导出CSV
表 5 各类立地因子影响因素方差比
Table 5. Variance ratio of the influencing factors of various site factors
第j类因素 1（坡位） 2（坡向） 3（土壤类型） 4（纬度） 5（海拔） 6（坡度） nσy2 方差（nσy2） 265.498 3 341.505 3 152.350 1 64.841 0 20.612 2 84.161 7 928.968 7 方差比$\left(\frac{\sigma_{j}^{2}}{\sigma_{y}^{2}}\right)$ 0.285 8 0.367 6 0.164 0 0.069 8 0.022 2 0.090 6
下载: 导出CSV
表 6 模型拟合误差分析表
Table 6. Error analysis table of the model-fitting
拟合统计量 平均绝对误差 均方根误差 相对均方根误差 临安地区样本 3.407 0.518 0.073 全部样本 2.125 1.421 0.062
下载: 导出CSV
-
[1] PRETZSCH H. Forest Dynamics, Growth and Yield: from Measurement to Model [M]. Berlin and Heidelberg: Springer-Verlag, 2009. [2] DHÔTE J F, HERVÉ J C. Productivity changes in four sessile oak forests since 1930: a stand-level approach[J]. Ann For Sci, 2000, 57(7): 651 - 680. [3] PHILLIPS O L, MALHI Y, BARBARA V, et al. Changes in growth of tropical forests: evaluating potential biases[J]. Ecol Appl, 2002, 12(2): 576 - 587. [4] VANCLAY J K, SKOVSGAARD J P, HANSEN C P. Assessing the quality of permanent sample plot databases for growth modelling in forest plantations[J]. For Ecol Manage, 1995, 71(3): 177 - 186. [5] 张雅梅, 何瑞珍, 安裕伦, 等.基于RS与GIS的森林立地分类研究[J].西北林学院学报, 2005, 20(4):147 - 152. ZHANG Yamei, HE Ruizhen, AN Yulun, et al. Forest site type classification and site quality evaluation based on technologies of RS and GIS[J]. J Northwest For Univ, 2005, 20(4):147 - 152. [6] 王永昌, 张金池.基于遥感技术的云台山立地分类及质量评价[J].南京林业大学学报(自然科学版), 2007, 31(1):85 - 89. WANG Yongchang, ZHANG Jinchi. The classification and evaluation of forest-sites in Yuntaishan Mountain area based on remote sensing technology[J]. J Nanjing For Univ Nat Sci Ed, 2007, 31(1):85 - 89. [7] 季碧勇, 陶吉兴, 王文武.基于连续清查固定样地生物量的立地质量评价[J].西南林业大学学报, 2012, 32(4):45 - 50. JI Biyong, TAO Jixing, WANG Wenwu. Forest site quality evaluation based on biomass data from continuous forest inventories[J]. J Southwest For Univ, 2012, 32(4):45 - 50. [8] 唐诚, 王春胜, 庞圣江, 等.广西大青山西南桦人工林立地指数表编制[J].浙江农林大学学报, 2019, 36(4):828 - 834. TANG Cheng, WANG Chunsheng, PANG Shengjiang. et al. Site index table construction for Betula alnoioles plantations in the Daqing Mountains, Guangxi[J]. J Zhejiang A&F Univ, 2019, 36(4):828 - 834. [9] 浦瑞良, 杨金中, 刘毓起, 等.紫金山彩红外片马尾松数量化立地指数表的编制[J].浙江林学院学报, 1994, 11(1):64 - 68. PU Ruiliang, YANG Jinzhong, LIU Yuqi, et al. Tabulation of quantifitation color infrared aerial photo site-index for Pinus massonianain Mt. Zijin[J]. J Zhejiang For Coll, 1994, 11(1):64 - 68. [10] 刘财富, 常新东, 庞家清, 等.刺槐人工林数量化立地指数表编制的研究[J].辽宁林业科技, 1997(6): 37 - 40. LIU Caifu, CHANG Xindong, PANG Jiaqing, et al. Study on the compilation of quantitative site index table of Robinia pseudoacacia plantation[J]. J Liaoning For Sci Technol, 1997(6): 37 - 40. [11] 汪笑安.旺业甸实验林场立地分类与质量评价研究[D].北京: 北京林业大学, 2013. WANG Xiaoan. The Classification and Evaluation of Site Conditions in Wangyedian Forest Farm[D]. Beijing: Beijing Forestry University, 2013. [12] 赖挺.四川巨按人工林立地分类研究[D].成都: 四川农业大学, 2005. LAI Ting. Study on Classification Eucalyptus grandis Platation in Sichuan[D]. Chengdu: Sichuan Agricultural University, 2005. [13] 杜健, 梁坤南, 周再知, 等.云南西双版纳柚木人工林立地类型划分及评价[J].林业科学, 2016, 52(9):1 - 10. DU Jian, LIANG Kunnan, ZHOU Zaizhi, et al. Site classification and evaluation of teak plantation in Xishuangbanna, Yunnan Province, China[J]. Sci Silv Sin, 2016, 52(9):1 - 10. [14] 倪志云, 陆贵巧, 张锁成, 等.河北省山区核桃立地质量评价研究[J].河北林果研究, 2006, 21(1):53 - 56. NI Zhiyun, LU Guiqiao, ZHANG Suocheng, et al. Study on site quality evaluation of walnut in the mountainous area of Hebei Province[J]. Hebei J For Orchard Res, 2006, 21(1):53 - 56. [15] 黄永平, 田小海.数量化理论Ⅰ在土地分级中的应用[J].湖北农学院学报, 2000, 20(3):236 - 240. HUANG Yongping, TIAN Xiaohai. Application of the quantification theory type Ⅰto land classification[J]. J Hubei Agric Coll, 2000, 20(3):236 - 240. [16] 于丽, 王明年, 房敦敏, 等.岩质围岩施工阶段亚级分级的数量化理论研究[J].岩土力学, 2009, 30(12):3846 - 3850. YU Li, WANG Mingnian, FANG Dunmin, et al. Study of quantification theory of rocky surrounding rock sub-classification during construction[J]. Rock Soil Mech, 2009, 30(12):3846 - 3850. [17] 秦利燕, 秦玉权, 邵春福, 等.应用改进的数量化理论预测道路断面事故率[J].中国安全科学学报, 2011, 21(11):87 - 91. QIN Liyan, QIN Yuquan, SHAO Chunfu, et al. Application of improved quantitative theory in prediction of road section accident rate[J]. China Saf Sci J, 2011, 21(11):87 - 91. [18] 孟宪宇.测树学[M]. 3版.北京:中国林业出版社, 2006. [19] 林狄显.杉木人工林立地质量评价方法的研究[J].内蒙古林业调查设计, 2011, 34(2):41 - 42. LIN Dixian. Study on site quality evaluation methods of Chinese fir plantation[J]. Inner Mongolia For Invest Des, 2011, 34(2):41 - 42. [20] 徐俊森.流沙海岸风口地段营建防护林的配套技术[J].浙江林学院学报, 2005, 22(2):139 - 143. XU Junsen. Techniques for constructing protection forest on the quicksand in drought costal area[J]. J Zhejiang For Coll, 2005, 22(2):139 - 143. [21] 郭明玲, 赵克昌.河西走廊农田防护林立地质量评价[J].防护林科技, 2004(5):5 - 7. GUO Mingling, ZHAO Kechang. Evaluation of site quality of farmland shelterbelt in Hexi corridor[J]. Prot For Sci Technol, 2004(5):5 - 7. [22] 张康健, 王蓝, 孙长忠.森林立地定量评价与分类[M].西安:陕西科学技术出版社, 1988:25 - 26. [23] 贾俊平.统计学基础[M].北京:中国人民大学出版社, 2010. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2020.01.014
点击查看大图
图(3) / 表(6)
计量
- 文章访问数: 1823
- HTML全文浏览量: 543
- PDF下载量: 38
- 被引次数: 0
