-
杉木（Cunninghamia lanceolata (Lamb.) Hook.）是我国南方主要造林树种。第九次全国森林资源清查结果显示，我国杉木人工林面积达到9.90 × 106hm2，蓄积量达7.55亿m3，分别占全国主要优势人工林树种的1/4及1/3[1]。然而，由于当前杉木人工林存在过纯、过密等，导致林分结构单一、林地生产力下降和林分结构不稳定[2-3]。同时，随着社会经济的发展，人们对杉木的木材结构需求发生了变化，逐渐由小径材转为大径材。因此，杉木纯林不能满足人们的需求，甚至影响到社会经济的可持续发展[4]。此外，有研究表明，营建混交林可以改善林内空间结构，增加林内光照的同时减少林内竞争[5]，有利于林木生长和提高生物多样性[6]，提高森林的经济和生态效益。因此，如何通过调整优化林分结构，提升杉木人工林的质量已成为亟待解决的问题。
在营造异龄复层混交林的研究上，国外主要集中在阔叶树种[7-10]，而国内重点是针叶树种。已有研究表明，间伐套种处理改变了林分密度以及生长空间，降低了林木竞争，有利于林分生长、林分生物量和生长力增加[11-13]。李婷婷[14]对杉木人工林进行间伐补植改造，发现林下补植阔叶树种可以显著提高林木单木材积和林分蓄积的年生长量，并且大叶栎、红椎、格木、灰木莲适合与杉木进行混交种植。孙冬婧等[15]在对红椎、米老排、大叶栎、润楠与杉木混交林的生长与生态效应研究中得出，混交林有利于林分生长，并形成丰富的物种多样性，无论在经济方面还是生态方面都比杉木纯林更加优化。欧建德[16]研究了福建南方红豆杉林下套种模式对人工林生长的影响，结果表明：林下套种后，林分冠幅显著大于纯林模式。陈莹莹[17]以福寿林场杉木人工中龄林为研究对象，通过在林下补植阔叶树，将杉木纯林改造为异龄针阔混交林，发现杉木的高径比、胸高形数和枝下高均随间伐强度的增强而减小。
目前，国内外对间伐补植经营的研究大多集中在间伐强度对林分结构和生长、生物多样性和林地土壤养分的影响等方面[18-21]，而对杉木人工林“间伐 + 冠下套种”后的杉木生长、干形形质特性和材种结构的量化评价不足。因此，本研究以安徽省青阳县杉木人工林为研究对象，调查林木的生长性状、空间利用能力、林分产量，通过主成分分析法综合评价“间伐 + 套种”对杉木林分质量的提升效果，为杉木人工林质量精准提升提供科学依据。
-
试验地位于安徽省池州市青阳县庙前镇三义村窑西林场（117°40′～118°05′ E，30°19′～30°50′ N），该地属于北亚热带湿润季风气候，温和多雨，年平均气温16.1℃，年平均降水量1 500 mm，年平均日照时数2 066 h，无霜期218 d；地形以丘陵为主，土壤为黄壤。试验林营建于2000年，平均林分密度为1 740株·hm−2，立地指数为18，无施肥和间伐措施。优势种为杉木；灌木植物以六月雪（Serissa japonica（Thunb.）Thunb.）、淡竹（Phyllostachys glauca McClure）、荚蒾（Viburnum dilatatum Thunb.）等为主；草本层主要有边缘鳞盖蕨（Microlepia marginata (Houtt.) C. Chr.）、点腺过路黄（Lysimachia hemsleyana Maxim. ex Oliv.）、山麦冬（Liriope spicata (Thunb.) Lour.）等。
-
2017年底，选取立地条件相近的杉木人工林地段，根据伐除木的数量设计47%（Ⅰ处理）、56%（Ⅱ处理）和65%（Ⅲ处理）3个间伐强度并设置样地。样地按照随机区组设计，Ⅰ处理、Ⅱ处理、Ⅲ处理、CK处理各设置20 m × 20 m的样地3个，同时Ⅰ处理和Ⅱ处理各增加1块对照样地，共计14块样地；为避免边界效应，实际间伐作业面积为30 m × 30 m。2018年初，在间伐处理后的样地中，用浙江楠（Phoebe chekiangensis C. B. Shang）和檫木（Sassafras tzumu（Hemsl.）Hemsl）1年生苗进行冠下套种，随机均匀排列种植，株行距为4 m × 5 m，2树种种植密度均为225株·hm−2。样地基本情况见表1。
表 1 样地基本概况
Table 1. General situation of sample plots
处理
Treatment杉木保留密度
The retained density/
(tree·hm−2)套种株数
The interplanting/
(tree·hm−2)坡向
Aspect坡度
Slope/(°)杉木平均胸径（2018）
Average DBH
(2018)/cm杉木平均树高（2018）
Average tree height
(2018)/mCK（对照组） 1 740 0 东南 15-20 13.1 10.7 Ⅰ（47%间伐 + 套种） 1 050 450 东南 15～20 15.1 11.4 Ⅱ（56%间伐 + 套种） 825 450 东 15～20 15.6 11.8 Ⅲ（65%间伐 + 套种） 675 450 东 15～20 17.2 12.8 分别在间伐套种第0年（2018年）、第2年（2020年）、第4年（2022年）对样地进行调查，对各样地的林木进行每木检尺，使用胸径尺测量杉木胸径和套种树种的地径，使用激光测高仪测量林木树高与枝下高，使用皮尺测量林木树冠的南北冠幅和东西冠幅，使用2022年的每木调查数据对林分生长进行统计分析。
-
（1）林木生长形质性状 采用胸径、树高、胸高形数、高径比、枝下高指标衡量杉木的生长形质性状。
（2）空间利用能力 空间利用能力从林木二维树冠和三维树冠两方面进行描述，二维树冠用冠幅、冠长、冠长率、冠形率来表示，三维树冠通过树冠表面积、树冠体积来表达。
（3）林分产量 林分产量从杉木单木材积、林分蓄积和材种结构3方面进行描述，材种结构用小条木、小径材、中径材、大径材出材率和用材、薪材、废材出材率来表示。
-
应用SPSS27.0中的单因素方差分析法对不同处理下杉木林分特征的差异进行分析，利用主成分分析法进行综合评价，步骤为[22]：
首先，将n种间伐套种处理的m项指标组成数据矩阵X。
$ X={\left[{X}_{ij}\right]}_{n\times m}\left(i=\mathrm{1,2},3\dots ,n;j=\mathrm{1,2},3\dots ,m\right) $
(1) 式中：Xij为第i种间伐套种处理的第j项指标的实测值。
其次，利用SPSS或Excel软件将数据标准化，Excel软件计算公式如下：
正向指标标准化：
$ {X}_{ij}^{*}={X}_{ij}/ \overline {{X}_{j}} $
(2) 逆向指标标准化：
$ {X}_{ij}^{*}= \overline {{X}_{j}}/{X}_{ij} $
(3) 式中：
$ {X}_{ij}^{*} $ 为Xij的标准化数据；$\overline {{X}_{j}}$ 为第j项指标的平均值。利用KMO检验法和Bartlett球体检验法对标准化处理后的数据进行适用性检验， KMO数值≥0.60且Sig数值＜0.05，表明各指标间关联程度较高，可以进行主成分分析。
选取主成分的特征值＞1且方差累计贡献率＞90%的前k个主成分，建立主成分与标准化指标间的关系，公式为：
$ {Y}_{p}={b}_{p1}\times {X}_{1}^{*} + {b}_{p2}\times {X}_{2}^{*} + \dots + {b}_{pm}\times {X}_{m}^{*}\left(p=\mathrm{1,2},3\dots ,k\right) $
(4) 式中：Yp为第p个主成分；bpm为第p个主成分的因子载荷。
用第p个主成分的方差贡献率和所选取的k个主成分的方差总贡献率的比值表示各个主成分的权重，将各个主成分的权重与k个主成分相结合，得到综合评价函数Y。Y的得分越高，则表明该间伐套种处理对杉木人工林的改造效果越好。公式如下：
$ Y=\sum _{k=1}^{k}\left({Y}_{p}\times \left({\lambda }_{p}/\sum _{k=1}^{k}{\lambda }_{p}\right)\right) $
(5) 式中：Y为不同间伐套种处理的综合得分；λp为第p个主成分的方差贡献率。
-
由图1可知，Ⅰ、Ⅱ、Ⅲ处理的胸径、树高均显著大于CK（P<0.05），与CK相比，Ⅰ、Ⅱ、Ⅲ处理的胸径分别增长了19.02%、30.01%、32.07%；树高则分别增长了10.40%、14.81%、21.60%。
图 1 不同间伐套种下林木胸径和树高生长变化
Figure 1. Changes of DBH and tree height growth under different thinning and interplanting
-
由图2可知，CK的高径比、胸高形数、枝下高均显著大于Ⅰ、Ⅱ、Ⅲ处理（P<0.05）。其中，Ⅱ处理的高径比最小，为0.73；Ⅲ处理的胸高形数最小，为0.52，显著低于Ⅰ处理；Ⅰ处理的枝下高最小，为3.85 m，显著低于Ⅱ与Ⅲ处理（P<0.05），Ⅱ与Ⅲ处理间的枝下高差异不显著。
图 2 不同间伐套种下林木干形生长变化
Figure 2. Changes of stem-form growth under different thinning and interplanting
-
对4个处理杉木的树冠结构进行分析（表2），发现与CK相比，Ⅰ、Ⅱ、Ⅲ处理的冠幅增长率分别为12.16%、31.00%、21.27%（P<0.05），其中Ⅱ处理的冠幅增长最快；冠长增长率依次为58.05%、55.58%、61.24%，Ⅲ处理的冠长增长最快（P<0.05）。
表 2 不同处理组树冠的生长变化
Table 2. Growth changes of crown of different treatments
处理
Treatment冠幅
Crown width/m冠长
Crown length/m冠长率
Crown ratio冠形率
Crown shapeCK 3.29 ± 0.80 D 5.65 ± 2.17 B 0.50 ± 0.17 C 1.78 ± 0.72 C Ⅰ 3.69 ± 0.85 C 8.93 ± 2.13 A 0.70 ± 0.13 A 2.58 ± 0.94 A Ⅱ 4.31 ± 0.63 B 8.79 ± 2.24 A 0.65 ± 0.13 B 2.07 ± 0.55 B Ⅲ 3.99 ± 0.67 A 9.11 ± 1.81 A 0.66 ± 0.11 B 2.35 ± 0.61 A 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).为充分反映不同间伐套种处理下杉木的树势和树冠立体状态的水平，对杉木的冠长率和冠形率进行了研究（表2）。间伐与阔叶树套种处理显著提高了冠长率与冠形率，其中Ⅰ、Ⅱ、Ⅲ处理的冠长率分别比CK增加了40.00%、30.00%、32.00%（P<0.05），Ⅱ与Ⅲ处理间的冠长率差异不显著；Ⅰ、Ⅱ、Ⅲ处理的冠形率显著高于CK（P<0.05），Ⅰ与Ⅲ处理间的冠形率差异不显著。
-
由表3可知，树冠表面积和树冠体积均随着试验处理强度的增强而先增大后减小。Ⅱ处理的值均最大，分别为61.99 m2、43.95 m3，其次依次为Ⅲ处理、Ⅰ处理、CK。其中，CK与Ⅰ、Ⅱ、Ⅲ处理间的树冠表面积存在显著性差异，Ⅰ与Ⅱ处理间的树冠表面积差异显著（P<0.05）；杉木人工林单木树冠体积在各处理间呈现显著性差异。杉木树冠表面积和树冠体积的变异系数表现为：CK>Ⅰ>Ⅱ>Ⅲ，可见 CK 中单木树冠表面积、树冠体积分化程度较大，而间伐与阔叶树套种处理可以降低树冠表面积和树冠体积的分化程度。
表 3 不同处理组三维树冠的生长变化
Table 3. Growth changes of three-dimensional canopy of different treatments
处理
Treatment树冠表面积均值
Average canopy surface
area/m2树冠表面积变异系数
Coefficients of variation for
canopy surface area/%树冠体积均值
Average canopy
volume/m3树冠体积变异系数
Coefficients of variation for
canopy volume/%CK 31.65C ± 15.73 49.70 17.82D ± 12.46 69.92 Ⅰ 52.91B ± 17.59 33.25 32.92C ± 17.09 51.91 Ⅱ 61.99A ± 18.74 30.23 43.95A ± 18.09 41.16 Ⅲ 58.94AB ± 15.71 26.65 39.15B ± 14.98 38.26 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05). -
由表4可知，随着试验处理强度的增强，杉木的单木材积呈现逐渐增长的趋势，杉木的林分蓄积则呈现逐渐下降的趋势。Ⅲ处理的单木材积最大，为0.22 m3；其次为Ⅱ处理，为0.21 m3；CK最小，为0.12 m3。杉木的林分蓄积则表现为CK>Ⅰ>Ⅱ>Ⅲ，与CK相比，Ⅰ、Ⅱ、Ⅲ处理的林分蓄积分别下降了27.50%、28.15%、36.94%（P<0.05）。
表 4 间伐套种后杉木单木材积和林分蓄积的生长变化
Table 4. Growth changes of single tree and stand volume of fir plantation after thinning and interplanting
指标
Index处理CK
Treatment CK处理Ⅰ
Treatment Ⅰ处理Ⅱ
Treatment Ⅱ处理Ⅲ
Treatment Ⅲ单木材积 Single tree volume/m3 0.12 ± 0.08 C 0.18 ± 0.10 B 0.21 ± 0.09 AB 0.22 ± 0.10 A 林分蓄积 Stand volume/(m3·hm−2) 226.12 ± 37.11 A 163.94 ± 22.84 BC 162.47 ± 33.27 C 142.58 ± 18.26 D 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05). -
由图3可知，随着试验处理强度的增强，林分的薪材、废材出材量以及总出材量逐渐减小，用材出材量则呈现先减少后增加再减少的趋势，并且CK的值均最大，Ⅲ处理的值均最小。
图 3 间伐套种后林分总出材量的结构
Figure 3. Structure of total stand output after thinning and interplanting
由图4可知，间伐套种处理降低了杉木林分内小条木、小径材的出材量，而中径材和大径材的出材量则在不同间伐套种处理下呈现不同的变化规律。其中，Ⅱ处理的中径材出材量最大，为80.87 m3·hm−2，Ⅰ处理的中径材出材量最小，为61.33 m3·hm−2；大径材出材量则表现为Ⅰ>Ⅲ>Ⅱ>CK。
图 4 间伐套种后林分用材出材量的结构
Figure 4. Structure of materials output after thinning and interplanting
对杉木用材、薪材、废材出材率的分析表明（表5），间伐与阔叶树套种处理可以提高杉木用材出材率，降低薪材和废材出材率。间伐套种处理后，杉木用材出材率表现为Ⅲ>Ⅱ>Ⅰ>CK，薪材和废材出材率则表现为CK>Ⅰ>Ⅱ>Ⅲ。其中，Ⅱ与Ⅰ、Ⅱ与Ⅲ处理间用材和薪材出材率差异不显著，Ⅱ与Ⅲ处理间废材出材率差异不显著。
表 5 间伐套种处理后不同材种出材率
Table 5. The timber assortments timber-produced rate after thinning and interplanting
处理
Treatment小条木
Small strips
of wood小径材
Small diameter
wood中径材
Medium diameter
wood大径材
Large diameter
wood用材
Materials薪材
Firewood废材
Waste woodCK 10.64 ± 19.05 A 42.87 ± 30.23 A 17.57 ± 29.38 C 0.73 ± 5.51 B 71.80 ± 6.85 C 6.66 ± 6.84 A 21.54 ± 1.93 A Ⅰ 4.23 ± 12.41 B 35.15 ± 32.08 B 29.16 ± 33.48 B 6.36 ± 20.05 A 74.90 ± 4.05 B 4.23 ± 2.82 BC 20.87 ± 1.29 C Ⅱ 1.24 ± 7.77 B 26.77 ± 33.55 BC 42.63 ± 34.42 A 5.53 ± 16.05 A 76.16 ± 4.05 AB 3.44 ± 4.00 CD 20.40 ± 0.73 BD Ⅲ 0.97 ± 6.92 B 25.09 ± 31.26 C 42.58 ± 32.99 A 8.02 ± 20.71 A 76.66 ± 2.68 A 3.03 ± 1.95 D 20.31 ± 1.68 B 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).间伐与阔叶树套种处理对各材种出材率的影响不同，其中，小条木、小径材的出材率呈下降趋势，与CK相比，Ⅰ、Ⅱ与Ⅲ处理的小条木出材率分别降低了60.24%、88.35%、90.88%（P<0.05）；中径材和大径材出材率则均呈升高趋势，间伐套种处理林分的中径材和大径材出材率显著高于未处理林分（P<0.05）。
-
为了消除各指标间单位和量纲的影响，对各指标（X1：胸径、X2：树高、X3：枝下高、X4：高径比、X5：胸高形数、X6：冠幅、X7：冠长率、X8：冠形率、X9：树冠表面积、X10：树冠体积、X11：小条木出材率、X12：小径材出材率、X13：中径材出材率、X14：大径材出材率、X15：薪材出材率、X16：废材出材率、X17：用材出材率、X18：单木材积、X19：林分蓄积）进行了标准化处理，并对标准化处理后的数据进行适用性检验，发现KMO数值为0.612，Sig数值为0.00，各指标间关联程度较高，可以进行主成分分析。
-
由表6可知，前6个主成分的累计贡献率高达91.228%，并且各主成分的特征值均大于1，说明这6个主成分可代表上述的19个指标。因此，提取这6个主成分，分别为Y1、Y2、Y3、Y4、Y5、Y6。
表 6 主成分提取结果
Table 6. Extraction results of principal compositions
成分
Composition初始特征值
Initial eigenvalue提取载荷平方和
Extract the sum of squares and load it总计
Total方差百分比
% of variance/%累积百分比
Accumulate/%总计
Total方差百分比
% of variance/%累积百分比
Accumulate/%1 8.831 46.479 46.479 8.831 46.479 46.479 2 3.110 16.370 62.849 3.110 16.370 62.849 3 1.791 9.428 72.276 1.791 9.428 72.276 4 1.342 7.065 79.341 1.342 7.065 79.341 5 1.199 6.311 85.652 1.199 6.311 85.652 6 1.059 5.576 91.228 1.059 5.576 91.228 -
由表7可知：在第一主成分的特征向量中，胸径、树高、胸高形数、单木材积的特征向量值较高。在第二成分的特征向量中，枝下高和冠长率的特征向量值较高；在第三主成分的特征向量中，小条木出材率和小径材出材率的特征向量值较高；在第四主成分的特征向量中，冠幅和冠形率的特征向量值较高；在第五主成分的特征向量中，大径材出材率、薪材出材率和废材出材率的特征向量值较高；在第六主成分的特征向量中，中径材出材率和大径材出材率的特征向量值较高。
表 7 主成分特征值与荷载结果
Table 7. Principal component characteristic values and load results
指标
Index成分 Composition Y1 Y2 Y3 Y4 Y5 Y6 X1 0.960 −0.147 0.125 0.136 0.095 −0.006 X2 0.890 −0.107 −0.201 0.135 0.044 −0.042 X3 0.091 −0.947 0.029 −0.011 0.013 0.009 X4 −0.663 0.204 −0.305 −0.022 −0.004 0.020 X5 −0.919 0.172 0.209 −0.135 0.009 0.039 X6 0.722 −0.008 0.103 −0.594 −0.120 0.189 X7 0.313 0.910 −0.140 0.078 −0.002 −0.044 X8 0.205 0.690 −0.302 0.562 0.101 −0.193 X9 0.826 0.459 −0.024 −0.226 −0.028 0.086 X10 0.801 0.376 0.063 −0.334 −0.042 0.135 X11 −0.709 0.296 0.520 0.106 −0.177 0.202 X12 −0.372 −0.180 −0.862 −0.162 0.046 0.168 X13 0.703 −0.021 0.432 0.029 −0.232 −0.503 X14 0.397 0.012 0.242 0.275 0.526 0.644 X15 −0.698 0.211 0.214 −0.278 0.530 −0.223 X16 −0.618 0.142 0.192 0.242 −0.592 0.371 X17 0.833 −0.239 −0.256 0.196 −0.338 0.109 X18 0.922 −0.069 0.231 0.178 0.189 0.077 X19 −0.362 −0.453 0.103 0.358 0.126 −0.120 -
由表8可知，各处理的综合得分排名由高到低依次为：Ⅲ处理（1.186）、Ⅱ处理（1.092）、Ⅰ处理（0.671）和CK（−0.709）。其中，Ⅲ、Ⅱ和Ⅰ处理的综合得分均大于0，林分生长均优于林分平均水平；CK的综合得分小于0，林分生长低于林分平均水平。
表 8 主成分得分及综合得分
Table 8. Principal component score and comprehensive score
处理
Treatment主成分
Principal component综合得分
Comprehensive score排名
RankingY1 Y2 Y3 Y4 Y5 Y6 CK −1.322 −0.658 0.063 0.143 0.008 −0.067 −0.709 4 Ⅰ 1.050 1.181 −0.309 0.161 0.046 0.084 0.671 3 Ⅱ 2.191 0.650 0.136 −0.632 −0.108 0.105 1.092 2 Ⅲ 2.308 0.745 −0.036 −0.159 0.028 0.079 1.186 1 -
杉木作为我国特有的用材林树种之一，间伐套种已经成为其营建针阔混交林的主要方式[23]。以往研究表明，间伐可以增加林木的可用资源，促使林木胸径和树高的增加[24-25]，本研究结果与其一致。
在本研究中，CK的高径比、胸高形数比间伐套种处理林分的高径比、胸高形数表现好[26]，但枝下高随着间伐强度的增强呈现先迅速减小后逐渐增大的趋势，与前人结论“枝下高随着林分密度的降低而逐渐增加[27]”不一致。原因与林分密度有关，林分密度大幅度降低后，林内光照和生长空间发生变化，单木下部枝条受到的自身和相邻木的遮阴减少，导致枝下高增长缓慢。但随着林分密度的降低，冠幅逐渐增大，下层冠层受光减少，加剧自然整枝的程度，导致枝下高逐渐增大[28]。
-
林木冠幅生长通常随着林分密度的增大而减小[29]，而在本研究中，随着林分密度的增大，杉木冠幅呈现先增大后减小的趋势。原因可能为Ⅲ处理的间伐强度为65%，伐除木过多，伐除过程中破坏了部分保留木冠层侧枝，导致冠幅因子出现误差。
在本研究中，冠长、冠长率、冠形率随着间伐强度的减弱整体呈现减弱的趋势，与公宁宁[30]对油松人工林树冠、段劼等[31]对侧柏树冠特征因子的研究结果一致。但冠长、冠长率、冠形率在Ⅰ处理出现上升，原因为：随着间伐强度的增强，杉木树高增大，但枝下高先减小后增大[28]并且在处理Ⅰ中最低，导致冠长在Ⅰ处理出现上升，从而间接影响了杉木冠长率、冠形率的生长；冠幅随着林分密度的增大而先增大后减小并且在Ⅱ处理中最高，导致树冠表面积、树冠体积在Ⅱ处理中出现上升。
-
杉木经过间伐套种处理之后，随着间伐强度的增强，杉木的单木材积呈现逐渐增长的趋势，这与张晓红等[32]的研究一致。但随着间伐强度的增强，林分蓄积量逐渐降低。原因在于林分的蓄积量受到单木材积和单位面积株数的共同影响，虽然各间伐套种处理间杉木的单木材积均有所增加，但是间伐之后保留的株数减少。如果间伐强度过大则会因为林分增长的蓄积量不足与伐除木的蓄积量抵消，进而出现单木材积增加但是林分蓄积量减少的情况，这与郑鸣鸣[33]对杉木中龄林间伐的研究、SULLIVAN[34]对美国黑松不同间伐强度的研究结果一致。
林分密度是影响林分材种结构、出材数量和质量的重要因素[35]，间伐套种处理后林分内保留株数减少，林木生长空间变大，促进林木径阶偏移，导致林分材种结构发生变化[36]。本研究发现，间伐套种有利于提高林分用材出材率，降低薪材和废材出材率，有效提高了林分材种效益，与前人的研究结果一致[37]。宋重升[38]以16年生杉木为研究对象，发现随着间伐强度的增强，小径材和小条木出材量和出材率呈现减弱的趋势，与本研究结果一致。此外，宋重升发现小条木仅出现在未间伐林分内，而在本研究中，各间伐套种处理中均出现了小条木，原因可能为初始间伐时林分林龄较大，并且间伐套种处理时间较短，导致林分胸径生长变化较小，仍有较多林木胸径位于14径阶下。王有良等[39]认为大径材蓄积量随保留密度的增大呈先增加而后减小的趋势，这与本研究结果不完全一致。可能是受各径阶活立木株数的影响，随着间伐强度的增大，林木胸径递增，径阶后移，但林木株数逐渐减少，导致林木大径材出材量呈现先增加后减小再增大的趋势[40]。
-
样地经过不同间伐套种模式改造之后，林内环境发生改变，导致林分生长形质性状、空间利用能力、林分产量均发生变化，本研究对其进行综合分析，发现得分由高到低依次为Ⅲ＞Ⅱ＞Ⅰ＞CK，表明间伐套种处理可以有效改善杉木人工林的生长特征。
-
间伐与阔叶树套种处理提高了23 a林龄的杉木人工林的规格材出材率，降低了非规格材出材率，有利于杉木大径材培育；CK是提升林木干形形质的最佳处理措施，提高了杉木的高径比、胸高形数和枝下高；Ⅱ处理（56%间伐 + 套种阔叶树）是提升林木空间利用能力的最佳处理措施，提高了杉木的冠幅、冠长、冠长率和冠形率、树冠表面积、树冠体积；Ⅲ处理（65%间伐 + 套种阔叶树）是快速提升杉木人工林生长特征的最佳处理措施，提高了杉木人工林林分特征的综合得分。
间伐套种对杉木人工林生长、干形形质和材种结构的影响
Effects of Thinning and Interplanting on the Tree Growth, Stem-form Quality and Timber Structure of Cunninghamia lanceolata
-
摘要:
目的 分析间伐套种措施对杉木林木生长形质性状、空间利用能力和林分产量的影响，明确林分经营方向，为杉木人工纯林改造提供科学依据。 方法 对安徽省青阳县23 a林龄的杉木人工林设置3种间伐强度（Ⅰ：47%；Ⅱ：56%；Ⅲ：65%），林下套种2种阔叶树浙江楠（Phoebe chekiangensis）和檫木（Sassafras tzumu），并以未间伐套种的杉木人工纯林为对照（CK：0%间伐，保留密度1 740株·hm−2），探究不同间伐与阔叶树套种处理对杉木人工林的影响，并通过主成分分析法综合评价间伐与阔叶树套种对杉木生长的影响。 结果 （1）间伐与阔叶树套种促进了杉木生长，Ⅲ处理的胸径、树高、单木材积、用材出材率均为最大，分别比CK增加了32.07%、21.60%、83.33%、6.77%（P<0.05）；但高径比和胸高形数、小条木、小径材、薪材和废材的出材率均低于CK，分别降低了8.64%、3.70%、90.88%、41.47%、54.50%、5.71%。（2）间伐与阔叶树套种改变了杉木的空间利用能力，冠幅、冠长、冠长率、冠形率、树冠表面积和树冠体积均显著高于CK（P<0.05），分别提高了12.16%～31.00%、55.58%～61.24%、30.00%～40.00%、16.11%～44.94%、67.17%～95.86%、84.74%～146.63%。（3）间伐套种处理可以有效改善杉木人工林的生长特征，各处理的综合得分由高到低依次为Ⅲ（1.186）>Ⅱ（1.092）>Ⅰ（0.671）>CK（−0.709）。 结论 杉木人工林经营若以快速提升杉木人工林生长特征为目的，本研究中Ⅲ处理（65%间伐 + 套种阔叶树）是最佳处理措施。 Abstract:Objective To analyze the effects of thinning and interplanting measures on the growth traits, spatial use capacity and stand yield for providing a scientific basis for managing Chinese fir plantation. Method The response of Cunninghamia lanceolata stand quality to thinning and interplanting was explored by comparing the pure forest of control treatment (CK: thinning intensity 0%, 1 740 trees left per hectare) with thinned forests (Ⅰ: 47%; Ⅱ: 56%; Ⅲ: 65%) and replanted Phoebe chekiangensis and Sassafras tzumu , and the improvement effect of thinning and interplanting on the growth of Chinese fir was comprehensively evaluated by principal component analysis. Result (1) Thinning and broadleaf trees interplanting promoted tree growth of Chinese fir. Forest had the highest DBH, tree height, individual tree volume, and timber yield under treatment III, which increased by 32.07%, 21.60%, 83.33%, and 6.77% respectively compared to CK (P<0.05) . However, the height-diameter ratio, breast-height form factor, the yield of small strips, small diameter wood, fuelwood and waste wood were lower than CK, with a decrease of 8.64%, 3.70%, 90.88%, 41.47%, 54.50% and 5.71% respectively compared with CK. (2) The space use capacity of Chinese fir plantations changed significantly after thinning and broadleaf trees interplanting. The crown width, crown length, crown ratio, crown shape, canopy surface and canopy volume were all significantly higher than those of CK, with the growth rates of 12.16%～31.00%, 55.58%～61.24%, 30.00%～40.00%, 16.11%～44.94%, 67.17%～95.86%, 84.74%～146.63%. (3) The thinning and interplanting treatments can effectively improve the growth of Chinese fir. The comprehensive score in treatment III (1.186) was the highest, followed by II (1.092), I (0.671), and CK (-0.709). Conclusion If the Chinese fir plantation management is to rapidly improve the growth, treatment III (65% thinning and interplanting broadleaf trees) is the best measure. -
图 1 不同间伐套种下林木胸径和树高生长变化
Figure 1. Changes of DBH and tree height growth under different thinning and interplanting
图 2 不同间伐套种下林木干形生长变化
Figure 2. Changes of stem-form growth under different thinning and interplanting
图 3 间伐套种后林分总出材量的结构
Figure 3. Structure of total stand output after thinning and interplanting
图 4 间伐套种后林分用材出材量的结构
Figure 4. Structure of materials output after thinning and interplanting
表 1 样地基本概况
Table 1. General situation of sample plots
处理
Treatment杉木保留密度
The retained density/
(tree·hm−2)套种株数
The interplanting/
(tree·hm−2)坡向
Aspect坡度
Slope/(°)杉木平均胸径（2018）
Average DBH
(2018)/cm杉木平均树高（2018）
Average tree height
(2018)/mCK（对照组） 1 740 0 东南 15-20 13.1 10.7 Ⅰ（47%间伐 + 套种） 1 050 450 东南 15～20 15.1 11.4 Ⅱ（56%间伐 + 套种） 825 450 东 15～20 15.6 11.8 Ⅲ（65%间伐 + 套种） 675 450 东 15～20 17.2 12.8 
下载: 导出CSV
表 2 不同处理组树冠的生长变化
Table 2. Growth changes of crown of different treatments
处理
Treatment冠幅
Crown width/m冠长
Crown length/m冠长率
Crown ratio冠形率
Crown shapeCK 3.29 ± 0.80 D 5.65 ± 2.17 B 0.50 ± 0.17 C 1.78 ± 0.72 C Ⅰ 3.69 ± 0.85 C 8.93 ± 2.13 A 0.70 ± 0.13 A 2.58 ± 0.94 A Ⅱ 4.31 ± 0.63 B 8.79 ± 2.24 A 0.65 ± 0.13 B 2.07 ± 0.55 B Ⅲ 3.99 ± 0.67 A 9.11 ± 1.81 A 0.66 ± 0.11 B 2.35 ± 0.61 A 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).
下载: 导出CSV
表 3 不同处理组三维树冠的生长变化
Table 3. Growth changes of three-dimensional canopy of different treatments
处理
Treatment树冠表面积均值
Average canopy surface
area/m2树冠表面积变异系数
Coefficients of variation for
canopy surface area/%树冠体积均值
Average canopy
volume/m3树冠体积变异系数
Coefficients of variation for
canopy volume/%CK 31.65C ± 15.73 49.70 17.82D ± 12.46 69.92 Ⅰ 52.91B ± 17.59 33.25 32.92C ± 17.09 51.91 Ⅱ 61.99A ± 18.74 30.23 43.95A ± 18.09 41.16 Ⅲ 58.94AB ± 15.71 26.65 39.15B ± 14.98 38.26 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).
下载: 导出CSV
表 4 间伐套种后杉木单木材积和林分蓄积的生长变化
Table 4. Growth changes of single tree and stand volume of fir plantation after thinning and interplanting
指标
Index处理CK
Treatment CK处理Ⅰ
Treatment Ⅰ处理Ⅱ
Treatment Ⅱ处理Ⅲ
Treatment Ⅲ单木材积 Single tree volume/m3 0.12 ± 0.08 C 0.18 ± 0.10 B 0.21 ± 0.09 AB 0.22 ± 0.10 A 林分蓄积 Stand volume/(m3·hm−2) 226.12 ± 37.11 A 163.94 ± 22.84 BC 162.47 ± 33.27 C 142.58 ± 18.26 D 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).
下载: 导出CSV
表 5 间伐套种处理后不同材种出材率
Table 5. The timber assortments timber-produced rate after thinning and interplanting
处理
Treatment小条木
Small strips
of wood小径材
Small diameter
wood中径材
Medium diameter
wood大径材
Large diameter
wood用材
Materials薪材
Firewood废材
Waste woodCK 10.64 ± 19.05 A 42.87 ± 30.23 A 17.57 ± 29.38 C 0.73 ± 5.51 B 71.80 ± 6.85 C 6.66 ± 6.84 A 21.54 ± 1.93 A Ⅰ 4.23 ± 12.41 B 35.15 ± 32.08 B 29.16 ± 33.48 B 6.36 ± 20.05 A 74.90 ± 4.05 B 4.23 ± 2.82 BC 20.87 ± 1.29 C Ⅱ 1.24 ± 7.77 B 26.77 ± 33.55 BC 42.63 ± 34.42 A 5.53 ± 16.05 A 76.16 ± 4.05 AB 3.44 ± 4.00 CD 20.40 ± 0.73 BD Ⅲ 0.97 ± 6.92 B 25.09 ± 31.26 C 42.58 ± 32.99 A 8.02 ± 20.71 A 76.66 ± 2.68 A 3.03 ± 1.95 D 20.31 ± 1.68 B 注：不同大写字母表示相同时间下不同间伐套种处理之间具有显著性差异（P<0.05）。
Note: Different capital letters indicate significant differences between different thinning and interplanting treatments at the same time (P<0.05).
下载: 导出CSV
表 6 主成分提取结果
Table 6. Extraction results of principal compositions
成分
Composition初始特征值
Initial eigenvalue提取载荷平方和
Extract the sum of squares and load it总计
Total方差百分比
% of variance/%累积百分比
Accumulate/%总计
Total方差百分比
% of variance/%累积百分比
Accumulate/%1 8.831 46.479 46.479 8.831 46.479 46.479 2 3.110 16.370 62.849 3.110 16.370 62.849 3 1.791 9.428 72.276 1.791 9.428 72.276 4 1.342 7.065 79.341 1.342 7.065 79.341 5 1.199 6.311 85.652 1.199 6.311 85.652 6 1.059 5.576 91.228 1.059 5.576 91.228
下载: 导出CSV
表 7 主成分特征值与荷载结果
Table 7. Principal component characteristic values and load results
指标
Index成分 Composition Y1 Y2 Y3 Y4 Y5 Y6 X1 0.960 −0.147 0.125 0.136 0.095 −0.006 X2 0.890 −0.107 −0.201 0.135 0.044 −0.042 X3 0.091 −0.947 0.029 −0.011 0.013 0.009 X4 −0.663 0.204 −0.305 −0.022 −0.004 0.020 X5 −0.919 0.172 0.209 −0.135 0.009 0.039 X6 0.722 −0.008 0.103 −0.594 −0.120 0.189 X7 0.313 0.910 −0.140 0.078 −0.002 −0.044 X8 0.205 0.690 −0.302 0.562 0.101 −0.193 X9 0.826 0.459 −0.024 −0.226 −0.028 0.086 X10 0.801 0.376 0.063 −0.334 −0.042 0.135 X11 −0.709 0.296 0.520 0.106 −0.177 0.202 X12 −0.372 −0.180 −0.862 −0.162 0.046 0.168 X13 0.703 −0.021 0.432 0.029 −0.232 −0.503 X14 0.397 0.012 0.242 0.275 0.526 0.644 X15 −0.698 0.211 0.214 −0.278 0.530 −0.223 X16 −0.618 0.142 0.192 0.242 −0.592 0.371 X17 0.833 −0.239 −0.256 0.196 −0.338 0.109 X18 0.922 −0.069 0.231 0.178 0.189 0.077 X19 −0.362 −0.453 0.103 0.358 0.126 −0.120
下载: 导出CSV
表 8 主成分得分及综合得分
Table 8. Principal component score and comprehensive score
处理
Treatment主成分
Principal component综合得分
Comprehensive score排名
RankingY1 Y2 Y3 Y4 Y5 Y6 CK −1.322 −0.658 0.063 0.143 0.008 −0.067 −0.709 4 Ⅰ 1.050 1.181 −0.309 0.161 0.046 0.084 0.671 3 Ⅱ 2.191 0.650 0.136 −0.632 −0.108 0.105 1.092 2 Ⅲ 2.308 0.745 −0.036 −0.159 0.028 0.079 1.186 1
下载: 导出CSV
-
[1] 国家林业和草原局. 中国森林资源报告(2014—2018)[M]. 北京: 中国林业出版社, 2019. [2] 赵铭臻, 刘 静, 邹显花, 等. 间伐施肥对杉木中龄林生长和材种结构的影响[J]. 南京林业大学学报(自然科学版), 2023, 47(2):70-78. [3] Lin N, Deng N, Lu D. et al. Short-Term Effects of Thinning on Tree Growth and Soil Nutrients in the Middle-Aged Chinese Fir (Cunninghamia lanceolata (Lamb. ) Hook. ) Plantations[J]. Forests, 2023, 14(1): 74. doi: 10.3390/f14010074 [4] 舒韦维, 卢立华, 李 华, 等. 林分密度对杉木人工林林下植被和土壤性质的影响[J]. 生态学报, 2021, 41(11):4521-4530. [5] 方碧江. 杉木灰木莲混交林生长情况及土壤理化性质[J]. 森林与环境学报, 2022, 42(1):46-52. [6] 李方兴, 张意苗, 易伟东, 等. 马尾松、木荷纯林及混交林的生长差异分析[J]. 南方林业科学, 2016, 44(5):17-20. [7] Bi H Q, Long Y S, Turner J, et al. Additive prediction of aboveground biomass for Pinus radiata (D. Don) plantations[J]. Forest Ecology and Management, 2010, 259(12): 2301-2314. doi: 10.1016/j.foreco.2010.03.003 [8] Jacobs D F, Selig M F, Severeid L R. Aboveground carbon biomass of plantation-grown American chestnut (Castanea dentata) in absence of blight[J]. Forest Ecology and Management, 2009, 258(3): 288-294. doi: 10.1016/j.foreco.2009.04.014 [9] Ming A G, Jia H Y, Zhao J L, et al. Above-and below-ground carbon stocks in an indigenous tree (Mytilaria laosensis) plantation chronosequence in subtropical China[J]. PLoS One, 2014, 9(10): e109730. doi: 10.1371/journal.pone.0109730 [10] Wang H, Liu S R, Wang J X, et al. Effects of tree species mixture on soil organic carbon stocks and greenhouse gas fluxes in subtropical plantations in China[J]. Forest Ecology and Management, 2013, 300: 4-13. doi: 10.1016/j.foreco.2012.04.005 [11] Li W H. Degradation and restoration of forest ecosystems in China[J]. Forest Ecology and Management, 2004, 201(1): 33-41. doi: 10.1016/j.foreco.2004.06.010 [12] Zhu J J, Yang K, Yan Q L, et al. Feasibility of implementing thinning in even aged Larix olgensis plantations to develop unevenaged larch-broadleaved mixed forests[J]. Journal of Forest Resource, 2010, 15: 71-80. doi: 10.1007/s10310-009-0152-6 [13] 欧建德. 不同林窗马尾松林对林下南方红豆杉形质生长的影响[J]. 西南林业大学学报(自然科学), 2020, 40(5):10-17. [14] 李婷婷, 陆元昌, 庞丽峰, 等. 杉木人工林近自然经营的初步效果[J]. 林业科学, 2014, 50(5):90-100. [15] 孙冬婧, 温远光, 罗应华, 等. 近自然化改造对杉木人工林物种多样性的影响[J]. 林业科学研究, 2015, 28(2):202-208. [16] 欧建德, 吴志庄. 林下套种对南方红豆杉树冠形态结构及干形变化的影响[J]. 西南林业大学学报, 2016, 36(5):106-110. [17] 陈莹莹. 间伐补植阔叶树大苗对杉木人工林生长的影响[D]. 长沙: 中南林业科技大学, 2017. [18] 郭丽玲, 潘 萍, 欧阳勋志, 等. 间伐补植对马尾松低效林生长及林分碳密度的短期影响[J]. 西南林业大学学报(自然科学), 2019, 39(3):48-54. [19] 彭文成, 黄士绮, 杨 佳, 等. 马占相思人工林近自然改造群落结构转型2年期效果初报[J]. 热带林业, 2023, 51(1):30-35. [20] 刘思泽, 尹海锋, 沈 逸, 等. 间伐强度对马尾松人工林间伐初期林下植被群落物种组成和多样性的影响[J]. 应用生态学报, 2020, 31(9):2866-2874. [21] 翟凯燕, 马婷瑶, 金雪梅, 等. 间伐对马尾松人工林土壤活性有机碳的影响[J]. 生态学杂志, 2017, 36(3):609-615. [22] 葛兆轩, 苑美艳, 黄选瑞, 等. 塞罕坝华北落叶松人工林不同经营模式效果评价[J]. 林业科学研究, 2020, 33(5):38-47. [23] 王俊男. 林下套种楠木对杉木人工林生态系统的影响[D]. 福州: 福建农林大学, 2018. [24] Ford K R, Jan K. Competition alters tree growth response to climate at individual and stand scales[J]. Canadian Journal of Forest Research, 2017, 47: 53-62. doi: 10.1139/cjfr-2016-0188 [25] 姜 丽, 张雄清, 段爱国, 等. 不同林分密度指标在杉木单木直径年生长模型的应用[J]. 林业科学研究, 2022, 35(4):123-129. [26] 韩明跃, 李莲芳, 郑 畹, 等. 间伐强度对云南松中龄低产林分结构的调整研究[J]. 中南林业科技大学学报, 2011, 31(2):27-33. [27] 熊光康, 厉月桥, 熊有强, 等. 低密度造林对杉木生长、形质和材种结构的影响[J]. 南京林业大学学报(自然科学版), 2021, 45(3):165-173. [28] 童书振, 盛炜彤, 张建国. 杉木林分密度效应研究[J]. 林业科学研究, 2002, 15(1):66-75. [29] 苏日娜. 不同经营措施对长白落叶松人工林生长和土壤性质的影响[D]. 哈尔滨: 东北林业大学, 2022. [30] 公宁宁, 马履一, 贾黎明, 等. 不同密度和立地条件对北京山区油松人工林树冠的影响[J]. 东北林业大学学报, 2010, 38(5):9-12. [31] 段 劼, 马履一, 贾黎明, 等. 北京地区油松人工林树冠竞争因子的测算与分析[J]. 东北林业大学学报, 2012, 40(3):14-18. [32] 张晓红, 张会儒, 卢 军, 等. 目标树抚育间伐对蒙古栎天然次生林生长的初期影响[J]. 林业科学, 2020, 56(10):83-92. [33] 郑鸣鸣, 任正标, 王友良, 等. 间伐强度对杉木中龄林生长和结构的影响[J]. 森林与环境学报, 2020, 40(4):369-376. [34] Sullivan T P, Sullivan D S. Acceleration of old-growth structural attributes in lodgepole pine forest: tree growth and stand structure 25 years after thinning[J]. For Ecol Manage, 2016, 365: 96-106. doi: 10.1016/j.foreco.2016.01.029 [35] 卢立华, 农 友, 李 华, 等. 保留密度对杉木人工林生长和生物量及经济效益的影响[J]. 应用生态学报, 2020, 31(3):717-724. [36] Zachara T. The influence of selective thinning on the social structure of the young Scots pine stand[J]. Prace Instytutu Badawczego Lesnictwa Seria A, 2000, 3: 35-61. [37] Nord-Larsen T, Cao Q V. A diameter distribution model for evenaged beech in Denmark[J]. Forest Ecology and Management, 2006, 231: 218-225. doi: 10.1016/j.foreco.2006.05.054 [38] 宋重升, 王有良, 张利荣, 等. 间伐强度对杉木人工林材种结构的影响[J]. 福建农林大学学报(自然科学版), 2022, 51(2):195-203. [39] 王有良, 宋重升, 何宗明, 等. 杉木材种结构及大径材生长保留密度与立地效应[J]. 北华大学学报(自然科学版), 2022, 23(2):246-252. [40] 程亚平, 吴伟华, 姚 俊, 等. 不同立地及间伐强度对杉木人工林生长的影响[J]. 南方林业科学, 2016, 44(2):24-2751. -
点击查看大图
图(4) / 表(8)
计量
- 文章访问数: 3669
- HTML全文浏览量: 1121
- PDF下载量: 94
- 被引次数: 0
