Processing math: 100%

内蒙古大兴安岭林区兴安落叶松人工林密度控制图研究

徐旭平, 吕延杰, 王建军

徐旭平, 吕延杰, 王建军. 内蒙古大兴安岭林区兴安落叶松人工林密度控制图研究[J]. 自然保护地,2023,3(4):79−88. DOI: 10.12335/2096-8981.2023070601
引用本文: 徐旭平, 吕延杰, 王建军. 内蒙古大兴安岭林区兴安落叶松人工林密度控制图研究[J]. 自然保护地,2023,3(4):79−88. DOI: 10.12335/2096-8981.2023070601
XU Xuping, LYU Yanjie, WANG Jianjun. Density Management Diagrams for Larix gmelinii Plantations in Daxing’an Mountains[J]. Natural Protected Areas, 2023, 3(4): 79−88. DOI: 10.12335/2096-8981.2023070601
Citation: XU Xuping, LYU Yanjie, WANG Jianjun. Density Management Diagrams for Larix gmelinii Plantations in Daxing’an Mountains[J]. Natural Protected Areas, 2023, 3(4): 79−88. DOI: 10.12335/2096-8981.2023070601

内蒙古大兴安岭林区兴安落叶松人工林密度控制图研究

详细信息
    通讯作者:

    王建军: E-mail:dreamjjwang@163.com

  • 中图分类号: S753

Density Management Diagrams for Larix gmelinii Plantations in Daxing’an Mountains

  • 摘要: 【目的】构建内蒙古大兴安岭地区兴安落叶松人工林收获预估密度控制图,为科学经营,合理调控林分密度、促进林木生长、提升蓄积收获量提供参考。【方法】以内蒙古大兴安岭林区的兴安落叶松人工林为研究对象,利用四期一类清查复测样地的142块样地数据,应用Reineke自然稀疏方程拟合生长曲线,构建林分密度控制图,并模拟无经营措施、合理间伐、重度疏伐3种模式下林分生长状态及收获量。【结果】Reineke自然稀疏方程拟合效果良好,其决定系数R2为0.911,均方根误差RMSE为0.008;优势高模型、蓄积模型2个辅助方程的拟合效果良好,其决定系数R2均在0.9以上;应用最优生长密度上限与下限、林分郁闭线以及2个辅助方程绘制了密度控制图。模拟结果表明,采用合理间伐、两次疏伐及达到目标胸径时,年平均收获量明显高于其他两种林分。同时,合理间伐林分达到目标胸径时的年龄、优势高均低于无经营措施林分,总收获量明显高于其他两种林分。【结论】大兴安岭林区兴安落叶松人工林密度控制图通过林分因子之间的关系直观展现了林分生长状态和抚育经营效果,并发现林分进入自然稀疏阶段后将密度控制在349~757株/hm2,能够提高达到目标胸径时林分蓄积总收获量。
    Abstract: 【Objective】A Density Management Diagram (DMD), which includes growth and yield estimation functions and is designed to facilitate reasonable silvicultural measures, has been developed for Larix gmelinii plantations in the Daxing’an Mountains. 【Method】The focus of this research was on the growth of Larix gmelinii plantations in Inner Mongolia Autonomous Region. Based on four times Chinese National Forest Inventory, 142 sample plots were collected. Reineke self-thinning equation was used to fit the growth curve and develop DMD. Based on DMD, the stand growth dynamic and yield was simulated under three different silvicultural measures: no management, reasonable thinning, and heavy thinning. 【Result】The Reineke self-thinning equation showed a good fitting effect with R2 of 0.911 and RMSE of 0.008. The dominant height, and volume models fitted as auxiliary equations, with R2 values above 0.9. The upper limit, lower limit, closed canopy lines, and two auxiliary equations were jointly drawn the DMD. Simulation results show that achieving the target stand average diameter at breast height (QMD) after two thinning resulted in a stand average yield of 19.25 m3/hm2 in the reasonably thinned stand, which was higher compared to the no-management stand, and the heavily thinned stand. The age and dominant height of the reasonably thinned stand in the target QMD were lower than that of the no-management stand, and the stand yield was significantly higher than the other two stands. 【Conclusion】The DMDs of Larix gmelinii plantations in Daxing’an Mountains can directly show the stand growth and the effectiveness of management practices through the relationship among stand parameters, and after self-thinning, stands with a density between 349 and 757 trees/hm2 could improve the yield of the stand under the target diameter.
  • 在生长过程中,林木对光、水分和空间等资源的竞争随着林木胸径的增加而增强,能力弱的自然枯死,林分进入自然稀疏状态(简称“自疏”)。该状态下,林木大小与单位面积林木株数之间表现出明显的线性关系[1-2]。 Yoda[3]和Reineke[4]认为无论树种、年龄和环境条件如何变化,林木大小−密度关系在对数轴上的斜率分别保持恒定值−1.500和−1.605。根据林木大小−密度关系进行林分密度调控,对实现特定经营目标具有一定的指导意义[5-6]

    林分密度控制图以图形方式展示林分密度与各项林分因子之间的数量关系,常用于抚育间伐[5,7-8]、防风防火[9]、资源调查[10]及造林设计等方面,对促进林分的生长、提高森林经营水平有重要意义[11-12]。田猛等[8]绘制了福建杉木(Cunninghamia lanceolata)人工林密度控制图,模拟了无措施、合理抚育间伐、对照抚育间伐3种不同经营措施下林分的生长情况。王雪等[5]利用林分密度控制图确定闽楠(Phoebe bournei)人工林各个发展阶段的目标树密度。Krefta等[13]利用邓恩桉(Eucalyptus dunnii)人工林密度控制图,提出达到林分胸径生长量进行间伐与达到自疏年龄进行间伐的2种森林经营方案。多项研究表明由于参数估计方法、树种、立地质量以及经营历史的不同,最大密度线的斜率值存在差异[14-16]。例如,孟京辉[15]通过普通最小二乘法回归(OLS)、简化主轴回归(RMA)、随机前沿函数(SFF)等回归方法拟合杉木人工林自疏方程,发现当随机误差项的方差趋于0时,SFF法的拟合结果最接近真实种群自疏过程。林分密度与林分生产力、林分结构、林下植被生长与分布、土壤水分与养分等紧密相关,且可以人为调控[17-20]。针对林分密度调控技术的研究,可以为森林经营提供合理林分密度,提高生产力,改善林分结构,实现特定经营目标。因此,构建树种区域性密度控制图,对当地林分生长研究、森林经营措施的制定有重要意义。

    兴安落叶松(Larix gmelinii)是内蒙古大兴安岭林区主要树种,水源涵养、空气净化等生态功能显著[21]。然而,兴安落叶松人工林经营较为粗放,相同的经营模式在不同经营目标、不同立地条件下同样采用的现象普遍存在。目前对兴安落叶松林的研究较多集中在生物量与碳储量[22-23]、林分结构[24-25]、土壤性质[26]等方面,针对兴安落叶松人工林密度调控技术的研究较少。本研究以内蒙古大兴安岭林区兴安落叶松人工林为研究对象,构建兴安落叶松人工林密度控制图,为合理调控密度以及制定科学经营措施提供参考。

    研究区域为内蒙古自治区东北部的大兴安岭林区,位于大兴安岭的西半部分,地理坐标为119°36′20″~125°20′50″E,46°08′40″~53°20′00″N,森林覆盖率为80.5%。寒温带大陆性季风气候,冬季寒冷漫长,夏季炎热多雨,年平均气温−2~4 ℃,极端高温37.5 ℃,极端低温−52 ℃,昼夜温差大,年均降雪厚度为20~40 cm,降雪期为每年10月到次年5月,年均降水量450 mm左右,全年降雨量较少。平均日照时长超过2 550 h,无霜期平均为80 d[27]。地带性森林植被是以兴安落叶松为主的针叶林,其他乔木树种有白桦(Betula platyphylla)、栎类(Quercus spp.)、杨树类(Populus spp.)等,土壤类型有棕色针叶林土、黑土、暗棕壤等。

    在第6次、第7次、第8次、第9次国家森林资源连续清查固定复样地中,筛选出兴安落叶松胸高断面积占比≥ 80%的人工纯林样地142块。数据包括树种、胸径、林木生长状态等样木信息,以及坡度、坡向、坡位、土壤厚度以及郁闭度等样地信息。样地面积为0.067 hm2,布设间距6 km × 6 km。样地数据特征量统计如表1所示。

    表  1  兴安落叶松人工林样地数据特征量统计
    Table  1.  Basic statistical information of sample plots of Larix gmelinii plantations
    林分变量最大值平均值最小值标准差
    海拔/m 1 160 839 390 194
    坡度/ (°) 25.00 4.92 1.00 5.26
    土壤厚度/cm 80.00 41.58 15.00 16.05
    腐殖质厚度/cm 15.00 5.39 2.00 3.05
    平均年龄/a 55 30 9 9
    平均树高/m 18.00 11.16 3.80 3.56
    郁闭度 95.00 67.86 20.00 22.37
    优势高/m 20.34 15.04 8.10 3.00
    林分平均胸径/cm 26.61 11.75 5.48 3.85
    公顷蓄积/ (m3·hm−2) 205.45 76.12 0.72 53.32
    林分密度/ (株·hm−2) 3657 1168 30 818
    公顷断面积/ (m2·hm−2) 35.89 12.74 0.18 8.57
    注:采用大兴安岭地区的兴安落叶松树高方程和材积方程分别计算兴安落叶松的树高和材积[28-29]。优势高表示为每公顷林分中胸径最大的前100株林木的平均高度[30]
    下载: 导出CSV 
    | 显示表格

    通过相对密度法筛选完满立木度样地,相对密度(relative density,缩写RD,式中记作RD)是指林分实际密度与其理论最大密度的比值(RDN/Nmax),RD > 0.7时林木出现枯死现象[31]。首先,根据Yoda[3]最大密度线方程,假设林分平均材积(记作VM)与林分密度(N)关系ln(VM)ln(N)的理论斜率为−1.5,找出ln(VM)+ln(N)最大的组合确定截距;然后,根据完满立木度最高的点代入求出每块样地的理论最大密度(Nmax),再计算相对密度(RD),相对密度RD > 0.7的样地为完满立木度样地。筛选得到完满立木度20块样地。

    选用Reineke自然稀疏方程[4]对完满立木度样地的林分密度与林分平均胸径(QMD,式中记作DM)进行拟合。为了降低参数估计产生的偏差,采用简化主轴回归(RMA)方法对方程进行参数估计,统计工具采用R语言的lmodel2包和smatr包。Reineke自然稀疏方程表达式如下:

    ln(N)=a1+1.605×ln(DM) (1)

    式中:N为林分密度,株/hm2DM为林分平均胸径,cm;a1为常数。

    密度控制图的绘制以林分株数密度为纵坐标、林分平均胸径为横坐标的双对数坐标系为基础,添加林分密度线和3个辅助方程等值线,绘图工具采用R语言的ggplot2包。

    林分密度线包括最大密度线、林分最优生长密度的上下限以及林分郁闭线。利用20块完满立木度样地数据拟合得到的Reineke方程反映了林木密度和林木大小的平均关系,将方程向上平移到截距的95%置信区间的上限,得到林分生长的最大密度线。林分最优密度指的是林木在竞争过程中能够充分地利用光照、水分、空间等资源,并且在竞争过程中不会因为资源缺乏而枯死的林分密度[8,31]。根据经验及相关研究结果,确定最大密度线的60%,30%和20%分别为最优密度的上限、下限和林分郁闭的限度值,分别保证了林木不会因为密度过高而出现激烈竞争,导致林木个体死亡,保证林木能够获得足够的生长空间,充分利用林内自然资源,保证林分的郁闭状态[8,32-33]

    此外,基于142块样地数据,采取非线性最小二乘法(NLS)回归方法拟合蓄积方程[式(2)]和优势高方程[式(3)][8],统计分析工具采用R语言的nls包。方程表达形式如下:

    M=a2×Nb2×DMc2; (2)
    HD=a3×Nb3×DMc3 (3)

    式中:M为林分蓄积;HD为优势木平均高;abc为模型参数。

    结合密度控制图,根据沈剑波等[34]构建的落叶松人工林多形地位指数方程,实现不同状态下林分年龄的预测。采用的立地指数方程如下:

    ISI=[HT×(1e0.0179×T0)1e0.0179×T]0.5550 (4)

    式中:ISI为地位级指数;T0为基准年龄,落叶松人工林的基准年龄为20 a;T为林分年龄;HT为在年龄T时的林分优势高。

    从142块兴安落叶松人工林样地筛选出20块RD > 0.70的完满立木度样地。通过RMA回归方法采用Reineke方程拟合林分密度,经计算方程的R2为0.911,RMSE为0.008。方程如下:

    ln(N)=11.9551.679ln(DM) (5)

    兴安落叶松人工林最大密度线[式(6)]、最优密度的上限[式(7)]、最优密度的下限[式(8)]和林分郁闭线[式(9)]4条平行线的方程如下:

    ln(N)=12.2291.679ln(DM); (6)
    ln(N)=11.7371.679ln(DM); (7)
    ln(N)=11.5361.679ln(DM); (8)
    ln(N)=11.3101.679ln(DM) (9)

    蓄积方程和优势高方程具体参数值和评价指标结果如表2所示,各方程的拟合效果良好,R2均在0.9以上。式(6)~(9)以及式(2)~(3)共同构建兴安落叶松人工林密度控制图(图1)。

    表  2  蓄积方程和优势高方程拟合参数结果和评价指标
    Table  2.  Parameter estimates and evaluation indicators of dominant height model and volume model
    方程参数估计值评价指标
    abcR2RMSE
    式(2)0.0001540.99662.44800.9991.837
    式(3)2.1909000.07620.57790.9320.773
    下载: 导出CSV 
    | 显示表格
    图  1  兴安落叶松人工林密度控制图
    Figure  1.  Density management diagram for Larix gmelinii plantations

    以密度控制图为基础,设计3种森林经营模式,通过模型计算各种不同经营措施下林分生长因子(图2表3),兴安落叶松人工林的经营目标胸径为40 cm,故以此为分析时限。图2中虚线a表示无经营措施的林分自然生长过程。经过计算,第22年开始出现自疏现象,此时林分平均胸径为19.3 cm,林分密度为1 424 株/hm2,随后超过了最优密度上限,蓄积量基本保持不变,林木株数减少。当林分平均胸径分别达到25,30,35和40 cm时,分别采伐504,747,901和1 006 株/hm2,蓄积量分别为367.00,423.00,476.19和528.09 m3/hm2,蓄积年均收获量分别为12.54,11.82,11.11和 10.40 m3/hm2

    图  2  兴安落叶松人工林密度控制图不同经营模式模拟情景
    Figure  2.  Simulations of different silvicultural measures of the density management diagram for Larix gmelinii plantations
    表  3  不同经营模式模拟结果计算
    Table  3.  Simulation results of different silvicultural measures
    模式疏伐方式时序林分平均
    胸径/cm
    密度/(株·hm−2)林龄
    /a
    优势高/m蓄积量/(m3·hm−2)总收获量/(m3·hm−2)年均收获量/ (m3·hm−2)
    无经营措施 郁闭到自疏点 之前 10.0 1 424 10 14.4 60.18
    之后 19.3 1 424 22 21.1 299.96
    自疏开始 之前 19.3 1 424 21.1 299.96
    皆伐(25 cm) 之后 25.0 920 29 23.7 367.00 366.92 12.54
    皆伐(30 cm) 之后 30.0 677 36 25.7 423.00 422.59 11.82
    皆伐(35 cm) 之后 35.0 523 43 27.5 476.19 476.19 11.11
    皆伐(40 cm) 之后 40.0 418 51 29.3 528.09 528.09 10.40
    合理间伐 郁闭到自疏点 之前 10 1 424 10 14.4 60.18
    之后 19.3 1424 22 21.1 299.96
    第一次疏伐 之前 19.3 1424 21.1 299.96
    之后 20.9 757 22 21.1 195.94 104.02
    皆伐(25 cm) 之后 25.0 757 28 23.3 302.10 406.12 14.37
    皆伐(30 cm) 之后 30.0 757 37 25.9 472.05 576.07 15.77
    第二次疏伐 之前 757 25.9 472.05
    之后 33.2 349 37 25.9 280.09 191.96
    皆伐(35 cm) 之后 35.0 349 39 26.7 318.21 614.19 15.56
    皆伐(40 cm) 之后 40.0 349 49 28.9 441.24 737.22 15.11
    重度疏伐 最优密度上限 之前 10.0 1424 10 14.4 60.18
    之后 19.3 1424 22 21.1 299.96
    重度疏伐 之后 12.4 285 10 14.4 20.34 39.84
    皆伐(25 cm) 之后 25.0 285 24 21.7 114.02 153.87 6.46
    皆伐(30 cm) 之后 30.0 285 30 24.1 178.17 218.01 7.17
    皆伐(35 cm) 之后 35.0 285 38 26.3 259.85 299.69 7.90
    皆伐(40 cm) 之后 40.0 285 47 28.4 360.32 400.16 8.57
    下载: 导出CSV 
    | 显示表格

    图2中虚线b表示采用合理间伐的林分生长过程,将林分密度控制在最优密度下限与上限之间,当林分生长轨迹达到最优密度线的上限时进行间伐,间伐后的林分密度在最优密度下限之上,间伐前林分密度在最优密度上限之下。第1次疏伐时的年龄为第22年,此时林分平均胸径为19.3 cm,林分密度为1424株/hm2,疏伐后保留密度为757株/hm2。当林分平均胸径为25和30 cm时,蓄积量分别为302.10和472.05 m3/hm2,蓄积年均收获量分别为14.37和15.77 m3/hm2。第2次疏伐在第37年,此时林分平均胸径为30 cm,密度为757株/hm2,疏伐后保留密度为349株/hm2。当林分平均胸径为35和40 cm时,蓄积量分别为318.21和441.24 m3,蓄积年均收获量分别为15.56和15.11 m3/hm2

    图2中虚线c是模拟重度疏伐的林木生长过程,为防止林木因生长竞争而自然枯死,采取一次性高强度采伐,经过估算采伐强度为80%。伐前林分平均胸径为10.0 cm,林分密度为1424株/hm2,伐后保留密度为285株/hm2。当林分平均胸径分别达到25,30,35和40 cm时,蓄积量分别为114.02,178.17,259.85和360.32 m3/hm2,蓄积年均收获量分别为6.46,7.17,7.90和8.57 m3/hm2

    表3可以看出:无经营措施林分平均胸径达到25 cm时具有最大的年均收获量,为12.54 m3/hm2,随后年均收获量越来越小;合理间伐林分的年均收获量持续增加;重度疏伐林分的年均收获量最小。当林分平均胸径达到30 cm时,合理间伐林分的年均收获量最大,为15.77 m3/hm2;其次是无经营措施林分,为11.82 m3/hm2;重度疏伐林分的最小,为7.17 m3/hm2

    应用密度控制图计算3种经营模式下林分达到目标胸径的林龄、林分平均胸径、优势高、密度、蓄积量和总收获量,结果如图3所示。图中CT1为第一次疏伐,CT2为第二次疏伐。采取合理间伐林分达到目标胸径25,35和40 cm时的林龄、优势高均低于无经营措施林分,达到目标胸径30 cm时的林龄、优势高与无经营措施相差不大。在达到目标胸径时,不同经营模式的蓄积量大小关系为:无经营措施 > 合理间伐 > 重度疏伐。达到胸径30 cm时,合理间伐林分的蓄积量和总收获量均高于无经营措施林分和重度疏伐林分;达到胸径35和40 cm时,合理间伐林分的蓄积量高于重度疏伐林分,其总收获量明显高于无经营措施林分和重度疏伐林分。

    图  3  不同经营模式林分生长因子对比图
    Figure  3.  Comparison of different silvicultural measures

    林分密度控制图能够直观展现不同发展阶段林分密度、林分平均胸径、林分年平均收获量等林分参数之间的数量关系,反映同龄林森林的抚育经营效果。通过合理间伐将林分密度控制在最优密度上限与下限之间,林分达到目标胸径时的年龄、优势高均低于无经营措施林分,总收获量明显高于无经营措施林分和重度疏伐林分。达到目标胸径30 cm时,合理间伐措施林分的年均收获量均高于无经营措施林分和重度疏伐林分,为15.77 m3/hm2。因此,在进入自然稀疏阶段将林分密度控制在349~757株/hm2,有助于林木生长,提高蓄积收获量,为内蒙古大兴安岭林区兴安落叶松人工林的合理林分密度调控与科学经营提供参考。

    林分密度控制图能够直观展示林分不同经营措施下林分参数状态,应用林分密度控制图可以发现间伐能够缩短林分达到目标胸径的时间。赵匡记等[35]认为林木胸径总生长量随间伐强度的增加而增大,张晓红等[36]认为影响单株目标树生长的决定因素是目标树的自由生长空间,这与应用兴安落叶松人工林密度控制图模拟的结果一致。本研究结果还表明间伐后林分达到目标胸径时的优势高低于无经营措施的林分,该优势高与采伐强度呈负相关。这一结果与陈哲等[37]研究结果表述不一致,其研究结果显示采取对照、22%、38%抚育间伐强度措施的落叶松人工林样地的树高年均生长量之间无显著性差异(P > 0.05)。这应该是由于达到目标胸径的时间缩短,间伐后树高生长时间短,优势高也就低,这与树高年均生长量无显著性差异是同向的,并不矛盾。

    不同经营模式模拟结果表明,基于最优密度的上限与下限,采用合理间伐措施有助于提高目标胸径下林分的收获量,这与多项基于林分密度控制图的应用研究结果一致。例如,田猛等[8]的研究结果表明,在最优密度范围内对福建杉木人工林进行合理抚育间伐,林分的年平均收获量较不采取任何经营措施、对照性抚育间伐表现出较高的趋势。Gyenge等[33]模拟了多枝桉(Eucalyptus viminalis)林分采取在林分达到自疏状态进行50%强度间伐作业、保持树冠覆盖度40%~50%的密度管理、无干预的林分生产力,结果表明采取50%强度间伐作业的林分能够获得较高的蓄积量和生物量。以最优密度上限、下限为基础,制定合理的采伐时间和采伐强度,林木能够充分利用空间、水、光等资源,从而提高林木总的收获量。

    兴安落叶松人工林在开始自疏后的林分最优生长密度区间为757~349株/hm2,经营性抚育采伐强度可达到50%。王雪等[5]利用闽楠人工林密度控制图确定目标树密度范围:质量选择阶段360~720株/hm2,目标树生长阶段165~330株/hm2,林分蓄积生长阶段65~130 株/hm2,经营性抚育采伐强度为50%。赵衍征等[38]的研究结果指出,35.5%左右的间伐强度有利于大兴安岭落叶松天然林林分结构稳定性和林分生产力的提高。李勃等[18]对比了1 312, 1 552和1 632 株/hm2兴安落叶松林分生长状况,发现胸径、树高、材积生长量随密度的增加而降低,兴安落叶松人工林的合理经营密度为1 312 株/hm2。该密度高于本研究的林分生长适宜密度,但是该研究表明林分生长随密度增加而减缓,说明林木未获得充足生长空间,林分的密度依然过高。因此,确定合理采伐强度,保留林分中生长良好的优势木,采伐林下小树或者干扰树,可以有效改善林分结构,促进林木生长[39]

    随着林业核心任务由木材生产到多功能经营的转变,单一的以提升蓄积量为目标的林分密度控制图已不能满足林业生产需求,因此建立以生态功能为主的多功能林分密度控制图是未来森林经营研究的方向和趋势之一。例如:马佳明[40]构建了油松(Pinus tabuliformis)人工林蓄积量、水源涵养、碳储量和生物量密度控制图,并应用密度控制图进行多功能分析;Yang等[41]将树种比例考虑到SDI方程中,构建了台湾地区松柏混交林密度控制图。为进一步提升我国森林多功能经营管理水平,构建多树种混交林的密度控制图也成为了我国森林经营的重点研究方向。

  • 图  1   兴安落叶松人工林密度控制图

    Figure  1.   Density management diagram for Larix gmelinii plantations

    图  2   兴安落叶松人工林密度控制图不同经营模式模拟情景

    Figure  2.   Simulations of different silvicultural measures of the density management diagram for Larix gmelinii plantations

    图  3   不同经营模式林分生长因子对比图

    Figure  3.   Comparison of different silvicultural measures

    表  1   兴安落叶松人工林样地数据特征量统计

    Table  1   Basic statistical information of sample plots of Larix gmelinii plantations

    林分变量最大值平均值最小值标准差
    海拔/m 1 160 839 390 194
    坡度/ (°) 25.00 4.92 1.00 5.26
    土壤厚度/cm 80.00 41.58 15.00 16.05
    腐殖质厚度/cm 15.00 5.39 2.00 3.05
    平均年龄/a 55 30 9 9
    平均树高/m 18.00 11.16 3.80 3.56
    郁闭度 95.00 67.86 20.00 22.37
    优势高/m 20.34 15.04 8.10 3.00
    林分平均胸径/cm 26.61 11.75 5.48 3.85
    公顷蓄积/ (m3·hm−2) 205.45 76.12 0.72 53.32
    林分密度/ (株·hm−2) 3657 1168 30 818
    公顷断面积/ (m2·hm−2) 35.89 12.74 0.18 8.57
    注:采用大兴安岭地区的兴安落叶松树高方程和材积方程分别计算兴安落叶松的树高和材积[28-29]。优势高表示为每公顷林分中胸径最大的前100株林木的平均高度[30]
    下载: 导出CSV

    表  2   蓄积方程和优势高方程拟合参数结果和评价指标

    Table  2   Parameter estimates and evaluation indicators of dominant height model and volume model

    方程参数估计值评价指标
    abcR2RMSE
    式(2)0.0001540.99662.44800.9991.837
    式(3)2.1909000.07620.57790.9320.773
    下载: 导出CSV

    表  3   不同经营模式模拟结果计算

    Table  3   Simulation results of different silvicultural measures

    模式疏伐方式时序林分平均
    胸径/cm
    密度/(株·hm−2)林龄
    /a
    优势高/m蓄积量/(m3·hm−2)总收获量/(m3·hm−2)年均收获量/ (m3·hm−2)
    无经营措施 郁闭到自疏点 之前 10.0 1 424 10 14.4 60.18
    之后 19.3 1 424 22 21.1 299.96
    自疏开始 之前 19.3 1 424 21.1 299.96
    皆伐(25 cm) 之后 25.0 920 29 23.7 367.00 366.92 12.54
    皆伐(30 cm) 之后 30.0 677 36 25.7 423.00 422.59 11.82
    皆伐(35 cm) 之后 35.0 523 43 27.5 476.19 476.19 11.11
    皆伐(40 cm) 之后 40.0 418 51 29.3 528.09 528.09 10.40
    合理间伐 郁闭到自疏点 之前 10 1 424 10 14.4 60.18
    之后 19.3 1424 22 21.1 299.96
    第一次疏伐 之前 19.3 1424 21.1 299.96
    之后 20.9 757 22 21.1 195.94 104.02
    皆伐(25 cm) 之后 25.0 757 28 23.3 302.10 406.12 14.37
    皆伐(30 cm) 之后 30.0 757 37 25.9 472.05 576.07 15.77
    第二次疏伐 之前 757 25.9 472.05
    之后 33.2 349 37 25.9 280.09 191.96
    皆伐(35 cm) 之后 35.0 349 39 26.7 318.21 614.19 15.56
    皆伐(40 cm) 之后 40.0 349 49 28.9 441.24 737.22 15.11
    重度疏伐 最优密度上限 之前 10.0 1424 10 14.4 60.18
    之后 19.3 1424 22 21.1 299.96
    重度疏伐 之后 12.4 285 10 14.4 20.34 39.84
    皆伐(25 cm) 之后 25.0 285 24 21.7 114.02 153.87 6.46
    皆伐(30 cm) 之后 30.0 285 30 24.1 178.17 218.01 7.17
    皆伐(35 cm) 之后 35.0 285 38 26.3 259.85 299.69 7.90
    皆伐(40 cm) 之后 40.0 285 47 28.4 360.32 400.16 8.57
    下载: 导出CSV
  • [1]

    Enquist B J, Brown J H, West G B. Allometric Scaling of Plant Energetics and Population Density[J]. Nature, 1998, 395(6698): 163−165. DOI: 10.1038/25977

    [2]

    Pretzsch H. A Unified Law of Spatial Allometry for Woody and Herbaceous Plants[J]. Plant Biology, 2002, 4(2): 159−166. DOI: 10.1055/s-2002-25732

    [3]

    Yoda K. Self-Thinning in Overcrowded Pure Stands under Cultivated and Natural Conditions (Intraspecific Competition among Higher Plants. XI)[J]. Journal of Biology Osaka City University, 1963, 14: 107−129.

    [4]

    Reineke L H. Perfecting a Stand-density Index for Even-aged forests[J]. Journal of Agricultural Research, 1963, 46(7): 627−638.

    [5] 王雪, 李际平, 曹小玉, 等. 基于林分密度控制图的闽楠人工林全周期经营目标树密度研究[J]. 中南林业科技大学学报, 2021, 41(9): 71−78. DOI: 10.14067/j.cnki.1673-923x.2021.09.008
    [6] 孙拥康, 汤景明, 王怡, 等. 鄂西山区日本落叶松人工林全林整体模型研究与应用[J]. 森林工程, 2023, 39(3): 57−63. DOI: 10.3969/j.issn.1006-8023.2023.03.007
    [7] 向玉国, 郑小贤, 刘波云, 等. 落叶松人工林生物量密度控制图的编制[J]. 中南林业科技大学学报, 2013, 33(10): 99−102. DOI: 10.3969/j.issn.1673-923X.2013.10.021
    [8] 田猛, 曾伟生, 孟京辉, 等. 福建杉木人工林密度控制图研制及应用[J]. 西北林学院学报, 2015, 30(3): 157−163. DOI: 10.3969/j.issn.1001-7461.2015.03.27
    [9] 王建军, 陆元昌, 赵秀海, 等. 马尾松人工林密度控制图及密度预警研究[J]. 北京林业大学学报, 2019, 41(5): 31−37. DOI: 10.13332/j.1000-1522.20190015
    [10] 向玉国, 郑小贤, 刘波云, 等. 福建将乐林场杉木碳储量密度控制图的编制[J]. 西北农林科技大学学报(自然科学版), 2014, 42(8): 99−104. DOI: 10.13207/j.cnki.jnwafu.2014.08.006
    [11]

    Newton P F. Stand Density Management Diagrams: Review of Their Development and Utility in Stand-Level Management Planning[J]. Forest Ecology and Management, 1997, 98(3): 251−265. DOI: 10.1016/S0378-1127(97)00086-8

    [12]

    Newton P F, Lei Y, Zhang S Y. Stand-Level Diameter Distribution Yield Model for Black Spruce Plantations[J]. Forest Ecology and Management, 2005, 209(3): 181−192. DOI: 10.1016/j.foreco.2005.01.020

    [13]

    Krefta S M, Hess A F, Atanazio K A, et al. A Density Management Diagram for Even-Aged Eucalyptus dunnii Stands in Southern Brazil[J]. Journal of Tropical Forest Science, 2021, 33(4): 406−413. DOI: 10.26525/jtfs2021.33.4.406

    [14]

    Ge F X, Zeng W S, Ma W, et al. Does the Slope of the Self-Thinning Line Remain a Constant Value across Different Site Qualities? —An Implication for Plantation Density Management[J]. Forests, 2017, 8(10): 355. DOI: 10.3390/f8100355

    [15] 孟京辉. 自然稀疏方程不同拟合方法的对比研究[J]. 北京林业大学学报, 2019, 41(12): 58−68. DOI: 10.12171/j.1000-1522.20190434
    [16]

    Madrigal-González J, Calatayud J, Ballesteros-Cánovas J A, et al. Global Patterns of Tree Density are Contingent Upon Local Determinants in the World’s Natural Forests[J]. Communications Biology, 2023, 6: 47. DOI: 10.1038/s42003-023-04419-8

    [17] 王冰, 张鹏杰, 张秋良. 不同林型兴安落叶松林土壤团聚体及其有机碳特征[J]. 南京林业大学学报(自然科学版), 2021, 45(3): 15−24.
    [18] 李勃, 德永军, 韦勤, 等. 3种密度对兴安落叶松人工林生长的影响[J]. 内蒙古农业大学学报(自然科学版), 2013, 34(2): 65−68. DOI: 10.16853/j.cnki.1009-3575.2013.02.013
    [19]

    Brunner A, Forrester D I. Tree Species Mixture Effects on Stem Growth Vary with Stand Density—an Analysis Based on Individual Tree Responses[J]. Forest Ecology and Management, 2020, 473: 118334. DOI: 10.1016/j.foreco.2020.118334

    [20] 董雪婷, 张静, 张志东, 等. 树种相互作用、林分密度和树木大小对华北落叶松生产力的影响[J]. 应用生态学报, 2021, 32(8): 2722−2728. DOI: 10.13287/j.1001-9332.202108.006
    [21] 李斯琴. 内蒙古大兴安岭森林资源问题分析及应对措施[J]. 广东蚕业, 2020, 54(8): 21−22. DOI: 10.3969/j.issn.2095-1205.2020.08.10
    [22] 鲍春生, 白艳, 青梅, 等. 兴安落叶松天然林生物生产力及碳储量研究[J]. 内蒙古农业大学学报(自然科学版), 2010, 31(2): 77−82.
    [23] 齐光, 王庆礼, 王新闯, 等. 大兴安岭林区兴安落叶松人工林植被碳贮量[J]. 应用生态学报, 2011, 22(2): 273−279. DOI: 10.13287/j.1001-9332.2011.0072
    [24] 王飞, 郭玉东, 张秋良. 兴安落叶松林生物量模型的构建[J]. 内蒙古农业大学学报(自然科学版), 2015, 36(4): 44−47. DOI: 10.16853/j.cnki.1009-3575.2015.04.008
    [25] 杨鑫, 王建军, 杜志, 等. 基于气候因子的兴安落叶松天然林单木直径生长模型构建[J]. 北京林业大学学报, 2022, 44(8): 1−11. DOI: 10.12171/j.1000-1522.20210353
    [26] 谷会岩, 金靖博, 陈祥伟, 等. 采伐干扰对大兴安岭北坡兴安落叶松林土壤化学性质的影响[J]. 土壤通报, 2009, 40(2): 272−275. DOI: 10.19336/j.cnki.trtb.2009.02.015
    [27] 王飞, 张秋良, 王冰, 等. 不同年龄杜香−兴安落叶松林粗木质残体贮量及特征[J]. 生态学杂志, 2012, 31(12): 2981−2989. DOI: 10.13292/j.1000-4890.2012.0413
    [28]

    Enzinga G Y, Jiang L C. Evaluation of Region and Subregion-Based Heightdiameter Models for Dahurian Larch (Larix gmelinii) in Daxing’an Mountains in China[J]. Applied Ecology and Environmental Research, 2019, 17(6): 13567−13591.

    [29] 中华人民共和国农林部. 标准立木材积表[S]. 北京: 技术标准出版社, 1978.
    [30]

    Sharma M, Amateis R L, Burkhart H E. Top Height Definition and Its Effect on Site Index Determination in Thinned and Unthinned Loblolly Pine Plantations[J]. Forest Ecology and Management, 2002, 168(1-3): 163−175. DOI: 10.1016/S0378-1127(01)00737-X

    [31]

    Solomon D S, Zhang L J. Maximum Size-Density Relationships for Mixed Softwoods in the Northeastern USA[J]. Forest Ecology and Management, 2002, 155(1/2/3): 163−170. DOI: 10.1016/S0378-1127(01)00556-4

    [32]

    Minoche D, Risio-Allione L, Herrero De Aza C, et al. Density Management Diagram for Teak Plantations in Tabasco, Mexico[J]. iForest-Biogeosciences and Forestry, 2017, 10(6): 909−915. DOI: 10.3832/ifor2247-010

    [33]

    Gyenge J, Lupi A N, Ferrere P, et al. Stand Density Management Diagrams of Eucalyptus viminalis: Predicting Stem Volume, Biomass and Canopy Cover for Different Production Purposes[J]. CERNE, 2019, 25(4): 463−472. DOI: 10.1590/01047760201925042666

    [34] 沈剑波, 雷相东, 雷渊才, 等. 长白落叶松人工林地位指数及立地形的比较研究[J]. 北京林业大学学报, 2018, 40(6): 1−8. DOI: 10.13332/j.1000-1522.20170400
    [35] 赵匡记, 纪福利, 刘延文, 等. 华北落叶松林分生长对间伐和修枝的响应[J]. 浙江农林大学学报, 2016, 33(4): 581−588. DOI: 10.11833/j.issn.2095-0756.2016.04.005
    [36] 张晓红, 张会儒, 卢军, 等. 目标树抚育间伐对蒙古栎天然次生林生长的初期影响[J]. 林业科学, 2020, 56(10): 83−92. DOI: 10.11707/j.1001-7488.20201009
    [37] 陈哲, 魏浩亮, 周庆营, 等. 抚育间伐对华北落叶松人工林林分结构的影响[J]. 中南林业科技大学学报, 2022, 42(5): 54−64. DOI: 10.14067/j.cnki.1673-923x.2022.05.006
    [38] 赵衍征, 李耀翔, 李春旭, 等. 间伐作业对天然次生林林分结构特征的影响[J]. 中南林业科技大学学报, 2021, 41(3): 83−95,110. DOI: 10.14067/j.cnki.1673-923x.2021.03.010
    [39] 林富成, 王维芳, 门秀莉, 等. 兴安落叶松人工林空间结构优化[J]. 北京林业大学学报, 2021, 43(4): 68−76.
    [40] 马佳明. 油松多功能人工林密度效应控制图的研制及应用[D]. 保定: 河北农业大学, 2021.
    [41]

    Yang T R, Lam T Y, Kershaw J A. Developing Relative Stand Density Index for Structurally Complex Mixed Species Cypress and Pine Forests[J]. Forest Ecology and Management, 2018, 409: 425−433. DOI: 10.1016/j.foreco.2017.11.043

  • 期刊类型引用(1)

    1. 刘林甫,洪光宇,弥宏卓,王志强,徐荣会,胡永宁. 大兴安岭不同林龄兴安落叶松生物量及固碳释氧功能. 林草资源研究. 2024(03): 88-95 . 百度学术

    其他类型引用(2)

图(3)  /  表(3)
计量
  • 文章访问数:  128
  • HTML全文浏览量:  37
  • PDF下载量:  30
  • 被引次数: 3
出版历程
  • 收稿日期:  2023-07-05
  • 修回日期:  2023-10-18
  • 网络出版日期:  2023-11-05
  • 刊出日期:  2023-12-17

目录

/

返回文章
返回