-
树冠结构属于森林结构的范畴,不仅决定林木形态、未来长势与活力[1],其变化还能降低林木复杂的竞争[2]。研究树冠结构能帮助理解林木利用生态空间的能力,揭示其生长对策和适应机制[3]。密集度是林分空间结构的重要属性,与林分密度有区别,通过判断林分空间结构单元中树冠连接来分析林木疏密程度,反映了树冠的密集程度与林冠层的连续遮盖程度和挤压状况[4]。密集度与树冠结构有关系。通常,密集度取值越大,林地微环境林木密集程度越高,树冠遮盖越连续,即林木竞争就越激烈;取值越小则说明林木密集程度越低,即竞争压力减小,说明此时林木生长稀疏,出现林隙的可能性增加。森林经营中,冠幅常用于表征树冠特征,要探究林木密集程度与树冠结构的关系,单一的树冠特征难以起到普遍性作用且意义不大。就整个林冠层挤压与林木竞争而言,林木密集程度很可能受到除冠幅外的其它树冠结构特征,如树冠面积、冠长和树冠体积等的影响。
针对树冠结构与林木密集程度关系的专门研究迄今极为少见,多数研究集中于树冠结构与林木竞争关系层面。早期研究表明,林分密度能够预测林冠层竞争[5],无论树种组成与林木年龄,树冠发育与胸径增长对整个林分的竞争都非常敏感[6]。尽管树冠结构与林木大小密切相关,但就混交林而言,冠幅因受林木径阶差异的影响,因此很难量化林分种间、种内竞争[7],这就说明树冠结构与林木竞争关系可能与林木大小差异有关。树冠结构的变化与林分密度有关。在分布较稠密的森林中,外力干扰(如大风侵袭)导致冠层出现林隙,树冠结构变的更加复杂[8]。因为风的频繁干扰,树冠脆弱部位容易折断,单木树冠完整性及林冠层均遭到不同程度的破坏,冠幅、冠长和冠幅面积等因此改变。一旦完整的树冠被破坏,树冠结构内部、树与树的挤压状态随之也被调整,林木密集程度可能相应发生变化,且这种变化很可能因林分密度的差异有所不同。对于未经干扰或干扰较轻的林分而言,树冠结构变化可能由林木自疏所致。Gause假说认为,对同一资源产生竞争的2个种无法长期共存,结果会导致1个种占优势,另1个种则被淘汰[9]。随着树的发育,个体对营养资源与生长空间的需求增加,导致林木竞争更加激烈[10],当光热资源不足以支持所有林木生长时,一些林木生长将减缓,而一些林木将最终死亡[11]。表明无论林木个体还是整个林分,树冠结构与密集度的关系都避免不了林木竞争的作用,它们的关系很可能与林分密度有关。为此,本研究将探究林分密度如何影响树冠结构与密集度的关系。
蒙古栎(Quercus mongolica Fisch. ex Ledeb.)林在小陇山林区北部属主要的次生混交林群落和典型的地带性顶级植被,对维持当地生态平衡、提高环境质量和保护生物多样性具有重要的生态价值。近年来,国内学者从优化森林结构的角度对小陇山蒙古栎林展开了大量研究,比如森林经营措施与林分空间结构[12-14]、群落结构与稳定性[15]等,而有关蒙古栎林树冠结构的研究相对缺乏。近期研究表明,小陇山蒙古栎林的密集度处于中等至密集状态,除与林木冠幅大小有关外,可能与林分密度也有关[15],尚需进一步论证。基于以上背景,本研究依托小陇山林区北部4个30 m × 30 m的蒙古栎次生混交林样地,分析不同密度林分树冠结构特征的差异,探究不同密度林分的林木胸径与树冠结构的关系,从而厘清林分密度与树冠结构及密集度的关系,以期深入掌握蒙古栎次生混交林结构特征,为蒙古栎林可持续经营和资源可持续利用提供更加充分的科学依据,也为发挥山地森林的多种生态效能提供理论基础。
-
样地位于甘肃省小陇山林业保护中心滩歌林场境内(104°40′~104°58′ E,34°27′~34°39′ N)。气候属温带大陆性气候,受季风影响较显著,具有冬冷无严寒、夏热无酷暑的特点。根据武山县气象站资料记载,本地年均温9.6 ℃,年降水量490 mm以上,相对湿度65%,降水集中于7—9月,约占全年降水量的70%,无霜期140~218 d。研究样地海拔在1 800~2 800 m,坡度为25°~40°,山地垂直带谱明显,土壤主要为棕色土,多石砾。研究区森林覆盖率达47.9%,除次生树种蒙古栎、复叶槭(Acer negundo L.)、椴树(Tilia tuan Szyszyl.)、华山松(Pinus armandii Franch.)、山杨(Populus davidiana Dode)和白桦(Betula platyphylla Suk.)外,人工栽培树种包括青海云杉(Picea crassifolia Kom.)、华北落叶松(Larix principis-rupprechtii Mayr)、油松(Pinus tabuliformis Carrière)等。
-
野外调查于2020年8月进行,选取立地条件典型且均一性较强、人为干扰较少的林分。检尺并记录胸径(DBH)>5 cm的活立木,用北斗手持导航仪GPS( 型号:G130BD)确定每株树的坐标。用皮尺、围尺、测高仪(型号:BERTHA DENTON-G系列)测定林木胸径、树高、活枝下高、冠幅(东西与南北方向)等林分因子,编号及记录每株树的健康状态(是否病害、倒伏、弯曲与直立)。按照5点取样法”在每个样地内取2 m × 2 m的小样方5个,调查幼树更新、灌木种类、盖度、高度、株数、基径等(表1)。
样地
Sample plot林分类型
Forest type林分密度/(株·hm−2)
Stand density平均胸径/cm
Mean DBH平均树高/m
Mean height蓄积量/m3
Stand volume胸高断面积/m2
Basal areaⅠ 蒙古栎-白桦混交林 1 675 10.3 7.2 51.67 0.092 Ⅱ 蒙古栎-白桦混交林 1 413 10.9 7.4 45.78 0.103 Ⅲ 蒙古栎-白桦混交林 1 300 10.9 7.7 42.56 0.103 Ⅳ 蒙古栎-白桦混交林 1 125 9.0 7.0 26.11 0.071 Table 1. Details of sample plot in different densities
-
采用冠幅(CD)、冠长(CL)、冠长率(CR)、冠形率(CS)、树冠面积(CPA)和树冠体积(CVOL)等6种指标表征树冠结构,其中,冠幅为树冠在东西与南北方向冠宽的均值;冠长为树高与活枝下高之差;冠长率为冠长与树高之比;冠形率为冠长与冠幅之比;树冠面积与树冠体积分别由冠幅和树冠投影计算得到。
-
利用重要值法[16]确定蒙古栎林的优势树种。依据惠刚盈等[4,14]对密集度(C)的表述与取值,从林木个体与林分2种尺度进行密集度分析,同时采用Hegyi指数[13]判别蒙古栎林个体与林分的竞争强度(表2)。
指标
Indicator尺度
Scale表达式
Formula意义
Explanation密集度 Crowding 林木个体 $C_{i}=\dfrac{1}{4} \sum\limits_{j=1}^{4} y_{ij}$ 式中,参照树i与相邻木j的树冠投影重叠时,yij=1;否则,yij=0;并包含5种取值与状态:0表示很稀疏;0.25表示稀疏;0.5表示中等密集;0.75表示较密集;1表示很密集,数值越大,说明树冠连接越密,树冠相互重叠挤压的程度也越高。 林分 $\bar{C}=\dfrac{1}{n} \sum\limits_{i=1}^n C_i $ 式中,Ci表示林木个体的密集程度,n为全林分株数。 竞争强度 Competition 林木个体 $C I_{t}=\sum\limits_{f=1}^{m} \dfrac{d_{j} }{d_{i} \cdot L_{ij} }$ 式中,CIi为林木i的竞争强度,Lij为林木i到相邻木j的距离,di为对象木i的胸径,dj为相邻木j的胸径,n为竞争木株数。 林分 $ C I=\dfrac{1}{N} \sum\limits_{i=1}^{N} C I_{i} $ 表示林分竞争强度,CIi为林木个体i的竞争强度,N为林分内林木总株数。 Table 2. Formula of crowding and competition
-
利用Office软件进行数据预处理,分别利用森林结构分析专业软件Winklemass与R studio计算密集度与竞争指数。采用最小显著性差异法(LSD)对树冠结构特征进行多重比较,采用Pearson相关系数评价树冠结构特征与林木胸径大小的相关性。采用结构方程模型(SEM)分析林分密度、树冠结构和密集度的关系,显著性水平均采用P=0.05,其中,依据拟合指数,包括不显著(p>0.05)卡方检验(χ2),SRMR<0.08,GFI>0.95,CHI>0.95,以及最小的AIC值选取最优的结构方程模型[17]。数据分析主要在R studio中完成,利用Origin9.4绘图。为避免边缘效应造成的误差,在样地周围设1.5 m的缓冲区域。
-
蒙古栎林中,冠幅与树冠面积随着林分密度的递增而递增,冠长与树冠体积的变化趋势相似,并随林分密度的变化呈单峰型,虽然冠长率与冠形率也呈单峰型,但其变化趋势相对平缓,说明不同密度林分的树冠结构特征存在差异(图1a)。分析样地间、样地内树冠结构的变异(图1b),结果表明树冠体积的变异最强,变异程度介于28%~45%;其次为树冠面积与冠长,分别为18%~28%与13%~20%;而冠幅、冠长和冠长率的变异相对较弱,分别介于9%~14%、7%~9%和13%~16%。总体上,蒙古栎林树冠结构的变异未与林分密度保持一致性,而是随着林分密度增大呈现出“先增后减”的单峰型特点。
-
如图2所示,样地Ⅰ中林木胸径与树冠结构特征的相关性均显著,胸径与树冠体积的相关性最强(r=0.823,P<0.01),而与冠形率(r=0.188,P<0.05)和冠长率(r=0.216,P<0.05)的相关性较弱。样地Ⅱ中树冠面积(r=0.68, P<0.01)、冠幅(r=0.67, P<0.01)分别与胸径的相关性最强,其次为树冠体积(r=0.65,P<0.01)。与样地Ⅱ相似,样地Ⅲ中冠幅、树冠面积分别与胸径的相关性最强(P<0.01),其次为树冠体积(P<0.01),但2个样地的冠形率、冠长率与胸径无显著相关性。样地Ⅳ中树冠面积、冠幅、树冠体积分别与胸径的相关性最强(P<0.01),而冠形率、冠长率与胸径的相关性较弱(P<0.05)。总体上,随林分密度的递减,冠幅、树冠面积、树冠体积分别与胸径的相关性趋势近似于“L”型,而冠长—胸径的相关性随林分密度的递减逐渐减弱,林分密度的改变对冠长率—胸径、冠形率—胸径相关性的影响不大。
-
蒙古栎林中,4.5%的林木处于很稀疏状态,15.9%为稀疏,28.8%为中等密集,20.9%为较密集,29.7%为很密集,林分密集程度处于中等密集至很密集状态。具体而言(图3),样地Ⅰ与样地Ⅲ中最大相对频率分别为45.2%、27.8%的林木属于很密集,而样地Ⅱ和Ⅳ中最大相对频率分别为32.7%、40.6%的林木属于较密集与中等密集。分析优势树种的密集程度可知(表3),其最大值为0.909,最小值为0.292,包括稀疏、中等密集、较密集和很密集4种状态,但多数优势树种处于中等密集、较密集、很密集状态。就各样地优势树种的密集程度而言,样地Ⅰ中蒙古栎、华山松和椴树的值最大,且至少54.1%的林木达到很密集,尽管处于很密集的白桦较少,但其中也有46.2%达到中等密集。样地Ⅱ中34.1%的蒙古栎达到较密集,40%的白桦很密集,而华山松与椴树分别均匀分布于中等密集至很密集状态。样地Ⅲ中35%的白桦与57.1%的椴树达到了很密集,而此时多数蒙古栎和华山松处于稀疏与较密集状态。样地Ⅳ的优势树处于很稀疏至中等密集状态,其中蒙古栎31.2%、华山松50%、白桦41.7%以及椴树50%分别为稀疏、很稀疏、稀疏和中等密集,说明优势树种对蒙古栎林密集程度起重要作用。
样地
Sample plot优势树
Dominant trees相对频率 Relative frequency 密集度值
Value很稀疏
Thinly scattered稀疏
Scattered中等密集
Medium较密集
Intensive很密集
Extremely intensiveⅠ 蒙古栎 Quercus mongolica 0.000 0.059 0.235 0.059 0.647 0.824 华山松 Pinus armandii 0.000 0.000 0.090 0.183 0.727 0.909 白桦 Betula platyphylla 0.000 0.307 0.462 0.076 0.155 0.520 椴树 Tilia tuan 0.000 0.043 0.333 0.083 0.541 0.781 Ⅱ 蒙古栎 Quercus mongolica 0.073 0.173 0.292 0.341 0.121 0.566 华山松 Pinus armandii 0.000 0.250 0.250 0.250 0.250 0.625 白桦 Betula platyphylla 0.000 0.000 0.266 0.333 0.400 0.783 椴树 Tilia tuan 0.000 0.000 0.333 0.333 0.333 0.749 Ⅲ 蒙古栎 Quercus mongolica 0.055 0.368 0.236 0.263 0.078 0.538 华山松 Pinus armandii 0.000 0.250 0.125 0.500 0.125 0.500 白桦 Betula platyphylla 0.000 0.200 0.250 0.200 0.350 0.675 椴树 Tilia tuan 0.000 0.000 0.142 0.285 0.571 0.713 Ⅳ 蒙古栎 Quercus mongolica 0.000 0.312 0.251 0.187 0.250 0.594 华山松 Pinus armandii 0.500 0.000 0.333 0.167 0.000 0.292 白桦 Betula platyphylla 0.083 0.417 0.333 0.084 0.083 0.417 椴树 Tilia tuan 0.167 0.167 0.500 0.083 0.083 0.437 Table 3. Crowding of dominant tree species in different densities
-
最优结构方程模型表明(图4,表4),林分密度对树冠结构与密集度关系的影响包含两种路径:一方面,林分密度通过林木竞争间接影响了树冠结构(p<0.05),通过树冠结构的介导作用从而间接影响密集度(p<0.05);另一方面,林分密度通过林木竞争先间接影响林木胸径(p<0.001),再通过林木胸径的介导作用影响树冠结构,进而对密集度产生影响(p<0.05);其中,就树冠结构特征而言,冠幅与冠长的显著作用不可忽视(p<0.001)。尽管如此,但林分密度直接通过林木竞争也可间接影响密集度(p<0.001)。另外,林分密度通过影响林木竞争与胸径大小进而对密集度产生间接影响,但此影响并不显著(p>0.05)。进一步判定蒙古栎林的林分密度、林木胸径与树冠结构对密集度的联合影响可知,林分密度能解释林木竞争47.1%的变异,林木竞争解释了胸径0.1%的变异,林木胸径与竞争共同解释了树冠结构62.8%的变异,林木竞争、胸径大小与树冠结构共同解释了密集度23.6%的变异。
响应变量
Response variable解释变量
Explanatory variable影响
EffectBeta P 值
P Value林木竞争 Tree competition 林分密度 Stand density 直接 Direct −0.686 <0.001 胸径 DBH 林木竞争 Tree competition 直接 Direct 0.024 0.612 树冠结构 Crown structure 林木竞争 Tree competition 直接 Direct 0.098 <0.05 密集度 Crowding 林木竞争 Tree competition 直接 Direct 0.416 <0.001 树冠结构 Crown structure 胸径 DBH 直接 Direct 0.784 <0.001 密集度 Crowding 胸径 DBH 直接 Direct −0.019 0.827 密集度 Crowding 树冠结构 Crown structure 直接 Direct 0.219 <0.05 胸径 DBH 林分密度 Stand density 间接 Indirect −0.017 — 树冠结构 Crown structure 林分密度 Stand density 间接 Indirect −0.08 — 密集度 Crowding 林分密度 Stand density 间接 Indirect −0.303 — Table 4. Direct and indirect effect of stand density, tree competition, crown structure on crowding
Relationships between Crown Structure and Crowding in Quercus mongolica Forest of Xiaolong Mountains under Different Stand Densities
- Received Date: 2023-04-06
- Accepted Date: 2023-06-05
- Available Online: 2024-02-20
Abstract: