生态学报  2017, Vol. 37 Issue (20): 6668-6678

文章信息

范春雨, 张春雨, 赵秀海
FAN Chunyu, ZHANG Chunyu, ZHAO Xiuhai.
择伐对吉林蛟河阔叶红松林群落结构及动态的影响
Effects of selective harvest on community structure and dynamics in a mixed broadleaved Korean pine forest in Jiaohe, Jilin Province
生态学报. 2017, 37(20): 6668-6678
Acta Ecologica Sinica. 2017, 37(20): 6668-6678
http://dx.doi.org/10.5846/stxb201607291550

文章历史

收稿日期: 2016-07-29
网络出版日期: 2017-06-01
择伐对吉林蛟河阔叶红松林群落结构及动态的影响
范春雨 , 张春雨 , 赵秀海     
北京林业大学森林资源与生态系统过程北京市重点实验室, 北京 100083
摘要: 科学的森林经营能够优化林分结构,是调控森林生产力和生物多样性的有效手段。择伐作为森林经营的重要方式之一,其对森林结构以及群落动态的影响一直未有定论,因此迫切需要利用更加全面的数据对择伐及伐后林分特征的变化进行长期监测。根据森林大样地建立规范,2010年在吉林蛟河建立了42hm2阔叶红松林动态监测样地,2011年冬季截取部分面积进行择伐经营,以经营样地为研究对象,运用数值变量描述采伐活动并分析择伐前后群落结构的变化;同时结合2015年的二次调查数据,以立地条件基本一致的对照样地为参照,比较林分水平和物种水平上死亡率、更新率的差异,并利用线性混合效应模型探究择伐活动对个体径向生长的影响。研究结果显示:经营样地的择伐强度为5.4%,受采伐干扰影响较大的物种主要包括色木槭、白牛槭、裂叶榆、胡桃楸、千金榆、水曲柳以及紫椴,采伐主要集中于林冠层树种,亚林层和灌木层个体很少涉及。择伐前后物种组成、径级结构等并未发生明显改变。5年间,经营样地和对照样地的林分密度都降低,对比对照样地,经营样地的死亡率较低,但其更新状况并未优于对照样地。从胸高断面积来看,经营样地整体的年平均生长量高于对照样地,表明择伐导致的稀疏对个体生长和存活起到了一定的促进作用。将采伐强度纳入线性混合效应模型中分析发现,胸径始终是影响个体生长的最重要因素,其次是树木个体之间的非对称竞争;采伐所涉及到的7个主要树种的年平均生长量均高于对照样地,但仅有紫椴的径向生长表现出对采伐干扰的显著响应。综合来看,低强度择伐对群落结构和动态的影响较小,不同物种的径向生长对择伐的响应存在一定差异。
关键词: 择伐     林分结构     采伐偏好     动态     线性混合效应模型    
Effects of selective harvest on community structure and dynamics in a mixed broadleaved Korean pine forest in Jiaohe, Jilin Province
FAN Chunyu , ZHANG Chunyu , ZHAO Xiuhai     
Key Laboratory for Forest Resources & Ecosystem Processes of Beijing, Beijing Forestry University, Beijing 100083, China
Abstract: Forest management can be used to optimize stand structure, and it is an effective measure to regulate forest productivity and species diversity. The influence of selective harvest, which is an important method of forest management influencing the structure and dynamics of forest stands, does, however, not achieve consistent results. Accordingly, there is an urgent need to use more comprehensive data to study selective harvest and the change in forest stands after harvest. Following the protocols of the Center for Tropical Forest Science forest dynamic plot, we constructed a 42 hm2 plot in a mixed broadleaved Korean pine forest, located in Jiaohe, Jilin Province. We selected 19 ha of the plot to manage in the winter of 2011. The study was based on the managed plot, describing the harvest events using simple numeric variables and analyzing the change in the forest structure before and after the harvest event. At the same time, using the data of the second inventory of the plot after 5 years, we compared the mortality and recruitment between the managed plot and the control plot, which had the same habitat conditions at the stand and species levels. A linear mixed-effects model was employed to explore the effects of harvest on the growth of individuals. The results showed that the harvest intensity of the managed plot was 5.4%, as computed by the basal area. The species most influenced by the harvest included Acer mono, Ulmus laciniata, Juglans mandshurica, Carpinus cordata, Fraxinus mandshurica, and Tilia amurensis. Canopy tree species were harvested almost exclusively, whereas only a very few shrubs and sub-canopy trees were harvested. The harvest events did not noticeably change the species composition and diameter at breast height (DBH) distribution. During the 5 years, the stand density was lower than that of 2010 in both plots. The mortality of the managed plot was lower, whereas the recruitment status was not better, than those of the control plot. The basal area increment per year of the managed plot was larger than that of the control plot, indicating that the thinning resulting from harvesting promoted the growth of trees. We entered the predictors representing harvest intensity into the model and found that the DBH was the most significant variable for the growth analysis, followed by the asymmetric competition factor. The predictor representing harvest showed a significant effect only for the growth model of T. amurensis. In general, low harvest intensity had little effect on the structure and dynamics of the community, whereas the radial growth of different species showed various responses to selective harvest.
Key words: selective harvest     stand structure     harvest preference     dynamics     linear mixed-effects model    

森林是陆地生态系统的重要组成部分, 一直为人类提供生存和发展的重要资源, 但是随着近代对于森林的过度开发利用, 森林的质量在不断下降, 森林生态功能的发挥持续受阻, 这些现状使得森林经营的意义得到显现[1]。如今的森林经营开始朝着维持和恢复林分结构, 保证木材产量, 同时注重森林多种生态功能发挥的方向发展[2]

相较于人工林, 异龄混交林的经营存在着很多挑战。国际上对异龄混交林已经有许多种经营方式, 诸如北美的异龄林管理和欧洲的近自然经营等, 这些经营方法在不同的地区被广泛实践着[3]。森林采伐是森林经营的重要手段, “完美”的森林经营方案是希望在采伐这种人为干扰下, 林分仍能保持一种自然的状态, 不影响其生态功能的发挥, 也就是所谓的“近自然经营”, 但是这种所谓“自然”的状态也是在不断变化的, 因此是难以定义和度量当前林分的“自然”程度的。于是科学可持续的森林经营除了需要考虑采伐木材的限制性利用, 还要对采伐活动对森林生态系统结构的影响进行长期细致地监测。

择伐是对经营林地的一种十分重要的干扰, 大量研究表明, 相对于皆伐, 择伐是最优的森林经营策略[4-5]。Archambault利用加拿大魁北克南部伐后50a的混交林长期监测数据分析得到:择伐在一定程度上增加了保留木的径向生长量[6]。小范围适当强度的择伐会有效地改善群落内的光照条件, 进而促进森林结构的优化, 对植物的生长和更新产生一定的有利影响[7-9], 而且这种影响会在针叶树种和阔叶树种之间、同一物种不同径级的个体之间存在差异[10-11]。但是, 择伐在增加光照的同时也提高了环境的温度, 导致水分的大量蒸发, 一些植物个体会因为干旱的生活条件而死亡或者生长受到抑制[12]。同时, 择伐对于森林群落动态的影响也会随着研究区域的变化而变化[13]。国内开展的相关研究往往是短期的结果, 缺少对择伐及其影响的长期定位监测, 并且较少涉及到择伐对死亡率的影响研究, 造成对伐后森林生态系统结构功能变化的估计偏差。

阔叶红松林(Pinus koraiensis)是我国东北地区的地带性森林植被, 是温带针阔混交林的典型代表。自20世纪50年代, 我国开始对阔叶红松林进行择伐, 目的在于逐步恢复其原生植被, 同时保证森林经营的最大收益。了解择伐对于东北地区阔叶红松林群落结构及动态的影响对该区域森林的可持续发展具有十分重要的意义。

2010年在吉林蛟河的阔叶红松林内, 建立了立地条件相对一致的经营样地和对照样地。对照样地主要用来研究自然生态系统的结构和动态, 在经营样地内进行一定强度的择伐, 施加规律性的人为干扰。利用两块样地的对比, 分析森林经营影响下的林分结构和动态。本研究希望解决以下几个问题:(1) 择伐活动对于物种组成、径级结构的影响如何;(2) 在择伐的影响下, 林分的死亡、更新及径向生长如何变化, 不同物种之间是否存在差异。

1 材料与方法 1.1 研究区概况

研究区位于吉林省蛟河林业实验区管理局, 该区域属于受季风影响的温带大陆性山地气候。年均温为3.8 ℃, 最热月7月日均温为21.7 ℃, 最冷月1月日均温为-18.6 ℃, 平均年降水量为695.9 mm。土壤类型为山地暗棕色森林土, 土层厚度20—100 cm。植被属于长白山植物区系, 植被类型属于典型的天然次生针阔混交林。组成林分的乔木树种主要包括红松(Pinus koraiensis)、沙松(Abies holophylla)等针叶树种和千金榆(Carpinus cordata)、胡桃楸(Juglans mandshurica)、水曲柳(Fraxinus mandshurica)、紫椴(Tilia amurensis)、白桦(Betula platyphylla)、白牛槭(Acer mandshuricum)、色木槭(A. mono)、蒙古栎(Quercus mongolica)、春榆(Ulmus japonica)、裂叶榆(U. laciniata)等阔叶树种, 灌木主要包括暴马丁香(Syringa reticulata)、毛榛(Corylus mandshurica)、瘤枝卫矛(Euonymus pauciflorus)和东北鼠李(Rhamnus davurica)等[14]

1.2 研究方法

2010年夏季在研究区阔叶红松林内建立42 hm2(500 m×840 m; 43°57.524′— 43°58.042′N, 127°44.111′— 127°44.667′E)固定监测样地。参照CTFS (Centre for Tropical Forest Science)样地建设技术规范, 在样地建立后, 将其划分为20 m×20 m的连续样方。为便于测量每个个体的相对位置, 将每个20 m×20 m样方进一步划分成16个5m×5m的亚样方, 记录样地内所有胸径DBH≥1 cm的全部个体物种名、胸径、树高、冠幅(东西冠幅长、南北冠幅长)及枝下高, 并挂牌标记。对每个个体定位的方式是测量植株到5 m×5 m亚样方4个顶点的距离, 换算为其在整个样地内的相对坐标。同时, 测量每个样方的林冠开阔度和土壤养分。其中, 林冠开阔度是于2012年8月阴天, 用冠层分析仪(WinSCANOPY, Quebec, Canada)在每个样方中心距离地面1.5 m处拍摄1张半球状照片, 再用对应的WinSCANOPY和XLScanopy软件对照片进行处理而得到的。土壤养分是于2012年在每个样方采集500 g的土壤样品, 带回实验室进行分析, 测定其全氮、全磷、全钾、有机碳、pH值和速效氮、磷、钾含量。

由于样地内的地形条件(海拔、坡度和坡向等)和物种组成基本一致, 在样地内截取19 hm2进行森林经营实验。2011年冬季, 对经营样地进行了择伐。与传统择伐方式不同, 择伐前样地不预设采伐强度, 而是制定了采伐木和保留木原则, 并对每一个采伐个体进行记录以实现长期监测。总的采伐原则是在保护生物多样性的前提下, 充分发挥立地生产力, 使林分经济产出最高。具体包括, 采伐木的确定:① 优先伐除病腐木、濒死木, 双杈多梢木、弯曲木、被压木、干形不良木, 以及影响目的树种生长的非目的树种;目的树种包括红松、云冷杉、水曲柳、椴树、胡桃楸、柞树、黄菠萝等经济价值较高的树种。② 生长在一起的多株林木应伐除非目的树种, 对于从根部分叉的林木应保留一枝干型好的。保留木的确定:① 对珍稀树种红松、黄菠萝、山槐等应全部保留;② 保留具有非木质高产出的经济树种, 如红松、紫椴等;③ 保留林分中的大树;④ 保留能够给野生动物提供冬季休眠场所的林木;⑤ 保留大的奇形怪状的、可能有潜在艺术价值的林木;⑥ 注意保护树种多样性, 对于林分中个体很少的非目的树种应予以保留;⑦ 保留特殊的经济树种等。

1.2.1 择伐活动分析

温带森林往往呈现出层次分明的垂直结构, 位于林冠层的树木个体对光照的利用最充分, 其次是亚林层和灌木层。本研究将调查到的所有木本植物个体划分为3个物种组:林冠层、亚林冠层和灌木层, 考察不同种组的择伐行为。

筛选样地中采伐最多的7个物种进行分析, 分别为:白牛槭、胡桃楸、裂叶榆、千金榆、色木槭、水曲柳和紫椴, 它们在择伐中所占的比例分别为12%, 6%, 6%, 10%, 40%, 4%和8%。每个物种的采伐强度rG用伐除该物种植株的胸高断面积比例表示, 即rG=伐除植株胸高断面积/该物种总胸高断面积, 用rN表示伐除植株比例, 则NG可以用来描述采伐活动对树木大小的选择偏好:

NG越大, 则表明采伐对象更倾向于小树。同时分析经营样地中采伐活动对这7个主要树种径级分布的影响。

1.2.2 择伐对群落动态的影响

2015年, 对样地内全部树木个体进行了重复调查, 复查记录个体的生长状态和胸径, 标记出死亡的个体, 对新进阶的树木(胸径≥1 cm)进行定位、鉴定和挂牌。以经营样地和对照样地2010年和2015年的两次测量数据资料, 对比两块样地群落动态的异同, 计算物种水平和林分水平的死亡率和更新率, 分析经营活动对群落动态的影响。

考虑到树木个体生长的随机效应, 利用线性混合效应模型(LMMs)分析经营样地短期内影响树木个体径向生长的主要因素。以往此方面的研究主要考虑的自变量包括个体的胸径(DBH), 个体所处的环境条件(包括生物和非生物环境)[15-16]。本研究中, 鉴于个体胸高断面积生长量与胸径之间的非线性关系, 加入DBH以及中心化(去除共线性)后DBH的平方项。由于样地内地形平坦, 故本研究中非生物因素只考虑样地内的土壤理化性质和林冠开阔度。

代表生物因素的自变量表征了个体所处的竞争环境, 非对称竞争在植物群落邻体关系中普遍存在, 植物在对空间、光照、水分和土壤营养等资源的竞争过程中, 资源在竞争植物之间进行不平等的分配。大植株相对于小植株始终保持着一定的竞争优势[17-18]。本研究中使用目标个体20 m邻域范围内(去除边缘效应)胸高断面积(BA)以及胸径大于目标个体的胸高断面积比例(relBAL)代表邻体竞争强度。

式中,i代表目标个体, j代表i个体周围20m范围内胸径大于i的个体。

为衡量采伐对树木生长的影响, 在模型中加入了2个表征采伐的指标, 分别为个体20 m邻域范围内采伐植株的平均胸径与全部个体的平均胸径之比(d/D)和采伐后保留的胸高断面积比例(residual basal area removed)。

利用条件系数检验法检验k值 < 100, 各个自变量之间不存在多重共线性。

模型的响应变量为个体每年的平均胸高断面积增长量(cm2):

为了更加准确反映物种之间的生长差异, 同时采伐活动所涉及到的物种也是有限的, 选取了7个主要树种进行物种水平的逐个模型分析。为了满足模型对响应变量正态性的要求, 对胸高断面积的增量进行Box-Cox转换。同时针对每个模型计算了决定系数R2的值以及AIC值。由于混合效应模型的特殊性, 选取Nakagawa & Schielzeth[19]在2013年提出的边际R2(marginal R2)值即固定效应所能解释的模型变异比例来计算。数据分析与作图通过国际通用软件R-3.2.4及Excel中完成。

2 结果 2.1 择伐对林分结构的影响

2010年经营样地内胸径≥1 cm的个体数为20696, 共计35个物种, 隶属于17科24属。对比采伐前后, 物种组成的变化很小, 但是择伐在一定程度上改变了部分树种的相对优势程度。例如色木槭、紫椴和胡桃楸等物种, 择伐使其相对优势度下降, 红松、水曲柳和千金榆等的相对优势度上升(表 1)。

表 1 经营样地择伐前后物种组成分析 Table1 Composition of speices in managed plot before and after harvest
物种
Species
2010(伐前)pre-harvest2011(伐后)Post-harvest
平均胸径/cm
Average DBH
相对优势度/%
Relative superiority
平均胸径/cm
Average DBH
相对优势度/%
Relative superiority
色木槭Acer mono14.4717.8813.3816.83
Pinus koraiensis21.8712.8621.7713.27
紫椴Tilia amurensis24.5412.1923.6811.98
曲柳Fraxinus mandshurica28.7911.7528.2811.85
胡桃楸Juglans mandshurica28.4210.6027.3010.35
白牛槭Acer mandshuricum10.897.1310.457.21
千金榆Carpinus cordata7.584.557.445.12
沙松Abies holophylla23.144.2923.104.44
裂叶榆Ulmus laciniata11.343.4210.853.46
蒙古栎Quercus mongolica24.603.1424.223.21
枫桦Betula costata25.482.8323.822.70
春榆Ulmus davidiana var. Japonica12.632.3812.072.35
大果榆Ulmus macrocarpa18.391.2518.111.28
糠椴Tilia mandshurica17.251.0816.961.11
水榆花楸Sorbus alnifolia13.331.0512.410.98
拧筋槭Acer triflorum19.411.0418.651.00
香杨Populus koreana54.800.7854.800.80
黄檗Phellodendron amurense19.630.4919.420.50
青楷槭Acer tegmentosum9.250.369.200.39
青杨Populus cathayana51.550.3051.550.30
白桦Betula platyphylla18.730.2018.730.21
暴马丁香Syring areticulata4.700.164.670.31
山杨Populus davidiana25.450.1322.910.13
稠李Padus racemose5.790.135.790.17
簇毛槭Acer barbinerve3.270.033.270.08
怀槐Maackia amurensis17.870.0217.870.02
鼠李Rhamnus davurica6.470.026.470.02
瘤枝卫矛Euonymus verrucosus3.250.013.130.01
毛榛Corylus mandshurica2.750.012.750.03
山丁子Malus baccata5.490.014.350.01
翅卫矛Euonymus macropterus3.420.003.370.01
花楷槭Acer ukurunduense5.350.005.350.00
茶条槭Acer ginnala3.500.003.500.00
东北山梅花Philadelphus schrenkii1.940.001.940.00
花曲柳Fraxinus rhynchophylla4.000.004.000.00
金银忍冬Lonicera maackii2.880.002.880.00
小花溲疏Deutzia parviflora5.000.005.000.00
总计Total-100-100
相对优势度:物种胸高断面积/总胸高断面积

以胸高断面积计算, 经营样地的择伐强度为5.4%, 属于较低强度的择伐。样地内几乎所有的采伐个体均属于林冠层物种, 仅有很少的灌木层和亚林层个体被采伐。采伐活动对林冠层物种的影响最大, 其中林冠层个体的胸高断面积减少量占全部减少量的99%。

在择伐的主要树种中, 色木槭的择伐强度最高(10%), 平均每hm2的采伐株树为7.1, 采伐个体的径级分布范围约为15—35 cm (表 2)。相对于rG值, 各个树种的rN值更低, 仅有0.01-0.04, 两个指标分别从胸高断面积和株树密度的角度计量了采伐强度。从水曲柳和紫椴的NG值可以看出, 被伐个体主要是中等径级(20—40 cm)个体。而千金榆相比较而言, 小径级个体占绝大多数, 采伐的胸高断面积比例较高, 少量千金榆的大径级个体被伐除, 故NG值较低。

表 2 经营样地中主要树种的采伐强度和大小选择偏好分析 Table2 The selective harvest events for the major tree species in manage plot
主要树种
Major species
N/
(stems/hm2)
NA/
(stems/hm2)
NH/
(stems/hm2)
GB/
(m2/hm2)
GA/
(m2/hm2)
GH/
(m2/hm2)
rG rNNG林层
Forest layers
水曲柳45.945.30.63.83.70.10.030.010.67林冠层
紫椴62.861.21.63.93.70.20.050.030.6林冠层
胡桃楸47.345.913.43.20.20.060.020.33林冠层
千金榆249.2247.221.61.50.10.060.010.17林冠层
白牛槭137.1134.52.62.42.20.20.080.020.25林冠层
裂叶榆56.154.71.41.110.10.090.020.22林冠层
色木槭188.9181.87.15.85.20.60.10.040.4林冠层
NB代表采伐前的株树密度;NA代表采伐后的株树密度;NH代表采伐的株树密度;GB代表采伐前每公顷的胸高断面积;GA代表采伐后每公顷的胸高断面积;GH代表采伐掉的每公顷胸高断面积; 表中仅列出采伐强度(胸高断面积)前7位的物种

总体来看, 择伐活动并未明显改变主要树种的径级分布结构, 胡桃楸的径级分布呈现单峰型, 在20—30 cm处出现峰值;此外, 除水曲柳和紫椴呈现近似单峰型分布, 其余物种均为倒J型分布。择伐时主要伐除了对应树种径级分布的峰值位置(图 1)。从平均胸径来看, 采伐后, 绝大多数树种的平均胸径改变并不明显(表 1)。

图 1 经营样地中采伐前后径级分布结构 Fig. 1 The DBH distributions before and after harvest in the managed plot 灰色柱状图代表采伐后的株树, 灰色柱体上部白色部分代表采伐株树
3.2 择伐对群落动态的影响 3.2.1 择伐对死亡率和更新率的影响

对照样地个体数共计30095, 隶属于17科25属37种。从物种组成来看, 除东北鼠李和黑樱桃以外, 对照样地和经营样地的物种组成基本一致。其中, 红松、胡桃楸、色木槭、水曲柳和紫椴的相对优势度均大于10%, 是阔叶红松林的主要组成物种。

在林分水平上, 对照样地和经营样地的更新率均小于死亡率(表 3), 林分密度五年间都减小。经营样地的死亡率要低于对照样地, 表明择伐导致的稀疏使个体更好地存活。但是经营样地的更新低于对照样地, 说明样地内低强度的择伐并未对个体的更新起到明显的促进作用。

表 3 5年间经营样地和对照样地动态变化 Table3 Dynamics of the control and managed plot in 5 years
样地
Plot
采伐前株树密度
Stem density
pre-harvest/
(株/hm2)
采伐后株树密度
Stem density
post-harvest/
(株/hm2)
2015年株树密度
Stem density
in 2015/
(株/hm2)
更新率/%
Recruitment
死亡率
Mortality/%
年平均生长量
Average basal
area increment/
(BAI cm2/a)
经营样地1089.31068.9971.66.717.57.89
对照样地1308.51186.08.918.36.24
更新率和死亡率均为更新/死亡株树与2010年株树之比

在物种水平上, 灌木层物种如簇毛槭、毛榛等死亡量较大, 如簇毛槭在经营样地和对照样地中死亡个体占全部个体的比例分别达到了50%和44%。将全部个体按照胸径大小划分为不同的径级可以发现, 死亡主要发生在小径级个体中。从更新情况看, 更新与物种的多度之间有一定的相关性, 多度大的物种更新量也较大。同时, 死亡量大的物种更新量也会相对更大, 二者之间存在显著的正相关关系(r=0.79***)。在两块样地中, 簇毛槭的死亡量均小于更新量, 种群密度在5年内实现增大, 经营样地中, 毛榛的更新量大于死亡量, 除此之外, 所有物种种群密度都下降。对比两块样地, 采伐的主要树种也并未出现更新数量的大幅度提升(表 4)。

表 4 经营样地和对照样地物种水平动态统计 Table4 Dynamics of species in control and managed plot
物种
Species
对照样地control plot经营样地manage plot
多度2010
Abundance
in 2010
多度2015
Abundance
in 2015
死亡量
Number of
dead stems
更新量
Number of
recruitment
多度2010
Abundance
in 2010
多度2015
Abundance
in 2015
死亡量
Number of
dead stems
更新量
Number of
recruitment
色木槭5459506378238635903138648196
千金榆4726465040633047354416582263
白牛槭3416335652246226052478388261
暴马丁香1989132578311978047632319
紫椴15521369205221194104217119
簇毛槭*14531518634699383473194284
水曲柳133212311065872795781
红松1228115598251287116414118
胡桃楸12161155621898811881
春榆116110061984364351415425
裂叶榆98291615387106695020791
拧筋槭75667110318140135116
毛榛*62060036534517519799121
蒙古栎585533564264242231
鼠李43740260253119120
水榆花楸4353766012131973216
稠李405282165422321301119
糠椴390393532215166512
青楷槭3182649036252160986
枫桦300259410261215460
白桦272249241352870
大果榆205184265175147357
怀槐1961445424220
沙松112101165385367224
黄檗104891618170110
黑樱桃7756276
山丁子6866108201370
瘤枝卫矛563230652402614
山杨5443110121020
翅卫矛40153165028319
花楷槭40291659641
香杨343310181530
青杨2521408710
茶条槭2217721010
花曲柳119311100
小花溲疏109101100
东北鼠李6543
金银忍冬22115702
东北山梅花14033131
合计300952727855242707206961846136121377
*该物种的更新量大于死亡量, 种群密度增大
3.2.2 择伐对个体生长的影响

分析胸高断面积的变化, 两块样地内胸高断面积(BA)均有所上升, 经营样地BA净增加1.12m2/hm2, 对照样地净增加1.37m2/hm2, 主要是由于原有个体的生长造成的, 仅有少部分是由于更新个体。但经营样地的5年内平均每株的胸高断面积增长量要高于对照样地, 这表明采伐对于经营样地内个体的生长起到了一定的促进作用。同时, 分析主要采伐树种的生长也发现了一致的结果(图 3), 色木槭、水曲柳和紫椴在经营样地中的生长量要显著高于对照样地, 但胡桃楸和裂叶榆在两块样地中的生长差异不大。

图 2 对照样地和经营样地的胸高断面积变化 Fig. 2 The change of the basal area of control and managed plot

图 3 经营样地和对照样地的胸高断面积生长量比较 Fig. 3 The comparison of basal area increment between the control and managed plot

植物个体的生长与个体大小及邻域环境条件密切相关, 运用线性混合效应模型分析经营样地的胸高断面积生长量与胸径、竞争、非生物环境因子和采伐活动的影响。模型中响应变量和自变量的取值范围见表 5

表 5 线性混合效应模型响应变量和部分自变量的取值范围 Table5 Summary statistics for response variables and a subset of growth predictors of the LMMs
物种
Species
BAI/(cm2/a)DBH/cmBA/(cm2/hm2)relBALd/D
meanSDmeanSDmeanSDmeanSDmeanSD
白牛槭4.679.729.849.7719.136.270.920.132.351.46
胡桃楸22.7218.7428.99.8519.056.780.630.231.991.32
裂叶榆5.3110.8410.6010.8418.766.760.900.172.171.41
千金榆3.486.997.265.0019.016.230.960.062.121.13
色木槭7.615.9013.1913.2919.146.490.850.222.161.29
水曲柳24.9722.1929.0813.6719.366.260.610.292.051.15
紫椴16.1419.3925.9513.1519.246.410.690.271.951.20

表 6 线性混合效应模型显著影响因子及模型评价统计量 Table6 The significant predictors and the model-fit evaluation statistics for the LMMs
物种Species显著影响因子Significant factorsR2AIC
白牛槭DBH DBH2 K relBAL0.468035.1
胡桃楸DBH DBH2 relBAL0.293183.3
裂叶榆DBH DBH2 TK relBAL0.581371.8
千金榆DBH DBH2 TN TK C P relBAL0.305913.9
色木槭DBH DBH2 relBAL TK K0.564350.8
水曲柳DBH DBH2 relBAL0.56997.8
紫椴DBH DBH2 d/D rel relBAL0.451424.9

针对每个树种逐个进行模型分析发现, 所选取的固定效应因子能够解释个体生长29%—58%的变异, 其中胸径是影响个体生长的主要因素。考虑到生长和胸径之间的非线性关系, 加入了去中心化的DBH2, 在一定程度上提高了模型的解释能力。在生物环境变量中, 竞争显著影响个体生长。目标个体周围相对更大的个体的胸高断面积和(relBAL)表现出与生长量之间显著的负相关关系, 并未检测到个体周围20 m范围内全部邻体的胸高断面积和(BA)和与生长的显著相关关系, 这代表群落中存在着明显的非对称竞争, 资源在不同径级个体之间的分配不均导致小个体的生长受到抑制。在模型所考虑的非生物环境变量中, K元素含量相对其他元素与样方内个体的生长表现出更强的相关性。但综合来看, 生物变量对个体生长的影响更大。尽管经营样地7个主要树种的年均胸高断面积生长量均高于对照样地, 但在物种水平的生长模型中加入与采伐相关的自变量后, 结果显示仅有紫椴对其周围的采伐行为有一定的生长响应, 邻域内采伐平均胸径增加, 采伐强度越高时, 紫椴的生长量越大, 其余物种均未呈现出个体生长与采伐变量的显著相关性。

4 讨论 4.1 择伐对林分结构的影响

由于近年来不合理的经营利用, 当前东北阔叶红松林普遍存在着过伐林林分组成简单、稳定性较差的问题, 无法充分利用林地的生产力, 因此迫切需要我们通过调整林分结构来改善阔叶红松林的现状[20]。择伐是优化林分结构的重要手段, 能够有效地改善由于过于稠密造成的生长抑制、病虫害蔓延等[21], 同时采伐作业时主要伐除处于劣势生态位的个体, 对于物种多样性和生态系统稳定性也会产生一定影响[21]。本研究中并未发现经营样地采伐后林分结构的明显改变, 因为低强度的择伐对林分的扰动较小, 群落结构能够较快恢复, 这与前人在低强度采伐样地中的结果相一致[3]。Kyle等在2012年研究美国俄勒冈州西部的幼龄林对于不同强度择伐的响应时发现, 短期内, 择伐并未加速林分结构朝着演替后期的方向发展[23]。雷相东等的研究表明20%—30%强度的择伐没有显著改变落叶松云冷杉林的物种多样性和结构[24]。由此看来, 采伐行为虽然可以使用强度进行量化, 但始终是模糊且受到多种标准的影响的, 难以避免采伐决策制定过程中的主观性[25, 3], 应该尽量利用全面准确的定位监测数据对采伐行为对林分的多方面影响进行综合评价, 以帮助我们制定合理的经营措施和政策。

4.2 择伐对群落动态的影响

在异龄混交林这样一个复杂的生态系统中, 存在着年龄、树种和径级都不同的树木个体[26, 8]。由于林内各活立木对生长空间和资源的相互竞争及对各种外界环境因素的不同反应, 异龄林内的动态是难以准确捕捉其规律的。研究者们逐渐开始倾向于使用新的分析手段挖掘动态数据, 诸如利用混合效应模型探究生长或死亡与个体大小、竞争以及环境因素等的关系[27];使用矩阵模型模拟群落动态变化等[28], 尽可能全面考虑森林群落中可能影响动态变化的干扰因素。

择伐作为一种重要的人为干扰, 在降低林分密度的同时, 改变着林内温度、光照、竞争、特定生境条件和土壤等生态系统结构要素[12]。采伐过后, 林分密度降低, 林下拥挤程度得到缓解, 树木之间的竞争减弱, 提高存活率[29];同时, 择伐增加了植物充分进行光合作用的时间, 使得喜光树种具备更大的生长潜力, 可利用资源的增多也会加快林下的更新。在我们的研究中, 经营样地内的死亡率低于对照样地, 但更新状况并未得到明显的改善, 表明尽管择伐能够在一定程度上优化林分结构, 但是天然林分内物种的更新对非生物环境条件和邻体竞争更为敏感[11], 低强度择伐不能改变乔木树种对于绝大部分资源的竞争优势, 这些树种的幼树仍然难以获得足够的资源和生活空间。大部分研究都证实择伐会加快林分的生长, 随着采伐强度的提高, 保留木的直径生长往往也随之增加[30-32]。蒋子涵和金光泽进一步利用树木年轮分析发现择伐后保留木个体的生长会出现加速现象且能够持续26—29a[7], 在本文的研究中也发现了择伐对于林分水平胸高断面积生长的促进作用。

当前很多研究开始依赖基于个体水平的模型分析择伐对于群落结构和动态的影响[33-34], 在分析经营样地种不同物种生长影响因素时首次加入了表征采伐强度的指标, 但是没有发现采伐活动对个体生长的显著影响, 表明择伐虽然改变了林下的小环境, 但是保留木的生长还与采伐地的物种组成和地形条件等因素密切相关[27-28], 同时, 树木生长对择伐活动的响应可能存在一定的时间延迟(time-lagged growth responses)[35-36], 因此采伐对群落动态的影响可能需要较长的时间才能显现出来, 仅用短期数据很难发现规律。

当前对于森林群落内树木的生长、死亡和更新的研究已经十分广泛[35], 但是结合采伐行为, 将采伐强度指标纳入动态模型分析中的相关结果还很少见。森林的动态变化是自然变化与森林经营措施等干扰综合作用的结果。不考虑采伐干扰的动态模型是不完整的。在今后的研究中, 还应加入包含林分密度变化量和结构性变化的估计量, 更加全面的反映采伐对林分机构和动态的影响。

参考文献
[1] 吴蒙蒙, 王立海, 侯红亚, 刘德洲. 采伐强度对阔叶红松林空间结构的影响. 东北林业大学学报, 2013, 41(9): 6–9.
[2] 包也, 于大炮, 周莉, 周旺明, 孟莹莹, 吴志军, 郭焱, 王晓雨, 代力民. 择伐对阔叶红松林资源可持续利用的影响. 生态学报, 2015, 35(1): 31–37.
[3] Zhang C Y, Zhao X H, Gadow K V. Analyzing selective harvest events in three large forest observational studies in North Eastern China. Forest Ecology and Management, 2014, 316: 100–109. DOI:10.1016/j.foreco.2013.07.018
[4] 克劳斯·冯佳多, 惠刚盈. 森林生长与干扰模拟. 北京: 中国林业科学研究院, 1998.
[5] 汤孟平, 唐守正, 雷相东, 李希菲. 林分择伐空间结构优化模型研究. 林业科学, 2004, 40(5): 25–31. DOI:10.11707/j.1001-7488.20040504
[6] Archambault L, Delisle C, Larocque G R. Forest regeneration 50 years following partial cutting in mixedwood ecosystems of southern Quebec, Canada. Forest Ecology and Management, 2009, 257(2): 703–711. DOI:10.1016/j.foreco.2008.09.056
[7] 蒋子涵, 金光泽. 择伐对阔叶红松林主要树种径向与纵向生长的影响. 生态学报, 2010, 30(21): 5843–5852.
[8] 赵俊卉, 亢新刚, 龚直文. 择伐对北方森林更新、生物多样性和生长的影响研究进展. 内蒙古农业大学学报, 2008, 29(4): 264–270.
[9] Pamerleau-Couture É, Krause C, Pothier D, Weiskittel A. Effect of three partial cutting practices on stand structure and growth of residual black spruce trees in north-eastern Quebec. Forestry, 2015, 88(4): 471–483. DOI:10.1093/forestry/cpv017
[10] Prévost M, Dumais D, Pothier D. Growth and mortality following partial cutting in a trembling aspen-conifer stand: results after 10 years. Canadian Journal of Forest Research, 2010, 40(5): 894–903. DOI:10.1139/X10-036
[11] Deal R L, Tappeiner J C. The effects of partial cutting on stand structure and growth of western hemlock-Sitka spruce stands in southeast Alaska. Forest Ecology and Management, 2000, 159(3): 173–186.
[12] Bataineh M, Kenefic L, Weiskittel A, Wagner R, Brissette J. Influence of partial harvesting and site factors on the abundance and composition of natural regeneration in the Acadian Forest of Maine, USA. Forest Ecology and Management, 2013, 306: 96–106. DOI:10.1016/j.foreco.2013.06.016
[13] Ruel J C, Fortin D, Pothier D. Partial cutting in old-growth boreal stands: an integrated experiment. The Forestry Chronicle, 2013, 89(3): 360–369. DOI:10.5558/tfc2013-066
[14] Zhang C Y, Zhao Y Z, Zhao X H, von Gadow K. Species-habitat associations in a northern temperate forest in China. Silva Fennica, 2012, 46(4): 501–519.
[15] Liang J J, Buongiorno J, Monserud R A, Kruger E L, Zhou M. Effects of diversity of tree species and size on forest basal area growth, recruitment, and mortality. Forest Ecology and Management, 2007, 243(1): 116–127. DOI:10.1016/j.foreco.2007.02.028
[16] Porté A, Bartelink H H. Modelling mixed forest growth: a review of models for forest management. Ecological Modelling, 2002, 150(1/2): 141–188.
[17] Hara T. Dynamics of size structure in plant populations. Trends in Ecology & Evolution, 1988, 3(6): 129–133.
[18] Zhang C Y, Jin W B, Gao L S, Zhao X H. Scale dependent structuring of spatial diversity in two temperate forest communities. Forest Ecology and Management, 2014, 316: 110–116. DOI:10.1016/j.foreco.2013.07.025
[19] Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods in Ecology and Evolution, 2013, 4(2): 133–142. DOI:10.1111/mee3.2013.4.issue-2
[20] 胡云云, 闵志强, 高延, 冯启祥. 择伐对天然云冷杉林林分生长和结构的影响. 林业科学, 2011, 47(2): 15–24. DOI:10.11707/j.1001-7488.20110203
[21] 惠刚盈, 胡艳波, 赵中华. 再论"结构化森林经营". 世界林业研究, 2009, 22(1): 14–19.
[22] 沈林, 杨华, 亢新刚, 岳刚. 择伐强度对天然云冷杉林空间分布格局的影响. 中南林业科技大学学报, 2013, 33(1): 68–74.
[23] Dodson E K, Ares A, Puettmann K J. Early responses to thinning treatments designed to accelerate late successional forest structure in young coniferous stands of western Oregon, USA. Canadian Journal of Forest Research, 2012, 42(2): 345–355. DOI:10.1139/x11-188
[24] 雷相东, 陆元昌, 张会儒, 张则路, 陈晓光. 抚育间伐对落叶松云冷杉混交林的影响. 林业科学, 2005, 41(4): 78–84. DOI:10.11707/j.1001-7488.20050414
[25] Zucchini W, von Gadow K. Two indices of agreement among foresters selecting trees for thinning. Forest and Landscape Research, 1995, 1: 199–206.
[26] 王飞, 代力民, 邵国凡, 谷会岩, 惠淑荣, 熊晓波, 王惠. 非线性状态方程模拟异龄林径阶动态——以长白山阔叶红松林为例. 生态学杂志, 2004, 23(5): 101–105.
[27] Coomes D A, Allen R B. Effects of size, competition and altitude on tree growth. Journal of Ecology, 2007, 95(5): 1084–1097. DOI:10.1111/jec.2007.95.issue-5
[28] Kariuki M, Kooyman R M, Brooks L, Smith R G B, Vanclay J K. Modelling growth, recruitments and mortality to describe and simulate dynamics of subtropical rainforests following different levels of disturbance. Forest Biometry Modelling and Information Sciences, 2006, 1: 22–47.
[29] Burgess D, Robinson C, Wetzel S. Eastern white pine response to release 30 years after partial harvesting in pine mixedwood forests. Forest Ecology and Management, 2005, 209(1/2): 117–129.
[30] Simard S W, Blenner-Hassett H, Cameron I R. Pre-commercial thinning effects on growth, yield and mortality in even-aged paper birch stands in British Columbia. Forest Ecology and Management, 2004, 190(2/3): 163–178.
[31] Forget E, Nolet P, Doyon F, Delagrange S, Jardon Y. Ten-year response of northern hardwood stands to commercial selection cutting in southern Quebec, Canada. Forest Ecology and Management, 2007, 242(2/3): 764–775.
[32] Mäkinen H, Isomäki A. Thinning intensity and long-term changes in increment and stem form of Scots pine trees. Forest Ecology and Management, 2004, 203(1/3): 21–34.
[33] Rüger N, Gutiérrez Á G, Kissling W D, Armesto J J, Huth A. Ecological impacts of different harvesting scenarios for temperate evergreen rain forest in southern Chile——a simulation experiment. Forest Ecology and Management, 2007, 252(1/3): 52–66.
[34] Temesgen H, Martin P J, Maguire D A, Tappeiner J C. Quantifying effects of different levels of dispersed canopy tree retention on stocking and yield of the regeneration cohort. Forest Ecology and Management, 2006, 235(1/3): 44–53.
[35] Thorpe H C, Thomas S C, Caspersen J P. Residual-tree growth responses to partial stand harvest in the black spruce (Picea mariana) boreal forest. Canadian Journal of Forest Research, 2007, 37(9): 1563–1571. DOI:10.1139/X07-148
[36] Thorpe H C, Vanderwel M C, Fuller M M, Thomas S C, Caspersen J P. Modelling stand development after partial harvests: an empirically based, spatially explicit analysis for lowland black spruce. Ecological Modelling, 2010, 221(2): 256–267. DOI:10.1016/j.ecolmodel.2009.10.005