徐永清,王承偉,劉春生,蔣慧亮,魏 磊

（1.黑龍江省氣候中心，黑龍江 哈爾濱 150030；2.黑龍江省氣象臺(tái)，黑龍江 哈爾濱 150030）

1 引言

中國(guó)暴雨洪澇災(zāi)害頻發(fā)，造成嚴(yán)重經(jīng)濟(jì)損失，對(duì)國(guó)家經(jīng)濟(jì)發(fā)展及建設(shè)帶來巨大挑戰(zhàn)。 觀測(cè)統(tǒng)計(jì)研究和氣候預(yù)估研究均表明， 全球變化背景下， 降水在時(shí)、空雙維度上變化均呈現(xiàn)顯著非均勻性，小時(shí)降水極值的增長(zhǎng)速度遠(yuǎn)超預(yù)期，氣候預(yù)估顯示，未來降水的時(shí)空分布非均勻性和極端性將更為突出， 致災(zāi)風(fēng)險(xiǎn)加?。?-3］。近年來，對(duì)短歷時(shí)強(qiáng)降水的時(shí)空分布特征研究逐漸增多［4-8］。 左璇等［9］選取鄭州氣象站不同短歷時(shí)強(qiáng)降水資料，采用頻率分析、小波分析和極值分析等方法，揭示了鄭州短歷時(shí)強(qiáng)降水的演變規(guī)律，并對(duì)各短歷時(shí)強(qiáng)降水進(jìn)行了分布擬合檢驗(yàn)和重現(xiàn)期估算。 郭渠等［10］利用重慶地區(qū)逐分鐘降水資料，對(duì)當(dāng)?shù)馗鞫虤v時(shí)強(qiáng)降水空間分布特征進(jìn)行了探討并分區(qū)。馬京津等［11］根據(jù)北京市強(qiáng)降雨的時(shí)空變化特征進(jìn)行分區(qū)，并進(jìn)行了不同歷時(shí)降雨量重現(xiàn)期的估算。 唐傳師等［12］通過分析江西省短歷時(shí)降水的時(shí)空分布特征，對(duì)當(dāng)?shù)囟虤v時(shí)降水與暴雨的關(guān)系進(jìn)行了探討。陸婷婷等［13］對(duì)北京2012 年“7.21” 和2016 年“7.20”特大暴雨過程開展觀測(cè)對(duì)比研究， 分析兩次過程的環(huán)流形勢(shì)、 降水分布與變化、 對(duì)流系統(tǒng)等方面的差異，并借助區(qū)域自動(dòng)站分鐘級(jí)降水觀測(cè)數(shù)據(jù)， 提取和對(duì)比分析兩次過程短歷時(shí)降水事件特征上的異同， 盡管“7·21”和“7·20”兩次特大暴雨過程降水總量相當(dāng)（均超過了400 mm），但災(zāi)害程度卻顯著不同，“7·21”過程歷時(shí)更短、雨勢(shì)更強(qiáng)，全市平均最大小時(shí)雨強(qiáng)明顯大于“7·20”過程，致災(zāi)性更強(qiáng)。

黑龍江省暴雨具有短歷時(shí)、突發(fā)性、局地性等特點(diǎn)［14-15］，暴雨落區(qū)與強(qiáng)度是預(yù)報(bào)中的重點(diǎn)和難點(diǎn)，本文使用黑龍江省國(guó)家氣象觀測(cè)站逐時(shí)降水資料，研究短歷時(shí)強(qiáng)降水的時(shí)空分布特征， 為防災(zāi)減災(zāi)和暴雨預(yù)報(bào)預(yù)警提供依據(jù)。

2 資料與方法

2.1 資料說明

使用國(guó)家氣象信息中心 “中國(guó)國(guó)家級(jí)地面氣象站逐小時(shí)降水?dāng)?shù)據(jù)集”中黑龍江省84 個(gè)國(guó)家級(jí)氣象觀測(cè)站(圖1)的數(shù)據(jù)，該數(shù)據(jù)統(tǒng)計(jì)了1991-2020 年共30 a 夏季1 h、3 h、6 h、12 h 所有≥0.1 mm 的降水及出現(xiàn)時(shí)間，均按升序排列，其中3 h、6 h、12 h 均為逐小時(shí)滾動(dòng)累計(jì)降水，時(shí)間精度為1 h。

圖1 黑龍江省84 個(gè)國(guó)家氣象觀測(cè)站分布

2.2 研究方法

2.2.1 百分位法

各站各時(shí)效降水量閾值按照百分位法進(jìn)行統(tǒng)計(jì)，具體如下：

按照升序排列的1 h、3 h、6 h、12 h 降水量序列X1，X2，…，Xm， …，Xn，某個(gè)值小于或等于Xm的概率

式中：m 為Xm的序號(hào)，n 為該時(shí)效降水序列中樣本個(gè)數(shù)。 如果序列中有300 個(gè)值（即n = 300），則第99 百分位上的值為排序后X297(P=98.8%)和X298(P=99.1%)的線性插值［16-17］。

由此得到各時(shí)效升序降水量序列中不同百分位等級(jí)（如90、95、99 等）對(duì)應(yīng)的降水量值。

2.2.2 空間插值使用反距離權(quán)重插值法對(duì)84 個(gè)氣象觀測(cè)站各短歷時(shí)（1 h、3 h、6 h、12 h）的極值、99%分位、95%分位和90%分位值在黑龍江省范圍內(nèi)進(jìn)行空間插值，按自然段點(diǎn)法將各歷時(shí)、分位降水量分為5 級(jí)，由低到高分別為低值區(qū)、較低值區(qū)、中間值區(qū)、較高值區(qū)和高值區(qū)。

3 短歷時(shí)降水空間分布特征

3.1 夏季平均降水空間分布特征

由圖2 可見， 黑龍江省夏季30 a 年平均降雨量呈現(xiàn)中部多，四周少的特點(diǎn)：中部地區(qū)的小興安嶺地區(qū)和張廣才嶺西部地區(qū)降雨量均在380 mm 以上，伊春、鶴崗、鐵力、尚志等站點(diǎn)400 mm 以上，為降雨中心［18］。 以此為中心雨量向周邊逐漸減少，其中北部和松嫩平原南部等個(gè)別站點(diǎn)雨量在270-300 mm，為黑龍江省夏季降水最少的地方， 其它地區(qū)基本在300-380 mm 之間。

圖2 黑龍江省30 a 年平均夏季降雨量空間分布

3.2 1 h 短歷時(shí)降水空間分布特征

由圖3 可見，黑龍江省1 h 降水極值空間差異較大，在34-88 mm 間，最大值和最小值間相差50 mm，60 mm 以上主要分布在松嫩平原西南部和北部、小興安嶺、 東南部山區(qū)和三江平原地區(qū)等地也有零散的分布；東南部山區(qū)和大興安嶺北部為低值區(qū)，均在44 mm 以下；較高值區(qū)域分布在高值區(qū)的外圍；較低值區(qū)域主要分布在大興安嶺和三江平原地區(qū)； 雨量分布呈現(xiàn)較明顯的西南高，向東向北逐漸降低。 這一特點(diǎn)在99%、95%、90%幾個(gè)分位上也表現(xiàn)得更加明顯和連續(xù)，但同時(shí)可以看到，在99%分位的小時(shí)雨量值區(qū)間在7.5-17.2 mm，95%、90%不足10 mm， 下降得非常迅速。

圖3 黑龍江省1 h（a)降水極值,(b)99%分位,(c)95%分位,(d)90%分位降水量分布圖

3.3 3 h 短歷時(shí)降水空間分布特征

黑龍江省3 h 降水極值在43-126 mm 之間 （圖4），松嫩平原大部、小興安嶺和三江平原東部較高，其它分布與1 h 接近，最低的區(qū)域在大興安嶺和張廣才嶺東北部、 三江平原西部， 基本在60 mm 以下。90%分位、95%分位、99%分位的降水高值區(qū)均在松嫩平原地區(qū)，向東向北逐漸減少，低值區(qū)域是大興安嶺和東部偏南地區(qū)。 與3 h 降水極值相比，各百分位區(qū)間降雨量同樣下降得非常迅速。

圖4 黑龍江省3 h（a)降水極值,(b)99%分位,(c)95%分位,(d)90%分位降水量分布圖

3.4 6 h 短歷時(shí)降水空間分布特征

由圖5 可見， 黑龍江省6 h 降水極值在52-170 mm，東南和西北低，其它地區(qū)相對(duì)較高，最高值區(qū)主要分布松嫩平原西部，在松嫩平原東部、小興安嶺等地呈不連續(xù)點(diǎn)狀分布，在100 mm 以上，較高值區(qū)分布在松嫩平原東部和小興安嶺南部， 降水在97-106 mm，低值區(qū)分布在大興安嶺北部和東部地區(qū)，降水在52-75 mm。 降水90%分位、95%分位、99%分位均是西部松嫩平原地區(qū)高， 北部大興安嶺地區(qū)和東南部山區(qū)低， 黑龍江省中部和東部地區(qū)降水量居于中間。

圖5 黑龍江省6 h（a)降水極值,(b)99%分位,(c)95%分位,(d)90%分位降水量分布圖

3.5 12 h 短歷時(shí)降水空間分布特征

12 h 降水極值在65-198 mm（圖6），高值區(qū)在松嫩平原的西部和小興安嶺東部， 降水在130-198 mm； 較高值區(qū)在松嫩平原南部和小興安嶺的南部；低值區(qū)在大興安嶺北部和東部大部， 降水量在90 mm 以下， 較低值區(qū)在東南山區(qū)西部和三江平原地區(qū)。降水90%分位、95%分位、99%分位均是西部松嫩平原地區(qū)高， 降水分別在14 mm、21 mm、41 mm 以上， 北部大興安嶺地區(qū)低， 降水分別在12 mm、15 mm、33 mm 以下,黑龍江省中部和東部地區(qū)降水量居于中間。

圖6 黑龍江省12 h（a)降水極值,(b)99%分位,(c)95%分位,(d)90%分位降水量分布圖

3.6 短歷時(shí)降水空間分布特征分析

由以上分析可見， 黑龍江省各短歷時(shí)降水極值及90%分位以上降水均呈現(xiàn)明顯的西南地區(qū)強(qiáng)、向東向北減弱的趨勢(shì)， 這與夏季年平均降水量分布并不一致，尤其是西南部地區(qū)，平均降雨量較小，但不同短歷時(shí)降雨的大值區(qū)均出現(xiàn)在這一區(qū)域， 說明該地區(qū)降雨較其它地區(qū)更具有短歷時(shí)的特征，時(shí)間短、突發(fā)性強(qiáng)。 從短歷時(shí)降水極值來看，在小興安嶺及張廣才嶺等地還有不連續(xù)的極值點(diǎn)出現(xiàn)， 說明在復(fù)雜地形下，也有出現(xiàn)更強(qiáng)短歷時(shí)降雨的可能。 另外，由以上分析還可以看到， 不同百分位的雨量下降得特別快，這也說明對(duì)于黑龍江而言，短歷時(shí)強(qiáng)降雨頻次相對(duì)較少。

4 短歷時(shí)降水時(shí)間分布特征

根據(jù)黑龍江省84 個(gè)氣象觀測(cè)站各短歷時(shí)降水年極值和年極值均值， 分析短歷時(shí)降水極端性的時(shí)間變化特征。 其中極值為當(dāng)年84 個(gè)站中該歷時(shí)降水最大值， 均值則為當(dāng)年各站該歷時(shí)最大值累計(jì)后的平均值。

4.1 1 h 短歷時(shí)降水時(shí)間變化特征

由圖7 可見， 黑龍江省1 h 降水極值和均值年際變化較大， 總體略有增加，極值增幅（0.5 mm/ a）＞平均值（0.03 mm/ a）。 從極值變化來看， 30 a 有6 a出現(xiàn)明顯峰值， 有3 次峰值出現(xiàn)在2011 年后， 共有4 a 降水極值＞80 mm，2011 年后3 a 達(dá)75%，而且是降水最大的3 a。 但均值在2001-2010 年間最大，30 a 間共有4 a 超過30 mm，均出現(xiàn)在這10 a 間。 因此，從1 h 降水來看，最大降水極值與均值并不同步，2011 年以后1 h 降水極端性有所增強(qiáng)。1 h 降水極值和均值的最小值均出現(xiàn)在2004 年， 分別為43.5 mm和27 mm。

圖7 黑龍江省1 h 降水極值和均值曲線圖

4.2 3 h 短歷時(shí)降水時(shí)間變化特征

黑龍江省3 h 降水極值和均值的年際變化趨勢(shì)基本一致（圖8），均值年際變化較小，極值年際變化加大,增加趨勢(shì)比較明顯（0.7 mm/ a）。3 h 降水極值和均值在2011-2020 年間均明顯高于其他年代， 極值最高出現(xiàn)在2013 年，達(dá)到126.7 mm，其次為2001 和2016 年，降水分別為121.9 mm 和119.3 mm。 3 h 降水極值和均值的最小值也出現(xiàn)在2004 年， 分別為60.2 mm 和43.9 mm。

圖8 黑龍江省3 h 降水極值和均值曲線圖

4.3 6 h 短歷時(shí)降水時(shí)間變化特征

由圖9 可見，黑龍江省6 h 降水極值和均值具有相同的年際變化趨勢(shì)， 峰值和谷值出現(xiàn)時(shí)間基本一致，均值與極值的年際波動(dòng)明顯較1 h 和3 h 小。 6 h降水均值年際變化較小， 極值年際增加趨勢(shì)比較明顯 （0.7 mm/ a)。 6 h 降水極值在2011-2020 年間最高，其中2020 年最高，達(dá)到170.3 mm，遠(yuǎn)大于次高年份的2010 年(132.2 mm）,均值在1991-2000 年最高。

圖9 黑龍江省6 h 降水極值和均值曲線圖

4.4 12 h 短歷時(shí)降水時(shí)間變化特征

黑龍江省12 h 降水極值和均值具有相同的年際變化趨勢(shì)（圖10），且與6 h 降水的變化趨勢(shì)基本相同，均值呈微弱地增加趨勢(shì)（0.06 mm/ a），極值增加相對(duì)明顯（0.5 mm/ a）。 12 h 降水極值最高值出現(xiàn)在2020 年， 達(dá)到198.7 mm， 遠(yuǎn)大于次高年份1998 年（167.5 mm），年際變化較大，最低值在85 mm 左右。各站12 h 年降水最大值的平均值， 在2011-2020 年最高，其它年代降水相差不大。

圖10 黑龍江省12 h 降水極值和均值曲線圖

4.5 分析

黑龍江省各短歷時(shí)降水的極值和均值年際變化趨勢(shì)基本一致，歷時(shí)越長(zhǎng)，一致性越高。 1 h、12 h 降水均值呈微弱的增加趨勢(shì)，3 h、6 h 均值略呈下降趨勢(shì)，各短歷時(shí)的極值均呈增加趨勢(shì)，相對(duì)來說歷時(shí)越短增加越明顯。 各短歷時(shí)降水均是在2011-2020 年間最高，極值的峰值均位于最近10 a，但最大降水均值并不同步，說明近10 a 的降水的極端性有所增強(qiáng)。另外，2020 年6 h、12 h 降水極值遠(yuǎn)超其它年份，分析認(rèn)為與臺(tái)風(fēng)“美莎克”影響有關(guān)，“美莎克”影響黑龍江省時(shí)，雖然1 h、3 h 降水未出現(xiàn)極值，但小時(shí)雨強(qiáng)較大的持續(xù)時(shí)間較長(zhǎng)，所以累計(jì)雨量大。

5 結(jié)論與討論

（1）黑龍江省短歷時(shí)降水分布趨勢(shì)與夏季降水量分布趨勢(shì)差異較大，基本上1 h、3 h、6 h、12 h 降水極值在西部松嫩平原地區(qū)最高，大、小興安嶺較低，東南山區(qū)最低，1 h、3 h、6 h、12 h 降水99%分位、95%和90%分位，在西部松嫩平原地區(qū)最高，東南山區(qū)較低，大、小興安嶺最低。 黑龍江省短歷時(shí)降水極值較高，但99%分位降水值迅速減小，較極值降低4倍左右，95%分位降水值較99%分位降低2 倍左右，90%分位降水值約占95%分位的60%， 黑龍江省短歷時(shí)降水高值在總降水樣本中出現(xiàn)比例較小。

（2）黑龍江省各短歷時(shí)降水的極值和均值年際變化趨勢(shì)基本一致， 1 h、12 h 降水均值呈微弱的增加趨勢(shì)，3 h、6 h 均值略呈下降趨勢(shì)， 各短歷時(shí)的極值均呈明顯增加趨勢(shì)，歷時(shí)越短增加越明顯。 各短歷時(shí)降水均是在2011-2020 年間最高， 近10 a 降水極值明顯增強(qiáng)。

考慮到資料的長(zhǎng)序列和可靠性， 本文沒有應(yīng)用區(qū)域自動(dòng)站觀測(cè)資料，僅用84 個(gè)國(guó)家氣象觀測(cè)站資料進(jìn)行分析。 另外，對(duì)于時(shí)間分布特征，僅考慮到降水最大值的極值與均值， 對(duì)于極值的時(shí)間分布特征未做詳細(xì)分析， 這些工作將在未來的研究中進(jìn)一步開展。