-
竞争是生物间相互作用的重要方面之一[1],可分为种内和种间竞争,对林木个体的生长发育、森林群落的结构和群落演替进程都有影响。林木竞争状况多通过林木竞争指数来反映,竞争指数对林木之间的竞争进行了定量描述,在形式上反映了林木个体生长与生存空间的关系,实质上反映了树木对环境资源需求及其在获取环境资源时要承受的竞争压力[2]。竞争随着林龄的增加以及经营措施和环境的改变而发生改变[3],目前国内学者对种内和种间竞争的研究主要体现在优势树种种内和种间竞争现状的分析[1, 4-6],有少量关于群落种内和种间竞争随时间变化的研究[7],但对林分内的竞争关系尤其是种内和种间竞争对经营措施的响应研究较少。研究林分的竞争关系及其对经营措施的响应对维护森林生态系统的稳定,更好实现森林的各项功能,科学的经营和管理森林至关重要[8]。
国内目标树经营相关的研究多集中在对林分生长、结构、更新、生物多样性等的影响[9-12],对林木竞争的影响研究较少。蒙古栎(Quercus mongolica Fisch. ex Ledeb)次生林在我国东北地区分布范围很广,多为阔叶红松(Pinus koraiensis Sieb. et Zucc.)林过伐后形成的处于次生演替阶段的森林群落,林地生产力较低且生态功能退化,如何通过经营来改善其林分状况,提高其经济价值和生态效益有很重要的意义[13]。
本研究以蒙古次生林为对象,分析了林木竞争关系以及不同抚育间伐方式对林分内竞争状态的影响差异。为了解蒙古栎次生林内部竞争关系及抚育间伐方式对竞争关系的影响效应和机制提供依据。
-
2013年在塔子沟林场选择具有代表性的蒙古栎次生林林分(林龄40~70 a),设置12块面积为1 hm2(100 m×100 m)的固定监测样地。调查样地内胸径5 cm及以上的林木,测量并记录坐标、树种、胸径、树高等信息。对12块样地采取完全随机区组设计,各组林木胸径、蓄积及断面积无显著差异。试验设置4种处理,3次重复。4种处理包括:(1)传统经营,即T1处理(1、7、9号样地),遵循《森林抚育规程》(GB/T15781—2009)的规定进行抚育间伐设计,采伐强度控制在25%以下,采伐对象主要是生长不良的林木,伐后郁闭度在0.6以上;(2)目标树经营密度1,即T2处理(2、4、11号样地),按照目标树单株木集约经营理念制定采伐设计,用材目标树选取密度约为70株·hm−2,目标树的树种选择上优先考虑乡土的、价值比较高的以及此类群落演替的顶级树种。另外考虑了林木生命力、质量和距离3种因素,将影响目标树生长的林木确定为采伐木;(3)目标树经营密度2,即T3处理(3、8、10号样地),用材目标树选取密度为100株·hm−2,其余措施与T2一致。(4)对照样地,即T0处理(5,6,12号样地)。为便于对比试验,T0和T1按照100株·hm−2选取了目标树。2018年9月对所有样地进行了复测,抚育经营前样地信息见表1。
样地号
Plot code处理
Treatments树种组成
Species
composition郁闭度
Canopy
density密度/
(株·hm−2)
Stand density胸径 Diameter/cm 树高 Height/m 最小
Min.最大
Max.平均
Average最小
Min.最大
Max.平均
AverageZH-01 T1 6蒙1杨1白1红1色 0.66 825 5.0 66.0 14.6 2.0 29.6 11.0 ZH-02 T2 4蒙2杨1椴1白1色1红 0.71 957 5.0 59.1 13.7 0.5 24.1 11.1 ZH-03 T3 4蒙2白1色1椴1杨1红 0.72 1 115 5.0 55.7 12.9 1.3 28.5 10.2 ZH-04 T2 4蒙2白1椴1红1色1杨 0.83 1 058 5.0 50.2 13.6 1.5 25.8 10.9 ZH-05 T0 5蒙1红1椴1色1白1杨 0.74 1 089 5.0 53.8 15.1 1.4 21.8 10.0 ZH-06 T0 4蒙3白2杂1落 0.46 617 5.0 50.7 15.0 1.8 27.1 13.3 ZH-07 T1 5蒙2杂1白1红1黑 0.85 1 071 5.0 50.9 13.1 1.5 21.8 8.9 ZH-08 T3 5蒙2白1红1杂1黑 0.88 1 162 5.0 58.5 12.8 3.2 31.1 10.9 ZH-09 T1 6蒙2杂1杨1水 0.90 857 5.0 52.2 14.3 1.4 18.9 9.8 ZH-10 T3 6蒙2白1杂1黑 0.77 1 011 5.0 70.1 14.3 3.1 29.4 11.7 ZH-11 T2 4蒙2白1杂1色1红1椴 0.87 1 054 5.0 58.2 14.1 2.2 26.5 10.9 ZH-12 T0 5蒙2落1杂1白1水 0.67 754 5.0 65.3 15.9 1.4 23.4 11.0 注:ZH代表蒙古栎次生林,蒙代表蒙古栎,白代表白桦,杨代表大青杨,红代表红松,色代表色木槭,椴代表糠椴,落代表长白落叶松,黑代表黑桦,胡代表胡桃楸,水代表水曲柳,杂代表其他树种。Notes: ZH stands for natural secondary Quercus mongolica forest, 蒙. Stands for Quercus mongolica, 白. Stands forBetula platyphylla, 杨. Stands for Populus ussuriensis,红. Stands for Pinus koraiensis, 色. Stands for Acer pictum, 椴. Stands for Tilia mandshurica,落. Stands for Larix olgensis, 黑. Stands for Betula dahurica.,水. Stands for Fraxinus mandshurica, 杂. Stands for other tree species. Table 1. General description of sample plots
样地边缘林木的竞争木可能处于边界之外,因此本研究将样地边界5 m内区域视为缓冲区,其中林木只作竞争木来参与计算,核心区(90 m×90 m)的林木既作为竞争木又作为对象木参与计算。
-
采用优势度分析法来确定优势树种[4, 14]。具体方法是计算蒙古栎次生林乔木树种的相对胸高断面积,作为优势度,然后通过公式(1)确定优势树种个数。
式中:d表示确定优势树种个数的阈值,N为总树种数,T为优势度排位在前的树种,即上位种,U为优势度排位在后的树种,及剩余种xi为上位种(T)的优势度,x为优势树种所占的理想百分比,xj为剩余种的优势度。如果乔木层只有一个优势树种,那么理想百分比为100%;如果有两个优势树种,则理想百分比为50%;如果有3个优势树种,理想百分比为33.3%,以此类推来分别计算d值。认为当d取值最小时的上位种数为优势树种个数M1,最后将优势度从大到小排序,选择前M1个树种作为该样地的优势树种。
主要树种确定标准是:12块蒙古栎次生林样地中,树种均出现且断面积占比达到1%。这样能更充分反映林分内的竞争关系,而种内和种间竞争关系仅着重分析了优势树种。
-
选择应用最广泛的Hegyi竞争指数[15],计算方法如下:
其中:Hegyii表示第i株对象木的竞争指数,Dj表示第j个竞争木的胸径,Di表示第i个对象木的胸径,dij表示第i个对象木与第j个竞争木之间的水平距离,n表示对象木的竞争木个数,Hegyi表示总竞争指数,Hegyi_m表示平均竞争指数。种内竞争指数(Hegyi_intra)和种间竞争指数(Hegyi_inter),定义为将竞争木中与对象木树种相同的和不同的分别作为竞争木计算得出的竞争指数。平均种内竞争指数(Hegyi_intra_m)和平均种间竞争指数(Hegyi_inter_m)分别是Hegyi_intra和Hegyi_inter的平均值。这些竞争指数值越大代表竞争越激烈。
-
目前竞争木选择的方法主要有胸高断面积法[16]、固定半径法[15]、影响域法等。一些学者[17]指出竞争木个数选择应为4,而另一些研究有不同见解[18-21]。本研究采用有序样本聚类分析的方法[22]来确定竞争木个数,即依次选取距离对象木最近的n株相邻木作为竞争木,分别计算样地的Hegyi、Hegyi_intra、Hegyi_inter、Hegyi_m、Hegyi_intra_m和Hegyi_inter_m。以这些竞争指数作为分类样本,采用有序样本聚类分析法分成2类,根据分类结果的一致性来确定竞争木个数。有序样本聚类算法采用最优分割法(或Fisher法)进行分类,其递推公式为:
式中
$\varphi $ 表示分类的损失函数,p0表示某一种分类方案,n表示竞争木个数,这里n=3,4,5,$\cdots $ ,14,D表示某类样本离差平方和。 -
根据抚育经营前样地调查的数据,利用公式(1)计算得出优势树种为蒙古栎、红松、白桦和杂木,这些树种胸高断面积和株数分别占样地总量的75.41%和70.26%(表2)。杂木主要是青楷槭(Acer tegmentosum Maxim.)等小乔木树种,多为丛状分布,因此确定样地的优势树种为蒙古栎、红松和白桦;主要树种为蒙古栎、白桦、红松、色木、糠椴、大青杨、水曲柳、落叶松和黑桦(表2)。
优势树种
Dominant tree species胸高断面积
Basal area/(m2·hm−2)比例
Percent/%株数
Tree number/(tree·hm−2)比例
Percent/%d值
d value蒙古栎Quercus mongolica 9.82 47.82 284 29.50 0.022 2 白桦Betula platyphylla 2.82 13.76 81 8.44 0.010 9 红松Pinus koraiensis 1.44 7.00 167 17.30 0.010 3 杂木Other tree species 1.40 6.83 145 15.02 0.010 1 大青杨Populus ussuriensis 1.21 5.92 64 6.59 0.010 2 色木槭Acer pictum 1.06 5.15 92 9.53 0.010 4 糠椴Tilia mandshurica 1.02 4.99 54 5.56 0.010 7 长白落叶松Larix olgensis 0.68 3.29 33 3.39 0.011 1 黑桦Betula dahurica 0.52 2.52 11 1.11 0.011 6 水曲柳Fraxinus mandshurica 0.27 1.29 8 0.85 0.012 2 Table 2. The number and basal area of tree species in Mongolian oak plots
-
利用抚育经营前样地调查数据,样地内主要树种的Hegyi和Hegyi_m随竞争木株数n增加的计算结果如图1所示。
Figure 1. The total (a) and average (b) competition index of the main species based on the different number of competing neighbors
随着竞争木增多,Hegyi和Hegyi_m呈上升趋势,且上升趋势逐步减缓。竞争强度较大树种的竞争排位并不随着选择竞争木数量的增多而改变。而竞争强度较低的树种的竞争排序因竞争木选择的数量不同而发生变化。各竞争指数Hegyi、Hegyi_intra、Hegyi_inter、Hegyi_m、Hegyi_intra_m和Hegyi_inter_m的分类结果一致:[3,4,5,6,7]和[8,9,10,11,12,13,14],这种划分下,组内残差平方和最小,也即7和8为最优分割点。考虑到竞争木株数越多,计算值越稳定,本研究选取8个竞争木来计算Hegyi,从而分析抚育经营前林分的竞争状态及竞争对间伐的响应。
-
单因子方差分析显示Hegyi,Hegyi_intra和Hegyi_inter在T0,T1,T2和T3下无显著差异。样地内主要树种的Hegyi排序为:蒙古栎>红松>色木槭>糠椴>白桦>大青杨>落叶松>水曲柳>黑桦,树种的竞争排序和株数有明显的正相关关系(图2a)。样地主要树种的Hegyi_m排序为:水曲柳>色木槭>糠椴>红松>白桦>蒙古栎>落叶松>大青杨>黑桦(图2b)。种内竞争大于种间竞争的树种为:蒙古栎、白桦、水曲柳,种间竞争大于种内竞争的为:红松、色木槭、糠椴、大青杨、落叶松和黑桦,样地内种内竞争(30132)略低于种间竞争(31937)(图2a)。
-
按照国际林联(IUFRO)的标准划分林层[23]。根据两期调查数据(2013年和2018年)分别对各个林层的Hegyi和Hegyi_m进行统计(表3)。
处理
Treatment林层
Forest storyHegyi Hegyi_m 2013 2018 2013 2018 T0 下林层
Downstory4 175.94 1 957.86 8.94 8.16 中林层
Middle story5 106.43 5 729.61 5.11 5.85 上林层
Up story1 211.92 1 384.12 2.33 2.15 T1 下林层
Downstory4 181.23 595.78 8.67 10.45 中林层
Middle story5 612.26 7 866.84 5.42 6.24 上林层
Up story2 013.44 2 203.51 2.94 2.62 T2 下林层
Downstory4 285.51 2 765.44 8.23 10.72 中林层
Middle story8 441.90 9 218.68 6.32 6.03 上林层
Up story1 801.41 1 453.73 2.58 2.13 T3 下林层
Downstory6 181.70 3 299.76 7.85 8.50 中林层
Middle story6 974.19 8 485.58 5.82 5.48 上林层
Up story1 754.21 1 090.48 2.48 1.95 Table 3. Hegyi and Hegyi_m of different forest story in different treatment
2013年各林层的Hegyi排序
为中林层>下林层>上林层,2018年T0,T2和T3各林层的竞争表现与2013年类似,只有T1的总竞争强度表现为中林层>上林层>下林层。2013年和2018年的Hegyi_m均呈现为下林层>中林层>上林层。Hegyi变化方面:T0下林层降低53%,而中、上林层分别提高12%和14%,T1下林层大幅降低(86%),但中上林层升高,T2,T3上林层分别下降了19%和38%,下林层也呈降低趋势,中林层均有所提高。就Hegyi_m的变化来说,T0下林层和上林层分别下降了9%和8%,中林层提高14%,T1下林层和中林层有所上升,上林层略有降低。T2和T3中林层和上林层下降,下林层上升。 -
2018年优势树种的竞争较2013年发生了明显变化(表4)。4种方式下白桦和蒙古栎的Hegyi和Hegyi_m均有呈下降趋势,T2红松的Hegyi_m和T1,T3红松的Hegyi和Hegyi_m均呈下降趋势。T2,T3白桦的Hegyi和Hegyi_m下降幅度最大,T0下降幅度最小,T3红松的Hegyi和Hegyi_m下降幅度较大,T1对蒙古栎的竞争影响最小。
处理
Treatment树种
SpeciesHegyi Hegyi_m 2013 2018 2013 2018 T0 白桦
Betula platyphylla399.04 295.54 3.47 2.90 红松
Pinus koraiensis1 843.18 2 160.02 6.85 7.25 蒙古栎
Quercus mongolica1 799.72 1 310.16 3.27 2.74 T1 白桦
Betula platyphylla341.03 213.97 3.67 2.82 红松
Pinus koraiensis2 074.83 1 928.72 5.84 5.76 蒙古栎
Quercus mongolica2 648.66 2 239.84 3.55 3.22 T2 白桦
Betula platyphylla1 273.16 783.27 4.22 3.29 红松
Pinus koraiensis2 535.81 2 759.87 6.25 6.20 蒙古栎
Quercus mongolica2 506.09 1 952.00 3.97 3.48 T3 白桦
Betula platyphylla1 330.06 693.11 4.35 3.19 红松
Pinus koraiensis3 693.93 3 111.12 6.36 6.25 蒙古栎
Quercus mongolica3 574.54 2 608.78 4.02 3.34 Table 4. Changes of Hegyi and Hegeyi_m under different tending and thinning measures
-
4种方式下Hegyi_intra_m均呈下降趋势(图3),分别下降18%,21%,13%和15%,Hegyi_inter_m基本呈上升趋势(图3),但上升幅度较小,其中T1上升最明显(8%)。
由图4可知,蒙古栎的种间竞争树种主要是红松和白桦。T0方式下,两期(2013,2018)蒙古栎种间竞争树种的竞争强度及排序变化不大。T1方式下,蒙古栎竞争树种排序未发生变化,来自红松的竞争有所提升,来自白桦和大青杨的竞争有所减少,其他树种变化不大。T2方式下,来自红松和糠椴的竞争提高,来自白桦和大青杨的竞争降低。T3方式下,来自白桦的竞争减少,而来自红松的竞争提高。
红松主要的竞争树种是蒙古栎,T0红松的竞争树种强度排序变化不大,但来自落叶松和水曲柳的竞争强度增大,来自蒙古栎、糠椴和大青杨的竞争强度有所下降。T1主要竞争树种的强度排序发生部分变化,但强度变化不大。T2红松的竞争树种强度排序无变化,但来自白桦的竞争有所下降,来自色木槭和糠椴的竞争强度提高。T3主要竞争树种的排序变化不大,来自糠椴的竞争提高,排序靠前。对于白桦来说,主要的竞争树种为蒙古栎和红松,其次为色木槭和大青杨。T0白桦竞争树种强度和竞争排序无明显变化。T1方式下,来自色木槭的竞争强度和排序提升。T2方式下,来自大青杨的竞争强度下降,而来自色木槭的竞争强度升高。T3竞争排序不变,但来自色木槭和糠椴的竞争强度提高,来自大青杨的竞争强度明显下降。
综上所述,T1,T0方式下优势树种的竞争树种组成和强度变化不大且无明显趋势,而T2和T3方式下,各优势树种来自糠椴、色木槭的竞争强度增大,来自大青杨的竞争强度减小。
-
4种方式下,目标树的Hegyi_m均呈下降趋势,T2和T3方式下目标树的Hegyi_m和Hegyi_intra_m的下降幅度(14.0%,15.3%)明显高于T0和T1(4.4%,7.8%),4种方式下目标树的Hegyi_inter_m下降程度接近(表5)。目标树的Hegyi变化趋势与Hegyi_m变化趋势和幅度一致。
处理
Treatment树木类别
Type2013 2018 Hegyi_m Hegyi_intra_m Hegyi_inter_m Hegyi_m Hegyi_intra_m Hegyi_inter_m T0 非目标树non crop trees 5.75 2.89 2.86 5.32 2.38 2.94 目标树crop trees 1.81 0.85 0.95 1.73 0.80 0.93 T1 非目标树non crop trees 5.81 3.21 2.60 5.37 2.53 2.84 目标树crop trees 2.19 1.05 1.14 2.02 0.96 1.06 T2 非目标树non crop trees 5.94 2.68 3.26 5.71 2.35 3.36 目标树crop trees 2.21 1.12 1.09 1.90 0.83 1.07 T3 非目标树non crop trees 5.87 2.69 3.19 5.53 2.32 3.21 目标树crop trees 2.36 1.38 0.98 2.00 1.07 0.93 Table 5. Average competition index of crop trees and other trees to different thinning treatment
目标树的Hegyi_m均明显小于非目标树,T2和T3方式下目标树的Hegyi_m下降幅度分别为14%和15%,而非目标树仅降低了3.8%和5.8%。T0和T1中非目标树的种内竞争下降幅度要明显大于目标树的种内竞争下降幅度,而T2和T3呈现出相反的趋势。目标树和非目标树的种间竞争变化不大(表5)。
Response of Tree Competition in Natural Secondary Quercus mongolica Forest to Thinning Treatment
- Received Date: 2020-01-11
- Accepted Date: 2020-03-10
- Available Online: 2021-02-20
Abstract: