文章信息
- 张瑞波, 袁玉江, 喻树龙, 陈峰, 张同文, 尚华明, 范子昂
- ZHANG Ruibo, YUAN Yujiang, YU Shulong, CHEN Feng, ZHANG Tongwen, SHANG Huaming, FAN Zi'ang.
- 树轮记录的青海过去300年5-6月平均最高气温时空变化
- The spatiotemporal variability of May-June maximum temperature in past 300 years on the Qinghai Plateau, according to tree ring records
- 生态学报[J]. 2016, 36(23): 7603-7613
- Acta Ecologica Sinica[J]. 2016, 36(23): 7603-7613
- http://dx.doi.org/10.5846/stxb201506241285
-
文章历史
- 收稿日期: 2015-06-24
- 修订日期: 2015-04-25
20世纪80年代开始, 全球变化科学引发了人类社会的广泛关注, 气候变化、形成原因及其影响已经成为地球科学中最为活跃的研究领域, 并取得重大进展。现代观测资料已经证实气候是不稳定的, 在不同时间尺度上存在各种变率。对人类社会具有迫切现实意义的是几十年到百年尺度上的变化, 关系到当代和今后几代人的生存环境[1-2]。20世纪气候变化的归因是当前全球关注的焦点之一, 而解决这一问题的途径之一是对过去千年的气候变化历史进行精确重建[3]。树木年轮由于具有空间分布广、时间序列长、分辨率高、定年准确、环境变化指示意义明确且可定量等优势而成为研究过去气候变化的首选代用资料之一, 在揭示气候变化规律及机理研究中发挥了重要作用[4-5]。近年来国际上著名的几条千年温度变化序列的建立都依靠了大量树轮资料[6-9]。树木年轮气候重建研究在国外已经从单点小区域的气候要素重建发展到空间上大尺度的气候要素变化的重建。特别是对于欧洲和北美地区, 由于树轮采样点较为密集, 资料共享程度高, 树轮资料丰富, 近年来用树轮重建大空间尺度气候变化的研究较多[10-14]。在亚洲, Ohyama等利用CRU格点资料和5个树轮年表重建了东亚地区春季温度气候场[15]。袁玉江等曾利用阿尔泰山和天山的多个树轮年表进行过较大地域范围的重建北疆500a干湿变化[16]。
现代气候学研究表明, 青藏高原是我国以及东亚气候变化的“敏感区”和“启动区”, 其气候变率和突变对我国气候变化有着十分重要的先兆意义。青海省位于青藏高原之东北部, 其气候变化具有明显的高原特色, 同时是我国西北干旱区的干旱中心, 是西部干旱区、东部季风区和青藏高原区三大区域的交汇地带, 是全球变化的敏感区和生态系统脆弱区[3]; 揭示青海省气候变化的规律, 对研究青藏高原及周边地区生态环境的演变及制定气候变化影响的对策等方面具有现实意义。过去, 对青海大量研究集中于利用祁连圆柏树轮重建过去千年气候序列[17-30], 科研人员已经在该地区建立了迄今为止中国最长的长达3585a的树轮年表[31-33], 重建了过去2485年温度状况并与历史朝代兴衰做了对比[34-35], 该曲线具有一定的全球代表性[36]。但是大量的研究都是基于单点的树轮资料, 很少有人利用整个青海不同气候区树轮多条序列重建气候场, 了解过去整个青海地区气候的时空变化特征。利用树木年轮资料重建大范围的气候场是当今国际树木年轮气候研究的新动向和热点问题, 也是我国树轮研究几近空白的薄弱环节。本研究利用青海地区不同地理单元的祁连圆柏树轮宽度年表重建该地区过去300年气候场, 分析青海过去300年气候变化的时间序列特征及变化规律。
1 资料与方法 1.1 气象资料的选取用于研究气候变化的资料不仅要有足够的序列长度, 还应有足够的精确度, 青海多数气象台站建立于20世纪50年代中后期, 建站至1960年, 不少台站存在较多的缺测、漏测和观测中断现象, 根据均一性、代表性和比较性原则, 本文选定青海省内30个站作为代表站(茶卡、达日、大柴旦、德令哈、都兰、刚察、格尔木、贵德、久治、冷湖、玛多、茫崖、门源、民和、囊谦、诺木洪、祁连、共和、清水河、曲麻莱、同仁、托勒、沱沱河、五道梁、西宁、小灶火、兴海、野牛沟、玉树、杂多), 取各站1961-2009年的平均气温、平均最高气温、平均最低气温和降水量资料, 并采用Mann-Kendall方法, 来检查站址迁移是否引起观测资料的不连续(图 1)。
空间代表性分析气温资料来自英国East Anglia大学的Climatic Research Unit (简称CRU) 1901年1月到2009年12月的高分辨率全球逐月格点数据集(CRU TS 3.1), 其网格距均为0.5°×0.5°。CRU数据集有以下优点:在资料的重建过程中包含了严格的时间均一性检验, 时间尺度更长, 空间分辨率更高, 该数据集被大量运用于气候变化研究中[37]。
1.2 树轮年表的建立根据青海祁连圆柏的分布和区域代表性, 研究组分别于2007年和2009年两次在青海柴达木盆地、祁连山、青南高原和青海东部共采集了17个采样点676棵树1367个祁连圆柏样芯(图 1, 表 1), 采样点海拔在3000-4300 m之间。将采集的树轮样本经过干燥、固定、打磨、交叉定年和轮宽测量; 用COFECHA定年质量控制程序进行交叉定年的检验; 采用winARSTAN年表研制程序完成树轮宽度年表的建立, 在年表研制过程中, 分别采用样条函数方法、负指数方法和区域曲线标准化方法进行去除生长趋势过程, 比较发现, 样条函数方法对拟合的序列对气候的响应较好, 因此本研究最终采用样条函数拟合的年表进行气候分析, 并剔除代表性较差或异常的序列, 最终建立青海12个树轮宽度年表(图 2), 其分布于青海不同重要地理单元, 如柴达木盆地(QRG、DUL、BSS、ELS、KKS)、青南高原(ASS、ASX)、祁连山地区(DLG、QYG)、青海东部(HBX、MAX、SIG)。
采点名 Sampling Location |
代号 Code |
北纬(N) Latitude |
东经(E) Longitude |
海拔/m Elevation |
坡向 Aspect |
坡度 Slope |
郁闭度 Canopy density |
大柳沟 | DLG | 38°16′ | 100°00′ | 3350 | S | 30° | 0.2 |
青羊沟 | QYG | 38°10′ | 100°26′ | 3200 | E | 20° | 0.2 |
柏树山 | BSS | 37°09′ | 97°20′ | 3500 | 0° | 0 | |
鄂拉山 | ELS | 36°23′ | 98°20′ | 3560 | NW | 40° | 0 |
曲日岗 | QRB | 36°01′ | 98°11′ | 3750 | S | 30° | 0.15 |
柯柯寺 | KKS | 37°19′ | 98°24′ | 3390 | S | 5° | 0 |
都兰 | DUL | 36°02′ | 98°13′ | 3724 | SE | 37° | 0.13 |
寺沟 | SIG | 37°08′ | 102°28′ | 3023 | SW | 36° | 0.43 |
麦秀 | MAX | 35°16′ | 101°53′ | 3214 | SW | 38° | 0.62 |
河北乡 | HBX | 34°47′ | 100°49′ | 3418 | NE | 38° | 0.57 |
昂赛上线 | ASS | 32°42′ | 95°37′ | 4249 | S | 35° | 0.47 |
昂赛下线 | ASX | 32°43′ | 95°37′ | 4047 | S | 43° | 0.64 |
采用DendroClim2002相关分析和响应分析解释树轮径向生长对气候的响应; 空间分析方法选取经验正交函数法(EOF)[38], 它能够把随时间变化的气象要素场分解为空间函数部分和时间函数(主分量)部分。空间函数部分概括场的地域分布特点, 这部分是不随时间变化的; 而时间函数部分则由空间点(变量)的线性组合而构成, 称为主分量, 这些主分量的头几个占有原空间点(变量)的总方差的很大部分。本研究采用着重表现空间的相关性分布结构的旋转经验正交函数(Rotated Empirical Orthogonal Function, REOF)方法分析树轮场对气温场的代表性, 采用了多元回归方法重建气候空间分布状况, 采用区域平均序列的逐步回归方法重建气候时间序列; 采用“逐一剔除法”[39]进行交叉检验; 利用功率谱分析方法对重建序列进行周期检验。
2 结果与分析 2.1 青海气温场和树轮宽度场特征将12个年表和30个气象站平均气温、平均最高气温、平均最低气温相关普查表明, 整个青海地区, 5-6月平均最高气温与12个年表相关均较好, 因此利用树轮重建的时段选取5-6月平均最高气温。
首先分析由30个气象站1961-2008年5-6月平均最高气温构成青海5-6月的平均最高气温场。REOF空间分析表明(表 2, 图 3), 青海气温的同步性较好, 反应其同步变化的第一特征向量占气温场方差总量的62.06%, 且均为正值。而第二特征向量仅占总方差的15.33%, 前3个特征向量对总方差的贡献率达到81.62%且特征值均大于1, 其后各特征向量贡献率很小, 所以青海地区5-6月平均最高气温具有较快的收敛速度, 前三个特征向量可以较好的说明温度场的空间分布特征。经验正交第一模态(图 3)表明, 整个青海地区的主成分载荷都为正值, 在0.13-0.22之间, 相差很小, 说明整个青海地区5-6月温度空间变化高度一致, 虽然青海地区下垫面地形复杂, 温差较大, 但是各地变化趋势是非常一致的。高载荷量地区分布在青海北部, 以祁连山和柴达木盆地为主。青南高原载荷量最低。
主要特征参数 Main parameters |
模态 Principal component |
特征值 Eigenvector value |
方差贡献率/% Variance contribution |
累计贡献率/% Accumulative variance contribution |
气温场Temperature field | 1 | 18.617 | 62.06 | 62.06 |
2 | 4.599 | 15.33 | 77.39 | |
3 | 1.268 | 4.23 | 81.62 | |
树轮场Tree-ring field | 1 | 4.737 | 39.47 | 39.47 |
2 | 1.720 | 14.33 | 53.81 | |
3 | 1.372 | 11.43 | 65.24 | |
4 | 1.134 | 9.45 | 74.70 |
将12个年表REOF展开结果显示, 第一特征向量均为正值, 树轮宽度场中各采点变化趋势基本一致, 但树轮宽度场收敛速度远小于气温场, 第一特征向量对总方差的贡献率只有39.47%, 第二、第三、第四特征向量分别占总方差的14.33%、11.43%和9.45%, 前4特征向量累计贡献率为74.70%, 这是因为树木年轮生长同时受到多个气候因子的影响, 虽然平均最高气温是最大的影响因子, 但其他气候因子仍然有较显著的影响, 从而在树轮宽度场中表现为较慢的收敛速度[40], 按照特征值大于1作为截留主分量标准, 本文用前4特征向量代表宽度场。
气温场和树轮宽度场第一特征向量相关系数为-0.465(P < 0.01), 两场的第一特征向量均表现为同步变化一致的特点。气温场和树轮宽度场第一特征向量高值中心位于青海北部的祁连山区和柴达木盆地, 而低值中心位于青南高原西南部和东南部(阿尼玛卿山)。
2.2 气温场的时间序列重建及检验将30个气象站5-6月最高气温进行平均, 得到青海的区域平均最高气温序列, 利用逐步剔除回归法分别重建各气象站和区域平均最高温度序列, 考虑到年表样本总体代表性(EPS), 以子样本信号强度大于0.85的标准选取年表分析长度, 最终共获得1700-2006年气温重建序列。其中5个年表进入青海地区平均最高气温序列的重建方程:
(1) |
式中,T56为青海地区5-6月平均最高气温的区域平均序列(图 4), ASS、ASX、BSS、ELS、QYG分别代表昂赛上线、昂赛下线、柏树山、额拉山、青羊沟树轮宽度年表。该方程复相关系数达到0.764, 重建值对实测值的方差解释量分别为58.3%, RAdj2=0.531, F检验值F5, 40=11.2, 超过0.0001的显著性水平。
重建效果较好的气象站分布在祁连山西部、柴达木盆地和青南高原地区, 大部分复相关系数在0.7以上, 较差的气象站分布在青海东部。这是因为青海东部海拔不高、并处于东南季风边缘区, 受暖湿气流影响, 因此该区域树木径向生长对气候的响应复杂, 温度不是限制性因子。一般认为, 高寒地区和森林上线树木径向生长对气温较为敏感, 而青南高原和祁连山西部温度低, 柴达木盆地虽然海拔较低, 但是其5-6月平均最高气温低于青海东部, 海拔较青海东部高。
R | r | S1 | S2 | RE | t |
0.684 | 0.611 | 17/45 | 12/46 | 0.463 | 4.748 |
R:复相关系数, r:逐一剔除后相关系数, S1:一阶差符号检验, S2:符号检验, RE:误差缩减值, t为乘积平均数 |
方程稳定性检验值如表, 剔除相关和一阶差相关均超过0.0001的显著水平, S1未通过检验, 而S2通过0.01的显著性水平, 说明相对于高频变化, 重建序列在低频变化上同实测序列吻合得更好些。t值通过0.001的显著性水平, RE值大于0.3, 可见重建方程是较为稳定的。区域平均序列重建值与近45年实测值一致性较好(图 5), 说明重建的过去300年青海地区5-6月平均最高气温是较为可靠。
2.3 过去300年青海5-6月平均最高气温变化及代表性分析为了解1700-2006年青海地区平均最高气温的冷暖阶段变化, 将青海地区5-6月平均最高气温区域平均重建序列进行11a滑动平均(图 4)。青海地区过去300年平均最高气温时间序列与喻树龙等利用川西14个树轮年表和14个气象站合成的川西高原过去6-7月平均最低气温场时间序列高度一致(图 4)[40], 表明整个青藏高原东北部过去300多年温度变化可能较为一致。过去300年青海气温大致可分为5冷5暖的变化阶段(表 4), 最冷的时段为19世纪30-40年代, 温度距平平均偏低0.22℃, 最长的偏冷期为19世纪末20世纪初的40多年, 最暖的时段都发生在20世纪30-50年代, 最长的偏暖期为18世纪末19世纪初的40a。20世纪60年代以来, 青海5-6月平均最高气温持续上升, 尤其是80年代到现在, 青海地区平均最高气温呈现急剧持续上升现象。器测资料研究表明, 柴达木盆地年平均气温在1985年前后发生了突变, 1987-2010年柴达木盆地各地的平均气温显著上升[41]。施雅风等研究表明[42], 我国西北地区1987-2000年年均温比1961-1986年高0.7℃。这些研究结果均与本文研究结果一致。20世纪初和60-80年代为偏冷期, 30-50年代和90年代至今均为偏暖期与秦宁生等[43]用树轮宽度重建的6月最高气温, Wu等[44]用最大密度年表重建的近百年6-9月平均气温的冷暖阶段完全一致。这些均不但表明本文重建结果具有一定的可靠性, 而且说明, 整个青藏高原东北部温度过去300年温度变化的一致性。
年份 Years |
年数 Number of years |
距平 Anomalies |
1706-1716 | 11 | -0.11 |
1755-1791 | 37 | -0.16 |
1832-1848 | 17 | -0.22 |
1885-1928 | 44 | -0.08 |
1960-1992 | 33 | -0.19 |
1717-1754 | 38 | +0.39 |
1792-1831 | 40 | +0.26 |
1849-1884 | 36 | +0.22 |
1929-1959 | 31 | +0.48 |
1993-2001 | 9 | +0.29 |
利用功率谱分析方法对青海5-6月平均最高气温时间序列进行周期分析, 以了解过去300年青海地区5-6月平均最高气温时间演变的周期性, 取最大滞后M=102, 相当于各序列长度的三分之一。分析发现(图 6), 过去300年青海地区5-6月平均最高气温在0.05显著性水平上具有2.1、3.1、8.5a、和68.0a的变化准周期。青海温度场重建序列具有2.1a和8.5a变化周期特征与“准两年脉动(QBO)”十分接近, 已有研究表明, QBO的影响存在于较大的范围, 其可能与海气间相互耦合振荡有关[45]。3.1a准周期与厄尔尼诺(ENSO)的周期3-5a一致, ENSO现象具有全球性, 说明气温的变化可能与海气间相互耦合振荡有关, 短周期振荡受大尺度海气相互作用的全球事件(ENSO)的影响[46-47]。25.5a对应太阳黑子活动22a准周期, 说明太阳黑子的活动对这一区域的温度变化也有影响[48]。68a的变化准周期与全球气候系统海洋-大气的50-88年周期震荡一致[49]。这些周期暗示了青海过去300年气温受当地气候变化影响外, 还受到更大范围气候变化的影响, 进一步证明了气候变化的全球性。
青海气候受多种因素的影响, 其中大气环流是控制该地区气候特征的主要因素。相关分析表明, 青海5-6月平均最高气温重建序列与NAO、AO及主要环流指数有较好的一致性(图 7), 与NAO (-0.198, P < 0.05)和AO (-0.184, P < 0.05)呈显著负相关关系, 与5-6月大西洋欧洲环流型W相关达到-0.533, 显著性水平超过0.0001, 还与5-6月北半球极涡强度指数、大西洋欧洲区极涡强度指数、北美区极涡强度指数、北半球极涡强度指数、北美大西洋副高北界、大西洋欧洲区极涡面积指数、西藏高原A相关通过了0.01显著性检验, 这些环流指数都与西风环流有着密切关系, 表明青海地区温度变化受西风影响显著。这是因为青海地区位于31°30-39°19′N, 根据500hPa高空等压面分析, 青海处在西风带之内, 10月至翌年6月, 高原面对西风环流产生分支、绕爬流、屏障等动力作用, 扩大了西风带的影响范围, 因此西风环流成为影响和控制青海气候的重要因子之一。另外, 青海重建序列还与4月份的南极涛动(-0.306, P < 0.05)、印缅槽指数(0.397, P < 0.01)以及南方涛动指数(-0.184, P < 0.05)呈负相关关系, 表明该地区温度变化还受到印度季风的影响。青海地区温度和太平洋和亚洲的环流指数相关并不显著, 说明太平洋气流及东亚季风对青海地区气候影响微弱。
空间相关分析表明(图 8), 1901-2006年100多年相关显示, 重建区域平均序列和青海CRU相关大部分在0.4以上, 而1950-2006年相关显示, 重建区域平均序列和青海CRU相关大部分在0.5以上。说明青海5-6月区域平均最该气温重建序列可较好的代表整个青海地区5-6月平均最高气温。另外, 无论是近100年还是近50年的空间相关均表明, 青海区域平均序列与青藏高原大部分地区(包括川西高原)及南亚次大陆5-6月平均最高气温显著相关, 青海气候场重建序列在一定程度上可代表西南季风区5-6月平均最高气温。
3 结论(1) 气温场和树轮宽度场第一特征向量相关系数为-0.465(P < 0.01), 两场的第一特征向量均表现为同步变化一致的特点:气温场和树轮宽度场第一特征向量高值中心位于青海北部的祁连山区和柴达木盆地, 而低值中心位于青南高原西南部和东南部(阿尼玛卿山)。
(2) 过去300年青海气温大致可分为5冷5暖的变化阶段, 最冷的时段为1830s-1840s年代, 最长的偏冷期为19世纪末20世纪初, 最暖的时段都发生在1930s-1950s年代, 最长的偏暖期为18世纪末19世纪初, 存在5个明显的持续增温时段和4个持续降温时段, 增温缓慢, 降温迅速。20世纪60年代以来, 青海5-6月平均最高气温持续上升, 尤其是80年代到现在, 青海地区平均最高气温呈现急剧持续上升。过去300年青海地区5-6月平均最高气温具有2.1、3.1、8.5、25.5a和68.0a的变化准周期。
(3) 青海5-6月平均最高气温受西风影响和印度季风显著, 太平洋及东亚季风对青海地区气候影响微弱; 青海气候场重建序列在一定程度上可代表青藏高原大部分地区甚至西南季风区5-6月平均最高气温。
致谢: 本研究使用的青海祁连圆柏树轮年表均是在中国科学院地理科学与资源研究所邵雪梅研究员和中国科学院青藏高原研究所朱海峰副研究员的指导下完成, 中国科学院地理科学与资源研究所王丽丽研究员参加了野外采样工作, 特此致谢。[1] | 叶笃正, 符淙斌. 全球变化的主要科学问题. 大气科学 , 1994, 18 (4) : 498–512. |
[2] | 秦大河. 中国西部环境演变评估(综合卷):中国西部环境演变评估综合报告. 北京: 科学出版社, 20028–9. |
[3] | 杨保, 史锋, SonechkinD M, 王章勇, 秦春. 过去千年气候变化重建研究新进展. 中国沙漠 , 2011, 31 (2) : 485–491. |
[4] | Graumlich L J. A 1000-year record of temperature and precipitation in the Sierra Nevada. Quaternary Research , 1993, 39 (2) : 249–255. DOI:10.1006/qres.1993.1029 |
[5] | Mann M E, Bradley R S, Hughes M K. Global-scale temperature patterns and climate forcing over the past six centuries. Nature , 1998, 392 (6678) : 779–787. DOI:10.1038/33859 |
[6] | Esper J, Cook E R, Schweingruber F H. Low-frequency signals in long tree-ring chronologies for reconstructing past temperature variability. Science , 2002, 295 (5563) : 2250–2253. DOI:10.1126/science.1066208 |
[7] | Mann M E, Bradley R S, Hughes M K. Northern hemisphere temperatures during the past millennium:inferences, uncertainties, and limitations. Geography Research Letter , 1999, 26 (6) : 759–762. DOI:10.1029/1999GL900070 |
[8] | Mann M E. Climate over the past two millennia. Annual Reviews of Earth and Planetary Sciences , 2007, 35 (1) : 111–136. DOI:10.1146/annurev.earth.35.031306.140042 |
[9] | Büntgen U, Tegel W, Nicolussi K, McCormick M, Frank D, Trouet V, Kaplan J O, Herzig F, Heussner K U, Wanner H, Luterbacher J, Esper J. 2500 Years of European climate variability and human susceptibility. Science , 2011, 331 (6017) : 578–582. DOI:10.1126/science.1197175 |
[10] | Briffa K R, Osborn T J, Schweingruber F H. Large-scale temperature inferences from tree rings:a review. Global and Planetary Change , 2004, 40 (1/2) : 11–26. |
[11] | Cook E R, Esper J, D'Arrigo R D. Extra-tropical Northern Hemisphere land temperature variability over the past 1000 years. Quaternary Science Reviews , 2004, 23 (20/22) : 2063–2074. |
[12] | Jones P D, Briffa K R, Barnett T P, Tett S F B. High-resolution Palaeoclimatic records for the last millennium:interpretation, integration and comparison with General Circulation Model control-run temperatures. The Holocene , 1998, 8 (4) : 455–471. DOI:10.1191/095968398667194956 |
[13] | Mann M E, Zhang Z H, Rutherford S, Bradley R S, Hughes M K, Shindell D, Ammann C, Faluvegi G, Ni F. Global signatures and dynamical origins of the little ice age and medieval climate anomaly. Science , 2009, 326 (5957) : 1256–1260. DOI:10.1126/science.1177303 |
[14] | Cook E R, Anchukaitis K J, Buckley B M, D'Arrigo R D, Jacoby G C, Wright W E. Asian monsoon failure and megadrought during the last millennium. Science , 2010, 328 (5977) : 486–489. DOI:10.1126/science.1185188 |
[15] | Ohyama M, Yonenobu H, Choi J N, Park W K, Hanzawa M, Suzuki M. Climate field reconstruction of East Asian spring temperature. Climate of the Past Discussions , 2012, 8 (4) : 3535–3549. DOI:10.5194/cpd-8-3535-2012 |
[16] | 袁玉江, 韩淑媞. 北疆500年干湿变化特征. 冰川冻土 , 1991, 13 (4) : 315–322. |
[17] | 杨保. 树轮记录的小冰期以来青藏高原气候变化的时空特征. 第四纪研究 , 2012, 32 (1) : 81–94. |
[18] | 朱海峰, 郑永宏, 邵雪梅, 刘晓宏, 徐岩, 梁尔源. 树木年轮记录的青海乌兰地区近千年温度变化. 科学通报 , 2008, 53 (15) : 1835–1841. |
[19] | Liu Y, An Z S, Ma H Z, Cai Q F, Liu Z Y, Kutzbach J K, Shi J F, Song H M, Sun J Y, Yi L, Li Q, Yang Y K, Wang L. Precipitation variation in the northeastern Tibetan Plateau recorded by the tree rings since 850 AD and its relevance to the Northern Hemisphere temperature. Science in China:Series D Earth Sciences , 2006, 49 (4) : 408–420. DOI:10.1007/s11430-006-0408-3 |
[20] | Bräuning A. Dendrochronology for the last 1400 years in eastern Tibet. GeoJournal , 1994, 34 (1) : 75–95. DOI:10.1007/BF00813972 |
[21] | Bräuning A. Climate history of the Tibetan Plateau during the last 1000 years derived from a network of Juniper chronologies. Dendrochronologia , 2001, 19 (1) : 127–137. |
[22] | Shao X M, Huang L, Liu H B, Liang E Y, Fang X Q, Wang L L. Reconstruction of precipitation variation from tree rings in recent 1000 years in Delingha, Qinghai. Science in China Series D:Earth Sciences , 2005, 48 (7) : 939–949. DOI:10.1360/03yd0146 |
[23] | Yang B, Bräuning A, Johnson K R, Shi Y F. General characteristics of temperature variation in China during the last two millennia. Geophysical Research Letters , 2002, 29 (9) : 38-1–38-4. DOI:10.1029/2001GL014485 |
[24] | Zhang Q B, Cheng G D, Yao T D, Kang X C, Huang J G. A 2, 326-year tree-ring record of climate variability on the northeastern Qinghai-Tibetan Plateau. Geophysical Research Letters , 2003, 30 (14) : 1739. |
[25] | Zhang Q B, Qiu H Y. A millennium-long tree-ring chronology of Sabina przewalskii on northeastern Qinghai-Tibetan Plateau. Dendrochronologia , 2007, 24 (2/3) : 91–95. |
[26] | Zhang Y, Gou X H, Chen F H, Tian Q H, Yang M L, Peng J F, Fang K Y. A 1232 years tree-ring record of climate variability in the Qilian Mountains, Northwestern China. IAWA Journal , 2009, 30 (4) : 407–420. DOI:10.1163/22941932-90000228 |
[27] | 康兴成, GraumlichL J, SheooardP. 青海都兰地区1835年来的气候变化-来自树轮资料. 第四纪研究 , 1997 (1) : 70–75. |
[28] | Sheppard P R, Tarasov P E, Graumlich L J, Heussner H U, Wagner M, Österle H, Thompson L G. Annual precipitation since 515 BC reconstructed from living and fossil juniper growth of northeastern Qinghai Province, China. Climate Dynamics , 2004, 23 (7/8) : 869–881. |
[29] | Yang B, Qin C, Wang J L, He M H, Melvin T M, Osborn T J, Briffa K R. A 3500-year tree-ring record of annual precipitation on the northeastern Tibetan Plateau. PNAS , 2014, 111 (8) : 2903–2908. DOI:10.1073/pnas.1319238111 |
[30] | Liu X H, Qin D H, Shao X M, Chen T, Ren J W. Temperature variations recovered from tree-rings in the middle Qilian Mountain over the last millennium. Science in China Series D:Earth Sciences , 2005, 48 (4) : 521–529. DOI:10.1360/03yd0063 |
[31] | 邵雪梅, 王树芝, 徐岩, 朱海峰, 许新国, 肖永民. 柴达木盆地东北部3500年树轮定年年表的初步建立. 第四纪研究 , 2007, 27 (4) : 477–485. |
[32] | Shao X M, Wang S Z, Zhu H F, Xu Y, Liang E Y, Yin Z Y, Xu X G, Xiao Y M. A 3585-year master ring-width dating chronology of Qilian Juniper from the Northeastern Qinghai-Tibetan Plateau. IAWA Journal , 2009, 30 (4) : 379–394. DOI:10.1163/22941932-90000226 |
[33] | Shao X M, Xu Y, Yin Z Y, Liang E, Zhu H, Wang S. Climatic implications of a 3585-year tree-ring width chronology from the northeastern Qinghai-Tibetan Plateau. Quaternary Science Reviews , 2010, 29 (17/18) : 2111–2122. |
[34] | Liu Y, An Z S, Linderholm H W, Chen D L, Song H M, Cai Q F, Sun J Y, Tian H. Annual temperatures during the last 2485 years in the mid-eastern Tibetan Plateau inferred from tree rings. Science in China Series D:Earth Sciences , 2009, 52 (3) : 348–359. DOI:10.1007/s11430-009-0025-z |
[35] | Liu Y, Cai Q F, Song H M, An Z S, Linderholm H W. Amplitudes, rates, periodicities and causes of temperature variations in the past 2485 years and future trends over the central-eastern Tibetan Plateau. Chinese Science Bulletin , 2011, 56 (28-29) : 2986–2994. DOI:10.1007/s11434-011-4713-7 |
[36] | 刘禹, 蔡秋芳, 宋慧明. 关于青藏高原2485年温度的季节和空间代表性问题. 第四纪研究 , 2013, 33 (1) : 108–114. |
[37] | Xu H M, Taylor R G, Kingston D G, Jiang T, Thompson J R, Todd M C. Hydrological modeling of river Xiangxi using SWAT2005:a comparison of model parameterizations using station and gridded meteorological observations. Quaternary International , 2010, 226 (1/2) : 54–59. |
[38] | Wilks D. Statistical Methods in the Atmospheric Sciences. 2nd ed. Amsterdam, The Netherlands:Elsevier, 2006:1857. |
[39] | 吴祥定. 树木年轮与气候变化. 北京: 气象出版社, 1990171–172. |
[40] | 喻树龙, 袁玉江, 魏文寿, 尚华明, 张同文, 陈峰, 张瑞波. 川西高原6-7月最低气温场重建研究. 中国沙漠 , 2012, 32 (4) : 1010–1016. |
[41] | 戴升, 申红艳, 李林, 王振宇, 肖建设. 柴达木盆地气候由暖干向暖湿转型的变化特征分析. 高原气象 , 2013, 32 (1) : 211–220. |
[42] | Shi Y F, Shen Y P, Kang E, Li D L, Ding Y J, Zhang G W, Hu R J. Recent and future climate change in northwest China. Climate Change , 2007, 80 (3/4) : 379–393. |
[43] | 秦宁生, 时兴合, 邵雪梅, 汪青春. 川西高原树木年轮所指示的平均最高气温变化. 高原山地气象研究 , 2008, 28 (4) : 18–24. |
[44] | Wu P, Wang L L, Shao X M. Reconstruction of summer temperature variation from maximum density of alpine pine during 1917-2002 for west Sichuan Plateau, China. Journal of Geographical Sciences , 2008, 18 (2) : 201–210. DOI:10.1007/s11442-008-0201-7 |
[45] | 钱维宏, 朱亚芬, 叶谦. 赤道东太平洋海温异常的年际和年代际变率. 科学通报 , 1998, 43 (10) : 1098–1102. |
[46] | Rittenour T M, Brigham-Grette J, Mann M E. El Niño-like climate teleconnections in new England during the late Pleistocene. Science , 2000, 288 (5468) : 1039–1042. DOI:10.1126/science.288.5468.1039 |
[47] | Huber M, Caballero R. Eocene El Niño:evidence for robust tropical dynamics in the "Hothouse". Science , 2003, 299 (5608) : 877–881. DOI:10.1126/science.1078766 |
[48] | Rasmusson E M, Wang X L, Ropelewski C F. The biennial component of ENSO variability. Journal of Marine Systems , 1990, 1 (1/2) : 71–96. |
[49] | Schlesinger M E, Ramankutty N. An oscillation in the global climate system of period 65-70 years. Nature , 1994, 367 (6465) : 723–726. DOI:10.1038/367723a0 |