高山松天然林林分直径大小多样性研究

熊河先, 徐建

凉山州林业调查规划设计院，四川 凉山 615000

林分结构多样性是反映森林稳定性的重要指标。科学、合理的量化林分结构不仅能准确的评价森林的经济、生态、社会效益，还能为编制森林经营方案、制定森林经营措施提供依据[1]。学者通常使用林木大小多样性、林木格局多样性及物种多样性3个指标来量化林分结构[2]。林分直径结构是最重要、最基本的林分结构，林木直径便于测定且林分内各种大小直径的树木的分配状态将直接影响树木的树高、干形、材积、树种及树冠等因子的变化[3]。林分直径大小多样性指数是对林分直径结构的量化，其在一定程度上反应林分直径结构是否合理，已有研究表明：林分直径结构不合理，不利于正常的经营管理，需在今后的森林经营中逐步进行结构调整[4]；林分直径结构不合理，林分密度也需进行适当的调整[5]。

为了量化的研究生物多样性，自20世纪20年代起，Gleason首次提出了格里森指数（Gleason index）用于衡量均匀生境下物种的数目。近代多样性指数常被应用于生态学物种丰富度及均匀性的量化研究中[6]。近年来也有学者将多样性指数应用于林分直径大小多样性的研究[7-8]。描述林木直径大小多样性的指数有香农-威纳指数（Shannon-Wiener Index）、辛普森指数（Simpson's diversity Index）指数等8个指标[9]。量化分析是生态学研究的重要手段，典型对应分析（Canonical Correspondence Analysis，CCA）可以结合多个环境因子一起分析物种与环境之间的相关关系，其结果能够客观的反映群落的生态关系[10]。目前，国内生态学学者对CCA的研究大多是用于研究植物群落多样性，对林分直径大小多样性与周围环境的研究较少。

高山松（Pinus densata）属松科（Pinaceae）松属（Pinus）高大乔木，能生于干旱瘠薄的环境，是高海拔地区的先锋树种[11]。木里藏族自治县高山松分布范围广，是主要的用材和造林树种。目前，高山松的研究大都集中于直径结构、蓄积量、生物量等方面[12-14]，对于高山松直径大小多样性的研究较少。

综上所述，本研究基于48块木里藏族自治县高山松天然林实测样地，选择香农-威纳指数、辛普森指数，分别描述样地总体以及高山松林分直径大小多样性；使用相关分析、一般回归分析分别分析高山松林分直径大小多样性随林分因子、地形因子及林下植被因子的变化规律；再利用CCA排序探讨高山松直径大小多样性与环境梯度间的关系，分析环境对高山松天然林林分直径大小多样性的影响，以期为高山松林的经营管理、保护与发展利用规划提供依据。

1 研究方法

1.1 研究区概况

研究区位于四川省凉山彝族自治州木里藏族自治县，木里县地处四川省西南边缘，东跨雅砻江，西抵贡嘎山，南临金沙江，北靠甘孜州，面积1.3万km2。全县平均海拔3 100 m，年平均气温14.0 ℃。木里有长江上游重要的水源涵养林，是我国仅存不多的成片原始林区。

1.2 样地、环境因子选择以及数据整理

本研究的样地数据为木里藏族自治县二调实测样地数据，二调实测样布点为按照不同总体随机布点，样地面积666.67 m2，本研究在全县实测样地数据中筛选高山松组成大于8成且起源为天然的样地48块（见图1），样地实测胸径（DBH）株数分布见图2，环境因子分别从林分因子、地形因子、林下植被因子中选取。其中，林分因子选用株数密度（DOTt）、林分平均胸径（Dgt）；地形因子选用海拔（ALT）、坡度（SLO）；林下植被因子选用下木盖度（PcXM）、地被物盖度（PcDBW）。样地基本信息见表1。

图1 研究区位置示意图Fig.1 Schematic diagram of the locations in the study area

多样性指数计算：分别计算样地整体、高山松的Shannon-Wiener Index、Simpson's diversity Index[15,16]。

式中：H为Shannon-Wiener Index，样地总体为Ht，高山松为Hs；D为Simpson's diversity Index，样地总体为Dt，高山松为Ds；为每个径级林木株数占总株数s的百分比，以2 cm标准整化径阶并分级。

图2 胸径株数分布图Fig.2 Distribution of DBH in number of plants

表1 高山松天然林样地基本信息Tab.1 Basic information of the sampled plot

1.3 分析方法

采用SPSS 19.0对林分直径大小多样性指数与环境因子进行相关性分析，分析环境因子与直径大小多样性指数的相关性；运用一般回归模型分析林分直径大小多样性指数随环境因子的变化规律；利用Cannoco 5.0对样地进行CCA排序分析，分析环境梯度上高山松天然林林分直径大小多样性的变化。

2 结果与分析

2.1 直径大小多样性变化

从 Shannon-Wiener Index（H）、 Simpson's diversity Index（D）进行分析：从表2可以看出：Ht的取值在 1.2～3.1 之间、Hs的取值在 1.2～2.8 之间，Dt的取值在 0.59～0.95之间，Ds的取值在0.59～0.93之间，说明样地总体的直径大小多样性变动幅度比高山松大，样地内其他树种对高山松的直径大小多样性产生了一定的影响；Ht与Hs、Dt与Ds均值差异均较小，表现为样地总体>高山松，说明样地总体的直径大小多样性与高山松具有一致性，二者径级数基本一致，这应该是由于样地选择时只选择了高山松组成大于8成的样地造成的。

2.2 林分直径大小多样性与环境的关系

从表3可以看出，林分因子（Dgt、DOTt）与H相关性为极显著，Dgt与D的相关性为极显著，DOTt与D的相关性为显著，其中Dgt与林分直径大小多样性指数呈正相关；而DOTt与林分直径大小多样性指数呈负相关；地形因子中ALT与林分直径大小多样性指数相关性为极显著，而SLO与林分直径大小多样性指数相关性不显著；林下植被因子（PcXM、PcDBW）与林分直径大小多样性指数相关性均为不显著。

表2 样地林分直径大小多样性指数基本情况Tab.2 Basic situation of stand diameter size diversity in the sampled plot

表3 林分直径大小多样性指数与环境因子相关关系表Tab.3 Correlation between stand diameter size diversity index and environment factors

2.2.1 林分直径大小多样性指数随林分因子的变化

林分直径大小多样性指数随林分因子变化见图3，从拟合曲线来看，直径大小多样性指数（H、D）随Dgt的增大而增大，随DOTt的增大而降低，这与相关分析中直径大小多样性指数（H、D）与Dgt呈正相关，与DOTt呈负相关是一致的；从拟合精度来看，Dgt与直径大小多样性指数的拟合精度均在0.4以上。其中，Hs与Dgt的拟合精度最高，R2达到0.6989。DOTt与直径大小多样性指数拟合效果劣于Dgt的拟合效果，说明Dgt对直径大小多样性指数的影响大于DOTt。

2.2.2 林分直径大小多样性指数随地形因子的变化

林分直径大小多样性指数随地形因子变化见图4，从拟合曲线来看，H随SLO的增大先减少后增大，D随地形因子（SLO、ALT）的增大而增大。从拟合效果来看，直径大小多样性指数（H、D）与ALT的拟合效果优于SLO。说明ALT对直径大小多样性指数（H、D）的影响比SLO大，这与相关分析中直径大小多样性指数与地形因子相关性表现一致。

图3 林分直径大小多样性在林分因子（Dgt、DOTt）上的变化趋势Fig.3 Variation trend of stand diameter size diversity on stand factors (Dgt、DOTt)

2.2.3 林分直径大小多样性指数随林下植被因子的变化

林分直径大小多样性指数随林下植被因子变化见图5，从拟合精度来看，林分直径大小多样性指数（H、D）与林下植被因子（PcXM、PcDBW）的拟合效果均不好，R2均小于0.06。林下植被因子（PcXM、PcDBW）对林分直径大小多样性指数（H、D）的影响极小，与相关分析结果一致。

2.3 CCA 排序分析

图4 林分直径大小多样性在地形因子（ALT、SLO）上的变化趋势Fig.4 Variation trend of stand diameter size diversity on topographical factors（ALT、SLO）

图5 林分直径大小多样性在林下植被因子（PcDBW、PcXM）上的变化趋势Fig.5 Variation trend of stand diameter size diversity on understory vegetation factors（PcDBW、PcXM）

从表4可以看出，所选6个环境因子进行CCA排序分析效果较好，仅第一排序轴多样性指数与环境关系的解释量就能达到98.4%。根据前两轴绘制二维排序图（见图6），从趋势看，由第一轴从左至右，DOTt、PcDBW、PcXM逐渐降低，ALT、SLO、Dgt逐渐升高；第二轴从下往上，ALT、DOTt逐渐降低，Dgt、PcXM、SLO、PcDBW逐渐升高。Ht、Hs、Dt及Ds聚集在一起，说明样地直径大小多样性与高山松直径大小多样性变化趋势相似。从排序轴长短来看，Dgt在第一轴最长，说明Dgt对直径大小多样性指数的影响最显著，这与相关分析、一般回归分析结果一致。PcDBW、PcXM、SLO在第一轴较短，说明这3个因子对林分直径大小多样性指数的影响最小，这与相关分析结果一致，相关分析中这3个因子与林分直径大小多样性指数的相关性均表现为不显著。

3 结论与讨论

本研究以四川省凉山彝族自治州木里藏族自治县48块高山松天然林样地为研究对象，使用Shannon-Wiener Index、Simpson's diversity Index量化高山松直径大小多样性，量化直径大小多样性后，样地总体的直径大小多样性与高山松具有一致性，由于样地选择时只选择了高山松组成大于8成的样地，说明当优势树种成分较大时，林分总体直径多样性与林分优势树种直径多样性表现一致。采用3种方法分别分析直径多样性指数与环境因子（林分因子、地形因子、林下植被因子）的相关关系，研究区内样地内林木直径大小多样性表现为样地总体>高山松；3种分析方法均显示，环境因子对林分直径大小多样性的影响：林分因子>地形因子>林下植被因子。其中，林分因子对林分直径大小多样性指数的影响最大，林下植被因子对直径大小多样性指数的影响最小，且两个林下植被因子均是表示覆盖率，说明上层乔木直径结构变化受下层植被数量的影响较小。

表4 CCA 排序重要参数Tab.4 Important parameters of canonical correspondence analysis

图6 环境因子 CCA 二维排序图Fig.6 CCA two-dimensional sequence diagram of environmental factors

本次研究中，直径大小多样性指数与林分株数密度、坡度、海拔在相关分析时相关性表现为显著，但在一般回归模型拟合时拟合精度较低。其原因可能是在模型拟合时仅使用了一般回归模型，而林木直径是森林生长的结果，森林的生长过程受周围环境的影响，使其生长过程不再是一个独立的过程，受空间效应的影响，经典统计学模型对其的估计是有偏的，没有充分考虑空间效应对直径大小多样性指数的影响使一般回归模型拟合效果较差[17]。

林分直径大小多样性指数一定程度上反映了林分直径分布规律，可以作为制定、检查营林技术措施效果的依据之一。特别是设计抚育间伐方案时要充分考虑林分直径结构和空间结构，在保持合理林分密度的同时充分考虑直径结构规律[3]。

本研究引入生态学Shannon-Wiener Index、Simpson's diversity Index，从林分总体、样地内高山松两个方面分析高山松天然林林分直径大小多样性随环境的变化规律，但本次研究仅考虑了与距离无关的多样性测度方法，考虑变化时也只考虑了固定效应对直径大小多样性的影响，以上两点有待进一步研究。