-
马尾松Pinus massoniana曾经是南方集体林区的主要造林树种,在荒山绿化、生态服务和发展经济中起过不可替代的重要作用。但长期大面积纯林营造,林地出现地力衰退、生产力降低、生态服务功能下降等不良现象,特别是马尾松林经营长期遭受松材线虫Bursaphelenchus xylophlius病的困扰。据国家林业和草原局统计,近30 a来,全国因松材线虫病损失的松树累计达数十亿株,造成的直接经济损失和生态服务价值损失高达上千亿元,松材线虫病已成为全球森林生态系统中最具危险性、毁灭性的病害之一[1]。对于松材线虫病发生的疫区疫点,相关办法都要求“加强疫情科学除治,全面实施以清理病死、枯死、濒死松树为核心措施”。各地据此采取了皆伐、择伐等干扰方式清除疫木,以防止松材线虫的进一步扩散。大多相关研究集中在松材线虫病入侵后皆伐形成的除治迹地上马尾松林群落演替、林分结构、植物多样性和种群生态位的评价上[1-4]。但采伐疫木对受害马尾松林综合林分结构长期影响的相关研究较少。林分结构不但包括树种、胸径、树高、年龄、林木个体大小差异等非空间结构,也包括林木分布格局、单木竞争指数等空间结构,是林分未来生长变化的决定性因素之一[5-6]。在生产实践中,采伐疫木后的马尾松群落结构将会发生哪些变化,有必要对受害马尾松林的健康经营开展研究。本研究在浙江省杭州市临安区选择松材线虫病入侵过的马尾松人工纯林,以未受松材线虫病危害的林分为对照,探讨疫木卫生伐后的天然更新对马尾松人工纯林林分结构的影响,比较2种类型的林分在树种组成、胸径结构、树高结构、林木个体大小差异、聚集指数和竞争指数等方面的差异,以期为松材线虫病入侵后的马尾松林分的经营管理提供参考。
-
在杭州市临安区河桥镇选择1994年同一批造林的马尾松纯林作为研究对象。21世纪初该马尾松人工纯林有部分遭受松材线虫病,为防止病害蔓延,2005年开始每年进行卫生伐,其他马尾松林木均予以保留,保护天然更新的非松科Pinaceae植株,持续7 a,共采伐约1 500株·hm−2。疫木伐除后,任其自然恢复更新。2017年8月,按照立地环境基本因子一致的原则,从中分别选择未入侵过的小块马尾松林(对照)和卫生伐过的马尾松林作为调查对象,2种类型的林分分别设置了5块30 m×30 m的标准样地。样地土壤为黄壤,土层厚度为0.8 m,坡向为西北坡,坡度35°~40°,立地条件差,pH 6.2。样地边界和林木坐标用全站仪测定,对样地进行了每木检尺,起测胸径为5 cm,调查因子包括树种、胸径、坐标、树高等。样地基本情况见表1。
表 1 卫生伐后的马尾松林分基本特征
Table 1. Characteristics of P. massoniana stands after sanitation thinning diseased trees
处理 密度/(株·hm−2) 胸径/cm 树高/m 郁闭度 马尾松 阔叶树 马尾松 阔叶树 马尾松 阔叶树 对照 1 800 13.0±2.4 12.5±1.7 0.9 卫生伐 300 1 000 17.6±2.8 7.1±1.8 12.5±2.2 5.6±1.1 0.6 -
物种丰富度是指样地内出现的物种数目。树种重要值是研究某个种在群落中的地位和作用的综合数量指标,物种丰富度和树种重要值反映了植物多样性。重要值计算采用的公式为:重要值
$= $ (相对多度+相对频度+相对显著度)/3。其中:相对多度$= $ (某个种的多度/所有物种的多度之和)×100%;相对频度$= $ (某个种的频度/所有物种的频度之和)×100%;相对显著度$= $ (某个种的胸高断面积之和/所有物种的胸高断面积之和)×100%。 -
林木个体大小差异在植物种群中普遍存在,相邻林木个体的大小会影响目标个体的胸径、高度的生长,在群落发展和维持群落物种多样性方面也具有深远意义,反映植物个体大小差异的指数采用库兹涅茨系数(KC)、基尼系数(GC)和洛伦茨不对称系数(LAC)[8]。
用二元材积公式[9]计算马尾松和阔叶树的单株林木材积,计算公式为[4]:
$$ {V_{{\text{马尾松}}}} = 0.000\;062\;341\;803{D^{1.855\;149\;7}}{H^{0.956\;824\;92}}\text{;} $$ (1) $$ {V_{{\text{阔叶树}}}} = 0.000\;068\;563\;4{D^{1.933\;221}}{H^{0.867\;885}}\text{。} $$ (2) 式(1)~(2)中:V为单株材积(m3);D为胸径(cm);H为树高(m)。
库兹涅茨系数(KC)是指单株材积为林分前20%的林木个体的材积之和在林分蓄积量中所占的比例,这个指数的数值越大,表示林木分化越严重,个体差异越大。基尼系数( GC)可用于研究林分内各林木之间差异性的大小[4, 8],计算公式为:
$$ {G_{\rm{C}}} = \frac{{\mathop \sum \limits_{i = 1}^n \left( {2i - n - 1} \right){x_i}}}{{{n^2}u}}\text{。} $$ (3) 式(3)中:n为林分内林木的株数;xi为林木个体按材积大小升序排列时对应的第i株林木的材积;u为林分平均材积。洛伦茨不对称系数(LAC)用于反映林木个体大小等级对不一致性程度的贡献,当LAC>1时,表示不一致性的产生主要在少数较大的林木个体中;当LAC<1时,表示不一致性的产生主要在大多数较小的林木个体中[4, 8]。将林分n株林木按单株材积从小到大排序(x1,x2,…,xn),LAC的计算公式为:
$$ {L_{\rm{AC}}} = \frac{{m + \theta }}{n} + \frac{{{L_m} + \theta {X_{m + 1}}}}{{{L_n}}}\text{;} $$ (4) $$ \theta = \frac{{u - {X_m}}}{{{X_{m + 1}} - {X_m}}}\text{。} $$ (5) 式(4)~(5)中:m为恰好小于林木平均材积μ的林木序号;Lm为前m株林木个体材积之和;Ln为n株林木个体材积之和;θ为参数;Xm为排序第m株林木的材积;Xm+1为排序第m+1株林木的材积。
-
聚集指数是由CLARK等[12]提出的,用于检验种群空间分布格局的常用指数,是最早采用的与距离有关的空间格局指数[13]。采用最近邻单株距离的平均值与随机分布下的期望平均距离之比来表示:
$$ R = \dfrac{{\dfrac{1}{N}\mathop \sum \limits_{i = 1}^N {r_i}}}{{\dfrac{1}{2}\sqrt {\dfrac{S}{N}} }}\text{。} $$ (6) 式(6)中:R为聚集指数;ri为第i株树木到其最近邻木的距离;N为样地内林木株数;S为样地面积。若R>1,林木呈均匀分布;若R<1,则林木趋于聚集分布;若R
$= $ 1,林木则表现为随机分布。 -
竞争指数采用HEGYI[14]的计算公式:
$$ {C_i} = \mathop \sum \limits_{j = 1}^n \frac{{{d_j}}}{{{d_i} {L_{ij}}}}\text{。} $$ (7) 式(7)中:Ci为对象木i的竞争指数;Lij为对象木i与竞争木j之间的距离;di为对象木i的胸径;dj为竞争木j的胸径;n为对象木i的竞争木株数,本研究采用4邻近树法,n
$= $ 4。 -
采用Excel 2003进行数据处理与图表制作,SPSS 20.0进行单因素方差分析(one-way ANOVA),LSD方法进行显著性检验(α
$= $ 0.05)。利用ArcGIS和ForStat软件对2种林分进行空间分析。数据为平均值±标准差。 -
由表2可知:马尾松人工纯林的重要值仍然为100.00%。在受害马尾松纯林疫木卫生伐后的林分中,马尾松的重要值下降至38.41%,仍然为群落的第1优势种;出现了大量天然更新的阔叶树种,短柄枹Quercus glandulifera的重要值为35.59%,是群落的第2优势种;檵木Loropetalum chinense的重要值为10.91%,是群落下层林分中主要的伴生树种;白栎Quercus fabri的重要值为10.54%,是继短柄枹之后的又一主要乔木树种;另外还有少量的无患子Sapindus mukorossi、枫香Liquidambar formosana、薄叶山矾Symplocos anomala和合欢Albizia julibrissin等乔木树种分布。可见,松材线虫入侵后的马尾松人工纯林在疫木卫生伐后经过长期恢复已演替为马尾松落叶阔叶混交林。而且该混交林分没有发现马尾松的小树和幼树,预示着未来该林分的阔叶树重要值将超过马尾松,形成阔叶树马尾松混交林。
表 2 疫木卫生伐后马尾松林分的树种组成
Table 2. Tree species composition of P. massoniana stands after sanitation thinning infected trees
处理 树种 重要值/% 对照 马尾松 100.00 卫生伐 马尾松 38.41 短柄枹 35.59 檵木 10.91 白栎 10.54 无患子 2.24 枫香 1.25 薄叶山矾 0.55 合欢 0.52 -
由图1可知:马尾松人工纯林的林分平均胸径为13.0 cm,其中,10、12和14 cm径阶的林木株数最多,分别为22.8%、29.6%和21.8%,呈单峰状左偏正态分布,小径阶和大径阶林木各占一小部分。疫木卫生伐后形成的马尾松阔叶混交林林分平均胸径为9.1 cm,林木大多数分布在6和8 cm径阶,株数分别占49.2%和23.4%,而且随着胸径的增大呈下降趋势,呈现反J型曲线,较大径阶22和24 cm的株树比例比马尾松人工纯林高。未遭受松材线虫入侵的马尾松林依然保持单层林;疫木卫生伐后经过长期恢复由单层纯林演变为由马尾松占据的主林层与阔叶树组成的次林层的复层混交林分(表1)。
-
由表3可知:马尾松人工纯林疫木卫生伐后,库兹涅茨系数(KC)增加,说明较大个体林木和较小个体林木分化更严重;基尼系数(GC)增加,说明林木个体大小不一致性增大;洛伦茨不对称系数(LAC)由略小于1增加为略大于1,说明林木个体差异的产生主要在于大多数较小的林木个体变为少数较大的林木个体。
表 3 卫生伐后马尾松林的林木个体大小差异
Table 3. Individual size inequality of P. massoniana stands after thinning diseased trees
处理 KC GC LAC 对照 0.35±0.08 a 0.26±0.05 a 0.97±0.15 a 卫生伐 0.74±0.13 b 0.66±0.09 b 1.04±0.23 a 说明:同列不同字母表示同一指标不同处理间差异显著 (P<0.05) -
由表4可知:马尾松人工纯林聚集指数(R)>1,表现为均匀分布,疫木卫生伐后形成的马尾松阔叶林的R
$= $ 0.95<1,呈现聚集分布。马尾松人工纯林的平均竞争指数(Ci)为8.16,疫木卫生伐后形成的马尾松阔叶林的平均竞争指数增加至10.12。表 4 卫生伐后马尾松林的聚集指数和竞争指数
Table 4. Aggregation index and competition index of P. massoniana stands after sanitation thinning infected trees
处理 聚集指数(R) 竞争指数(Ci) 对照 1.27±0.21 a 8.16±3.52 a 卫生伐 0.95±0.14 b 10.12±5.60 b 说明:同列不同字母表示同一指标不同处理间差异显著 (P<0.05) 由图2可知:马尾松人工纯林各径阶处受到的竞争压力基本一致,未发生很大的变化。疫木卫生伐后形成的马尾松阔叶林于6 cm径阶处林木受到的竞争压力最大,其次是8 cm径阶林木,各径阶处受到的竞争压力差异较大。
Effects of sanitation cutting pine wilt diseased trees on the stand structure of pure Pinus massoniana plantation
-
摘要:
目的 探讨松材线虫Bursaphelenchus xylophlius病疫木采伐对马尾松人工纯林林分结构的影响,对松材线虫病入侵后的马尾松林分的经营管理有理论和实践价值。 方法 在浙江省杭州市临安区,分别选择松材线虫病危害及周边未受害的马尾松Pinus massoniana人工纯林设置样地并进行调查。 结果 疫木卫生伐12 a后的结果表明:与马尾松纯林相比较,卫生伐之后马尾松重要值显著下降,从100.00%降至38.41%;胸径结构由正态分布变为反J型曲线;树高结构由单层林变为复层林;林木大小个体差异显著增加,两级分化显著加剧;林木分布从均匀分布变为随机分布;林木竞争压力增加,6 cm径阶林木竞争压力最大。 结论 卫生伐促使受害马尾松纯林向混交林演替的速度加快,多个林分结构指标发生了重大变化,单层马尾松同龄纯林演替为复层马尾松落叶阔叶树混交异龄林。图2表4参18 Abstract:Objective This study aims to explore the effects of sanitation cutting of the diseased pine trees infected by Bursaphelenchus xylophlius on the stand structure of pure Pinus massoniana plantation and provide a reference for the management of P. massoniana invaded by B. xylophlius. Method Field plots were selected of both a pure P. massoniana plantation infected by B. xylophlius and its neighbouring plantation uninfected to be investigated in Lin’an, Zhejiang, China. Result After 12 years of sanitation cutting, compared with the uninfected pure P. massoniana plantation, it is shown that the importance values of P. massoniana significantly decreased from 100.00% to 38.41%. The stand diameter distribution followed an inverse J-shaped curve. In terms of height stricture, single storied stand changed into multi-storied stand. The individual size inequality significantly increased. The distribution pattern of trees changed from uniform distribution to random distribution. The competitive index of trees significantly increased, and that of trees at 6 cm diameter class was the highest. Conclusion The sanitation cutting accelerated the succession speed of the infected pure P. massoniana, with significant changes in several stand structure indicators and the single-storied even-aged pure P. massoniana stand succeeded by the mixed multi-storied uneven-aged P. massoniana and broad-leaved trees stand. [Ch, 2 fig. 4 tab. 18 ref.] -
表 1 卫生伐后的马尾松林分基本特征
Table 1. Characteristics of P. massoniana stands after sanitation thinning diseased trees
处理 密度/(株·hm−2) 胸径/cm 树高/m 郁闭度 马尾松 阔叶树 马尾松 阔叶树 马尾松 阔叶树 对照 1 800 13.0±2.4 12.5±1.7 0.9 卫生伐 300 1 000 17.6±2.8 7.1±1.8 12.5±2.2 5.6±1.1 0.6 表 2 疫木卫生伐后马尾松林分的树种组成
Table 2. Tree species composition of P. massoniana stands after sanitation thinning infected trees
处理 树种 重要值/% 对照 马尾松 100.00 卫生伐 马尾松 38.41 短柄枹 35.59 檵木 10.91 白栎 10.54 无患子 2.24 枫香 1.25 薄叶山矾 0.55 合欢 0.52 表 3 卫生伐后马尾松林的林木个体大小差异
Table 3. Individual size inequality of P. massoniana stands after thinning diseased trees
处理 KC GC LAC 对照 0.35±0.08 a 0.26±0.05 a 0.97±0.15 a 卫生伐 0.74±0.13 b 0.66±0.09 b 1.04±0.23 a 说明:同列不同字母表示同一指标不同处理间差异显著 (P<0.05) 表 4 卫生伐后马尾松林的聚集指数和竞争指数
Table 4. Aggregation index and competition index of P. massoniana stands after sanitation thinning infected trees
处理 聚集指数(R) 竞争指数(Ci) 对照 1.27±0.21 a 8.16±3.52 a 卫生伐 0.95±0.14 b 10.12±5.60 b 说明:同列不同字母表示同一指标不同处理间差异显著 (P<0.05) -
[1] 吴敏娟, 尤誉杰, 张晓红, 等. 不同干扰模式对受害马尾松人工纯林林分结构的影响[J]. 应用生态学报, 2019, 30(1): 58 − 66. WU Minjuan, YOU Yujie, ZHANG Xiaohong, et al. Effects of different interference modes on the stand structure of artificial pure forest of Pinus massoniana [J]. Chin J Appl Ecol, 2019, 30(1): 58 − 66. [2] 吕文艳, 王柏泉, 曾德山, 等. 松材线虫病疫木不同伐除方式对森林植物群落演替的影响[J]. 湖北林业科技, 2008, 37(3): 4 − 8. LÜ Wenyan, WANG Baiquan, ZENG Deshan, et al. Effects of different cutting methods of pine wood nematode disease on forest plant community succession [J]. Hubei For Sci Technol, 2008, 37(3): 4 − 8. [3] 王国明, 赵颖, 陈斌, 等. 松材线虫病除治迹地自然恢复过程中物种多样性的动态变化[J]. 浙江林学院学报, 2010, 27(2): 170 − 177. WANG Guoming, ZHAO Ying, CHEN Bin, et al. Dynamic changes of species diversity during the natural restoration process of pine wood nematode disease [J]. J Zhejiang For Coll, 2010, 27(2): 170 − 177. [4] 石娟, 骆有庆, 宋冀莹, 等. 松材线虫入侵后不同伐倒干扰强度对马尾松林植物多样性的影响[J]. 应用生态学报, 2006, 17(7): 1157 − 1163. SHI Juan, LUO Youqing, SONG Jiying, et al. Effects of different cutting intensity on the plant diversity of Pinus massoniana forest after invasion of pine wood nematode [J]. Chin J Appl Ecol, 2006, 17(7): 1157 − 1163. [5] 孟宪宇. 测树学[M]. 北京: 中国林业出版社, 2006. [6] 吴建强, 王懿祥, 杨一, 等. 干扰树间伐对杉木人工林林分生长和林分结构的影响[J]. 应用生态学报, 2015, 26(2): 340 − 348. WU Jianqiang, WANG Yixiang, YANG Yi, et al. Effects of interfering tree thinning on the growth and stand structure of Chinese fir plantation [J]. Chin J Appl Ecol, 2015, 26(2): 340 − 348. [7] 朱婷婷, 王懿祥, 朱旭丹, 等. 遮光对木荷和枫香光合特性的影响[J]. 浙江农林大学学报, 2017, 34(1): 28 − 35. ZHU Tingting, WANG Yixiang, ZHU Xudan, et al. Effect of shading on photosynthetic characteristics of Schima superba and Liquidambar formosana [J]. J Zhejiang A&F Univ, 2017, 34(1): 28 − 35. [8] 王懿祥, 张守攻, 陆元昌, 等. 林木个体大小不一致性指标对人工林间伐方式的即时性响应[J]. 应用生态学报, 2014, 25(6): 1645 − 1651. WANG Yixiang, ZHANG Shougong, LU Yuanchang, et al. The immediate response of individual size inconsistency indexes to the mode of artificial forest thinning [J]. Chin J Appl Ecol, 2014, 25(6): 1645 − 1651. [9] 王懿祥, 张守攻, 陆元昌, 等. 干扰树间伐对马尾松人工林目标树生长的初期效应[J]. 林业科学, 2014, 50(10): 67 − 73. WANG Yixiang, ZHANG Shougong, LU Yuanchang, et al. The initial effect of interfering tree thinning on the growth of target trees in Pinus massoniana plantation [J]. Sci Silv Sin, 2014, 50(10): 67 − 73. [10] 周红敏, 惠刚盈, 赵中华, 等. 林分空间结构分析中样地边界木的处理方法[J]. 林业科学, 2009, 45(2): 1 − 5. ZHOU Hongmin, HUI Gangying, ZHAO Zhonghua, et al. The treatment method of boundary wood in the spatial structure analysis of stands [J]. Sci Silv Sin, 2009, 45(2): 1 − 5. [11] 仇建习, 汤孟平, 沈利芬, 等. 近自然毛竹林空间结构动态变化[J]. 生态学报, 2014, 34(6): 1444 − 1450. QIU Jianxi, TANG Mengping, SHEN Lifen, et al. Dynamic changes of spatial structure of near-natural bamboo forests [J]. Acta Ecol Sin, 2014, 34(6): 1444 − 1450. [12] CLARK P J, EVANS F C. Distance to nearest neighbor as a measure of spatial relationships in populations [J]. Ecology, 1954, 35(4): 445 − 453. [13] 安慧君, 张韬. 聚集指数边界效应的校正方法与应用[J]. 南京林业大学学报(自然科学版), 2005, 29(3): 57 − 60. AN Huijun, ZHANG Tao. Correction method and application of boundary effect of clustering index [J]. J Nanjing For Univ Nat Sci Ed, 2005, 29(3): 57 − 60. [14] HEGYI F. A simulation model for managing jack-pine stands [J]. Growth Models Tree Stand Simu, 1974, 30: 74 − 90. [15] 王懿祥. 人工马尾松和杉木林目标树经营理论与实践[D]. 北京: 中国林业科学研究院, 2012. WANG Yixiang. Theory and Practice of Target Tree Management in Pinus massoniana and Cunninghamia lanceolata Plantation[D]. Beijing: Chinese Academy of Forestry, 2012. [16] 姜俊, 谢阳生, 陆元昌, 等. 不同林龄阶段马尾松人工林群落结构特征及经营策略[J]. 西北林学院学报, 2015, 30(6): 1 − 7. JIANG Jun, XIE Yangsheng, LU Yuanchang, et al. Community structure and management strategy of Pinus massoniana plantation in different ages [J]. J Northwest For Univ, 2015, 30(6): 1 − 7. [17] 柴希民, 蒋平. 松材线虫的发生和防治[M]. 北京: 中国农业出版社, 2003. [18] 张慧, 周国模, 白尚斌, 等. 目标树抚育对亚热带天然次生灌丛群落结构和多样性的影响[J]. 应用生态学报, 2017, 28(5): 1414 − 1420. ZHANG Hui, ZHOU Guomo, BAI Shangbin, et al. Effects of target tree rearing on the structure and diversity of subtropical natural secondary shrub communities [J]. Chin J Appl Ecol, 2017, 28(5): 1414 − 1420. -
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20190487