苗婷+黃曉梅+蘇軍鋒+魏邦憲

摘要 利用甘肅省康縣1971—2014年月報表中日降水?dāng)?shù)據(jù)和日照數(shù)據(jù)，運用線性傾向估計、累積距平法、Mann-Kendall突變檢驗法、滑動T檢驗法和小波分析法，分析了康縣連陰雨時間變化特征、突變情況和變化周期。結(jié)果表明：康縣連陰雨在時間尺度上呈“兩降兩升”型變化。康縣連陰雨從1992年開始突變減少，到2007年左右達(dá)到顯著性水平。連陰雨在變化過程中存在著3～7年、8～16年、17～26年以及27年以上4類尺度的周期變化規(guī)律，21年為第一主周期。

關(guān)鍵詞 連陰雨；累積距平法；突變檢驗；小波分析；甘肅康縣；1971—2014年

中圖分類號 P426.61 文獻(xiàn)標(biāo)識碼 A 文章編號 1007-5739（2017）06-0238-03

康縣地處甘肅省東南部，屬亞熱帶向暖溫帶過渡地區(qū)，氣候溫和，四季分明，雨水充沛，日照較少，災(zāi)害性天氣頻繁發(fā)生[1]。主要災(zāi)害性天氣有暴雨、干旱、雷暴、寒潮、大風(fēng)、倒春寒及連陰雨等。根據(jù)季節(jié)的不同連陰雨會給農(nóng)業(yè)生產(chǎn)造成不同的影響。1987年夏季連陰雨對康縣冬小麥的成熟和收割打碾造成嚴(yán)重影響，使冬小麥出現(xiàn)不同程度的發(fā)芽、霉變、黑穗等現(xiàn)象，導(dǎo)致冬小麥產(chǎn)量和品質(zhì)下降；1989年春季低溫連陰雨造成冬小麥生長發(fā)育期比往年推遲10 d左右，后期氣溫回升慢，誘發(fā)小麥病蟲害；2009年秋季連陰雨天氣影響康縣經(jīng)濟(jì)作物秋收，導(dǎo)致核桃等經(jīng)濟(jì)林果不能按時采收，品質(zhì)下降，減產(chǎn)明顯。此外，雨量大、持續(xù)時間長的連陰雨也會造成洪澇災(zāi)害、山體滑坡，對交通運輸業(yè)及人們生命財產(chǎn)安全產(chǎn)生極大影響。

國內(nèi)外氣象工作者對連陰雨的研究主要包括2個方向：一是基于預(yù)報的連陰雨產(chǎn)生的環(huán)流背景及影響因素，二是基于測報的連陰雨時空分布特征及其對農(nóng)作物和交通運輸?shù)鹊挠绊慬2-25]。方建剛等[26]利用NCEP/NCAR再分析資料，分析了陜西省2011年秋季強(qiáng)連陰雨期間歐亞大氣環(huán)流的異常特征；王建兵等[27]分析了甘南高原秋季連陰雨的氣候特征及主要環(huán)流形勢，得出副高偏強(qiáng)有利西南氣流把孟加拉灣水汽輸送到甘南，對連陰雨天氣的形成有重要作用；潘 旸等[28]分析了2009年2—3月我國南方一次連陰雨過程的環(huán)流背景，得出冷空氣和南方暖濕氣流形成的穩(wěn)定準(zhǔn)靜止鋒，是影響這次連陰雨過程的主要天氣系統(tǒng)；張 智等[29]運用多種方法分析了寧夏連陰雨氣候變化特征，得出連陰雨發(fā)生次數(shù)的氣候變化特征及突變、周期變化；王 棟等[30]分析了1960年以來山西秋季連陰雨的變化趨勢和時空結(jié)構(gòu)特征，并建立了秋季連陰雨強(qiáng)度指數(shù)模型；魏 鋒等[31]分析甘肅省近35年連陰雨天氣氣候特征，得出甘肅省連陰雨的分布特點是自西北向東南呈階梯狀增加。

本文以甘肅省康縣為例，利用近44年康縣氣象月報表中的日降水?dāng)?shù)據(jù)和日照數(shù)據(jù)，運用多種方法分析了年際、季節(jié)、月際尺度上連陰雨的時間分布特征、突變和周期變化。利用線性傾向估計和累積距平法得到了多種尺度上連陰雨的時間分布特征，利用Mann-Kendall突變檢驗和滑動T檢驗法得到連陰雨的突變時間，避免了一種方法導(dǎo)致的虛假突變，利用Morlet小波分析法得到了康縣年際連陰雨的周期變化。通過以上方法，希望在精細(xì)化的大趨勢下為康縣長期預(yù)報災(zāi)害性天氣和防災(zāi)減災(zāi)提供依據(jù)，減少連陰雨造成的農(nóng)業(yè)、交通業(yè)和人民生產(chǎn)生活的損失。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)來源

利用甘肅省康縣1971—2014年氣象月報表中的日降水?dāng)?shù)據(jù)和日照數(shù)據(jù)，統(tǒng)計出近44年逐月連陰雨次數(shù)，并得到年代際逐月連陰雨總次數(shù)（表1）。根據(jù)本地天氣氣候特點及對農(nóng)業(yè)生產(chǎn)影響的情況，各地對連陰雨天氣的規(guī)定略有差異。本文所用標(biāo)準(zhǔn)為甘肅省氣象局所定，即連陰雨需同時滿足以下條件：連續(xù)陰雨日數(shù)在5 d或者5 d以上，總降水過程雨量≥15 mm，期間單天日照≤10 h，過程的開始、結(jié)束日雨量≥0.1 mm（允許過程中有1 d微量或無降水），即稱為一次連陰雨天氣過程。將連陰雨按季節(jié)劃分，12月至次年2月為冬季連陰雨，3—5月為春季連陰雨，6—8月為夏季連陰雨，9—11月為秋季連陰雨。

1.2 研究方法

1.2.1 線性傾向估計法。線性傾向估計是一種特殊的線性回歸形式，用一條合理的直線yi=a+bti（i=1，2，…，n）表示連陰雨次數(shù)y與其時間t之間的關(guān)系。其中，回歸系數(shù)b為傾向值，傾向值正（負(fù)）表明連陰雨隨時間呈上升（下降）趨勢[32]。

1.2.2 累積距平法。累積距平是氣象上一種常用的、由曲線直觀判斷變化趨勢的方法。累積距平曲線呈上升趨勢，表示距平值增加；呈下降趨勢，則表示距平值減小。從曲線明顯上下起伏可以判斷，其長期顯著的變化趨勢及持續(xù)性變化，從曲線小的波動變化考察其短期的距平值變化[32]。

1.2.3 Mann-Kendall突變檢驗法。Mann-Kendall突變檢驗法是一種非參數(shù)統(tǒng)計檢驗方法，也稱無分布檢驗，其優(yōu)點是不需要序列遵循一定的分布，也不受少數(shù)異常值的干擾，更適用于順序變量的分析，計算也較簡單，而且可以明確突變開始的時間，并指出突變區(qū)域，是一種常見的突變檢驗方法[32]。

1.2.4 滑動T檢驗?；瑒覶檢驗是一種均值突變，通過考察2組樣本平均值的差異是否顯著來檢驗突變，這種方法子序列的選擇具有一定的人為性，需要反復(fù)變動子序列的長度以提高計算結(jié)果的可靠性[32]。

1.2.5 Morlet小波變換。Morlet小波變換是連續(xù)復(fù)小波變換，可以同時給出時間序列變化的位相和振幅兩方面的信息，并消除用實小波變換系數(shù)作為判據(jù)而產(chǎn)生的虛假振蕩，揭示出隱藏在時間序列中的多種變化周期，充分反映時間序列在不同時間尺度中的變化趨勢，并能對時間序列未來發(fā)展趨勢進(jìn)行定性估計[33-34]。

2 結(jié)果與分析

2.1 時間變化特征

利用康縣1971—2014年年際連陰雨數(shù)據(jù)進(jìn)行線性傾向估計，并求得近44年康縣年際連陰雨距平值（圖1）。可以看出，康縣年際連陰雨發(fā)生次數(shù)呈減少趨勢，但未通過T檢驗，說明這種減少趨勢不明顯。連陰雨累積距平曲線呈“兩降兩升”型變化，存在明顯的偏多年和偏少年。其中，20世紀(jì)70年代到80年代初期連陰雨呈增加趨勢，1977年以前呈上升趨勢，趨勢明顯，1978—1984年呈波動上升，趨勢較為緩慢，總體上康縣進(jìn)入連陰雨多發(fā)階段；80年代中期有一個明顯的短暫下降趨勢，康縣進(jìn)入連陰雨少發(fā)階段；80年代末期到90年代末期，康縣連陰雨次數(shù)呈直線上升趨勢，進(jìn)入多發(fā)階段；90年代末至今呈波動減少趨勢，其中90年代末呈直線下降趨勢，21世紀(jì)初開始呈波動減少趨勢，康縣連陰雨進(jìn)入少發(fā)階段。

2.2 突變檢驗分析

本文采用Mann-Kendall法和滑動T檢驗共同分析康縣連陰雨的突變情況。用Mann-Kendall法檢驗1971—2014年康縣年際、季節(jié)、月際連陰雨序列的突變，給定顯著性水平α=0.01（即μ0.01=±2.57）。由圖2可知，近44年康縣連陰雨呈波動下降趨勢，在20世紀(jì)70年代存在突變點，但交點過多，不能確定突變點的正確年份；1992年接近突變點，但未突變。經(jīng)反復(fù)變化步長得出，取步長為5的滑動T檢驗法驗證突變點的位置最為理想。從圖3可以看出，通過0.01顯著性水平的突變點為1992年和1996年2處，可以證明Mann-Kendall法檢驗出的20世紀(jì)70年代的幾個突變點為虛假突變?；瑒覶檢驗出的1996年在Mann-Kendall法中沒有體現(xiàn)出來，故真正的突變點為1992年。從這一年開始，康縣連陰雨突變減少，到2007年左右達(dá)到顯著性水平。

2.3 小波周期分析

本文利用Morlet小波系數(shù)實部等值線圖反映連陰雨不同時間尺度的周期變化及其在時間域中的分布，進(jìn)而判斷在不同時間尺度上連陰雨的周期和未來變化趨勢。去除邊界效應(yīng)，由圖4可看出，連陰雨年際變化過程中存在著3～7年、8～16年、17～26年以及27年以上4類時間尺度的周期變化規(guī)律。其中，在17～26年尺度上出現(xiàn)了多發(fā)—少發(fā)交替的準(zhǔn)3次振蕩，并在整個分析時段表現(xiàn)穩(wěn)定，具有全域性；8～16年尺度上2000年以前振蕩比較穩(wěn)定，2001年開始波動較大；3～7年和37年以上尺度上規(guī)律不明顯。

小波方差圖能反映連陰雨時間序列的波動能量隨時間尺度的分布情況，可用來確定連陰雨變化過程中存在的主周期。從圖5可以看出，康縣年際連陰雨存在3個較為明顯的峰值，它們依次對應(yīng)21年、10年和7年的時間尺度。其中，最大峰值對應(yīng)21年的時間尺度，說明21年左右的周期振蕩最強(qiáng)，為年際連陰雨變化的第一主周期；10年時間尺度對應(yīng)著第二峰值，為年際連陰雨變化的第二主周期；連陰雨的第三主周期對應(yīng)著7年的時間尺度。這說明上述3個周期的波動控制著年際連陰雨在整個時間域內(nèi)的變化特征。圖中在28年以后有上升趨勢，但未出現(xiàn)拐點。由于本文數(shù)據(jù)僅為近44年，能否得出其他更大尺度的周期，有待后續(xù)研究。

3 結(jié)論

研究結(jié)果表明：康縣年際連陰雨總體呈現(xiàn)減少趨勢，累積距平曲線呈“兩降兩升”型變化，存在明顯的偏多年和偏少年；通過Mann-Kendall法和滑動T檢驗共同檢驗康縣連陰雨序列的突變可知，康縣年際連陰雨從1992年開始突變減少，到2007年左右達(dá)到顯著性水平。通過Molet小波分析，康縣年際連陰雨變化過程中存在3～7年、8～16年、17～26年以及27年以上4類時間尺度的周期變化規(guī)律，并存在3個較為明顯的峰值，連陰雨在28年以后有上升趨勢，但由于本文數(shù)據(jù)僅為近44年，能否得出其他更大尺度的周期，有待后續(xù)研究。

4 參考文獻(xiàn)

[1] 韓建熙，苗婷，王大偉.2015年康縣氣候影響評價[J].安徽農(nóng)業(yè)科學(xué)，2016（11）：206-208.

[2] 張熙，李芳，趙京峰，等.2014年9月濟(jì)寧地區(qū)連陰雨降水特征分析[J].現(xiàn)代農(nóng)業(yè)科技，2015（23）：250-251.

[3] 李恩莉，邊赟.寶雞地區(qū)秋季連陰雨的氣候特征及強(qiáng)連陰雨年的環(huán)流特征[J].現(xiàn)代農(nóng)業(yè)科技，2016（2）：238-240.

[4] 馬玉坤，李菲，韋伯龍，等.臨夏州一次連陰雨天氣過程診斷分析[J].現(xiàn)代農(nóng)業(yè)科技，2016（3）：254-256.

[5] 魏秀蘭，趙熙，黃玉芳.2010年和2011年魯西南初秋罕見連陰雨天氣對比分析[J].干旱氣象，2015（4）：626-634.

[6] 周鵬康，秦金梅，孫秀芬，等.云南1981—2010年5—10月連陰雨時空分布特征[J].干旱氣象，2015（5）：796-801.

[7] 席林華.春季低溫連陰雨天氣的熱帶環(huán)流和總云量分布[J].氣象，1989（8）：10-15.

[8] 李琳，馮艷霞，王弢.春播期連陰雨天氣成因分析[J].現(xiàn)代農(nóng)業(yè)科技，2013（7）：268.

[9] 孫美蓮，楊佳，張憲冬，等.2007年甘肅省連陰雨天氣過程分析[J].現(xiàn)代農(nóng)業(yè)科技，2013（9）：255.

[10] 佘錄青，高宕鳳.成縣連陰雨發(fā)生規(guī)律及其對農(nóng)作物的影響[J].現(xiàn)代農(nóng)業(yè)科技，2013（11）：260-261.

[11] 孫昕，王丹.2011年秋季連陰雨對商水縣農(nóng)業(yè)生產(chǎn)的影響[J].現(xiàn)代農(nóng)業(yè)科技，2012（8）：299.

[12] 李娜.2014年濟(jì)寧市主要農(nóng)業(yè)氣象災(zāi)害及其影響分析[J].現(xiàn)代農(nóng)業(yè)科技，2015（4）：236-237.

[13] 常靜，羅兵，常自祎，等.2014年4月14—21日臨夏州春季連陰雨過程分析[J].現(xiàn)代農(nóng)業(yè)科技，2015（22）：232-233.

[14] 吳洪顏，高蘋，趙凱.春季連陰雨對江蘇省夏收作物產(chǎn)量的影響[J].災(zāi)害學(xué)，2003（3）：48-51.

[15] 姜愛軍，田心如，王冰梅.連陰雨災(zāi)害評估模型的研究[J].災(zāi)害學(xué)，1997（2）：49-53.

[16] 姚靜，高慶九，俞小鼎，等.一次農(nóng)業(yè)災(zāi)害嚴(yán)重的連陰雨機(jī)理分析[J].現(xiàn)代農(nóng)業(yè)科技，2016（4）：247.

[17] 黃德珍，任淑華，杜長林.秋季連陰雨對農(nóng)業(yè)生產(chǎn)的影響[J].現(xiàn)代農(nóng)業(yè)科技，2010（4）：322.

[18] 胡江波，蘇俊輝.漢中市2009年8月16—22日連陰雨、暴雨天氣過程分析[J].現(xiàn)代農(nóng)業(yè)科技，2010（6）：282-283.

[19] ZHANG Jingzhen，ZHANG Li，LI Jianming，et al.Analysis on the Variat-ion Characteristics of Precipitation in North China during Recent 30 Years[J].Meteorological and Environmental Research，2010，10：20-24.

[20] JING Anhua，ZHANG Zonghao，KONG Fanzhong.Climatic Analysis of Continuous Rain Southwest of Shandong[J].Meteorological and Environ-mental Research，2010（9）：15-17.

[21] ZHAI P M，SUN A J，REN F M，et al.Changes of climate extremes in China[J].Climate Change，1999，42（1）：203-218.

[22] 張秀雯.上海地區(qū)早秋連陰雨的中期預(yù)報方法[J].氣象，1979（12）：26-28.

[23] 鄒浩，陸善焌，樓世博，等.用模糊聚類分析方法作春季連陰雨預(yù)報[J].氣象學(xué)報，1983（2）：231-235.

[24] 項瑛，程婷，王可法，等.江蘇省連陰雨過程時空分布特征分析[J].氣象科學(xué)，2011，31（增刊1）：36-39.

[25] 宋連春，楊興國，韓永翔，等.甘肅氣象災(zāi)害與氣候變化問題的初步研究[J].干旱氣象，2006，24（2）：63-69.

[26] 方建剛，易俊蓮.2011年陜西秋季強(qiáng)連陰雨天氣的大氣環(huán)流異常特征分析[J].災(zāi)害學(xué)，2013（1）：6-10.

[27] 王建兵，安華銀，汪治桂，等.甘南高原秋季連陰雨的氣候特征及主要環(huán)流形勢[J].干旱氣象，2013（1）：70-77.

[28] 潘旸，李建，廖捷，等.2009年2—3月我國南方連陰雨天氣過程分析[J].氣象，2010（3）：39-46.

[29] 張智，梁培，陳玉華，等.寧夏連陰雨氣候變化特征分析研究[J].災(zāi)害學(xué)，2010（1）：69-72.

[30] 王棟，譚桂容，耿新.1960年以來山西秋季連陰雨的氣候特征分析[J].災(zāi)害學(xué)，2015（1）：75-81.

[31] 魏鋒，白虎志，孫秉強(qiáng).甘肅省近35年連陰雨天氣氣候特征分析[J].成都信息工程學(xué)院學(xué)報，2005（4）：479-482.

[32] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)[M].北京.氣象出版社，2007：37-65.

[33] 廖雪萍，覃衛(wèi)堅，唐炳莉，等.廣西近50年暴雨日數(shù)變化的小波分析[J].氣象，2007（12）：39-45.

[34] 張家誠.中國氣候[M].上海：上?？萍汲霭嫔?，1985.