3.1 年内变化特征
结合张力样条函数法对省31个气象站点降水倾向率进行插值,在年降水量变化趋势空间分布(图2A)中可以看出,1960—2019年,除仁怀、余庆、都匀、黎平、平柱、惠水和榕江站点外,其余气象台站降水的气候倾向率均为负值,说明研究区整个区域存在普遍降水减少的趋势。降水量呈现增加趋势最快的站点主要有仁怀、平柱和榕江,其中平柱降水倾向率最大(30.84 mm/10 a),其次是榕江(20.01 mm/10 a); 降水呈减小趋势最快的站点主要有晋安、盘县、安顺和兴仁,其中盘县站点减少最多(-45.72 mm/10 a),其次是晋安(-42.70 mm/10 a)和安顺(-36.45 mm/10 a); 降水变化趋势不显著的站点主要是罗甸站点(-0.05 mm/10 a)。从省降水倾向率的空间分布可以看出,从西到东,年降水量逐渐增加; 其中六盘水、安顺及黔西南北部为降水减少的集中区,黔东南为降水增加的集中区; 省南部和北部降水呈减少趋势,中部黔南为降水较多的区域。
对省近60 a的年内降水量的气候倾向率进行统计(图2B),结果表明省区域范围内降水量呈增加趋势的月份主要有1月、3月、6月,趋势率分别为4.76 mm/10 a,3.59 mm/10 a和0.82 mm/10 a,其余月份均呈减小趋势,减少最快的主要是4月(-1.66 mm/10 a)和11月(-1.30 mm/10 a)。可见,省区域内降水减少主要是由于各月降水减少引起。
3.2 年际变化特征
图3A表明,1960—2019年省多年降水量波动幅度很大,但整体变化趋势不显著,趋势率仅为-9.25 mm/10 a(p>0.1)。近60 a来,最小年降水量出现在2011年,其年降水量比多年降水量少349 mm; 1966年、1989年、2009年的降水量也较少。年降水量最大值出现在1967年,比多年降水量多211 mm; 1976年、1977年、2015年降水也相对较多。期间出现3个较明显的波峰,即1969年(1 290 mm),1978年(1 300 mm)和1995年(1 282 mm); 出现2个明显的波谷,1987年(1 081 mm)和2011年(1 027 mm)。
对全省1960—2019年的降水距平百分率进行统计(图3B),结果表明,整体上降水量的年际变化分布并不均匀。降水负距平百分率显著偏大的年份有1966年、1981年、1989年、2005年、2009年、2011年、2013年,其中2011年降水负距平百分率最大(-29%),属于轻度干旱; 此外,1966年降水距平百分率为-17%,轻旱,而1967年降水距平百分率为18%,轻涝,可见1960s降水变化幅度很大,此间1966年发生了轻度干旱,即1966—1967年发生了旱涝交替的自然灾害现象。1984—1990年、2009—2013年出现连续降水偏少的现象、而1967—1971年、2014—2019年出现连续降水偏多的现象; 结合表1,发现近60 a来,贵州共发生7次轻度干旱事件,分别为1966年、1981年、1989年、2005年、2009年、2011年、2013年; 两次洪涝事件1967年和1977年。
3.3 季节时间变化
图4为1960—2019年贵州省季节降水量距平序列变化。总体上省4个季节的逐年降水量在近60 a来的变化趋势均不显著,其中春季和秋季降水量呈微弱减少趋势,下降速率为-5.41 mm/10 a(p>0.1)和-6.55 mm/10 a(p>0.1); 夏季和冬季降水量呈微弱增加的趋势,上升速率分别为1.60 mm/10 a(p>0.1)和2.07 mm/10 a(p>0.1)。此外,春季降水量在2011年达到最低,其数值比近60 a来的平均值低129.12 mm; 最高值出现在2016年,比多年均值高91.08 mm,可见近60 a来春季降水量的波动很大,振幅约为220 mm(图4A)。夏季降水量在1972年达到最低(比平均值低282.88 mm); 而最高值出现在1979年,比平均值高210.47 mm; 近60 a来夏季降水量振幅约为493 mm。秋季降水量振幅约为296 mm(图4C); 冬季降水量震荡上升(趋势不显著),极小值降水量为34.94 mm(2009年),极大值降水量为131.31 mm(2019年)(图4D)。
3.4 季节空间变化
对研究区31个气象站点四季降水量气候倾向率进行张力样条函数插值,并将站点降水倾向率以不同等级颜色表现出来(图5)。春季降水倾向率中,除仁怀、余庆、都匀、平柱和榕江外,其余站点降水倾向率均为负值,说明1960—2019年近60 a来,研究区春季存在普遍的降水减少趋势。其中从图5A中可以看出,负高值变化率以安顺站点(-18.85 mm/10 a)为中心,向外降水量呈增加趋势; 其次是以松桃站点(-15.29 mm/10 a)为负高值中心,向外降水量呈增加趋势; 降水量增加趋势最大的站点是榕江(7.00 mm/10 a),其次是余庆(6.18 mm/10 a),表明这两个区域降水增加多。春季降水量形成西南部及东北部降水较少,而东南部降水较丰富的空间分布格局。夏季降水趋势率空间分异较小,有11个气象站点降水趋势率为负值,其余站点均为正值,说明全省范围内,夏季降水整体以增加的趋势为主。其中,负高值以盘县站点为中心(-18.90 mm/10 a),正高值以平柱站点为中心(26.14 mm/10 a)。从图5B中可以看出夏季西南部降水呈减少趋势,东部降水呈现增加趋势,地区降水趋势差异不大。
图5 1960-2019年贵州省降水倾向率空间变化
秋季气象站点除仁怀、望漠、罗甸和榕江外,其余站点降水倾向率均为负值(图5C); 其中呈增加趋势的站点倾向率最大的为仁怀站点(7.69 mm/10 a),形成以仁怀为中心向外降水逐渐减少的空间分布,其余站点均在1.00 mm/10 a以下; 呈下降趋势的站点中趋势最大的站点是晋安站点(-22.28 mm/10 a),其次是盘县(-21.15 mm/10 a); 从整体分布格局上看,秋季西南部降水少,其余地区较多。冬季除威宁、晋安、桐梓、毕节和兴仁外,其余站点均为正值,说明研究区冬季降水存在普遍的增加趋势,图5D中可以看出低值以晋安站点(-2.45 mm/10 a)为中心,向外降水呈逐渐增加的趋势; 较为特殊的是冬季省区域范围内还出现以都匀(5.79 mm/10 a)和平柱(6.57 mm/10 a)为高值的降水中心点,向外呈环状延伸递减。
3.5 不同地貌区年际降水变化
喀斯特地貌是由于地下地表水对可溶性岩石进行侵蚀和堆积而成,由于地下地表水作用的时限长短及程度不同,不同时期形成不同的地貌类型; 对1960—2019年全省区域逐年降水分别进行张力样条函数空间插值,获得面状年降水空间分布。借鉴前人基于不同地貌类型区[32]对降水变化的影响分析,对贵州省不同地貌成因类型及构成进行统计分析(表2),并对不同地貌类型区降水量进行统计以了解降水变化趋势(表3)。从地貌类型特征分区统计结果可知,喀斯特低中山型区域成因以溶蚀—侵蚀类型为主,占比达81.0%; 峰丛谷地型区域以溶蚀为主型居多,占比81.3%; 峰林盆地(溶原)型区域以溶蚀为主型居多,占比89.0%; 峰丛洼地型以溶蚀为主型居多,占比64.9%; 峰林谷地型以溶蚀—侵蚀类型为主,占比53.6%; 混合型地貌区以侵蚀—剥蚀类型为主,占比64.7%。
从表3可知,近60 a,六大地貌类型区年降水均呈下降趋势。参照中国根据成因、形态组合及发育程度来划分中国岩溶七大类型的依据,从地貌成因类型的角度结合所占百分比来说:(1)溶蚀作用强烈的地区,地表岩溶形态组合以峰林地形及溶蚀洼地为主要特征; 从表2可知,以溶蚀作用为主的地貌类型区主要有峰丛谷地型、峰林盆地(溶原)型和峰丛洼地型,从表3年降水量统计可以看出,3种地貌类型中以峰丛洼地型降水最多(1 237.9 mm/a)且以-19.30 mm/10 a(p>0.01)的下降速率显著下降,峰林盆地(溶原)型地貌区降水量比峰丛谷地型区域多约27 mm/a。(2)溶蚀—侵蚀作用强烈的地区,以侵蚀作用形成的地貌为主,从表2及表3可以看出研究区地貌类型主要有喀斯特低中山型及峰林谷地型; 喀斯特低中山型比峰林谷地型地貌区年降水量少约24 mm/a,且降水倾向率为-12.13 mm/10 a(p>0.1),但不显著。此外,降水下降趋势较慢的为峰林谷地型地貌区,下降速率仅为-2.84 mm/10 a(p>0.1)。(3)以侵蚀—剥蚀为主的地貌区是混合型,年降水量较多(1 178.7 mm),年降水变化以-3.25 mm/10 a不显著下降。
表3 1960-2019年不同地貌类型年降水趋势及降水量
3.6 不同地貌区四季降水变化
由表4可知,1960—2019年,贵州省区域四季降水在不同的地貌类型差异较大。在6种地貌类型中,春季降水变化趋势均为负值; 其中峰丛洼地型地貌区降水量以-6.99 mm/10 a(p>0.01)的下降速率在显著减少,其次是喀斯特低中山地貌区降水最少(-5.79 mm/10 a),但下降速率不显著。与春季相似,秋季降水在6个地貌类型区均呈下降趋势,且峰丛洼地地貌区下降趋势明显且最大(-9.23 mm/10 a(p>0.01)),其次是峰林盆地(溶原)地貌类型区(-7.22 mm/10 a(p>0.1))和喀斯特低中山(-7.00 mm/10 a(p>0.01))区域降水呈下降趋势。夏季除峰丛洼地型地貌区降水呈下降趋势外,其余地貌类型区均呈上升趋势,其中以峰林谷地型降水上升趋势率最大(4.96 mm/10 a(p>0.1)),但显著性均不明显。冬季时,峰丛洼地型地貌区降水也以-0.94 mm/10 a的速率呈不显著下降趋势,降水增加地貌区中,以混合型地貌区为最大,上升趋势率为0.93 mm/10 a(p>0.1)。综合而言,自1960—2019年,全省峰丛洼地型地貌区四季降水均呈最快下降趋势,其中春季和秋季下降趋势显著,夏季和冬季下降趋势不显著。其余地貌类型区均有变化,但变化趋势不显著。
表4 1960-2019年不同地貌类型季节降水倾向率mm/10 a
从表4可知,春季降水在六大地貌区均呈减少趋势,除峰丛洼地减少趋势显著外,其余均不显著。从图6可以看出,六大地貌类型区春季降水累积距平值大致呈先上升后波动下降的变化趋势。其中,喀斯特低中山地貌区多年降水均值为294 mm,在1985年达到峰值,累积距平为575 mm,在2011年达到谷值,累积距平为-107 mm。峰丛谷地型地貌区多年降水均值为298 mm,在1978年达到峰值,累积距平519 mm,在2011年达到谷值,累积距平为-181 mm。峰林盆地(溶原)型地貌区,多年降水均值为303 mm,在1984年达到峰值,累积距平为572 mm,2011年达到谷值,累积距平为-170 mm。峰丛洼地型地貌区多年降水均值275 mm,在1963年达到谷值,累积距平为-83 mm,在1985年达到峰值,累积距平为609 mm。峰林谷地型地貌区多年降水均值为335 mm,在1984年达到峰值,累积距平为486 mm,在2011年达到谷值,累积距平为-144 mm。混合型地貌区多年降水均值为320 mm,在1985年达到峰值,累积距平为434 mm,在2011年达到谷值,累积距平为-128 mm。
图6 1960-2019年贵州省不同地貌类型春季降水距平变化
从表4可知,夏季降水在六大地貌区除峰丛洼地地貌区呈显著减少趋势外,其余地貌区均呈增加趋势。从图7可以看出,六大地貌类型区夏季降水累积距平值大致呈先波动变化后上升的变化趋势。其中,喀斯特低中山地貌区多年降水均值为560 mm,在1990年达到谷值,累积距平为-307 mm,在2002年达到峰值,累积距平为507 mm。峰丛谷地型地貌区多年降水均值为519 mm,在1990年达到谷值,累积距平为-358 mm,在2002年达到峰值,累积距平为488 mm。峰林盆地(溶原)型地貌区多年降水均值为530 mm,在1990年达到谷值,累积距平为-391 mm,在2002年达到峰值,累积距平为524 mm。峰丛洼地型地貌区多年降水均值为638 mm,在1975年达到谷值,累积距平为-272 mm,在2001年达到峰值,累积距平为617 mm。峰林谷地型地貌区多年累积降水均值为522 mm,在1990年达到谷值,累积距平为-549 mm,在2002年达到峰值,累积距平为269 mm。混合型地貌区多年降水均值为540 mm,在1990年达到谷值,累积距平为-385 mm,在2002年达到峰值,累积距平为290 mm。
图7 1960-2019年贵州省不同地貌类型夏季降水距平变化
从表4可知,秋季降水在六大地貌区除峰丛洼地地貌区呈显著减少趋势外,其余地貌区均呈不显著减少趋势。从图8可以看出,六大地貌类型区秋季降水累积距平值大致呈先上升后下降的变化趋势。其中,喀斯特低中山型地貌区多年降水均值为248 mm,在1983年达到峰值,累积距平为551 mm,在2013年达到谷值,累积距平为-130 mm。峰丛谷地型地貌区多年降水均值为235 mm,在1983年达到峰值,累积距平为488 mm,在2013年达到谷值,累积距平为-88 mm。峰林盆地(溶原)型地貌区多年降水均值为241 mm,在1983年达到峰值,累积距平为584 mm,在2009年达到谷值,累积距平为-76 mm。峰丛洼地型地貌区多年降水均值为258 mm,在1983年达到峰值,累积距平为670 mm,在2013年达到谷值,累积距平为-187 mm。峰林谷地型地貌区多年降水均值为245 mm,在1983年达到峰值,累积距平为501 mm,在2009年达到谷值,累积距平为-202 mm。混合型地貌区多年降水均值为234 mm,在1983年达到峰值,累积距平为432 mm,在2013年达到谷值,累积距平为-205 mm。
图8 1960-2019年贵州省不同地貌类型秋季降水距平变化
从表4可知,冬季降水在六大地貌区有增有减,且变化趋势均不显著。从图9可以看出,六大地貌类型区冬季降水累积距平值大致呈振荡变化的趋势。其中,喀斯特低中山型地貌区多年降水均值为68 mm,在1968年达到谷值,累积距平为-75 m,在1970年达到峰值,累积距平为113 mm。峰丛谷地型地貌区多年降水均值为70 mm,在1968年达到谷值,累积距平为-82 mm,在1970年达到峰值,累积距平为112 mm。峰林盆地(溶原)型地貌区多年降水均值为75 mm,在1970年达到峰值,累积距平为105 mm,在2014年达到谷值,累积距平为-89 mm。峰丛洼地型地貌区多年降水均值为70 mm,在1968年达到谷值,累积距平为-67 mm,在1994年达到峰值,累积距平为99 mm。峰林谷地型地貌区多年降水均值为93 mm,在1972年达到峰值,累积距平为79 mm,在2014年达到谷值,累积距平为101 mm。混合型地貌区多年降水均值为88 mm,在1988年达到谷值,累积距平为-110 mm,在1997年达到峰值,累积距平为58 mm。
3.7 季节降水周期性
从1960—2019年贵州省四季降水周期图可看出近60 a四季均有多个正负相位交错的振荡中心,说明季节降水量存在着明显的周期振荡。其中,春季降水量存在5个主振荡周期,分别是28 a,18 a,13 a,9 a和5 a。其中,28 a左右的振荡周期最强(图 10),是春季降水的第一主周期,整个研究时段内振荡最为明显且一直存在,经历了“偏多—偏少—偏多—偏少—偏多—偏少—偏多”7个交替变化; 其次,13 a左右的振荡周期仅在20世纪80—90年代较为明显,为春季降水的第三周期; 9 a左右的降水周期主要出现在21世纪,其余周期变化均不明显。
夏季降水主要存在28 a,15 a,7 a,4 a共4个周期(图 11)。其中,观察小波方差发现4个较明显的峰值,分别是28 a对应的第1峰值,15 a对应的第2峰值; 其中28 a为最强振荡主周期,且经历了“偏多—偏少—偏多—偏少—偏多—偏少—偏多”7个交替变化。
图9 1960-2019年贵州省不同地貌类型冬季降水距平变化
图 10 1960-2019年贵州省春季降水量变化
1960—2019年贵州省秋季降水主要存在28 a,18 a,9 a,5 a和3 a共5个周期(图 12)。其中28 a的振荡周期最强,是秋季降水的第一主周期,整个研究时段内振荡最为明显,且经历了“偏多—偏少—偏多—偏少—偏多—偏少—偏多”7个交替变化。
1960—2019年贵州省冬季降水主要存在28 a,18 a,8 a,3 a共4个周期(图 13)。其中,以18 a的周期振荡最强,为第一主周期,经历了“偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多—偏少—偏多”11个交替变化。
图 11 1960-2019年贵州省夏季降水量变化
图 12 1960-2019年贵州省秋季降水量变化
图 13 1960-2019年贵州省冬季降水量变化
3.8 季节降水突变性
对1960—2019年贵州省近60 a四季降水时间序列进行Mann-Kendall非参数突变检验[33](图 14); 其中春季降水量的检验统计值UFk与UBk曲线在临界值±1.96之间于1984年附近有一个交点(图 14A); 结合图4A可以看出,1984年附近累积距平达到峰值,且在1984年两侧累积距平曲线变化趋势由上升变为下降,由此可得,春季降水量的时间突变点是1984年。夏季降水量在临界值±1.96两侧交点分别在1963年、1966年、2004年、2006年、2007年、2014年、2017年和2018年附近; 结合图4B知在2014年两侧,降水变化趋势由下降变为上升,因此夏季降水量变化趋势时间突变点为2014年。秋季降水量的检验统计量UFk与UBk曲线交于1984年附近; 结合图2C知,1984年之前降水以上升为主,1984年之后降水量以下降为主,所以秋季降水量的突变时间为1984年。冬季降水量在临界值±1.96之间交点分别在1980年、1982年、1983年、1984年、1987年、1997年、1998年、2005年、2014年、2015年和2018年附近; 结合图4D知,冬季累积降水距平值波动较大没有显著的增减趋势变化点,因此冬季降水量变化趋势没有发生突变。
图 14 1960-2019年贵州省季节降水量突变检验
