桂林三月岭马尾松树轮宽度对气候变化的响应

段丙闯，蔡秋芳

桂林三月岭马尾松树轮宽度对气候变化的响应

段丙闯1,2，蔡秋芳1

1.中国科学院地球环境研究所 黄土与第四纪地质国家重点实验室，西安 710061
2. 中国科学院大学，北京 100049

利用采自广西桂林三月岭湘桂古道的马尾松树轮样本，建立研究区的标准树轮宽度年表。树轮气候响应分析发现：马尾松径向生长与当年6 — 9月平均温度（r = −0.51，p<0.01）和上年2 — 11月平均温度（r = − 0.53，p<0.01）显著负相关，与上年10月降水量显著正相关，与当年6月、7月、8月的平均相对湿度和帕尔默干旱指数（PDSI）显著正相关。说明生长季当年6 — 9月伴随高温产生的干旱胁迫抑制了马尾松的生长，上年气候对当年树木生长存在滞后影响。空间分析显示，研究区树轮年表与我国大范围季风区、西北干旱半干旱地区上年2 — 11月平均温度变化相关关系显著，与中南半岛北部地区当年6 — 9月平均温度显著负相关。1976 — 1977年前后树轮年表与海表温度空间相关型的变化，表明研究区夏季的温度气候要素与大尺度的气候驱动太平洋年代际震荡（PDO）有一定联系。

桂林；马尾松（Pinus massoniana）；中国东部亚热带；树轮宽度；气候响应；温度

树木年轮以其定年准确、连续性强、分辨率高且样本地域分布广等特点成为全球气候变化研究的主要代用指标之一，在国际地圈生物圈计划（IGBP）和过去全球变化（PAGES）研究中扮演着重要角色。近几十年来，我国学者开展了大量的树木年轮气候学研究，其工作主要集中在北方及西北的干旱半干旱气候敏感区（邵雪梅等，2004；刘禹等，2006，2009，2011，2013；蔡秋芳等，2012；蔡秋芳和刘禹，2013a，2015；Cai et al，2015；Li et al，2015），树轮宽度响应模式的相关研究也已十分成熟与深入，并借此进行了大量气候重建工作（勾晓华等，2001；方克艳等，2008；朱海峰等，2008）。然而在深受亚洲夏季风影响的热带、亚热带地区，由于气候湿润、树木与气候的响应关系相对复杂等因素，树木生长的限制因子具有更多的不确定性，使得树轮定年工作较为困难，相关的树轮气候研究在这一区域发展较为缓慢（Cai and Liu，2016；Zheng et al，2016）。

近十多年来，陆续有学者在中国东部亚热带地区进行了数百年来温度指标的重建并取得了很多有价值的历史气候信息（Shi et al，2010；Chen et al，2012；Duan et al，2013；Chen et al，2015）。但想要深入探寻中国南方广大地区历史气候变化规律，仍需要开展大量的树轮气候学基础性研究。本文选取广西桂林地区的马尾松（Pinus massoniana）作为研究对象，运用树木年轮学研究方法，建立马尾松树轮宽度年表，探讨该地区马尾松树木年轮宽度变化与气候因子的关系，为深入开展中国亚热带地区相关树轮重建工作提供基础性研究支持。

1 研究区概况

本组树芯于2015年9月采自广西省桂林市的三月岭湘桂古道，编号为XGGD。采样点位于25°22′ — 25°23′N，110°31′ — 110°32′E（图1），海拔约480 m。研究区属亚热带季风气候，光照充足，热量丰富，雨量充沛，无霜期长，具有夏长冬短，四季分明且雨热基本同季的特点。年平均气温19.3℃。7月最热，月平均气温为28℃；1月最冷，月平均气温7.9℃（图2）。年平均无霜期309 d，年平均降雨量1949.5 mm。平均蒸发量1490 — 1905 mm。年平均相对湿度为73% — 79%，年平均日照时数为1670 h（况雪源等，2007）。

图1 树轮采样点与气象站分布图Fig.1 Map showing the tree-ring sampling sites and the meteorological station

图2 桂林、融安、道县气象站平均月均气温和月降水量年内分配（1959 — 2014年）Fig.2 Mean monthly temperature and total precipitation at the Guilin station, Rong’an station, Daoxian station for the period 1959 to 2014

采样点植被种类丰富，马尾松为主要优势乔木，林相整齐，林下灌木茂盛，部分马尾松被采过松脂，树干挺拔，高~30 m，胸径50 — 110 cm。马尾松树群呈条带状分布在三月岭湘桂古道两侧。研究区属红壤土带，土壤类型为风化程度较高的红壤土。

马尾松广泛分布于中国的亚热带地区。它能够在贫瘠、恶劣的自然条件下良好生长，因而成为中国亚热带地区绿化造林和利用荒山废地恢复植被的主要先锋树种。马尾松也是南方开展区域气候重建的重要树种，但仍需进一步探讨该树种在本研究区进行气候重建的潜力与前景。

2 数据与方法

2.1 样本采集与年表建立

按照树木年轮气候学的采样方法选取生长状况良好的树体，使用口径为5.15 mm的生长锥在呈对角或90º角的两个方向上在胸高部位钻取树轮样芯，依次编号并放置在纸样管内。在采集马尾松树芯样本过程中，尽量钻取目测树龄比较大的树木。一般每棵树钻取两根样芯，个别树钻取第三芯用于稳定同位素分析。最终在采样点钻取42棵树90根样芯。

将样芯放置于干燥处自然风干，然后将其粘贴在特制的木槽内固定，进行磨光处理（Cook and Kairiukstis，1990）后，进行目测定年，然后在精度为0.01 mm的LINTAB树轮宽度测量仪上逐年测量树轮宽度。采用COFECHA程序进行交叉定年的质量控制（Fritts，1976；Holmes，1983），去掉同大多数序列变化趋势相异的序列。采用负指数函数或直线拟合去除马尾松序列中由于树龄导致的幼龄效应（Fritts，1976）。用原始测量序列除以拟合序列即得到标准化树轮宽度指数序列。

利用ARSTAN程序进行树轮宽度年表的制作（Cook，1985），最终得到了马尾松树轮宽度标准化年表（STD）（图3）。考虑到样本年代长短不一，采用子样本信号强度（sub-sample signal strength，SSS）（Wigley et al，1984）来确定可靠年表的起始时间。本文采用SSS>0.85的标准，并对此时段年表的统计特征进行分析。所分析的统计特征有：标准差（SD）、所有序列间相关系数（r1）、树间相关系数（r2）、树内相关系数（r3）、信噪比（SNR）、第一主成分所解释方差量（PC1）、平均敏感度（MS）、一阶自相关系数（AR1）、样本总体代表性（EPS）（表1）。

图3 树轮宽度STD年表及相应时段样本量Fig.3 Standardized tree-ring width chronology and sample depth

表1 桂林三月岭马尾松年轮标准年表统计（1942 — 2014）Tab.1 Statistics of the STD chronologies for Pinus massoniana (1942 — 2014)

2.2 气候资料

本文选取距采样点较近的广西省桂林市、广西融安以及湖南道县三个气象站点的数据。其记录时段分别为1951 — 2014年、1957 — 2014年和1959 — 2014年。由于广西北部地处南岭山地，地理环境复杂，为了尽量减小局地气候对数据分析的影响，故选取三地观测记录的等权重平均值与树轮宽度进行相关分析（1959 — 2014年，其中1959年数据只有8月及之后数据，故舍去不用）。所选用的气候要素包括月平均温度、月平均最高（低）温度、月降水量和月相对湿度数据。此外，美国国家大气研究中心（NCAR）Dr. Aiguo Dai的帕尔默干旱指数（Palmer Drought Severity Index，PDSI）（108.75° — 111.25°E，23.75° — 26.25°N，1942 — 2005）也被用来进行分析（http://climexp. knmi.nl）。

3 结果与讨论

3.1 树轮年表

共有40棵树的63根树芯通过COFECHA程序的质量检验。利用这63个序列建立了树轮宽度STD年表，时段为1939 — 2014年，可靠年表区间为1942 — 2014年。年表的各个统计值如表1所示。样本原始序列的序列间相关系数为0.59，标准化年表的r1、r2和r3也均较高，表明各序列年轮宽窄变化比较同步，可能受一致的外界气候环境变化的影响。MS是度量相邻年轮宽度变化情况的一个指标，平均敏感度越大，表明对气候变化越敏感，样本中保留的气候信息相对就越多（蔡秋芳和刘禹，2013b）。本研究中MS值为0.22，符合亚热带地区树轮气候学研究的标准（Speer，2010）。马尾松标准年表AR1为0.43，说明上年气候对下年树木生长滞后效应明显，气候响应分析时需将上年的气候因子考虑进去。

3.2 气候响应分析

选取上一年和当年共24个月的气候因子与建立的树轮年表进行Pearson相关分析，确定影响研究区树木径向生长的主要气候因子。考虑到本区气候特点和马尾松生长习性，又将气候要素进行了不同月份的组合，进一步分析马尾松径向生长与气候要素不同月份组合之间的关系。马尾松树轮宽度年表与气候要素的相关分析结果见图4。

如图4a所示，树轮宽度与上年及当年大部分月份的温度呈负相关关系，且与上年2月、3月、7 — 10月及当年7 — 9月的单月温度呈不同程度的显著负相关。如果按照月平均温度高于10 — 10.5℃马尾松进入生长期（封晓辉，2011），则研究区2 — 11月马尾松均处于生长期。然而，当年2 — 5月平均温度与树轮宽度指数均无显著相关，是因为马尾松自2月进入抽芽期至4月的径向生长期，温度适宜，水分充足，对气候变化的响应较差；而5月该地区平均温度虽然较高，但此时出现的降水高峰期（图2）避免了高温可能带来的干旱胁迫。树轮宽度年表与当年7月、8月、9月的平均温度呈显著负相关，与当年6 — 9月4个月份组合的平均温度相关系数为−0.51（p< 0.01），考虑到是因为该时段持续的高温造成了土壤水分蒸发减少，土壤和马尾松自身的蒸腾作用也加强，每月逐渐减少的降水已无法缓解高温带来的生理干旱，从而限制了马尾松的生长（Alley，1984；卫捷等，2003）。此外，树轮宽度年表与上年2 — 11月组合的月平均温度相关系数最高，相关系数为−0.53（p< 0.01）。指示上年生长季的气温对马尾松的生长有显著的滞后效应。Cai and Liu（2016）在江西三清山黄山松的研究中也发现了类似的滞后效应，即黄山松树轮宽度对上年3 — 10月气候要素变化的响应十分显著（r= 0.75，p < 0.01），高出了对当年气候要素的响应。可见，此处树轮宽度与气候要素的响应分析结果是合理的。

桂林三月岭马尾松树轮宽度除与上年10月降水量显著正相关外（r = 0.32，n= 64，p< 0.01），与其他月降水量均呈较低相关（图4b），表明降水量对研究区马尾松生长的影响表现在季风区的干旱季节。10月份属于马尾松的生长季晚期，上年10月份较多的降水有利于树体存储较多的营养物质或改良土壤水热状况，使其在下一年生长季到来时更好地生长。

此外，将树轮年表与相对湿度以及PDSI也进行了相关分析。相对湿度指空气中水汽压与饱和水汽压的百分比（ 段若溪和姜会飞，2002）。对于植物-土壤-大气统一体而言，相对湿度和土壤湿度直接相关。温度较高的夏季，强烈的蒸腾作用使得马尾松的生长受到干旱胁迫的影响。相对湿度较高的环境能有效地降低叶片内外的蒸汽压差，有效地减弱其蒸腾作用，利于植物体水分的保持。树轮年表与当年6月、7月、8月份平均相对湿度显著正相关，与6 — 7月份相关系数最高（r = 0.55，p < 0.01），表明此时段较高的相对湿度，能够有效促进植物的生长，缓解高温带来的生理干旱。

帕尔默干旱指数是一个综合考虑了当时的水分条件、前期水分状况及持续时间的干旱指标，可作为水分供需累积效应的一个近似度量。相关分析发现，树轮宽度指数与当年6 — 11月的单月PDSI指数正相关显著，当年7月份单月最高（r= 0.44，n= 64，p< 0.01）（图4d），说明了水热供需对马尾松生长的限制作用主要发生在夏季及夏季后的干旱季节；与上年7月至当年1月的单月PDSI指数均呈显著正相关，这与陈峰等（2015）在闽中北柳杉树轮宽度序列的研究结果大体一致。陈峰等（2015）在其文中指出该宽度年表和上年7月至当年2月的PDSI有显著的正相关关系。

图4 桂林三月岭马尾松树轮宽度标准化年表与（a）月平均、最高、最低温度、（b）降水量、（c）相对湿度及（d）PDSI的相关系数（横坐标的P代表上年，C代表当年）Fig.4 Correlation coeffi cients for tree-ring chronology with (a) monthly temperatures, (b) precipitation, (c) relative humidity, (d) PDSI (P and C represent previous and current year, respectively)

3.3 空间相关分析

研究区气象站2 — 11月平均温度（图5a）及6 — 9月平均温度（图5b）与同期CRU温度格点数据的空间相关分析结果如图5，表明所使用气象站2 — 11月及6 — 9月的平均温度数据与大范围地区的温度变化同步性较好，具有广泛的空间代表性。树轮宽度指数与CRU格点数据空间相关分析表明，桂林三月岭树轮宽度指数与我国较大范围地区上年2 — 11月平均气温存在显著负相关（1960 — 2014）（图5c），主要包括东南沿海地区、西北干旱区以及青藏高原南缘一带；树轮宽度指数与当年6 — 9月平均气温的空间相关结果显示，相关显著的地区包括东南沿海、台湾地区、中南半岛北部等地区。综上可知，研究区马尾松对周边亚热带地区的夏季温度具有较强的指示意义，而对广大东亚季风区及西北干旱区的温度变化也有良好的相关关系，其指示意义存在于上年生长季平均温度对当年树木生长的滞后作用中。

陶诗言和张庆云（1998）指出，ENSO事件中，El Niño年印度夏季风会相对减弱，La Niña年则使得印度夏季风增强。为进一步证实研究区马尾松对夏季主要受西南季风影响的亚热带地区的温度确实存在较强的指示意义，本文进一步将树轮年表与南方涛动指数（SOI）（http://climexp. knmi.nl）及海温进行相关分析。南方涛动指数是应用数理统计的方法将塔希提站与达尔文气象观测站的海平面气压差值进行处理后得到的一个用于衡量南方涛动强弱的指数，很好地反映了El Niño现象的活跃程度（杨梅学和姚檀栋，1999）。

将1942 — 2014年的树轮年表与同年6 — 9月SOI进行相关分析，相关系数未能达到显著水平。但有学者指出，1976 — 1977年气候突变，太平洋年代际震荡（PDO）由冷位相进入暖位相年后，El Niño发生更加频繁（Trenberth et al，2002；丁一汇和孙颖，2006；姜超等，2011）。因此将树轮年表与SOI指数做分段相关，时段分别为1942 — 1976年和1977 — 2014年，发现在前一时段中两者并无相关关系，后一时段年表与SOI相关系数显著（r = 0.44，p < 0.01）。进一步将树轮年表与相关海域的海温数据进行分段空间相关，海温数据选用HadlSST1格点数据（http://climexp. knmi.nl），结果如图6。1942 — 1976年树轮年表与印度洋及孟加拉湾海域及SST显著相关，而在1977 — 2014时段中并无相关关系，后一时段中树轮年表与Niño 3.4及临近海域温度显示出正相关关系。产生该结果的原因可能如下：1942 — 1976年，西南季风相对较强，研究区受西南季风携带的水热影响较大，因而树木的生长与季风水热源区下垫面的海温状况相关显著；1977 — 2014年，El Niño活跃期，印度季风相对较弱的时期，西南季风水热源区对树木的生长的影响有限，受东亚夏季风水热状况影响较大，表现为与El Niño 3.4区的海温呈正相关。

图5 （a）气象站2—11月平均温度观测记录（1959 — 2013）与同期CRU TS3.23格点数据的空间相关分析，（b）气象站6—9月平均温度观测记录（1960 — 2014）与同期CRU TS3.23格点数据的空间相关分析，（c）树轮宽度年表（1960 — 2014）与上年2 — 11月CRU TS3.23格点数据的空间相关分析，（d）树轮宽度年表（1960 — 2014）与当年6 — 9月CRU TS3.23格点数据的空间相关分析（图中黑点为采样点）Fig.5 Spatial Pearson correlation plots for February—November monthly mean temperature (a) and June—September monthly mean temperature (b) with CRU TS3.23 grid data of the corresponding period, Spatial Pearson correlation plots for chronology with (c) February—November monthly mean temperature, (d) June—September monthly mean temperature

值得一提的是，在XGGD的树轮年表中，2007年的树轮宽度属正常偏高值，在2008年急剧下降为记录时段的年轮宽度最低值，其后树轮宽度也处于1939年以来最低值阶段。考虑到是受2008年的雪灾影响。2008年1月10日起发生在中国南方大范围的雪灾冻害事件也波及到广西地区（李晓靖等，2011）。取样时据当地人介绍，2008年雪灾严重影响到了当地树木生长，许多高大的马尾松被雪压垮，主要受害情况有断梢断杆等，且半数左右马尾松因灾死亡。结合2008年之后树木年轮宽度一直处于极低值这一结果，说明了极端气候事件会记录在树轮中，且其滞后影响可能会持续很多年。李金建等（2014）在松潘的研究中，其年表及重建结果也与雪灾发生的年份很好地对应了起来。由此可见，研究区马尾松受当年及上年生长季的水热状况的影响，但是生长季前期冬季的极端冷事件同样在树木年轮中留下了印记。以往的研究主要根据相关系数显著的响应时段季节时段进行气候重建，检测不同年份重建季节的极端气候事件。但是这些研究并未考虑响应季节以外的突发极端气候的影响，由此可能导致重建工作中存在对极端气候事件的误读；考虑到突发气候事件而忽略其滞后影响也会使得极端或异常年份之后时段的重建工作与历史事实有所偏差。

图6 树轮年表与HadlSST1格点数据不同时段的空间相关：（a）1942 — 1976年；（b）1977 — 2014年Fig.6 Spatial Pearson correlation plots for grid data HadlSST1 with chronology in (a)1942 — 1976 and (b)1977 — 2014

4 结论

利用采自桂林三月岭地区马尾松树轮样本建立了1939 — 2014年的树轮宽度年表，其可靠时段为1942 — 2014年。通过对树轮宽度STD年表与气候因子的相关分析发现，桂林三月岭马尾松树轮年表主要受当年6 — 9月和上年2 — 11月平均气温的影响（负相关），树轮宽度仅与上年10月的降水量显著负相关。此外，树轮宽度与6月、7月、8月相对湿度及PDSI均显著正相关，与6 — 7月的相对湿度及7月PDSI相关最高，表明此时段内温度变化引起的水分亏缺是限制研究区马尾松径向生长的主要因素。空间分析发现该树轮宽度年表对夏季主要受西南季风控制的较大范围地区的夏季温度具有指示意义，而对东亚广大季风区温度的指示意义体现在温度变化对树木生长的滞后效应中。并对全年树木生长状况做出总结：（1）桂林三月岭地区在相对偏暖的年份，极可能2月便已进入水分充足、温度适宜的生长期，并在5月及5月之前保持着良好的生长状态；（2）6 — 9月，研究区马尾松受到高温带来的干旱胁迫影响，温度成为影响研究区马尾松生长的主要气候因子；（3）10 — 11月树木依旧处于生长季，此时树木的生长受水热因子的综合作用影响较大；（4）研究区上年气候条件（2 — 11月平均温度及10月降水）对次年马尾松的径向生长有明显的滞后影响。

本文所建的树轮年表对研究区及周边较大范围的6 — 9月平均温度及上年2 — 11月平均温度均有较高相关，但相关性尚未达到气候重建的要求，且本文所建的树轮宽度年表长度较短，未来有必要在周边区域寻找更多序列更长的树木年轮资料，以获取研究区较长尺度上的温度信号。

致谢：感谢王艳超野外协助采样；感谢沈宝印在样品处理过程中给予的帮助。

蔡秋芳, 刘 禹, 王艳超. 2012. 陕西太白山树轮气候学研究[J]. 地球环境学报, 3(3): 702 – 708. [Cai Q F, Liu Y, Wang Y C. 2012. Dendroclimatic investigation of Chinese Pine in Taibai Mountain, Shaanxi Province [J]. Journal of Earth Environment, 3(3): 702 – 708.]

蔡秋芳, 刘 禹. 2013a. 山西五鹿山油松树轮宽度年表的建立及过去百余年5—6月平均气温变化[J]. 第四纪研究, 33(3): 511 – 517. [Cai Q F, Liu Y. 2013. The development of a tree-ring width chronology and the May—June mean temperature variability in Wulu Mountain, Shanxi Province of north-center China [J]. Quaternary Sciences,33(3): 511 – 517.]

蔡秋芳, 刘 禹. 2013b. 湖北麻城马尾松树轮宽度对气候的响应及1879年以来6—9月平均最高气温重建[J]. 科学通报, 58(1): 169 – 177. [Cai Q F, Liu Y. 2013. The June—September maximum mean temperature reconstruction from Masson pine tree-rings in Macheng, southeast China since 1879 AD [J]. Chinese Science Bulletin, 58(1): 169 – 177.]

蔡秋芳, 刘 禹. 2015. 山西中条山白皮松和华山松径向生长对气候变化的响应及气候意义[J]. 地球环境学报, 6(4): 208 – 218. [Cai Q F, Liu Y. 2015. Climatic responses of Pinus bungeana Zucc. and Pinus armandii Franch. in Mt Zhongtiao, Shanxi Province, and May—July temperature variations since 1920 [J]. Journal of Earth Environment, 33(6): 690 – 695.]

陈 峰, 袁玉江, 喻树龙. 2015. 闽中北柳杉树轮指示的气候信号与季风区不同地域干湿变化关系[J]. 山地学报, 33(6): 690 – 695. [Chen F, Yuan Y J, Yu S L. 2015. Drought signals in the tree-ring width record of Cedar (Cryptomeria fortunei) trees from north central Fujian: Linkages to the Monsoonal Regions [J]. Mountain Research, 33(6): 690 – 695.]

丁一汇, 孙 颖. 2006. 国际气候变化研究新进展[J]. 气候变化研究进展, 2(4): 161 – 167. [Ding Y H, Sun Y. 2006. Recent advances in climate change science [J]. Advances in Climate Change Research, 2(4): 161 – 167.]

段若溪, 姜会飞. 2002. 农业气象学[M]. 北京: 气象出版社: 1 – 275. [Duan R X, Jiang H F. 2002. Agrometeorology [M]. Beijing: China Meteorological Press: 1 – 275.]

方克艳, 勾晓华, 陈发虎, 等. 2008. 树轮生态学研究进展[J].

冰川冻土, 30(5): 825 – 834. [Fang K Y, Gou X H, Chen F H, et al. 2008. The advance of dendroecology [J]. Journal of Glaciology and Geocryology, 30(5): 825 – 834.]

封晓辉. 2011. 气候变化对中国马尾松径向生长的影响[D].北京: 中国林业科学研究院. [Feng X H. 2011. The effect of climate change on radial growth of Pinus massoniana in China [D]. Beijing: Chinese Academy of Forestry.]

勾晓华, 陈发虎, 王亚军, 等. 2001. 利用树轮宽度重建近280年来祁连山东部地区的春季降水[J]. 冰川冻土, 23(3): 292 – 296. [Gou X H, Chen F H, Wang Y J, et al. 2001. Spring precipitation reconstructed in eastern Qilian Mountains during the last 280 a by tree-ring widths [J]. Journal of Glaciology and Geocryology, 23(3): 292 – 296.]

姜 超, 徐永福, 季劲钧, 等. 2011. ENSO年代际变化对全球陆地生态系统碳通量的影响 [J]. 地学前缘, 18(6): 107 – 116. [Jiang C, Xu Y F, Ji J J, et al. 2011. Infl uences of the decadal variation of ENSO on the carbon fl ux in the terrestrial ecosystems [J]. Earth Science Frontiers, 18(6): 107 – 116.]

况雪源, 苏 志, 涂方旭. 2007. 广西气候区划[J]. 广西科学, 14(3): 278 – 283. [Kuang X Y, Su Z, Tu F X. 2007. Climate regionalization of Guangxi [J]. Guangxi Sciences, 14(3): 278 – 283.]

李金建, 邵雪梅, 李媛媛, 等. 2014. 树轮宽度记录的松潘地区年平均气温变化[J]. 科学通报, 59: 1446 – 1458. [Li J J, Shao X M, Li Y Y, et al. 2014. Annual temperature recorded in tree-ring from Songpan region [J]. Chinese Science Bulletin, 59: 1446 – 1458.]

李晓靖, 洪奕丰, 王宗星, 等. 2011. 2008年南方冰雪灾害对森林生态系统影响研究[J]. 安徽农业科学, 39(36): 22374 – 22377. [Li X J, Hong Y F, Wang Z X, et al. 2011. Effects of ice and snow storm on forest eco-system in the southern of China in 2008 [J]. Journal of Anhui Agricultural Sciences, 39(36): 22374 – 22377.]

刘 禹, 安芷生, Linderholm H W, 等. 2009. 青藏高原中东部过去2485年以来温度变化的树轮记录[J]. 中国科学D辑:地球科学, 39(2): 166 – 176. [Liu Y, An Z S, Linderholm H W, et al. 2009. Annual temperatures during the last 2485 years in the mid-eastern Tibetan Plateau inferred from tree rings [J]. Science China Earth Sciences, 52(3): 348 – 359.]

刘 禹, 安芷生, 马海州, 等. 2006. 青海都兰地区公元850年以来树轮记录的降水变化及其与北半球气温的联系[J]. 中国科学: D辑, 36(5): 461 – 471. [Liu Y, An Z S, Ma H Z, et al. 2006. Precipitation variation in the northeastern Tibetan Plateau recorded by the treerings since 850 AD and its relevance to the Northern Hemisphere temperature [J]. Science in China (Series D), 36 (5): 461 – 471.]

刘 禹, 蔡秋芳, 宋慧明, 等. 2011. 青藏高原中东部2485年温度变化幅度、速率、周期、原因及未来趋势[J].科学通报, 56(28/29): 2986 – 2994. [Liu Y, Cai Q F, Song H M, et al. 2011. Amplitudes, rates, periodicities and causes of temperature variations in the past 2485 years and future trends over the central-eastern Tibetan Plateau [J]. Chinese Science Bulletin, 56(28/29): 2986 – 2994.]

刘 禹, 蔡秋芳, 宋慧明. 2013. 关于青藏高原2485年温度的季节和空间代表性问题[J]. 第四纪研究, 33(1): 108 – 114. [Liu Y, Cai Q F, Song H M. 2013. Seasonal and spatial representativeness of the tree-ring based2485-year annual mean temperature reconstruction in the northeastern Tibetan Plateau [J]. Quaternary Sciences, 33(1): 108 – 114.]

邵雪梅, 黄 磊, 刘洪滨, 等. 2004. 树轮记录的青海德令哈地区千年降水变化[J]. 中国科学: D辑, 34(2): 145 – 153. [Shao X M, Huang L, Liu H B, et al. 2004. Past precipitation variations during the last 1000 a in Delingha area, Qinghai province, as recorded in tree-ring index [J]. Science in China (Series D), 34(2): 145 – 153.]

陶诗言, 张庆云. 1998. 亚洲冬夏季风对ENSO事件的响应[J]. 大气科学, 22(4): 399 – 407. [Tao S Y, Zhang Q Y. 1998. Response of the Asian winter and summer monsoon to ENSO events [J]. Scientia Atmospherica Sinica, 22(4): 399 – 407.]

卫 捷, 陶诗言, 张庆云. 2003. Palmer干旱指数在华北干旱分析中的应用[J]. 地理学报, 58(S1): 91 – 99. [Wei J, Tao S Y, Zhang Q Y. 2003. Analysis of drought in northern China based on the Palmer Severity Drought Index [J]. Acta Geographica Sinica, 58(S1): 91 – 99.]

杨梅学, 姚檀栋. 1999. 季风现象与ENSO事件研究 [J].

环境科学进展, 7(4): 96 – 101. [Yang M X, Yao T D. 1999. The studies on the monsoon Phenomena and enso events [J]. Advances in Environmental Science, 7(4): 96 – 101.]

朱海峰, 郑永宏, 邵雪梅, 等. 2008. 树木年轮记录的青海乌兰地区近千年温度变化[J]. 科学通报, 53(15): 1835 – 1841. [Zhu H F, Zheng Y H, Shao X M, et al. 2008. Millennial temperature reconstruction based on tree ring widths of Qilian juniper from Wulan, Qinghai Province, China [J]. Chinese Science Bulletin, 53(24): 3914 – 3920.] Alley W M. 1984. The Palmer drought severity index: limitations and assumptions [J]. Journal of Climate and Applied Meteorology, 23(7): 1100 – 1109.

Cai Q F, Liu Y, Liu H, et al. 2015. Reconstruction of drought variability in North China and its association with sea surface temperature in the joining area of Asia and Indian-Pacific Ocean [J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 417: 554 – 560.

Cai Q F, Liu Y. 2016. Two centuries temperature variations over subtropical southeast China inferred from Pinus taiwanensis Hayata tree-ring width [J]. Climate Dynamics, 48(5): 1813 – 1825.

Chen F, Yuan Y J, Wei W S, et al. 2012. Reconstructed temperature for Yong’an, Fujian, Southeast China: linkages to the Pacific Ocean climate variability [J]. Global and Planetary Change, 86: 11 – 19.

Chen F, Yuan Y J, Wei W S, et al. 2015. Tree-ring response of subtropical tree species in southeast China on regional climate and sea-surface temperature variations [J]. Trees, 29(1): 17 – 24.

Cook E R. 1985. A time series analysis approach to tree ring standardization (dendrochronology, forestry, dendroclimatology, autoregressive process) [D]. Arizona: The University of Arizona: 1 – 171.

Cook E T, Kairiukstis L A. 1990. Methods of dendrochronology [M]. Dordrecht: Kluwer Academic Publishers: 1 – 394.

Duan J P, Zhang Q B, Lyu L X. 2013. Increased variability in cold-season temperature since the 1930s in subtropical China [J]. Journal of Climate, 26(13): 4749 – 4757.

Fritts H C. 1976. Tree rings and climate [M]. New York: Academic Press: 1 – 579.

Holmes R L. 1983. Computer-assisted quality control in treering dating and measurement [J]. Tree-Ring Bulletin, 43(1): 69 – 78.

Li Q, Liu Y, Nakatsuka T, et al. 2015. The 225-year precipitation variability inferred from tree-ring records in Shanxi Province, the North China, and its teleconnection with Indian summer monsoon [J]. Global and Planetary Change, 132: 11 – 19.

Shi J F, Cook E R, Lu H Y, et al. 2010. Tree-ring based winter temperature reconstruction for the lower reaches of the Yangtze River in southeast China [J]. Climate Research, 41(2): 169 – 175.

Speer J H. 2010. Fundamentals of tree-ring research [M]. Arizona: University of Arizona Press: 1 – 509.

Trenberth K E, Stepaniak D P, Caron J M. 2002. Interannual variations in the atmospheric heat budget [J]. Journal of Geophysical Research: Atmospheres, 107. DOI: 10.1029/2000JD000297.

Wigley T M L, Briffa K R, Jones P D. 1984. On the average value of correlated time series, with applications in dendroclimatology and hydrometeorology [J]. Journal of climate and Applied Meteorology, 23(2): 201 – 213.

Zheng Y H, Shao X M, Lu F, et al. 2016. February—May temperature reconstruction based on tree-ring widths of Abies fargesii from the Shennongjia area in Central China [J]. International Journal of Biometeorology, 60: 1175 – 1181.

Responses of tree-ring width of Pinus massoniana to climate change in Guilin, Guangxi Province

DUAN Bingchuang1,2, CAI Qiufang1
1. State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi’an 710061, China
2. University of Chinese Academy of Sciences, Beijing 100049, China

Background, aim, and scope It is essential for researchers to know more about past climate variations from the study of climatic proxy data, such as ice core, loess, tree-ring and stalagmite. Tree-ring, due to its accurate dating, continuous, high resolution and sensitive to climate change, has been widely used to reconstruct the past climatic index. Many paleoclimate series has been rebuilt based on tree-ring width by researchers. In China, tree-ring studies in the subtropical regions are much rarer than extratropical regions, which are mainly caused by the complex relationship between environmental factors and the tree-ring growth in subtropical area and the diffi culty in cross dating. Our study aimed to reveal the responding mechanism of tree-rings to climatefactors. Materials and methods In our study, the standardized tree-ring width chronology was developed based on 90 tree-ring cores from 42 healthy Pinus massoniana, collected from Guilin, in subtropical south China (25°22′ — 25°23′N, 110°31′ — 110°32′E, 480 m a.s.l) during September, 2015. All the tree-ring cores were dealt with standard dendrochronological method. The cores were visually dated and measured using the LINTAB measurement machine with resolution of 0.01 mm and cross dated by the COFECHA program. Then we use the series that successfully passed the COFECHA program to develop the standardized tree-ring chronology via the ARSTAN program. Considering the complex environment of sampling site, we use the mean meteorological data of the nearest three weather stations (Guilin station, Rong’an station and Daoxian station, records from 1959 to 2014) to reduce the possible infl uence of micro-climate. To determine whether the climate had impact on the radial growth of Pinus massoniana and which climate factor limited tree growth most, we used Pearson correlation analyses to calculate the correlation coefficients between ring-width STD chronology and climate factors. Global Monthly Dai Palmer Drought Severity Index (PDSI) (108.75° — 111.25°E, 23.75° — 26.25°N, 1942 — 2005) were also analysed against the STD chronology. Spatial correlation patterns of the observed temperature and the STD chronology and the gridded data (land temperatures and SST) were calculated via the website http://climexp.knmi.nl. Results Growth-climate response analyses showed that the growth of Pinus massoniana negatively correlated with the monthly mean temperature variability of current June — September (r = −0.51, p<0.01) and previous February — November (r = −0.53, p<0.01), positively correlated with the previous October precipitation (r = 0.32, p<0.05). Tree rings also signifi cantly correlated with the monthly relative humidity and PDSI of current June, July and August. Result of spatial correlation showed that the XGGD tree-ring STD chronology significantly correlated with the temperatures in the monsoon fringe area and north Indochina Peninsula area. Spatial correlation pattern between tree-ring STD chronology and SST also showed a shift before and after 1977. Discussion In the study, no signifi cant correlation was found between the meteorological data and chronology in the early stage of tree growth. The meteorological data showed an abundant precipitation during the early stage of tree growth (April — June) and continuous high temperature from May to September. It is not possible that growth rate of trees are limited by precipitation. But when it comes to summer, higher temperature and less precipitation may lead to stronger evapotranspiration which can cause severe water stress for plants and thus limit the growth of trees. Negative correlations between tree-ring STD chronology and temperature indicate that in study area the growth of Pinus massoniana was greatly limited by severe water stress caused by high temperature in current June — September. And there exists a quite significant lag-effect of temperature on the growth of Pinus massoniana, which has been widely reported for different species and was particularly signifi cant in southeast China. Moreover, tree-ring STD chronology showed different spatial correlation pattern while correlating with the monthly mean temperature of current June — September and previous February — November. Tree growth in the study area has a well correlationship with current June — September’s temperature of north Indochina Peninsula area, which may mean June — September’s hydrothermal condition of study area were directly infl uenced by the southwest monsoon. Moreover, extreme climate events such as the snow disaster in 2008 also impacted on the growth of tree rings in and after 2008. The error caused by extreme climate events should be taken into account when dealing with climate reconstruction. Conclusions Result showed that the ringwidth STD chronology from the study area had the potential to obtain climate signal especially the summer’s temperature. In the study area, the tree-ring width is sensitive to the temperature of current June — September and previous February — November. Moreover, extreme climate events also have long-term effects on tree growth, which may give rise to some undesirable infl uence on climate reconstruction. Change of the spatial correlation patterns for SST with chronology before and after 1977 indicated that there was a certain link between PDO and the summer’s temperature in the research area. Recommendations and perspectives Tree-ring research in southChina is relatively rare. More tree-ring studies should be carried on in this subtropical regions, which will help us to know more about the climatic response of tree growth in subtropical regions and untie the shackles that limit our fully understanding of critical climate issues.

Guilin; Pinus massoniana; subtropical southeast China; tree-ring width; climate response; temperature

Date: 2016-12-15; Accepted Date: 2017-03-20

CAS “Light of West China” Program; National Natural Science Foundation of China (41671212, 41171170, 41371221); Project of State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences

CAI Qiufang, E-mail: caiqf@ieecas.cn

2016-12-15；录用日期：2017-03-20

中国科学院“西部之光”项目；国家自然科学基金项目（41671212，41171170，41371221）；中国科学院黄土与第四纪地质国家重点实验室开放基金项目

蔡秋芳，E-mail: caiqf@ieecas.cn

段丙闯, 蔡秋芳. 2017. 桂林三月岭马尾松树轮宽度对气候变化的响应[J]. 地球环境学报, 8(3): 243 – 252.

: Duan B C, Cai Q F. 2017. Responses of tree-ring width of Pinus massoniana to climate change in Guilin, Guangxi Province [J]. Journal of Earth Environment, 8(3): 243 – 252.