立地因子 数量化值 1 2 3 4 坡向 阳坡 阴坡 半阴半阳 坡位 上坡 中坡 下坡 谷地
3结果与分析 3.1 不同立地条件下毛竹生长状况分析 通过40个毛竹标准地的调查,可以发现不同的立地条件对毛竹生长具有一定的影响,但不同的立地因子影响的程度不一样。每个因子对毛竹生长影响的具体情况是:半阴半阳坡的生长较好,阴坡生长较差;坡度较小的生长较好,坡度太大的生长较差;中下坡的生长较好,上坡的生长较差;腐质层偏厚的生长较好,腐质层偏薄的生长较差;郁闭度适中的生长较好,郁闭度过大或过小的生长较差;而在调查范围内,立株度对毛竹生长的影响不明显,但立株度偏小的毛竹胸径大一些,立株度偏大的毛竹胸径小一些。 3.2逐步回归分析 用坡向(X1)、坡度(X2)、坡位(X3)、腐质层厚度(X4)、郁闭度(X5)、立株度(X6)6个立地因子对林分胸径生长量(y)进行逐步回归,回归方程为y 5.80+0.03X2+0.42X3+0.13X4,F=7.55**回归方程达到极显著水平,因此可以用该方程对毛竹林分的胸径生长量进行预测。 对回归方程的截距和每个自变量进行偏回归检验,结果见表2。表2 偏回归的假设检验变量 回归系数 标准误 F Pr>F 截距 5.79556 0.77597 55.78 <0.0001 X2 0.02865 0.00816 12.32 0.0012 X3 0.41537 0.18023 5.31 0.0271 X4 0.12844 0.04629 7.70 0.0087 复相关系数R=0.621 4**(r0.05=0.440 r0.01=0.517)
从表2可知:截距、X2与X4的偏回归都达到极显著,而X3的偏回归也达到显著水平。这也说明,坡度、坡位和腐质层厚度分别对毛竹生长有显著影响,它们即为影响毛竹生长的主导因子,且影响大小顺序依次为:坡度>腐质层厚度>坡位。而另外3个立地因子对毛竹生长的影响不大。 3.3 立地因子间相关分析 计算6个立地因子的平均数、标准差及变异系数,结果见表3。表3 各立地因子的平均数、标准差及变异系数名称 X1 X2 X3 X4 X5 X6 平均数 1.090 2.750 20.050 7.200 28.450 5.550 标准差 0.884 2.468 13.077 14.133 26.668 6.677 变异系数 0.811 0.897 0.652 1.963 0.937 1.203
从表3可知:在这40块标准地中,腐质层厚度的差异程度最大,其次是立株度,而坡位的变异程度最小,其次是坡向。 计算6个立地因子两两间的相关系数,结果见表4。表4 立地因子间的相关系数表X1 X2 X3 X4 X5 X6 X1 1.000 X2 0.550** 1.000 X3 -0.752** -0.608** 1.000 X4 0.518** 0.266 -0.458** 1.000 X5 -0.543** -0.500** 0.718** -0.326** 1.000 X6 0.745** 0.887** -0.707** 0.489** -0.505** 1.000
注:r0.05=0.312 r0.01=0.403。由表4可以看出:只有腐质层厚度与坡度的相关性是不显著的,说明腐质层厚度与坡度2个因子之间的相互作用较小,它们能分别对毛竹胸径生长量产生影响。另外,郁闭度与腐质层厚度两因子之间的相关性达到显著,而其它两两因子之间的相关性达到极显著。 3.4主成分分析 在表3与表4的基础上,对6个立地因子进行主成分分析,结果见表5。表5 各主成分的特征根、特征向量、贡献率和累积贡献率主成分特征向量因子 Y1 Y2 Y3 Y4 Y5 Y6 X1 0.434 0.165 0.009 0.728 0.450 0.227 X2 0.408 -0.497 0.379 -0.320 0.004 0.584 X3 -0.449 -0.003 0.319 -0.225 0.804 -0.025 X4 0.304 0.826 0.229 -0.395 -0.025 0.127 X5 -0.375 0.101 0.730 0.399 -0.384 0.095 X6 0.460 -0.181 0.411 -0.036 0.059 -0.763 特征根 3.914 0.798 0.661 0.375 0.195 0.058 贡献率 0.652 0.133 0.110 0.063 0.033 0.010 累积贡献率 0.652 0.785 0.895 0.958 0.990 1.000
从表5可知:前2个主成分的累积贡献率已达到78.5%,本文选取前2个主成分。第1主成分的贡献率为65.2%,即该成分已表达了总信息的65.2%,而在第1主成分的特征向量中,X6和X3的绝对值最大,说明立竹度和坡位对Y1的影响最大,再根据表4可知:立竹度与坡位的相关系数达-0.707,可以认为立竹度是随坡位变化而变化的1个因子,因此,有理由把Y1看作是反映坡位的因子;第2主成分的贡献率为13.3%,可表达总信息量的13.3%,而在第2主成分的特征向量中,X4和 X2的绝对值最大,且两者的相关系数只有0.266,说明腐质层厚度与坡度对Y2的影响都很大,有理由把Y2看作是反映腐质层厚度与坡度的因子。 3.5 主成分聚类分析 在主成分分析结果的基础上,本文只取累积贡献率已达到78.5%的前2个主成分作为立地因子主成分模型,即: Y1=0.434X1+0.408X2-0.449X3+0.304 X4-0.375X5+0.460X6 Y2=0.165X1-0.497X2-0.003X3+0.826X4+0.101X5-0.181X6 将各立地因子指标数据分别代入上面2个公式,即可求出以Y1和Y2为坐标轴的坐标值,以Y1为纵轴,Y2为横轴,可作标准地的散点分布图,见图1。每一类中的标准地在分布图中都相距较近,说明它们具有相似的立地条件,并结合各标准地资料可知,其中第工类的特点是:大多位于谷地,处于阳坡,且坡度很小,一般在5°以下,腐质层厚度较厚,接近20 cm;第Ⅱ类的特点是:大多处于半阴半阳坡,坡度适中,一般在15°左右,腐质层较第1类薄,在12 cm左右;第Ⅲ类的特点是:大多也处于半阴半阳坡,但坡度较大,一般在25°以上,且腐质层较薄,在8 cm左右。根据各类不同的特点,可以采用相应的经营管理措施。但每一类中的标准地,可以采用相似的经营方式和管理措施。而不包含在这三大类中的标准地,因各具特点,在今后的生产中因地制宜即可。 4 小结 (1)通过采用数量化的方法,先将坡向与坡位2个立地因子数量化,然后与另外4个立地因子一并进行逐步回归分析,筛选得出影响毛竹生长的3个主导因子,分别为坡度、坡位与腐质层厚度,且影响大小的顺序依次为:坡度>腐质层厚度>坡位。 (2)通过主成分分析,将原来的6个变量转变成2个主成分,其中1个可代表坡位因子,另1个可代表腐质层与坡度因子。这2个主成分对原来所有成分的贡献已达到78.5%,完全可以利用这2个主成分来描述其对毛竹林分生长的影响。与逐步回归分析的结果相比具有较好的一致性,2种分析方法都认为坡度、腐质层厚度与坡位对毛竹生长的影响较大。 (3)在主成分分析结果的基础上,通过聚类分析,对40块标准地进行了聚类,将立地条件相近的标准地划分为一类,主要分为三大类。在今后的生产中,可以对每一类采用相似的经营方式和管理措施,为实践生产提供了很大的帮助。 (崔鸿侠 熊德礼 陈双田 湖北省林业科学研究院 张维佗 湖北生态工程职业技术学院徐敏门 湖北省嘉鱼县林业局 蔡建国 湖北省国有林场工作站)