-
在全球气候变化背景下,中国政府在第七十五届联合国大会上提出二氧化碳排放力争于2030年前达到峰值、2060年前实现碳中和的重大战略目标,人工林碳汇是其中的重要组成部分[1]。影响林分年内生长的因素复杂,包括温度、降水等气候因子以及林分密度、竞争等结构特征。林分密度是最重要的林分结构特征之一[2],与径向生长有直接关系,如Jordan Bello等[3]发现,低密度林中无柄橡树生长高于中密度林。在其他许多关于温带针叶林径向生长影响因子的研究中,太阳辐射、降水、饱和水汽压差和日最低气温等因子则表现出不同规律[4-9]。
华北落叶松(Larix principis-rupprechtii Mayr.)作为宁夏六盘山地区的主要造林树种,发挥着水土保持、水源涵养、固碳释氧等多种生态功能。目前关于华北落叶松林分的径向生长的年内变化特征,特别是对不同密度人工林如何响应环境因子变化,尽管已有较为丰富的研究成果[10-12],但受限于不同观测方法和分析手段及研究时段,当前研究结果中对同一因子影响的正负性和主导因子判定存在不同结论,如熊伟[10]等发现日最低温度与径向生长量有显著正相关关系,王亚蕊等[13]的研究则表明,该变量与生长量呈不显著负相关;管伟等[11]研究结果表明,空气相对湿度与生长存在正相关,Liu[12]的研究结果则表现为不显著的负相关关系。本研究通过线性混合效应模型分析了各环境因子对华北落叶松径向生长的影响,以期从另一角度探讨环境因子的驱动作用,并与相关分析结果对比验证。本研究拟解决的主要问题为:(1)分析不同密度下华北落叶松人工林年内径生长的动态变化;(2)确定影响其生长的主要环境因子,并量化各因素影响华北落叶松人工林径向生长的相对贡献。为此,本研究利用树木生长尺监测了六盘山华北落叶松树木个体的径向生长动态,并分析其树干径向生长的季节变化格局,以明确不同密度下树木年内径向生长的主要影响因素,为准确评估和预测未来华北落叶松人工林树木径向生长和森林动态、实现可持续管理与多功能维持提供科学依据。
-
研究区位于宁夏六盘山南坡的香水河小流域(106°12′~106°16′ E,35°27′~35°33′ N),海拔2 070~2 931 m,小流域总面积43.7 km²。研究区属于温带半湿润气候,雨热不同期,年均气温6.8 ℃,年均降水量710.3 mm,降水主要集中在6—9月份,全年无霜期90~130 d。土壤类型为森林灰褐土,土壤质地为壤土,土壤厚度为80 cm。石砾含量为8%~18%。研究区内乔木以华北落叶松人工林为主,天然次生林有白桦(Betula platyphlla Suk.)、华山松(Pinus armandii Franch.)等。灌木林主要有虎榛子(Ostyopsis davidiana Decne.)、灰栒子(Cotongaster acutifolius Turcz.)等。林下草本主要有苔草(Cyperaceae carex L.)、东方草莓(Fragaria orientali Lozinsk.)等。
-
在香水河小流域内选择4块不同林分密度的华北落叶松人工纯林样地,样地面积40 m × 40 m。为消除小地形影响,样地坡向相同(东南坡向)、坡度接近、坡位相同(均为下坡位)且海拔高度基本一致的相邻样地(表1)。对于样地间坡度差异,取林分密度接近、坡度有差异的2个高密度样地(样地3、4)进行了方差分析,结果表明其平均累积径向生长增量差异不显著,本研究样地间坡度差异不影响生长量分析。
样地编号
Sample NO.海拔
Elevation/m坡度
Slope/(°)林龄
Age/a平均胸径
DBH/cm平均树高
Height/m优势木平均高
Average height of
dominant trees/m样本量
Sample
size林分密度
Stand density/
(株·hm−2)1 2 336 18.0 39 21.4 ± 4.6 17.6 ± 2.3 20.9 ± 1.8 76 544 2 2 351 28.8 39 20.8 ± 3.0 19.3 ± 1.5 21.2 ± 1.1 91 742 3 2 300 29.0 39 19.0 ± 3.9 17.7 ± 2.4 20.5 ± 1.5 151 999 4 2 346 16.2 39 19.7 ± 4.0 19.5 ± 2.9 22.4 ± 1.5 144 1 020 Table 1. General characteristics of the sample plot
-
将4个样地内的462株华北落叶松均选为样树,在4个样地内每棵树干距地1.3 m处刮去干枯树皮,安装QT-W304树木生长尺。观测时间从2022年5月开始,记录开始生长时的初始刻度值,每隔5~7 d观测生长尺刻度变化,观测至10月底生长季结束。最终获取到438株样木的生长尺有效数据。
-
在距离研究样地100 m处的空旷草地内设有标准气象观测场,使用 Weatherhawk232自动气象站监测气象因子,监测指标包括空气温度、降水量、空气相对湿度(以下简称空气湿度)、太阳辐射、风速(m·s−1),每隔1 min采集1次数据,每隔5 min记录1次平均值。利用空气温度和相对湿度计算饱和水汽压差[14]。
在样地内按 0~10、10~20、20~40 、40~60 和60~80 cm土层深度布设 EC-5 土壤水分传感器(Decagon,美国),并利用 EM 50 数据采集器(Decagon,美国)每 5 min 记录 1 次数据。本研究环境因子分析提到的土壤含水量均为土壤体积含水量,为0~80 cm土层实测土壤体积含水量的平均值(图1)。
-
Gompertz 函数在树木径向生长拟合中有广泛的应用[15],可以用来确定树木的生长季。使用Gompertz函数对不同密度华北落叶松生长季内平均累积径向生长量随时间变化进行拟合,然后对拟合函数求导得到华北落叶松生长速率。Gompertz函数模拟公式如下式[16]:
式中,Y为树干累积径向生长量;A为上渐近线;β表示与Y初始值有关的参数;k表示决定生长速率的参数;t表示时间;将拟合方程求导,得到华北落叶松的径向生长速率。
采用树径变化速率法判定快速增大阶段[13,17-18],即认为树干直径变化速率开始大于4 µm·d−1时为快速生长期的开始时间,开始小于4 µm·d−1时为结束时间[17]。
利用线性混合效应模型分析2022年华北落叶松日径向生长量与环境因子之间的关系,将胸径(表征树木个体间的竞争能力)、土壤水分、平均气温、饱和蒸汽压差、降水量、太阳辐射和风速作为固定效应,将林分密度(表征林分尺度的树木竞争)、生长时期(以月份表征,作为生长阶段的次一级分期)作为随机效应。用R语言中“lmerTest”包的“lmer”函数进行线性混合效应模型的分析。利用glmm.hp包[19]分析影响树木径向生长各因素的相对重要性,使用Origin 2018进行绘图。
-
图2为生长季节内华北落叶松的径向变化图,华北落叶松径向生长的季节变化呈现“S”型。2022年DOY 152—163(6月1日—6月13日)起,华北落叶松径向生长量开始大幅度增加;参考王亚蕊等[13]对六盘山华北落叶松研究结果,可将本研究生长季内树木直径变化过程分为3个阶段:1)相对稳定阶段(5月及9月上旬以后,DOY138—150及DOY245以后),此时树木直径净增量基本为零;2)快速增大阶段(6月初—7月下旬,DOY152—210),此时树木直径快速增大,但因树木蒸腾失水和土壤体积含水量变化等影响而呈波动增长趋势;3)涨缩波动阶段(8月上旬—9月上旬,DOY212—245),此时树木直径虽波动增长,但生长较慢(图2)。
Figure 2. Seasonal pattern of stem diameter variation of Larix principis-rupprechtii at different stand densities
图3为2022年测量和模拟的不同林分密度下华北落叶松平均累积径向生长量和日平均径向生长速率变化图,由图3a可知,Gompertz函数模拟结果显示,华北落叶松平均累积径向生长量拟合曲线的决定系数(R2)均≥0.94,模型模拟与实测数据趋势基本一致。华北落叶松年内径向生长呈“S”型。不同林分密度下,华北落叶松的生长开始于5月—6月初(DOY121—163),7月中下旬(DOY198—209)达到高峰,9月上旬左右(DOY245)停止生长。
Figure 3. Gompertz function-modeled curves of the cumulative stem growth(a) and daily growth rate(b) at different stand densities
图3b为模拟得到的华北落叶松日平均径向生长速率,本研究中分析华北落叶松树干径向生长的快速增大阶段为DOY 152—212(6月至7月),即生长速率>4 µm·d−1的阶段。4个林分密度的华北落叶松快速生长期分别为58,57,44和36 d,表现出林分密度越小快速增大阶段越长的趋势,随着林分密度的增大,快速增大阶段持续时间缩短。从生长速率来看,2022年华北落叶松径向生长速率曲线表现为单峰型;4个林分密度下最大生长速率依次为32.6、23.6、17.6 和22.5 µm·d−1,林分密度越低,最大生长速率越高。最大生长速率出现在DOY 171—174(6月20日—6月23日)左右。
不同林分密度下生长的季节变化趋势不尽相同(图3b)。随着林分密度增大,华北落叶松快速增大阶段开始的时间较晚,结束时间较早,导致径向生长快速增大阶段时间少于低林分密度树木。但快速生长期开始时间均为6月初,快速生长期结束时间为7月末。就整个生长季而言,华北落叶松5月中旬开始生长,但生长速率较低,8月末生长结束,9月停止生长。由图4可以看出,2022年各密度下华北落叶松径向生长主要发生在阶段B快速增大阶段(6月初—7月下旬,DOY152—210),在544、742、999和1 020 株·hm−2的密度下,华北落叶松快速增大阶段相对径向生长量分别为71.5%、67.5%、66.7%和64.1%,除1 020 株·hm−2密度外,不同密度华北落叶松径向生长在快速增大阶段(阶段B)均超过整个生长季的2/3;相对稳定阶段和涨缩波动阶段华北落叶松径向生长均较低。
-
表2为混合效应模型分析的气象、树木大小、林分密度和生长时期对华北落叶松日径向生长的影响结果。由表2可见,树木胸径、日平均气温、土壤含水量、降水量、饱和水汽压差、太阳辐射和风速均对华北落叶松日径向生长量有显著影响。作为表征树木之间竞争强弱的因子,树木胸径与华北落叶松的径向生长呈正相关性,即胸径越大,树木生长越快。对气候因素来说,树木日径向生长随平均气温、降水量和风速增加而显著增加,但随土壤体积含水量、饱和水汽压差和太阳辐射增加而显著下降,这可能与影响华北落叶松树木生长直接相关的冠层蒸腾作用有关。对比不同因素在决定华北落叶松日径向生长中的作用发现,土壤水分和平均气温是影响日径向生长的最主要因子,其相对贡献率分别达到30.30%和27.24%。太阳辐射、风速和饱和蒸气压差作为影响树木蒸腾作用的主要因素,它们的相对贡献率合计为38.03%,个体竞争因素(胸径大小)对树木日径向生长的影响作用非常小。然而,林分密度(表征林分尺度的竞争)和生长时期(即表征其生长时期)也显著影响华北落叶松的日径向生长。模型模拟结果表明,同时考虑林分密度、生长季月份、树木大小以及气象因子的模型可解释72.0%的方差,但排除林分密度和生长时期因素后,树木大小和气象因子共同仅能解释46.4%的方差。
模型变量
Predictors估计值
Estimate95%置信区间95%
Confidence intervals标准误
SEp值
p value偏相关系数
Partial R2相对贡献率
Relative importance/%固定效应 Fixed effects 截距 Intercept 0.06 0.05~0.06 0.003 7 <0.001 胸径 DBH 0.00 0.00~0.00 0.000 0 <0.001 0.004 5 0.97 土壤含水量 Soil content −0.20 −0.22~−0.19 0.006 6 <0.001 0.140 7 30.30 平均气温 Mean temperature 0.00 0.00~0.00 0.000 0 <0.001 0.126 5 27.24 饱和水汽压差 VPD −0.01 −0.01~−0.01 0.000 9 <0.001 0.027 9 6.01 降水量 Rainfall 0.00 0.00~0.01 0.000 6 <0.001 0.016 1 3.47 太阳辐射 Solar radiation −0.00 −0.00~−0.00 0.000 0 <0.001 0.063 3 13.63 风速 Wind speed 0.02 0.01~0.02 0.000 6 <0.001 0.085 4 18.39 随机效应 Random effects σ2 0.00 τ00 月 0.00 τ00 林分密度 0.00 ICC 0.48 N林分密度 4 N月 6 观测数量 Observations 6 568 边际R2/条件R2
Marginal R2/Conditional R20.464/0.720 注:σ2为残差的方差,τ为随机效应引起的方差,ICC为类内相关系数,N是月或林分密度分类的数量,边际R2是固定效应的变异系数,条件R2是固定效应和随机效应影响的变异系数。
Note: σ2 is the variance of residuals, τ is the variance caused by random effects, ICC is intra-class correlation coefficient, N is the number of plots, Marginal R2 is the coefficient of variation for fixed effects, and Conditional R2 is the coefficient of variation for fixed and random effects.Table 2. Effects of daily meteorological variables, tree size, stand density and month in growth season on daily radial increment of Larix principis-rupprechtii plantation during growing season
Intra-annual Radial Growth Dynamics of Larix principis-rupprechtii Plantations with Different Stand Densities and Its Response to Environmental Factors
- Received Date: 2023-05-06
- Accepted Date: 2023-06-26
- Available Online: 2024-04-27
Abstract: