地理科学进展, 2021, 40(2): 343-356 doi: 10.18306/dlkxjz.2021.02.014

研究综述

树轮密度对气候的响应及重建研究进展

邓国富1,2, 李明启,1,2,*

1. 中国科学院地理科学与资源研究所,中国科学院陆地表层格局与模拟重点实验室,北京 100101

2. 中国科学院大学,北京 100049

Advances of study on the relationship between tree-ring density and climate and climate reconstruction

DENG Guofu1,2, LI Mingqi,1,2,*

1. Key Laboratory of Land Surface Pattern and Simulation, Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China

2. University of Chinese Academy of Sciences, Beijing 100049, China

通讯作者: *李明启(1979— ),男,山东东明人,博士,副研究员,主要研究方向为树轮气候学。E-mail: limq@igsnrr.ac.cn

收稿日期: 2020-04-20   修回日期: 2020-08-21   网络出版日期: 2021-02-28

基金资助: 国家自然科学基金项目.  41977391
国家自然科学基金项目.  41630529
国家自然科学基金项目.  41571194

Received: 2020-04-20   Revised: 2020-08-21   Online: 2021-02-28

Fund supported: National Natural Science Foundation of China, No.  41977391
National Natural Science Foundation of China, No.  41630529
National Natural Science Foundation of China, No.  41571194

作者简介 About authors

邓国富(1996— ),女,湖南株洲人,硕士生,主要研究方向为气候与环境变化。E-mail: dengguofu18@mails.ucas.edu.cn

摘要

树轮密度是重建过去气候变化的重要代用指标。过去几十年,国内外学者针对树轮密度开展了大量的研究工作,取得了丰硕的成果。论文收集了国内外有关树轮密度的研究文献,从树轮密度对气候变化的响应和气候重建等方面进行了概述。尽管树轮密度对气候变化的响应可能受到样点的地理位置、地形(如海拔、坡向等)、树种等因素的影响,但湿冷地区树轮晚材最大密度能较好地反映生长季或生长季末期温度变化;而干旱地区的早材平均密度能揭示生长季早期的降水变化。因此,目前大多数研究以重建生长季或生长季末期温度为主,也有重建降水量、海平面气压变化的案例。此外,重建方程的方差解释量存在空间和树种间的差异,但这些差异需要更多的研究工作加以验证。最长的重建序列来自于欧洲,长达2018 a;国内基于树轮密度重建的最长温度序列来自于西南地区,长达449 a。另外,树轮密度序列在对温度响应的分异性问题以及评估强火山喷发的降温效应方面还存在争议。其争议源于树轮密度与温度关系还受到其他环境变量与人为因素的影响,因此今后研究还需关注树轮密度与降水、光照等其他气候变量的关系,同时也要注重树种、海拔、实验方法等非气候因子的影响。

关键词: 树轮密度 ; 气候变化 ; 温度重建 ; 火山喷发

Abstract

Tree-ring density is one of the important proxies for climate reconstruction. In the past few decades, international and Chinese scholars have carried out a large number of studies on the relationship between tree-ring density and climate, and many climate series were reconstructed for different regions. In this study, we collected publications about the responses to climate change and climate reconstruction based on tree ring density, and reviewed some progresses in these areas. Although the relationship between tree-ring density and the climate may be affected by tree species, terrain (altitude, slope direction, and so on), and other factors, the maximum latewood density is a well-known proxy of summer or early fall temperature variation in cold and moist areas, while earlywood density can be used to reconstruct precipitation variation in arid areas. Therefore, most of the scholars reconstructed the temperature of the summer or the late fall, and only few precipitation and sea level pressure series were reconstructed based on tree-ring density. In addition, the variance that the reconstruction can explain varies from region to region, and from tree species to species. To date, the reconstruction series spanned 2018 years in Europe, which is the longest in the world, and 449 years in Southwest China, which is the longest series in China. There still exist controversies around the divergence in tree-ring density response to temperature and the cooling effect of volcanic eruptions, because of the unstable relationship between tree-ring density and temperature. Future research should pay more attention to the relationship between tree-ring density and other climate variables such as precipitation and light, and the influence of nonclimatic factors such as altitude, tree species, and experimental methods.

Keywords: tree-ring density ; climate change ; temperature reconstruction ; volcanic eruption

PDF (3272KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

邓国富, 李明启. 树轮密度对气候的响应及重建研究进展. 地理科学进展[J], 2021, 40(2): 343-356 doi:10.18306/dlkxjz.2021.02.014

DENG Guofu, LI Mingqi. Advances of study on the relationship between tree-ring density and climate and climate reconstruction. Progress in Geography[J], 2021, 40(2): 343-356 doi:10.18306/dlkxjz.2021.02.014

树木年轮具有分辨率高、定年准确、分布广泛等特点,是重建过去气候变化的重要代用资料[1,2,3],常用的代用指标有宽度、密度和同位素[4,5,6,7]。树轮密度指木材的组织密度,包括早材最小密度(minimum earlywood density, MID)、早材平均密度(mean earlywood density, EWD)、晚材平均密度(mean latewood density, LWD)和晚材最大密度(maximum latewood density, MXD),其大小与早晚材细胞壁的厚度及细胞的大小有关。相对于宽度,树轮密度受生长季前期气候条件影响较小,能够更好地揭示生长季气候变化[8,9,10]。20世纪60年代,国际上开启了基于树轮密度的气候重建工作[11]。之后,MXD被广泛应用于区域气候重建[12,13,14,15,16]。国内起步相对较晚,始于20世纪90年代[17,18]。目前,国内外学者一方面力求利用密度年表重建更长时间尺度的气候序列[19,20,21];另一方面则是建立高分辨率的树轮密度网络,以此提高重建结果可信度,更好地揭示区域气候变化过程[21,22,23,24]

本文将从研究区分布、树种、使用的密度指标、重建的长度等几个角度出发,综述国内外近几十年来的树轮密度对气候的响应及其重建工作,进一步探讨树轮密度研究的潜在问题,发现树轮密度的研究潜力,为今后树轮密度研究提供背景参考。

1 树轮密度对气候的响应研究

1.1 树轮密度与气候关系

揭示树轮密度指标与气候要素关系是重建过去气候变化的前提。近几十年来,国内外学者开展了大量树轮密度对气候变化的响应研究,取得了丰硕的研究成果。树轮密度的大小取决于树轮细胞的大小及其细胞壁的厚度。生长季前期的形成层活动较强,细胞分裂和扩大,树木以径向生长为主;生长季后期形成层活动减弱,细胞中合成细胞壁的物质沉积,细胞壁加厚,因此树轮早材密度较小而晚材密度较大[25,26,27]。生长季或生长季末期较高的温度能加强植物的光合作用,有利于细胞壁加厚,增加树轮晚材密度[25],尤其是在树木生长受温度限制的地区,树轮晚材密度对温度的响应更敏感。如在Pyrenees、Alps、Carpathians、Tatras、Altai、Norikura等高山树线地区[27,28,29,30,31,32,33,34](表1),以及在北美、西伯利亚北部等高纬树木生长北界地区,其MXD与生长季温度有较好的正相关关系[35,36,37,38]。中国天山、长白山、崆峒山等山区的MXD也是对夏季温度响应最强,但相关性偏低[39,40,41,42,43,44,45](表2)。此外,MXD在反映极端温度和相应的大气环流形势上也具有一定的潜力,如芬诺斯坎迪亚的34个MXD年表的第一主成分能反映1901—1978年夏季极端高温天数的变化,方差解释量为50%,并通过对高度场、风场的分析发现:2个序列的正极值年份,反气旋多分布于芬诺斯坎迪亚的中东部,常伴有高空急流的经向移动[46]

表1   国际上树轮密度对气候响应研究

Tab.1  International studies on responses of tree-ring density to climate change

研究区树种拉丁名指标海拔/m气候变量相关系数/时段文献来源
Pyrenees欧洲赤松
山赤松
欧洲冷杉
Pinus sylvestris
Pinus uncinata
Abies alba
MXD2010T(5—9月)0.80/1952—2002年[27]
Mount Norikura富士山冷杉Abies veitchiiMXD1600T(7、9月)0.30/1979—2005年[28]
Mount Norikura富士山冷杉Abies veitchiiMXD1900T(7—9月)0.35/1979—2005年[28]
Mount Norikura大白叶冷杉Abies mariesiiMXD1900T(7—9月)0.34/1979—2005年[28]
Mount Norikura大白叶冷杉Abies mariesiiMXD2400T(7—9月)0.44/1979—2005年[28]
Rhaetian Alps, Italy瑞士五针松Pinus cembraMXD2075~2350T(5—9月)0.60/1901—2015年[29]
Eastern Carpathians欧洲赤松Pinus sylvestrisMXD750~1200Tmax(7月30日—9月24日)0.63/1961—2013年[30]
Tatras欧洲云杉Picea abiesMXD1480T(4—9月)0.62/1901—2004年[31]
Dolina Mengusovska欧洲落叶松Larix deciduaMXD1450T(5—6月)0.68/1951—2012年[32]
Village Vernar欧洲落叶松Larix deciduaMXD850P(7月)0.43/1951—2012年[32]
Altai西伯利亚落叶松Larix sibiricaMXD2380T(6—7月)0.57/1963—2000年[33]
Yakutia卡氏落叶松Larix cajanderiMXD22T(6—7月)0.55/1950—2000年[33]
Alps欧洲落叶松Larix deciduaMXD>1500T(6—9月)0.73/1911—2003年[34]
Franklin Mountains白云杉Picea glaucaMXD653T(5—8月)0.45/1944—1977年[35]
Hornby Cabin白云杉Picea glaucaMXD143T(5—8月)0.58/1944—1977年[35]
Cri Lake白云杉Picea glaucaMXD108T(5—9月)0.50/1944—1977年[35]
Coppermine白云杉Picea glaucaMXD200T(5—8月)0.48/1944—1977年[35]
North-eastern Finland欧洲赤松Pinus sylvestrisMXD200T(4—9月)0.68/1876—2013年[36]
Laanila欧洲赤松Pinus sylvestrisMXD220T(6—8月)0.61/1958—2002年[37]
Suntar Khayata卡氏落叶松Larix cajanderiMXD900Tmin(6—8月)0.67/1929—2000年[38]
Iberian Peninsula西班牙刺柏Juniperus thuriferaMID530~1375P(4—6月)-0.67/1951—2000年[47]
Southern Urals欧洲赤松Pinus sylvestrisMID740P(5月)-0.58/1950—2002年[48]
Khangai西伯利亚落叶松Larix sibiricaMID1920P(6月)-0.57/1950—2002年[48]
Sierra de Gúdar欧洲黑松Pinus nigraMID1090P(5月)-0.65/1950—2002年[48]

注:T表示平均温度,Tmax表示平均最高温度,Tmin表示平均最低温度;P表示降水。下同。

新窗口打开| 下载CSV


表2   中国树轮密度对气候响应研究

Tab.2  Studies on responses of tree-ring density to climate change in China

样点树种海拔/m指标响应变量相关系数/时段文献来源
吉普克雪岭云杉2555MXDT(8月)0.35/1959—2004年[39]
阿乌里亚乔克山雪岭云杉2690MXDT(4—8月)0.42/1961—2005年[39]
小白代沟雪岭云杉2682MXDT(6—8月)0.33/1956—2006年[39]
长白山长白落叶松1848MXDTmax(6—8月)0.38/1958—2008年[40]
长白山长白落叶松1585MXDTmax(4—9月)0.44/1958—2008年[40]
长白山长白落叶松983MXDTmin(8月)-0.26/1958—2008年[40]
漠河兴安落叶松650MXDTmax(7月)0.45/1922—1998年[41]
漠河兴安落叶松650MIDTmin(3月)-0.323/1922—1998年[41]
漠河樟子松650MXDTmax(8月)0.49/1922—1998年[41]
漠河樟子松650MIDP(4月)0.343/1922—1998年[41]
崆峒山油松1950MIDP(6—8月)-0.4/1951—2006年[42]
崆峒山油松1950EWDP(6—8月)-0.36/1951—2006年[42]
崆峒山油松1950MXDT(6—8月)-0.26/1951—2006年[42]
崆峒山油松1950LWDT(6—9月)-0.28/1951—2006年[42]
艾肯达坂天山云杉2455MXDT(7—8月)0.30/1968—2005年[43]
巩乃斯林场北天山云杉1982EWDP(6—7月)-0.50/1968—2005年[43]
艾肯达坂雪岭云杉2450MXDTmax(5—8月)0.304/1958—2008年[44]
小五台山青杨(雌)1600MXDTmax(8月)0.348/1982—2011年[45]
小五台山青杨(雄)1600MXDTmax(4月)-0.429/1982—2011年[45]
祁连山中段青海云杉3140~3390EWDP(5月)-0.623/1956—2009年[49]

新窗口打开| 下载CSV


在树木生长受温度限制的地区,早材密度对温度的响应不同于晚材密度,温度的升高更有利于早材细胞的径向增大,而细胞壁加厚是次要的,因此高温导致早材密度更小[26]。但干旱地区的树轮早材密度对生长季早期的降水更为敏感,当春季干旱稍明显时,早材细胞径向扩张减弱,引起密度增大[47,48];当干旱较严重时,早材细胞脱水,原生质与细胞壁收缩,在细胞壁上形成很多锐利的折叠,撕破原生质结构,导致早材密度增大[42-43,49]。如Iberian Peninsula、Southern Urals Khangai、Sierra de Gúdar等气候干冷地区的MID与5—6月的降水呈显著的负相关关系[47,48](表1)。中国崆峒山的早材密度与6—8月的降水呈负相关关系[42],祁连山中段的早材密度与5月降水负相关[49],巩乃斯林场北的早材密度与6—7月的降水负相关[43](表2)。

此外,树轮密度年内变化受到气候季节变化的影响。地中海地区生长季早期的干旱使树轮密度突然增大,秋季的多雨使晚材密度突然变小,从而形成伪轮[50,51,52]。但由于早材树轮细胞的生长要消耗植物体内存储的养料,即早材密度还受到上年冬季气候的影响,所以早材中出现伪轮的情况较少,只有在极端干旱、寒冷的年份下才出现。而晚材密度主要受当年生长季气候的影响,所以出现变异的频次较多[50]

1.2 影响树轮密度与气候关系的因素

树木生长环境和自身遗传因素均能影响树轮密度与气候的关系,不同地形条件下决定树轮密度大小的主要环境因子可能不同,而遗传因素控制了树木应对环境变化的生理机制[3,5]。此外,树轮密度的大小还可能受人为干扰、虫害、获取树轮密度方法的影响,掩盖了原本气候信号。因此,影响树轮密度与气候关系的要素和过程相对复杂。

海拔对树轮密度与气候关系的影响主要体现在气候条件和土壤持水能力差异2个方面[32]。其一,高海拔地区气候寒冷,其MXD的增长主要受生长季温度的控制[32],如Tatras、Norikura及长白山地区高海拔的MXD与生长季温度的相关明显高于低海拔[28,32,40]。其二,不同海拔地区土壤持水能力的差异可能影响树轮密度对降水的响应,如Tatras低海拔山脉的土壤由石灰岩风化而成,储水能力差,树木易遭受干旱胁迫,MXD对7月干旱响应敏感;而高海拔地区土层的母质为花岗岩,土壤的蓄水能力较低海拔强,降水对MXD影响较弱[32]。此外,不同海拔MXD大小不同,并且在去除非气候因素影响时并不能完全消除此差异,如果不同海拔的样本量在时间上分布不均匀,会扭曲树轮密度年表的趋势,干扰气候信息的获取[53,54]

坡向影响树轮密度对季节气候变化的敏感程度,如在地中海地区,南坡光周期更长,蒸发量更大,水分条件相对较差,其树轮密度更易因水分条件的变化而产生变异,因此南坡树轮比北坡树轮更易形成伪轮[51]

由于遗传因素的影响,不同树种或同树种雌雄株应对环境变化的生理机制也有差异,从而影响树轮密度对气候变化的响应[41,45]。树种差异会影响树轮密度对气候响应的敏感程度和时间。如地中海地区栎树(Quercus ilex)扎根深,在干旱条件下仍可以发生蒸发;而松树扎根浅,易因干旱引起气孔关闭,减弱蒸发。因此,栎树的树轮细胞更易丢失水分而发生变形,易发生密度突变,对气候变化更敏感[52]。再如漠河兴安落叶松(Larix gmelinii)的晚材细胞加厚的活跃高峰可能比樟子松(Pinus sylvestris var. mongolica)早,其MXD对温度的最敏感时段比樟子松有所提前[41]。但也存在同一地区不同树种树轮密度与气候关系基本一致的情况,如日本Norikura山1900 m处富士山冷杉与大白叶冷杉,其MXD对气候的响应基本一致[28](表1)。因此,树种差异对树轮密度与气候关系的影响存在区域差异,仍需更多研究来验证。树木性别可能影响树轮密度对气候响应的时间,如小五台山青杨(Populus cathayana)雄株MXD对初春温度有负的响应,而雌株没有,其原因是合适的低温下,雌株萌发先于雄株,雄株在养分竞争中处于劣势,使得生长季节后期的形成层活动减弱,径向生长减缓,MXD增大[45]。但目前只见青杨雌雄株树轮密度对气候响应差异的报道,其他树种雌雄异株的影响尚未可知。

树轮密度的大小还可能受到人为干扰、虫害、树轮密度获取方法等的影响。Franceschini等[55]认为,间伐可使树木获取的光照和水分增加,抵抗环境胁迫的能力增强,其密度大小受气候变化的影响可能变小,但由于采样地生境较好,研究结果并未证实其假设。他指出需进一步在环境相对恶劣条件下采样来验证该假设。此外,树木因虫害而落叶,光合作用减弱,用于细胞径向生长和细胞壁增厚的光合产物减少,导致早材密度偏大而晚材密度偏小[56]。然而在出现气候极端的年份树轮密度也会出现类似的情况,并且极端气候也易引发虫害,因此树轮密度、虫害、极端气候之间的关系较为复杂,目前主要是结合同一地区非宿主植物的树轮年表以及虫害记录来剔除虫害噪声[6]。此外,近年相关研究表明,相比于X射线法,树轮径向解剖技术获得的参数能精确地反映密度变化,其获取的LWD与温度关系在高低频域上均较好[36,57]

2 气候变化重建研究进展

过去几十年,基于树轮密度与气候要素之间的关系,研究人员做了大量的气候重建工作,揭示了过去几百年,甚至千年的区域气候变化过程,为模式预测未来气候变化提供了基础数据。

2.1 国外基于树轮密度的气候重建研究

国际上的重建工作多分布在Scandinavia[57,58,59,60]、Pyrenees[61,62,63]、Pindus[64]、Scotland[65]、Alaska[66,67,68,69]等高山高纬地区,且基本上是以MXD来重建区域生长季(4—9月)或生长季晚期(7—9月)的温度历史,其在器测时期温度重建的方差解释量均超过了50%,表明这些地区的MXD能较好地反映研究区温度变化。而在美洲西部、澳大利亚东南部等地区生长季或夏季温度重建的方差解释量相对较低,如Washington地区1905—1983年生长季温度的重建方差解释量为39%[68],Tasmania中部地区1960—2009年1月最高温度重建的方差解释量为34%[70],这可能是由于树轮密度对温度变化敏感度降低所引起,也有可能是复本量较少、标准化方法不同等因素造成的。如Pyrenees中部地区复本量由216个增至414个,重建的方差解释量由28.09%(1944—2005年)增至51.84%(1950—2014年)[61,63]。虽然在Scandinavia地区研究中温度重建的方差解释量较高,但在AD 850—950年和AD 1100—1300年期间的复本量较低,导致重建不确定性较高。Linderholm等[57]在Zhang等[54]研究的基础上,增加了45个来自当地的干木材样本,删除了来自古建筑的样本,所重建的中世纪暖期的温度要比之前重建的低,但由于中世纪暖期的复本量还是偏少,目前该时期的温度变化还是难以确定。Dorado Liñán等[62]指出,当研究区地理范围较广时,直接使用所有测量序列建立区域标准化年表时,如没有满足区域标准化方法中复本量和年龄分类的要求,结果的可靠性可能会较差;而如果先建立不同样点的小区域年表、再将小区域年表合成大区域年表,则可能提高方差解释量。

除用树轮密度单一指标重建外,还有研究认为MXD与EWD均受春季气候影响,两者存在一定相关性[71],为剔除生长季前期气候影响,利用MXD与EWD差值重建夏季平均温度变化过程,方差解释量高达70%(1855—2006年)[59]。此外,也有研究尝试重建温度以外的其他气候变量。如Briffa等[72]用北半球341个MXD序列重建了海平面气压变化,但由于校准和验证时段过短、大气环流与树轮密度是间接关系等因素影响,重建结果并不理想。他指出,通过利用更长时间的器测温度序列,以及对降水敏感的代用资料,或许能改善重建结果的信度。

大多数针叶树种树轮密度均可用来进行气候重建(图1),但由于环境差异,导致重建所用树种也存在区域差异。欧洲Scandinavia山脉主要使用欧洲赤松[54,57,59-60,65],Pyrenees山脉采用山赤松(Pinus uncinata)[61,63],Pindus山脉采用波斯尼亚松(Pinus heldreichii)[64]。美洲Alaska地区采用白云杉[66,69,73],Maine地区采用红云杉(Picea rubens)[15],中亚地区主要为雪岭云杉(Picea schrenkiana)[74]和西伯利亚落叶松[75],日本中部为库页云杉(Picea glehnii)[76]。当进行大空间尺度重建时,研究者多利用不同树种的树轮密度资料,如Briffa等[16]采用云杉、落叶松、冷杉的MXD重建了北美北部地区夏季温度变化;Wiles等[68]用白云杉、北美云杉(Picea sitchensis)、花旗松(Pseudotsuga menziesii)和大果铁杉(Tsuga mertensiana)的MXD重建了阿拉斯加和太平洋西北部海岸4—9月的温度变化。

图1

图1   国际上树轮密度重建工作

Fig.1   International climate reconstruction based on tree ring density


就年表长度与重建序列长度来看,欧洲多千年长度的MXD年表,因Scandinavia地区的欧洲赤松活树的年龄可达500 a以上,死树年龄可达1000 a以上并且保存得较为完整[57],因此目前该区拥有世界上最长的MXD年表,长达2226 a(216 BC—AD 2010),由瑞典北部Torneträsk及Fennoscandia北部2条序列合成,并重建了过去2018 a(8 BC—AD 2010)温度变化[60]。北美洲地区在21世纪10年代也将MXD年表延长到了千年尺度,目前最长MXD年表及重建序列为930 a(AD 1073—2002),其研究区位于Alaska的Firth河流域[69]。中亚和日本地区目前的树轮密度年表与重建序列主要在三四百年间[74,75,76],相对偏短。

2.2 国内基于树轮密度的气候重建研究

中国重建工作多是基于几个样点树轮密度资料开展的区域气候重建,分布在西南、西北和东北地区(图2),其中西南主要分布在青藏高原东南缘[77,78]、横断山脉及周边地区[9,79-80];西北在天山[81]、祁连山[70]、贺兰山[82]等地;东北在长白山[40]及大兴安岭[83]。近几年来大空间尺度的重建也在开展,如Duan等[23]用了青藏高原地区17个采样点的树轮密度资料重建研究区自1572年来8—9月的温度变化,陈峰等[84]用新疆北部地区5个序列重建了研究区近353 a的5—8月温度变化。方差解释量同样存在空间差异(图2a),西南地区重建的结果相对较好,进一步验证了气候高寒区域的树轮晚材密度能够更好地揭示其生长季温度的历史变化。

图2

图2   中国树轮密度重建工作

注:本图基于自然资源部标准地图服务网站下载的审图号为GS(2019)1823号的标准地图制作,底图无修改。

Fig.2   Climate reconstruction based on tree-ring density in China


重建气候要素存在空间差异(图2b)。东北和西南地区夏季降水充沛,树木生长的主要限制因子为温度,多数工作采用MXD重建区域生长季温度[9,77,79,83,85]。而西北地区气候干旱,温度和降水均有可能限制树木生长[17],使用的密度指标与重建要素相对复杂,如石门山低海拔地区生长季早期气候暖干,温度升高加剧土壤水分的蒸发,树轮密度受温度影响,故其EWD被用于重建5—6月的平均最高温度变化[86];吉木萨尔、陕西黄陵生长季早期的树轮生长主要受降水影响,因而其MID被用于重建生长季降水变化[17,18];鄂尔多斯地处中国北方干旱半干旱地区,其树轮晚材形成受到夏末到冬初的降水的影响,降水增多使树木生长季延长,树轮晚材多形成内径小、壁质厚的分生细胞,使得树轮晚材密度偏大,因此其MXD被用于重建6—10月的降水变化[87]

国内用于重建的树种也较多。中国西南地区针叶树种资源丰富,先后有研究使用川西云杉(Picea likiangensis)[19,78,85,88-89]、紫果云杉(Picea purpurea)[9]、高山松(Pinus densata)[90]及峨眉冷杉(Abies fabri)[91]、急尖长苞冷杉(Abies georgei)[77]、大果云杉(Picea asperata)[79]、麦吊云杉(Picea brachytyla)[92]等树种MXD进行气候重建。但目前该研究区各树种重建的方差解释量存在差异,如峨眉冷杉和高山松的重建方差解释量在50%以上,而急尖长苞冷杉和怒江红杉的方差解释量只有39%。该差异可能来源于遗传因素、环境或实验误差的影响,也可能是由于树轮密度揭示气候变化的能力存在种间差异,今后还需采取同一地点不同树种的树轮密度资料进行验证。西北和东北地区用于重建的树种相对偏少,但均为各研究区内树轮生长对气候响应较为敏感的树种[6-7,84]:天山地区为雪岭云杉[18,81,93-94]、阿尔泰山为西伯利亚落叶松[84,95]、祁连山地区青海云杉(Picea crassifolia)[70],石门山、贺兰山、陕西黄陵等地为油松(Pinus tabuliformis)[17,82,86-87],东北地区主要是兴安落叶松[83]、长白落叶松(Larix olgensis)[96](图2c)。

相对于欧美地区,国内密度年表与重建的温度序列相对较短。在西南地区,青藏高原东部马尔康的树轮密度年表最长,为640 a (AD 1375—2014)[19];基于MXD重建的最长温度序列为449 a (AD 1563—2011)[89](图2d),来源于青藏高原东南部。西北地区阿尔泰山的树轮密度最长年表为521 a (AD 1488—2008)[84],但重建的最长温度序列来源于乌鲁木齐河源地区,为384 a (AD 1623—2006)[94]。东北长白山密度年表最长为393 a (AD 1615—2008),其重建序列为238 a (AD 1762—1999)[96]。由此可见,树轮密度年表最长的地方,其重建序列不一定最长,这表明需要更多年龄更大的样芯,通过增加更早历史时期的复本量来延长重建长度。

3 其他方面树轮密度研究工作

3.1 辅助火山喷发事件的定年

历史或地质时期火山喷发时间主要基于文献和冰芯记录数据来确定的,但有时候其定年存在误差[97]。由于树轮MXD与温度相关性高,其能为火山喷发事件提供准确定年。如Schneider等[98]根据树轮MXD重建的北半球温度序列,发现过去千年有14个显著降温事件,并与冰芯中的火山喷发记录相吻合,为火山喷发提供了确切年代;Esper等[22]用20个MXD年表重建了北半球过去600 a温度变化序列,结果显示1453年降温程度最大,表明15世纪中期的火山喷发事件发生在1452年,而不是1458年;Briffa等[99]分析了北半球区域树轮密度年表,推测Billy Mitchell火山爆发可能在1495年或1587年。

3.2 评估区域气候对火山爆发的响应

研究表明:赤道区域和北半球中低纬度的强火山喷发后1~2 a内,全球各地出现不同程度的降温,导致树轮密度偏低[16,97,100],因此常用树轮密度及其重建序列来研究火山喷发的降温效应。不同地区气候对火山喷发的响应存在差异,如高纬度地区对强火山喷发的响应要比为中纬度地区明显,延迟时间在1~2 a左右[101]。北美西北部地区对火山喷发的响应较强,而北美的中部地区则响应偏弱[102]。此外,在青藏高原地区树轮MXD中也发现了强的火山信号,青藏高原东南部在Tambora火山喷发之后有显著降温[79]。青藏高原对中高纬度的火山爆发响应很弱,对热带火山爆发有显著的降温响应,但该降温响应有可能被厄尔尼诺的影响所抵消[103]。总之,火山喷发对区域气候的影响与研究区地理位置、火山喷发的位置和强度有关[33]

关于用树轮密度资料估计火山喷发带来的降温程度存在争议。Mann等[104]认为,火山喷发后的降温年份可能发生丢轮,从而低估了火山喷发带来的降温效应。但Esper等[105]假设欧洲北部和中部的MXD区域年表在1816年丢年,发现1722—1815年的MXD与夏季温度的相关系数分别从0.53、0.41降至0.02、0.03,从而否定了Mann等[104]的理论。也有研究认为光强减弱也能影响光合产物的合成,从而使密度偏低,故用密度重建的温度高估了火山喷发的降温效应[106]。但Stoffel等[107]用树轮密度重建的温度变化与气候模式模拟的结果基本一致。Lea Schneider等[108]也指出树轮MXD受上个生长季树木生理条件的影响小,主要受当年生长季气候的影响,因此MXD能够较好地反映其区域在火山喷发后的降温程度。总之,丢轮、光照对树轮密度的影响,增加了重建温度的不确定性,降低了树轮密度对温度的代表性。

3.3 树轮密度的分异性问题

“分异性”问题是指基于树轮代用资料重建的温度序列在近几十年没有显著上升的趋势,与器测记录产生了显著差异[69,109]。最早在阿拉斯加地区白云杉树轮宽度与密度对气候响应的研究中发现“分异性”问题[67]。随后,在白令海峡与楚科奇地区[66]、芬诺斯坎迪亚北部[58]、比利牛斯山脉[63]等高纬度地区的树轮密度研究中也出现了“分异性”问题。中低纬度地区出现“分异性”问题的研究相对较少。在中国青藏高原东南部也存在该问题,Li等[79]发现在2005年后MXD与温度相关性降低。

关于造成“分异性”问题的原因,目前众说不一,主要的解释有3种。第一种观点认为“分异性”问题是由于低估降水影响造成的[27]。近几十年来,高纬度地区的增暖大于低纬度地区,可能使高纬地区树轮生长的主要限制因子发生了改变,从而使温度与树木生长之间相关性降低[109]。研究表明,近几十年来高纬地区树轮与降水的关系加强,干旱胁迫对树轮宽度与密度产生了很大的影响[63,67]。第二种观点认为由于近些年来北半球高纬度地区植物可利用的光合有效辐射减少,树木径向生长缓慢,影响了树轮对升温的响应[106,110]。第三种观点认为,标准化过程、空间分辨率、树种、树龄等因素影响了重建结果[62]。不同树种、不同年龄阶段的树木生长趋势存在差异,而在建年表的过程中采用同一生长曲线拟合,再进行去趋势处理,可能导致重建温度与器测温度之间的偏差[73],如Grudd等[58]指出,当样本中加入了很多比较小树龄的样本时会影响生长趋势的拟合,最终影响重建结果;Schneider等[108]重建的北半球温度序列显示在近几十年没有显著的上升趋势,这可能是由于树轮密度空间分辨率过小引起的。

4 总结与展望

本文通过对收集的国内外树轮密度研究成果进行归纳总结,试图得出一定的规律,挖掘树轮密度研究的潜力,以期为未来树轮密度研究提供一定借鉴与帮助。

综合国内外树轮密度响应研究工作,发现采样地多分布在高海拔或高纬地区,少数研究分布在干旱的山区。由于采样地地理位置、生长环境的差异,树轮密度对气候响应的关系错综复杂。在高大山体树木的生长上限及高纬地区树线北界,MXD能更好地揭示温度变化;而在干旱的地区,树轮密度序列中含有降水信号。然而目前鲜少有研究分析干热河谷地区的树轮密度对气候的响应,今后树轮密度研究可选取该区进行采样,由于该区气候高温低湿,有可能从其树轮密度年表中提取出降水信号[32]。此外,树轮密度与气候的关系还受到树种、海拔、坡向、人类活动、虫害等因子的影响,但各方面研究的数量偏少,且目前主要是根据相关系数或方差解释量大小来判断不同树种或不同海拔树轮密度与气候变量关系的差异,而该差异有可能是由样本量、实验误差等造成的,因此今后还需更多研究工作进一步验证。另外,间伐对树轮密度的影响本质是树木对光照、养分竞争带来的影响,但要寻找合适的间伐采样地还相对困难,今后研究可从种群密度的角度来分析竞争对树轮密度与气候关系的影响[60]

国内外树轮密度重建工作主要分布于高纬度或高海拔地区,以重建生长季或生长季末期的温度为主,但重建的方差解释量存在空间和树种差异,而该差异可能源于环境、实验误差、复本量、标准化方法的影响。亚洲和北美洲的树轮密度年表多在1000 a以下,比欧洲短。重建工作涉及的树种包括云杉、落叶松、冷杉、松树、黄杉、铁杉等多个属的树种,但不同树种树轮密度揭示气候变化的能力可能存在差异,因此找出适合于重建区域气候变化的最优树种是今后树轮密度研究关注的重点之一,而对于利用多个树种密度重建空间大尺度的气候变化时,则需要考虑树种差异的影响。此外,国际上各研究区基本上都是以晚材密度重建生长季温度历史为主,但其实早材密度指标也包含有丰富的气候信息[26,32],今后重建工作可基于其他密度指标,或利用不同的数学方法将不同密度指标合成新的指标来开展,从而更好地揭示过去气候变化。在重建的气候变量方面,可以考虑其他要素,如海平面气压变化等。

综合树轮密度在其他方面的应用来看,树轮密度用于评估火山爆发的降温程度还存在争议,原因是树轮密度与温度关系受到了光照、伪轮等其他因子的影响,因此,在严格控制数据质量的同时,今后需关注光照对树轮密度与温度关系的影响。对于树轮密度研究中存在的分异性问题,多数解释是由于温度与树轮密度之间的相关性降低,与降水的相关性升高,因此树轮密度与气候变量的关系不一定是线性关系,今后可考虑用非线性模型解决该问题[24]。此外,树轮密度与气候变量的关系并非一成不变,重建方程未必能代表所有时间段的树轮密度与气候因子的关系,因此今后的研究可比较不同校准方法重建同一个气候变量所产生的偏差[60]

总之,为提高树轮密度对温度等气候变量的代表性,确定树轮密度与各个环境变量和人类活动的关系、建立一个能同时表征树轮密度与各气候变量关系的数理统计模型是今后树轮密度研究要发展的方向之一。

参考文献

Cook E R, Seager R, Kushnir Y , et al.

Old World megadroughts and pluvials during the Common Era

[J]. Science Advances, 2015,1(10):e1500561. doi: 10.1126/sciadv.1500561.

URL     PMID:26601136      [本文引用: 1]

Climate model projections suggest widespread drying in the Mediterranean Basin and wetting in Fennoscandia in the coming decades largely as a consequence of greenhouse gas forcing of climate. To place these and other

Büntgen U, Myglan V S, Ljungqvist F C , et al.

Cooling and societal change during the Late Antique Little Ice Age from 536 to around 660 AD

[J]. Nature Geoscience, 2016,9(3):231-236.

DOI:10.1038/ngeo2652      URL     [本文引用: 1]

Fritts H C. Tree rings and climate [M]. London, UK: Academic Press, 1976: 434-505.

[本文引用: 2]

张芳芳, 郑永宏, 潘国艳 , .

神农架地区树轮δ 18O序列的气候指示意义

[J]. 地理科学进展, 2018,37(7):946-953.

DOI:10.18306/dlkxjz.2018.07.008      URL     [本文引用: 1]

Abies forrestii)树轮纤维素氧同位素的分析,建立了1790-2011时段的年分辨率的δ18O年表。将树轮δ18O序列与邻近的巴东气象站记录的气候数据、大气环流数据以及南海海温进行对比分析,研究了树轮δ18O序列对气候要素的响应。结果表明,神农架地区树轮δ18O序列与该地区气候要素密切相关:与当年6-7月降水量呈负相关、与当年3-4月温度呈正相关,与当年6-7月相对湿度呈负相关。同时,树轮δ18O序列与当年2-3月北半球副高面积指数(5E-360)、9月南海副高面积指数(100E-120E)、1月太平洋极涡面积指数(Zone 2,150E-120W)及4-5月南海海温(SST)呈显著正相关;与当年1月南方涛动指数(SOI)呈负相关。表明树轮δ18O序列在一定程度上记录了大尺度的大气环流信息。该研究有助于进一步认知神农架地区树轮稳定氧同位素在过去气候变化研究中的潜在价值,为进一步气候重建研究提供科学依据。]]>

[ Zhang Fangfang, Zheng Yonghong, Pan Guoyan , et al.

Climatic significance of tree-ring δ 18O in Shennongjia Mountain

Progress in Geography, 2018,37(7):946-953. ]

[本文引用: 1]

郑永宏, 张永, 邵雪梅 , .

大别山地区黄山松和油松树轮宽度的气候意义

[J]. 地理科学进展, 2012,31(1):72-77.

[本文引用: 2]

[ Zheng Yonghong, Zhang Yong, Shao Xuemei , et al.

Climate significance of tree ring width of Huangshan Pine and Chinese Pine in the Dabie Mountains

Progress in Geography, 2012,31(1):72-77. ]

[本文引用: 2]

孙宇, 王丽丽 .

全球落叶松属树轮气候学研究进展

[J]. 地理科学进展, 2013,32(12):1760-1770.

DOI:10.11820/dlkxjz.2013.12.005      URL     [本文引用: 3]

落叶松属常分布于高纬高海拔地区,年际生长节律明显,在全球气候变化重建研究中具有重要的地位。本文对欧洲、北美、亚洲地区主要的落叶松属树种近年来取得的成果进行了综合分析,以地区为单位,梳理研究思路,概述研究重点及特色。全球落叶松属研究以对气温敏感性较高的轮宽为主,约占72%;其次是最大晚材密度,约占16%。其中,尤以欧洲落叶松和西伯利亚落叶松研究成果最多。在各地区中,欧洲的研究成果约占61%,倾向于大空间尺度多树种综合研究。北美研究常结合特殊的地形条件及生态事件进行,落叶松叶蜂虫害信息的提取是研究的重点。亚洲地区以西伯利亚落叶松为主,在20 世纪90 年代就已出现大量的重建成果。中国落叶松属有10 种1 变种,但研究相对落后于其他地区,近年来才在一些树种上取得突破,并出现了树轮密度的研究成果,未来应重点转向空间“场”的研究和多树种综合研究。

[ Sun Yu, Wang Lili .

Global research progresses in dendroclimatology of Larix miller

Progress in Geography, 2013,32(12):1760-1770. ]

[本文引用: 3]

王亚军, 李明启 .

中国利用树轮资料重建干湿变化研究进展

[J]. 地理科学进展, 2016,35(11):1397-1410.

DOI:10.18306/dlkxjz.2016.11.010      URL     [本文引用: 2]

树木年轮资料定年准确、分辨率高、连续性强,已成为研究过去全球变化的重要信息来源之一。中国树轮与干湿变化研究,样点分布广泛,研究树种较多,重建尺度较长。基于树轮宽度、密度和同位素等各种代用指标,在分析树木径向生长对干湿变化响应、重建历史时期干湿变化等方面取得很大进展。一般来说,在干旱半干旱区的低海拔区域,水分条件,尤其是春季、秋季和年湿润状况,是树木径向生长的重要影响因素。森林上限、相对冷湿区域以及亚热带气候区,树木生长与降水的相关性偏弱。在时间尺度上,基于树轮重建的长于1000年甚至超过2000年的干湿变化序列已有多条,主要分布在青藏高原地区。目前最长年表已达4500年,最长降水序列为3500年。本文通过对部分树轮成果的简要概括,以期为树轮采样和干湿变化重建提供借鉴与帮助。

[ Wang Yajun, Li Mingqi .

Research progress of dry-wet climate reconstruction by tree-ring in China

Progress in Geography, 2016,35(11):1397-1410. ]

[本文引用: 2]

Liang E Y, Lu X M, Ren P , et al.

Annual increments of juniper dwarf shrubs above the tree line on the central Tibetan Plateau: A useful climatic proxy

[J]. Annals of Botany, 2012,109(4):721-728.

DOI:10.1093/aob/mcr315      URL     [本文引用: 1]

Forty-three discs from 33 stems of Wilson juniper were sampled near the north-eastern shore of the Nam Co (Heavenly Lake). Cross-dating was performed along two directions of each stem, avoiding the compression-wood side as far as possible. A ring-width chronology was developed after a negative exponential function or a straight line of any slope had been fit to the raw measurements. Then, correlations were calculated between the standard ring-width chronology and monthly climate data recorded by a weather station around 100 km away.Our study has shown high dendrochronological potential of Wilson juniper, based on its longevity (one individual was 324 years old), well-defined growth rings, reliable cross-dating between individuals and distinct climatic signals reflected by the ring-width variability. Unlike dwarf shrubs in the circum-arctic tundra ecosystem which positively responded to above-average temperature in the growing season, moisture turned out to be growth limiting for Wilson juniper, particularly the loss of moisture caused by high maximum temperatures in MayJune.Because of the wide distribution of shrub and dwarf shrub species on the central Tibetan Plateau, an exciting prospect was opened up to extend the presently existing tree-ring networks far up into one of the largest tundra regions of the world.]]>

Yin H, Liu H B, Linderholm H W , et al.

Tree ring density-based warm-season temperature reconstruction since A.D. 1610 in the eastern Tibetan Plateau

[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2015,426:112-120.

DOI:10.1016/j.palaeo.2015.03.003      URL     [本文引用: 4]

Jones P D, Melvin T M, Harpham C , et al.

Cool North European summers and possible links to explosive volcanic eruptions

[J]. Journal of Geophysical Research: Atmospheres, 2013,118(12):6259-6265.

DOI:10.1002/jgrd.50513      URL     [本文引用: 1]

Fritts H C, Blasing T J, Hayden B P , et al.

Multivariate techniques for specifying tree-growth and climate relationships and for reconstructing anomalies in paleoclimate

[J]. Journal of Applied Meteorology, 1971,10(5):845-864.

DOI:10.1175/1520-0450(1971)010<0845:MTFSTG>2.0.CO;2      URL     [本文引用: 1]

Jacoby G C, Ivanciu I S, Ulan L D .

A 263-year record of summer temperature for northern Quebec reconstructed from tree-ring data and evidence of a major climatic shift in the early 1800's

[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 1988,64:69-78.

DOI:10.1016/0031-0182(88)90143-5      URL     [本文引用: 1]

Schweingruber F H, Bartholin T, Schaur E , et al.

Radiodensitometric-dendroclimatological conifer chronologies from Lapland (Scandinavia) and the Alps (Switzerland)

[J]. Boreas, 1988,17(4):559-566.

DOI:10.1111/bor.1988.17.issue-4      URL     [本文引用: 1]

Hughes M K, Schweingruber F H, Cartwright D , et al.

July-August temperature at Edinburgh between 1721 and 1975 from tree-ring density and width data

[J]. Nature, 1984,308:341-344.

DOI:10.1038/308341a0      URL     [本文引用: 1]

Conkey L E .

Red spruce tree-ring widths and densities in Eastern North America as indicators of past climate

[J]. Quaternary Research, 1986,26(2):232-243.

[本文引用: 2]

Briffa K R, Jones P D, Schweingruber F H.

Summer temperatures across northern North America: Regional reconstructions from 1760 using tree-ring densities

[J]. Journal of Geophysical Research: Atmospheres, 1994,99(D12):25835-25844.

[本文引用: 3]

刘禹, 吴祥定, 安芷生 , .

树轮密度、稳定C同位素对过去近100 a陕西黄陵季节气温与降水的恢复

[J]. 中国科学(地球科学), 1997,27(3):271-276.

[本文引用: 4]

[ Liu Yu, Wu Xiangding, An Zhisheng, et a1.

Seasonal precipitation and temperature reconstruction based on tree-ring density and stable carbon isotope

Scientia Sinica Terrae, 1997,27(3):271-276. ]

[本文引用: 4]

张志华, 李骥 .

用树轮密度及宽度资料重建新疆吉木萨尔县的季节降水和最高温度

[J]. 气象学报, 1998,56(1):77-86.

[本文引用: 3]

[ Zhang Zhihua, Li Ji .

Precipitation and average monthly high temperature in the Jimusare, Xinjiang as reconstructed from tree density and tree widths

Acta Meteorologica Sinica, 1998,56(1):77-86. ]

[本文引用: 3]

Duan J P, Li L, Ma Z G , et al.

Post-industrial late summer warming recorded in tree-ring density in the eastern Tibetan Plateau

[J]. International Journal of Climatology, 2020,40(2):795-804.

[本文引用: 3]

Esper J, Klippel L, Krusic P J , et al.

Eastern Mediterranean summer temperatures since 730 CE from Mt. Smolikas tree-ring densities

[J]. Climate Dynamics, 2020,54:1367-1382.

[本文引用: 1]

St. George S, Esper J.

Concord and discord among Northern Hemisphere paleotemperature reconstructions from tree rings

[J]. Quaternary Science Reviews, 2019,203:278-281.

[本文引用: 2]

Esper J, Büntgen U, Hartl-Meier C , et al.

Northern Hemisphere temperature anomalies during the 1450 s period of ambiguous volcanic forcing

[J]. Bulletin of Volcanology, 2017,79(6):41. doi: 10.1007/s00445-017-1125-9.

[本文引用: 2]

Duan J P, Ma Z G, Li L , et al.

August-September temperature variability on the Tibetan Plateau: Past, present and future

[J]. Journal of Geophysical Research: Atmospheres, 2019,124(12):6057-6068.

[本文引用: 2]

Ljungqvist F C, Thejll P, Björklund J , et al.

Assessing non-linearity in European temperature-sensitive tree-ring data

[J]. Dendrochronologia, 2020,59:125652. doi: 10.1016/j.dendro.2019.125652.

[本文引用: 2]

Björklund J, von Arx G, Nievergelt D, et al.

Scientific merits and analytical challenges of tree‐ring densitometry

[J]. Reviews of Geophysics, 2019,57(4):1224-1264.

[本文引用: 2]

Björklund J, Seftigen K, Schweingruber F , et al.

Cell size and wall dimensions drive distinct variability of earlywood and latewood density in Northern Hemisphere conifers

[J]. New Phytologist, 2017,216(3):728-740.

[本文引用: 3]

Büntgen U, Frank D, Trouet V , et al.

Diverse climate sensitivity of Mediterranean tree-ring width and density

[J]. Trees, 2010,24(2):261-273.

[本文引用: 4]

Takahashi K, Okuhara I, Tokumitsu Y , et al.

Responses to climate by tree-ring widths and maximum latewood densities of two Abies species at upper and lower altitudinal distribution limits in central Japan

[J]. Trees, 2011,25(4):745-753.

[本文引用: 8]

Cerrato R, Salvatore M C, Gunnarson B E , et al.

A Pinus cembra L. tree-ring record for late spring to late summer temperature in the Rhaetian Alps, Italy

[J]. Dendrochronologia, 2019,53:22-31.

[本文引用: 1]

Nagavciuc V, Roibu C-C, Ionita M , et al.

Different climate response of three tree ring proxies of Pinus sylvestris from the Eastern Carpathians, Romania

[J]. Dendrochronologia, 2019,54:56-63.

[本文引用: 2]

Kaczka R J, Spyt B, Janecka K , et al.

Different maximum latewood density and blue intensity measurements techniques reveal similar results

[J]. Dendrochronologia, 2018,49:94-101.

DOI:10.1016/j.dendro.2018.03.005      URL     [本文引用: 2]

Klippel L, Büntgen U, Konter O , et al.

Climate sensitivity of high- and low-elevation Larix decidua MXD chronologies from the Tatra Mountains

[J]. Dendrochronologia, 2020,60:125674. doi: 10.1016/j.dendro.2020.125674.

[本文引用: 9]

Churakova O V, Fonti M V, Saurer M , et al.

Siberian tree-ring and stable isotope proxies as indicators of temperature and moisture changes after major stratospheric volcanic eruptions

[J]. Climate of the Past, 2019,15(2):685-700.

[本文引用: 4]

Büntgen U, Frank D C, Nievergelt D , et al.

Summer temperature variations in the European Alps, A.D. 755-2004

[J]. Journal of Climate, 2006,19(21):5606-5623.

[本文引用: 2]

D'Arrigo R D, Jacoby G C, Free R M.

Tree-ring width and maximum latewood density at the North American tree line: Parameters of climatic-change

[J]. Canadian Journal of Forest Research, 1992,22(9):1290-1296.

[本文引用: 5]

Björklund J, Seftigen K, Fonti P , et al.

Dendroclimatic potential of dendroanatomy in temperature-sensitive Pinus sylvestris

[J]. Dendrochronologia, 2020,60:125673. doi: 10.1016/j.dendro.2020.125673.

[本文引用: 3]

Campbell R, McCarroll D, Loader N J, et al.

Blue intensity in Pinus sylvestris tree-rings: Developing a new palaeoclimate proxy

[J]. The Holocene, 2007,17(6):821-828.

[本文引用: 2]

Kirdyanov A V, Treydte K S, Nikolaev A , et al.

Climate signals in tree-ring width, density and δ 13C from larches in Eastern Siberia (Russia)

[J]. Chemical Geology, 2008,252(1/2):31-41.

[本文引用: 2]

Yuan Y J, Zhang T W, Wei W S , et al.

Development of tree-ring maximum latewood density chronologies for the western Tien Shan Mountains, China: Influence of detrending: Method and climate response

[J]. Dendrochronologia, 2013,31(3):192-197.

[本文引用: 4]

Shen C C, Wang L L, Li M Y .

The altitudinal variability and temporal instability of the climate-tree-ring growth relationships for Changbai larch (Larix olgensis Henry) in the Changbai mountains area, Jilin, Northeastern China

[J]. Trees, 2016,30(3):901-912.

[本文引用: 6]

王丽丽, 邵雪梅, 黄磊 , .

黑龙江漠河兴安落叶松与樟子松树轮生长特性及其对气候的响应

[J]. 植物生态学报, 2005,29(3):380-385.

[本文引用: 7]

[ Wang Lili, Shao Xuemei, Huang Lei , et al.

Tree-ring characteristics of Larix gmelinii and Pinus sylvestris var. mongolica and their response to climate in Mohe, China

Acta Phytoecologica Sinica, 2005,29(3):380-385. ]

[本文引用: 7]

侯迎, 王乃昂, 张学敏 , .

高频光密度测量法在崆峒山树轮年表建立中的应用

[J]. 干旱区地理, 2010,33(2):236-242.

[本文引用: 7]

[ Hou Ying, Wang Nai'ang, Zhang Xuemin, et a1.

Application of high- frequency densitometry in establishment of chronology from Mt. Kongtong, Gansu Province

Arid Land Geography, 2010,33(2):236-242. ]

[本文引用: 7]

喻树龙, 袁玉江, 陈峰 , .

巩乃斯河源树木年轮密度年表特征分析

[J]. 沙漠与绿洲气象, 2010,4(4):6-11.

[本文引用: 5]

[ Yu Shulong, Yuan Yujiang, Chen Feng , et al.

Tree-ring density chronology features of Gongnaisi riverhead area in western Tianshan Mountain

Desert and Oasis Meteorology, 2010,4(4):6-11. ]

[本文引用: 5]

张同文, 袁玉江, 喻树龙 , .

树轮灰度与树轮密度的对比分析及其对气候要素的响应

[J]. 生态学报, 2011,31(22):6743-6752.

[本文引用: 2]

[ Zhang Tongwen, Yuan Yujiang, Yu Shulong , et al.

Contrastive analysis and climatic response of tree-ring gray values and tree-ring densities

Acta Ecologica Sinica, 201l, 31(22):6743-6752. ]

[本文引用: 2]

黄科朝, 胥晓, 李霄峰 , .

小五台山青杨雌雄植株树轮生长特性及其对气候变化的响应差异

[J]. 植物生态学报, 2014,38(3):270-280.

[本文引用: 5]

[ Huang Kechao , Xu xiao, Li Xiaofeng, et al.

Gender-specific characteristics of tree-ring growth and differential responses to climate change in the dioecious tree Populus cathayana in Xiaowutai Mountains, China

Chinese Journal of Plant Ecology, 2014,38(3):270-280. ]

[本文引用: 5]

Zhang P, Ionita M, Lohmann G , et al.

Can tree-ring density data reflect summer temperature extremes and associated circulation patterns over Fennoscandia?

[J]. Climate Dynamics, 2017,49(7/8):2721-2736.

[本文引用: 1]

Camarero J J, Rozas V, Olano J M .

Minimum wood density of Juniperus thurifera is a robust proxy of spring water availability in a continental Mediterranean climate

[J]. Journal of Biogeography, 2014,41(6):1105-1114.

[本文引用: 3]

Camarero J J, Fernández-Pérez L, Kirdyanov A V , et al.

Minimum wood density of conifers portrays changes in early season precipitation at dry and cold Eurasian regions

[J]. Trees, 2017,31(5):1423-1437.

[本文引用: 5]

杨银科, 黄强, 刘禹 , .

云杉树轮生长密度对气候要素的响应分析

[J]. 西安理工大学学报, 2012,28(4):432-438.

[本文引用: 3]

[ Yang Yinke, Huang Qiang, Liu Yu , et al.

Response analysis between climate factors and the density of wood growth of Picea crassifolia

Journal of Xi'an University of Technology, 2012,28(4):432-438. ]

[本文引用: 3]

Kurz-Besson C B, Lousada J L, Gaspar M J, et al.

Effects of recent minimum temperature and water deficit increases on Pinus pinaster radial growth and wood density in Southern Portugal

[J]. Frontiers in Plant Science, 2016,7:1170. doi: 10.3389/fpls.2016.01170.

URL     PMID:27570527      [本文引用: 2]

Zalloni E, Battipaglia G, Cherubini P , et al.

Site conditions influence the climate signal of intra-annual density fluctuations in tree rings of Q-ilex L

.[J]. Annals of Forest Science, 2018,75(3):68. doi: 10.1007/s13595-018-0748-0.

[本文引用: 2]

Zalloni E, Battipaglia G, Cherubini P , et al.

Contrasting physiological responses to Mediterranean climate variability are revealed by intra-annual density fluctuations in tree rings of Quercus ilex L. and Pinus pinea L

[J]. Tree Physiology, 2018,38(8):1213-1224.

DOI:10.1093/treephys/tpy061      URL     PMID:29920596      [本文引用: 2]

Mixed forests of Quercus ilex L. and Pinus pinea L. are widely found throughout the Mediterranean Basin, being representative of two co-existing functional types: evergreen-sclerophyllous drought-resistant species and Mediterranean-adapted drought-avoidant conifers. Their contrasting physiological strategies to cope with water deficit influence all the processes regulating their growth such as wood formation, leading to peculiar tree-ring anatomical features such as intra-annual density fluctuations (IADFs). Intra-annual density fluctuations are abrupt changes in wood anatomical traits within a tree ring, appearing as latewood-like cells within earlywood or earlywood-like cells within latewood, and are frequently found in Mediterranean species as a response to seasonal climate changes. In this study, we characterized the anatomical traits and composition of carbon and oxygen stable isotopes in IADFs occurring in tree rings of Q. ilex and P. pinea trees co-existing at a same site in Southern Italy, in order to link their xylem hydraulic properties with the related physiological mechanisms. The relationships between IADF occurrence and seasonal mean temperature and total precipitation were investigated, with the aim of assessing whether they can be used as indicators of species-specific responses to intra-annual climate fluctuations. Results show that IADF period of formation is during autumn months for both species. The influence of climate on IADF occurrence was found to be an indicator of species-specific response to climate: an increased stomatal conductance associated to the formation of a wood safer against embolism was found in Q. ilex, while a tighter stomatal control associated to a more efficient wood with regard to hydraulic conductivity occurred in P. pinea. Moreover, the assessment of the influence of climate on IADF occurrence indicates that, with rising temperatures, Q. ilex would form fewer IADFs compared with P. pinea. Other study cases are desirable to assess the suggested forecasts and to link the plasticity of the species to form IADFs with their effective adaptive capability to compete for resources, and to explain how it may influence future population development.

Zhang P, Björklund J, Linderholm H W .

The influence of elevational differences in absolute maximum density values on regional climate reconstructions

[J]. Trees, 2015,29(4):1259-1271.

[本文引用: 1]

Zhang P, Linderholm H W, Gunnarson B E , et al. 1200

years of warm-season temperature variability in central Scandinavia inferred from tree-ring density

[J]. Climate of the Past, 2016,12(6):1297-1312.

[本文引用: 3]

Franceschini T, Gauthray-Guyénet V, Schneider R , et al.

Effect of thinning on the relationship between mean ring density and climate in black spruce (Picea mariana (Mill.) B.S.P.)

[J]. Forestry: An International Journal of Forest Research, 2018,91(3):366-381.

[本文引用: 1]

Koprowski M, Duncker P .

Tree ring width and wood density as the indicators of climatic factors and insect outbreaks affecting spruce growth

[J]. Ecological Indicators, 2012,23:332-337.

[本文引用: 1]

Linderholm H W, Gunnarson B E .

Were medieval warm-season temperatures in Jämtland, central Scandinavian Mountains, lower than previously estimated?

[J]. Dendrochronologia, 2019,57:125607. doi: 10.1016/j.dendro.2019. 125607.

[本文引用: 5]

Grudd H.

Torneträsk tree-ring width and density AD 500-2004: A test of climatic sensitivity and a new 1500-year reconstruction of north Fennoscandian summers

[J]. Climate Dynamics, 2008,31(7/8):843-857.

[本文引用: 3]

Linderholm H W, Björklund J, Seftigen K , et al.

Fennoscandia revisited: A spatially improved tree-ring reconstruction of summer temperatures for the last 900 years

[J]. Climate Dynamics, 2015,45(3/4):933-947.

[本文引用: 3]

Matskovsky V V, Helama S .

Testing long-term summer temperature reconstruction based on maximum density chronologies obtained by reanalysis of tree-ring data sets from northernmost Sweden and Finland

[J]. Climate of the Past, 2014,10(4):1473-1487.

[本文引用: 5]

Büntgen U, Frank D, Grudd H , et al.

Long-term summer temperature variations in the Pyrenees

[J]. Climate Dynamics, 2008,31(6):615-631.

[本文引用: 3]

Dorado Liñán I, Büntgen U, González-Rouco F , et al.

Estimating 750 years of temperature variations and uncertainties in the Pyrenees by tree-ring reconstructions and climate simulations

[J]. Climate of the Past, 2012,8(3):919-933.

[本文引用: 3]

Büntgen U, Krusic P J, Verstege A , et al.

New tree-ring evidence from the Pyrenees reveals Western Mediterranean climate variability since Medieval Times

[J]. Journal of Climate, 2017,30(14):5295-5318.

[本文引用: 5]

Klippel L, Krusic P J, Konter O , et al.

A 1200 + year reconstruction of temperature extremes for the northeastern Mediterranean region

[J]. International Journal of Climatology, 2019,39(4):2336-2350.

[本文引用: 2]

Rydval M, Gunnarson B E, Loader N J , et al.

Spatial reconstruction of Scottish summer temperatures from tree rings

[J]. International Journal of Climatology, 2017,37(3):1540-1556.

[本文引用: 2]

D'Arrigo R, Mashig E, Frank D, et al.

Reconstructed warm season temperatures for Nome, Seward Peninsula, Alaska

[J]. Geophysical Research Letters, 2004,31(9):L09202. doi: 10.1029/2004gl019756.

[本文引用: 3]

Jacoby G C, D'Arrigo R D.

Tree-ring width and density evidence of climatic and potential forest change in Alaska

[J]. Global Biogeochemical Cycles, 1995,9(2):227-234.

[本文引用: 3]

Wiles G C, D'Arrigo R D, Jacoby G C.

Temperature changes along the gulf of Alaska and the Pacific Northwest coast modeled from coastal tree rings

[J]. Canadian Journal of Forest Research, 1996,26(3):474-481.

[本文引用: 3]

Anchukaitis K J, D'Arrigo R D, Andreu-Hayles L, et al.

Tree-ring-reconstructed summer temperatures from Northwestern North America during the last nine centuries

[J]. Journal of Climate, 2013,26(10):3001-3012.

[本文引用: 4]

O'Donnell A J, Allen K J, Evans R M, et al.

Wood density provides new opportunities for reconstructing past temperature variability from southeastern Australian trees

[J]. Global and Planetary Change, 2016,141:1-11.

[本文引用: 3]

Björklund J A, Gunnarson B E, Seftigen K , et al.

Blue intensity and density from northern Fennoscandian tree rings, exploring the potential to improve summer temperature reconstructions with earlywood information

[J]. Climate of the Past, 2014,10(2):877-885.

[本文引用: 1]

Briffa K R, Osborn T J, Schweingruber F H , et al.

Tree-ring width and density data around the Northern Hemisphere: Part 2, spatio-temporal variability and associated climate patterns

[J]. The Holocene, 2002,12(6):759-789.

[本文引用: 1]

Davi N K, Jacoby G C, Wiles G C .

Boreal temperature variability inferred from maximum latewood density and tree-ring width data, Wrangell Mountain region, Alaska

[J]. Quaternary Research, 2003,60(3):252-262.

[本文引用: 2]

陈峰, 王慧琴, 袁玉江 , .

树轮最大密度记录的吉尔吉斯斯坦天山山区公元1650年以来的7—8月温度变化

[J]. 沙漠与绿洲气象, 2014,8(4):1-7.

[本文引用: 2]

[ Chen Feng, Wang Huiqin, Yuan Yujiang , et al.

July-August temperature variability since A.D. 1650 inferred from tree-ring maximum density of Picea schrenkiana in the Tianshan Mountains, Kyrgyzstan

Desert and Oasis Meteorology, 2014,8(4):1-7. ]

[本文引用: 2]

Chen F, Yuan Y J, Wei W S , et al.

Tree ring density-based summer temperature reconstruction for Zajsan Lake area, East Kazakhstan

[J]. International Journal of Climatology, 2012,32(7):1089-1097.

[本文引用: 2]

Davi N, D'Arrigo R, Jacoby G , et al.

Warm-season annual to decadal temperature variability for Hokkaido, Japan, inferred from maximum latewood density (AD 1557-1990) and ring width data (AD 1532-1990)

[J]. Climatic Change, 2002,52(1/2):201-217.

[本文引用: 2]

Li M Y, Duan J P, Wang L , et al.

Late summer temperature reconstruction based on tree-ring density for Sygera Mountain, southeastern Tibetan Plateau

[J]. Global and Planetary Change, 2018,163:10-17.

[本文引用: 3]

Liang H X, Lyu L X, Wahab M .

A 382-year reconstruction of August mean minimum temperature from tree-ring maximum latewood density on the southeastern Tibetan Plateau, China

[J]. Dendrochronologia, 2016,37:1-8.

[本文引用: 2]

Li M Q, Huang L, Yin Z Y , et al.

Temperature reconstruction and volcanic eruption signal from tree-ring width and maximum latewood density over the past 304 years in the southeastern Tibetan Plateau

[J]. International Journal of Biometeorology, 2017,61(11):2021-2032.

DOI:10.1007/s00484-017-1395-0      URL     PMID:28685206      [本文引用: 5]

This study presents a 304-year mean July-October maximum temperature reconstruction for the southeastern Tibetan Plateau based on both tree-ring width and maximum latewood density data. The reconstruction explained 58% of the variance in July-October maximum temperature during the calibration period (1958-2005). On the decadal scale, we identified two prominent cold periods during AD 1801-1833 and 1961-2003 and two prominent warm periods during AD 1730-1800 and 1928-1960, which are consistent with other reconstructions from the nearby region. Based on the reconstructed temperature series and volcanic eruption chronology, we found that most extreme cold years were in good agreement with major volcanic eruptions, such as 1816 after the Tambora eruption in 1815. Also, clusters of volcanic eruptions probably made the 1810s the coldest decade in the past 300 years. Our results indicated that fingerprints of major volcanic eruptions can be found in the reconstructed temperature records, while the responses of regional climate to these eruption events varied in space and time in the southeastern Tibetan Plateau.

Li M Y, Wang L, Fan Z X , et al.

Tree-ring density inferred late summer temperature variability over the past three centuries in the Gaoligong Mountains, southeastern Tibetan Plateau

[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2015,422:57-64.

[本文引用: 1]

Yu S L, Yuan Y J, Wei W S , et al.

A 352-year record of summer temperature reconstruction in the western Tianshan Mountains, China, as deduced from tree-ring density

[J]. Quaternary Research, 2013,80(2):158-166.

[本文引用: 2]

陈峰, 袁玉江, 魏文寿 , .

树轮记录的贺兰山北部近208年5—7月温度变化

[J]. 应用气象学报, 2011,22(4):463-471.

[本文引用: 2]

[ Chen Feng, Yuan Yujiang, Wei Wenshou , et al.

May-July temperature variability since 1801 inferred from tree rings of Pinus tabulaeformis of Helan Mountains in China

Journal of Applied Meteorological Science, 2011,22(4):463-471. ]

[本文引用: 2]

Sun Y, Wang L L, Chen J , et al.

Reconstructing mean maximum temperatures of May-August from tree-ring maximum density in North Da Hinggan Mountains, China

[J]. Chinese Science Bulletin, 2012,57(16):2007-2014.

[本文引用: 3]

陈峰, 袁玉江, 魏文寿 , .

利用树轮密度重建新疆北部5—8月温度变化

[J]. 冰川冻土, 2017,39(1):43-53.

[本文引用: 4]

[ Chen Feng, Yuan Yujiang, Wei Wenshou , et al.

Air temperature from May through August in northern Xinjiang reconstructed from multi-site tree-ring density

Journal of Glaciology and Geocryology, 2017,39(1):43-53. ]

[本文引用: 4]

Xing P, Zhang Q B, Lv L X .

Absence of late-summer warming trend over the past two and half centuries on the eastern Tibetan Plateau

[J]. Global and Planetary Change, 2014,123:27-35.

[本文引用: 2]

Chen F, Yuan Y .

May-June maximum temperature reconstruction from mean earlywood density in north Central China and its linkages to the summer monsoon activities

[J]. PLoS ONE, 2014,9(9):e107501. doi: 10.1371/journal.pone.0107501.

URL     PMID:25207554      [本文引用: 2]

杨银科, 黄强, 刘禹 , .

以树木年轮密度资料重建鄂尔多斯中部地区6至10月降水量的变化

[J]. 西北农林科技大学学报(自然科学版), 2013,41(8):96-102, 109.

[本文引用: 2]

[ Yang Yinke, Huang Qiang, Liu Yu , et al.

Tree-ring density based precipitation reconstruction from June to October in Central Ordos. Journal of Northwest A&

F University (Natural Science Edition), 2013,41(8):96-102, 109. ]

[本文引用: 2]

Wang L, Duan J P, Chen J , et al.

Temperature reconstruction from tree-ring maximum density of Balfour spruce in eastern Tibet, China

[J]. International Journal of Climatology, 2010,30(7):972-979.

[本文引用: 1]

Duan J P, Zhang Q B .

A 449 year warm season temperature reconstruction in the southeastern Tibetan Plateau and its relation to solar activity

[J]. Journal of Geophysical Research: Atmospheres, 2014,119(20):11578-11592.

[本文引用: 2]

吴普, 王丽丽, 邵雪梅 .

采用高山松最大密度重建川西高原近百年夏季气温

[J]. 地理学报, 2005,60(6):998-1006.

[本文引用: 1]

[ Wu Pu, Wang Lili, Shao Xuemei .

Reconstruction of summer temperature from maximum latewood density of Pinus densata in west Sichuan

Acta Geographica Sinica, 2005,60(6):998-1006. ]

[本文引用: 1]

段建平, 王丽丽, 李论 , .

树轮最大密度记录的贡嘎山区公元1837年以来的温度变化

[J]. 科学通报, 2010,55(11):1036-1042.

[本文引用: 1]

[ Duan Jianping, Wang Lili, Li Lun , et a1.

Temperature variability since A.D. 1837 inferred from tree-ring maximum density of Abies fabric in Gongga Mountains, China

Chinese Science Bulletin, 2010,55(11):1036-1042. ]

[本文引用: 1]

Fan Z X, Bräuning A, Yang B , et al.

Tree ring density-based summer temperature reconstruction for the central Hengduan Mountains in Southern China

[J]. Global and Planetary Change, 2009,65(1/2):1-11.

[本文引用: 1]

陈津, 王丽丽, 朱海峰 , .

用天山雪岭云杉年轮最大密度重建新疆伊犁地区春夏季平均最高温度变化

[J]. 科学通报, 2009,54(9):1295-1302.

[本文引用: 1]

[ Chen Jin, Wang Lili, Zhu Haifeng, et a1.

Reconstructing mean maximum temperature of growing season from the maximum density of the Sehrenk Spruce in Yili, Xinjiang, China

Chinese Science Bulletin, 2009,54(9):1295-1302. ]

[本文引用: 1]

陈峰, 袁玉江, 魏文寿 , .

树轮记录的过去384a乌鲁木齐河源7月温度变化

[J]. 冰川冻土, 2011,33(1):55-63.

[本文引用: 2]

[ Chen Feng, Yuan Yujiang, Wei Wenshou , et al.

July temperature at the upper treeline recorded in the tree-ring in the headwaters of the Urumqi River

Journal of Glaciology and Geocryology, 2011,33(1):55-63. ]

[本文引用: 2]

陈峰, 袁玉江, 魏文寿 , .

用西伯利亚落叶松年轮最大密度重建和布克赛尔5—8月份平均最高温度

[J]. 生态学报, 2010,30(17):4652-4658.

[本文引用: 1]

[ Chen Feng, Yuan Yujiang, Wei Wenshou , et al.

Dendroclimatic reconstruction of mean maximum May-August temperature from the maximum density of the Larix sibirica in Hoboksar, China

Acta Ecologica Sinica, 2010,30(17):4652-4658. ]

[本文引用: 1]

孙毓, 王丽丽, 陈津 .

长白落叶松树轮生长对气候变化的响应以及对夏季气温的重建

[J]. 地球环境学报, 2012,3(3):889-899.

[本文引用: 2]

[ Sun Yu, Wang Lili, Chen Jin .

Response of tree growth to climate change and reconstruction of summer temperature based on Korean larch

Journal of Earth Environment, 2012,3(3):889-899. ]

[本文引用: 2]

Jones P D, Briffa K R, Schweingruber F H .

Tree-ring evidence of the widespread effects of explosive volcanic-eruptions

[J]. Geophysical Research Letters, 1995,22(11):1333-1336.

[本文引用: 2]

Schneider L, Smerdon J E, Pretis F , et al.

Corrigendum: A new archive of large volcanic events over the past millennium derived from reconstructed summer temperatures

[J]. Environmental Research Letters, 2017,12(11):119501. doi: 10.1088/1748-9326/aa9426.

[本文引用: 1]

Briffa K R, Jones P D, Schweingruber F H , et al.

Influence of volcanic eruptions on Northern Hemisphere summer temperature over the past 600 years

[J]. Nature, 1998,393:450-455.

[本文引用: 1]

李明启, 邵雪梅 .

基于树轮资料初探过去千年强火山喷发与青藏高原东部温度变化关系

[J]. 地球科学进展, 2016,31(6):634-642.

[本文引用: 1]

[ Li Mingqi, Shao Xuemei .

Study on the relationship between large volcanic eruptions and temperature variation based on tree-ring data in the eastern Tibetan Plateau during the past millennium

Advances in Earth Science, 2016,31(6):634-642. ]

[本文引用: 1]

Esper J, Schneider L, Smerdon J E , et al.

Signals and memory in tree-ring width and density data

[J]. Dendrochronologia, 2015,35:62-70.

[本文引用: 1]

Briffa K R, Jones P D, Schweingruber F H. Tree-ring density reconstructions of summer temperature patterns across Western North America since 1600[J].

Journal of Climate

1992,5(7):735-754.

[本文引用: 1]

Duan J P, Li L, Ma Z G , et al.

Summer cooling driven by large volcanic eruptions over the Tibetan Plateau

[J]. Journal of Climate, 2018,31(24):9869-9879.

[本文引用: 1]

Mann M E, Fuentes J D, Rutherford S .

Underestimation of volcanic cooling in tree-ring-based reconstructions of hemispheric temperatures

[J]. Nature Geoscience, 2012,5(3):202-205.

[本文引用: 2]

Esper J, Büntgen U, Luterbacher J , et al.

Testing the hypojournal of post-volcanic missing rings in temperature sensitive dendrochronological data

[J]. Dendrochronologia, 2013,31(3):216-222.

[本文引用: 1]

Tingley M P, Stine A R, Huybers P .

Temperature reconstructions from tree-ring densities overestimate volcanic cooling

[J]. Geophysical Research Letters, 2014,41(22):7838-7845.

[本文引用: 2]

Stoffel M, Khodri M, Corona C , et al.

Estimates of volcanic-induced cooling in the Northern Hemisphere over the past 1500 years

[J]. Nature Geoscience, 2015,8(10):784-788.

[本文引用: 1]

Schneider L, Smerdon J E, Büntgen U , et al.

Revising midlatitude summer temperatures back to A.D. 600 based on a wood density network

[J]. Geophysical Research Letters, 2015,42(11):4556-4562.

[本文引用: 2]

D'Arrigo R, Jacoby G, Buckley B, et al.

Tree growth and inferred temperature variability at the North American Arctic treeline

[J]. Global and Planetary Change, 2009,65(1/2):71-82.

[本文引用: 2]

Briffa K R, Schweingruber F H, Jones P D , et al.

Reduced sensitivity of recent tree-growth to temperature at high northern latitudes

[J]. Nature, 1998,391:678-682.

[本文引用: 1]

/