1. 引言
1960年以来，化石燃料燃烧加剧、碳排放量逐年增长带来近百年甚至近千年来前所未有的相关气候变化。出于对潜在的全球气候变化问题的重视，由世界气象组织(WMO)和联合国环境规划署(UNEP)于1988年成立了联合国政府间气候变化专门委员会(Intergovernmental Panel on Climate Change, IPCC)，并在1990年发布了《第一次评估报告》 [1]。自此气候变化的成因，其潜在的环境、社会、经济影响，以及缓解气候变化的可能性对策就成为了各国政府和各类国际研究机构的热点议题。
根据IPCC第五次评估报告结果，全球气候变暖趋势不可避免，其具体表现在全球气温上升、海平面上升、极地冰盖消融、极端气候事件频增等方面 [2]。气候变化已经对人类健康、人类安全、生计与贫困等人类福祉方面产生了负面影响，增加了额外负担。气候变化将加剧现有问题，放大各类冲突，为国家安全带来新的隐患 [3]。在全球变暖气候的威胁下，应对极端天气的可能性对策成为全球共同关注的焦点。其中各类极端降水将导致城市内涝危害城市交通安全 [4]、农田积水危害粮食安全 [5]、地表水量短时间激增影响堤坝泄洪 [6]，与人民的生产生活安全息息相关，是国内外气象工作者们共同关注的重点问题。
四川东部即四川盆地及周围山地属中亚热带湿润气候区，又兼有海洋性气候特征，全年温暖湿润；盆地地区云量多，晴天少 [7]。成都是四川的省会城市，位于四川盆地西部，地势平坦、物产丰富、农业发达，自古有“天府之国”的美誉，也是国家重点高新技术产业基地、商贸物流中心和综合交通枢纽 [8]。温江地处成都，温江站的气象观测数据的相关研究分析结论可以较好代表成都气象要素变化。对该站点的近60年降水数据的研究分析，能够了解成都地区的极端降水变化特征，确保人民生产生活正常运行和社会经济可持续发展。
据此，本文通过百分位法确定温江站的极端降水事件阈值，通过一元线性方程进行极端降水事件变化趋势的表示，同时由相关系数检验判断变化趋势是否显著，通过累加距平进行变化趋势的突变检验。以温江站作为代表探究成都地区极端降水的变化特征。
2. 资料与方法
2.1. 资料概况
本文重点对四川省成都市近60 a的极端降水变化特征进行分析，所用的数据为成都市温江站1960~2019年的逐日降水资料。气象数据来源于中国气象数据网。
2.2. 研究方法
本文采用百分位法 [9] [10] 进行极端降水阈值的确定，根据确定后的极端降水阈值统计年际与年代的极端降水量、极端降水日数、极端降水强度；用线性倾向估计法 [11] 对各个气象要素的线性变化趋势及其显著性进行分析；用累积距平方法 [11] 对年际变化进行突变分析。
3. 研究结果
3.1. 极端降水量的变化特征
3.1.1. 极端降水量的年际变化特征
通过百分位法的计算，成都市温江站极端降水事件的阈值为27.1 mm。根据传统的气象学名词解释，24 h内降雨量大于等于25 mm且小于50 mm的天气现象叫做大雨。由此来看，对于夏季少暴雨少急雨、短时雨量较小的成都地区来说，24 h内降雨量达到大雨级别即可粗略认定为发生了极端降水天气过程。
从由线性倾向估计方法得出的线性直线(图1)来看，1960~2019年的极端降水量年际变化呈现随年份增加的负相关趋势，其线性相关系数r值为−0.181，未通过α = 0.05的显著性检验，表明极端降水量的线性相关性较弱，即成都地区的极端降水量减少的年际变化不显著。近60 a年平均极端降水量为425.6 mm，最大值出现在1961年为886 mm，最小值出现在1965年为105.9 mm。
Figure 1. 1960~2019 interannual variation trend of annual extreme precipitation in Wenjiang
图1. 1960~2019年温江站年极端降水量变化趋势
根据累积距平所绘制的曲线(图2)来看，1982年前累积距平曲线呈波动上升趋势，并在1982年达到最大值1208.7 mm，成为最大值点；从1982~2012年累积距平曲线呈波动下降趋势，距平值迅速减小，在2012年达到最小值−502.4 mm；在2012~2016年期间，累积距平值波动平缓，并在2016年再次取得最小值。即可以初步判断1982年和2012年为极端降水量变化趋势的突变时间点。1982年前年极端降水量平均值为478.5 mm，1982~2012年极端降水量平均值为371.9 mm，2012年后年极端降水量平均值为497.4 mm。
Figure 2. 1960~2019 Accumulated anomaly of annual extreme precipitation in Wenjiang
图2. 1960~2019年温江站年极端降水量累积距平
3.1.2. 极端降水量的年代变化特征
根据极端降水量的年代变化特征表(表1)，可以看出80 s和00 s分别存在一个极端降水量变化的低值区，印证了前文极端降水量年际变化的累积距平分析结果的基本正确。年代平均极端降水量在60 s达到最大为531.9 mm，其后逐渐减小，近10 a又开始增多。
Table 1. 1960~2019 decadal variation trend of extreme precipitation in Wenjiang
表1. 1960~2019年温江站极端降水量年代变化趋势
3.2. 极端降水日数的变化特征
3.2.1. 极端降水日数的年际变化特征
本文根据国际组织推荐标准的极端降水统计方法得出极端降水事件阈值，即第95个百分位上的值定义为极端降水事件发生。由此可得，一般来说一年内的极端降水事件日数应当在19日以内。成都地区1960~2019年的年均极端降水日数为7.9 d，60 a内无年极端降水日数在19日以上。说明成都地区的极端降水日数较少，受到由于极端降水导致的危险伤害的日数较少。
从由线性倾向估计方法得出的线性直线(图3)来看，1960~2019年的极端降水日数的年际变化同极端降水量的年际变化趋势相似，均呈现随年份增加的负相关趋势，在60 a的时间序列上呈波动下降。但是其得出的线性相关系数r值为−0.179，未通过α = 0.05的显著性检验，即极端降水日数下降的趋势不显著。近60 a年平均极端降水日数为7.9 d，最大值出现在1965年和1983年为14 d，最小值出现在1964年为2 d。
根据累积距所绘制的曲线(图4)来看，1990年前累积距平曲线呈波动上升趋势，并在1990年达到最大值23.7 d；从1990~2012年累积距平曲线呈波动下降趋势，距平值迅速减小，在2012年达到最小值−5.1 d；在2012~2019年期间，波动较为平缓，即可以初步判断1990年和2012年为极端降水日数变化趋势的突变时间点。1990年之前年极端降水日数平均值为8.6 d，1991~2012年极端降水日数平均值为3.6 d，2012年之后年极端降水日数平均值为8.6 d。
3.2.2. 极端降水日数的年代变化特征
根据极端降水日数的年代变化特征表(表2)，可以看出60 s存在一个极端降水日数的最高值为9.4 d，00 s存在一个极端降水日数的最低值为6.3 d。
Figure 3. 1960~2019 interannual variation trend of days of extreme precipitation in Wenjiang
图3. 1960~2019年温江站年极端降水日数变化趋势
Figure 4. 1960~2019 accumulated anomaly of days of extreme precipitation in Wenjiang
图4. 1960~2019年温江站年极端降水日数累积距平
Table 2. 1960~2019 decadal variation trend of days of extreme precipitation in Wenjiang
表2. 1960~2019年温江站极端降水日数年代变化趋势
3.3. 极端降水强度的变化特征
3.3.1. 极端降水强度的年际变化特征
极端降水强度由年极端降水量和年极端降水日数之比得出，表明了极端降水事件对于当地生活生产影响的可能程度。和极端降水量以及极端降水日数相同，是极端降水研究中重要的极端降水指数。
从由线性倾向估计方法得出的线性直线(图5)来看，直线与x轴基本平行，表示了1960~2019年的极端降水强度的无显著的增加或减少，即成都地区极端降水强度在60 a的时间序列上呈上下波动变化。近60 a年平均极端降水强度为54.0 mm/d，最高值出现在1997年为87.2 mm/d，最低值出现在2011年为36.1 mm/d。
Figure 5. 1960~2019 interannual variation trend of annual extreme precipitation intensity in Wenjiang
图5. 1960~2019年温江站年极端降水强度变化趋势
根据累积距平所绘制的曲线来看(图6)，在1973年前累积距平曲线达到最大值，在1997年累积距平曲线达到最小值，曲线整体呈现波动变化趋势，可初步判断1973年和1997年分别出现了极端降水强度的突变。1973年前年极端降水强度平均值为57.0 mm/d，1974~1997年极端降水日数平均值为49.9 mm/d，1997年之后年极端降水强度平均值为56.2 mm/d。
Figure 6. 1960~2019 accumulated anomaly of annual extreme precipitation intensity in Wenjiang
图6. 1960~2019年温江站年极端降水强度累积距平
3.3.2. 极端降水强度的年代变化特征
根据极端降水强度的年代变化特征表(表3)，可以看出80s存在极端降水强度最低值为49.4 mm/d；90 s存在极端降水强度最大值为56.9 mm/d，主要原因为1997年极端降水强度为87.2 mm/d，强度最大。年代降水强度多年平均值为54.2 mm/d，同各个年代的极端降水强度基本相近，变化了较小，和前文年际变化趋势分析相同。
Table 3. 1960~2019 decadal variation trend of extreme precipitation intensity in Wenjiang
表3. 1960~2019年温江站极端降水强度年代变化趋势
4. 结论
本文利用1960~2019年温江站的逐日降水量数据，在经过对于缺测值的初步处理后，通过百分位法计算得出了极端降水事件阈值。从极端降水量、极端降水日数、极端降水强度三种不同极端降水指数方面，在年际与年代两种时间尺度上进行了分析，结论如下：
1) 极端降水量在年际时间序列上呈现波动下降趋势，但下降趋势不显著。1982年与2012年为极端降水量变化趋势的突变时间点。60 s出现了极端降水量的最高值，最低值出现在00 s。
2) 极端降水日数在年际时间序列上同样呈波动下降趋势，且下降趋势不显著。极端降水日数的突变时间点同极端降水量基本保持一致，且最高值和最低值同样出现在60 s和00 s。
3) 极端降水强度呈现出波动变化，但无线性趋势。年极端降水强度值在多年平均值线上下小幅波动，年代极端降水强度变化不大。