Numerical classification and ordination of vegetation communities in Jinzifeng National Forest Park, Qingyuan, Zhejiang

巾子峰国家森林公园位于浙江省庆元县西部，面积5 752 hm²，27°27′~27°39′N，118°50′~119°06′E，最高峰巾子峰海拔1 563 m。该区气候属中亚热带季风气候区，温暖潮湿，雨水充沛，四季分明。年平均气温为17.2 ℃，最高月平均气温26.9 ℃，最低月平均气温7.0 ℃，年平均降水量1 689.0 mm，平均相对湿度76%；年平均日照时数1 828.2 h，日照率41%；年平均无霜期256.0 d。土壤类型为红壤和黄壤2个土类，均为地带性土壤。红壤广泛分布于海拔800 m以下的低山丘陵，黄壤分布于800 m以上的中低山。主要植被类型有常绿阔叶林、常绿-落叶阔叶混交林、针阔叶混交林、针叶林、毛竹Phyllostachys edulis林等，地带性典型植被为常绿阔叶林。

2019年4−5月，在巾子峰国家森林公园海拔500~1 500 m内，隔100 m设置1个取样区间，每个取样区间布设2~8个20 m×20 m样方，共布设10个取样区间48个样方。记录样方内胸径(DBH)≥2 cm的木本植物个体种名、高度、胸径、盖度，每种灌木及草本的种名、高度、盖度、株数或丛数，同时记录每个样方的海拔、坡向、坡度、土壤类型、土壤湿度、土壤温度。其中，海拔高度用全球定位系统测量，坡向和坡度用罗盘仪测定，土壤温度和土壤湿度用Veinasa-TS土壤温湿速测仪测量。48个样方的各个生境因子如表1所示。

用重要值表示每个物种在森林群落中的优势程度，按照乔木树种、灌木及草本植物不同的生长特征，选择不同的重要值计算方法。计算公式如下：乔木树种重要值=(相对密度+相对优势度+相对频度)/3。灌木及草本植物重要值=(相对盖度+相对高度)/2。

海拔高度、坡度、土壤湿度、土壤温度直接测量数值。用数字表示不同的土壤类型，其中：0代表红壤，1代表黄壤。为了分析森林群落的向阳程度，将坡向因子进行转换，以数字表示不同的坡向等级，其中：1表示北坡(0°~22.5°和337.5°~360°)，2表示东北坡(22.5°~67.5°)，3表示西北坡(292.5°~337.5°)，4表示东坡(67.5°~112.5°)，5表示西坡(247.5°~292.5°)，6表示东南坡(112.5°~157.5°)，7表示西南坡(202.5°~247.5°)，8表示南坡(157.5°~202.5°)。数值越大说明越向阳。

使用Wintwins 2.3软件，采用双向指示种分析分类方法，对由48个样方及样方内194个种的重要值组成的48×194维矩阵进行分析，得出数量分类结果。

使用CANOCO 4.5软件^[11]，采用去趋势对应分析排序方法，对由48个样方及样方内194个种的重要值组成的48×194维矩阵进行分析，得出排序结果。采用典范对应分析对由48个样方及样方内194个种的重要值组成的48×194维矩阵和48个样方及6个环境因子组成的48×6维矩阵进行分析，得出排序结果。

野外调查共记录种子植物193种，隶属66科129属。科物种组成占总物种数比例分别为壳斗科Fagaceae 28.79% (19种)、樟科Lauraceae 27.27% (18种)、山茶科 Theaceae 13.64% (9种)、杜鹃花科Ericaceae 13.64% (9种)、蔷薇科Rosaceae 13.64% (9种)、冬青科Aquifoliaceae 9.09% (6种)、禾本科Gramineae 9.09% (6种)、茜草科Rubiaceae 9.09% (6种)、金粟兰科Chloranthaceae 7.58% (5种)、百合科 Liliaceae 7.58% (5种)。根据表2中群落不同层次物种重要值的排序可知(排名前20)，乔木层优势物种主要为杉木、木荷、黄山松、青冈、甜槠、红楠、少叶黄杞和拟赤杨。灌木层优势种为杜茎山、草珊瑚、隔药柃和映山红。草本层优势种为蕨、芒萁、黑足鳞毛蕨、中华薹草、白茅、江南卷柏和荩草。

根据中国植被群落的分类原则^[12]和双向指示种分析分类结果，将48个样方分成10个群丛(图1)，采用生境指示种及群丛优势种的方式对群丛进行命名。各群丛环境因子描述见表3。

Figure 1.  Dendrogram of TWINSPAN classified result in Jinzifeng Forest Park of Qingyuan, Zhejiang

群丛Ⅰ：黄山松-映山红-东方古柯-少花万寿竹群丛Association Pinus taiwanensis-Rhododendron simsii-Erythroxylum sinense-Disporum uniflorum，包括38、39、43号样方。乔木层优势种为黄山松，主要伴生种为木荷。灌木层优势种为映山红、东方古柯，主要伴生种为宜昌荚蒾、麂角杜鹃Rhododendron latoucheae。草本层优势种为少花万寿竹、藜芦。

群丛Ⅱ：黄山松-木荷-映山红-多花黄精群丛Association Pinus taiwanensis-Schima superba-Rhododendron simsii-Polygonatum cyrtonema，包括34、41、42、44、45号样方。乔木层优势种为黄山松、木荷、青冈、甜槠，主要伴生种为浙闽樱Prunus schneideriana、云和新木姜子Neolitsea aurata var. paraciculata、苦茶槭、湖北鹅耳枥Carpinus hupeana、多脉鹅耳枥Carpinus polyneura。灌木层优势种为映山红、麂角杜鹃，主要伴生种为映山红、江南越橘、中华野海棠。草本层优势种为多花黄精、少花万寿竹。

群丛Ⅲ：甜槠-木荷-黄丹木姜子-少花万寿竹群丛Association Castanopsis eyrei-Schima superba-Litsea elongata-Disporum uniflorum，包括32、33号样方。乔木层优势种为甜槠、木荷、青冈，主要伴生种为马银花、多脉鹅耳枥、江南油杉Keteleeria fortunei var. cyclolepis、浙江樟Cinnamomum chekiangense、红楠、四照花Cornus kousa subsp. chinensis。灌木层优势种为黄丹木姜子Litsea elongata，主要伴生种为毛花连蕊茶、隔药柃、江南越橘。草本层优势种为少花万寿竹、藜芦、蕨。

群丛Ⅳ：杉木-毛花连蕊茶-江南卷柏群丛Association Cunninghamia lanceolata-Camellia fraterna-Selaginella moellendorffii，包括15、19、22、24、31、35、36、37、40号样方。乔木层优势种为杉木，主要伴生种为木荷、拟赤杨、虎皮楠、刨花楠、光皮桦、青冈、薯豆、少叶黄杞。灌木层优势种为毛花连蕊茶，主要伴生种为细枝柃、杜茎山、隔药柃、矩形叶鼠刺、毛冬青Ilex pubescens、黄绒润楠、乌药、赤楠、山鸡椒Litsea cubeba、刺毛杜鹃、豹皮樟、山莓、映山红。草本层优势种为江南卷柏、蕨、赤车Pellionia radicans，主要伴生种为贯众、黑足鳞毛蕨、芒萁、日本蛇根草。

群丛Ⅴ：杉木-木荷-杜茎山-芒萁群丛Association Cunninghamia lanceolata-Schima superba-Maesa japonica-Dicranopteris pedata，包括1、2、3、5、6、7号样方。乔木层优势种为杉木、木荷、薯豆、拟赤杨，主要伴生种为少叶黄杞、细叶青冈、树参、青冈、刨花楠、米槠、红楠。灌木层优势种为杜茎山、刺毛杜鹃，主要伴生种为江南越橘、草珊瑚、矩形叶鼠刺、檵木、黄绒润楠。草本层优势种为芒萁、黑足鳞毛蕨，主要伴生种为蕨、赤车、江南卷柏、中华薹草。

群丛Ⅵ：青冈-木荷-细枝柃-黑足鳞毛蕨群丛Association Cyclobalanopsis glauca-Schima superba-Eurya loquaiana-Dryopteris fuscipes，包括4、11、12、18、21号样方。乔木层优势种为青冈、木荷、锥栗，主要伴生种为冬青、虎皮楠、少叶黄杞、米槠、马银花。灌木层优势种为细枝柃，主要伴生种为杜茎山、石楠 Photinia serrulata、草珊瑚、乌药、狗骨柴Diplospora dubia、赤楠、矩形叶鼠刺、刺毛杜鹃、江南越橘。草本层优势种为黑足鳞毛蕨、芒萁、中华薹草。

群丛Ⅶ：拟赤杨-虎皮楠-山矾-日本蛇根草群丛Association Alniphyllum fortunei-Daphniphyllum oldhami-Symplocos sumuntia-Ophiorrhiza japonica，包括8、10、28号样方。乔木层优势种为拟赤杨、虎皮楠，主要伴生种为米槠、树参、光皮桦、少叶黄杞、木荷、马银花。灌木层优势种为山矾Symplocos sumuntia，主要伴生种为杜茎山、矩形叶鼠刺、黄绒润楠、草珊瑚、乌药、隔药柃。草本层优势种为日本蛇根草、黑足鳞毛蕨，主要伴生种为中华薹草、蕨、江南卷柏。

群丛Ⅷ：刨花楠-隔药柃-蕨群丛Association Machilus pauhoi-Eurya muricata-Pteridium aquilinum var. latiusculum，包括16、20、25、26、27、29、46号样方。乔木层优势种为刨花楠、甜槠、栲，主要伴生种为少叶黄杞、深山含笑Michelia maudiae、青冈、大叶青冈、冬青、虎皮楠、锥栗。灌木层优势种为隔药柃，主要伴生种为毛花连蕊茶、细枝柃、山矾、乌药、黄绒润楠、石楠、矩形叶鼠刺、赤楠、豹皮樟、乌冈栎Quercus phillyreoides。草本层优势种为蕨，主要伴生种为中华薹草、江南卷柏、贯众、地菍。

群丛Ⅸ：马尾松-木荷-草珊瑚-中华薹草群丛Association Pinus massoniana-Schima superba-Sarcandra glabra-Carex chinensis，包括9、13、17、30号样方。乔木层优势种为马尾松、木荷、青冈，主要伴生种为少叶黄杞、拟赤杨、虎皮楠、锥栗、细叶青冈、红楠、蓝果树。灌木层优势种为草珊瑚，主要伴生种为檵木、黄绒润楠、朱砂根、映山红、江南越橘、隔药柃、山莓、老鼠矢。草本层优势种为中华薹草、芒萁、日本蛇根草。

群丛Ⅹ：毛竹-山鸡椒-地菍群丛Association Phyllostachys edulis-Litsea cubeba-Melastoma dodecandrum，包括14、23、47、48号样方。乔木层优势种为毛竹，主要伴生种为杉木、拟赤杨、马尾松。灌木层优势种为山鸡椒，主要伴生种为隔药柃、乌药、黄绒润楠、檵木、胡枝子Lespedeza bicolor、老鼠矢。草本层优势种为地菍，主要伴生种为芒萁、杏香兔耳风Ainsliaea fragrans、蕨、中华薹草。

48个样方的去趋势对应分析结果表明：4个排序轴的特征值分别为0.549、0.400、0.295、0.225，因前2个排序轴的特征值比较大，能显示出一定的生态意义，所以采用前2个排序轴的数据分别作二维排序图(图2)。从图2可以看出：去趋势对应分析排序的横轴和纵轴基本反映了不同的生态意义。第1排序轴(横轴)反映了各群丛所在环境的海拔和土壤类型梯度，即沿第1排序轴从左往右，随着海拔高度逐渐升高，土壤类型由红壤变为黄壤；第2排序轴(纵轴)反映了各群丛所在的坡度和土壤情况变化，沿第2排序轴从下往上，坡度由陡坡转为缓坡，土壤温度逐渐降低，土壤湿度逐渐升高。

Figure 2.  DCA two-dimensional ordination diagram for samples of plant communities in Jinzifeng Forest Park of Qingyuan, Zhejiang

结合双向指标种分析的分类结果对排序图进行分析，发现不同的森林群落类型在排序图中呈现出一定规律的分布。黄山松、黄山松-木荷群丛分布在海拔较高的中山带阴坡、半阴坡，土壤类型为黄壤，土壤温度和湿度均较低，位于排序图的右下方。甜槠-木荷群丛分布在中山半阳坡、半阴坡，坡度较陡，土壤类型为黄壤，土壤温度较高，土壤湿度较小，位于排序图的左下方。杉木-毛花连蕊茶群丛分布在中山带半阳坡、半阴坡，土壤类型为红壤或黄壤，位于排序图的中心。杉木-木荷、马尾松-木荷群丛分布在低山半阳坡、半阴坡，土壤类型为红壤，位于排序图的中心靠左位置。青冈-木荷群丛分布在低山半阳坡、半阴坡，土壤类型为红壤，位于排序图的左方。拟赤杨-虎皮楠群丛分布在低山半阳坡，土壤类型为红壤，位于排序图的左上方。刨花楠群丛分布在低山阳坡、半阳坡，土壤类型为红壤，位于排序图的左下方。毛竹群丛分布在低山半阳坡，土壤类型为红壤，位于排序图的上方。

浙江庆元巾子峰国家森林公园植被群落物种重要值和样方矩阵的典范对应分析排序所有轴的梯度长度最大为4.115，所以选择单峰模型。蒙特卡罗检验结果显示：所有典范排序轴通过统计检验(F=1.629, P=0.002)，说明所选择的环境因子对植被群落物种分布具有显著影响。由表4可知：典范对应分析排序轴中前2轴累计解释了物种-环境因子关系的55.2%，说明典范对应分析前2轴能较好地反映森林植被与环境因子之间的关系。根据环境因子与排序轴的相关性大小可知，海拔、土壤类型与第1排序轴呈极显著正相关(P＜0.01)，土壤温度与第1排序轴呈显著负相关(P＜0.05)，相关性从大到小依次为海拔、土壤类型、土壤温度，表明第1排序轴反映的是海拔和土壤类型的梯度变化。土壤温度与第2排序轴呈极显著负相关(P＜0.01)，与坡度呈显著正相关(P＜0.05)，表明第2排序轴反映了土壤温度的梯度变化。

样方-环境因子的典范对应分析排序图较为直观地呈现了植被群落各群丛类型分布与环境因子的关系(图3)，带箭头射线代表环境因子，射线长度表示植被群落与环境因子相关性的大小，射线所在象限表示环境因子与排序轴相关性的正负。典范对应分析排序第1轴代表了海拔和土壤类型，从左往右海拔逐渐升高，土壤由红壤转变为黄壤；第2排序轴代表了土壤温度，从下往上土壤温度逐渐增加。群丛Ⅰ和群丛Ⅱ的分布与海拔和土壤类型关系比较密切，分布在海拔最高的区域，且均是黄壤。群丛Ⅲ分布在较高海拔。群丛Ⅳ分布比较宽泛，在海拔600~1 300 m均有分布，土壤有红壤还有黄壤。群丛Ⅴ和群丛Ⅶ与土壤湿度关系较密切。群丛Ⅷ、群丛Ⅸ和群丛Ⅹ则与土壤温度和土壤湿度关系密切。

Figure 3.  CCA diagram for samples and environmental variables of the plant communities in Jinzifeng Forest Park of Qingyuan, Zhejiang

优势乔木树种对群落结构和群落环境的形成有明显控制作用，个体数量不一定很多，但却能决定群落结构和内部环境条件，是植被群落的建造者。对森林群落中20个优势乔木树种进行典范对应分析排序，排序结果如图4所示。海拔、土壤类型、土壤湿度和坡度对物种分布起决定作用。黄山松分布在中山带阴坡，土壤类型为黄壤，位于排序图的右方。甜槠、锥栗分布在低山半阳坡，土壤类型为黄壤，位于排序图的上方。杉木分布在低山、中山带半阳坡、半阴坡，土壤类型为红壤或黄壤，位于排序图的中间靠上位置。青冈、木荷分布在低山、中山带半阳坡、半阴坡，土壤类型为红壤或黄壤，位于排序图的中间靠下位置。拟赤杨、石栎分布在低山半阳坡，土壤类型为红壤，位于排序图的中间靠左位置。马尾松分布在低山半阳坡，土壤类型为红壤，位于排序图的左方。刨花楠、虎皮楠、栲、米槠、红楠、冬青、薯豆、少叶黄杞、细叶青冈、黄绒润楠分布在低山半阳坡、半阴坡，土壤类型为红壤，位于排序图的左下方。毛竹分布在低山阳坡，土壤类型为红壤，位于排序图的左上方。

Figure 4.  CCA diagram for trees and environmental variables in Jinzifeng Forest Park of Qingyuan, Zhejiang

本研究根据中国植被群落的分类原则和双向指标种分析方法，将浙江庆元巾子峰国家森林公园内48个样方分成10个群丛。对10个群丛及其相应样方进行去趋势对应分析排序。结果表明：各群丛在排序图中的分布均较为集中，说明去趋势对应分析排序结果与双向指标种分析分类结果较为一致，综合应用双向指标种分析分类和去趋势对应分析排序能够更好地解释森林群落之间的差异性和连续性，这与以往的研究结果相同^[13-17]。双向指标种分析数量分类是依据指示种和地区分布差异划分的，而去趋势对应分析是依据排序轴的综合信息，将样方或植物物种排列在一定的空间^[2]。去趋势对应分析排序轴能够反映一定的生态梯度。本研究横轴主要反映了样方的海拔高度和土壤类型，纵轴主要反映了样方的土壤温度和坡度。各群丛在排序图中的位置能有效反映出各群丛的环境特征及群丛间的连续分布关系，可以看出植被类型间出现了交错现象，表明植被类型与环境因子间关系比较复杂，植被类型的分布受海拔、土壤环境、坡度、坡向等多种环境因子的共同影响。

典范对应分析排序与去趋势对应分析排序相比，去趋势对应分析排序仅使用物种数据，在反映样方之间或物种之间的关系上具有较大优势，但不能较好地反映与环境因子之间的相关性；典范对应分析同时使用了物种数据和环境因子数据，能较好地表达群落的环境梯度，排序轴能更好地对环境因子进行解释^{[2, 4, 18]}。本研究划分的10种群丛类型与样方在典范对应分析排序图上的分布格局基本吻合，能较好地揭示植被类型与环境因子之间的关系，反映特定的植被类型对不同的环境因子差异具有指示作用。从典范对应分析排序图上可以看出：海拔、土壤类型和土壤温度这3个主要环境因子共同影响了植被的分布。物种丰富度海拔梯度格局及其形成机制一直是宏观生态学研究的热点。研究表明：海拔梯度是影响植被群落物种分布格局的主要控制因子^[4,19]。因为通过海拔梯度的变化可以改变局部水分状况、土壤理化性质、干扰程度等其他生态环境，直接或间接地对群落分布格局产生影响^[20-22]。本研究也证明海拔对山地植被分布具有决定性作用，海拔是影响巾子峰国家森林公园植被群落分布的最重要环境因子。

综上所述，本研究中影响植被群落分布最重要的环境因子是海拔，其次为土壤类型、土壤温度和坡度，土壤湿度和坡向对植被群落分布的影响相对较小。对20个优势树种进行的典范对应分析排序发现，优势树种的分布特征与各群丛及相应样方的分布特征具有较好的相似性，影响各优势树种分布最重要的环境因子也是海拔，其次为土壤类型、土壤温度和坡度。这与各群丛及相应样方典范对应分析排序结果一致，说明各优势树种的分布在一定程度上决定了各群丛类型的分布。