陳魯鳳 劉 杰 蔡斌斌

(山東科技大學(xué) 測繪與空間信息學(xué)院，山東 青島 266590)

0 引言

極端降水事件對人類生產(chǎn)生活有著相當(dāng)重要的影響,近年來國外科學(xué)家對全球各地極端降水的趨勢和極值等變化進(jìn)行了大量研究[1]。1998年,KARL等[2]指出20世紀(jì)相比于19世紀(jì)的美國年降水量增加了8%,而增加的降水量大部分歸因于暴雨量和暴雨日數(shù)的增加。ROY等[3]發(fā)現(xiàn)整個19世紀(jì)期間印度的極端降水事件發(fā)生的頻率處在增加狀態(tài)中,并且增加趨勢明顯。全球陸地各氣象站降水觀測數(shù)據(jù)顯示,在過去20世紀(jì)期間,中低緯度地區(qū)強降水事件發(fā)生的頻率普遍增加[4]。

中國極端降水事件的變化趨勢存在明顯的區(qū)域性差異[5]。楊金虎等[6]通過對全國各個氣象站的降水?dāng)?shù)據(jù)研究,發(fā)現(xiàn)中國的西北部、青藏高原、長江中下游地區(qū)和東南沿海地區(qū)的極端降水事件呈現(xiàn)明顯的增加的趨勢。閔屾等[7]通過分析中國20世紀(jì)的極端降水事件,研究發(fā)現(xiàn)了中國極端降水事件具有時間維度上的持續(xù)性和空間分布的區(qū)域性等特點。

本次研究在借鑒其他學(xué)者研究理論及方法的基礎(chǔ)上[8],運用山東省20個氣象站點的近44年的日降水?dāng)?shù)據(jù)來分析山東省暴雨的時空分布及變化特征。希望通過這一研究,能在一定基礎(chǔ)上揭示山東省暴雨的時空變化特征,對人們的生產(chǎn)生活提供一定幫助。

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

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

降水?dāng)?shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享網(wǎng),由于研究時間區(qū)間跨度較大,數(shù)據(jù)精度到日降水量且要有連續(xù)性,所以經(jīng)過篩選選取了山東省境內(nèi)具有代表性的 20個氣象站點,這20個站點除萊蕪和棗莊外均勻分布在山東省的每個地級市內(nèi),并且在沿海多降水地區(qū)如煙臺威海等地境內(nèi)平均有兩個氣象站點,這樣的選擇能對山東省暴雨有更全面的分析。確定站點后,獲取了各站點1970—2013年的月降水量和日降水量數(shù)據(jù)。

在識別降水量的基礎(chǔ)上確定暴雨事件,其關(guān)鍵是確定暴雨的閾值。本次研究采用絕對閾值的方法定義降水日和暴雨日,即根據(jù)國家降水等級劃分標(biāo)準(zhǔn),日降水量≥ 0.1 mm確定為一個降水日,日降水量≥ 50 mm為一個暴雨日[9]。

1.2 研究方法

1.2.1 反距離權(quán)重空間插值法

反距離權(quán)重法是一種常用而簡便的空間插值方法,它以插值點與樣本點間的距離為權(quán)重進(jìn)行加權(quán)平均,離插值點越近的樣本點賦予的權(quán)重越大[10]。相比于其他插值方法,反距離權(quán)重法(Inverse Distance Weight，IDW)插值的自然鄰近關(guān)系具有良好的自適應(yīng)分布特性,其結(jié)合自然鄰近關(guān)系,可有效地解決參考點分布不均勻問題,是一種自適應(yīng)的反距離權(quán)重的插值方法[11]。

1.2.2氣候傾向率

用xi表示樣本量為n的某一氣候變量,用ti表示xi所對應(yīng)的時間,建立xi與ti之間的一元線性回歸方程

xi=a+bti(i=1,2,…,n)

(1)

式(1)的含義是用一條合理的直線表示氣候變量x與其時間t之間的關(guān)系,其中a為常數(shù),b為回歸系數(shù)[12]。對觀測數(shù)據(jù)xi及相應(yīng)的時間ti,回歸系數(shù)b和常數(shù)a的最小二乘法估計為

(2)

(3)

其中,b× 10即氣候傾向率,單位為mm/10a?；貧w系數(shù)b的符號表示氣候變量x的趨勢傾向,b>0,說明隨時間t的增加x呈上升趨勢,b值越大,上升趨勢越明顯,b< 0,說明隨時間t的增加x呈下降趨勢,b值越小,下降趨勢越明顯[13]。

1.2.3小波分析

小波分析能清晰地揭示出隱藏在時間序列中的多種變化周期,充分反映系統(tǒng)在不同時間尺度中的變化趨勢,并能對系統(tǒng)未來發(fā)展趨勢進(jìn)行定性估計[14]。當(dāng)選擇好合適的小波函數(shù)后,通過小波變換獲得小波系數(shù),然后利用相關(guān)軟件繪制小波系數(shù)等值線圖[15]。

2 山東省降水量時空變化特征

2.1 降水量年際變化特征

通過對山東省近44年來的平均降水量趨勢分析(圖1)發(fā)現(xiàn)，1970—2013年山東省平均最低降水量出現(xiàn)在2002年,為416.0 mm;最高出現(xiàn)在2003年,為924.4 mm,相差達(dá)508.4 mm。在44年的時間里,降水量的總體趨勢是上升的,氣候傾向率達(dá)到0.91 mm/10a。

圖1 年降水量趨勢圖

對山東省近44年來的平均降水日數(shù)趨勢分析(圖2)發(fā)現(xiàn)。1970—2013年山東省最少降水日最少出現(xiàn)在1999年,為59.1 d;最多出現(xiàn)在2003年,為93.5 d,相差近34 d。在44年的時間里,降水日數(shù)的總體趨勢是減少的,降水日數(shù)傾向率為-1.73 d/10a,且減少較明顯。

圖2 年降水日數(shù)趨勢圖

2.2 年降水量空間分布特征

通過對各站點的年均降水量和年均降水日數(shù)空間特征分析,發(fā)現(xiàn)年均降水量和年均降水日數(shù)的空間分布特征具有很強的一致性。其受緯度因素、海陸因素以及地形因素的影響,由東南向西北逐漸減少,呈明顯的帶狀分布,并且有明顯的高低值中心。如圖3～4所示，審圖號為GS(2020)4814號。

圖3 年降水量空間分布圖

圖4 年降水日數(shù)空間分布圖

3 山東省暴雨時空變化特征

3.1 暴雨時間變化特征

通過對山東省近44年來的年均暴雨量趨勢分析(圖5)發(fā)現(xiàn)。1970—2013年山東省平均最低暴雨量出現(xiàn)在2002年,為95.87 mm;最高出現(xiàn)在2007年,為317.62 mm。在44年的時間里,降水量的總體趨勢是增加的,氣候傾向率達(dá)到5.34 mm/10a。

圖5 暴雨量趨勢圖

對山東省近44年來的年均暴雨日數(shù)趨勢分析(圖6)發(fā)現(xiàn)。1970—2013年山東省最少暴雨日最少出現(xiàn)在2002年,為0.67 d;最多出現(xiàn)在1990年,為4.07 d,相差近3.4 d。在44年里,暴雨日數(shù)傾向率為0.05 d/10a,暴雨日數(shù)的總體趨勢是增加的。

圖6 暴雨日數(shù)趨勢圖

對近44年來山東省不同年代年暴雨降水量和年降水量、年暴雨日數(shù)和年暴雨日數(shù)的貢獻(xiàn)率進(jìn)行分析(表1)。結(jié)果顯示，暴雨降水量對年降水量的貢獻(xiàn)很大,且暴雨降水量年代貢獻(xiàn)率之和一般穩(wěn)定在25%以上,表明暴雨降水量是山東省降水量的重要組成部分。

表1 暴雨貢獻(xiàn)率 單位：%

3.2 暴雨空間分布特征

分析山東省年均暴雨降水量空間分布如圖7～9所示,審圖號為GS(2020)4814號。年均暴雨量、暴雨日數(shù)和暴雨強度存在三個明顯的高值中心,分別在泰山山地地區(qū)、魯東南山地丘陵地區(qū)和山東半島上的威海地區(qū);而明顯的低值中心有兩個,分別是山東西部的聊城市和山東省北部的濱州市。其都受緯度海陸和地形因素影響:(1)受到緯度地帶性的影響,整體上呈現(xiàn)出自南向北遞減的趨勢;(2)受到海陸因素的影響,在整體上呈現(xiàn)出自沿海向內(nèi)陸遞減的規(guī)律,因此,在靠近海岸的處于山東半島的煙臺威海地區(qū)形成高值地區(qū);而山東省北部的德州和西部的聊城靠近內(nèi)陸,形成明顯的低值地區(qū);(3)受到受地形因素影響,表現(xiàn)在山地地區(qū)降水較多,這也是泰山山地地區(qū)和魯中魯南的山地丘陵地區(qū)形成明顯的高值中心的原因。

圖7 暴雨量分布圖

圖8 暴雨日數(shù)分布圖

圖9 暴雨強度分布圖

3.3 暴雨空間變化特征

在山東省年均暴雨降水量空間變化趨勢圖中(圖10),年均暴雨降水量變化率在空間分布上存在明顯差異。從圖10中可以看出，暴雨降水傾向率值存在3個明顯的高值中心,即威海市、淄博中部和菏澤西部,說明這些地方暴雨降水量的增加趨勢顯著;而暴雨降水量傾向率值存在3個明顯的低值中心,包括濰坊北部、臨沂南部棗莊和德州西北部,說明這些地區(qū)暴雨降水量減少趨勢較明顯;而在山東的其他地區(qū),暴雨降水量傾向率在0上下浮動不大,說明高地區(qū)的暴雨降水量變化趨勢不明顯。

圖10 暴雨量變化圖

在山東省年均暴雨日數(shù)空間變化趨勢圖中(圖11),年均暴雨日數(shù)變化趨勢在空間分布上差異明顯。從圖11中可以看出，暴雨降水傾向率值存在3個明顯的高值中心,即濟(jì)南中部地區(qū)、淄博中部地區(qū)和菏澤西部地區(qū),說明這些地方暴雨日數(shù)的變化趨勢處在增加狀態(tài)中,并且增加趨勢顯著;而暴雨降水量傾向率值存在1個明顯的低值中心,位于魯南部的臨沂棗莊地區(qū),說明此地區(qū)暴雨日數(shù)變化趨勢處在明顯地減少趨勢;而在山東的其他地區(qū),暴雨日數(shù)傾向率在0上下,說明這些地區(qū)的暴雨日數(shù)變化趨勢不明顯。

圖11 暴雨日數(shù)變化圖

在山東省年均暴雨強度空間變化趨勢圖中(圖12),年均暴雨強度變化趨勢在空間分布上差異明顯。從圖12中可以看出，暴雨強度傾向率數(shù)值存在1個明顯的高值中心,即濱州地區(qū),說明此地暴雨強度的變化呈現(xiàn)增加趨勢,并且增加趨勢顯著;而暴雨降水強度傾向率值存在1個明顯的低值中心,位于聊城西部,說明此地區(qū)暴雨強度變化呈現(xiàn)明顯的下降趨勢;而在山東的其他地區(qū),暴雨強度傾向率在0上下,說明這些地區(qū)的暴雨強度變化趨勢不明顯。

圖12 暴雨強度變化圖

暴雨空間變化特征圖10～12，審圖號為GS(2020)4814號。

4 結(jié)束語

本研究通過對山東省近44年來山東省暴雨降水的時間和空間的變化特征分析得出以下結(jié)論:

(1)山東省年降水量呈現(xiàn)整體緩慢增加趨勢,年降水日數(shù)呈現(xiàn)不明顯的下降趨勢。山東省的降水量、降水日數(shù)的空間分布特征受緯度和海陸因素的影響,呈現(xiàn)出由南向北遞減和由東部沿海向西北內(nèi)陸遞減的分布規(guī)律。

(2)山東省近44年來暴雨降水量和暴雨日數(shù)呈現(xiàn)明顯的上升趨勢,并且在時間變化上具有明顯的一致性,即暴雨降水量隨暴雨水日數(shù)的增多而增多,也表明山東省的暴雨強度在時間變化上也呈現(xiàn)上升趨勢。

(3)山東省的暴雨降水量、暴雨日數(shù)和暴雨強度的空間分布特征都受緯度、海陸和地形因素的影響,呈現(xiàn)出由南向北遞減和由東部沿海向西北內(nèi)陸遞減的分布規(guī)律,并且在泰山地區(qū)形成暴雨發(fā)生的高頻地帶。

(4)山東省暴雨降水量、暴雨日數(shù)和暴雨強度的空間變化趨勢各有不同,其中,暴雨量變化趨勢與暴雨強度變化趨勢在空間分布上具有明顯的一致性。

整體看來,山東暴雨特征主要與位置地形有很大關(guān)系:受海陸因素影響,沿海地區(qū)暴雨強度較大,而山東西北部地區(qū)暴雨強度較小;收地形因素影響,泰安山地地區(qū)一直是暴雨發(fā)生的高頻地帶。通過分析確定了暴雨事件高發(fā)地區(qū),可以通過采取相關(guān)措施,降低暴雨對人們生產(chǎn)生活帶來的不利影響,包括城市定時清理河道和下水道,防止暴雨時發(fā)生堵塞造成城市內(nèi)澇;居民增加暴雨危機(jī)意識,養(yǎng)成關(guān)注天氣提前防范的良好意識等。

本研究分析了1970—2013年中國山東省地區(qū)的暴雨氣候變化特征,為以后相關(guān)研究提依據(jù)。但強降水變化過程和機(jī)理,以及人類活動對強降水的影響等較為復(fù)雜,有待于在今后研究中進(jìn)行進(jìn)一步探討。