1960-2014年祁连山南坡及其附近地区降水时空变化特征

付建新， 曹广超， 李玲琴， 曹生奎， 唐仲霞，杨晓敏， 蒋 刚， 虞 敏， 袁 杰， 刁二龙

(1.青海师范大学 青海省自然地理与环境过程重点实验室, 西宁 810008; 2.青海师范大学 地理科学学院, 西宁 810008)

IPCC报告表明全球变暖已经是不容争辩的事实，全球变暖对降水有可能产生影响[1]。秦大河等[2]认为全球变暖可能造成中纬度地区降水增加，前人对于祁连山区的降水研究主要从时间尺度与空间尺度进行了分析。时间尺度包括小时、季节、年等角度。小时尺度：殷雪莲等[3]对祁连山区2003年7月25日与2005年6月30日两次降水过程从环流背景、动力机制以及地形特征等方面做了研究。季节尺度：刘雪梅等[4]对祁连山区夏季降水的日变化特征进行了分析；陈乾等[5]对祁连山区2006年夏季降水过程从3类气流型进行了分析。年尺度：强芳等[6]基于格点数据对祁连山区1961—2012年面雨量的特征进行了分析；魏锋等[7]使用小波变换方法分析了祁连山地区17个测站1960—2004年5—9月逐日降水量的特征；张小明等[8]对祁连山近45 a异常降水的时空特征进行了研究；张存杰等[9]采用EOF和REOF等方法,研究了近40 a来祁连山附近气温和降水的时空分布特征。空间尺度包括流域、分段、整个区域。流域尺度：蓝永超等[10]研究了黑河山区汇流区降水对全球变暖的响应；牛赟等[11]对祁连山大野口流域降水特征进行了分析。分段尺度：王宁练等[12]对祁连山中段北坡最大降水高度带进行了分析。整个区域：贾文雄[13]对近50 a来祁连山及河西走廊降水的时空变化进行了研究；汤懋苍[14]对祁连山区降水的地理分布特征进行了阐述。对降水成因的研究主要从经纬度、地形、大气环流等角度。陈少勇等[15]认为祁连山主体的降水多于周围地区，祁连山东部降水比西部稳定；李岩瑛等[16]分析了不同降水强度的时空分布特征及其与海拔的关系；贾文雄[17]认为祁连山的气温与地理位置的关系存在相关性关系，降水与地理位置的关系较为复杂；程鹏等[18]分析了高空西风急流对祁连山区降水的影响。也有学者从降水与祁连山区水源涵养的关系角度进行了阐述。常学向等[19]研究了祁连山林区森林对降水的截留作用；张学龙等[20]分析了祁连山青海云杉林截留对降水的分配效应。前人对于祁连山区降水主要集中在时空视角、成因分析与对水源涵养的影响等方面，在方法上使用传统与信息技术结合的较少，本文主要对祁连山南坡附近的降水特征进行分析，使用传统气象分析方法与GIS与SPPSS软件制图分析方法并用进行研究。

祁连山区的生态是国家当前与今后重点保护与建设的对象，该区域是“丝绸之路经济带”重要的水源涵养功能区[21]，所以在新形势下，研究祁连山南坡及其附近地区降水的变化具有重要的理论与现实意义。

1 研究区概况

本文研究区概况见图1，核心区位于青海省境内，但是核心区气象站点较少，所以增加了外围区域的气象站点，外围区为核心区的南北两侧，分别位于青、甘两省。祁连山南坡，海拔2 257～5 235 m，是我国重要的西北—东南走向的一系列山脉，是我国重要的地理分界线之一[22]，具有高原大陆性气候特征[23]。祁连山南坡是青海省重要的水源涵养保护区，气候的变化直接会对当地植被等生态环境造成影响，通过研究祁连山南坡及其附近地区降水的变化为进一步研究祁连山南坡水源涵养区的生态系统等研究做铺垫。

图1祁连山南坡及其附近地区气象站点分布

2 数据来源与研究方法

本文数据来源于中国气象数据网(http:∥data.cma.cn/)，选取1960—2014年祁连山南坡及其附近地区19个气象站点的降水日数、降水量、降水强度数据，为了保证研究结果的准确性，所以对数据进行一致性检验的验证，对区域内平均降水序列与各站点降水的原始序列的相关性进行了检测，相关系数均在0.793以上，说明降水数据表现出较高的相关性，所以数据不仅能反映各站点的气候变化，同时也能反映出区域内的气候变化；降水日数为日降水量≥0.1 mm日数的总和，降水强度为各等级降水总量与其对应降水日数之比。研究方法包括线性趋势法、相关分析法、多项式趋势法、5 a滑动平均，Mann-Kendall与滑动t检验[24-25]法，利用ArcGIS 10.0对降水日数与降水强度的空间变化做制图分析，对于文中显著性检验，使用SPSS统计学软件中的相关性分析模块，对时间序列与原序列变量之间的相关系数进行不同程度的显著性检验(p<0.05，p<0.01)。

3 降水的时间序列变化

3.1 多年降水日数与降水强度年际变化趋势

祁连山区的降水受到3种气流的影响：其一来自太平洋的东南季风带来的湿润水汽，其二受到来自大西洋冷湿气流的影响，夏季还要受到翻越高大青藏高原之后的印度洋暖湿气团的影响[26]。

祁连山南坡及其附近地区1960—2014年多年降水日数(图2A)为缓慢增长趋势(0.030 d/a)，其平均值为77.67 d，降水日数大于平均值的年份占49.09%，小于平均值的年份占50.91%；降水日数最多的为93.42 d(1988年)，最少的为64.89 d(1997年)，二者相差28.53 d；90 d以上占5.45%，80 d以上占32.73%；从5 a滑动平均来看，多年降水日数20世纪90年代以后较以前变化大，大致分为2个大周期(20世纪80年代中后期—21世纪初期、21世纪初期以来)与3个小周期(20世纪60年代中期—70年代中期、70年代中期—80年代初期、80年代初期—后期)。降水强度(图2B)整体上为缓慢波动增长趋势(0.009 mm/a)，降水强度的平均值为3.247 mm/a，52.73%的年份降水强度超过了均值，47.27%的年份降水强度低于均值，最大值为3.79mm/a(2014年)，最小值为2.63 mm/a(1962年)二者相差1.16 mm/a；从5 a滑动平均来看，其波动较降水日数小，期间存在2个大的周期(20世纪60年代中期—70年代中期、90年代初期—21世纪初期)，同时还并存若干小的周期。

图2祁连山南坡及其附近地区多年平均降水日数与降水强度时间变化

3.2 降水日数与降水强度的季节变化趋势

从四季降水的线性回归图(图3)得出，春季(3月—5月)、夏季(6月—8月)、秋季(9月—11月)、冬季(12月—翌年2月)的降水日数倾斜率分别为：-0.000 2，0.002 3，-0.000 04，0.008 2 d/a，除了冬季降水日数为增加趋势，其余季节变化趋势不大，四季的降水日数的增加趋势均小于多年平均降水日数(0.034 9 d/a)；降水强度倾斜率分别为：0.005 5 mm/a(p<0.01)，0.013 6 mm/a(p<0.01)，0.007 2 mm/a(p<0.05)，0.005 mm/a(p<0.01)，均通过了显著性检验，四季降水强度均为增加趋势，上升的速率由大到小依次为：夏、秋、春、冬，夏季降水强度的增长率高于多年平均降水强度的增长率(0.008 8 mm/a)，说明夏季降水强度对全年降水贡献最大。

春季降水日数(图3A)中最大值为8.404 d(1967年)，最小值为3.737 d(2013年)，二者相差4.667 d，其中大于平均值(5.949 d)的占到47.27%，小于平均值的占到52.73%；5 a滑动平均显示3个升降周期：20世纪60年代中期—70年代末、70年代末—21世纪初期、21世纪初期—研究期末，其中第2个周期最长，包括了1个显著上升区(20世纪70年代末—90年代初)。夏季降水日数(图3C)整体变化与春季相似，变化不明显，不同的是夏季略微上升，春季略微下降；最大值14.754 d(1983年)与最小值9.702 d(1962年)相差5.052 d，其中大于平均值(11.908 d)的占到45.45%，；20世纪80年代中期为1个降水日数高峰期。秋季降水日数(图3E)变化趋势基本不变，平均值为6.383 d，大于与小于平均值所占比重大致相当，基本上在平均值上下波动，20世纪70年代中期—90年代中期为1个较长下降区间，20世纪90年代中期—21世纪10年代初期为1个较长的上升区间。冬季降水日数(图3G)的变化趋势较其他季节的变化大，呈现增长趋势(0.008 2 d/a)，2008年降水日数最多(5.035 d)，比最少降水日数(1.456 d)多3.579 d，小于平均值的比重比大于平均值的比重多16.36%，5 a滑动平均表现4个升降区间。

春季降水强度(图3B)为上升趋势(0.005 5 mm/a)，春季降水强度的平均值为2.521 mm/a，41.82%的年份变化趋势大于均值，58.18%的年份变化趋势小于均值，最大值3.799 mm/a(1985年)与最小值差1.440 mm/a(1995年)差值较大，为2.359 mm/a；5 a平均滑动显示两个大的周期，即20世纪60年代中期—80年代中期与20世纪80年代中期—21世纪初。夏季降水强度(图3D)的倾斜率(0.013 6 mm/a)比春季、秋季、冬季降水强度的倾斜率分别高0.008，0.006，0.009 mm/a；夏季降水强度范围为3.343～5.014 mm/a，平均值为4.298 mm/a，大于均值的比重为52.73%，小于均值的比重为47.27%，大致从20世纪90年代降水强度开始大于均值；5 a滑动平均表现出5个周期，其中包括了多个小的波动。秋季降水强度(图3F)的增长率为0.072 mm/a，平均值为3.196 mm/a，期间包括1个大的周期与若干个小的周期，大的周期为20世纪80年代—90年代中期。冬季降水强度(图3H)的增长率(0.005 mm/a)是四季中最低的，其均值为0.899 mm，大于均值的年份与小于均值的年份之比为5∶6，20世纪60年代中期—70年代中期为低于均值的1个较长周期，21世纪开始了新的1个较长周期。

图3祁连山南坡及其附近地区降水日数与降水强度的季节变化

3.3 季节降水量的突变检验

春季降水量(图4A)：正序列UF与反序列UB的交点是1982年，交点且位于±1.96临界线之间，所以1982年为突变开始的年份，这与贾文雄等[27]对祁连山东部季节降水突变检验的结论一致，1986年降水量增加趋势明显，但是之后UF与UB基本没有超过临界线，UF只在1993年超过了临界线，所以突变并不明显；夏季降水量(图4B)：正序列曲线UF的值于1973年大于0，之后降水量波动上升，期间UF与UB有有效的交点为3个，1972年、2008年、2010年，结合滑动t检验，当n1=n2=5，10时，在a=0.01的显著性检验，可以确定1972年为突变年份开始年份，1987—1991年、1993—2001年、2006—2009年、2010—2014年为突变时间段，一直到2010年突变明显加快；秋季降水量(图4C)：突变开始的年份为2003年，但是正序列曲线UF一直没有超过临界线，突变不明显；冬季降水量(图4D)：正序列UF与反序列UB的交点共有4个，为1989年、1993年、1996年、1998年，结合滑动t检验，当n1=n2=5，9时，在a=0.05的显著性检验，确定为1998年为突变开始的年份，从1990年降水增多，尤其在2005年之后，增加趋势迅速。

图4祁连山南坡及其附近地区四季降水量突变检验

3.4 降水日数与降水强度的年代际变化趋势

年代际年平均降水日数(图5)由大到小依次为：81.03 d(20世纪80年代)，79.97 d(2000—2014年)，77.01 d(20世纪70年代)，76.89 d(20世纪60年代)，74.70 d(20世纪90年代)，20世纪80年代与2000—2014年高于年降水日数平均值(77.67 d)；20世纪80年代增长最快，为1.071 d/a，20世纪90年代为负增长(-0.326 d/a)，二者相差1.397 d/a，年代际增长率波动较大。20世纪60，70，80，90年代、2000-2014年的平均降水强度分别3.055 mm，3.219 mm，3.144 mm，3.308 mm，3.448 mm，可以看出60—80年代降水强度较低，90年代之后降水强度增加明显，尤其进入21世纪降水强度增加愈加突出，这与中国60—80年代中期降水偏少，80年代中后期增多的趋势相似[28]；降水强度增长率最高的是0.381 mm/10 a(20世纪90年代)，增长率最低的是0.074 mm/10 a(20世纪70年代)，大致在0～0.4 mm/10 a。

4 降水的空间变化

4.1 多年平均降水日数与降水量空间分布

从图6中可以看出，多年平均降水量与平均降水日数均呈现带状分布，大致为西北—东南走向，由东南向西北递减，主要因为夏季风东南季风对此区域影响由东南向西北减弱；研究区多年平均降水日数与降水量的均值分别为80.76 d与267.79 mm，降水日数排在前4位即降水日数大于100 d的站点为乌鞘岭、门源、野牛沟、祁连、刚察，这些站点均位于东南部，降水量大于400 mm的站点为门源、野牛沟、祁连、乌鞘岭，这4个站点同样均位于东南部；位于西北部的高台、酒泉降水日数大约分别为45，40 d，降水量大约分别为109，88 mm；祁连山南坡与北坡的降水量之比为4∶1，南坡多年降水量最多与最少的站点分别为门源(520.86 mm)、德令哈(180.39 mm)，北坡多年降水量最多与最少的站点分别为永昌(204.29 mm)、酒泉(87.51 mm)；南坡多年降水日数最多与最少的站点分别为乌鞘岭(136.53 d/a)、德令哈(52.53 d/a)，北坡多年降水日数最多与最少的站点分别为永昌(7.82 d/a)、民勤(39.46 d/a)，南坡降水明显多于北坡，因为祁连山南坡为东南季风的迎风坡，来自西北太平洋的暖湿气流顺着湟水谷地沿着祁连山南坡不断爬升，受到地形的抬升作用，在一定海拔高度降水量增加，南坡形成了多雨区，北坡处于东南季风的背风坡，水汽受到山地阻挡，成为雨影区，气候类型为温带大陆性气候，海洋水汽影响弱，降水量少。

图5祁连山南坡及其附近地区降水日数与降水强度的年代际变化

图6祁连山南坡及其附近地区多年平均降水量与平均降水日数空间分布

4.2 降水日数与降水强度年际变化趋势空间分布

从(图7A)可知，1960—2014年各站点的年降水日数为正增长的站点数量为10个，占到全部站点52.63%，为负增长的站点数量为9个，占到全部站点47.37%，正负增长之比大体相当；增长最快与最慢的站点分别为德令哈(0.365 d/a)与托勒(0.004 d/a)，下降最快与最慢的站点分别为门源(-0.218 d/a)与西宁(-0.012 d/a)；核心研究区倾斜率的平均值为-0.09 d/a，呈下降趋势，外围研究区倾斜率的平均值为0.051 d/a，呈上升趋势；祁连山南坡站点倾斜率平均值为0.026 d/a，北坡站点倾斜率平均值为0.020 d/a，南北坡降水日数基本相当；只有3个站点通过了显著性检验，德令哈为p<0.01，门源、乌鞘岭为p<0.05。

由图7B可知，1960—2014年各站点的年降水强度除了民和(-0.023 mm/10 a)与恰卜恰(-0.002 mm/10 a)为负增长外，其余站点均为正增长，占全部站点89.47%；增长最快与最慢的站点分别为都兰(0.207 mm/10 a)与张掖(0.037 mm/10 a)，二者相差0.170 mm/10 a；倾斜率>0.1 mm/10 a的站点数量为8个，占到全部站点42.11%，增长率<0.1 mm/10 a的站点数量为11个，占到全部站点57.89%；核心研究区3个站点倾斜率>0.1 mm/10 a，野牛沟(0.155 mm/10 a)、托勒(0.154 mm/10 a)、祁连(0.117 mm/10 a)，倾斜率平均值为0.120 mm/10 a，外围研究区倾斜率平均值为0.079 mm/10 a，核心研究区的增长速度是外围研究区的1.524倍，增长较快；祁连山南坡站点倾斜率平均值为0.100 1 mm/10 a，北坡站点倾斜率平均值为0.066 mm/10 a，所以南坡增长速度大约是北坡的1.513倍；63.16%的站点通过显著性检验，其中p<0.01的站点为5个，占全部站点26.32%，占通过检验站点41.67%；p<0.05的站点为7个，占全部站点36.84%，占通过检验站点58.33%。

图7祁连山南坡及其附近地区降水日数与降水强度年际变化空间分布

4.3 降水日数与降水强度的季节变化空间分布

各站点的春季降水日数(图8A)正负增长之比为11∶8，分别占到全部站点57.89%与42.11%；正增长的范围为0.002～0.062 d/a，负增长的范围为-0.105 8～-0.000 4 d/a，增长最快与最慢的站点分别为德令哈与高台，下降最快与最慢的站点分别为门源与张掖，增长最快与下降最快的站点与年际变化的站点一致；核心研究区倾斜率的平均值为-0.049 d/a，呈下降趋势，外围研究区倾斜率的平均值为0.008 d/a，呈略微上升趋势；祁连山南坡站点倾斜率平均值为-0.008 d/a，北坡站点倾斜率平均值为0.003 d/a，南北坡降水日数变化趋势差别不大；门源、都兰、刚察3个站点通过了显著性检验(p<0.05)。夏季降水日数(图8C)正负增长的站点数量基本相当，增长最快的站点0.203 d/a(德令哈)与下降最快的站点-0.112 d/a(门源)相差0.315 d/a；核心研究区野牛沟、祁连、门源呈下降趋势，托勒呈增长趋势，整体上呈下降趋势(-0.025 d/a)；外围研究区整体上呈上升趋势(0.023 d/a)；祁连山南坡站点整体上呈上升趋势(0.029 d/a)，北坡整体上呈下降趋势(-0.016 d/a)；通过显著性检验的站点占到全部站点21.05%，其中德令哈通过(p<0.01)的检验。秋季降水日数(图8E)为正增长的站点数量为10个，负增长的为9个，整体上为略微下降趋势(-0.005 d/a)；变化最大的两个站点为乌鞘岭(0.079 d/a)与野牛沟(-0.072 d/a)，二者相差0.151 d/a；只有乌鞘岭通过了(p<0.05)相关性检验。冬季降水日数(图8G)的正增长、零增长、负增长三者之比为15∶1∶3，整体上表现出正增长(0.024 d/a)，增长最快的是乌鞘岭(0.108 d/a)，降低最快的是-0.034 d/a，二者相差0.074 d/a；核心研究区整体上为增长趋势(0.014 d/a)，托勒是唯一零增长的站点；3个站点通过了显著性检验。

各站点的春季(图8B)、夏季(图8D)、秋季(图8F)、冬季(图8H)(以下皆为此顺序)降水强度增长率之比分别为15∶4，18∶1，13∶6，18∶1，增长率最高的站点分别为都兰(0.251 mm/10 a)、永昌(0.252 mm/10 a)、乌鞘岭(0.629 mm/10 a)、高台(0.155 mm/10 a)，下降最快的站点分别为茶卡(-0.074 mm/10 a)、民和(-0.102 mm/10 a)、酒泉(-0.095 mm/10 a)、托勒(-0.003 mm/10 a)；核心研究区的平均增长率分别为0.069，0.171，0.096，0.027 mm/10 a；通过显著性检验的站点数量依次为0，3，2，10；其中通过p<0.01检验的站点数量为1(秋季)，6(冬季)，春、夏均未通过此检验，通过p<0.05检验的站点数量为3(夏季)，1(秋季)，4(冬季)，春季未通过此检验。

4.4 降水日数与降水强度的年代际变化空间分布

20世纪60年代年均降水日数(图9A)84.21%的站点为正增长(0.188～1.546 d/a)，15.79%的站点为负增长(-0.418～-0.273 d/a)，整体上表现出正增长(均值为0.708 d/a)；武威与门源成为增长最快与下降最快的的站点，二者相差1.127 d/a；核心研究区整体上为下降趋势(均值为-0.197 d/a)；其中武威与高台通过了p<0.05的相关性检验。20世纪60年代年均降水强度(图9B)13个站点为正增长(0.127～1.238 mm/10 a)，6个站点为负增长(-0.288～-0.234 mm/10 a)，整体呈现正增长(0.300 mm/10 a)；永昌与民和分别为增长最快与最慢的站点，门源与德令哈为下降最快与最慢的站点；核心研究区整体上呈上升趋势(均值为0.171 mm/10 a)；通过显著性检验的站点为祁连(p<0.01)、永昌(p<0.01)、都兰(p<0.05)。20世纪70年代年均降水日数(图9C)正负增长之比为14∶5，整体上呈正增长(均值为0.578 d/a)；核心研究区整体上呈略微下降趋势(均值为-0.001 5 d/a)，外围研究区呈增长趋势(均值为0.733 d/a)。20世纪70年代年均降水强度(图9D)整体上表现为上升趋势(均值为0.066 mm/10 a)，核心研究区的倾斜率为0.326 mm/10 a，比全部站点的倾斜率大0.260 mm/10 a，外围区的倾斜率为-0.003 1 mm/10 a，呈现略微下降的趋势；通过相关性检验的比例为21.05%。

图8祁连山南坡及其附近地区降水日数与降水强度季节变化空间分布

20世纪80年代年均降水日数(图9E)变化范围为-0.709～3.794 d/a，二者相差3.085 d/a，变化差异较大，整体上呈现增长趋势(均值为1.071 d/a)；核心研究区的倾斜率为0.641 d/a，外围研究区的倾斜率为1.186 d/a；正增长的站点数量为13，负增长的站点数量为6；刚察与托勒为增长最快与最慢的站点；刚察与野牛沟、茶卡、民和分别通过了p<0.01与p<0.05相关性检验。20世纪80年代年均降水强度(图9F)整体上呈增长趋势(0.249 mm/10 a)，南坡倾斜率为0.304 mm/10 a，北坡倾斜率为0.130 mm/10 a，南坡较北坡增加的快；正增长的比例为63.16%，负增长的比例为36.84%；只有门源通过了相关性检验(p<0.05)。20世纪90年代年均降水日数(图9G)68.42%的站点为负增长，下降最快的为民和(-1.606 d/a)，核心研究区整体上为负增长(均值为-0.406 d/a)，外围研究区也为负增长(均值为-0.253 d/a)。20世纪90年代年均降水强度(图9H)正负增长之比为12∶7，正增长的范围为0.015～1.383 mm/10 a，增长最快与最慢的站点分别为刚察与门源，负增长的范围为-0.910～-0.081 mm/10 a，下降最快与最慢的站点分别为酒泉与都兰；其中6个站点通过了相关性检验(p<0.05)。2000—2014年年均降水日数(图9I)的平均值为0.24 d/a，正增长范围为0.246～1.021 d/a，增长最快的为茶卡，民勤为零增长，正增长范围为-0.496～-0.068 d/a，门源下降最快；茶卡与恰卜恰通过了p<0.05的相关性检验。2000—2014年年均降水强度(图9J)正负增长之比为13∶6，平均值为0.133 mm/10 a，正增长的范围为0.037～0.812 mm/10 a，负增长的范围为-0.564～-0.050 mm/10 a；核心研究区的倾斜率为0.310 mm/10 a，外围研究区的倾斜率为0.086 mm/10 a；南坡的倾斜率为0.202 mm/10 a，北坡的倾斜率为0.016 mm/10 a；其中野牛沟通过了p<0.01的相关性检验。

图9祁连山南坡及其附近地区降水日数与降水强度年代际变化空间分布

5 讨 论

(1) 时间序列变化。多年降水日数与降水强度整体上表现为波动缓慢增长趋势；祁连山区的降水主要受到太平洋暖湿气流的影响，同时也受到来自大西洋冷湿气流的影响，降水年际变化大，夏季降水强度对全年贡献最大；20世纪80年代中后期之后较60—80年代中期的降水强度有明显增加，这种趋势与全国降水特征一致；季节突变减压发现除了秋季突变不明显，春、夏、冬发生突变的年份分别为1982年、1972年、1998年。

(2) 空间分布规律。多年平均降水量与平均降水日数受到东南季风的影响，均呈现出西北—东南走向的带状分布，各个站点降水分布呈现的走向与山脉走向一致，南坡为东南季风的迎风坡降水明显多于处在背风坡的北坡；影响降水的因素诸多，海拔、地形、迎风坡与背风坡、大气环流等；1960—2014年各站点的年降水日数为正负增长的站点数量之比大体相当，为10∶9，89.47%的站点年降水强度为正增长，这与张存杰等[29]研究结论基本一致；夏季降水强度增长最明显。

(3) 20世纪80年代中后期降水强度有所增加，尤其进入21世纪增加更为明显，这会导致河流径流量增加，比如黑河干流的径流量从20世纪偏枯状态到21世纪年径流量出现回升的现象[30]。

参考文献：

[1] Wentz F J, Ricciardulli L, Hilburn K, et al. How much more rain will global warming bring[J]. Science, 2007,317(5835):233-235.

[2] 秦大河, Thomas Stocker. IPCC第五次评估报告第一工作组报告的亮点结论[J].气候变化研究进展,2014,10(1):1-6.

[3] 殷雪莲,郭建华,董安祥,等.沿祁连山两次典型强降水天气个例对比分析[J].高原气象,2008,27(1):184-192.

[4] 刘雪梅,张明军,王圣杰,等.2008—2014年祁连山区夏季降水的日变化特征及其影响因素[J].地理学报,2016,71(5):754-767.

[5] 陈乾,张久林,陈添宇.祁连山区2006年夏季降水过程的天气分析[J].干旱气象,2008,26(4):22-28.

[6] 强芳,张明军,王圣杰,等.基于格点数据的1961—2012年祁连山面雨量特征分析[J].地理学报,2015,70(7):1125-1136.

[7] 魏锋,王劲松,李宝梓,等.祁连山近45 a5—9月日降水气候特征[J].干旱气象,2010,28(3):285-290.

[8] 张小明,魏锋,陆燕.祁连山近45 a年降水异常的气候特征[J].干旱气象,2006,24(3):35-41.

[9] 张存杰,郭妮.祁连山区近40年气候变化特征[J].气象,2002,28(12):33-39.

[10] 蓝永超,丁永建,康尔泗.近50年来黑河山区汇流区温度及降水变化趋势[J].高原气象,2004,23(5):723-727.

[11] 牛赟,刘贤德,敬文茂,等.祁连山大野口流域气温、降水、河川径流特征分析[J].干旱区地理:汉文版,2014,37(5):931-938.

[12] 王宁练,贺建桥,蒋熹,等.祁连山中段北坡最大降水高度带观测与研究[J].冰川冻土,2009,31(3):395-403.

[13] 贾文雄.近50年来祁连山及河西走廊降水的时空变化[J].地理学报,2012,67(5):631-644.

[14] 汤懋苍.祁连山区降水的地理分布特征[J].地理学报,1985(4):323-332.

[15] 陈少勇,董安祥,韩通.祁连山东、西部夏季降水量时空分布的差异及其成因研究[J].南京气象学院学报,2007,30(5):715-719.

[16] 李岩瑛,张强,许霞,等.祁连山及周边地区降水与地形的关系[J].冰川冻土,2010,32(1):52-61.

[17] 贾文雄.祁连山气候的空间差异与地理位置和地形的关系[J].干旱区研究,2010,27(4):607-615.

[18] 程鹏,陶健红,张新荣,等.高空西风急流对祁连山区降水的影响[J].气象科技,2007,35(4):489-494.

[19] 常学向,赵爱芬,王金叶,等.祁连山林区大气降水特征与森林对降水的截留作用[J].高原气象,2002,21(3):274-280.

[20] 张学龙,罗龙发,敬文茂,等.祁连山青海云杉林截留对降水的分配效应[J].山地学报,2007,25(6):678-683.

[21] 张强,杜志成.丝绸之路经济带上区域生态安全评价研究:以祁连山冰川与水涵养生态功能区为例[J].生态经济:中文版,2016,32(10):169-173.

[22] 王海军,张勃,靳晓华,等.基于GIS的祁连山区气温和降水的时空变化分析[J].中国沙漠,2009,29(6):1196-1202.

[23] 戴声佩,张勃.基于GIS的祁连山植被NDVI对气温降水的旬响应分析[J].生态环境学报,2010,19(1):140-145.

[24] 符淙斌,王强.气候突变的定义和检测方法[J].大气科学,1992,16(4):482-493.

[25] 魏凤英.现代气候统计诊断与预测技术[M].2版.北京:气象出版社,2007.

[26] 丁永建,叶柏生,周文娟.40年来西北干旱区黑河流域降水时空分布特征[J].冰川冻土,1999,21(1):42-48.

[27] 贾文雄,何元庆,李宗省,等.祁连山区气候变化的区域差异特征及突变分析[J].地理学报,2008,63(3):257-269.

[28] 秦大河,丁一汇,苏纪兰,等.中国气候与环境演变(上卷)[M].北京:科学出版社,2005.

[29] 张存杰,郭妮.祁连山区近40年气候变化特征[J].气象,2002, 28(12):33-39.

[30] 蓝永超,康尔泗,金会军,等.黑河出山径流量年际变化特征和趋势研究[J].冰川冻土,1999,21(1):49-53.