李金建， 劉希勝， 王銳婷， 楊 濤

(1.四川省農(nóng)業(yè)氣象中心，四川成都610072;2.青海省水文水資源勘測局，青海 西寧810001;3.四川省大氣探測中心，四川成都610072)

0 引言

降水過程是全球水分循環(huán)的重要環(huán)節(jié)，為眾多的江河提供了水分補給，其變化特征一直以來都是氣候變化、水資源等領(lǐng)域的重要研究對象［1-5］，眾多學(xué)者利用不同的方法、指標(biāo)，針對區(qū)域降水量、雨日、降水強度等要素的變化做了大量研究［6-10］。

岷江作為長江流域流量最大的支流，發(fā)源于四川省岷山南麓，整個流域均在四川境內(nèi)，也是四川省重要的河流之一，其上游5縣位于青藏高原向四川盆地的過渡地帶，一方面，較強的降水是徑流補給的重要來源;另一方面，也容易造成山洪、滑坡、泥石流等災(zāi)害，尤其是2008年汶川地震以來，該區(qū)域為地震重災(zāi)區(qū)，地震導(dǎo)致生態(tài)環(huán)境更加脆弱［11］，一遇降雨，尤其是較強的降雨過程(中雨及其以上)極易發(fā)生滑坡、泥石流等次生地質(zhì)災(zāi)害，因此，中雨以上降水具有怎樣的分配特征，變化趨勢如何等，對于區(qū)域水資源的開發(fā)利用，生態(tài)環(huán)境保護等都有重要意義，有必要進行針對性的研究。雖然，已有學(xué)者探討了岷江流域年降水以及汛期降水的變化特征［12-13］，但對于中雨以上量級的降水特征的研究不多，基于此，采用1966～2012年47年的氣象觀測資料，分析岷江上游地區(qū)中雨以上降水變化特征，旨在為該地區(qū)防澇抗旱、水資源合理利用等提供科學(xué)參考依據(jù)。

1 數(shù)據(jù)來源與研究方法

1．1 數(shù)據(jù)來源

研究采用四川省氣象局信息中心提供的岷江上游5縣氣象站1966～2012年逐日降水觀測數(shù)據(jù)，在逐日降水資料基礎(chǔ)上，通過計算分析得到每個站點建站以來的中雨以上降水序列(含降水量及降水日數(shù))。

1．2 研究方法

首先利用Mann-Kendall(MK)法［13］對岷江上游地區(qū)中雨以上雨量、雨日進行突變檢驗，尋求可能的突變點，然后應(yīng)用滑動 t檢驗(MTT)法［14］和 Yamamoto法［14］，對可能的突變點進行驗證，增強突變分析結(jié)果的可信度。

采用Morlet小波變換分析岷江上游地區(qū)中雨以上雨量、雨日的周期變化特征。Morlet小波變換原理見文獻［15-16］。在進行小波變換時，為了減弱“頭部影響”在序列兩端產(chǎn)生的誤差，利用對稱外延法［15］對時間序列進行擴展。

2 結(jié)果與分析

2．1 年內(nèi)分配特征

對岷江上游5縣(茂縣、松潘、汶川、理縣和黑水)氣象站1966～2012年年總降水量及年中雨以上降水量進行統(tǒng)計發(fā)現(xiàn)，各站降水量及中雨以上降水量均有較大的差異，從多年平均狀況來看，5站多年平均降水量分別為471.7 mm、719.3 mm、502.9 mm、611.6 mm和822.3 mm;中雨以上降水量則分別為185 mm、330 mm、198.5 mm、244 mm和445.9 mm。盡管數(shù)值上有較大差異，但是從中雨以上降水量占全年總降水量的比重來看，各均達(dá)到較高的比重，其中茂縣最小，也達(dá)到39.2%，黑水最大，可達(dá)54.2%。由此可見，中雨以上降水量對于全年降水變化具有重要的作用，研究其變化特征無論對降水，還是河流徑流變化都具有重要意義。

岷江上游5縣多年平均的中雨以上降水量及雨日數(shù)的月變化趨勢圖如圖1所示。從圖1可以看出，降水量和雨日具有相似的分配特征，1～3月及11～12月降水量和雨日數(shù)都比較少，從4月份開始增加，但每個站之間又略有差異，其中茂縣、理縣和汶川均呈單峰型，最大平均降雨量分別出現(xiàn)在8月、6月和8月，而最大平均雨日數(shù)則分別出現(xiàn)在6月、6月和8月;松潘和黑水均表現(xiàn)為雙峰型，松潘最大平均降雨量和最大平均雨日數(shù)均在9月達(dá)到最大值，次大值出現(xiàn)在6月(其中5、6月份的平均雨日數(shù)相同)，黑水最大平均降雨量和最大平均雨日數(shù)均在6月達(dá)到最大值，9月次之。

圖1 岷江上游中雨以上降水量及雨日數(shù)月變化

2．2 趨勢變化特征

岷江上游5縣中雨以上雨量、雨日的逐年變化、5年滑動平均以及線性變化趨勢如圖2所示。從圖2可以看出，就整體變化趨勢而言，5個站的雨量和雨日的年際變化趨勢都不顯著(除理縣的雨量、雨日以及黑水的雨日呈微弱上升趨勢外，其余均表現(xiàn)為微弱的下降趨勢，均未通過0.05顯著度檢驗)，但是，在不同的時間段則呈現(xiàn)不同的變化趨勢。從5年滑動平均曲線可以看出，茂縣雨量(圖2a)、雨日(圖2b)均在20世紀(jì)70年代末以前、80年代中期至90年代初及21世紀(jì)初以來呈現(xiàn)上升趨勢，減少趨勢則主要出現(xiàn)在20世紀(jì)80年代初及90年代;松潘雨量(圖2c)、雨日(圖2d)在20世紀(jì)60年代末至70年代初、80年代中期至90年代初以及21世紀(jì)初以來呈現(xiàn)上升趨勢，減少趨勢主要出現(xiàn)在70年代初至80年代中期及90年代初至90年代中期，而90年代末至21世紀(jì)初則變化不明顯;汶川雨量(圖2e)、雨日(圖2f)則在20世紀(jì)60年代末70年代初、80年代初中期以及90年代中期至21世紀(jì)10年代末表現(xiàn)為較明顯的下降趨勢，而在20世紀(jì)70年代中期至80年代初、80年代后期至90年代初及21世紀(jì)10年代末以來呈現(xiàn)明顯的上升趨勢;理縣雨量(圖2g)、雨日(圖2h)則在20世紀(jì)60年代末至70年代初、70年代后期、80年代中至90年代初及21世紀(jì)10年代中期以來呈上升趨勢，而在70年代中期、80年代初以及90年代呈現(xiàn)下降趨勢;黑水雨量(圖2i)、雨日(圖2j)上升趨勢主要出現(xiàn)在20世紀(jì)60年代末至70年代初、80年代中期至90年代初及21世紀(jì)以來，下降趨勢則主要出現(xiàn)在20世紀(jì)80年代初以及90年代。

岷江上游5縣中雨以上降水年際年代際變化特征與川西高原［17］乃至盆地西部［18］夏季降水變化具有較好的一致性，而其特殊的地理位置，使降水變化受到多種環(huán)流系統(tǒng)的綜合影響［19］，研究表明，當(dāng)研究區(qū)及周邊地區(qū)降水處于偏多階段時，東亞中高緯地區(qū)多呈現(xiàn)較強的兩脊一槽環(huán)流型，同時西太副高整體偏北;而當(dāng)降水處于偏少階段時，東亞中高緯地區(qū)很少出現(xiàn)或者僅表現(xiàn)為較弱的兩脊一槽環(huán)流型，貝湖以南地區(qū)的槽區(qū)減弱甚至為脊區(qū)，同時西太副高整體偏南［17-18］;20世紀(jì)80年代末以來，夏季副熱帶高壓的強度偏強，西界偏西，北界偏南，與西南大部分地區(qū)的降水減少有重要聯(lián)系［20］。

圖2 各站雨量、雨日逐年變化曲線

為進一步探明各個站點降水的共性變化特征，對岷江上游5縣中雨以上雨量、雨日進行EOF分析，其前3個模態(tài)的方差貢獻及累積方差貢獻率如表1所示，結(jié)果表明兩者的第1主成分的方差貢獻率分別為49.53%和46.78%，且相應(yīng)的荷載值均為正值，同時，無論雨量還是雨日的第1主成分時間系數(shù)與各站的相關(guān)系數(shù)也均通過0.001的信度檢驗(表2)，表明第1主成分時間系數(shù)的演變能夠反映區(qū)域雨量、雨日的共性變化特征。

表1 中雨以上雨量及雨日前3個特征向量方差貢獻率

表2 第1主成分時間系數(shù)與各站點雨量、雨日相關(guān)系數(shù)

圖3給出雨量、雨日的第1主成分時間系數(shù)(PC1)變化特征，由圖3可以發(fā)現(xiàn)，與各站點的情況相似，雨量PC1與雨日PC1變化具有非常高的一致性，從長期變化趨勢看，區(qū)域整體呈微弱的下降趨勢，但不同階段，仍表現(xiàn)比較明顯的波動，從圖3的5年滑動平均曲線可以看出，區(qū)域雨量、雨日均在20世紀(jì)60年代末至70年代初、90年代中至21世紀(jì)初處于偏少階段，而在70年代中至80年代初、80年代末至90年代中表現(xiàn)為偏多階段。這些特征也是各個站點序列表現(xiàn)出來共性特點，這種結(jié)果進一步顯示了利用第1主成分時間系數(shù)表征區(qū)域變化特征的可行性。

圖3 區(qū)域雨量PC1及雨日PC1變化曲線

2．3 突變特征

前面利用EOF分析，得到代表整個區(qū)域雨量、雨日變化曲線，即雨量PC1及雨日PC1，下面將主要利用這兩條序列分析區(qū)域變化特征。

圖4 區(qū)域雨量、雨日MK檢驗

首先利用MK法對雨量PC1及雨日PC1進行突變檢測，結(jié)果如圖4所示，兩者具有十分相似的變化特征，兩者均在20世紀(jì)70年代中至90年代的大部分時段UF值＞0，表明雨量、雨日有增加的趨勢，而21世紀(jì)以來UF值基本＜0，雨量、雨日有減少的趨勢;UF和UB曲線在±1.96信度線之間有3個交點，分別出現(xiàn)在1995年、2008年及2011年，但后2個交點(2008年和2011年)，只能視為同一個時期，必然只能保留一個突變點，綜合考慮氣候突變意義、序列長度等因素，保留2008年的突變點。

為進一步的驗證MK方法檢測的突變點，對區(qū)域可能突變年前后不同時段作滑動t檢驗(MTT)以及信噪比檢驗(Yamamoto)，見表3。分別取 n1=n2=2、3、4，對區(qū)域雨量而言，1995年的突變點，利用MTT檢驗，當(dāng) n1=n2=3時通過0.05的信度檢驗;利用Yamamot檢驗，則當(dāng)n1=n2=2、3時通過0.01的信度檢驗。2008年的突變點，利用MTT檢驗，當(dāng)n1=n2=2、4時通過0.01的信度檢驗，當(dāng)n1=n2=3時通過0.05的信度檢驗;利用Yamamot檢驗，則當(dāng)n1=n2=2、3時通過0.01的信度檢驗，當(dāng)n1=n2=4時通過了0.05的信度檢驗。對區(qū)域雨日而言，1995年的突變點，利用MTT檢驗，當(dāng)n1=n2=2、3時通過了0.05的信度檢驗;利用Yamamot檢驗，則當(dāng)n1=n2=2、3時通過了0.01的信度檢驗。2008年的突變點，利用MTT檢驗，當(dāng)n1=n2=3、4時通過0.05的信度檢驗，當(dāng)n1=n2=3時通過0.05的信度檢驗;利用Yamamot檢驗，則當(dāng)n1=n2=3時通過0.01的信度檢驗，當(dāng)n1=n2=3、4時通過0.05的信度檢驗。綜上，通過兩種方法的佐證，這兩個突變點具有一定的可信性。

表3 突變點的信度檢驗(MTT、Yamamot)

大量研究表明，20世紀(jì)以來，氣候存在明顯的突變現(xiàn)象，在20世紀(jì)初、20世紀(jì)20年代、60年代和90年代是突變現(xiàn)象出現(xiàn)的主要時段［21-25］，但是，由于局地性差異，不同地區(qū)突變點出現(xiàn)的時間也存在一定的差異［26］，檢測的1995年突變點與90年代出現(xiàn)的大范圍的氣候突變相符，而2008年的突變點還有待更多研究進一步驗證。

2．4 周期變化分析

將岷江上游區(qū)域雨量PC1、區(qū)域雨日PC1進行Morlet小波變換，得到小波變換實部，如圖5所示。由圖可以看出，兩者具有一致的周期變化特征，周期變化明顯，干濕交替顯著。10年以上的年代際變化主要存在準(zhǔn)18年周期振蕩和準(zhǔn)12年周期振蕩，其中準(zhǔn)18年周期表現(xiàn)顯著且穩(wěn)定，貫穿整個時域，在過去47年間始終存在，且主周期幾乎不變，強度較強，主要經(jīng)歷5次干濕更替，其轉(zhuǎn)換點分別為1970年、1980年、1989年、1998年和 2008年，其中 1966～1970年、1981～1989年及1999～2007年為偏干(或雨日偏少)階段，1971～1980年、1990～1998年及2008～2012年為偏濕(或雨日偏多)階段。從這個時間尺度來看，研究區(qū)域目前正處于降水偏多時期，且未來一段時間內(nèi)仍將呈偏多的狀態(tài)。而準(zhǔn)12年周期振蕩則在20世紀(jì)80年代中后期逐漸轉(zhuǎn)化為準(zhǔn)10年周期振蕩。10年以下尺度周期變化特征則已2年為中心的周期與4年為中心的周期的轉(zhuǎn)化，20世紀(jì)80年代初以前主要表現(xiàn)為2年為中心的周期，80年代中其分化為2年周期和4年周期這兩個主振蕩周期，自20世紀(jì)80年代中期開始4年為中心的周期振蕩顯著且穩(wěn)定的存在，到目前為止，該周期振蕩由7個偏干階段和6個偏濕階段組成;而自80年代中期分化以來，2年周期則表現(xiàn)出時強時弱的態(tài)勢，其在20世紀(jì)90年代至21世紀(jì)初信號較強，其他時段信號較弱，甚至與4年周期相融合。

圖5 區(qū)域雨量、雨日小波變換實部圖

3 結(jié)論

利用統(tǒng)計學(xué)方法，從年內(nèi)分配特征、趨勢分析、突變及周期變化等角度，對岷江上游地區(qū)中雨以上降水量及雨日進行分析，主要得到以下結(jié)論:

(1)岷江上游5縣多年平均的中雨以上降水量以及雨日具有相似的年內(nèi)分配特征，均表現(xiàn)為1～3月及11～12月降水量和雨日數(shù)較少，而4～9月為較為集中期。

(2)利用EOF分析，獲取代表區(qū)域中雨以上降水量以及雨日的序列，即第1主成分時間系數(shù)(雨量PC1與雨日PC1)，表明區(qū)域中雨以上降水量以及雨日均呈微弱的下降趨勢，其5年滑動平均曲線表明20世紀(jì)60年代末至70年代初、90年代中至21世紀(jì)初處于降水(雨日)偏少階段，而在70年代中至80年代初、80年代末至90年代中表現(xiàn)為降水(雨日)偏多階段。

(3)通過MK法，并以滑動t檢驗法及YAMAMOT法為佐證，發(fā)現(xiàn)區(qū)域中雨以上降水量以及雨日在1995年和2008年存在突變現(xiàn)象。

(4)區(qū)域雨量、雨日具有一致的周期變化特征，干濕交替顯著。10年以上尺度存在2個主要振蕩模態(tài)，中心周期分別為12年、10年;10年以下的年際變化也存在2個主要振蕩模態(tài)，中心周期分別為4年、2年。

［1］ 牛存穩(wěn)，張利平，夏軍.華北地區(qū)降水量的小波分析［J］.干旱區(qū)地理，2004，27(1):66-70.

［2］ 匡正，李仲貞，林一驊.華北降水時間序列資料的小波分析［J］.氣候與環(huán)境研究，2000，5(3):312-317.

［3］ 左洪超，呂世華，胡隱樵.中國近50a氣溫及降水量的變化趨勢分析［J］.高原氣象，2004，23(2):238-244.

［4］ 李少娟，何大明，張一平.縱向嶺谷區(qū)降水量時空變化及其地域分異規(guī)律［J］.科學(xué)通報，2007，52:51-63.

［5］ 楊輝，宋正山.華北地區(qū)水資源多時間尺度分析［J］.高原氣象，1999，18(4):496-508.

［6］ 張琪，李躍清.近48年西南地區(qū)降水量和雨日的氣候變化特征［J］.高原氣象，2014，33(2):372-383.

［7］ 陳少勇，任燕，喬立，等.中國西北地區(qū)大雨以上降水日數(shù)的時空分布特征［J］.資源科學(xué)，2011，33(5):958-965.

［8］ 孫善磊，周鎖銓，金博，等.淮海地區(qū)降水周期及突變特征分析［J］.氣象科學(xué)，2010，30(2):221-227.

［9］ Huang J，Sun S，Xue Y，et al.Spatial and Temporal Variability of Precipitation and Dryness/Wetness During 1961-2008 in Sichuan Province［J］.West China.Water Resources Management，2014，28(6):1655-1670.

［10］ 蘇布達(dá)，姜彤，任國玉，等.長江流域極端降水時空分布和趨勢［J］.氣候變化研究進展，2006，2(1):9-14.

［11］ 徐新良，江東，莊大方，等.汶川地震災(zāi)害核心區(qū)生態(tài)環(huán)境影響評估［J］.生態(tài)學(xué)報，2008，28(12):5899-5908.

［12］ 周大杰，張明珠，俞火亙，等.岷江流域汛期降水時空特征研究［J］.水文，2009，29(4):26-33.

［13］ 張明珠，俞火亙，宋新山，等.岷江上游及青衣江流域近50年來降水變化特征分析［J］.南水北調(diào)與水利科技，2008，6(3):46-53.

［14］ 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)［M］.北京:氣象出版社，2007:10-116.

［15］ 吳洪寶，吳蕾.氣候變率診斷和預(yù)測方法［M］.北京:氣象出版社，2005:33-245.

［16］ 林振山，鄧自旺.子波氣候診斷技術(shù)的研究［M］.北京:氣象出版社，1999:1-38.

［17］ 陳權(quán)亮，劉曉冉，范廣洲，等.川西高原夏季降水變化特征及其異常年環(huán)流形勢［J］.中國沙漠，2010，30(3):706-711.

［18］ 周長艷，李躍清，房靜，等.高原東側(cè)川渝盆地東西部夏季降水及其大尺度環(huán)流特征［J］.高原山地氣象研究，2008，28(2):1-9.

［19］ 徐裕華.西南氣候［M］.北京:氣象出版社，1991:201-312.

［20］ 慕巧珍，王紹武，朱錦紅，等.近百年夏季西太平洋副熱帶高壓的變化［J］.大氣科學(xué)，2001，25(6):787-797.

［21］ 符淙斌.氣候突變現(xiàn)象的研究［J］.大氣科學(xué)，1994，18(3):373-384.

［22］ 孫秀忠，羅勇，張霞，等.近46年來我國降雪變化特征分析［J］.高原氣象，2010，29(6):1594-1601.

［23］ 趙芳芳，徐宗學(xué).黃河蘭州以上氣候要素長期變化趨勢和突變特征分析［J］.氣象學(xué)報，2006，64(2):246-255.

［24］ 李金建，李媛媛，秦寧生，等.利用樹木年輪重建川西松潘高原5月降水變化特征［J］.高原氣象，2013，32(4):1195-1202.

［25］ 王素萍，宋連春，韓永翔，等.近50年河西走廊地區(qū)降水均值突變特征分析［J］.高原氣象，2011，30(5):1286-1293.

［26］ 馬振鋒，彭駿，高文良，等.近40年西南地區(qū)的氣候變化事實［J］.高原氣象，2006，25(4):633-642.