三峡库区马尾松林林分结构特征对灌木层物种多样性的影响

《三峡库区马尾松林林分结构特征对灌木层物种多样性的影响》由会员分享，可在线阅读，更多相关《三峡库区马尾松林林分结构特征对灌木层物种多样性的影响（19页珍藏版）》请在装配图网上搜索。

1、三峡库区马尾松(Pinus massoniana)林林分结构特征对灌木层物种多样性的影响胡文杰;潘磊;雷静品;唐万鹏;庞宏东;崔鸿侠;王晓荣【摘要】林分结构特征对灌木层物种多样性有重要影响,对其之间关系进行研究,可 为提高生态系统功能稳定性、进行合理的林分经营管理提供理论依据.以马尾松林 为研究对象,采用多元逐步回归、因子分析和典型相关分析等方法分上林层、中林 层和下林层研究林分结构对灌木层多样性的影响结果表明:Pearson相关分析结 果表明,上林层是影响灌木层物种多样性的主要层次,而林分总体结构特征与灌木层 物种多样性无显著性关系;(2)用13 个林分结构参数和灌木层物种多样性各参数进 行

2、多元逐步回归拟合方程都通过显著性检验,R2为0.246-0.630,且仅剔除了树高 和大小比数2 个变量而保留了其它11 个林分结构变量,说明林分结构的综合作用对 灌木层物种多样性有较大影响;(3)林分结构可分为 3 类,分别为林分竞争状况因子、 水平结构及冠幅因子和垂直结构及生长因子;(4)上林层林分结构总体上与灌木层物 种多样性之间存在着较强的相关性，典型相关系数为0.822(P=0.014),其中垂直结构 及生长因子为主要影响因子,权重为 0.962.在马尾松林中,上林层林分结构中的垂直 结构与生长因子和灌木层物种多样性表现为正相关关系,在经营管理中,可通过调控 上林层林木的胸径、冠长及

3、冠层垂直结构等以达到提高林分灌木层物种多样性的目 的.【期刊名称】生态环境学报【年(卷),期】2019(028)007【总页数】9页(P1332-1340)关键词】马尾松;林分结构;灌木层物种多样性;典型相关分析 【作 者】胡文杰;潘磊;雷静品;唐万鹏;庞宏东;崔鸿侠;王晓荣【作者单位】湖北省林业科学研究院,湖北武汉 430075;湖北省林业科学研究院,湖 北武汉 430075;中国林业科学研究院林业研究所,北京 100091;湖北省林业科学研 究院,湖北武汉 430075;湖北省林业科学研究院,湖北武汉 430075;湖北省林业科 学研究院,湖北武汉 430075;湖北省林业科学研究院,湖北

4、武汉 430075【正文语种】中文【中图分类】S718.5;X17林分结构特征和物种多样性是评估森林生态系统功能和林分经营质量的关键指标（Burrascano et al. , 2011 ;方升佐等，2012 ）。灌木层作为森林生态系统的组 成部分，对保持生态系统功能稳定性具有重要作用，其物种多样性在提高林分水源 涵养能力、改善土壤养分状况等方面（褚建民等， 2007 ）发挥着不可替代的作用， 甚至对乔木层的自然更新也会产生一定的影响（杜忠等，2016）。乔木层林分结 构在一定程度上决定了林内个体间的空间生态位（袁士云等，2010），对林内透 光条件、土壤基质活性及林内微环境等（燕亚飞等，20

5、14;黎芳等，2016;周晓 果等， 2017 ）有明显的影响，进而导致林下植被生长情况发生变化。因此，研究 林分结构特征对灌木层物种多样性的影响对实现生态系统健康发展具有重要意义。 林分结构包括空间结构和非空间结构。林下物种多样性影响因素的研究主要集中在 海拔、坡度、环境、土壤等（唐志尧等，2004;王顺忠等，2005;王家鸣等， 2018），关于林分结构对林下物种多样性的研究相对较少。这些研究对于了解不 同地理气候区域的物种多样性状况具有重要作用，但在特定环境条件下难以通过经 营调控这些因子达到提高多样性的目的。所以，如何在特定的环境条件和立地条件 下通过改善林分结构以提高物种多样性、促进

6、林分健康稳定生长是值得研究的重要 课题。林分结构对林下植被多样性具有明显影响(张建宇等， 2017；冯琦雅等，2018)。 朱光玉等(2018)对栎类林的研究结果表明，混交度对灌木层物种多样性的影响 力最大，但其他林分类型的林分结构对灌木层物种多样性的影响有待进一步开展。 马尾松林是三峡库区重要的防护林之一，本研究以该地区马尾松林为研究对象，在 全面分析林分结构特征的基础上，通过因子分析、典型相关分析等耦合林分结构综 合特征与灌木层物种多样性的关系，以期在经营管理中，为提高库区马尾松林灌木 层物种多样性、维持生态系统稳定性提供理论依据。1 研究区概况秭归县(1100-11018E,3038-3

7、11TN )位于湖北省宜昌市，地形整体上 西南部高于东北部，属于典型的长江三峡山地地貌。气候条件属于亚热带大陆季风 气候，年均温18 C左右，雨量充沛，海拔高度差达1400 m以上，气温随海拔梯 度变化明显。本研究中的马尾松林为1970年代的飞播林，平均海拔1200 m左右， 无人为经营，极少受到人为干扰，植被都以自然演替为主，土壤类型为黄壤。乔木 树种主要有马尾松、亮叶桦(Betula luminifera )、檫木(Sassafras tzumu )、 木姜子(Litsea pungens )、漆(Toxicodendron vernicifluum )等，林分平均 胸径15.23 cm，

8、平均树高11.51 m，林分密度950 - 1600 planthm-2。灌木层 物种主要有宜昌荚蒾(Viburnum erosum )、高粱泡(Rubus lambertianus )、 野蔷薇(Rosa multiflora )等。2 研究方法2.1 样地设置与调查 样地设置在秭归县国有林场，在同一坡面以样带形式选取立地条件相对一致的 15 块马尾松标准样地，面积20 mx20 m,样地都位于中坡位。对样地内所有胸径 5 cm的树挂牌，对所有挂牌的树进行每木检尺，测量指标包括胸径（DBH）、 树高（H）、枝下高、东南西北方向冠幅等，设置5 m缓冲区，用全站仪定位， 获取结构单元中对象木与邻

9、近木的信息。在每个样地的4个顶角及中心处设置5 个2mx2 m的灌木样方，记录样方内属于灌木层的所有物种的名称、盖度、频度 等信息。2.2 林层划分林层划分采用树冠光竞争高度（CCH ）法，该方法考虑了不同高度层树叶在光合 作用中贡献差异的生态学现象，根据树高和冠长确定林层的高度。本研究将林分分 为上林层、中林层和下林层，计算公式为CCH二aCL+UB,其中，CL为冠长，即 树高与枝下高的差值，UB为枝下高，a为截止系数。以高度低于下林层最低高度 的林木个体数不超过林木个体总数的10%为准，结合本文数据实际情况，截止系 数取0.4。具体方法参照郑景明等（2007）和陈科屹等（2017）。2.3

10、 结构特征参数及多样性指数的计算本研究选取胸径、树高、枝下高、冠长、冠幅、林层指数、开敞度、竞争指数、空 间分布密度指数、角尺度、大小比数、混交度及林分密度等 13个参数反映林分结 构综合特征。部分参数计算方法及含义如下：（1）冠幅（CW） 本研究中的冠幅取东、南、西、北4个方向的平均值（符利勇等， 2013） ，计算 公式为：式中,CW东、CW南、CW西和CW北分别代表东、南、西、北方向冠幅长度。（2 ）林层指数（S ）林层指数表达的是林层的垂直结构，其反映了空间结构单元中林层的垂直分布格局及林层多样性（Lv et al. , 2012 ;曹小玉等，2015 ），计算公式为：式中，zi表示以

11、对象木i为空间结构单元中的林层个数；当对象木i与相邻木j属 于同一层时，Sij=0 ；当对象木i与相邻木j不属于同一层时，Sij=1。其值越大代 表林层的垂直结构越复杂。（3）开敞度（K）。开敞度是用来反映林内的透光条件的指标（汪平等，2013），计算公式为：式中，Dij为对象木i与相邻木j之间的水平距离，Hij为对应相邻木j的树高。（4 ）林木竞争指数（CI） 竞争指数反映了林木个体所承受的来自相邻木的竞争压力大小（李建军等， 2010） ，计算公式为：式中，Cli为对象木i的竞争指数，di和dj分别为对象木i和相邻木j的胸径，Lij 代表对象木i与相邻木j之间的距离，LD代表空间结构单元中

12、4株相邻木与对象 木的距离平均值，n为竞争木株树。（5）水平分布格局林分水平分布格局选取林木空间分布密度（D）（Hui et al.，2004 ）和角尺度（W）（惠刚盈等，2004 ）两个指标，其计算公式分别为：式中，ri为对象木i包含相邻n株树时的最小半径，rmax为林分中两株相邻树木 之间的最大距离；式中，n为相邻木株数，当对象木i与相邻木j的夹角小于标准角时，则Zij=1, 否则Zij=O。（6 ）大小比数（U ） 大小比数反映的是林木大小的分化程度（惠刚盈等，1999），本研究以胸径大小 来表征，计算公式为：式中，n为相邻木株数，当对象木i比相邻木j小时，则Kij = 1 ,否则Kij

13、=O。（7）混交度（M） 混交度用来反映树种隔离程度的指标（唐万鹏等，2015），其值越大表明树种混 交程度越高。计算公式为：式中，n为相邻木株数，若对象木i与相邻木j的树种不同，则mij=1,否则 mij=0。（8）灌木层多样性指标 在耦合林分结构与物种多样性的关系时，应综合选取与多样性相关的各种指标体系（张建宇等，2017 ），本研究选取Margalef丰富度指数（R）（袁士云等， 2010）、丰富度（Richness）、Shannon-Wiener 指数（H，）、Pielo u 均匀度 指数（J ）和Simpson多样性指数（S）（Hu et al. , 2016 ）等5个指标反映灌 木

14、层多样性情况。计算公式如下：式中，S为物种数，N为所有物种的个体总数，Pi为物种i的重要值。 不同林层林分结构特征见表1。2.4 因子分析及典型相关分析在因子分析中，采用KMO （ Kaiser MeyerOlkin ）量度对主成分分析的结果进行 检验，一般认为KMO值要大于0.5，如果小于0.5则需要剔除异常变量或引入与 异常变量有关联的其它变量(高慧璇，1997)。主成分分析结果中，特征值接近 0的主成分具有多重共线性，此时可将特征值接近0的主成分中权重较大的变量剔 除，再用剩余变量重新进行主成分分析，直到KMO值满足一般条件(袁志发等， 2002)。表1不同林层林分结构特征统计Table

15、 1 Statistical of forest structure characteristics for different layersUpper layer Middle layer Lower layer Variables Mean value Range Standard deviation Variable coefficient/%Mean value Range Standard deviation Variable coefficient/%Mean value Range Standard deviation Variable coefficient/%DBH/cm 2

16、4.42 18.18-30.20 2.90 11.86 14.538.84-19.05 3.26 22.46 8.03 5.87-10.75 1.40 17.41 Tree height (H)/m 16.12 13.56-17.97 1.36 8.42 11.178.79-13.81 1.40 12.56 7.52 5.83-9.38 0.88 11.76 Under crown height(UB)/m 7.74 6.25-8.87 0.78 10.10 6.293.93-8.38 1.22 19.35 4.25 2.775.82 0.90 21.07 Crown length (CL)/

17、m 8.41 6.49-10.12 1.09 12.96 4.904.106.11 0.61 12.37 3.44 2.25-4.79 0.67 19.54 Crown width (CW)/m 2.77 2.263.27 0.27 9.75 2.021.38-2.55 0.33 16.42 1.75 1.16-2.25 0.32 18.28 Canopy layer index(S) 0.39 0.19-0.53 0.11 29.11 0.080.00-0.17 0.06 70.74 0.01 0.000.04 0.01 171.27 Open degree (K) 0.24 0.18-0.

18、32 0.04 16.78 0.260.20-0.33 0.05 18.71 0.25 0.18-0.30 0.04 14.67 Competition index(CI) 0.75 0.33-0.98 0.15 20.61 1.541.10-2.20 0.35 22.38 2.64 2.04-5.54 0.83 31.23 Spatial distribution density index (D) 0.80 0.76-0.85 0.03 3.47 0.810.70-0.87 0.05 6.18 0.82 0.78-0.90 0.03 3.80 Uniform angle index(W)

19、0.52 0.39-0.67 0.07 12.83 0.540.35-0.75 0.10 18.10 0.53 0.45-0.67 0.07 12.42 Neighborhood comparison (U) 0.20 0.04-0.33 0.07 35.30 0.450.27-0.65 0.10 22.65 0.660.47-0.83 0.08 12.70 Mixed degree (M) 0.62 0.33-0.79 0.12 19.94 0.660.330.81 0.14 20.95 0.69 0.50-0.84 0.11 15.73 Stem density (SD)/(treehm-

20、2) 410 150-675 114.49 27.92 345 100-725 174.80 50.67 433 75-675 127 29 典型相关分析可用来判断两组变量之间的相关性，即使变量组 ui=ai1x1+ai2x2+ai3x3和vi=bi1y1+bi2y2在线性组合ui和vi的方差为1的约 束条件下，让其相关系数达到最大(朱光玉等，2018)。本研究通过因子分析构 建反映林分结构特征和灌木层多样性的因子，提取具有最大相关性的线性组合并得 到典型相关系数。Pearson简单相关分析、主成分分析和典型相关分析在SPSS 22.0中完成。3 结果与分析3.1 灌木层群落特征灌木层Shan

21、non-Wiener指数、Margalef指数、Pielou扌旨数、Simpson指数和 丰富度的平均值分别为 2.04、 1.67、 0.83、 5.93和 12.07。共有 44个物种，其 中，耐阴物种、阳性物种和中性物种的个数分别为6个、13个和25个。在耐阴 物种中的乔木种有1个，为檫木，在阳性物种中的乔木种有化香树(Platycarya strobilacea )、亮叶桦和山樱花(Cerasus serrulata )，在中性物种中的乔木种 有梾木(Swida macrophylla )、茅栗(Castanea seguinii)、木姜子、漆、青 麸杨(Rhus potaninii)

22、和香椿(Toona sinensis )。各样地灌木层物种重要值排 名前三的物种主要有木姜子、高粱泡(Rubus lambertianus )、宜昌荚蒾(Viburnum erosum )、桦叶荚蒾(Viburnum betulifolium )、麻核枸子(Cotoneaster foveolatus )、檫木、亮叶桦、漆、湖北算盘子(Glochidion wilsonii)、蜡莲绣球(Hydrangea strigosa )和中国旌节花(Stachyurus chinensis)等。3.2 林分结构与多样性指数相关性分析Pearson相关性检验结果表明(表2),林分总体结构特征各参数与多样性

23、指数均 无显著相关性(P 0.05 )。在中林层和下林层，Simpson指数分别与树高和胸径 呈显著正相关(P v 0.05 ), Pielou指数则在下林层中与林木竞争指数负相关(P v 0.05 )，与林分密度正相关(P v 0.01)，而其他多样性指数与结构特征则无显 著性关系。但在上林层，各多样性指数均与林分结构表现出显著或极显著关系， Shannon-Wiener指数随着上林层树高及枝下高的增大而增大(P v 0.01), Pielou指数和Simpson指数与上林层胸径、树高、枝下高、冠长呈显著或极显著 正相关；上林层的垂直结构越复杂，对应的Pielou指数越高；林木混交度越大，

24、对应的Simpson指数越小；林木大小分化程度越高，混交程度越低，Margalef 指数与丰富度越大，Margalef指数与枝下高表现出显著正相关，丰富度R与其的 相关性则未通过显著性检验。综上，上林层是影响灌木层多样性的关键层次，因此， 在后续分析中都以上林层的结构特征参数为基础进行统计分析。表2林分结构特征与多样性指数相关分析Table 2 Correlation analysis between forest structure characteristics and diversity index*: Pv0.01; *: Pv 0.05Variables All layer Uppe

25、r layer Middle layer Lower layer ShannonWiener Margalef Pielou Simpson Richness Shannon-Wiener Margalef Pielou Simpson Richness Shannon-Wiener MargalefPielouSimpson Richness Shannon-Wiener Margalef PielouSimpson Richness DBH/cm 0.38 0.27 0.38 0.50 0.27 0.31 -0.08 0.61* 0.52* -0.08 0.28 0.07 0.39 0.3

26、9 0.07 0.49 0.42 0.44 0.56* 0.42 Tree height(H) 0.44 0.45 0.29 0.45 0.45 0.72* 0.37 0.81* 0.81* 0.37 0.51 0.47 0.42 0.53* 0.47 0.30 0.22 0.34 0.42 0.22 Under crown height (UB) 0.39 0.39 0.30 0.40 0.40 0.66* 0.51* 0.57* 0.60* 0.51 0.42 0.33 0.44 0.42 0.33 0.13 0.12 0.18 0.28 0.12 Crown length(CL) 0.3

27、4 0.27 0.26 0.41 0.28 0.42 0.07 0.62* 0.59* 0.07 0.33 0.42 0.09 0.38 0.42 0.24 0.060.31 0.25 0.06 Crown width(CW) -0.19 -0.05 -0.31 -0.12 -0.05 -0.36 -0.43 - 0.13 -0.13 -0.43 -0.18 -0.05 -0.32 -0.12 -0.05 0.00 0.11 -0.19 -0.03 0.11 Canopy layer index (S) 0.36 0.07 0.50 0.41 0.07 0.29 -0.16 0.60* 0.4

28、0 -0.16 0.43 0.28 0.44 0.39 0.28 0.29 0.11 0.38 0.31 0.11 Open degree(K)-0.06 -0.24 0.12 -0.10 -0.24 -0.30 -0.30 -0.20 -0.32 -0.30 0.03 -0.15 0.16 -0.06 -0.15 - 0.13 -0.22 -0.03 -0.16 -0.22 Competition index (CI) -0.11 -0.25 0.08 0.06 - 0.25 0.20 0.34 0.04 0.17 0.35 -0.14 -0.20 -0.02 -0.02 -0.19 -0.

29、40 -0.12 -0.58* -0.38 -0.12 Spatial distribution density index(D)-0.10 -0.24 0.04 0.02 -0.240.05 -0.43 0.35 0.14 -0.44 -0.14 0.00 -0.24 -0.10 0.000.00 0.03 -0.11 -0.03 0.03 Uniform angle index (W) -0.33 -0.48 -0.08 -0.20 -0.48 -0.08 -0.28 0.09 -0.02 -0.28 -0.15 -0.27 0.06 0.00 -0.27 -0.39 -0.32 -0.3

30、6 -0.38 -0.32 Neighborhood comparison(U)0.25 0.38 0.07 0.18 0.380.48 0.55* 0.25 0.33 0.55* 0.06 0.38 -0.24 -0.01 0.38-0.21 -0.04 -0.33 -0.21 -0.04 Mixed degree(M) -0.23 -0.38 0.02 -0.29 -0.38 -0.49 -0.53* -0.33 -0.53*-0.53* -0.09 -0.17 0.06 -0.22 -0.17 0.05 -0.08 0.22 0.02 -0.08 Stem density(SD) 0.0

31、8 - 0.04 0.20 0.15 -0.04 0.11 0.16 0.03 0.07 0.16 -0.16 -0.04-0.23 -0.14 -0.04 0.26 -0.17 0.65* 0.39 -0.173.3 多元逐步回归结果 在相关分析的基础上，以上林层结构特征参数为自变量，以多样性指数为因变量进 行多元逐步回归，从而保留影响灌木层物种多样性的主要解释变量，由于 Margalef丰富度指数和丰富度代表的意义较为接近，故多元逐步回归中仅以 Margalef代表。从表3可以看出，Shannon-Wiener指数、Pielou指数及 Simpson指数的回归方程中，都剔除了树高而保留了其

32、它所有变量，Margalef指 数回归方程则仅剔除了大小比数。这说明灌木层物种多样性受到了多种林分结构特 征的综合影响。3.4 因子分析对林分结构变量和多样性指数分别做主成分分析。采用KMO量度筛选变量，剔除 具有多重共线性的变量，以达到多指标降维的目的。对 13个林分结构变量做因子 分析，其KMO=0.282，此时删除特征值为0.000的成分中权重最大的树高变量， 再以相同方式对剩下的 12 个变量做因子分析，此时删除特征值=0.010（接近于 0） 的成分中权重最大的枝下高变量，再做因子分析，此时KMO=0.505，大于一般 基本要求的 0.5 ，且最小的主成分特征值为0.19 ，不可视为

33、接近于0。因此，采用 保留的11个变量进行主成分分析，选取特征值大于1的主成分，最终保留3个主 成分，其旋转后因子载荷表见表 4。对反映多样性的 5个指数做因子分析， KMO=0.721 ，符合基本要求，保留了2个主成分，其旋转后因子载荷表见表5。 从表4可以看出，上林层林分结构特征变量可分为3类因子，在第1类因子中， 开敞度、林木竞争指数、大小比数、混交度及林分密度的载荷绝对值较大，这些变 量分别反映了林木的透光条件、林木承受的竞争压力、林木大小分化程度、树种隔 离程度和林分的稀疏状况，这些变量基本都是林木通过自然竞争后的结果，反映了 林分中林木的竞争状态，因此，可将第1个类因子定义为林分竞

34、争状况因子。第2 个主成分中则是林木空间密度、角尺度以及冠幅的载荷值较大，由于前两个变量反 映了林分水平空间分布格局，因此可将其定义为水平结构及冠幅因子。第3个主 成分中，胸径、冠长及反映林层垂直结构的林层指数载荷值较大，因此，将该主成 分定义为垂直结构及生长因子。从表5可以看出，多样性变量指标可分为2类因 子，即丰富度因子和多样性因子。表3多样性指数与上林层结构特征的多元逐步回归Table 3 The multiplestepwise regression model with diversity index and the structural characteristics of the

35、 upper layerDiversity index Regression equationShannon-Wiener y=-0.231DBH+0.340UB-0.481CL-0.314CW-0.145S- 0.055K+0.134CI-0.139D-0.135W+0.274U-0.259M+0.042SD (R2=0.474, P=0.003)Margalef y=0.087DBH+0.209H+0.333UB+0.020CL- 0.297CW+0.056S-0.022K-0.304CI-0.308D-0.052W-0.305M-0.370SD (R2=0.246, P=0.035)Pi

36、elou y=0.168DBH+0.097UB-0.122CL- 0.077CW+0.234S+0.093K-0.038CI+0.258D+0.031W-0.020U-0.034M- 0.041SD (R2=0.623, P=0.000)Simpson y=0.014DBH+0.155UB-0.222CL- 0.081CW-0.062S-0.044K+0.090CI+0.045D-0.075W+0.068U- 0.272M+0.001SD (R2=0.630, P=0.000)表 4 旋转后林分结构因子载荷 Table 4 The loading values of forest struct

37、ure factors after rotationMeaning of factors Forest competition status factor Horizontal structure and canopy width factor Vertical structure and growth factor DBH 0.05 0.53 0.80 Crown length -0.19 -0.02 0.91 Crown width 0.04 0.71 -0.05 Canopy layer index 0.55 0.08 0.74 Open degree 0.80 -0.29 -0.09

38、Competition index -0.78 -0.40 -0.20 Spatial distribution density index 0.08 0.77 0.37 Uniform angle index 0.38 0.72 0.10 Neighborhood comparison - 0.82 -0.34 -0.02 Mixed degree 0.86 0.14 -0.07 Stem density -0.62 -0.46 - 0.23表5旋转后多样性指数因子载荷Table 5 The loading values of diversity index factors after ro

39、tationMeaning of factors Richness factor Diversity factor Shannon-Wiener 0.69 0.72 Margalef 0.98 0.21 Pielou 0.08 0.99 Simpson 0.37 0.91 Richness 0.98 0.213.5 典型相关分析 根据以上因子分析可得到每个样地的3个林分结构因子和2个灌木层多样性因子 的标准化因子得分，并以此为基础进行典型相关分析，结果共得到2组典型相关结构（表6）,其中，第1典型相关结构通过显著性检验（P=0.014 ），典型相关系数为 0.822，这表明上林层林分结构特征与

40、灌木层多样性总体上有较强的相关性 从表6可以看出，反映林分结构特征的3个因子中，垂直结构及生长因子的载荷 最大，为 0.962，表明该因子是影响灌木层多样性的主要因子，且灌木层多样性随 着该因子的增大而变丰富。其次是林分竞争因子，其载荷为-0.265，说明灌木层 多样性随着该因子的增大而变小。而灌木层多样性对水平结构及冠幅因子的敏感性 则较小。4 讨论 林下植被生长情况是多种因素共同作用的结果，其多样性很大程度上受到了林分演 替阶段、过程及空间异质性等对生境改造的影响（鱼腾飞等，2011）。乔木层作 为森林生态系统的重要组成部分，其生长状况可对林分空间异质性（黎芳等，2016 ）、土壤条件及微

41、气候（Uliassi et al. , 2002 ;Rodriguezcalcerrada et al.,2008 ）等产生影响，从而导致林下物种多样性发生变化（ Lhotka et al.，2008；Parker et al. , 2008 ）。Pearson相关分析表明，林分总体结构特征各变量均与灌木层多样性无显著性关系,分林层后发现中林层和下林层的结构特征与灌木层多 样性有较微弱的关系,而上林层中的结构变量则与灌木层多样性有较丰富的显著相 关性,这在一定程度上表明了林下植被物种多样性影响因素的复杂性（刘冠成等, 2018 ） 。上林层对整个林分的树种组成有关键作用（玉宝等, 2015 ）

42、 ,在本研究 中,上林层的蓄积量占整个林分蓄积的74.73% ,该层次多由林分中径级较大的林 木组成,在自然状态下,大径级林木（喻理飞等,2002）的生物量与凋落物量相 对较多（王莉莉等,2014）,对土壤及整个生态系统的养分循环有重要影响,而 土壤条件又与物种多样性有密切关系（王顺忠等,2005）,且上层林木对林内光 照分布及物种生态位有明显影响（欧芷阳等,2017）,这在一定程度上解释了灌 木层物种多样性与上林层的结构特征的关系更为紧密的原因。采用多元逐步回归的主要目的是为了排除干扰变量而保留解释变量（Hu et al., 2018 ），本研究中，保留的12个结构特征变量分别解释了 Pie

43、lou指数、 Simpson 扌旨数和 Shannon-Wiener 指数的 62.3%、63.0%和 47.4%，而仅解释 了 Margalef指数的24.6%。朱光玉等（2018 ）采用混交度、聚集指数和开敞度 等3个林分结构特征变量与多样性指数作多元逐步回归的结果显示，这3个变量 仅解释了 Auclair&Goff多样性指数的23.3%和Pielou指数的16.3% , Margalef 则未能通过回归方程的显著性检验。该结果与本研究结果对比表明，灌木层物种多 样性是林分结构多种特征共同作用的结果。为探索林分结构的变量组合与灌木层多样性之间的关系，本文采用了典型相关分析。 结果表明上林层

44、林分结构特征和灌木层多样性之间存在显著的相关关系，其中垂直 结构及生长因子的影响最大，即胸径越大、冠长越长，林层的垂直结构越丰富，灌 木层多样性就越高。这主要是因为处于上林层的大树能够获得更好的光照条件，并 且能够产生更多的凋落物，与土壤间的养分流动更多，从而提高了土壤有机质含量， 而土壤孔隙度、蓄水量及养分等会随着土壤有机质的提高而增大（王莉莉等， 2014），这些都是影响灌木层植被生长的潜在因素。垂直结构越复杂，林下的光 照可能会更强（Pelt et al. , 2000 ），从而提高了物种多样性（蒋政等，2015 ）。 林分竞争因子整体与灌木层多样性呈负相关，这可能是由于林分竞争越激烈，

45、对于 有限资源的利用就越充分，灌木的生长由此受到影响。有研究认为，林分混交度越 大，其对应的灌木层多样性越高（朱光玉等，2018），然而在本研究中，多元回 归方程和典型相关分析结果都显示灌木层多样性随混交度的增大而减小，这与林分 类型有较大的关系，马尾松林混交度越低，其透光性越好，越利于林下层植被的生 长（刘延滨等，2012），这与张建宇等（2017）对落叶松林的研究结果类似。此 外，还与灌木层树种特性有关（何友均等，2013）。表6典型相关分析结果及各因子的典型载荷Table 6 The results of canonical correlation analysis and typica

46、l loads of each factorTypical structure Canonical correlation coefficient Significance level Forest competition status factor Horizontal structure and canopy width factor Vertical structure and growth factor Richness factor Diversity factor No.1 canonical correlation structure 0.822 0.014 -0.265 0.0

47、61 0.962 0.013 1.000 No.2 canonical correlation structure 0.523 0.173 -0.666 -0.733 -0.137 1.000 0.012 曹梦等（2018）在研究马尾松飞播林林下物种多样性与环境因子关系时，选取了 包括海拔、株数密度、石砾含量、平均胸径等 20个生物因子和非生物因子，以前 向选择和蒙特卡洛检验法进行关系耦合，最后得到了9个影响因子，解释了物种 多样性变异的70.9%，并得出生物因子是影响多样性的主要因子的结论，这与本 研究的结果具有一定的相似性。典型相关分析对样本数量和变量个数有一定的要求， 如果变量过多而样本

48、数量不足时，其分析效果会受到较大的影响。因此本研究通过 因子分析将变量进行了降维分类，取得了较好的效果。5 结论 在马尾松林中，上林层是影响灌木层物种多样性的主要层次，灌木层物种多样性受 到该层次多种林分结构特征的共同影响。上林层林分结构可分为3类，其中，垂 直结构及生长因子对灌木层物种多样性的影响最大，表现为正相关关系。在经营管 理中，可通过调控上林层的胸径、冠长及冠层垂直结构以达到提高林分灌木层物种 多样性的目的。相关文献】 BURRASCANO S, SABATINI F M, BLASI C, 2011. Testing indicators of sustainable forest

49、 management on understorey composition and diversity in southern Italy through variation partitioning J. Plant Ecology, 212(5): 829-841.HU W J, WANG P C, ZHANG Y, et al., 2016. Effect of gap disturbances on soil properties and understory plant diversity in a Pinus massoniana plantaion in Hubei, cent

50、ral China J. The Journal of Animal & Plant Sciences, 26(4): 988-1001.HU W J, WANG X Y, CUI H X, et al., 2018. Distribution characteristics and factors influencing the ecosystem in western Hubei J. Acta Ecologica Sinica, 38(2): 165-172.HUI G Y, VON G K, HU Y B, et al., 2004. Characterizing forest spa

51、tial distribution pattern with the mean value of uniform angle index J.Acta Ecologica Sinica, 24(6): 1225- 1229.LHOTKA J M, LOEWENSTEIN E F, 2008. Influence of canopy structure on the survival and growth of underplanted seedlings J. New Forests,35(1): 89-104.LV Y, ZANG H, WAN X J, et al., 2012. Stor

52、ey structure study of Cyclobalanopsis myrsinaefolia mixed stand based on storey index J.Forest Resources Management, 3: 8184.PARKER W C, DEY D C, 2008. Influence of overstory density on ecophysiology of red oak (Quercus rubra) and sugar maple (Acer saccharum) seedlings in central Ontario shelter woo

53、ds J. Tree Physiology, 28(5): 797-804.PELT R V, FRANKLIN J F, 2000. Influence of canopy structure on the understory environment in tall, old-growth, conifer forests J.Canadian Journal of Forest Research, 30: 1231-1245.RODRfGUEZCALCERRADA J, REGNERI S M, ALONSO J, et al., 2008.1nflue nee of overstory

54、 density on understory light, soil moisture, and survival of two underplanted oak species in a Mediterra nean montane Scots pine forest J. I nvestigaci6 n Agraria Sistemas Y Reeursos Forestales, 17(1): 31-38.ULIASSI D D, RUESS R W, 2002. Limitations to symbiotic nitrogen fixation in primary successi

55、on on the Tanana River floodplain J.Ecology, 83(1): 88-103.曹梦，潘萍，欧阳勋志,等,2018.飞播马尾松林林下植被组成、多样性及其与环境因子的关系J.生 态学杂志, 37(1): 1-8.CAO M, PAN P, OUYANG X Z, et al., 2018. Relationships between the composition and diversity of understory vegetation and environmental factors in aerially seeded Pinus massonia

56、na plantationsJ. Chinese Journal of Ecology, 37(1): 1-8.曹小玉，李际平,封尧,等,2015.杉木生态公益林林分空间结构分析及评价J.林业科学,51(7): 37- 48.CAO X Y, LI J P, FENG R, et al., 2015. Analysis and evaluation of the stand spatial structure of Cunninghamia lanceolata ecological forestJ. Scientia Silvae Sinicae, 51(7): 37-48. 陈科屹，张会儒

57、，雷相东,等,2017.云冷杉过伐林垂直结构特征分析J.林业科学研究,30(3): 450-459.CHEN K Y, ZHANG H R, LEI X D, et al., 2017. Analysis of vertical structure characteristics for Spruce-Fir over-cutting forest J. Forest Research, 30(3): 450-459.褚建民卢琦，崔向慧,等,2007.人工林林下植被多样性研究进展J.世界林业研究,20(3): 9-13.CHU J M, LU Q, CUI X H, et al., 2007.

58、 Review on species diversity of undergrowth vegetation in plantation ecosystem J. World Forestry Research, 20(3): 9-13.杜忠,蔡小虎，包维楷,等,2016.林下层植被对上层乔木的影响研究综述几应用生态学报,27(3): 963-972.DU Z, CAI X H, BAO W K, et al., 2016. Understory effects on overstory trees: A review J. Chinese Journal of Applied Ecology

59、,27(3): 963-972.方升佐，田野,2012.人工林生态系统生物多样性与生产力的关系J.南京林业大学学报：自然科学版, 36(4): 1-6.FANG S Z, TIAN Y, 2012. The relationship between biodiversity and productivity in the artificial plantation ecosystem J. Journal of Nanjing Forestry University (Natural Science Edition), 36(4): 1-6.冯琦雅，陈超凡，覃林，等,2018.不同经营模式对蒙

60、古栋天然次生林林分结构和植物多样性的影响J.林业科学, 54(1): 12-21.FENG Q Y, CHEN C F, QIN L, et al., 2018. Effects of different management models on stand structure and plant diversity of natural secondary forests of Quercus mongolica J. Scientia Silvae Sinicae,54(1): 12-21.符利勇，孙华，张会儒,等,2013.不同郁闭度下胸高直径对杉木冠幅特征因子的影响J.生态学报， 33

61、(8): 2434-2443.FU L Y, SUN H, ZHANG H R, et al., 2013. Effects of diameter at breast height on crown characteristics of Chinese Fir under different canopy density conditions J. Acta Ecologica Sinica, 33(8):2434-2443.高慧璇,1977. SAS系统SAS/STAT软件使用手册M.北京:中国统计出版社:160-524.GAO H X, 1997. SAS system SAS/STAT

62、 software user manual M.Beijing: China Statistics Press: 160-524.何友均, 梁星云, 覃林, 等, 2013. 南亚热带马尾松红椎人工林群落结构、物种多样性及基于自然的森 林经营J林业科学,49(4): 24-33.HE Y J, LIANG X Y, QIN L, et al., 2013. Community Structure，Species diversity of Pinus massoniana and Castanopsis hystrix plantation and the nature-based forest

63、 management in the southern subtropical China J. Scientia Silvae Sinicae, 49(4): 24-33.惠刚盈,GADOW K V,胡艳波,2004.林分空间结构参数角尺度的标准角选择J.林业科学研究， 17(6): 687-692.HUI G Y, GADOW K V, HU Y B, 2004. The optimum standard angle of the uniform angle index J. Forest Research, 17(6): 687-692.惠刚盈,GADOW K V, MATTHIAS A

64、, 1999. 一个新的林分空间结构参数:大小比数J.林业科学研 究, 12(1): 1-6.HUI G Y, GADOW K V, MATTHIAS A, 1999. A new parameter for stand spatial structureneighbourhood comparison J. Forest Research,12(1): 1-6.蒋政,丁雨龙，宋东杰,等,2015.安徽横山植物群落的组成、结构与多样性分析J.植物资源与环境 学报, 24(3): 99-106.JIANG Z, DING Y L, SONG D J, et al., 2015. Analyses

65、 on composition, structure and diversity of plant community of Hengshan Mountain in Anhui province J. Journal of Plant Resources and Environment, 24(3): 99-106.黎芳,潘萍,宁金魁，等,2016.飞播马尾松林林分空间结构对林下植被多样性的影响J东北林业大 学学报, 44(11): 31-35.LI F, PAN P, NING J K, et al., 2016. Effects of stand spatial structure on understory vegetait on diversity of aerial seeding Pinus massoniana plantations J. Journal of Northeast Forestry University, 44(11): 31-35.李建军，李际平，刘素青,等