珠海市香洲城區(qū)降水變化特征及成因分析

曹永強(qiáng)，李玲慧，邵薇薇，李元菲

(1.遼寧師范大學(xué)地理科學(xué)學(xué)院，遼寧 大連 116029；2.中國水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室，北京 100038)

在全球變暖以及城市化進(jìn)程加快的雙重背景下，城市建成區(qū)面積迅速擴(kuò)張、建筑密度變大，人類活動加劇使得空氣中凝結(jié)核增多以及對流層空氣層結(jié)不穩(wěn)定性增強(qiáng)，致使產(chǎn)生熱島效應(yīng)和雨島效應(yīng)[1-2]，進(jìn)而影響城市降水特征發(fā)生變化。另外，地表不透水率增加過程伴隨的填洼、下滲量的減少和滯洪量的增加也使得防洪減災(zāi)問題更加突出[3]。關(guān)于降水特征的研究一直是我國氣候相關(guān)領(lǐng)域的研究熱點(diǎn)。如，陳隆勛等[4]對1951—1989年的中國降水的研究表明，全國年降水量呈減少趨勢，但西部降水量增長顯著；張建云等[5]研究發(fā)現(xiàn)盡管我國北方地區(qū)近年來降水量有所增加，但仍低于多年平均值；王澄海等[6]對我國西北地區(qū)降水周期變化特征進(jìn)行了研究，結(jié)果表明，西北地區(qū)降水呈現(xiàn)出非周期性振蕩和多平衡態(tài)的特點(diǎn)；林紓等[7]也對我國西北地區(qū)1960年以來的暴雨氣候變化特征進(jìn)行了分析，結(jié)果顯示西北地區(qū)暴雨發(fā)生最為頻繁是其東部地區(qū)，且以6—9月發(fā)生暴雨的概率最大。2000年以來，國際減災(zāi)戰(zhàn)略的實(shí)施表明，在自然災(zāi)害預(yù)防、防備和減災(zāi)3項(xiàng)工作中，災(zāi)害預(yù)防工作最為重要，且災(zāi)情信息獲取是災(zāi)害預(yù)防的重要工具[8-9]。暴雨是內(nèi)澇之源，暴雨時程分配特征對于暴雨徑流和城市防洪排澇有著重要影響[10-12]，以往的研究在中國各地區(qū)降水特征方面已取得了豐碩成果，但針對城區(qū)降水特征及暴雨時程分配的研究還相對較少。

由于特殊的地理位置和地形條件，珠海市香洲城區(qū)常遭受臺風(fēng)暴潮、外江洪水和當(dāng)?shù)亟涤陰淼暮闈碁?zāi)害。2017年第13號臺風(fēng)“天鴿”為珠海建市記載以來遭遇的最強(qiáng)臺風(fēng)，其引發(fā)的風(fēng)暴潮疊加天文大潮的波峰共同影響了潮位[13]，后續(xù)帶來的超強(qiáng)降雨、風(fēng)暴潮和巨浪嚴(yán)重威脅了沿海低洼地區(qū)人員的生命安全，全市受災(zāi)人員高達(dá)64.14萬人，直接經(jīng)濟(jì)損失高達(dá)204.5億元。鑒于此，本文分析氣候變暖背景下珠海市香洲城區(qū)不同等級降水的分布頻率、分布周期、變化趨勢以及暴雨的時程分配特征，以期為預(yù)防內(nèi)澇災(zāi)害發(fā)生及降低災(zāi)害損失等提供參考。

1 研究區(qū)概況

珠海市是我國首批4個經(jīng)濟(jì)特區(qū)之一，位于南海之濱、珠江口西岸(22°21′N～22°33′N，113°46′E～113°62′E)。香洲城區(qū)地勢相對低平，自西北向東南傾斜，是由板障山四周的平地構(gòu)成的中心城區(qū)(圖1)。氣候?qū)賮啛釒ШＱ笮约撅L(fēng)氣候，多年平均降水量約為2 065 mm，降水季節(jié)長，年降水日數(shù)為103天左右。近年來，香洲城區(qū)城鎮(zhèn)化發(fā)展迅速，2018年土地面積達(dá)559.29 km2，常住人口為109.81萬人，城鎮(zhèn)人口比為90.08%[14]。在旅游業(yè)等產(chǎn)業(yè)發(fā)展的影響下，城區(qū)內(nèi)土地利用類型變化顯著，加之氣候變化，使得該地區(qū)降水特征有所改變，洪澇災(zāi)害的發(fā)生越來越頻繁，給城市排水系統(tǒng)增添了極大壓力。

圖1 香洲城區(qū)地形

2 資料與方法

本文所使用的數(shù)據(jù)包括珠海市1981—2010年逐日降水?dāng)?shù)據(jù)，來源于珠海市規(guī)劃設(shè)計院，以及2011—2018年CMADS中國大氣同化驅(qū)動數(shù)據(jù)集V1.0及V1.1。由于研究區(qū)相關(guān)數(shù)據(jù)比較缺乏，可獲得的氣象數(shù)據(jù)年份受限，因此選用不同種類數(shù)據(jù)進(jìn)行補(bǔ)充研究。0.5 h間隔降水?dāng)?shù)據(jù)來源于NASA GPM IMERG，時間為2000—2018年，用于分析暴雨及以上等級降水的時程分配特征；下墊面數(shù)據(jù)中林地、水域、綠地等來源于Open Street Map影像；建筑物輪廓來源于91衛(wèi)圖助手；DEM數(shù)據(jù)來源于ALOS對地觀測衛(wèi)星，分辨率為10 m。

本文利用Manner-Kendall檢驗(yàn)法及小波分析法從多種角度分析香洲城區(qū)不同等級降水量變化趨勢及周期特征，具體方法參考文獻(xiàn)[15-17]。

3 結(jié)果與分析

3.1 不同等級降水分布頻率及變化趨勢

通過對收集的數(shù)據(jù)進(jìn)行計算可得到，1981—2018年香洲城區(qū)降水日數(shù)共5 487 d，平均每年141 d，其中，降水日數(shù)最多的是2016年，為191 d，最少的是2004年，為101 d[18-19]。圖2為香洲城區(qū)降水量年際變化及顯著性檢驗(yàn)，可見香洲城區(qū)降水量變化呈波動減小趨勢，下降率每10 a為65.70 mm，且下降趨勢通過了0.05的顯著性檢驗(yàn)。

(a) 降水量年際變化

根據(jù)GB/T 28592—2012《降水量等級標(biāo)準(zhǔn)》，以日雨量計可將降雨劃分為小雨、中雨、大雨、暴雨、大暴雨及特大暴雨6種不同等級，劃分標(biāo)準(zhǔn)為：24 h 降水總量在0.1～10 mm之間為小雨，10～25 mm 之間為中雨，25～50 mm之間為大雨，50～100 mm之間為暴雨，100～250 mm之間為大暴雨，250 mm以上為特大暴雨，本文將暴雨、大暴雨和特大暴雨統(tǒng)稱為總暴雨。根據(jù)統(tǒng)計，可得到不同等級降水日數(shù)頻率分布與降水量占比情況：香洲城區(qū)小雨只占總降水量的11.95%，但降水日數(shù)占總降水日數(shù)的69.55%；中雨占總降水量的17.21%，降水日數(shù)占總降水日數(shù)的14.53%；大雨占總降水量的22.98%，降水日數(shù)占總降水日數(shù)的9.02%；總暴雨量占總降水量的47.86%，降水日數(shù)占總降水日數(shù)的6.91%。結(jié)果表明，不同等級降水對總降水量的貢獻(xiàn)以大雨以上為主。春、秋季以暴雨為主，其降水量分別占該季總降水量的25.35%和28.70%；夏季以大雨為主，其降水量占該季總降水量的24.97%；冬季以中雨為主，其降水量占該季總降水量的33.03%。

小雨、中雨、大雨、暴雨、大暴雨和特大暴雨等不同等級降雨的變化率每10 a分別為-2.45 mm、0.52 mm、-16.62 mm、-46.37 mm、16.79 mm和28.23 mm。進(jìn)一步分析可知，1981—2018年，香洲城區(qū)大暴雨呈顯著增加趨勢(每10 a為16.79 mm)，中雨和特大暴雨每10 a增加幅度分別為0.52 mm和28.23 mm，但增加趨勢未經(jīng)過0.05顯著性檢驗(yàn)。小雨、大雨及暴雨呈減少趨勢，每10 a減少幅度分別為2.45 mm、16.62 mm和46.37 mm。從不同等級降水量占比情況來看，由大到小排序?yàn)椋罕┯辍⒋笥?、大暴雨、中雨、小雨、特大暴雨，且大雨和暴雨兩個等級降雨占比均超過20%，因此，1981—2018年，香洲城區(qū)降水量的減少主要是由于大雨及暴雨減少造成的。

3.2 不同等級降水的周期變化

3.2.1 降水日數(shù)的周期變化

圖3為不同等級降水日數(shù)的小波變換的實(shí)部等值線圖。分析香洲城區(qū)總降水和小雨量級的年降水日數(shù)的小波變換實(shí)部等值線圖(圖3(a)(b))可知，1981—2018年，在年降水日數(shù)的變化過程中均存在著1～7 a、7～18 a及18～38 a的3種尺度的周期變化規(guī)律。對應(yīng)7～18 a周期，降水日數(shù)經(jīng)歷了增多—減少—增多—減少—增多的轉(zhuǎn)換趨勢，正負(fù)交替出現(xiàn)4次振蕩；對應(yīng)18～38 a周期，降水日數(shù)經(jīng)歷了增多—減少—增多的變化趨勢。中雨的年降水日數(shù)變化過程(圖3(c))中存在著1～5 a、5～9 a、9～12 a、12～18 a、18～36 a和36～38 a多尺度的周期變化規(guī)律；大雨的年降水日數(shù)變化過程(圖3(d))中存在著1～5 a、5～11 a、11～19 a、19～36 a和36～38 a的多尺度的周期變化規(guī)律；總暴雨的年降水日數(shù)變化過程(圖3(e))中存在著1～7 a、7～18 a、18～31 a和31～38 a的周期變化規(guī)律。總體來說，香洲城區(qū)不同等級的降水日數(shù)呈現(xiàn)多時間尺度的復(fù)雜嵌套，小尺度嵌套在較大尺度的周期性變化特征。

(a) 總降水

(c) 中雨

(d) 大雨

(e) 總暴雨

圖4為不同等級降水日數(shù)的小波方差圖，可以看出，總降水日數(shù)及小雨日數(shù)的小波方差圖中存在3個較為明顯的峰值，降水日數(shù)序列的能量中心分別位于5 a、14 a和25 a的特征時間尺度，說明在這4年左右的振蕩最為強(qiáng)烈，其中，最大峰值對應(yīng)25 a的時間尺度，為降水日數(shù)變化的第一主周期，14 a和 5 a的時間尺度對應(yīng)著第二、第三峰值，分別為降水日數(shù)變化的第二、第三主周期。中雨存在5個較為明顯的峰值，分別為4 a、7 a、10 a、14 a和25 a，其中，降水日數(shù)變化的第一主周期同樣為25 a，14 a、4 a、10 a 和 7 a 分別為中雨日數(shù)變化的第二至第五主周期。大雨日數(shù)存在4個較為明顯的峰值，能量中心分別位于 14 a、25 a、10 a和2 a特征時間尺度，其中，14 a左右的周期震蕩為大雨日數(shù)變化的第一主周期，25 a、10 a 和2 a分別為大雨日數(shù)變化的第二至第四主周期?？偙┯甏嬖?個較為明顯的峰值，分別為4 a、14 a和24 a，其中，14 a左右的周期震蕩為總暴雨日數(shù)變化的第一主周期，24 a和4 a分別為總暴雨日數(shù)變化的第二、第三主周期。此外，從主振蕩周期上看，香洲城區(qū)總暴雨日數(shù)減少的等值線還未閉合，說明在未來一段時間內(nèi)該地區(qū)總暴雨日數(shù)將處于偏少的時期。

圖4 不同等級降水日數(shù)小波方差

3.2.2 降水量的周期變化

圖5為不同等級降水量小波變換的實(shí)部等值線圖。由圖5(a)可見，1981—2018年，年總降水量的變化過程中均存在著1～2 a、2～7 a、7～18 a、18～33 a 和33～38 a的5種尺度的周期變化規(guī)律。對應(yīng)7～18 a周期，降水量經(jīng)歷了增多—減少—增多—減少—增多的正負(fù)交替4次豐枯振蕩；對應(yīng)18～33 a周期，降水日數(shù)經(jīng)歷了增多—減少—增多2次振蕩變化。其中，小雨(圖5(b))存在1～9 a、9～18 a和18～38 a的周期變化規(guī)律，中雨(圖5(c))存在1～5 a、5～9 a、9～18 a、18～35 a和35～38 a的多尺度周期變化規(guī)律，大雨(圖5(d))存在1～6 a、6～11 a、11～19 a、19～36 a和36～38 a 的變化規(guī)律，總暴雨(圖5(e))存在1～2 a、2～9 a、9～20 a、20～31 a和31～38 a的周期變化規(guī)律?？傮w來說，與降水日數(shù)變化周期相似，香洲城區(qū)不同降水量級的年降水量同樣呈現(xiàn)出多時間尺度的復(fù)雜嵌套，小尺度嵌套在較大尺度的周期性變化特征，且不同等級降水的日數(shù)與降水量的振蕩周期大致相似。

(a) 總降水

圖6為不同等級降水量小波方差圖，可以看出，總降水量及小雨量級的小波方差均存在3個較為明顯的峰值，在5 a、14 a和25 a的時間尺度振蕩最為強(qiáng)烈。其中，最大峰值對應(yīng)著25 a的時間尺度，為降水量的變化的第一主周期，14 a和5 a時間尺度分別為第二、第三主周期。中雨存在4個較為明顯的峰值，為4 a、7 a、14 a和25 a，其中，最大峰值同樣對應(yīng)25 a的時間尺度，為降水量變化的第一主周期，14 a、7 a、4 a分別為中雨降水量變化的第二至第五主周期。大雨同樣存在4個較為明顯的峰值，為2 a、10 a、14 a和25 a，其中，25 a左右的周期震蕩為大雨日數(shù)變化的第一主周期，14 a、10 a和2 a分別為大雨日數(shù)變化的第二至第四主周期。總暴雨存在3個較為明顯的峰值，分別為5 a、15 a和25 a，其中，15 a左右的周期震蕩為總暴雨量變化的第一主周期，5 a和 25 a分別為總暴雨量變化的第二、第三主周期。此外，從主振蕩周期上看，香洲城區(qū)總降水量和總暴雨量減少的等值線還未閉合，說明在未來一段時間內(nèi)該地區(qū)總降水量及總暴雨量將處于偏少的時期。

圖6 不同等級降水量小波方差

3.3 總暴雨場次年際及年內(nèi)變化特征

1981—2018年，珠海市香洲城區(qū)降水達(dá)暴雨以上等級共計412場，總暴雨出現(xiàn)頻次總體呈波動上升趨勢，速率每10 a為0.36次，2013年總暴雨頻次達(dá)到峰值，2017年總暴雨頻次最少。以旬為單位統(tǒng)計暴雨在年內(nèi)發(fā)生的次數(shù)用于識別場次暴雨年內(nèi)分布特征，結(jié)果表明，研究期內(nèi)珠海市香洲城區(qū)暴雨最早出現(xiàn)于1月中旬，最遲出現(xiàn)于12月上旬，集中發(fā)生于4—10月，暴雨發(fā)生頻次合計363場次，占場次暴雨總數(shù)的88.11%，其中，在4—6月前汛期雨季和7—10月后汛期雨季暴雨發(fā)生頻次分別占場次暴雨總數(shù)的43.69%和44.42%。

3.4 場次大暴雨時程分配及成因分析

通過珠海市2000—2018年逐日降水?dāng)?shù)據(jù)篩選出降水等級達(dá)大暴雨及以上等級的日期，以日為單位，對應(yīng)NASA GPM IMERG時間間隔0.5 h的降水?dāng)?shù)據(jù)，利用Origin 8.0繪制統(tǒng)計出的研究期間51場大暴雨及特大暴雨的時程分配特征，結(jié)果見圖7。由圖7可見，場次暴雨持續(xù)時間在4.5～23 h 之間，且降雨時長有增加趨勢。為進(jìn)一步探究其原因，對香洲城區(qū)城市化發(fā)展帶來的土地利用類型變化情況進(jìn)行統(tǒng)計分析，結(jié)果如圖8所示。可見，1981—2018年林地面積占比由39.29%上升至41.41%，這是由于珠海市設(shè)立經(jīng)濟(jì)特區(qū)伊始，大力發(fā)展旅游業(yè)，將保護(hù)生態(tài)環(huán)境作為城市發(fā)展的戰(zhàn)略目標(biāo)[20]。農(nóng)田面積占比由23.34%下降到4%，農(nóng)村居民點(diǎn)從有到無，大量農(nóng)業(yè)人口轉(zhuǎn)變?yōu)槌擎?zhèn)人口[21]，產(chǎn)業(yè)和經(jīng)濟(jì)快速發(fā)展，使得城市建設(shè)用地面積占比由35.11%上升到51.52%，從而導(dǎo)致香洲城區(qū)產(chǎn)生“城市熱島效應(yīng)”，使城區(qū)降雨時長增長[22-23]。此外，51場暴雨中有39場峰值區(qū)間出現(xiàn)在夜間，8場出現(xiàn)在上午，6場出現(xiàn)在午間，5場出現(xiàn)在午后。究其原因，與空氣中氣溶膠粒子的日變化規(guī)律有關(guān)。據(jù)統(tǒng)計，近年來珠海市空氣中PM2.5和PM10質(zhì)量濃度分別為28.12 μg/m3和45.94 μg/m3，SO2和NO2質(zhì)量濃度分別為7.62 μg/m3和30.21 μg/m3，CO和O3質(zhì)量濃度分別為 0.77 μg/m3和58.83 μug/m3。且空氣中PM2.5質(zhì)量濃度在5：00—7：00達(dá)到高值，顆粒物濃度隨太陽輻射的增強(qiáng)和近地面大氣混合層高度的增加而下降。由于東南沿海的地理位置，在午后多出現(xiàn)局部對流天氣，PM2.5質(zhì)量濃度降低，隨著晚高峰的出現(xiàn)，氣溶膠濃度有所回升，因此，暴雨峰值多出現(xiàn)在夜間。

(a) 2000-06-18至2008-04-19

(a) 1981年

4 結(jié) 論

a.1981—2018年，珠海市香洲城區(qū)總降水量、小雨、大雨及暴雨等級降水量呈波動減小趨勢，下降率每10 a分別為65.70 mm、2.45 mm、16.62 mm和 46.37 mm；中雨、大暴雨和特大暴雨呈增加趨勢，每 10 a 增加幅度分別為0.52 mm、16.79 mm和28.23 mm。

b.總降水量的減少主要是由于大雨及暴雨減少所造成。年降水日數(shù)及總降水量的周期變化均經(jīng)歷了增多—減少—增多—減少—增多的4次正負(fù)振蕩變換，并以25 a的周期振蕩為第一周期。

c.場次暴雨持續(xù)時間在4.5～23 h之間，且降水時長有增加趨勢，主要原因?yàn)槌鞘谢M(jìn)程引起地面粗糙度增大對移動的降雨系統(tǒng)產(chǎn)生了阻礙效應(yīng)，使城區(qū)雨時增長；暴雨峰值的出現(xiàn)時間頻率由大到小排序?yàn)橐归g、上午、午間、午后，這與空氣中氣溶膠粒子的日變化規(guī)律相關(guān)。