To reveal the radial growth trends of Pinus koraiensis and the differences in their responses to climate factors among different latitudes, and to assess the dynamic characteristics, adaptability and sensitivity of P. koraiensis under the influence of climate change, dendrochronological techniques were used to study the relationships between the climatic variables and the radial growth of P. koraiensis in broadleaved Korean pine forests with different latitudes. The results showed that there were differences in the responses of the radial growth of P. koraiensis in four different latitudes to local climatic factors. In Baishilazi Nature Reserve, the southernmost sample plot, the radial growth of P. koraiensis was significantly positively correlated to the average relative humidity but significantly negatively correlated to the average maximum temperature of growing season. In low altitude of Changbai Mountain Nature Reserve, the intermediate sample plot, it was significantly positively correlated with precipitation, average relative humidity and PDSI, and significantly negatively correlated with the average maximum temperature in growing season. In Liangshui Nature Reserve, the intermediate sample plot, significantly positive correlations occurred for average relative humidity and PDSI in growing season, while significantly negative correlations occurred for the average temperature and the average maximum temperature in growing season. Nevertheless, in Shengshan Nature Reserve, the northernmost sample plot, it was positively correlated to the temperature factors of most months. The climatic variables in June of current year were the main factors limiting the radial growth of P. koraiensis in all latitudes, and the average maximum temperature in June of current year was highly negatively correlated with radial growth for P. koraiensis measured at all sites. In recent four decades, with the rising of temperature, the radial growth of P. koraiensis decreased significantly in the southernmost point, increased significantly in the northernmost point, and did not change significantly in middle latitudes. The distribution area of P. koraiensis would be reduced if the temperature increased and precipitation remained unchanged in the future.
全 文 :不同纬度阔叶红松林红松径向生长
对气候因子的响应
刘 敏1,2 毛子军1∗ 厉 悦2 孙 涛1 李兴欢1 黄 唯1 刘瑞鹏1 李元昊2
( 1东北林业大学森林植物生态学教育部重点实验室, 哈尔滨 150040; 2齐齐哈尔大学, 黑龙江齐齐哈尔 161006)
摘 要 运用树木年轮气候学方法,研究了不同纬度阔叶红松林红松径向生长趋势及其与气
候因子的关系,以期揭示不同纬度红松径向生长对气候因子响应的差异以及气候变化影响下
红松的动态特征、适应性及敏感性.结果表明: 不同纬度红松径向生长对当地气候因子的响应
存在差异,最南部的白石砬子自然保护区红松径向生长与生长季的平均相对湿度呈显著正相
关,与平均最高气温呈显著负相关.中部的长白山自然保护区的低海拔区域红松径向生长与
生长季的水分因子(降水、相对湿度和帕尔默干旱指数 PDSI)呈显著正相关,而与平均最高气
温呈显著负相关;凉水自然保护区红松径向生长与生长季的水分因子(相对湿度和 PDSI)呈
显著正相关,与气温因子(平均气温和平均最高气温)呈显著负相关.而最北部的胜山自然保
护区红松径向生长则与大部分月份的气温因子呈显著正相关.当年 6 月气候因子是影响所有
纬度红松径向生长的关键气候因子,4个样地都与当年 6月平均最高气温呈显著负相关.在气
温不断上升的近 40年,最南部的红松径向生长呈显著下降趋势,最北部呈显著上升趋势,中
部的变化不显著.如果未来气温升高而降水不变,红松分布区可能缩小.
关键词 红松; 树木年轮; 纬度; 气候变化
本文由国家自然科学基金面上项目(31270494)、国家重点基础研究发展计划项目(2010CB951301)、齐齐哈尔大学青年教师科研启动支持计划
项目(2012k⁃M23)和黑龙江省大学生创新创业训练计划项目(201410221030)资助 This work was supported by the National Natural Science Foun⁃
dation of China (31270494), the National Basic Research Program of China (2010CB951301), the Young Teachers’ Scientific Research Support Pro⁃
gram of Qiqihar University ( 2012k⁃M23) and the Innovation and Entrepreneurship Training Program for College Students in Heilongjiang Province
(201410221030).
2015⁃10⁃27 Received, 2016⁃02⁃23 Accepted.
∗通讯作者 Corresponding author. E⁃mail: zijunm@ nefu.edu.com
Response of radial growth of Pinus koraiensis in broad⁃leaved Korean pine forests with
different latitudes to climatical factors. LIU Min1,2, MAO Zi⁃jun1∗, LI Yue2, SUN Tao1, LI
Xing⁃huan1, HUANG Wei1, LIU Rui⁃peng1, LI Yuan⁃hao2 ( 1Ministry of Education Key Laboratory
of Forest Plant Ecology, Northeast Forestry University, Harbin 150040, China; 2Qiqihar University,
Qiqihar 161006, Heilongjiang, China) .
Abstract: To reveal the radial growth trends of Pinus koraiensis and the differences in their respon⁃
ses to climate factors among different latitudes, and to assess the dynamic characteristics, adaptabi⁃
lity and sensitivity of P. koraiensis under the influence of climate change, dendrochronological tech⁃
niques were used to study the relationships between the climatic variables and the radial growth of
P. koraiensis in broad⁃leaved Korean pine forests with different latitudes. The results showed that
there were differences in the responses of the radial growth of P. koraiensis in four different latitudes
to local climatic factors. In Baishilazi Nature Reserve, the southernmost sample plot, the radial
growth of P. koraiensis was significantly positively correlated to the average relative humidity but sig⁃
nificantly negatively correlated to the average maximum temperature of growing season. In low alti⁃
tude of Changbai Mountain Nature Reserve, the intermediate sample plot, it was significantly posi⁃
tively correlated with precipitation, average relative humidity and PDSI, and significantly negatively
correlated with the average maximum temperature in growing season. In Liangshui Nature Reserve,
the intermediate sample plot, significantly positive correlations occurred for average relative humidi⁃
ty and PDSI in growing season, while significantly negative correlations occurred for the average
应 用 生 态 学 报 2016年 5月 第 27卷 第 5期 http: / / www.cjae.net
Chinese Journal of Applied Ecology, May 2016, 27(5): 1341-1352 DOI: 10.13287 / j.1001-9332.201605.020
temperature and the average maximum temperature in growing season. Nevertheless, in Shengshan
Nature Reserve, the northernmost sample plot, it was positively correlated to the tempe⁃rature fac⁃
tors of most months. The climatic variables in June of current year were the main factors limiting the
radial growth of P. koraiensis in all latitudes, and the average maximum temperature in June of cur⁃
rent year was highly negatively correlated with radial growth for P. koraiensis measured at all sites.
In recent four decades, with the rising of temperature, the radial growth of P. koraiensis decreased
significantly in the southernmost point, increased significantly in the northernmost point, and did
not change significantly in middle latitudes. The distribution area of P. koraiensis would be reduced
if the temperature increased and precipitation remained unchanged in the future.
Key words: Pinus koraiensis; tree ring; latitude; climate change.
红松(Pinus koraiensis)是我国名贵又珍稀的树
种,以该种为主要建群种或优势种的阔叶红松林是
我国东北地区的典型地带性植被,阔叶红松林的动
态变化关系到东北地区森林植被的稳定.由于红松
属典型温带湿润型山地大乔木树种,对温湿状况适
应的生态幅较窄[1] .气候变化情景下,红松的动态变
化引起了很多学者的关注[2-5] .一些学者对长白山等
地红松与气候的响应进行了研究[2,5-13],结果显示,
不同海拔红松对气候因子的响应存在明显差
异[9-11] .但是,气候因子对不同纬度红松径向生长的
影响是否也存在差异? 每个纬度影响红松径向生长
的关键气候因子是什么? 在气候变暖的情景下,不
同纬度红松的径向生长的变化趋势是否一致? 这些
问题目前尚未明了.
气候变化对树木的影响是直接作用于树木生理
过程,树木每一个年轮的形成都是多个生态因子综
合作用的结果,每个年轮的宽窄都记录着众多气象
因子和生态因子的信息.气候敏感区域的树木,运用
树木年轮学方法在去除生长趋势和其他非气候因素
的影响后,树轮宽度年表可保留非常强的气候信
号[14-15],因此可利用树轮宽度年表与气候因子的相
关性来揭示影响树木径向生长的关键气候因子[16] .
树木年代学由于具备定年准确、重复性高等特点,已
经成为分析气候敏感区典型树种径向生长与气候变
化关系的重要手段与主要途径[17] .
本文运用树木年轮学方法,以我国温带针阔混
交林内 4个不同纬度原始阔叶红松林的红松为研究
对象,研究不同纬度红松径向生长对气候因子响应
的异同,不同纬度影响红松径向生长的关键气候因
子及气候变暖背景下不同纬度红松的径向生长动态
特征、适应性及敏感性,以期为预测全球气候变化背
景下我国原始阔叶红松林的动态提供数据,为阔叶
红松林的合理经营提供科学参考.
1 研究地区与研究方法
1 1 研究区概况
在我国原始阔叶红松林分布区从南到北 4个以
阔叶红松林为保护对象的国家级自然保护区内选取
样地(图 1),自然保护区受人为干扰小,研究结果主
要反映自然因素对植物生长的影响.
白石砬子国家级自然保护区(40°50′—40°57′ N,
124°44′—124°57′ E)位于辽宁省东部宽甸县大川头
镇,保护区内有较完整的大面积东北亚地区地带性
原生型阔叶红松混交林[19] .此区域属温带季风气
候,冬季比较寒冷,夏季温暖湿润,冬夏、春秋昼夜的
温差变化较大,1958—2012年(宽甸气象站)的年平
均气温为 7.1 ℃,年总降水量 1097.1 mm,雨热同季.
长白山国家级自然保护区 ( 42° 04′—42° 23′ N,
127°55′—128°08′ E)位于吉林省延吉地区安图县、
抚松县、长白县境内,是我国目前面积最大、自然环
境和生态系统保存比较完整的森林保护区[20],红松
在保护区内广泛分布[21] .该区属温带大陆性山地气
图 1 采样地和气象站点分布
Fig.1 Distribution of sampling sites and weather stations.
Ⅰ: 样地 Sampling site; Ⅱ: 气象站 Weather station; Ⅲ: 阔叶红松林
分布区[18]Broad⁃leaved Korean pine forest[18] . A: 白石砬子自然保护
区 Baishilazi Nature Reserve; B: 长白山自然保护区 Changbai Moun⁃
tain Nature Reserve; C:凉水自然保护区 Liangshui Nature Reserve; D:
胜山自然保护区 Shengshan Nature Reserve. 下同 The same below.
2431 应 用 生 态 学 报 27卷
候,冬季漫长寒冷,夏季短暂且温暖多雨,春季风大
干燥,秋季凉爽多雾. 1958—2012 年(松江气象站)
的年平均气温 2.8 ℃,年降雨量 674 mm.凉水国家级
自然保护区(47°7′—47°14′ N,128°48′—128°55′ E)
位于我国小兴安岭南坡达里带岭支脉东坡,属于黑
龙江省伊春市,是目前我国保存下来的最典型、最完
整的北温带针阔叶混交林系统[22-23],属于现代红松
分布中心[24] .此区域属大陆性季风气候, 1958—
2012年(伊春气象站)的年平均气温 1.1 ℃,年降水
量 629. 5 mm.胜山国家级自然保护区 ( 49° 25′—
49°40′ N,126°27′—127°02′ E)地处小兴安岭西北
坡,处于大、小兴安岭交错过渡地带,是我国原始红
松林分布的最北区域,也是阔叶红松林分布的最西
北区域[25] .此区域属寒温带大陆性气候,冬季寒冷
漫长,夏季温暖多雨,秋季短促温湿.1958—2012 年
(孙吴气象站)的年平均气温 - 0. 5 ℃,年降水量
537.7 mm.
1 2 采样与年表建立
1 2 1样本采集与处理 2012 年 9—10 月,在 4 个
国家级自然保护区内设置 4个样地(表 1).
按照国际树轮库(ITRDB)的标准,分别在 4 个
样地内选择达到林冠层的生长健康的成熟红松进行
样芯采集,每个样地选取 25 株以上样木,测定每株
样木的胸径,用树木生长锥在每株样木胸高处南北
相对方向钻取 2个样芯,确保通过髓心.树芯取出后
放入塑料管内,用标签纸封口并注明采样地、树号、
树芯编号和样地的生境条件及各采样木的特征,带
回室内固定、自然风干、打磨.
1 2 2年轮宽度测定及年表建立 用骨架示意图方
法[26]进行初步交叉定年[14],初步定年后按照树木
年代学方法对树芯进行标年,整 10、50 和 100 年处
各自做出标注.利用 LinTab 5 年轮分析仪在 0.001
mm水平上测量树轮宽度.由于采样时间某些样地红
松还处于生长期末期,因此舍弃当年(2012 年)的年
轮,从 2011 年开始测定年轮宽度.应用 TSAP 和
COFECHA[27]软件对样芯宽度序列进行交叉定年,
确保定年与测量的准确.去除腐心、断裂、敏感度低
和奇异点较多的序列.
运用 ARSTAN程序[28]对序列宽度进行标准化
并建立年轮宽度指数年表.由于本研究区属于较湿
润的森林地区,研究中采用样芯序列长度的 2 / 3 作
为步长的样条函数法(该方法无需假定年轮样本生
长趋势的变化形式,直接采用连续光滑插值方法对
具有持续性生长以及种间竞争产生的非同步扰动的
树木进行生长趋势拟合,适于湿润地区去生长趋
势[29-31])去除红松生长趋势和非气候因素,并对去
趋势的序列以双重平均法合成标准年表(STD).4 个
样地中白石砬子自然保护区红松的年表较短,对其
进行 1932—2011年(80年)公共区间分析,对其他 3
个样地红松标准年表进行 1912—2011 年(100 年)
公共区间分析,得出相关年表统计指标,如样本总体
代表性、树间相关系数、一阶自相关.
1 2 3气象资料 分别选取与 4 个样地最近的宽甸
气象站(海拔 260. 1 m)、松江气象站(海拔 591. 4
m)、伊春气象站(海拔 240.9 m)和孙吴气象站(海
拔 234.5 m)的月降水量、月均相对湿度、月平均气
温、月平均最高和最低气温,此数据来自中国气象科
学数据共享网服务网(http: / / cdc.cma.gov.cn).鉴于
气象站海拔与样地海拔存在差异,根据海拔每上升
100 m 气温下降 0. 6 ℃ 的一般规律计算样地气
温[32],得出样地月平均气温、月平均最高和最低气
温.月降水量、月平均相对湿度数据采用气象站数
据.4个气象站的位置和样地气象因子分别见图 1和
图 2.
表 1 采样地基本情况
Table 1 Basic status of sampling sites
样地
Site
经度
Longitude
纬度
Latitude
海拔
Altitude (m)
坡向
Aspect
郁闭度
Crown density
样本量
Sample size
主要植物
Main plant species
A 126°46′ 49°27′ 560~590 北坡 North 0.7 31 1、2、3、4、5、6、7、8、9
B 128°54′ 47°11′ 390~410 北坡 North 0.7 32 1、3、4、8、10、11、12
C 128°05′ 42°24′ 740~750 北坡 North 0.6 44 1、2、4、6、13、14、15
D 124°47′ 40°55′ 790~820 北坡 North 0.8 37 1、2、4、10、12、16、17、18、19
A: 胜山自然保护区 Shengshan Nature Reserve; B: 凉水自然保护区 Liangshui Nature Reserve; C: 长白山自然保护区 Changbai Mountain Nature
Reserve; D: 白石砬子自然保护区 Baishilazi Nature Reserve. 下同 The same below. 1) 红松 Pinus koraiensis; 2) 鱼鳞云杉 Picea jezoensis; 3) 臭冷
杉 Abies nephrolepis; 4) 紫锻 Tilia amurensis; 5) 白桦 Betula platyphylla; 6) 糠椴 Tilia mandshurica; 7) 狗枣猕猴桃 Actinidia kolomikta; 8) 毛榛
Corylus mandshurica; 9) 山葡萄 Vitis amurensis; 10) 花楷槭 Acer ukurunduense; 11) 五角枫 Acer mono; 12) 刺五加 Eleutherococcus senticosus; 13)
蒙古栎 Quercus mongolica; 14) 水曲柳 Fraxinus mandshurica; 15) 山杨 Populus davidiana; 16) 裂叶榆 Ulmus laciniata; 17) 核桃楸 Juglans mand⁃
shurica; 18) 枫桦 Betula costata; 19) 忍冬 Lonicera japonica.
34315期 刘 敏等: 不同纬度阔叶红松林红松径向生长对气候因子的响应
图 2 4个样地 1958—2012年月平均降水量、月平均气温、月平均最高和最低气温以及季节划分
Fig.2 Monthly average precipitation, average temperature, average maximum and minimum temperature and season division in the
four sampling sites in 1958—2012.
P: 平均降水量 Average precipitation; Tm:平均气温 Average temperature; Tmax:平均最高气温 Average maximum temperature; Tmin:平均最低气温
Average minimum temperature. WD:冬季Winter dormant period; BG:生长季前期 Early growing season; GS:生长季 Growing season; EG:生长季末
期 End of growing season. 下同 The same below.
从荷兰皇家气象研究所的数据共享网站( ht⁃
tp: / / climexp.knmi.nl)获取帕尔默干旱指数(PDSI)
(2.5° × 2. 5°的格点数据) [33] . PDSI 数据时间段为
1981—2005年. 4 个样地 PDSI 数据的经纬度分别
为:白石砬子国家级自然保护区(39. 0—41. 5° N,
123. 5—126. 0° E )、长白山国家级自然保护区
(41.5—44.0° N, 126.5—129.0° E)、凉水国家级自然
保护区(46.5—49.0° N ,128.5—131.0° E)、胜山国家
级自然保护区(48.5—51.0° N ,126.0—128.5° E).
考虑到红松生长节律和气候因子对红松生长影
响的滞后作用,选取上一年 4 月到当年 10 月(共 19
个月)的月气象数据进行分析.另外,考虑到气候因
子具有累积和长期影响效应,根据月平均气温(Tm)
和月平均最低气温(Tmin),将年气候资料分为以下
几个季节阶段[10,23]:冬季(WD,Tm≤0 ℃)、生长季
前期(BG,Tm >0 ℃且 Tmin <5 ℃)、生长季(GS,Tmin
≥5 ℃)、生长季末期(EG,Tm > 0 ℃且 Tmin < 5 ℃)
(图 2).
1 3 数据处理
采用 SPSS 19.0 软件 Pearson 相关法进行红松
树轮年表与月气象数据和季节数据的相关分析[34]
(α= 0.05).
2 结果与分析
2 1 不同纬度红松年表特征
表 2 显示,白石砬子自然保护区红松年表序列
最短,只有 102年;而凉水自然保护区红松年表序列
最长,达 238 年. 4 个样地红松的平均敏感度均
>0.170,其中,白石砬子自然保护区最高,为 0.197.
各采样地的标准差为 0.172 ~ 0.537,白石砬子自然
保护区最大,而长白山自然保护区最低.4 个样地的
信噪比都较高,为 8.894~21.505.各样地的样本总体
表 2 4个样地红松年表特征及公共区间分析
Table 2 Chronology characteristics and common analysis of STD chronology of Pinus koraiensis in four sites
样地
Site
样芯长度
Length of series
(a)
公共区间
Common period
(a)
平均敏感度
Mean
sensitivity
标准差
Standard
deviation
信噪比
Signal to
noise ratio
样本总体
代表性
EPS
树间相关
系数
Coefficient
between
trees
一阶
自相关
Autocorrelation
order 1
EPS
>0.85
A 160 (1852—2011) 100 (1912—2011) 0.173 0.302 8.894 0.899 0.270 0.807 1895—2011
B 238 (1774—2011) 100 (1912—2011) 0.182 0.318 13.051 0.929 0.273 0.472 1829—2011
C 154 (1858—2011) 100 (1912—2011) 0.170 0.172 11.978 0.923 0.295 0.253 1895—2011
D 102 (1910—2011) 80 (1932—2011) 0.197 0.537 21.505 0.956 0.577 0.625 1927—2011
EPS: General representative of samples.
4431 应 用 生 态 学 报 27卷
代表性都超过了 0.85 的标准值.4 个样地树间相关
系数为 0.270~0.577,胜山自然保护区最低,白石砬
子自然保护区的最高.胜山自然保护区红松的一阶
自相关最高,为 0.807,表明胜山自然保护区红松径
向生长受上一年气候因素的影响最大.
2 2 红松径向生长与气候因子的关系
图 3 显示,不同纬度红松径向生长对当地的气
候因子响应情况存在差异,白石砬子自然保护区主
要对水分(相对湿度和降水)和高温比较敏感;长白
山自然保护区对水分比气温更敏感;凉水自然保护
区主要对当年 6 月的气候因子敏感;胜山自然保护
区主要对气温因子敏感.具体表现为:白石砬子自然
保护区红松径向生长与上一年 5月的降水量和月均
相对湿度、当年 5—7月的月均相对湿度呈显著正相
关,与当年 1月和 9月的降水量、当年 6月的月均气
温和月均最高气温呈显著负相关;长白山自然保护
区红松径向生长与当年 4月和 9月的月平均气温、6
月的降水量和月均相对湿度及 6—7、9—10 月的
PDSI 呈显著正相关,与上一年 7月的降水量、当年 6
月的月最高气温和 4月的月均相对湿度呈显著负相
图 3 红松年轮宽度与月份气象因子的相关系数
Fig.3 Correlation coefficients between STD chronology of Pinus koraiensis and monthly meteorological data.
p4~ p12: 上一年 4—12 月 Previous April-December; 1 ~ 10: 当年 1—10 月 Current year January-October. RH: 平均相对湿度 Average relative
humidity; PDSI: 平均帕尔默干旱指数 Average Palmer drought index. ∗P<0.05; ∗∗P<0.01. 下同 The same below.
54315期 刘 敏等: 不同纬度阔叶红松林红松径向生长对气候因子的响应
关;凉水自然保护区红松径向生长与当年 6 月的月
平均气温、月均最高和最低气温及当年 2 月的月均
相对湿度呈显著负相关,与当年 6 月的降水和月均
相对湿度、当年 6—7 月的 PDSI 及上一年 4 月的月
均气温和月均最高气温呈显著正相关;胜山自然保
护区红松径向生长与上一年 4 月、7—8 月、10—12
月和当年 2—5 月、7—9 月的月均气温,当年 2—3
月的月均最高气温,上一年 4 月—当年 5 月、当年
7—9月的月均最低气温,上一年5月的月均相对湿
图 4 红松年表与季节气象因子的相关系数
Fig.4 Correlation coefficients for STD chronology of Pinus ko⁃
raiensis and seasonal meteorological data.
PBG: 上一年生长季前期 Previous early growing season; PGS: 上一年
生长季 Previous growing season; PEG: 上一年生长季末期 Previous
end of growing season.
度和当年 6月的 PDSI呈显著正相关,与当年 6月的
月均最高气温呈显著负相关.
与月气候变量的关系相比,红松径向生长与季
节气候变量的关系要更为清晰.从图 4 可以看出,白
石砬子自然保护区的红松径向生长主要与当年生长
季的季平均气温和平均最高气温呈显著负相关,与
季均相对湿度呈显著正相关;长白山自然保护区主
要与当年生长季的降水和相对湿度以及生长季和生
长季末期的 PDSI 呈显著正相关,与生长季平均最
高气温呈显著负相关;凉水自然保护区主要与生长
季的平均相对湿度和 PDSI 呈显著正相关,与生长
季的季均气温和季均最高气温呈显著负相关;胜山
图 5 1958—2011年 4个气象站的年平均气温、年均最高气
温、年均最低气温和年降水量动态
Fig.5 Dynamics of annual average temperature, annual average
maximum temperature, annual average minimum temperature,
and annual precipitation from 1958 to 2011.
6431 应 用 生 态 学 报 27卷
自然保护区的红松径向生长对每个季节的大部分气
温因子呈显著正相关.
2 3 红松径向生长对气候变暖的响应
由图 5可以看出,4个样地年平均气温、年平均
最高和最低气温在 1970 年之后都出现了显著上升
(长白山自然保护区年平均最高气温除上升未达显
著水平).其中,年均最低气温上升幅度最大,胜山自
然保护区和长白山自然保护区上升速率分别为 1.20
和 0.80 ℃·(10 a) -1,凉水自然保护区上升幅度在
4个样地中较小,但也达到了 0.37 ℃·(10 a) -1 .4
个样地的年平均气温上升幅度也较明显,每 10年最
少上升 0.36 ℃,其中,胜山自然保护区每十年上升
0.73 ℃.除长白山自然保护区外,其他 3 个样地的年
平均最高气温上升显著,上升幅度为 0. 13 ~ 0. 32
℃·(10 a) -1 .相对于显著上升的气温,4 个样地的
降水变化都不显著.
从不同纬度样地红松径向生长对气候因子的响
应情况可以看出,4 个纬度红松每年年轮宽度的变
化与气温因子关系密切.在气温不断上升的近 40
年,不同纬度红松径向生长的变化趋势各异.从图 6
可以看出,1970—2011年,最南端的白石砬子自然保
护区的红松径向生长呈显著下降趋势,年轮指数下降
图 6 各样地红松标准年表和样本量
Fig.6 Standard chronologies and sample size of Pinus koraiensis in each site.
Ⅰ: 年轮宽度指数 Ring width index; Ⅱ: 平均值线 Mean line; Ⅲ: 样本数量 Sample size.
74315期 刘 敏等: 不同纬度阔叶红松林红松径向生长对气候因子的响应
速率为 0.38·(10 a) -1;最北端的胜山自然保护区呈
显著上升趋势,年轮指数上升速率为 0.07·(10 a) -1;
中间 2个纬度样地红松年表在近 40年变化不显著.
3 讨 论
3 1 不同纬度红松径向生长对气候因子响应的
异同
植物群落的分布与自然地理环境条件有着极为
密切的关系,地带性植被的形成是与该地区气候条
件长期适应的结果,所以植被分布区的中心地段为
最适分布区,而分布区的边缘区域则是受某些环境
因子制约的生态胁迫区.本研究中,4 个样地包含了
原始阔叶红松林分布区的中心区域和南、北缘,跨越
近 10个纬度,不同纬度的生态环境条件存在一定差
异.不同纬度红松径向生长对气候因子的响应有一
致的地方,但也存在很多差异.
不同之处主要为:最南端的白石砬子自然保护
区红松径向生长主要对生长季的高温和降水敏感,
而最北端的胜山自然保护区对每个月的气温因子都
很敏感;中间纬度的长白山低海拔地区主要对生长
季水分因子敏感,凉水自然保护区主要对当年 6 月
的气温和水分因子敏感.处于最低纬度的白石砬子
自然保护区,虽然降水量最多,但其气温最高,高温
使土壤蒸发量增大,其 PDSI 在 4 个纬度样地中最
低,说明土壤干旱最严重.在气温较高但降水量偏少
的月份,高温将促进红松蒸腾作用和呼吸作用,增强
气温的胁迫,加之高温造成的土壤持水能力下降,进
一步抑制了红松的生长.因此,此区域生长季高温及
高温造成的干旱是红松径向生长的主要限制因子.
处于最高纬度的胜山自然保护区,气温低、降水量
少,月均相对湿度和 PDSI 较低.对于高纬度寒冷地
区来说,树木的径向生长更大程度上受有机物合成
的影响,生长期间过低的气温会降低叶片的光合作
用,从而抑制红松的径向生长.因此,本地区红松的
径向生长与大部分月份(季节)的平均气温和平均
最低气温均呈显著正相关,此结果与及莹[35]对黑河
红松的研究结果一致.处于中间纬度的长白山自然
保护区和凉水自然保护区,气温和降水量都处于中
间水平.长白山自然保护区低海拔地区生长季气温
较适合红松生长,降水相对较少,造成红松对生长季
的降水量、平均相对湿度和 PDSI 都呈显著正相关,
与最高气温呈显著负相关,此研究结果与 Yu 等[10]
和高露双等[31]对本地区红松的研究结果一致.凉水
自然保护区 6 月的平均气温和平均最高气温都很
高,但是 6月的降水量相对于红松快速生长期的需
求量来说不是很充足,因此此区域红松径向生长与
6月的降水量呈显著正相关;另外,6 月的高温会加
大土壤水分蒸发,降低土壤湿度,在降水量不是特别
充足的 6 月,高温会直接限制红松径向生长,因此
红松径向生长与 6 月的平均气温、平均最高气温呈
显著负相关.此研究结果与及莹[35]对凉水红松的研
究结果一致.通常在寒冷地区,气温与不同纬度的同
种植物的径向生长的相关性存在差异,主要是因为
在其分布区的北端(或上限),较高的气温有利于物
种度过寒冷的冬天[36],高温还能促进融雪和生长季
的提前[37],并提高生长季的植物生长速率[38];而在
物种分布的南端(或下限),生长季的最高温已接近
物种可以容忍的上限,因此,更高的气温只会增强胁
迫和阻碍生长速率[39] .
共同之处主要表现为 2 个方面:一是 4 个不同
纬度红松径向生长都与当年 6月的平均最高气温呈
显著负相关.将 4 个样地当年 6 月的平均最高气温
的 10年滑动平均序列与当地红松标准年表指数的
10年滑动平均序列进行比较(图 7),发现它们之间
存在明显的负相关性,凉水自然保护区和白石砬子
自然保护区呈显著负相关,长白山自然保护区呈显
著负相关,胜山自然保护区的相关性值偏小,说明当
年 6月的平均最高气温对 4个样地红松的年轮宽度
起重要作用.1970年以来,只有白石砬子自然保护区
的 6月的平均最高气温显著上升[上升速率为 0.36
℃·(10 a) -1],而其他 3个样地的上升幅度均未达
到显著水平.另外,近 40 年来白石砬子自然保护区
的年表显著下降,而长白山和凉水自然保护区变化
不显著(图 6),这也从另一个角度印证了 6 月最高
气温对红松径向生长的影响.
二是,4 个样地红松径向生长除了与 6 月的平
均最高气温响应敏感外,对 6 月的其他气象因子也
都比较敏感:与 6月的平均气温均呈负相关,与 6 月
的降水量、月均相对湿度和 PDSI 均呈正相关,并且
相关系数值都较大,表明 6 月的气象因子对红松的
径向生长影响很大.4 个样地 6 月的气温都较高,是
红松径向生长迅速时期,但 6 月各地的降水量相对
生长需要来说还不充足,6 月的高温加重了土壤水
分的匮乏.土壤干旱胁迫会降低树木叶水势,使叶肉
细胞对水分的吸附力增大,限制了叶内水分的运动,
从而使树木的光合作用速率减弱,影响到树木的生
长[40] .此外,红松是浅根性树种,调查样地内大部分
红松都是生长于坡地、浅土层,红松对土壤水分非常
8431 应 用 生 态 学 报 27卷
图 7 4个样地红松年表与 6 月平均最高气温的 10 年滑动
平均序列
Fig.7 Ten years moving average sequence of STD chronologies
of Pinus koraiensis and the average maximum temperature in
June in the four sampling sites.
Ⅰ: 年轮宽度指数 Ring width index; Ⅱ: 6 月平均最高气温 Mean
maximum temperature of June.
敏感.因此,红松与 6 月的气温因子呈负相关,与水
分因子呈正相关.
3 2 不同纬度红松径向生长对气候变暖的响应
有研究发现,不同地区红松在 1970年后随着气
候变暖年轮宽度指数发生了明显变化.如 Yu 等[10]
研究发现,长白山低海拔红松从 1970年后年轮宽度
指数显著下降,而中高海拔显著上升.李广起等[2]也
发现,长白山高海拔红松年轮指数在 1980年后显著
上升;及莹[35]研究显示,黑河的红松宽度指数在
1980年后表现为显著下降趋势.
本研究样地贯穿我国原始阔叶红松林分布区从
南到北不同纬度,4个样地的气温在 1970 年之后显
著上升(降水变化都不显著),而红松年轮宽度指数
在 1970年之后变化趋势各异,南端显著下降,中间
变化不显著,北端显著上升.处于最南端的白石砬子
自然保护区近 40年气温上升幅度显著,而红松年表
指数呈显著下降趋势.表明在红松分布区的南缘高
温地区,气温上升显著阻碍了红松的径向生长.近 40
年,白石砬子自然保护区的年均气温、年均最高和最
低气温分别以 0.26、0.33、0.43 ℃·(10 a) -1的速率
上升(图 5).其中,6月的气温上升尤其显著,2000—
2009年 6月的月均和月均最高气温分别比 1970—
1979年上升了 1.2和 1.5 ℃ .而相关分析显示,此区
域红松的径向生长与当年 6月的月均气温和月均最
高气温呈显著负相关.白天最高气温的增加促进相
对湿度降低,同时造成蒸腾加强引起植物体内水份
失调,从而导致生长量降低[41] .较高的 6 月的高温
和不太充足的降水,使本来处于径向生长旺季的红
松的径向生长减速甚至停止,造成年轮宽度指数随
着气温的升高而逐渐下降.可以预测,如果此区域的
气温持续上升,尤其是 6月的气温上升,将会使红松
的径向生长进一步下降,进而影响到红松的分布,使
红松分布区南缘发生北移.分布于中部纬度的长白
山自然保护区低海拔地区和凉水自然保护区的红松
年轮宽度在近 40 年变化不显著.虽然近 40 年这 2
个区域年均气温、年均最高和最低气温都上升显著,
但与红松径向生长呈显著负相关的关键月份的气温
(长白山自然保护区的 6 月平均最高气温;凉水自
然保护区的 6 月平均气温、平均最高和最低气温)
的升高幅度并未达到显著水平,因此 1970年之后的
气候变暖现象对此区域的红松径向生长影响不显
著.另外,长白山和小兴安岭是红松分布的核心区
域,尤其凉水自然保护区是现代红松分布的中心区
域[24],这 2个地区的气候环境应该是最适合红松生
长的.红松径向生长动态显示至今为止本区域的气
候变化幅度还在红松对各因子的适宜生态幅范围
内,未对其径向生长造成明显的扰动.凉水自然保护
区的红松年轮宽度呈不显著的下降,这与及莹[35]对
凉水红松的研究结果一致,与延晓冬[42]、程肖侠
等[43]模型预测的快速衰退存在一定差异.长白山自
然保护区的结果与 Yu 等[10]的研究结果有差异,他
们的研究显示,长白山海拔 750~800 m地区红松径
94315期 刘 敏等: 不同纬度阔叶红松林红松径向生长对气候因子的响应
向生长在 1970年后显著下降,而本文研究结果显示
无显著变化.结果差异的原因可能与采样的具体样
地、生境等有关.处于最北端的胜山自然保护区,近
40年气温上升幅度最大,红松年轮指数也在 1970
年后表现为显著上升趋势.此区域地处严寒地区,红
松径向生长与大部分月份的气温因子都呈显著正相
关,气温显著升高,有利于红松的光合作用和有机物
的积累,促进了红松的径向生长.本研究结果与及
莹[35]的研究结果相反,及莹[35]认为黑河红松的年
轮宽度在 1980年之后呈下降趋势.但长白山红松垂
直分布的树轮研究表明,最近几十年,随着气温不断
升高,相对于低海拔分布的红松,高海拔地段红松年
轮宽度指数显著上升[2,10] .可见,气温升高有利于寒
冷区域红松的生长,这与本研究结果一致.
很多模型模拟显示,在气候持续变化的情况下
东北阔叶红松林可能有一个快速衰退的过程,红松
的分布区可能会改变[42-45] .鉴于以上 4 个纬度红松
的径向生长在气候变暖后的动态可以预测,如果气
温持续升高,尤其是 6月平均最高气温持续升高,低
纬度的红松生长将首先受到影响,生长速率将减慢,
会造成原始红松林的退化,阔叶红松林分布区的南
缘可能发生北移;中部纬度地区红松受到的影响暂
时不明显;最北地区红松径向生长会有所提高.但
是,是否能够北移,目前还不能定论,需要深入研究.
因为阔叶红松林适宜于空气湿润、土壤相对肥沃、排
水良好的生境,在我国境内如果北移,则红松将进入
大兴安岭,而大兴安岭地区受地形地貌影响属于寒
温带大陆性季风气候,具有大气干旱、土壤干旱瘠薄
的特点,未来气候变暖情景下,大兴安岭地区的降水
和空气湿度如何尚需进一步的科学推断.尤其是土
壤特性的改变需要相当长的时间.因此,在未来气温
升高、降水不变的情景下,大兴安岭的气候、土壤条
件可能还是不利于阔叶红松林的迁移.当然,在大、
小兴安岭交错地带以及大兴安岭东坡的低洼背风土
壤水热条件相对较好的局部小生境有可能适合红松
生长.综上,未来气温升高、降水不变的气候条件下
可能导致红松分布区缩小.
参考文献
[1] Wei L (卫 林), Wang H⁃M (王辉民), Wang Q⁃D
(王其冬), et al. The influence of climate changes on
Korean pine forest in China. Geographical Research (地
理研究), 1995, 14(3): 17-26 (in Chinese)
[2] Li G⁃Q (李广起), Bai F (白 帆), Sang W⁃G (桑卫
国). Different responses of radial growth to climate war⁃
ming in Pinus koraiensis and Picea jezoensis var. komaro⁃
viiat their upper elevational limits in Changbai Mountain.
Chinese Journal of Plant Ecology (植物生态学报),
2011, 35(5): 500-511 (in Chinese)
[3] Yin L (尹 璐), An N (安 宁), Long L⁃P (龙良
平), et al. Response of carbon isotope composition of
Chinese red pine annual growth ring cellulose to east
China temperature variations. Acta Mineralogica Sinica
(矿物学报), 2005, 25(2): 103-106 (in Chinese)
[4] Wang X⁃C (王晓春), Zhao Y⁃F (赵玉芳). Growth re⁃
lease determination and interpretation of Korean pine
and Koyama spruce in Shengshan National Nature Re⁃
serve, Heilongjiang Province, China. Acta Ecologica
Sinica (生态学报), 2011, 31(5): 1230- 1239 ( in
Chinese)
[5] Wang H (王 辉), Shao X⁃M (邵雪梅), Fang X⁃Q
(方修琦), et al. Responses of Pinus koraiensis tree ring
cell scale parameters to climate elements in Changbai
Mountains. Chinese Journal of Applied Ecology (应用生
态学报), 2011, 22(10): 2643-2652 (in Chinese)
[6] Chen L, Wu SH, Pan T. Variability of climate⁃growth
relationships along an elevation gradient in the Changbai
Mountain, Northeastern China. Trees, 2011, 25: 1133-
1139
[7] Gao L⁃S (高露双), Wang X⁃M (王晓明), Zhao X⁃H
(赵秀海). Response of Pinus koraiensis and Picea
jezoensis var. komaroviito climate in the transition zone of
Changbai Mountain. Chinese Journal of Plant Ecology
(植物生态学报), 2011, 35(1): 27-34 (in Chinese)
[8] Hao Z⁃Q (郝占庆), Dai L⁃M (代力民), He H⁃S (贺
红士), et al. Potential response of major tree species to
climate warming in Changbai Mountain, Northeast Chi⁃
na. Chinese Journal of Applied Ecology (应用生态学
报), 2001, 12(5): 312-319 (in Chinese)
[9] Chen L (陈 力), Yin Y⁃H (尹云鹤), Pan T (潘
韬), et al. Growth of Pinus koraiensis and its response
to temperature change at Changbai Mountain, Northeast
China. Resources Science (资源科学), 2012, 34(11):
2139-2145 (in Chinese)
[10] Yu DP, Liu JQ, Lewis BJ, et al. Spatial variation and
temporal instability in the climate⁃growth relationship of
Korean pine in the Changbai Mountain region of Nor⁃
theast China. Forest Ecology and Management, 2013,
300: 96-105
[11] Chen L (陈 力), Wu S⁃H (吴绍洪), Dai E⁃F (戴
尔阜 ). Analysis of the tree⁃ring width chronologies
of Pinus koraiensis and Larix olgensison Changbai
Mountains, Northeast China. Geographical Research (地
理研究), 2011, 30(6): 1147-1155 (in Chinese)
[12] Wang M (王 淼), Bai S⁃J (白淑菊), Tao D⁃L (陶
大立), et al. Effect of rise in air⁃temperature on tree
ring growth of forest on Changbai Mountain. Chinese
Journal of Applied Ecology (应用生态学报), 1995, 6
(2): 128-132 (in Chinese)
[13] Chen L (陈 列), Gao L⁃S (高露双), Zhang Y (张
赟), et al. Characteristics of tree⁃ring chronology of Pi⁃
nus koraierrsis and its relationship with climate factors on
0531 应 用 生 态 学 报 27卷
the northern slope of Changbai Mountain. Acta Ecologica
Sinica (生态学报), 2013, 33(4): 1285- 1291 ( in
Chinese)
[14] Fritts HC. Tree Rings and Climate. London: Academic
Press, 1976
[15] Cook ER. A Time Series Analysis Approach to the Tree
Ring Standardization. PhD Thesis. Tucson, AZ: Univer⁃
sity of Arizona, 1985
[16] Zhang X⁃L (张先亮), Cui M⁃X (崔明星), Ma Y⁃J
(马艳军), et al. Larix gmelinii tree⁃ring width chrono⁃
logy and its responses to climate change in Kuduer,
Great Xing’ an Mountains. Chinese Journal of Applied
Ecology (应用生态学报), 2010, 21(10): 2501-2507
(in Chinese)
[17] Zhao Z⁃J (赵志江). The Responses of Abies faxoniana
and Picea purpurea to Climate Factors in Subalpine of
Western Sichuan Province, China. PhD Thesis. Beijing:
Beijing Forestry University, 2013 (in Chinese)
[18] Zhou Y⁃L (周以良). Vegetation Geography in Northeast
China. Beijing: Science Press, 1997 (in Chinese)
[19] Guo Y⁃T (郭元涛), Liu Z⁃P (刘忠平), Song G⁃L
(宋桂莲), et al. Climate characteristics of Baishilazi
Nature Reserve and its relationship with forest plants.
Journal of Liaoning Forestry Science & Technology (辽
宁林业科技), 2011(4): 11-12 (in Chinese)
[20] Wang M (王 淼), Ji L⁃Z (姬兰柱), Li Q⁃R (李秋
荣), et al. Effects of soil temperature and moisture on
soil respiration in different forest types in Changbai
Mountain. Chinese Journal of Applied Ecology (应用生
态学报), 2003, 14(8): 1234-1238 (in Chinese)
[21] Hao Z⁃Q (郝占庆), Guo S⁃L (郭水良), Cao T (曹
同). Plant Diversity and Distribution Patterns in Chang⁃
bai Mountains. Shenyang: Liaoning Science and Tech⁃
nology Press, 2002 (in Chinese)
[22] Ma Y⁃X (马玉心), Cai T⁃J (蔡体久), Tan X⁃J (谭
晓京), et al. A research on spermatophyte flora of pri⁃
mary broad⁃leaved and Korean pine forest in Liangshui
Natural Reserve. Journal of Northeast Normal University
(Natural Science) (东北师大学报: 自然科学版),
2007, 39(1): 84-90 (in Chinese)
[23] Liu M (刘 敏), Mao Z⁃J (毛子军), Li Y (厉
悦), et al. Climatic effects on radial growth of Korean
pines with different bark forms in Liangshui Natural Re⁃
serve, Northeast China. Chinese Journal of Applied Eco⁃
logy (应用生态学报), 2014, 25(9): 2511-2520 (in
Chinese)
[24] Zhang H⁃Q (张恒庆), An L⁃J (安利佳), Zu Y⁃G
(祖元刚). RAPD analysis on genetic variation of Pinus
koraiensis in Liangshui National Reserve. Bulletin of Bo⁃
tanical Research (植物研究), 2000, 20(2): 201-206
(in Chinese)
[25] Wang X⁃C (王晓春), Zhao Y⁃F (赵玉芳). Growth re⁃
lease determination and interpretation of Korean pine
and Koyama spruce in Shengshan National Nature Re⁃
serve, Heilongjiang Province, China. Acta Ecologica
Sinica (生态学报), 2011, 31(5): 1230- 1239 ( in
Chinese)
[26] Stokes MA, Smiley TL. An Introduction to Tree⁃ring
Dating. Chicago: University of Chicago Press, 1968
[27] Holmes RL. Dendrochronology Program Library, Version
1992-1. Tucson, AZ: Laboratory of Tree⁃ring Re⁃
search, University of Arizona, 1992
[28] Cook ER, Holmes RL. Users Manual for Program AR⁃
STAN. Tucson, AZ: Laboratory of Tree⁃ring Research,
University of Arizona, 1986
[29] Cook ER, Kairiukstis LA. Methods of Dendrochronolo⁃
gy: Applications in the Environmental Sciences. Dor⁃
drecht, the Netherlands: Kluwer Academic Publishers,
1990
[30] Wu X⁃D (吴祥定). Tree Rings and Climate Change.
Beijing: China Meteorological Press, 1990 ( in Chi⁃
nese)
[31] Gao L⁃S (高露双). The Radial Growth in Relation to
Climate Factors for Typical Species in Changbai Moun⁃
tain. PhD Thesis. Beijing: Beijing Forestry University,
2011 (in Chinese)
[32] Li G⁃Q (李广起), Bai F (白 帆), Sang W⁃G (桑卫
国). Different responses of radial growth to climate war⁃
ming in Pinus koraiensis and Picea jezoensis var. komaro⁃
vii at their upper elevational limits in Changbai Moun⁃
tain, China. Chinese Journal of Plant Ecology (植物生
态学报), 2011, 35(5): 500-511 (in Chinese)
[33] Dai A, Trenberth KE, Qian TT. A global dataset of
Palmer drought severity index for 1870-2002: Relation⁃
ship with soil moisture and effects of surface warming.
Journal of Hydrometeorology, 2004, 5: 1117-1130
[34] Yu D⁃P (于大炮), Wang S⁃Z (王顺忠), Tang L⁃N
(唐立娜), et al. Relationship between tree⁃ring chro⁃
nology of Larix olgensisin Changbai Mountains and the
climatic change. Chinese Journal of Applied Ecology (应
用生态学报), 2005, 16(1): 14-20 (in Chinese)
[35] Ji Y (及 莹). Climate⁃growth Relationships of Korean
Pine in Heilongjiang and Their Potential for Global War⁃
ming. Master Thesis. Harbin: Northeast Forestry Univer⁃
sity, 2010 (in Chinese)
[36] Liang EY, Shao XM, Xu Y. Tree⁃ring evidence of re⁃
cent abnormal warming on the southeast Tibetan Plateau.
Theoretical and Applied Climatology, 2009, 98: 9-18
[37] Dang HS, Jiang MX, Zhang QF, et al. Growth respon⁃
ses of subalpine fir (Abies fargesii) to climate variability
in the Qinling Mountain, China. Forest Ecology and
Management, 2007, 240: 143-150
[38] Frank D, Esper J. Characterization and climate response
patterns of a high⁃elevation, multi⁃species tree⁃ring net⁃
work in the European Alps. Dendrochronologia, 2005,
22: 107-121
[39] Wang T, Ren HB, Ma KP. Climatic signals in tree ring
of Picea schrenkiana along an altitudinal gradient in the
central Tianshan Mountains, northwestern China. Trees,
2005, 19: 736-742
[40] Guo J⁃P (郭建平), Gao S⁃H (高素华), Wang L⁃M
(王连敏), et al. Experimental study on the impacts of
CO2 concentration and soil water stress on nursery stocks
of Pinus koraiensis Sibe. et. Zucc and Picea asperata
15315期 刘 敏等: 不同纬度阔叶红松林红松径向生长对气候因子的响应
Mast. Acta Meteorologica Sinica (气象学报), 2004, 65
(4): 493-497 (in Chinese)
[41] Rolland C. Tree⁃ring and climate relationships for Abies
alba in the internal Alps. Tree⁃Ring Bulletin, 1993, 53:
1-11
[42] Yan X⁃D (延晓冬), Zhao S⁃D (赵士洞), Fu C⁃B
(符淙斌), et al. How should the Xiao Xinggan Mt.
forests change with potential climate change: A simula⁃
tion study. Journal of Natural Resources (自然资源学
报), 1999, 14(4): 372-376 (in Chinese)
[43] Cheng X⁃X (程肖侠), Yan X⁃D (延晓冬). Effects of
climate change on typical forest in the northeast of Chi⁃
na. Acta Ecologica Sinica (生态学报), 2008, 28(2):
534-543 (in Chinese)
[44] Zhou D⁃H (周丹卉), He H⁃S (贺红士), Li X⁃Z (李
秀珍), et al. Potential responses of different stand age
classes to climate changes in the Xiaoxinganling
Mountains, Northeastern China. Journal of Beijing
Forestry University (北京林业大学学报), 2007, 29
(4): 110-117 (in Chinese)
[45] Guo Q⁃S (郭泉水), Yan H (阎 洪), Xu D⁃Y (徐德
应), et al. Effects of climate changes on geographical
distribution of Pinus koraiensis in China. Acta Ecologica
Sinica (生态学报), 1998, 18(5): 484-488 ( in Chi⁃
nese)
作者简介 刘 敏,女,1979 年生,副教授,博士研究生. 主
要从事植物生态学和树木年轮生态学研究. E⁃mail: liumi⁃
naa168@ 126.com
责任编辑 孙 菊
刘敏, 毛子军, 厉悦, 等. 不同纬度阔叶红松林红松径向生长对气候因子的响应. 应用生态学报, 2016, 27(5): 1341-1352
Liu M, Mao Z⁃J, Li Y, et al. Response of radial growth of Pinus koraiensis in broad⁃leaved Korean pine forests with different latitudes
to climatical factors. Chinese Journal of Applied Ecology, 2016, 27(5): 1341-1352 (in Chinese)
2531 应 用 生 态 学 报 27卷