间伐与气候对长白落叶松树轮宽度的影响*

崔诗梦 向 玮

(北京林业大学林学院 北京 100083)

全球气候变化对森林生态系统动态有着重要影响(Henttonenetal., 2014; Guillemotetal., 2015)。森林生长对气候变化如何响应，而人类对此应采取怎样的应对措施受到广泛关注。树木生长受树种、年龄等内在因素(Gerendiainetal., 2012)及外部环境条件如林分密度、立地因子和气候因子等影响。雷相东等(2009)发现林分密度、立地因子和单木竞争对树木生长均有显著影响。气候因子对树木生长的影响较为复杂，不仅生长季内的温度和降水等气候因子会影响树木生长(Wimmeretal., 1997; 李宗善等, 2010; Chenetal., 2012)，生长季以前的气候因子也会对树木生长产生影响(Lebourgeois, 2000; Linderholm, 2001; Mäkinenetal., 2003)。这种现象与在生长季之前的气候因子对土壤水分、树木的生理活动等方面的影响有关(Rolland, 1993; Jyskeetal., 2012)。树木生长与气候因子的响应关系会受树种(Franketal., 2005; 王丽丽等, 2005)、年龄(Szeiczetal., 1994; Carrer, 2004)、地理位置(Savvaetal., 2006; Henttonenetal., 2014)等因素的影响。为保证树木生长对气候因子变化的敏感性，研究地多选在寒冷、干旱地区的林线处(kirdyanovetal., 2003; 于健等, 2016)。在林分内部树木生长的研究中，通常认为在一个区域内气候是较为稳定的，因而不考虑气候因子对树木生长的影响。然而在全球气候变化的背景下，气候的稳定性正在发生改变，有必要考虑到气候因子对森林内部树木生长的影响。

目前，人类活动干扰也已成为影响树木生长的因素。间伐是常用的营林措施，可以降低林内竞争(Millaretal., 2007; 贾忠奎等, 2012a; Magruderetal., 2013)，改善林内土壤、光照等条件(Barbouretal., 1994; 李国雷等, 2008; 贾芳等, 2009;贾忠奎等, 2012b)，提高保留木对气候因子的适应能力(Missonetal., 2003; Kerhoulasetal., 2013; Lebourgeoisetal., 2014; Oliveretal., 2014)，促进保留木生长(Kogaetal., 1997; 雷相东等, 2005; Guller, 2007)。不同强度的间伐对保留木的促进作用也有所差异，通常，中、高强度间伐对树木生长的促进作用较强，而低强度间伐对树轮宽度的影响较小或不显著(马履一等, 2007; 莫日根等, 2013; 龚固堂等, 2015)。间伐对树木生长的影响会在伐后数年间逐渐增强后减弱(Missonetal., 2002; Kogaetal., 2002; Mehtataloetal., 2014)。研究不同间伐强度对林木生长的影响，有助于为特定林分找出合适的经营方案。

树木年轮记录了树木逐年的径向生长，具有分辨率高和易于获取等特点(吴祥定等, 1996)。利用年轮宽度数据研究树木生长时，由于同一样地或同一个体间存在相关性，因此不能满足独立同分布的假设。混合模型可通过规定不同的协方差结构来表示相关的误差，允许数据间存在相关性和方差异质性，可提高预测精度和解释随机误差来源，因此被广泛应用于树木生长的研究中(雷相东等, 2009; 李春明等, 2010)。长白落叶松(Larixolgensis)以其生长快、适应性强的特点，作为先锋树种，大量用于我国东北地区人工林的营造，其生长与收获受到广泛关注。本研究以长白落叶松为对象，建立树轮宽度指标的非线性混合模型，分析间伐与气候因子对树木生长的影响，为全球气候变化条件下的林业生产等营林活动、森林生态系统管理提供参考依据。

1 研究区概况

研究地位于吉林省汪清林业局金沟岭林场(130°05′—130°20′E，43°17′—43°25′N)，属长白山系老爷岭山脉雪岭支脉。低山丘陵，海拔300～1 200 m，平均坡度10°～25°。温带大陆性季风气候，年均气温3.9 ℃，1月份气温最低，平均最低气温-32 ℃，7月份气温最高，平均最高气温32 ℃； 全年无霜期138天； 年均降水量600～700 mm，多集中在7月。土壤种类以暗棕壤为主，平均厚度40 cm左右。该地区植被属长白山植物区系。

2 研究方法

2.1 样地设置和调查

本研究样地数据采自金沟岭林场的20块固定样地。其起源为1964—1967年营造的长白落叶松人工林，经多年演变，大部分成为长白落叶松云冷杉针阔混交林。林内除长白落叶松外，还有鱼鳞云杉(Piceajezoensis)、臭冷杉(Abiesnephrolepis)、红松(Pinuskoraiensis)、色木槭(Acermono)、白桦(Betulaplatyphylla)、水曲柳(Fraxinusmandschurica)、椴树(Tiliatuan)和榆树(Ulmuspumila)等。20块样地按照区组试验设计，分为5个区，每区包含的4块样地在1987年分别按胸高断面积的百分比进行轻度(20%)、中度(30%)、重度(40%)下层疏伐及对照处理，分别标记为A1, A2, A3, A4, B1, B2, B3, B4, C1, C2, C3, C4, D1, D2, D3, D4, D5, E1, E2, E3, E4； 其中D1和D2样地于1993年再次分别进行了强度为30%和20%的间伐。各样地概况与间伐情况见表1。样地在1987—2012年间进行过13次调查，间隔期为2～3年，内容包括树种、胸径和立地因子等。对于其间未调查年份的数据通过线性插值补全(Filipescuetal., 2014)。

2.2 数据来源

年轮宽度数据来源为2016年夏季采集的长白落叶松生长锥年轮条。年轮条采样前，依据样地调查数据中的胸径转化为径阶将林木划分为3个等级： 1级为样地内径阶不低于四分位点的林木，2级为径阶不足1级木但高于林分平均胸径所在径阶的林木，3级为低于林分平均胸径所在径阶的林木。通过简单随机抽样的方式从各样地内1，2，3级林木中分别抽取4株长白落叶松钻取年轮条。共采集243株长白落叶松林木的年轮条，其中，4条样品(1条采自2级木，3条采自3级木)断裂难以恢复。对剩余的年轮条通过Lintab 6树轮宽度测量仪测量年轮宽度，并将年轮宽度数据通过COFECHA程序进行交叉定年。交叉定年结果显示8条样品年轮宽度数据与主序列相关性不显著，确认测量无误后舍弃。最终231条样品年轮宽度数据通过交叉定年。保留的各样地年轮条样品信息见表1。年轮宽度数据的80%用于建模，20%用于检验。

气候数据为样地附近气象站提供数据的线性插值结果(Shenetal., 2015)，包括1964—2010年间每月最高气温、最低气温、平均气温、降水量及生长季大于5 ℃积温等。

本研究使用的1987—2010年的样地、气候与树轮宽度数据。

2.3 模型建立

(1)

对于树轮宽度随形成层年龄变化的趋势，选用三参数Logistic方程模拟，形式如下：

f(α,A)=α0+α1/[1+e(α2×A)]。

(2)

式中：α0表示生长稳定时的树轮宽度指标；α1为靠近髓心的树轮宽度指标；α2为尺度参数。

间伐对树轮宽度的影响通常表现为在伐后数年间先增大后逐渐减弱的趋势(Missonetal., 2002; Mehtataloetal., 2014)，因此在比较不同方程形式后选择公式(3)拟合树轮宽度指标受间伐影响并随伐后时间变化的趋势。不同间伐强度对树轮宽度指标的影响以哑变量的形式表现。

g(β,T)=(β111T11+β112T12+β113T13)t1/

(3)

式中：T11为第1次轻度间伐的哑变量，T12为第1次中度间伐的哑变量，T13为第1次重度间伐的哑变量，T21为第2次轻度间伐的哑变量，T22为第2次中度间伐的哑变量，所有哑变量均是为1，否为0。β111为第1次轻度间伐对树轮宽度指标影响程度的固定效应参数；β112为第1次中度间伐对树轮宽度指标影响程度的固定效应参数；β113为第1次重度间伐对树轮宽度指标影响程度的固定效应参数；β121为第2次轻度间伐对树轮宽度指标影响程度的固定效应参数；β122为第2次中度间伐对树轮宽度指标影响程度的固定效应参数；β21为第1次间伐对树轮宽度指标影响程度与持续时间的固定效应参数；β22为第2次间伐对树轮宽度指标影响程度与持续时间的固定效应参数；t1为第1次间伐后年数；t2为第2次间伐后年数。

2.4 模型自变量选择与参数估计

影响树木生长的因子较多，为研究方便，将除形成层年龄与间伐因素以外的自变量分为4组： 立地因子组、林分因子组、单木竞争组和气候因子组。各变量组包含变量及说明见表2。考虑到应用，单木竞争组变量采用较为方便获取的与距离无关的指数。林分密度指数(SDI)通过公式(4)计算。

SDI=Ns(Dg/D0)a，

(4)

式中：Ns为林分每公顷株树；Dg为林分平均胸径；D0为标准胸径；a为自然稀疏率； 本研究中长白落叶松标准直径与自然稀疏率的取值与李春明等(2010)的研究一致：D0=20 cm，a=1.58。

本研究自变量较多，且气候因子间存在相互作用，各自变量间可能存在共线性，影响参数估计的精确性。研究中利用方差膨胀因子(VIF)判断自变量间的多重共线性。仅保留VIF<5且回归系数显著的自变量(雷相东等, 2009)。本研究中各自变量量纲不一致，为通过模型参数比较各自变量对因变量影响的大小，将表2中自变量标准化，此时模型参数代表在自变量变化1个标准差时因变量的变化情况。

对于可能存在的误差的方差异质性和序列相关性，将优选出合适的方差函数(幂函数、常数幂函数)及协方差结构(自回归、滑动平均、自回归滑动平均)。模型的选择基于Akaike’s information criterion (AIC值)(Akaike, 1974)。对于嵌套模型，AIC能够辨识出较好的模型，其值小者为优。模型参数估计采用极大似然法。所有计算通过R软件的nlme包完成。

①生长季为4—9月Growth season refers to the period between April and September; 夏季为6—8月Summer refers to the period between June and August; 冬季为前一年12月至该年3月Winter refers to the period from December in the previous year to March of that year.

2.5 模型评价及检验

通过以下指标对模型拟合结果进行评价： 决定系数R2、绝对偏差(Bias)、相对偏差(Bias%)、均方根误差(RMSE)和相对均方根误差(RMSE%)，见公式(5)-(9)(雷相东等, 2009)。

；

(5)

(6)

(7)

(8)

(9)

3 结果与分析

3.1 模型结果

除形成层年龄和间伐因素外，将9个自变量经标准化后引入模型，9个自变量分别为林分密度指数(SDI)、对象木胸径与林分平均胸径的比值(Dr)、生长季大于5 ℃积温(Dgs)、夏季热湿指数(Hs)、冬季降水量(Pw)、7月最高气温(T7max)、9月最低气温(T9 min)、5月降水量(P5)和4月平均温(T4mean)。最终模型见式(10)。模型的参数估计结果见表3。除表示第1次强度为轻度的间伐造成影响的参数β111不显著外，其他参数均在95%置信水平上显著(P<0.05)。模型中加入样地与单木水平随机效应后，AIC值降低，因此将样地与单木水平的随机效应引入模型。协方差结构选择一阶自回归结构AR(1)，此时模型AIC值最小。

；

其中：

f(α,A)=α0+α1/[1+e(α2×A)]，

g(β,T)=(β111T11+β112T12+β113T13)t1/

γX=γ1SDI+γ2Dr+γ3Dgs+γ4Hs+

γ5Pw+γ6T7max+γ7T9min+γ8P5+γ9T4mean，

εijk～(0,σ2Rij)，

σ2Rij=σ2Γij(φ)，

Γij(φ)=AR(1)。

(10)

式中：Dr为对象木胸径与林分平均胸径的比值；Dgs为生长季大于5 ℃积温；Hs为夏季热湿指数；Pw为冬季降水量；T7max为7月最高气温；T9 min为9月最低气温；P5为5月降水量；T4mean为4月平均气温；γ1,γ2,γ3,γ4,γ5,γ6,γ7,γ8和γ9均为待估参数；Γij(φ)为序列相关函数；AR(1)为一阶自回归结构。

表3 模型参数估计、标准误差和P值①Tab.3 Parameter estimates, standard errors, and P-values for the model

①SD(bi)为bi的标准差 SD(bi) is the standard deviation ofbi； SD(bij)为bij的标准差 SD(bij) is the standard deviation ofbij.

图1为模型参数估计结果代入公式(1)所绘制曲线，显示形成层年龄在15～45年的时间段内，树轮宽度指标逐渐减小并趋于稳定。

图1 树轮宽度指标随形成层年龄变化趋势Fig.1 Trend of tree-ring width index variation according to cambial age

项目ItemBiasBias%RMSERMSE%R2固定效应R2R2offixedeffect建模结果Modelingdata00003-05592006131669660730103680检验结果Testingdata00015-01294006131684180728703294

在间伐效应方面，模型中第1次轻度间伐的参数没达到95%置信水平显著，表明该次间伐对树轮宽度的影响并不显著，其余间伐措施均对长白落叶松树轮宽度指标表现出显著影响。图2为将模型参数估计结果代入公式(2)绘制的曲线，图中显示树轮宽度指标在不同间伐时间与间伐强度间有所差异： 在第1次间伐中，重度间伐引起树轮宽度指标增量最多，中度间伐次之，轻度间伐最少，在伐后1～2年达到峰值，分别为0.038，0.029和0.014个单位，后逐渐减弱并持续至约伐后4年； 第2次间伐中，轻度间伐引起的树轮宽度指标增量比中度间伐大，且在伐后4～5年达到峰值，分别为0.069和0.057个单位，后逐渐减弱约持续至伐后15年。

图2 不同间伐时间及间伐强度对树轮宽度指标影响Fig.2 Effects of different thinning time and tinning intensity on tree-ring width index

描述立地因子的备选变量均未对长白落叶松树轮宽度指标表现出显著影响； 林分因子组中SDI与单木竞争组中Dr对长白落叶松树轮宽度指标表现出显著影响。气候因子组中有7个变量对长白落叶松树轮宽度指标表现出显著影响。其中生长季大于5 ℃ 积温和7月最高气温每升高1个标准差所引起树轮宽度指标的变化量最多，分别为-0.023 6和0.022 0个单位； 夏季热湿指数的影响次之，每升高1个标准差引起树轮宽度指标降低0.010 0个单位； 4月平均温、5月降水量、9月最低气温和冬季降水量对树轮宽度指标的影响相对较小。表4给出了模型中的自变量标准化计算时所用到的均值和标准差。

表4 模型自变量均值和标准差Tab.4 Mean value and standard deviation of variables

3.2 模型评价与检验

将检验数据代入模型，计算出评价统计量，与建模数据拟合结果的评价统计量进行对比(表5)，显示模型检验与拟合的各评价统计量基本一致。引入随机效应参数后，模型的决定系数(R2)有大幅提高，对于建模和检验数据而言，分别从0.37和0.33提高至0.73和0.73，表明样地及树木个体间的差异是随机误差的重要来源。绝对偏差(Bias)、相对偏差(Bias%)、均方根误差(RMSE)和相对均方根误差(RMSE%)在建模数据和检验数据的拟合结果中相差不大，认为模型稳定可靠。

图3显示模型残差分布均匀，未见明显趋势，因此不考虑方差异质性。

图3 模型残差与拟合值关系Fig.3 Standardized residuals versus fitted values of model

4 讨论

间伐对长白落叶松的生长表现出先增强后减弱的促进作用(图2)，这样的趋势也见于Misson等(2002)、 Koga等(2002)和Mehtatalo等(2014)的研究中。这种促进效果受间伐强度影响，并且在2次间伐中表现出差异。第1次间伐时林龄为20年，间伐后树木生长加快，重度间伐下林木径向生长最快，中度间伐次之，轻度间伐与对照组差异不显著，与其他学者的研究结果一致(马履一等, 2007; 莫日根等, 2013; 龚固堂等, 2015)。第2次间伐中轻度间伐比中度间伐对树木生长表现出了更强的促进作用，而其他学者关于间伐的研究中通常为轻度间伐对树木生长的促进效果弱于中度间伐(马履一等, 2007; 莫日根等, 2013; 龚固堂等, 2015)，原因可能是本研究中2次间伐相隔时间较短，仅有6年，树木生长还没有受到林内资源的强烈限制。在这样的情况下，30%强度的间伐与20%强度相比，对林内环境的破坏更多，从而抵消了部分间伐对林木生长带来的好处，因此表现出中强度间伐比低强度间伐对树木生长的促进作用弱。第2次间伐对长白落叶松生长的影响比第1次间伐更强且持续时间更长。原因可能是第1次间伐时林龄较小，树木生长旺盛，间伐对林内竞争的缓解由于林木的快速生长而迅速减弱。随着林龄增长，第2次间伐时林木的生长速度减缓(图1)，因此该次间伐对林木竞争的缓解可以持续更长时间。然而，本研究中仅2块样地进行了第2次间伐，样本量小，有偶然因素造成这些现象的可能。

本研究中对长白落叶松生长表现出显著影响的气候因子与其他学者的相关研究结果有所差异(于大炮等, 2005; 于健等, 2016; 陈力等, 2014)。原因可能是不同研究地生长季的起止时间不同，各时间段的气候因子对树木生长的生物学意义也不一致(Savvaetal., 2006)。同时不同研究地的气候条件存在差异，树木生长的限制因子不尽相同(Mäkinenetal., 2003; Wilmkingetal., 2004)。

模型结果显示温度对长白落叶松树木生长影响较大，其中影响程度最大的是7月最高气温和生长季大于5 ℃积温。7月是该地区一年中温度最高的月份，较高的最高气温会促进林内的蒸发散，加剧干旱胁迫，高温又会促进呼吸作用，加速营养物质消耗，不利于树木生长(Lebourgeois, 2000; Wilmkingetal., 2004)。Wang等(2002)在研究中发现树木早材宽度对生长季早期积温敏感，且早材比晚材在树轮宽度的变化中贡献更大，因此，生长季积温的提高对树木生长有促进作用。此外，夏季热湿指数对长白落叶松生长也表现出较大影响，较高的夏季热湿指数意味着在树木生长最旺盛的夏季，温度相对较高而降水不充足，干旱条件会限制树木生长(Wilmkingetal., 2004; Henttonenetal., 2014)。4月平均温、5月降水量、9月最低气温和冬季降水量的变化所引起的树轮宽度指标的变化量相对较小，但也表现出显著影响。温暖的4月有助于土壤解冻及激活形成层，影响木质部生长的开始和持续时间，延长生长季，促进树木生长(Savvaetal., 2006; Jyskeetal., 2012)。5月降水对长白落叶松生长表现出抑制作用，推测原因是此时正处于春季，已消融的积雪补充了土壤水分，因而5月的降水不是树木生长的限制因子。而降水的增多意味着阴天与太阳辐射减少，不利于树木新枝条的生长与光合作用，并且寒冷潮湿的环境会抑制落叶松的生长(Chenetal., 2012)。9月进入了生长季末期，温度较低，树木蒸腾作用减缓，水分对树木生长的限制减少，温度成为树木生长的限制因子，因此表现出9月最低气温升高会促进树木生长(史江峰等, 2006)。冬季降水会在林内形成积雪，可以隔离树根并降低冰冻深度，且积雪融化可以补充土壤水分，有助于树木抵御干旱，有利于其生长(Linderholm, 2001)。

5 结论

本研究表明对长白落叶松云冷杉混交林内的长白落叶松生长有显著影响的因素有： 形成层年龄、林分密度、单木竞争、间伐和气候因子。在林龄20年左右，中高强度的间伐对长白落叶松树木生长有较强的促进作用，促进效果可持续4年。低海拔长白落叶松生长受大量气候因子综合作用，其生长的主要限制因子是7月最高气温和生长季大于5 ℃积温。生长于林分内部的树木对气候因子变化敏感。在全球气候变化的背景下，在对树木生长进行研究时，为使研究结果更精确，应当同时考虑气候因子的影响。

关于气候因子及间伐的交互作用，包括不同间伐措施对树轮宽度与气候因子的响应关系的影响和在不同气候条件下间伐对树轮宽度影响的变化在本研究中并没有很好地体现，将在以后的研究中加以完善。

陈 力, 尹云鹤, 赵东升, 等. 2014. 长白山不同海拔树木生长对气候变化的响应差异. 生态学报, 34(6): 1568-1574.

(Chen L, Yin Y H, Zhao D S,etal. 2014. Climate response of tree growth along altitudinal gradient in the Changbai Mountains, Northeast China. Acta Ecologica Sinica, 34(6): 1568-1574. [in Chinese])

龚固堂, 牛 牧, 慕长龙, 等. 2015. 间伐强度对柏木人工林生长及林下植物的影响. 林业科学， 51(4): 8-15.

(Gong G T, Niu M, Mu C L,etal. 2015. Impact of different tinning intensities on growth ofCupressusfunebrisplantation and understory plants. Scientia Silvea Sinicea, 51(4): 8-15. [in Chinese])

贾 芳, 贾忠奎, 马履一, 等. 2009. 抚育间伐对北京山区油松幼龄人工林水源涵养功能的影响. 水土保持学报, 23(6): 235-239.

(Jia F, Jia Z K, Ma L Y,etal. 2009. Effects of thinning on water conservation of youngPinustabulaeformisplantation in Beijing mountain area. Journal of Soil and Water Conservation, 23(6): 235-239. [in Chinese])

贾忠奎, 公宁宁, 姚 凯, 等. 2012a. 间伐强度对塞罕坝华北落叶松人工林生长进程和生物量的影响. 东北林业大学学报, 40(3): 5-31.

(Jia Z K, Gong N N, Yao K,etal. 2012a. Effects of thinning intensity on the growth and biomass ofLarixprincipis-rupprechtiiplantation in Saihanba, Hebei province. Journal of Northeast Forestry University, 40(3): 5-31. [in Chinese])

贾忠奎, 温志勇, 贾 芳, 等. 2012b. 北京山区侧柏人工林水源涵养功能对抚育间伐的响应. 水土保持学报, 26(1): 62-71.

(Jia Z K, Wen Z Y, Jia F,etal. 2012b. Effects of thinning on water conservation ofPlatycladusorientalisplantation in Beijing mountain area. Journal of Soil and Water Conservation, 26(1): 62-71. [in Chinese])

雷相东, 李永慈, 向 玮. 2009. 基于混合模型的单木断面积生长模型. 林业科学, 45(1): 74-80.

(Lei X D, Li Y C, Xiang W. 2009. Indicidual basal area growth model using multi-level linear mixed model with repeated measures. Scientia Silvea Sinicea, 45(1): 74-80. [in Chinese])

雷相东, 陆元昌, 张会儒, 等. 2005. 抚育间伐对落叶松云冷杉混交林的影响. 林业科学, 41(4): 78-85.

(Lei X D, Lu Y C, Zhang H R,etal. 2005. Effects of thinning on mixed stands ofLarixolgensis,AlbiesnephrolepisandPiceajazoensis. Scientia Silvea Sinicea, 41(4): 78-85. [in Chinese])

李春明, 唐守正. 2010. 基于非线性混合模型的落叶松云冷杉林分断面积模型. 林业科学, 46(7): 106-113.

(Li C M, Tang S Z. 2010. The basal area model of mixed stands ofLarixolgensis,AbiesnephrolepisandPiceajezoensisbased on nonlinear mixed model. Scientia Silvea Sinicea, 46(7): 106-113. [in Chinese])

李国雷, 刘 勇, 甘 敬, 等. 2008. 飞播油松林地土壤酶活性对间伐强度的季节响应. 北京林业大学学报, 30(2): 82-88.

(Li G L, Liu Y, Gan J,etal. 2008. Seasonal response of soil enzyme activity to thinning intensity in aerialPinustabulaeformisstands. Journal of Beijing Forestry University, 30(2): 82-88. [in Chinese])

李宗善, 刘国华, 张齐兵, 等. 2010. 利用树木年轮宽度资料重建川西卧龙地区过去159年夏季温度的变化. 植物生态学报, 34(6): 628-641.

(Li Z S, Liu G H, Zhang Q B,etal. 2010. Tree ring reconstruction of summer temperature variations over the past 159 years in Wolong National Natural Reserve, western Sichuan, China. Chinese Journal of Plant Ecology, 34(6): 628-641. [in Chinese])

马履一, 李春义, 王希群, 等. 2007. 不同强度间伐对北京山区油松生长及其林下植物多样性的影响. 林业科学, 43(5): 1-9.

(Ma L Y, Li C Y, Wang X Q,etal. 2007. Effects of thinning on the growth and diversity of undergrowth ofPinustabulaeformisplantation in Beijing mountain areas. Scientia Silvea Sinicea, 43(5): 1-9. [in Chinese])

莫日根, 张秋良, 吕竟斌, 等. 2013. 抚育间伐对油松及华北落叶松人工林生长量的影响. 内蒙古林业科技, 39(2): 28-31.

(Mo R G, Zhang Q L, Lü J B,etal. 2013. Influence of tending thinning on increment ofPinustabulaeformisandLarixprincipis-rupprechtiiplantation. Journal of Inner Mongolia Forestry Science & Technology, 39(2): 28-31. [in Chinese])

史江峰, 刘 禹, 蔡秋芳, 等. 2006. 油松(Pinustabulaeformis)树轮宽度与气候因子统计相关的生理机制——以贺兰山地区为例. 生态学报, 26(3): 697-705.

(Shi J F, Liu Y, Cai Q F,etal. 2006. A case study of physiological characteristics of statistical correlation betweenPinustabulaeformistree-ring widths and climatic factors. Acta Ecologica Sinica, 26(3): 697-705. [in Chinese])

王丽丽, 邵雪梅, 黄 磊, 等. 2005. 黑龙江漠河兴安落叶松与樟子松树轮生长特性及其对气候的响应. 植物生态学报, 29(3): 380-385.

(Wang L L, Shao X M, Huang L,etal. 2005. Tree-ring characteristics ofLarixgmeliniiandPinussylvestrisvar.mongolicaand their response to climate in Mohe, China. Acta Phytoecologica Sinica, 29(3): 380-385. [in Chinese])

吴祥定, 邵雪梅. 1996. 采用树轮宽度资料分析气候变化对树木生长量影响的尝试. 地理学报, 51(增刊): 92-101.

(Wu X D, Shao X M. 1996. Analysis of the impact of climate change on tree growth by tree-ring width data. Acta Geographica Sinica, 51(supplement): 92-101. [in Chinese])

于大炮, 王顺忠, 唐丽娜, 等. 2005. 长白山北坡落叶松年轮年表及其与气候变化的关系. 应用生态学报, 16(1): 14-20.

(Yu D P, Wang S Z, Tang L N,etal. 2005. Relationship between tree-ring chronology ofLarixolgensisin Changbai Mountains and the climate change. Chinese Journal of Applied Ecology, 16(1): 14-20. [in Chinese])

于 健, 徐倩倩, 刘文慧, 等. 2016. 长白山东坡不同海拔长白落叶松径向生长对气候变化的响应. 植物生态学报, 40(1): 24-35.

(Yu J, Xu Q Q, Liu W H,etal. 2016. Response of radial growth to climate change forLarixolgensisalong an altitudinal gradient on the eastern slope of Changbai Mountain, Northeast China. Chinese Journal of Plant Ecology, 40(1): 24-35. [in Chinese])

Akaike H. 1974. A new look at the statistical model identification. Automatic Control IEEE Transactions on, 19(6): 716-723.

Barbour R J, Fayle D C F, Chauret G. 1994. Breast-height relative density and radial growth in mature jack pine (Pinusbanksiana) for 38 years after thinning. Canadian Journal of Forest Research, 24(12): 2439-2447.

Carrer M, Urbinati C. 2004. Age-dependent tree-ring growth responses to climate inLarixdeciduaandPinuscembra. Ecology, 85(3): 730-740.

Chen F, Yuan Y J, Wei W S,etal. 2012. Climatic response of ring width and maximum latewood density ofLarixsibiricain the Altay Mountains, reveals recent warming trends. Annals of Forest Science, 69(6):723-733.

Filipescu C N, Lowell E C, Koppenaal R,etal. 2014. Modeling regional and climatic variation of wood density and ring width in intensively managed Douglas-fir. Canadian Journal of Forest Research, 44(3): 220-229.

Frank D, Esper J. 2005. Characterization and climate response patterns of a high-elevation, multi-species tree-ring network in the European Alps. Dendrochronologia, 22(2): 107-121.

Gerendiain A Z, Oromi J G, Mehtatalo L,etal. 2012. Effects of cambial age, clone and climatic factors on ring width and ring density in Norway spruse (Piceaabies) in southeastern Finland. Forest Ecology and Management, 263(1): 9-16.

Guillemot J, Klein E K, Davi H,etal. 2015. The effects of thinning intensity and tree size on the growth response to annual climate inCedrusatlantica: a linear mixed modeling approach. Annals of Forest Science, 72(5):651-663.

Guller B. 2007. The effect of thinning treatments on density, MOE, MOR, and maxium crushing strength ofPinusbrutiaTen. wood. Annals of Forest Science, 64(4): 467-475.

Henttonen H M, Mäkinen H, Heiskanen J,etal. 2014. Response of radial increment variation of Scots pine to temperature, precipitation and soil water content along a latitudinal gradient across Finland and Estonia. Agricultural & Forest Meteorology, 198-199:294-308.

Jyske T, Manner M, Mäkinen H,etal. 2012. The effects of artificial soil frost on cambial activity and xylem formation in Norway spruce. Trees, 26(2):405-419.

Kerhoulas L P, Kolb T E, Hurteau M D,etal. 2013. Managing climate change adaptation in forests: a case study from the U.S. Southwest. Journal of Applied Ecology, 50(6): 1311-1320.

Kirdyanov A, Hughes M, Vaganov E,etal. 2003. The importance of early summer temperature and date of snow melt for tree growth in the Siberian Subarctic. Trees, 17(1):61-69.

Koga S, Oda K, Tsutsumi J,etal. 1997. Effect of thinning on the wood structure in annual growth rings of Japanese larch (LarixLeptolepis). Iawa Journal, 18(3): 281-290.

Koga S, Zhang S Y, Begin J. 2002. Effects of precommercial thinning on annual radial growth and wood density in Balsam fir (Abiesbalsamea). Wood and Fiber Science, 34(4): 625-642.

Lebourgeois F. 2000. Climatic signals in earlywood, latewood and total ring width of Corsican pine from western France. Annals of Forest Science, 57(2): 155-164.

Lebourgeois F, Eberle P, Merian P,etal. 2014. Social status-mediated tree-ring responses to climate ofAbiesalbaandFagussylvaticashift in importance with increasing stand basal area. Forest Ecology and Management, 328: 209-218.

Linderholm H W. 2001. Climatic influence on Scots pine growth on dry and wet soils in the central Scandinavian mountains, interpreted from tree-ring widths. Silva Fennica, 35(4): 415-424.

Magruder M, Chhin S, Palik B,etal. 2013. Thinning increases climatic resilience of red pine. Canadian Journal of Forest Research, 43(9): 878-889.

Mäkinen H, Nöjd P, Kahle H P,etal. 2003. Large-scale climatic variability and radial increment variation ofPiceaabies(L.) Karst. in central and northern Europe. Trees, 17(2):173-184.

Mehtatalo L, Peltola H, Kilpelainen A,etal. 2014. The response of basal area growth of Scots pine to thinning: a longitudinal analysis of tree-specific series using a nonlinear mixed-effects model. Forest Science, 60(4): 636-644.

Millar C I, Stephenson N L, Stephens S L. 2007. Climate change and forests of the future: managing in the face of uncertainty. Ecological Applications A Publication of the Ecological Society of America, 17(8): 2145-2151.

Misson L, Guiot J, Nicault A. 2003. Effects of different thinning intensities on drought response in Norway spruce (Piceaabies(L.) Karst.). Forest Ecology and Management, 183(1): 47-60.

Misson L, Vincke C, Devillez F. 2002. Frequency responses of radial growth series after defferent thinning intensities in Norway spruce (Piceaabies(L.) Karst.) stands. Forest Ecology and Management, 177(1): 51-63.

Oliver J, Bogino S, Rathgeber C,etal. 2014. Thinning has a positive effect on growth dynamics and growth-climate relationships in Aleppo pine (Pinushalepensis) trees of different crown classes. Annals of Forest Science, 71(3): 395-404.

Rolland C. 1993. Tree-ring and climate relationships forAbiesalbain the internal Alps. Tree-ring Bulletin, 53: 1-11.

Savva Y, Oleksyn J, Reich P B,etal. 2006. Interannual growth response of Norway spruce to climate along an altitudinal gradient in the Tatra Mountains, Poland. Trees, 20(6): 735-746.

Shen C, Lei X, Liu H,etal. 2015. Potential impacts of regional climate change on site productivity ofLarixolgensisplantations in northeast China. iForest-Biogeosciences and Forestry, 8: 642-651.

Szeicz J M, Macdonald G M. 1994. Age-dependent tree-ring growth responses of subarctic white spruce to climate. Canadian Journal of Forest Research, 24(1):120-132.

Wang L, Payette S, Bégin Y. 2002. Relationships between anatomical and densitometric characteristics of black spruce and summer temperature at tree line in northern Quebec. Canadian Journal of Forest Research, 32(3): 477-486

Wilmking M, Juday G P, Barber V A,etal. 2004. Recent climate warming forces contrasting growth responses of white spruces of white spruce at treeline in Alaska through temperature thresholds. Global Change Biology, 10(10): 1724-1736.

Wimmer R, Grabner M. 1997. Effects of climate on vertical resin duct density and radial growth of Norway spruce (Piceaabies(L.) Karst.). Trees, 11(5): 271-276.