Classification, ordination and correlation analysis of dominant species of ecological non-commercial forests in Songyang, Zhejiang

松阳县(28°14′~28°36′N，119°10′~119°42′E)位于浙江省西南部，全县总面积为14.06 万hm²，公益林面积为6.38 万hm²，占林业用地面积的55.0%，森林覆盖率75.0%，林木绿化率75.9%^[15-16]。气候属亚热带季风气候，温暖湿润，四季分明，无霜期长，具有一定的山地立体气候。林地土壤有红壤、黄壤、岩性土、潮土等4类，其中红壤和黄壤为主要的土壤类型。松阳县地带性植被属于中亚热带常绿阔叶林，现有森林的乔木林林分面积中马尾松Pinus massoniana林、杉木Cunninghamia lanceolata林、阔叶林、经济林的面积比为39∶36∶17∶8^[17]。群落类型以阔叶林、针阔混交林和针叶林群落为主，优势树种有马尾松、杉木、枫香Liquidambar formosana、木荷Schima superba、青冈Cyclobalanopsis glauca、苦槠Castanopsis sclerophylla、麻栎Quercus acutissima等^[18]。

以浙江省松阳县的公益林作为总体样本，按一定比例随机抽取监测小班，在每个小班中选取1个最能代表该小班状况的区域设置20 m×20 m的固定监测样地，共设128个。使用哈尔滨牌DQL-12Z型12倍正像森林罗盘仪进行打样地，闭合差控制在1/200以内。对样地内乔木树种(胸径≥5.0 cm)采用每木调查，记录树种名称、胸径、树高、冠幅、枝下高等数据。记录样地基本地理信息和环境因子，包括海拔、坡度、坡位、坡向、土壤类型、腐殖质厚度和凋落物厚度^[19]，并利用手持GPS仪记录样地经纬度坐标。该县被抽中的128个样地中包含马尾松林、杉木林、针阔混交林、阔叶林、竹林等林分类型。样地分布涵盖了中山、低山、丘陵、盆地等地形，海拔为115~1 278 m。

各个种群的重要值反映种群在群落中的作用和地位，计算公式^[18]如下：乔木层重要值=(相对显著度+相对频度+相对多度)/3×100%；相对显著度=某一树高胸高断面积之和/所有树种胸高断面积之和×100%；相对多度=某一树种株数之和/所有树种株数之和×100%；相对频度=某一树种频度/所有树种的频度和×100%。

TWINSPAN是在指示种分析(indicator species analysis)的基础上改进而来的。它以二岐式的分割法来划分类型，其划分根据“指示种”(或称为指示指标)将样地类型与物种组成依次划分为各个等级的类型单位，完成物种分类和样地分类^[20]。本研究依据物种重要值排序，选取重要值大于0.1的乔木层树种，利用WinTWINS 2.3软件分析^[20]。重要值采用5个等级作为样地划分的指示，分别为0~0.05、0.05~0.10、0.10~0.30、0.30~0.50、＞0.50。

CCA是一种非线性多元直接梯度分析方法，它把对应分析与多元回归结合起来，每一步计算结果都与环境因子进行回归，研究物种与环境的关系^[21]。本研究选取坡位、坡度、坡向、海拔高度、土壤类型、腐殖质厚度、凋落物厚度为环境因子。将这些环境因子建立矩阵，在CANOCO 4.5中进行分析。此外，样地分布环境涵盖较广，也可以使得CCA排序的结果更加可靠。

采用Shannon-Wiener多样性指数、Simpson指数、Pielou均匀指数和物种丰富度作为衡量群落生物多样性的指标^[22]。用SPSS单因素方差分析比较不同群丛的物种多样性差异。使用R 3.6.1的boxplert包绘制箱型图^[23]。由于χ²检验可判断种间关联性程度的强弱，但是无法定量关联程度大小，因此采用联结系数和共同出现百分率对研究区11个优势种(重要值＞1.1%)进行种间联结分析^[24-27]。利用Excel 2019进行数据整理，使用R 3.6.1的spaa包进行种间关联分析，corrplot包进行种间关联半矩阵图输出。

表1列出了128个样地中重要值排名前20的物种，主要有杉木、马尾松、毛竹、木荷、青冈、甜槠、苦槠、石栎等树种。

采用研究区128个样地重要值矩阵数据在TWINSPAN中进行群丛划分，得到6个群丛类型，每个群丛的样地及环境信息如图1和表2所示。其中：群丛Ⅰ为板栗+杉木+杨梅群丛，群丛Ⅱ为毛竹+杉木+马尾松群丛，群丛Ⅲ为马尾松+杉木+油茶Camellia oleifera群丛，群丛Ⅳ为马尾松+木荷+杉木群丛，群丛Ⅴ为杉木+木荷+马尾松群丛，群丛Ⅵ为杉木+马尾松+木荷群丛。

Figure 1.  Dendrogram illustrating the presence of six vegetation associations using TWINSPAN analysis of 128 plots in Songyang non-commercial forest

第1轴和第2轴与环境因子的特征值分别为0.302和0.167，物种环境相关性分别为0.718和0.555，包含了65.7%的生态信息。因此本研究使用第1轴和第2轴的数据来分析样地中物种和周围环境之间的关系(表3)。CCA排序图的每个箭头代表1种环境因子，箭头的长度表示该环境因子对物种的影响大小^[28-29]。从图2可看出：CCA排序第1轴与海拔高度呈反比，相关系数为−0.583，与坡向呈正相关0.271(表3)。第2轴与土壤类型、坡位和凋落物厚度相关性较高。

Figure 2.  CCA ordination diagram of 128 plots in Songyang non-commercial forest

从图3可见：6个群丛在物种多样性指数分析中存在显著差异(P＜0.001)。其中多样性指数(Shannon-Wiener指数和Simpson指数)变化趋势基本一致，群丛Ⅴ的Shannon-Wiener指数和Simpson指数显著高于其他群丛(P＜0.001)，且群丛Ⅰ和Ⅱ的Simpson指数与群丛Ⅲ差异显著(P＜0.001)；群丛Ⅲ的Pielou指数与除群丛Ⅴ之外的群丛差异显著(P＜0.001)，群丛Ⅴ与群丛Ⅰ、Ⅱ存在显著差异(P＜0.001)；物种丰富度显示群丛Ⅴ显著高于其他群丛(P＜0.001)。

Figure 3.  Changes of species diversity index among six associations

选择重要值排名前11位的树种，分别计算其两两之间的χ²值、联结系数和共同出现百分率，并分别做半矩阵图。将χ²结果对数化，lg3.841=0.584，lg6.635=0.822。如图4A所示：4个种对lgχ²处于[0.584，0.822]内，6个种对lgχ²＞0.822，表明有4个(7.2%)种对呈显著关联，6个(10.9%)种对间关联极显著。表现极显著联结的种对有：马尾松对赤杨叶、木荷对甜槠、木荷对青冈、木荷对赤杨叶、青冈对苦槠、青冈对甜槠。表现显著联结的种对有：杉木对青冈、杉木对苦槠、青冈对石栎、青冈对白栎。图4B表明：优势物种种间关系有36对(65.5%)呈正关联，19对(34.5%)呈负关联，正负关联比为1.9。其中，联结系数处于[0.5，1.0]值域的有0对，处于[0.2，0.5]值域的有6对(10.9%)：木荷对青冈、木荷对甜槠、青冈对甜槠、青冈对苦槠、甜槠对白栎、甜槠对赤杨叶。而处于[−1.0，−0.5]值域的有7对(12.7%)：马尾松对黄山松、毛竹对石栎、毛竹对白栎、青冈对黄山松、苦槠对黄山松、石栎对黄山松、黄山松对赤杨叶，处于[−0.5，−0.2]值域的有6对(10.9%)：马尾松对毛竹、马尾松对甜槠、马尾松对赤杨叶、毛竹对木荷、毛竹对青冈、毛竹对甜槠。根据共同出现百分率矩阵图(图4C)可知：共同出现百分率处于[0.4，0.7]的种对有3对，占总种对数的5.5%，这些种对呈显著正联结。共同出现百分率处于[0.2，0.4]的种对有17对，占总种对数的30.9%，这些种对为不显著正联结。处于[0，0.2]的种对有35对，占总种对数的63.6%，这些种对几乎无联结。

Figure 4.  Semi-matrices graphs of interspecific association of 11 dominant species

CCA第1轴排序主要与海拔、坡向以及土壤有关，本质上反映植物对水热和阳光需求的变化，这与钟娇娇等^[12]对秦岭山地天然次生林进行排序分类的结论一致。群落的形成与环境因子密不可分，通常是气候、地形和土壤等因子的综合作用，而这些环境因子在不同地理空间上对于植物群落格局形成的影响也不同^[30]。本研究结果中影响群落分布的第1主导因子是海拔，第2主导因子是土壤类型。陈云等^[31]在小秦岭的研究表明：森林群落物种组成和分布格局的前2个主要环境因子是分别是海拔和坡度，其中第2主导因子与本研究结果存在差异，这种差异在与叶诺楠等^[32]的研究比较中也得到了体现。可见在不同的地理和气候区划中，环境因子对植被群落分布和类型的影响是不同的。

在不同的立地条件下松阳县公益林的群落组成存在显著差异，同时物种多样性在环境梯度上也呈现一定分异规律，这与早期的研究也较为接近^[32]。本研究中，杉木+木荷+马尾松群丛(V)的Shannon-Wiener指数高于其他群丛，可能由于群丛Ⅴ在各海拔均有分布且林分为混交林，因此群落物种结构复杂、多样性高。而其他群丛有些在单一环境下分布，群落组成简单、结构单一，故群落物种多样性低。此外，毛竹入侵造成群落中的光照和土壤等环境条件变化^[33]，这可能将造成乔木和灌木多样性降低^[34]，这也是毛竹为主导的群丛多样性较低的主要原因之一。

研究区大多数乔木树种呈现不显著联结，说明乔木层的物种间联系不明显，趋于独立分布，这也有可能是生境异质性导致的^[35]。同时，χ²检验难以定量确定种间联结程度的大小，因此再采用AC值和PC值同时进行种间联结强度测定能够从不同角度阐明物种种间关联特性。样地内呈现正联结的种对数多于负联结，说明乔木层群落结构尚不稳定^[36]。物种入侵将影响群落间的种间关联性，一些入侵的物种会降低群落的正关联性，增加负联结性^[37]。

CCA结果显示群丛Ⅰ的样地数量不多，主要位于低海拔的阳坡，早期为杉木林内伐除杉木后人为种植的板栗和杨梅等经济林，为便于经营管理，近年来被划入公益林内。该类林分结构单一，生物多样性低，但立地条件尚可，可以较大强度地间伐板栗和杉木，补植木荷、青冈等珍贵树种，逐步改造为常绿落叶阔叶混交林。需要引起高度重视的是群丛Ⅱ，即毛竹群丛，样地数量不少，但生物多样性是研究区所有类型中最低的。以毛竹为代表的鞭生竹克隆生长使其能够入侵相邻林分，并逐渐发展为新的“优势种”，成为一种“局部入侵”物种，影响周围森林群落结构、土壤质量、幼苗更新及幼树生长，从而影响整个森林生态系统。因此，加强公益林内毛竹群丛的林分改造，防止竹林的无序扩张是今后公益林管理的重要问题。在CCA排序图中，样地众多的群丛Ⅴ和群丛Ⅵ主要分布在海拔500~1 000 m且土层厚度较好的生境中，林内物种较为丰富，以木荷为代表的阔叶树种在乔木层占有一定的优势，今后可以适当间伐杉木、马尾松等针叶树，对阔叶树种实施定株抚育或目标树经营，促进阔叶树的生长，加速形成以阔叶树为主的林分类型。

本研究对松阳县生态公益林群落进行数量分类与排序，分析了多样性指数和种间关系，提出了公益林的差别管理措施建议。在日后公益林管理工作中，可考虑对不同立地条件不同林分的公益林进行不同的补植措施，对毛竹林进行边缘控制，对一些能适应各类环境的群落组合进行适当推广等措施。