錢(qián)姝妮

(廣州市水務(wù)規(guī)劃勘測(cè)設(shè)計(jì)研究院有限公司,廣東 廣州 510641)

強(qiáng)降雨易造成洪水、泥石流、滑坡等自然災(zāi)害,往往造成嚴(yán)重的經(jīng)濟(jì)和生命財(cái)產(chǎn)損失。中國(guó)每年因強(qiáng)降雨導(dǎo)致的損失巨大,例如2012年7月發(fā)生在北京的特大暴雨,導(dǎo)致79人死亡,直接經(jīng)濟(jì)損失高達(dá)116億元[1];再如2018年登陸廣東省臺(tái)風(fēng)“山竹”,誘發(fā)了珠江流域下游連續(xù)多天的暴雨,大量農(nóng)田、城市交通、民房受損,直接經(jīng)濟(jì)損失逾30億元[2]。研究強(qiáng)降雨時(shí)空分布特征,對(duì)提高變化環(huán)境下強(qiáng)降雨演變規(guī)律認(rèn)識(shí)、增大強(qiáng)降雨防災(zāi)減災(zāi)能力具有重要意義。

Wood等[3]分析了歐洲日尺度和次日尺度極端降水變化,結(jié)果表明極端降水的變化程度隨著季節(jié)、歷時(shí)變化而異,次日尺度極端降水比日尺度極端降水增加更快。Ruiz-Alvarez等[4]研究了墨西哥極端日降水指數(shù)的變化趨勢(shì),發(fā)現(xiàn)日降水量超過(guò)50 mm的天數(shù)并沒(méi)有顯著增加,但年最大1天降水事件增加趨勢(shì)明顯。Wu等[5]分析了中國(guó)極端降水時(shí)程分配,發(fā)現(xiàn)雨峰靠后占主導(dǎo),且西北地區(qū)極端降水顯著增加。吳孝情等[6]解析了珠江流域降雨極值非平穩(wěn)性,結(jié)果表明珠江三角洲和東江流域南部、柳江流域東北部地區(qū)50/100年一遇的日降雨量級(jí)較大,而南盤(pán)江西部地區(qū)則較小。

強(qiáng)降雨往往由異常天氣系統(tǒng)所致,其發(fā)生前后通常伴隨著多種氣象變量(如風(fēng)速)的異常波動(dòng)或變化。有關(guān)研究表明,風(fēng)速與強(qiáng)降雨的強(qiáng)度、量級(jí)等特征存在一定的關(guān)聯(lián)性[7]。Zscheischler等[8]研究發(fā)現(xiàn),歐洲大部分地區(qū)極端降水與風(fēng)速呈正相關(guān)關(guān)系,特別當(dāng)極端日降水量超過(guò)100 mm時(shí),這種正相關(guān)關(guān)系尤為明顯。Um等[9]研究亦發(fā)現(xiàn),韓國(guó)南部沿海地區(qū)臺(tái)風(fēng)期間降雨量隨著風(fēng)速的增加而增大。Feng等[10]以中國(guó)臺(tái)灣西南地區(qū)為例,分析臺(tái)風(fēng)過(guò)程中上游風(fēng)速與地形強(qiáng)降雨的關(guān)系,發(fā)現(xiàn)低層風(fēng)速與小時(shí)尺度強(qiáng)降雨有很好的相關(guān)性,上游風(fēng)速?zèng)Q定了下游山區(qū)強(qiáng)降雨量。王柯等[11]分析了塔中雷暴降水天氣過(guò)程風(fēng)場(chǎng)變化,發(fā)現(xiàn)雷暴降水過(guò)程低層風(fēng)速增大而高層風(fēng)速減小。蔡奇霖等[12-13]以廣東省大隆洞水庫(kù)和東湖水庫(kù)為例,探討了強(qiáng)風(fēng)暴雨遭遇規(guī)律,結(jié)果顯示庫(kù)區(qū)年最大風(fēng)速和年最大暴雨存在一定的遭遇概率。

總結(jié)以往文獻(xiàn)發(fā)現(xiàn),強(qiáng)降雨具有高度時(shí)空異質(zhì)性,且由于強(qiáng)降雨定義不同、數(shù)據(jù)質(zhì)量差異等原因,同一地區(qū)的強(qiáng)降雨研究成果也有所差異。此外,在強(qiáng)降雨與風(fēng)速關(guān)系上,現(xiàn)有研究主要分析強(qiáng)降雨發(fā)生當(dāng)日的風(fēng)速變化(或兩者遭遇問(wèn)題),較少探討強(qiáng)降雨發(fā)生前后風(fēng)速的變化規(guī)律,對(duì)強(qiáng)降雨與前后期風(fēng)速的響應(yīng)關(guān)系認(rèn)識(shí)不足。基于此,以珠江流域?yàn)槔?定義95%百分位值為日強(qiáng)降雨閾值,運(yùn)用Mann-Kendall等方法分析1961—2018年期間強(qiáng)降雨時(shí)空變化特征,解析強(qiáng)降雨與前后期風(fēng)速響應(yīng)關(guān)系,以期為流域暴雨洪水防災(zāi)減災(zāi)提供參考依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

珠江流域面積約45萬(wàn)km2(圖1),地形主要以山地、丘陵為主,地勢(shì)總體上東低西高。流域氣候?qū)賮啛釒Ъ撅L(fēng)氣候,溫和多雨,年平均氣溫在14～22℃之間,年平均降雨量約1 200～2 200 mm,年內(nèi)雨量分配不均,且年際變化大。下游珠江三角洲地區(qū)城市密集,人口眾多,城鎮(zhèn)化程度高;由于面臨南海,易受臺(tái)風(fēng)、強(qiáng)對(duì)流天氣影響,強(qiáng)降雨事件頻發(fā),并引發(fā)區(qū)域甚至流域性大洪水。因此,開(kāi)展流域強(qiáng)降雨時(shí)空演變研究,對(duì)流域防洪減災(zāi)有重要意義。

圖1 珠江流域上游、中游、下游及183個(gè)氣象站點(diǎn)分布

1.2 基礎(chǔ)數(shù)據(jù)

采用分布于流域內(nèi)的183個(gè)氣象站日尺度降水、風(fēng)速數(shù)據(jù),時(shí)間跨度為1961—2018年,站點(diǎn)空間分布見(jiàn)圖1,來(lái)源于中國(guó)氣象局。將流域劃分為上游、中游和下游,其中位于上游區(qū)域的氣象站點(diǎn)41個(gè),中游94個(gè),下游48個(gè)。

1.3 研究方法

1.3.1日尺度強(qiáng)降雨事件定義

定義日降雨量超過(guò)0.1 mm為有雨天,將1961—2018年有雨天的日降雨量從小到大排序,并取95%百分位值為強(qiáng)降雨閾值,日降雨量超過(guò)閾值視為發(fā)生一次強(qiáng)降雨事件。采用頻率、降雨量2個(gè)基本指標(biāo)刻畫(huà)強(qiáng)降雨事件特征。

1.3.2binning方法

采用binning方法研究不同風(fēng)速條件下強(qiáng)降雨閾值變化情況[14],具體如下:基于風(fēng)速觀測(cè)記錄最小間隔(0.1 m/s),將風(fēng)速按0.1 m/s單位間隔劃分成多組(即0.1、0.2、0.3 m/s…),對(duì)每組風(fēng)速的日降雨量(如1 m/s風(fēng)速下的所有日降雨觀測(cè)記錄)從小到大排序,與強(qiáng)降雨定義相對(duì)應(yīng),取95%百分位值為當(dāng)前風(fēng)速條件下的強(qiáng)降雨閾值。

1.3.3Mann-Kendall趨勢(shì)分析方法

采用Mann-Kendall方法檢驗(yàn)強(qiáng)降雨和風(fēng)速的變化趨勢(shì)。Mann-Kendall為非參數(shù)統(tǒng)計(jì)方法,在時(shí)間序列趨勢(shì)診斷、突變分析方面應(yīng)用廣泛[15]。Mann-Kendall簡(jiǎn)要原理過(guò)程如下。

給定時(shí)間序列xi(i=1,2,…,n),計(jì)算統(tǒng)計(jì)量z和s:

(1)

(2)

當(dāng)統(tǒng)計(jì)量z大于零時(shí),序列呈增加趨勢(shì),小于零為減少趨勢(shì)。通過(guò)計(jì)算傾斜度β確定趨勢(shì)值:

(3)

2 結(jié)果與分析

2.1 日尺度強(qiáng)降雨與風(fēng)速時(shí)空變化特征

為分析強(qiáng)降雨事件及風(fēng)速年際變化及空間分布,統(tǒng)計(jì)每個(gè)站點(diǎn)、每一年發(fā)生強(qiáng)降雨事件的總降雨量、總頻率和風(fēng)速的年平均值,并總結(jié)上游、中游和下游區(qū)域內(nèi)的總體情況。表1為流域上游、中游和下游地區(qū)強(qiáng)降雨事件的總降雨量變化趨勢(shì)?？梢园l(fā)現(xiàn),強(qiáng)降雨總量最大的為下游地區(qū),最小為上游地區(qū),分別為532.72、321.72 mm,其分布規(guī)律與現(xiàn)有研究基本一致[16]。流域上游地區(qū)強(qiáng)降雨總量呈顯著下降趨勢(shì),變化率為-0.5 mm/a,中游和下游為不顯著增加趨勢(shì)。圖2展示了流域1961—2018年期間強(qiáng)降雨總量與頻率年際變化趨勢(shì)空間分布。由圖可知,年強(qiáng)降雨總量趨勢(shì)空間分布不均,空間上無(wú)明顯的分布規(guī)律;年頻率呈增加趨勢(shì)的站點(diǎn)主要分布在流域中游地區(qū),其中北部部分地區(qū)增加趨勢(shì)明顯。

表1 珠江流域上、中、下游1961—2018年面平均強(qiáng)降雨量及變化趨勢(shì)

a)降雨量

b)頻率圖2 年強(qiáng)降雨量與頻率變化趨勢(shì)空間分布

表2和圖3分別給出了流域上游、中游和下游地區(qū)的年平均風(fēng)速變化趨勢(shì)及空間分布情況。由表2可知,上游平均風(fēng)速最大,為2.7 m/s,中游風(fēng)速最小,為1.9 m/s。上游年平均風(fēng)速均呈上升趨勢(shì),且通過(guò)95%顯著性水平,變化率為0.01 m/s/a;與之相反,中游和下游平均風(fēng)速呈顯著下降趨勢(shì),變化率分別為-0.007 、-0.003 m/s/a。從圖3可知,平均風(fēng)速呈顯著上升趨勢(shì)的站點(diǎn)基本分布在上游地區(qū),而中游大部分地區(qū)站點(diǎn)的風(fēng)速呈現(xiàn)不顯著下降趨勢(shì)。

表2 珠江流域上、中、下游1961—2018年區(qū)域平均風(fēng)速及變化趨勢(shì)

圖3 珠江流域年平均風(fēng)速變化趨勢(shì)空間分布

2.2 不同風(fēng)速條件下日尺度強(qiáng)降雨閾值變化規(guī)律

為分析珠江流域上游、中游、下游地區(qū)不同風(fēng)速下強(qiáng)降雨閾值的變化,將區(qū)域內(nèi)(如上游)所有站點(diǎn)的日平均風(fēng)速、日降雨量數(shù)據(jù)利用binning方法進(jìn)行分類計(jì)算,結(jié)果見(jiàn)圖4。由圖可知,上游強(qiáng)降雨閾值隨著風(fēng)速的增加而減小,說(shuō)明在風(fēng)速較大條件下,強(qiáng)降雨量較小。中游和下游強(qiáng)降雨閾值與風(fēng)速的對(duì)應(yīng)關(guān)系在風(fēng)速5 m/s以下較為相似,強(qiáng)降雨閾值基本上隨著風(fēng)速遞增先增大后降低,在大約2 m/s的風(fēng)速條件下,強(qiáng)降雨閾值達(dá)到最大,在60～80 mm之間,這說(shuō)明中游和下游地區(qū)強(qiáng)降雨量并非隨風(fēng)速單調(diào)變化。值得注意的是,風(fēng)速在5 m/s以上時(shí),下游強(qiáng)降雨閾值出現(xiàn)個(gè)別較大值,當(dāng)風(fēng)速為5.4 m/s時(shí),強(qiáng)降雨閾值約90 mm,說(shuō)明個(gè)別情況下風(fēng)速大時(shí)強(qiáng)降雨量級(jí)也較大。

注:綠色、藍(lán)色和紅色曲線分別為上游、中游和下游強(qiáng)降雨閾值滑動(dòng)平均線。圖4 上游、中游和下游強(qiáng)降雨閾值隨風(fēng)速變化關(guān)系

2.3 日尺度強(qiáng)降雨事件前后期風(fēng)速變化特征

為研究珠江流域上游、中游、下游地區(qū)強(qiáng)降雨事件發(fā)生前后風(fēng)速的變化特征,對(duì)區(qū)域內(nèi)(如上游)1961—2018年發(fā)生的所有強(qiáng)降雨事件發(fā)生當(dāng)天、發(fā)生前后2 d的風(fēng)速進(jìn)行統(tǒng)計(jì)平均,結(jié)果見(jiàn)圖5。由圖可知,上游、中游、下游地區(qū)強(qiáng)降雨發(fā)生當(dāng)天平均風(fēng)速分別為1.97、1.76、2.16 m/s。上游地區(qū)發(fā)生強(qiáng)降雨前2 d風(fēng)速逐漸降低,下降平均速率約-0.065 m/s/d;強(qiáng)降雨發(fā)生后風(fēng)速繼續(xù)下降至約1.95 m/s,之后開(kāi)始回升,且風(fēng)速在強(qiáng)降雨發(fā)生后1 d出現(xiàn)顯著性變化(降低)。與之相反,中游和下游地區(qū)強(qiáng)降雨發(fā)生前風(fēng)速逐漸增大,發(fā)生前2 d至發(fā)生當(dāng)天的風(fēng)速平均增加速率分別為0.11、0.1 m/s/d,強(qiáng)降雨發(fā)生當(dāng)天風(fēng)速最大,且當(dāng)天風(fēng)速出現(xiàn)顯著性變化(增大);強(qiáng)降雨發(fā)生后風(fēng)速逐漸回落,至發(fā)生后2 d的風(fēng)速平均下降速率分別為-0.12、-0.18 m/s/d。此外,3個(gè)地區(qū)強(qiáng)降雨發(fā)生后1 d風(fēng)速均有不同程度的下降。綜合可知,珠江流域強(qiáng)降雨發(fā)生前后風(fēng)速均有明顯的變化,但不同地區(qū)風(fēng)速對(duì)強(qiáng)降雨的響應(yīng)(減少或增加)并不一致。

注:橫坐標(biāo)“-2”、“-1”、“1”、“2”分別表示強(qiáng)降雨發(fā)生前2 d、前1 d、后1 d和后2 d,帶“×”為風(fēng)速變化通過(guò)學(xué)生t檢驗(yàn)95%顯著性水平。圖5 上游、中游和下游強(qiáng)降雨事件發(fā)生當(dāng)天及前后2 d風(fēng)速變化

3 討論

百分位值法是目前國(guó)內(nèi)外定義極端降水事件的主要方法之一,其原理是將降水的逐年有雨天序列進(jìn)行頻率分析,一定程度上考慮了降水的空間差異性,但目前尚無(wú)流域統(tǒng)一的百分位選擇標(biāo)準(zhǔn),常用的有90%、95%和99%百分位等。研究表明,同一區(qū)域不同百分位值標(biāo)準(zhǔn)對(duì)應(yīng)的研究成果并不相同,甚至差異明顯。伍麗麗等[17]利用珠江流域51 a降水資料,通過(guò)比較90%、95%和97%百分位對(duì)應(yīng)的極端降水閾值空間分布,發(fā)現(xiàn)百分位值標(biāo)準(zhǔn)越高,流域東西部地區(qū)極端降水閾值差別越明顯,采用偏低(如90%)或偏高百分位均對(duì)結(jié)論合理性造成影響?；诖?本研究選用的95%百分位可認(rèn)為合理考慮了流域東西部降水差異及研究結(jié)果的合理性。

強(qiáng)降雨和風(fēng)速研究不可避免地受地形、數(shù)據(jù)質(zhì)量等因素的制約,相同地區(qū)的研究結(jié)果可能存在差異。為規(guī)避地形因素的影響,以往研究者采用網(wǎng)格化數(shù)據(jù)進(jìn)行研究,如Wu等基于網(wǎng)格化降水?dāng)?shù)據(jù),發(fā)現(xiàn)珠江流域1961—2013年上游極端降雨呈下降趨勢(shì),中下游為上升趨勢(shì),且中游增加趨勢(shì)較上游、下游明顯[18]。本研究中,上游強(qiáng)降雨量減少而中下游增加(表1),與該研究結(jié)論基本一致,可認(rèn)為本研究與基于網(wǎng)格化數(shù)據(jù)的研究結(jié)果差別不大。同樣在風(fēng)速變化研究上,有關(guān)研究表明珠江流域上游風(fēng)速表現(xiàn)為增強(qiáng)趨勢(shì)而中下游為減弱趨勢(shì)[19],這與本研究相關(guān)結(jié)論吻合(表2),說(shuō)明本研究中站點(diǎn)風(fēng)速與網(wǎng)格化風(fēng)速數(shù)據(jù)差異不明顯。另外,鑒于本研究采用站點(diǎn)較多(180多個(gè)測(cè)站),本研究中強(qiáng)降雨和風(fēng)速空間分布特征較以往有關(guān)基于站點(diǎn)數(shù)據(jù)的研究結(jié)果要復(fù)雜、規(guī)律性也較低。

在強(qiáng)降雨閾值與風(fēng)速關(guān)系分析中,上游強(qiáng)降雨閾值隨風(fēng)速增加而降低,中下游隨風(fēng)速先遞增后下降,但下游風(fēng)速較大時(shí)強(qiáng)降雨閾值出現(xiàn)極大值,特別當(dāng)風(fēng)速為5.4 m/s時(shí),強(qiáng)降雨閾值約90 mm,說(shuō)明風(fēng)速大時(shí)強(qiáng)降雨量級(jí)也較大,這可能與下游地處沿海、易發(fā)生臺(tái)風(fēng)、強(qiáng)對(duì)流等極端天氣有關(guān)[20]。除此特殊情況外,中下游強(qiáng)降雨閾值在2 m/s風(fēng)速條件下達(dá)到極大值。這些均說(shuō)明風(fēng)速越大,強(qiáng)降雨量級(jí)不一定越高,因此,現(xiàn)實(shí)中需要結(jié)合其他天氣條件綜合確定強(qiáng)降雨情況。

4 總結(jié)

以珠江流域?yàn)檠芯繀^(qū),通過(guò)Mann-Kendall等方法分析了流域內(nèi)183個(gè)氣象站1961—2018年日降雨、風(fēng)速資料,總結(jié)了日尺度強(qiáng)降雨與風(fēng)速時(shí)空變化特征以及強(qiáng)降雨發(fā)生前后風(fēng)速的響應(yīng)規(guī)律,得到如下結(jié)論。

a)流域上游年強(qiáng)降雨總量減少趨勢(shì)明顯,但中下游趨勢(shì)不明顯,上游年平均風(fēng)速呈顯著上升趨勢(shì),中游年平均風(fēng)速變化趨勢(shì)與上游相反,下游年平均風(fēng)速變化不明顯。

b)流域上游強(qiáng)降雨閾值隨著風(fēng)速的增加而減小,中下游強(qiáng)降雨閾值隨著風(fēng)速先遞增后下降。

c)珠江流域強(qiáng)降雨事件發(fā)生前后風(fēng)速出現(xiàn)明顯變化,上游強(qiáng)降雨發(fā)生前風(fēng)速逐漸下降,強(qiáng)降雨發(fā)生第2天開(kāi)始回升,中下游強(qiáng)降雨發(fā)生前風(fēng)速增加,強(qiáng)降雨發(fā)生后風(fēng)速逐漸下降。

本研究結(jié)果將豐富珠江流域極端降雨時(shí)空變化規(guī)律理論體系,提高強(qiáng)降雨與風(fēng)速的響應(yīng)關(guān)系認(rèn)識(shí),同時(shí)有助于強(qiáng)降雨早期癥候快速診斷、強(qiáng)降雨預(yù)報(bào)、預(yù)防規(guī)避強(qiáng)降雨導(dǎo)致的洪澇等災(zāi)害具有重要的參考意義。