四川省強降水分布及變化特征研究

次旺頓珠 白瑪多吉 羅桑丹增

摘要 利用1980—2012年四川省36個站點逐日降水資料，運用數理統計、一元線性回歸、M-K突變檢驗、EEMD經驗模態分解等方法，研究了四川省強降水分布以及時間演變規律。結果表明：四川省總體年平均降水量為1 039.6 mm，且年降水量總體呈遞減趨勢，每10年減少22.9 mm左右;中部地區年平均降水量相對較多，東部地區年平均降水量其次，西部地區年平均降水量相對較少。四川省降水量季節變化春季呈上升趨勢，夏季、秋季和冬季呈現下降趨勢。四川省大雨日數呈下降趨勢，大約每10年大雨日數隨時間會減少4.0 d。暴雨日數也呈下降趨勢，大約每10年暴雨日數隨時間會減少1.5 d。中部地區大雨日數、暴雨日數均較多，東部地區其次，西部地區相對較少。1992年為四川省年降水量開始突變年;降水序列在年際演變過程中，存在2～3年小尺度、5～10年中尺度、20年左右大尺度這3類時間尺度變化規律。根據趨勢項得出四川省年降水量呈遞減趨勢。

關鍵詞 降水分布;變化特征;一元線性回歸;氣候傾向率;M-K突變檢驗;EEMD經驗模態分解;四川省

中圖分類號 S161.6文獻標識碼 A

文章編號 0517-6611（2019）04-0224-07

Autonomous Region，Bangor，Tibet 852500;3.Tibet Autonomous Region Atmospheric Detection Technology and Equipment Support Center，Lhasa，Tivet 850001）

Abstract Using the daily precipitation data of 36 stations in Sichuan Province from 1980 to 2012，the mathematical statistics，linear regression，MK mutation test，EEMD empirical mode decomposition were used to study the distribution and time evolution of strong precipitation in Sichuan Province.The results showed that the overall average annual precipitation was 1 039.6 mm in Sichuan Province，and the mean annual precipitation showed a decreasing trend，every 10 years reduced about 22.9 mm.The annual precipitation in central area of relatively large，the eastern area rainfall secondly，western region the average annual precipitation was less.Based on the variation of seasonal precipitation of Sichuan Province，the spring was on the rise，summer，autumn and winter showed a downward trend.Heavy rain days showed a decreasing trend，approximately every 10 years the number of heavy rain days time would be reduced by 4 days.Rainstorm day number also declined，approximately every 10 years the number of heavy rain days time would be reduced by 1.5 days.The central area of heavy rain days，rainstorm days were more，the eastern region secondly and western region was less.The results of MK test show that，Sichuan Province in 1992 for the annual precipitation to the mutation years.Based on the annual precipitation was EEMD empirical mode decomposition obtains，precipitation sequence in the interannual evolution，variation of these 3 time scales are small scale，2-3 years 5-10 years，20 years or so large scale.

Key words Precipitation distribution;Variation characteristics;Linear regression;Climate tendency rate;M-K mutation test;EEMD empirical mode decomposition;Sichuan Province

近年來，氣候變化形勢愈發嚴峻，從20世紀70年代初，全球氣候就開始發生顯著變化，降水特征在一定程度上也發生了變化，氣候的嚴峻變化形勢已引起各個國家和人民的極大關注[1-5]。

目前，越來越多的學者開始研究氣候變化問題。周長艷等[6]對四川近50年降水氣候特征進行分析，結果表明年平均暴雨日數從西到東呈“增—減—增”的變化趨勢。孫銀鳳等[7]指出，南京市年降水序列存在2、6～7、14～15和20年的4類準變化周期，且降水量在20世紀70年代較少，而后降水量開始增多。尹承美等[8]分析了濟南市降水年際、月、時以及強度變化特征。王燕等[9]指出，甘肅省降水自西北向東南呈增加趨勢，夏季降水量較多，降水量的年際變化不顯著。李玲萍等[10]認為，石羊河流域降水呈現逐年遞增的趨勢，夏季降水較多，冬季降水最少，且降水分布不均。馬京津等[11]研究表明，華北地區春季和夏季降水量年代際變化有一定的差別。郝立生等[12]從華北地區降水年代際變化特征出發，研究海溫、東亞夏季風、積雪等因素變化對夏季降水量的影響機制。安文芝等[13]研究表明，西北地區降水呈現從東南和西北兩側向中部減少的分布特征，東部地區降水呈逐年遞減趨勢，西部則呈增多趨勢。還有很多學者對我國近幾十年來降水氣候特征作了一系列研究得出，我國降水分布不均勻，西部、東南部以及東北地區降水較為集中，華北地區降水較少;東部地區夏季降水量呈逐年遞增趨勢，秋季和冬季降水逐年減少，西部地區冬季降水逐年增加[14-17]。

New等[18]研究表明，全球降水量除了南極以外各個區域年降水量均呈遞增趨勢。Taschetto等[19]研究發現，澳大利亞年降水量東西部差異較為顯著，西部年降水量呈遞增趨勢，而東部年降水量呈遞減趨勢。

四川省處于中國西南腹地，區域氣候表現差異較大，氣候災害較多，主要有干旱、暴雨、洪澇等。因此，研究四川省降水變化規律對防范各類氣象災害有重大意義。筆者將利用四川省36個站點逐日降水資料，采用一元線性回歸、滑動平均、M-K突變檢驗、EEMD經驗模態分解等方法，研究四川省降水演變規律。

1 資料與方法

1.1 資料

利用1980—2012年四川省對36個基準地面氣象觀測站實測逐日降水資料，分別計算上述每個站點年、月降水量平均值，用來代表四川省總體年、月降水量變化情況。季節范圍劃分如下：春季為3～5 月，夏季為6～8月， 秋季為9～11月，冬季為12月～次年2月。把日降水量超過50 mm的降水事件稱為暴雨[20]，日降水量超過25 mm的降水事件稱為大雨。

1.2 方法

1.2.1 一元線性回歸。

采用氣象傾向率反映降水量的變化趨勢，以此反映降水活動變化情況。氣候傾向率即用最小二乘法擬合線性趨勢變化：y=a+bx。式中，b為氣候傾向率，y表示降水量（mm），x表示降水量對應的時間，a為常數。當b>0時，表示降水量隨時間年限的增加呈上升趨勢;當b<0時，表示降水量隨時間年限的增加呈下降趨勢。b值的大小反映了降水量上升或下降的速率[21]。

筆者使用滑動平均對近33年四川省降水量進行趨勢擬合，用來確定降水量變化趨勢。樣本量為n的降水量序列x其滑動平均序列表示為：

1.2.2 M-K突變檢驗。

取顯著性水平a=0.05，則臨界值曲線為U0.05=±1.96。對年降水量時間序列進行M-K突變檢驗，即將UF年降水順序統計曲線、UB年降水逆序統計曲線以及臨界線繪制于同一幅圖中。如果UF順序統計曲線大于0，則表明降水量呈遞增趨勢;相反，呈遞減趨勢。若UF順序統計曲線超出0.05臨界線水平，則表明降水遞增或遞減的趨勢非常顯著。當UF順序統計曲線與UB逆序統計曲線交點處于0.05臨界線水平之內，那么交點對應的時刻便是突變開始的時間[22]。

1.2.3 EEMD經驗模態分解。

Wu等[23]首先提出EEMD分解方法。EEMD是一種新的噪聲輔助數據分析方法，該分解方法彌補了原始EMD方法的不足。算法具體步驟如下：①對原始降水資料疊加強度很低的白噪聲;

②對所疊加噪聲的降水資料進行EMD分解，得到IMF分量;

③每次疊加不同強度的白噪聲，并重復①、②N次;④將疊加N次經EMD分解得到的IMF分量求和，并取其平均值，用來代表最后的IMF分量。公式如下：

式中，cj（t）表示對原始降水資料經過EEMD經驗模態分解后得到的第j個IMF分量。

2 結果與分析

2.1 四川省降水時間變化特征

2.1.1 降水量年際變化特征。四川省總體年平均降水量為1 039.6 mm，最大年降水量為1 180.2 mm（出現于1998年），最小年降水量為878.1 mm（出現于2006年）。四川省各個縣市年平均降水量統計結果見表1。從表1可以看出，年平均降水量最多的是峨眉山地區，為1 664.3 mm;其次是雅安地區，為1 652.4 mm。四川省中部地區年降水量相對較多，東部地區年降水量其次，西部地區年平均降水量相對較少。

圖1為四川省近33年來年降水量及距平變化趨勢。圖1a中虛星線為年降水量實測值，實線為一元線性回歸值，一元回歸趨勢線擬合效果通過了0.05的顯著水平，說明擬合效果較好。近33年來四川省年降水量波動幅度較大，總體呈遞減趨勢，其氣候傾向率為-22.9 mm/10 a，表示年降水量每10年減少22.9 mm左右。圖1b為四川省近33年的年降水量距平以及滑動平均演變，柱狀圖為年降水量實測距平，虛線為5年滑動平均曲線，直線為一元線性回歸擬合值。由圖1可見，從5年滑動平均曲線可以看出，1993年是四川省年降水量的一個轉折點。1980—1993年曲線值均大于0，說明在此期間年降水量較多;1993—2012年，曲線值以小于0為主，表明在此期間四川省年降水量較少。1980—1985年年降水量變化處于相對穩定狀態，距平值以正數為主（僅1983年距平值為負），高于平均值水平，表明在此期間四川省年降水量較多。1987—1991年年降水量變化也處于相對穩定狀態，距平值均為正數，說明在此期間四川省年降水量也較多。

從四川省36個站點年降水量氣候傾向率（表2）可以看出，色達、馬爾康、稻城、康定、西昌這5個地區氣候傾向率大于0，表明每10年降水量總體呈遞增趨勢;其余站點年降水量氣候傾向率均小于0。其中，峨眉山地區年降水量遞減幅度最大，氣候傾向率為-88.9 mm/10 a，表明峨眉山地區年降水量每10年減少88.9 mm;其次是宜賓地區，年降水量每10年減少65.9 mm。

2.1.2 降水量季節變化特征。

以四川省1980—2012年的四季降水量距平為變量，以年份為自變量，依次對各個序列做趨勢，結果見圖2。圖2中柱形為季節降水量實測距平值，實線為5 a滑動平均值，虛線為一元回歸方程擬合值，擬合值均通過0.05顯著性水平檢驗。4個季節中春季降水量呈上升趨勢，夏季、秋季和冬季呈現下降趨勢。其中，夏季降水量下降趨勢最為顯著，冬季降水量下降趨勢最弱。

由圖1a可見，1980—1995年春季降水量以負距平為主，表明在此段時間內四川省境內降水量較低，且低于平均降水量。1995—2012年春季降水量以正距平為主，表明在此期間四川省境內降水量較高。降水量氣候傾向率為3.804 mm/10 a，表示春季降水量每10 a約上升3.8 mm。從5年滑動平均看，四川省1980—1995年降水量呈波動變化，沒有明顯趨勢，1987年、1994年前后降水量趨勢達到谷值，1983年、1990年前后降水量趨勢達到峰值，1995年后降水量呈現上升趨勢，一直到2006年上升趨勢消失，之后降水量呈下降趨勢。

由圖2b可見，夏季降水量在1980—1990年以正距平為主，說明在此段時間內四川省境內降水量較大，大于平均降水量。1990—2012年夏季降水量以負距平為主，說明此期間四川省境內降水量較低。降水量氣候傾向率為-18 mm/10 a，表示夏季降水量每10年下降18 mm。從5年滑動平均看，四川省1980—1995年降水量一直呈下降趨勢，在1995年下降趨勢達到谷值，1995—2000年降水量呈上升趨勢，并于2000年達到最大值，2000—2005年降水量呈下降趨勢，2005年后降水量出現上升趨勢。

由圖2c可見，秋季降水量在1980—1995年以正距平為主，說明在此段時間內四川省境內降水量較大，大于平均降水量。1990—2012年秋季降水量以負距平為主，說明此期間四川省境內降水量較低。降水量氣候傾向率為-8 mm/10 a，表示秋季降水量每10年下降8 mm。從5年滑動平均看，1980—1988年，降水量趨勢呈現先下降后上升、再上升的趨勢，其中1986年降水量達到最小值，1984年、1988年達到最大值;而1988—1997年降水量一直呈下降趨勢，在1997年達到谷值，1997—2002年降水量呈上升趨勢，2002年后呈下降趨勢，并在2008年結束。

由圖2d可見，冬季降水量在1980—1988年和1997—2003年以負距平為主，說明在此段時間內四川省境內降水量較小，1988—1997年和2003—2012年，冬季降水量以正距平為主，說明此期間四川省境內降水量較大。降水量氣候傾向率為-0.5 mm/10 a，表示冬季降水量每10年下降0.5 mm，降水量下降趨勢不明顯。從5年滑動平均看，1980—1995年，降水量整體呈上升趨勢，其中在1985年前后出現短暫的下降趨勢;1995—2001年四川省境內降水量呈下降趨勢，并在2001年達最小值;2001—2007年降水量呈上升趨勢，2007年后呈下降趨勢。

2.1.3 大雨日數變化特征。

根據全國一般定義，將25 mm≤日降水量≤50 mm的降水事件稱為大雨。筆者以時間為橫坐標，以對應年份中出現大雨的日數為縱坐標，繪制四川省大雨日數變化圖（圖3），圖中星直線為逐年統計值，實線表示大雨日數的趨勢變化，虛線表示大雨日數的擬合趨勢。

圖3表明，四川省大雨日數呈下降趨勢，下降率為4 d/10 a，大約每10年大雨日數減少4 d。四川省1980—1995年大雨日數呈下降趨勢，在1995年達到最小值;1995—2000年大雨日數呈上升趨勢，并在2000年達到最大值;2000—2012年從四川省各地區大雨日數下降率統計結果（表3）可以看出，峨眉山和雅安地區大雨日數下降率分別為9.94、9.30 d/10 a。四川省中部地區大雨日數下降率較高，東部地區其次，西部地區大雨日數下降率相對較低。

2.1.4 暴雨日數變化特征。

根據日降水量超過50 mm的降水事件為暴雨的統計標準，以時間為橫坐標，以對應年份出現暴雨的日數為縱坐標，繪制四川省暴雨日數變化趨勢圖（圖4）。圖4中星直線為逐年統計值，實線表示暴雨日數的趨勢變化，虛線表示暴雨日數的擬合趨勢。

由圖4可知，四川省暴雨日數呈下降趨勢，下降率為1.5 d/10 a，大約每10年暴雨日數隨時間會減少1.5 d。四川省1980—1995年暴雨日數呈下降趨勢，在1995年達到最小值，1995—2012年四川省暴雨日數呈上升趨勢。1980—2012年四川省暴雨日數每年平均有66 d，最多為90 d，最少為47 d。

從四川省各地區暴雨日數下降率統計結果（表4）可以看出，峨眉山和雅安地區暴雨日數下降率分別為5.55、6.42 d/10 a。四川省中部地區暴雨日數下降率較高，東部地區其次，西部地區暴雨日數下降率相對較低，幾乎不會出現暴雨情況。

2.2 四川省降水M-K突變檢驗

氣候突變是指從一個平均值狀態到另一個平均值狀態的急劇變化過程。M-K突變檢驗具有樣本不需要遵循一定的分布、也不受少數異常點干擾等優點。為了研究四川省年降水量突變情況，筆者用36個站逐年降水平均值來代表四川省年降水情況，對逐年降水資料進行M-K非參數檢驗。

從圖5a四川省全省年降水量M-K統計曲線可以看出，從1992年開始UF年降水順序統計曲線小于0，表明從1992年開始四川省年降水量開始呈遞減趨勢;從2006年開始UF年降水順序統計曲線開始超出0.05顯著性水平線，說明從2006年開始四川省年降水量呈顯著性減少趨勢;UF順序曲線和UB逆序曲線交點約為1992年，且交點處于0.05顯著性水平之內，說明交點能夠通過顯著性檢驗。因此可以得出，1992年為四川省年降水量開始突變年。

分別選取甘孜、峨眉山、巴中這3個地區的年降水量來研究四川省西部、中部、東部降水突變年情況。從圖5b甘孜年降水量M-K統計曲線可以看出，1980—1989年UF年降水順序統計曲線小于0，表明在此期間甘孜年降水量呈遞減趨勢;1990—2006年UF年降水順序統計曲線大于0，表明在此期間甘孜年降水量呈遞增趨勢;UF順序曲線和UB逆序曲線交于點1997年、2005年，且2個交點均處于0.05顯著性水平之內，說明交點能夠通過顯著性檢驗。因此，可以得出1997年、2005年為甘孜年降水量開始突變年。

從圖5c峨眉山年降水量M-K統計曲線可以看出，從1992年開始UF年降水順序統計曲線小于0，表明從1992年開始峨眉山年降水量開始呈遞減趨勢;UF順序曲線和UB逆序曲線交于點1991年，且交點均處于0.05顯著性水平之內，說明交點能夠通過顯著性檢驗。因此，可以得出1991年為峨眉山年降水量開始突變年。從圖5d巴中年降水量M-K統計曲線可以看出，從1984年開始UF年降水順序統計曲線小于0，表明從1984年開始巴中年降水量開始呈遞減趨勢;UF順序曲線和UB逆序曲線交于點1984年左右，且交點均處于0.05顯著性水平之內，說明交點能夠通過顯著性檢驗。因此，可以得出1984年為巴中年降水量開始突變年。

2.3 四川省降水EEMD經驗模態分解

使用EEMD方法，四川省1980—2012年逐年降水資料序列能夠分解成4個固有模態IMF分量和一個趨勢項分量r，表明四川省近33年來年降水量具有多個時間尺度特征。由圖6可見，固有模態第1個IMF1分量波動幅度最大。其中，20世紀80年代初期至90年代中期波動幅度相對較小，90年代中期至21世紀00年代中期波動幅度相對較大，而后波動幅度又變小，但整體表現出2～3年的準波動周期。第2個IMF2分量波動幅度相對于IMF1分量較小，但波動較長，且在20世紀90年代初期至90年代末期波動幅度較大，其余時間內波動幅度均較小，整體表現出5～10年的準波動周期。第3個IMF3分量波動幅度相對于IMF2較小，整個波動階段也存在5～10年的準周期，波動幅度在整個21世紀00年代相對較大，其余年代波動幅度較小。第4個IMF4分量波動幅度相對于其他3個IMF分量均較小，但波動最長，整個波動階段存在20年左右的準周期。由圖6i可以看出，趨勢項r在時間尺度上呈遞減趨勢，表明年降水量呈逐年遞減趨勢。這一結論與一元線性回歸方程值所求得降水氣候傾向率小于0的結論相一致。

3 結論

利用1980—2012年四川省36個站點逐日降水資料，采用一元線性回歸、M-K突變、EEMD經驗模態分解等方法，研究了四川省降水氣候時間演變規律，得出以下結論：

（1）四川省總體年平均降水量為1 039.6 mm，最大年降水量為1 180.2 mm（1998年），最小年降水量為878.1 mm（2006年）。年降水量總體呈遞減趨勢，每10年減少22.9 mm 左右。四川省中部地區年降水量相對較多，東部地區年降水量其次，西部地區年平均降水量相對較少。

（2）四川省降水量春季呈上升趨勢，降水量每10年上升3.8 mm，夏季、秋季和冬季呈現下降趨勢，降水量每10年分別減少18.0、8.0、0.5 mm。可見，四川省夏季降水量下降趨勢最為顯著，冬季降水量下降趨勢最小。

（3）四川省大雨日數呈下降趨勢，大約每10年會減少4 d，且中部地區大雨日數較多，東部地區其次，西部地區大雨日數相對較少;四川省暴雨日數統計分析得出，暴雨日數呈下降趨勢，大約每10年會減少1.5 d。

（4） 四川省年降水量從1992年開始呈遞減趨勢，其中從2006年開始年降水量呈顯著性減少趨勢，且1992年為四川省年降水量開始突變年。降水序列在年際演變過程中，存在2～3年小尺度、5～10年中尺度、20年左右大尺度這3類時間尺度變化規律。根據趨勢項得出四川省年降水量呈遞減趨勢。

參考文獻

[1] 王紹武.現代氣候學研究進展[M].北京：氣象出版社，2001：7-20.

[2] KRISHNAKUMAR K N，PRASADA RAO G S L H V，GOPAKUMAR C S.Rainfall trends in twentieth century over Kerala，India [J].Atmospheric environment，2009，43：1940-1944.

[3] DING Y H，WANG Z Y.A study of rainy seasons in China [J].Meteorology and atmospheric physics，2008，100：121-138.

[4] ZHANG Q，XU C Y，ZHANG Z，et al.Spatial and temporal variability of precipitation over China，1951-2005[J].Theoretical and applied climatology，2009，95：53-68.

[5] YUAN Y J，HE Q，ZHANG J B，et al.The primary research of climate change features and their causes for the recent 40 years in Xinjiang [J].Proceedings of SPIE，2003，4890：373-379.

[6] 周長艷，岑思弦，李躍清，等.四川省近50年降水的變化特征及影響[J].地理學報，2011，66（5）：619-630.

[7] 孫銀鳳，陸寶宏.基于EEMD的南京市降水特征分析[J].中國農村水利水電，2013（3）：5-9.

[8] 尹承美，梁永禮，冉桂平，等.濟南市區短時強降水特征分析[J].氣象科學，2010，30（2）：262-267.

[9] 王燕，王潤元，王毅榮，等.近37年甘肅省降水特征分析[J].干旱區資源與環境，2009，23（4）：94-99.

[10] 李玲萍，楊永龍，錢莉.石羊河流域近45年氣溫和降水特征分析[J].干旱區研究，2008，25（5）：705-710.

[11] 馬京津，高曉清，曲迎樂.華北地區春季和夏季降水特征及與氣候相關的分析[J].氣候與環境研究，2006，11（3）：321-329.

[12] 郝立生，丁一匯.華北降水變化研究進展[J].地理科學進展，2012，31（5）：593-601.

[13] 安文芝，祝玲敏，周文強，等.西北地區降水特征及變化規律分析[J].干旱地區農業研究，2006，24（4）：194-199.

[14] 王紹武，蔡靜寧，朱錦紅，等.中國氣候變化的研究[J].氣候與環境研究，2002，7（2）：137-145.

[15] 葉柏生，李羽中，楊大慶，等.我國過去50a來降水變化趨勢及其對水資源的影響（Ⅰ）：年系列[J].冰川凍土，2004，26（5）：587-594.

[16] 葉柏生，李羽中，楊大慶，等.我國過去50a來降水變化趨勢及其對水資源的影響（Ⅱ）：月系列[J].冰川凍土，2005，27（1）：100-105.

[17] 施雅風.中國西北氣候由暖干向暖濕轉型問題評估[M].北京：氣象出版社，2003.

[18] NEW M G，TODD M C，HULME M，et al.Precipitation measurements and trends in the twentieth century[J].International journal of climatology，2001，21（15）：1899-1922.

[19] TASCHETTO A S，ENGLAND M H.An analysis of late twentieth century trends in Australian rainfall[J].International journal of climatology，2009，29（6）：791-807.

[20] 郭元喜，王建立.山東極端降水的氣候變化[J].中國水運（學術版），2008，8（1）：175-177.

[21] 紀曉玲.寧夏雷暴天氣氣候和環流特征及典型過程分析[D].蘭州：蘭州大學，2009.

[22] 劉正源，姜蘇，曹洪亮.呼和浩特市雷暴氣候特征及其變化分析[J].江西農業學報，2013，25（2）：79-82.

[23] WU Z H，HUANG N E.Ensemble empirical mode decomposition：A noise-assited data analysis method [J].Advances in adaptive data analysis，2009，1（1）：1-41.