陸琴琴 汪翔
蚌埠市汛期降水多時(shí)間尺度分析
陸琴琴 汪翔
安徽省蚌埠市地處淮河沿岸,屬暖溫帶向亞熱帶過(guò)度的氣候帶。蚌埠市汛期為5~9月,1952~2015年,汛期平均降水量為647.0mm,占年均降水的69.6%,最多的年份占全年降水的81.8%(1995年)。蚌埠市汛期降水主要受梅雨鋒擺動(dòng)、低渦、切變線或西北太平洋副熱帶高壓影響,多連陰雨或區(qū)域性大—暴雨過(guò)程。蚌埠汛期是暴雨洪澇、高溫干旱等災(zāi)害的高發(fā)時(shí)段,因此深入研究汛期降水的氣候演變規(guī)律,對(duì)做好汛期的短期氣候預(yù)測(cè)以及防災(zāi)減災(zāi)有著重要意義。目前,針對(duì)蚌埠汛期降水的長(zhǎng)期氣候演變規(guī)律的研究還較少。本文根據(jù)蚌埠市地面觀測(cè)站近64年實(shí)測(cè)資料,采用氣候傾向率、Morlet小波等方法從汛期降水線性變化特征、旱澇分布結(jié)構(gòu)、周期演變規(guī)律以及干濕期轉(zhuǎn)折突變等方面著手,探究蚌埠汛期降水多時(shí)間尺度的氣候演變規(guī)律。
本文以蚌埠市1952~2015年汛期月平均降水量為研究對(duì)象,資料序列長(zhǎng)度為64年,共320個(gè)月。首先分析蚌埠汛期降水年變化、月變化基本特征。在此基礎(chǔ)上,通過(guò)Morlet小波變換分析汛期降水的多時(shí)間尺度周期性分布,然后利用功率譜檢驗(yàn)找出最顯著的時(shí)頻變化周期。
采用最小二乘法求取線性回歸趨勢(shì)方程,其回歸方程表達(dá)式為:

式中,a為回歸系數(shù),b為回歸常數(shù)。
采用復(fù)值Morlet小波作為母小波函數(shù),其表達(dá)式為:

式中:c為常數(shù),i表示虛部。其小波變換式為:

式中:Wf(a,b)稱為小波變換系數(shù),f(t)為1個(gè)分析信號(hào)函數(shù),a為分辨尺度(也叫放大因子),b為平移因子為小波函數(shù)ψ(t)的復(fù)共軛函數(shù)。
將時(shí)間域上關(guān)于a的所有小波變換系數(shù)的平方進(jìn)行積分,即為小波方差。對(duì)于離散水文時(shí)間序列,小波方差可以采用下式計(jì)算:

功率譜分子中使用平滑功率譜,進(jìn)行漢寧平滑,公式為:

式中:St是功率譜密度估計(jì)值,r(τ)為時(shí)間序列的自相關(guān)系數(shù)序列,m是最大時(shí)間后延,n是時(shí)間序列長(zhǎng)度,周期Tl=2m/l,l為波數(shù)。根據(jù)相關(guān)文獻(xiàn),m一般選在n/3~n/10之間,文中取m=8。
圖1給出了蚌埠1952~2015年汛期降水距平圖。蚌埠汛期降水分布年代間分布不均。5年滑動(dòng)平均可以看出明顯的年際和年代際變化特征:20世紀(jì)90年代中期前,蚌埠汛期降水存在約4~6年的豐(水期)枯(水期)交替特征,即有8~12年左右的旱澇交替周期,分別表現(xiàn)為四個(gè)豐水期及四個(gè)枯水期(1959年以前、1963~1966年、1970~1975年、1989~1994年四個(gè)豐水期及1959~1963年、1966~1970年、1976~1983年、1983~1989年四個(gè)枯水期)。1993年后交替周期延長(zhǎng)為16~20年左右,表現(xiàn)為1993~2001年的枯水期及2002~2012年的豐水期。一元線性回歸分析顯示,整個(gè)序列呈上升趨勢(shì),但上升趨勢(shì)不顯著,降水傾斜率為1.4mm/a,每10年增加14mm。
文中定義降水量距平百分率≤-25%為旱,降水量距平百分率≥25%為澇。從圖1可看出,蚌埠干旱年份明顯多于洪澇,但洪澇強(qiáng)度明顯大于干旱。1968年之前,蚌埠旱澇交替頻繁,振幅也較大;20世紀(jì)60年代末期起,降水距平偏旱程度明顯,這種干旱程度一直持續(xù)至80年代末期;90年代起旱澇又再次頻繁交替,但幅度略低于60年代較低。近64年,共出現(xiàn)14個(gè)偏旱年,降水距平百分率小于-40%的有5年,其中1966年和1978年偏旱程度最顯著,降水距平百分率高達(dá)-65%和-62%;偏澇年工出現(xiàn)12年,超過(guò)40%的有4年,其中1956年、1991年、2007年偏澇程度顯著,降水距平百分率分別達(dá)89%、69%、69%。統(tǒng)計(jì)分析表明,蚌埠旱澇大多由西太平洋副高的位置和強(qiáng)度異常所致,如60年代末~80年代末偏旱程度明顯與當(dāng)年年汛期副高持續(xù)偏強(qiáng)有關(guān),西太平洋副高的位置較常年偏西、偏北,受副高控制,蚌埠出現(xiàn)持續(xù)的干旱天氣。

圖1 蚌埠市1952~2015年汛期降水距平圖
蚌埠市64年汛期各月降水量統(tǒng)計(jì)情況見表1。可見,蚌埠汛期降水各月分布不均勻,7月最多,8月次之,5月最少,7月降水量約達(dá)5月的3倍。月均降水日數(shù)也表現(xiàn)為7月、8月最多,9月、5月較少;總體而言,蚌埠汛期降水多集中在7~8月,7、8兩個(gè)月的降水量占汛期降水總量的55.6%。主要原因可能是因?yàn)?、8月份江淮切變線多維持,同時(shí)也多受到熱帶氣旋外圍影響,雨量及雨日均最多;而5月份雨日多于9月份雨量卻較少,則可能是由于前期梅雨期多連陰雨而雨量不大,而后期副高增強(qiáng)多對(duì)流及熱帶氣旋影響,雨日少但降水強(qiáng)度大。
對(duì)64年來(lái)汛期各月的降水量做一元線性回歸分析,得到5~9月各月降水量的傾斜率(表1)。由表1可知,5~9月降水傾斜率均為正,其中6月降水增加最明顯,7月份次之,但均未通過(guò)0.05顯著性檢驗(yàn),表明各月的增長(zhǎng)趨勢(shì)并不顯著。從增幅來(lái)看,64年來(lái),蚌埠市6月份和7月份降水量分別上升39.7mm和28.2mm,平均增幅22.3%,整個(gè)汛期上升87.7mm,上升幅度為13.6%。
多時(shí)間尺度小波分析的關(guān)鍵在于繪制小波系數(shù)圖和小波方差圖。小波系數(shù)圖可以得到關(guān)于時(shí)間序列變化的小波變化特征。在尺度一定的情況下,小波系數(shù)隨時(shí)間的變化過(guò)程反映了系統(tǒng)在該尺度下的變化特征:正的小波變換系數(shù)對(duì)應(yīng)于偏多期,負(fù)的小波變換系數(shù)對(duì)應(yīng)于偏少期,小波變換系數(shù)為零對(duì)應(yīng)著突變點(diǎn)。圖2為蚌埠市1952~2015年汛期月降水距平序列Morlet小波變換系數(shù)圖,從圖上可以看出不同時(shí)間尺度波動(dòng)情況。
對(duì)于1~10年年際尺度的周期變化來(lái)看,信號(hào)的主要變化特征體現(xiàn)在2年、5年及6~8年的周期上。首先,2年左右的周期貫穿始終,在80年代及近10年有所減弱。其次,5年左右的周期在70年代中期至2005年左右表現(xiàn)得比較明顯;第三,6~8年左右的周期在全時(shí)域都有明顯的響應(yīng),正負(fù)中心有規(guī)律地交替出現(xiàn)。1970年前,振蕩強(qiáng)度較強(qiáng),閉合中心值最大為1.5,等值線密集,70年代后,強(qiáng)度有所減弱。
對(duì)于10年以上的年代際尺度周期變化,首先存在16~18年周期,整個(gè)時(shí)域經(jīng)歷了四個(gè)完整周期,即1952~1960年、1968~1978年、1985~1993年、2002~2010年的降水較多期和1960~1968年、1978~1985年、1993~2002年的降水較少期,旱澇分界的轉(zhuǎn)折點(diǎn)出現(xiàn)在1960年、1968年、1978年、1985年、1993年和2002年。且振蕩強(qiáng)度逐漸增強(qiáng),1998年及2005年達(dá)到最強(qiáng)。汛期降水最多的時(shí)段出現(xiàn)在2002~2010年,對(duì)應(yīng)正的小波系數(shù)達(dá)1.5的閉合曲線,中心點(diǎn)在2005年,表明這一時(shí)段的汛期降水是近64年來(lái)最多的。對(duì)于32年尺度來(lái)說(shuō),蚌埠市汛期降水并沒有表現(xiàn)出明顯的周期性。
由小波系數(shù)的變化趨勢(shì)可以預(yù)測(cè)未來(lái)幾年氣候旱澇變化趨勢(shì)。根據(jù)年際尺度上強(qiáng)16年的周期振蕩規(guī)律,2015年后蚌埠市汛期降水演變趨勢(shì)仍處于負(fù)值中心附近,表明蚌埠市汛期降水仍呈減少的趨勢(shì)。
小波方差隨尺度變化的過(guò)程稱為小波方差變化圖,該圖能反映時(shí)間序列中所包含的各種尺度(周期)的波動(dòng)及其強(qiáng)弱(能量大小)隨尺度變化的特性。在一定尺度下,小波方差表示時(shí)間序列中該尺度周期波動(dòng)的強(qiáng)弱(能量大小)。因此,通過(guò)小波方差圖,可以確定序列存在的主要時(shí)間尺度,即主周期。
圖3為蚌埠市1952~2015年汛期月降水距平序列小波方差圖。從圖上可以看出8年尺度的變化周期較為突出,即蚌埠市汛期降水存在8年左右的主周期;峰值為18年,說(shuō)明蚌埠市汛期降水還存在18年左右的次周期;還存在著2年尺度周期。
通過(guò)對(duì)蚌埠市汛期降水量時(shí)間序列進(jìn)行功率譜分析,可以從譜曲線中的譜值最大來(lái)確認(rèn)主要振動(dòng)及其對(duì)應(yīng)的周期,進(jìn)而揭示和驗(yàn)證其變化規(guī)律。圖4為蚌埠市1952~2015年汛期降水量功率譜估計(jì)圖,對(duì)Morlet小波分析的結(jié)果進(jìn)行檢驗(yàn)。圖中實(shí)線為功率譜值,虛線為95%置信水平的紅色噪音譜。由圖4可見,在周期長(zhǎng)度為18年處,功率譜估計(jì)值為一峰值且大大超過(guò)標(biāo)準(zhǔn)譜,因此18年是第一顯著周期,并可確定12~18年的周期振蕩;另外在2.4~2.25年及5.14年處也超過(guò)了標(biāo)準(zhǔn)譜,所以還存在2年、5年的顯著周期。所以蚌埠市汛期降水序列,在年際尺度上2年、5年左右的周期通過(guò)了檢驗(yàn),在年代際尺度上存在明顯的18年周期振蕩。

圖2 蚌埠市1952~2015年汛期降水距平序列Morlet小波變換系數(shù)圖

表1 蚌埠市1952~2015年汛期各月降水量統(tǒng)計(jì)表

圖3 蚌埠市1952~2015年汛期降水距平序列小波方差圖

圖4 蚌埠市1953~2015年汛期降水量功率譜曲線圖
利用蚌埠市1952~2015年的月降水量資料,采用線性傾向估計(jì)和Morlet連續(xù)小波變換、功率譜分析等方法對(duì)近64年蚌埠市汛期(5~9月)降水多時(shí)間尺度演變特征進(jìn)行了分析。得出以下結(jié)論:
(1)近64年來(lái)蚌埠市市汛期降水主要呈波動(dòng)性變化;20世紀(jì)80年代前,汛期旱澇交替周期為8~12年,90年代后交替周期延長(zhǎng)為16~20年。
(2)汛期降水多集中在7~9月,7月降水量及降水日均最多。
(3)近64年來(lái),蚌埠汛期各月降水量及總降水量呈上升趨勢(shì),但上升趨勢(shì)并不顯著,均未通過(guò)0.05顯著性檢驗(yàn)。
(4)Morlet連續(xù)小波變換表明蚌埠市汛期降水存在6~8年左右的主周期,18年左右的次周期以及2年、5年周期,其中2年、5年、18年周期振蕩通過(guò)功率譜檢驗(yàn)
(作者單位:安徽省蚌埠市氣象局233040)