1. 引言
大气中的液态水或固态水在重力的作用下,克服空气阻力,从高空降落到地面的现象称为降水 [1]。降水是水循环过程中最基本的环节,是地表径流的本源,是地下水的主要补给来源,也是清洁水的最终来源,而降水在时空分布上的不均匀是引起洪涝和旱灾的直接原因 [1],所以降水是连接天气、气候和水循环等很多方面的关键物理过程,是关系到国计民生的最重要的气象要素之一。IPCC第五次评估报告指出,近百年来全球气候变暖毋庸置疑,气温升高使得大气中可承载的最大水汽量含量增多,造成水循环明显加快,降水的时空分布特征也随之发生了很大变化,高纬度和一些中纬度地区降水量呈显著增加的趋势,以至于全球范围内频繁发生极端降水事件 [2] [3] [4]。
本文的研究区域为我国四川地区。四川盆地是年降雨量比较充沛的地区,四川盆地内的降水对其下游长江的水位及三峡工程的建设等有很大的影响,由于每年的降雨量在空间和时间上分配并不均匀,盆地内常会有降雨过强导致的暴雨天气,给盆地带来洪涝、泥石流、山地滑坡等严重灾害,给农业、交通、建筑等带来重大影响 [5]。
就四川地区强降水分布特征的研究而言,四川地区的极端降水值、极端降水强度以及极端降水频数都在增多 [6]。张顺谦表明:近50年来四川区域性暴雨日数、暴雨次数和极端降水日数呈略微减少的趋势,极端降水强度、极端降水量占年降水量的比例呈微弱增加的趋势,但极为不显著,与我国极端降水趋多趋强的变化特征不一致 [7]。胡豪然等利用四川盆地1961~2006年降水资料表明在汛期该地区极端降水事件发生频次的分布与降水量分布的差异比较大,自西向东有阶梯状递减趋势,汛期内极端降水事件发生的频次除了在重庆东部区和四川盆地西部区分别呈较弱的增长和减少趋势以外,在其余各区的线性趋势都较为明显 [8]。
进一步了解全球变暖背景下四川地区强降水的时空变化特征,不仅对认识该地区气候背景状况和该地水资源分布有重要作用,而且对于提高四川短期气候预测能力具有借鉴意义。
2. 数据和方法
2.1. 数据
本文使用的是四川省1980~2018年的欧洲中期数值预报中心(ECMWF) ERA5降水再分析数据,气象要素包括平均总降水率、总降水量、平均降雪率、降雪量等。在本文中,只使用总降水量这一个要素来分析四川地区强降水的时空变化特征。时间分辨率:小时(h);时间跨度:1980~2018。
2.2. 方法
2.2.1. 降水强度的定义
四川有些地区的年最大降雨量不能够达到气象部门规定的暴雨标准,如果以50 mm日降水量来定义极端降水事件的绝对阈值,那么有些地区将是空白区并且无法分析该地区极端降水事件。基于这种影响,赵衍斌 [9] 采用非参数法(百分位法)在研究中给出了四川地区强降水的定义:根据每个测站的日降水量,定义不同地区强降水事件的阈值,并按从小到大的顺序对1980~2019年的日降水量资料(降水量大于0)进行排序,并将某站点降水量的95%的数值定义为强降水事件的阈值,并且当该站的日降水量超过该站的强降雨阈值时,就称该日发生了极端降水事件。
根据强降水事件发生的次数,对四川地区的强降水强度进行定义:
(1)
其中S代表强降水强度,t代表某一年强降水事件发生的次数,sr代表该年强降水逐日降水量 [9]。
2.2.2. Mann-Kendall检验
Mann和Kendall首先提出了Mann-Kendall突变检验的原理并发展了这一方法,它是一种非参数检验方法,也称为无分布检验,一开始仅仅用来检测已知序列的变化趋势,后来经过在实践中不断的完善和改进,形成了成了一种检验序列突变的方法 [10]。Mann-Kendall突变检验方法的优势是样本无需遵从确定的分布,就能够大致确定变化趋势的起始位置。
Mann-Kendall突变检验的计算方法参照《现代气候统计诊断与预测技术》 [11]:
若原始时间序列为
,且
表示的是第i个样本
大于
的累计数
,则有:
(2)
假设原始的时间序列随机独立,那么
的均值和方差为:
(3)
(4)
将(2)、(3)、(4)式标准化:
(5)
经过上述计算得到正向样本序列的统计量UF,再利用类似计算,得到反向样本序列的统计量UF,得到正、反向样本序列的统计量UF和UB之后,在所给定的显著性水平下(本文给定显著性水平α = 0.05),分别画出UF和UB曲线,然后利用两者曲线的变化来分析研究中物理量样本的变化趋势。若所得到的曲线图上UF > 0那么表明所检测序列呈现上升趋势,UF < 0则表明所检测序列呈现下降趋势,当U = 0时,则所检测序列变化持平,且如果UF值超过临界线,那么可以判断为所检测序列变化显著上升或者显著下降,如果UF和UB的交点出现在临界值之间,那么UF与UB交点所在位置对应的时刻就是突变发生的时刻。
2.2.3. 其他方法
本文用到的其它方法有:1) 通过设立统计阈值,挑选出具有统计意义的极端降水事件,利用区域平均方法获得降水量的空间分布、降水频次的空间分布以及降水强度的空间分布;2) 采用经验正交函数(EOF)分析方法来获得主要空间变化模态,来提供降水变化的主要空间特征;3) 采用线性回归来计算降水次数的长期变化趋势;4) 采用功率谱方法获得降水次数的主要变化周期及演变规律。
3. 四川地区强降水时空分布特征
3.1. 四川地区强降水空间分布特征
3.1.1. 阈值
图1为1980~2018年来四川地区极端降水量的阈值分布图。分析图1可以发现,四川地区强降水的阈值分布自西向东逐渐增大,自南向北逐渐增大,整体表现为西小东大、南小北大的特征。由于四川地区西高东低的特殊地形影响,川西高原的阿坝、甘孜等地阈值在20 mm以下,普遍偏低;四川盆地以及盆地以东地区阈值在20 mm以上,普遍较高;广元、巴中、南充、广安、达州一带,成都、雅安一带以及西昌、攀枝花一带是四川盆地及以东地区主要存在的3个阈值的大值中心,强降水阈值均可达到30 mm以上,可以看出四川不同地区的强降水阈值差异明显,阈值最大差值可以达到20 mm。
![](//html.hanspub.org/file/9-2420349x20_hanspub.png)
Figure 1. Threshold distribution map of extreme precipitation in Sichuan
图1. 四川地区极端降水量的阈值分布图
3.1.2. 强降水量的空间分布
图2为1980~2018年来四川地区季节平均强降水量空间分布图。可以看出四川地区强降水量自西向东月平均降雨量逐渐增多,川西高原的阿坝、甘孜等地四季降水都最少的地方,四川西部和西南地区四季降水比较多。春季,整个四川地区降水大致呈西南–东北带状分布,自西向东呈“少–多–少”分布规律,大部分地区强降水量不超过120 mm,也出现了凉山、雅安、广安三个强降雨量大值中心,凉山、雅安等地局部降水量可达120 mm以上,广安局部降水量高达160 mm以上。夏季,整个四川地区降水依旧大致呈西南–东北带状分布,与春季有所不同的是,夏季强降水自西向东呈“少–多”分布规律,而且夏季强降雨量的大值中心相对于春季而言逐渐东移,移到四川中东部德阳、广元等地,四川西部地区强降水量基本不超过200 mm,德阳局部降水量可达500 mm以上,广安局部降水量高达600 mm以上。秋季,整个四川地区降水分布比较散,自西向东呈“多–少–多”分布规律,降雨量整体减少,大部分地区强降水量不超过80 mm,最大值中心移到四川东北部广元东北边缘地区,出现了凉山、广元两个强降雨量大值中心,凉山局部降水量可达120 mm以上,广安局部降水量高达240 mm以上。冬季,没有强降水的发生。
![](//html.hanspub.org/file/9-2420349x21_hanspub.png)
Figure 2. Spatial distribution map of seasonal average heavy precipitation in Sichuan from 1980 to 2018
图2. 1980~2018年来四川地区季节平均强降水量空间分布图
图3为1980~2018年来四川地区年平均强降水量空间分布图。可以看出1980~2018年四川地区多年年平均强降水量总体上表现为西少东多的趋势,即自西向东年平均强降雨量逐渐增多。而且存在明显的区域变化特征,在各地区分布差异很大,多年年平均强降水量最小值162.533 mm,最大值为2323.33 mm。阿坝和甘孜等地多年年平均强降水量较少,尤其是甘孜西北部地区,多年年平均强降水量值接近最小值162.533 mm,除此之外凉山西北部和东北部地区多年年平均强降水量也较少,在200~300 mm之间。广元、巴中和达州北部等地多年平均强降水量较多,多年年平均强降水量值达到700~800 mm之间,其中广元和巴中交界处的北部达到四川境内多年平均强降水量的最大值,另外,在雅安和甘孜的交界处以及凉山南部局部也有多年年平均强降水量大值中心,多年年平均强降水量在600 mm以上。
将图3与图1进行对比,可以发现四川地区极端降水量阈值的与四川地区多年平均强降水量的分布趋势基本上是保持一致的,阈值低的地区,极端降水事件的多年平均强降水量也较低,反之也成立。
![](//html.hanspub.org/file/9-2420349x22_hanspub.png)
Figure 3. Spatial distribution map of annual average heavy precipitation in Sichuan from 1980 to 2018
图3. 1980~2018年来四川地区年平均强降水量空间分布图
3.1.3. 强降水发生频率的空间分布
图4为1980~2018年间四川地区各个季节累计强降水次数。可以看出四川地区强降水次数的分布存在明显的区域变化特征,各个季节都有不同的趋势。春季,整个四川地区强降水次数大致呈西南-东北带状分布,自西向东呈“少–多–少”分布规律,即自西到东强降水次数先增大后减少,大部分地区强降水次数不超过10次,出现了阿坝、甘孜、凉山一带为主轴的强降水次数大值中心,甘孜、凉山等地局部强降水次数可达10次以上,广安局部强降水次数高达15次以上。夏季,整个四川地区降水依旧大致呈西南–东北带状分布,与春季有所不同的是,夏季强降水分布规律与春季相反,自西向东呈“多–少–多”分布规律,即自西到东强降水次数先减少后增多,也可以说是夏季强降水次数的大值中心相对于春季而言逐渐东移,移到四川东部,这样夏季四川地区强降水次数有两个大值中心,一个大值中心位于甘孜西部地区,强降水次数可达44次以上,另一个大值中心位于巴中、广元、绵阳、德阳、南充、遂宁、成都、资阳、眉山、雅安、宜宾、内江、乐山以及泸州等地,强降水次数达到44次以上。秋季,整个四川地区强降水分布比较散,自西向东呈“多–少–多”分布规律,强降水次数整体减少,大部分地区强降水次数18次,强降水次数最大值中心移到四川东部和阿坝北部,出现了阿坝、广元两个强降雨次数大值中心,阿坝局部强降水次数可达22次以上,广安局部强降水次数可达18次以上。冬季,没有强降水的发生。
图5为1980~2018年间四川地区累计强降水次数。可以看出1980~2018年四川地区39年累计强降水次数在各地区分布差异很大,整个四川地区强降水次数大致呈西南–东北带状分布,自西向东呈“少–多–少”分布规律,即自西到东强降水次数先增大后减少,多年累计强降水量次数最多为713次,最少为488次。大部分地区强降水次数都在60次以上,出现了阿坝、甘孜北部为中心、甘孜中部为中心、攀枝花为中心以及绵阳、德阳为中心的四个强降水次数大值中心即强降水高频区,绵阳、德阳强降水高频区、甘孜中部降水高频区强降水次数可达61次以上,阿坝、甘孜北部强降水高频区、攀枝花强降水高频区强降水次数可达62次以上;凉山东北部、乐山、阿坝南部局部地区形成了强降水低频区,强降水次数在57次以下;其余地区强降水次数在58~60次之间。
![](//html.hanspub.org/file/9-2420349x23_hanspub.png)
Figure 4. The cumulative number of heavy precipitation in each season in Sichuan from 1980 to 2018
图4. 1980~2018年间四川地区各个季节累计强降水次数
![](//html.hanspub.org/file/9-2420349x24_hanspub.png)
Figure 5. The cumulative number of heavy precipitation in Sichuan from 1980 to 2018
图5. 1980~2018年间四川地区累计强降水次数
将图5和图1进行对比,可以发现极端降水发生次数同极端降水阈值的分布情况并不一致,阈值大的地区,其强降水发生次数未必就少,反之亦然,二者之间并无必然的联系。进一步对比发现,极端降水发生次数较多的以阿坝、甘孜北部为中心、甘孜中部地区为中心、攀枝花为中心以及绵阳、德阳为中心的四个强降水高频区的极端降水阈值都较小。
3.1.4. 强降水降水强度的空间分布
图6为1980~2018年四川地区强降水强度的空间分布图。可以看出海拔较高的地区比如川西高原等地,强降水强度较弱,而海拔较低的川东平原,强降水强度较强,即四川地区强降水强度空间分布主要表现为自西向东逐渐增大的趋势,这种趋势与四川地区西高东低的地形存在一致性。四川地区强降水强度最大值和最小值差值较大,四川地区强降水强度存在三个强度高值中心,雅安南部与阿坝的交界处有一高值中心,强度高达175 mm/d,达州有一高值中心,强度可达150 mm/d,凉山有一高值中心,强度高达150 mm/d;四川地区强降水强度存在三个强度低值中心,内江和甘孜南部分别有一个低值中心,强度分别为低于75 mm/d、50 mm/d,甘孜北部和阿坝北部则有一个强度低值带,强降水强度基本在50 mm/d以下,四川地区强降水强度最小值为28.089 mm/d。
将图6与图1进行对比,可以发现四川地区极端降水量阈值的分布与强降水强度的分布趋势基本一致,阈值低的地区,强降水强度也较低,反之亦成立。将图6与图3进行对比,可以发现四川地区年平均强降水量的分布与强降水强度的分布趋势基本一致,年平均强降水量少的地区,强降水强度也较低,反之亦成立。对于强降水事件,阈值越低、强降水量愈少,强降水事件强度越小,反之亦成立。
![](//html.hanspub.org/file/9-2420349x25_hanspub.png)
Figure 6. Spatial distribution map of heavy precipitation intensity in Sichuan from 1980 to 2018
图6. 1980~2018年四川地区强降水强度的空间分布图
3.1.5. 主要空间变化模态分析
图7给出了四川地区多年逐月平均强降水量的EOF前四个模态的空间分布特征,来获得主要空间变化模态,且分离出显著的主要模态,以此来提供降水变化的主要空间特征。第一模态的特征向量占方差总贡献的58.147%,是四川地区多年逐月平均强降水量的最主要分布型。四川西部大部分地区均为负值,负值中心位于甘孜、阿坝和凉山北部地区,中心值比较集中,正值中心在四川西北部地区,主要分布在广元、南充、达州等地,中心值为0.03~0.04之间。四川地区多年逐月平均强降水量空间分布的第一类型具有多年逐月平均强降水量“东多西少”的分布特征。第二模态占方差总贡献的23.0249%,也是一个较重要的分布型。较大的正值区出现在四川西部大部分地区,正值中心广元、南充、巴中等地的中心值大于0.06,负值中心在四川南部的凉山地区。四川地区多年逐月平均强降水量空间分布的第二类型主要表现出四川西北部强降水量偏多,其余地区强降水量一致性偏小的分布特征。第三模态占方差总贡献的11.6517%。正值区出现在甘孜、阿坝大部分地区以及凉山、攀枝花等地,正值中心中心值大于0.03,负值中心在巴中、南充、达州等地。四川地区多年逐月平均强降水量空间分布的第三类型从四川西北到东南主要表现出“− → + → −”的格局。第四模态占方差总贡献的6.42744%。正值区出现在绵阳、德阳、成都等地的部分区域,正值中心中心值大于0.05,负值中心在四川西北部广元、南充、绵阳等地。四川地区多年逐月平均强降水量空间分布的第四类型从四川西部到东部主要表现出“+ → −”的格局。
![](//html.hanspub.org/file/9-2420349x26_hanspub.png)
Figure 7. Spatial distribution of EOF based on monthly average heavy precipitation in Sichuan for many years
图7. 基于四川地区多年逐月平均强降水量的EOF空间分布
3.2. 四川地区强降水时间分布特征
3.2.1. 极端降水事件发生次数的年际变化
图8给出了逐年夏季极端降水事件频数时间序列及其拟合的折线图来分析强降水的长期变化趋势。从图中可以看出,1998年夏季四川地区极端降水事件发生的次数最多,2006年夏季四川地区极端降水事件发生的次数最少,且近39年来总体呈明显的下降趋势,这种倾斜倾向通过了信度为90%的显著性检验。从年际变化来看,20世纪90年代以前四川地区极夏季端降水事件发生的次数变化趋势基本平稳,从1993年以后开始呈波动变化。
3.2.2. 极端降水事件发生次数的周期特征
图9给出了四川地区逐年夏季极端降水事件频数功率谱分布来分析四川地区夏季强降水的主要变化周期及演变规律。可以看出四川地区夏季极端降水的功率谱年际变化存在明显的2.2年、3.5年和13年峰值,但只有2.2年的峰值通过了95%的显著性检验,而其他峰值没有通过显著性检验,说明其很可能是由随机性因素所引起的,而并非真正的波动周期。综上所述,四川地区39年夏季极端降水存在约2.2年的准周期。
![](//html.hanspub.org/file/9-2420349x27_hanspub.png)
Figure 8. The time series and fitting of the frequency of summer extreme rainfall events in Sichuan
图8. 四川地区逐年夏季极端降水事件频数时间序列及其拟合
![](//html.hanspub.org/file/9-2420349x28_hanspub.png)
Figure 9. Frequency power spectrum of extreme summer precipitation events in Sichuan
图9. 四川地区逐年夏季极端降水事件频数功率谱
3.2.3. 极端降水事件发生次数的突变年份
利用1980年到2018年间四川地区逐年夏季极端降水事件频数时间序列进行M-K显著性检验对其气候态的突变点进行分析。从图10中曲线的变化可以看出,四川地区逐年夏季极端降水事件频数的变化整体表现为20世纪90年代呈现先上升后下降的波动变化,但这种短暂的振荡现象并无显著趋势,在21世纪初,四川地区逐年夏季极端降水事件频数开始呈现出稳定下降的趋势,并在2003年达到了0.05的显著性水平,在其后有一个较为显著的突变区域,UF和UB的交点表示四川地区逐年夏季极端降水事件频数突变的年份在1981、1985年与1991年,由于交点同时位于虚线之间,表明两次突变都通过了0.05显著性水平检验。综上可知,近39年来四川地区逐年夏季极端降水事件频数突变发生在1981、1985年和1991年。
![](//html.hanspub.org/file/9-2420349x29_hanspub.png)
Figure 10. MK statistics of the frequency of extreme summer precipitation events in Sichuan
图10. 四川地区逐年夏季极端降水事件频数MK统计量
4. 结论
通过对1980~2018近39年来四川地区强降水时空特征分析,主要得出以下几点结论:
1) 四川地区强降水的阈值整体表现为西小东大、南小北大的趋势,存在明显的区域变化特征。
2) 四川不同地区的多年年平均强降水量总体上表现为西少东多的趋势;四川地区39年累计强降水次数自东向西呈“少–多–少”分布规律;四川地区强降水强度空间分布自西向东逐渐增大,这种趋势与四川地区西高东低的地形存在一致性。
3) 四川地区强降水量和强降水次数在各个季节都有不同的趋势,降水主要集中于雨季,强降水量大值中心和强降水高频区都会随着季节变化逐渐东移。
4) 四川地区极端降水量阈值的分布与四川地区多年平均强降水量的分布、强降水强度的分布趋势基本一致,阈值低的地区,极端降水事件的多年平均强降水量也较低,强降水强度也较低,反之也成立。四川地区极端降水量阈值的分布与极端降水发生次数的分布情况并不一致,阈值大的地区,其强降水发生次数未必就少,反之亦然,二者之间并无必然的联系。
5) 四川地区多年逐月平均强降水量空间表现为:第一模态具有多年逐月平均强降水量“东多西少”的分布特征,第二模态主要表现出四川西北部强降水量偏多,其余地区强降水量一致性偏小的分布特征。第三模态自北向南呈现“− → + → −”,第四模态从四川西部到东部主要表现出“+ → −”的格局。
6) 四川地区逐年夏季极端降水事件频数的变化整体表现为20世纪90年代呈现波动变化,在年际变化上总体呈明显的下降趋势;四川地区1980~2018年夏季极端降水次数存在约2.2年的准周期;四川地区逐年夏季极端降水事件频数突变发生在1981、1985年和1991年。