尤立, 齊天杰,2, 胡春明*

(1.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心工業(yè)廢水無(wú)害化與資源化國(guó)家工程研究中心, 北京 100085; 2.河北先進(jìn)環(huán)保產(chǎn)業(yè)創(chuàng)新中心有限公司, 石家莊 050035)

“一帶一路”是中國(guó)提出的新時(shí)期國(guó)際和地區(qū)合作倡議,旨在形成一個(gè)完整的歐、亞、非經(jīng)貿(mào)合作共同體,進(jìn)而密切沿線各國(guó)間的國(guó)際關(guān)系[1]。“一帶一路”合作倡議實(shí)施面臨的重要問(wèn)題之一是水資源問(wèn)題[2-3],特別是沿線國(guó)家基本上以發(fā)展中國(guó)家為主,受人類活動(dòng)和氣候變化的雙重影響,部分區(qū)域極端水文氣象事件頻發(fā)、水資源不均衡問(wèn)題突出,給“一帶一路”合作倡議的整體布局和實(shí)施帶來(lái)了不確定性影響[4]。據(jù)統(tǒng)計(jì),1995—2015年,對(duì)“一帶一路”合作倡議沿線國(guó)家而言,極端降水是危害最大的自然災(zāi)害致災(zāi)因子[5]。

斯里蘭卡位于印度洋東西交通的“十字路口”,是中國(guó)建設(shè)“21世紀(jì)海上絲綢之路”的重要伙伴國(guó)[6]。受熱帶季風(fēng)氣候影響,全年高溫,旱季和雨季分明。島內(nèi)水資源時(shí)空分布不均衡,降水期相對(duì)集中,地表水體面積變化波動(dòng)性較大,當(dāng)?shù)貥O端降水受地形特征的影響十分復(fù)雜,暴雨洪澇頻繁[7]。如2021年6月和11月,斯里蘭卡兩次遭遇極端降水災(zāi)害,引發(fā)洪水和山體滑坡,超過(guò)65萬(wàn)人日常生活受到影響。由于極端降水的突發(fā)性較強(qiáng),對(duì)區(qū)域社會(huì)經(jīng)濟(jì)發(fā)展的威脅較大,使之成為當(dāng)?shù)刎叫杞鉀Q的民生問(wèn)題[8]。

識(shí)別“絲綢之路經(jīng)濟(jì)帶”氣候變化特征是“一帶一路”合作倡議與建設(shè)的迫切需求[9]。受限于長(zhǎng)序列的降水資料缺失和觀測(cè)站的空間代表性不足,“一帶一路”相關(guān)地區(qū)降水時(shí)空變化的研究較少。Hurst 指數(shù)提供了對(duì)時(shí)間序列數(shù)據(jù)特征記憶的長(zhǎng)期存儲(chǔ)真實(shí)度量方法,在所有領(lǐng)域的時(shí)間序列分析中具有廣泛性和普適性,同樣也適用于較短時(shí)間序列的數(shù)據(jù)分析。基于斯里蘭卡2008—2020年降水?dāng)?shù)據(jù),通過(guò)經(jīng)驗(yàn)正交函數(shù)分析、RClimDex極端降水模型和Morlet小波分析法對(duì)斯里蘭卡極端降水演變特征開(kāi)展系統(tǒng)研究和評(píng)估,在揭示區(qū)域極端降水時(shí)空變化規(guī)律的基礎(chǔ)上,采用Hurst指數(shù)預(yù)測(cè)未來(lái)趨勢(shì),以期為當(dāng)?shù)胤罏?zāi)減災(zāi)和水資源管理提供重要借鑒,也可為研究“一帶一路”沿線的東南亞和南亞國(guó)家極端降水演變規(guī)律提供參考。

1 研究區(qū)概況及數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

斯里蘭卡位于印度洋中北部,連接著馬六甲海峽與波斯灣,是東西方海運(yùn)的十字路口,被譽(yù)為“印度洋上的明珠”。地理位置在5°55′N～9°50′N、79°42′E～81°53′E,全島總面積為65 610 km2,陸地面積約62 700 km2,水體面積(包括河流、湖泊、水庫(kù)和坑塘)約為2 910 km2[10]。海島地勢(shì)差異明顯,北部沿海平原寬闊,中南部為地形崎嶇的高原內(nèi)陸。氣候類型為熱帶季風(fēng)氣候,年內(nèi)只分旱季和雨季,降水具有明顯的季節(jié)特征。1998—2020年斯里蘭卡多年平均降水量為1 759.37 mm,當(dāng)?shù)亟邓磕觌H變化幅度較大?？臻g分布上主要表現(xiàn)為西南多東北少的特點(diǎn),西南部年降水量2 000～3 000 mm,東北部年降水量約1 000 mm。

1.2 數(shù)據(jù)來(lái)源

所使用到的斯里蘭卡氣象站點(diǎn)實(shí)測(cè)數(shù)據(jù)來(lái)源于2008—2020年的數(shù)據(jù)統(tǒng)計(jì),主要包括了斯里蘭卡海島的10個(gè)氣象站的逐日降水?dāng)?shù)據(jù),為降低因海島地勢(shì)差和站點(diǎn)分布不均問(wèn)題對(duì)雨量計(jì)算的影響,采用泰森多邊形法來(lái)計(jì)算海島的年均降水量[11],各站點(diǎn)的空間分布如圖1所示。

圖1 斯里蘭卡海島氣象站空間分布Fig.1 Spatial distribution of weather stations in Sri Lanka

為確保氣候變化研究的可靠性,在進(jìn)行極端降水計(jì)算之前,采用RClimDex軟件對(duì)斯里蘭卡原始?xì)庀髷?shù)據(jù)進(jìn)行嚴(yán)格的質(zhì)量控制及均一化檢驗(yàn),并結(jié)合廣泛推薦應(yīng)用的非參數(shù)統(tǒng)計(jì)方法M-K(Mann-Kendall)突變檢驗(yàn)法來(lái)揭示整體時(shí)間序列的變化趨勢(shì),從而判斷是否存在氣候突變[12-13]。圖2為給定顯著性水平P=0.05下的M-K檢驗(yàn)結(jié)果,其中P=0.05時(shí),統(tǒng)計(jì)量臨界值U0.05=±1.96。結(jié)果表明:在2008—2020年,降水量正序列(UF)和反序列(UB)均處于臨界值線內(nèi)(U0.05=±1.96),降水序列未發(fā)生明顯突變。兩曲線分別在2009—2010年和2015—2016年存在交點(diǎn),2008年后UF統(tǒng)計(jì)量在0線以下,表明斯里蘭卡本階段年降水呈下降趨勢(shì),從2010年之后,UF統(tǒng)計(jì)量均處于0線以上,斯里蘭卡降水趨勢(shì)呈穩(wěn)步上升趨勢(shì),其突變并未有顯著趨勢(shì)。

圖2 降水量M-K統(tǒng)計(jì)檢驗(yàn)Fig.2 M-K statistical test of precipitation

2 資料與方法

2.1 經(jīng)驗(yàn)正交函數(shù)分析

經(jīng)驗(yàn)正交函數(shù)(empirical orthogonal function,EOF)分析,因其極強(qiáng)的實(shí)用意義成為統(tǒng)計(jì)氣象要素變化的重要方法。原理是將要素場(chǎng)分解得到兩個(gè)不同的函數(shù):空間函數(shù)和時(shí)間函數(shù),將體現(xiàn)要素強(qiáng)烈顯著性的數(shù)據(jù)集中于幾個(gè)空間分布模態(tài)中,用極少的兩個(gè)或幾個(gè)基本涵蓋要素場(chǎng)信息的模態(tài)來(lái)客觀描述其空間結(jié)構(gòu)和時(shí)間變化趨勢(shì)[14]。

將氣象降水在m個(gè)格點(diǎn)上的n次連續(xù)監(jiān)測(cè)資料看作一個(gè)m維隨機(jī)向量X,經(jīng)驗(yàn)正交函數(shù)EOF分解是將降水場(chǎng)X分解為特征向量矩陣S和對(duì)應(yīng)的時(shí)間系數(shù)矩陣T,可表示為

(1)

降水?dāng)?shù)據(jù)經(jīng)預(yù)處理后,得到距平矩陣X,將X矩陣與其轉(zhuǎn)置矩陣XT乘積得到協(xié)方差矩陣C,將協(xié)方差矩陣C的特征向量λ由大到小順次排列,分解所對(duì)應(yīng)的特征向量,從而得到降水場(chǎng)的第1模態(tài),第2模態(tài),…,第m模態(tài)。對(duì)于任一模態(tài),空間特征向量分量值越大,則降水變化的振蕩強(qiáng)度越顯著;時(shí)間系數(shù)的正負(fù)代表與空間變化的異同性,時(shí)間系數(shù)為正值,表明與當(dāng)前空間分布同向,為變化主趨勢(shì);時(shí)間系數(shù)為負(fù)值,表明當(dāng)前空間分布狀態(tài)與降水變化趨勢(shì)相反。

2.2 RClimDex極端降水模型

基于R語(yǔ)言的RClimDex模型,經(jīng)模型本身的數(shù)據(jù)質(zhì)量控制模塊對(duì)逐日降水和氣溫?cái)?shù)據(jù)進(jìn)行嚴(yán)格篩查后,可對(duì)不同種類的極端氣候指數(shù)進(jìn)行分析討論[15-16]?？紤]到斯里蘭卡熱帶氣候下不存在小于0 ℃的溫度條件,對(duì)極端氣溫指數(shù)研究意義不大,故選取6項(xiàng)極端降水指標(biāo)來(lái)探究斯里蘭卡極端降水的時(shí)空變化特征。各極端降水指數(shù)定義如表1所示。

表1 極端降水指數(shù)定義Table 1 Definition of extreme precipitation index

2.3 Morlet小波分析

Morlet小波分析是將某一時(shí)空域信號(hào)經(jīng)平移、尺度變換后分解為系列小波函數(shù)的疊狀態(tài)的分析方法,能直觀反映變量在整個(gè)時(shí)段內(nèi)的周期震蕩情況,同時(shí)定性分析預(yù)測(cè)系統(tǒng)的未來(lái)變化趨勢(shì),常用于水文氣象的周期分析中[17-18]。定義小波函數(shù)為

(2)

2.4 R/S分析

R/S分析法常被用于探索時(shí)間的分形特征及長(zhǎng)期記憶過(guò)程中,依據(jù)計(jì)算所得Hurst指數(shù)H來(lái)判斷在時(shí)間序列上的可持續(xù)性[19]。Hurst指數(shù)取值在[0,1],當(dāng)H=0.5時(shí),表明該降水指數(shù)的變化具隨機(jī)性,現(xiàn)有變化與未來(lái)趨勢(shì)無(wú)相關(guān);0≤H<0.5時(shí),過(guò)去與未來(lái)的變化趨勢(shì)呈反持續(xù)性,H越接近0,反持續(xù)性越高;0.5

3 結(jié)果與討論

3.1 年際降水量EOF時(shí)空分解

3.1.1 空間分布特征

應(yīng)用EOF法對(duì)斯里蘭卡2008—2020年降水場(chǎng)進(jìn)行時(shí)空分解,按照North顯著性檢驗(yàn)方法,只有當(dāng)各模態(tài)下的特征根誤差上下限互相無(wú)重疊時(shí),才能判定此模態(tài)通過(guò)了North顯著性檢驗(yàn)。由表2可知,前6個(gè)模態(tài)的累計(jì)方差貢獻(xiàn)率達(dá)98.60%。其中第1、第2模態(tài)下特征根誤差上下限無(wú)重疊,因此通過(guò)North檢驗(yàn),其累計(jì)方差貢獻(xiàn)率達(dá)到90.13%,分布特征能夠很好地表征斯里蘭卡降水場(chǎng)的空間結(jié)構(gòu)特征。

表2 斯里蘭卡年降水量EOF分解的前6個(gè)特征向量貢獻(xiàn)率Table 2 Contribution of the first six eigenvectors of the EOF decomposition of annual precipitation in Sri Lanka

第1特征向量場(chǎng)(EOF1)的方差貢獻(xiàn)率為83.38%,遠(yuǎn)高于其他向量場(chǎng)的方差貢獻(xiàn)率,故斯里蘭卡降水場(chǎng)的空間分布類型以EOF1為主。由圖3(a)可知,所有站點(diǎn)的特征值均落在0.1～0.5,特征值為正則表明斯里蘭卡降水在空間變化上呈高度一致性,降水高值中心處在海島西南部,低值中心處在海島東南部。降水空間分布形勢(shì)與季風(fēng)氣候和地形密切相關(guān):斯里蘭卡受由印度洋吹向陸地的濕潤(rùn)西南季風(fēng)影響,其中強(qiáng)烈濕潤(rùn)的西南季風(fēng)在海島中南部高原地帶爬升遇冷,引起西南局部強(qiáng)降水,從而形成了斯里蘭卡降水以西南高—東南低為主的空間分布態(tài)勢(shì)。

圖3 斯里蘭卡年降水量EOF分析中的前2個(gè)特征向量場(chǎng)Fig.3 The first 2 eigenvector fields in the EOF analysis of annual precipitation in Sri Lanka

第2特征向量場(chǎng)(EOF2)是海島降水空間分布的次表現(xiàn)形式,方差貢獻(xiàn)率為6.75%。從圖3(b)可以看出,EOF2下站點(diǎn)特征值范圍處在-0.3～0.6,海島東西空間分布差異明顯。海島東部為正值區(qū),西部為負(fù)值區(qū),呈現(xiàn)明顯的東高西低的反向分布模式,反映出斯里蘭卡降水空間分布的次變化特征為由東向西依次遞減,該表現(xiàn)形式主要受海島東北干燥季風(fēng)降水影響。

3.1.2 時(shí)間分布特征

對(duì)于任意年,降水場(chǎng)的模態(tài)類型是由該年時(shí)間系數(shù)的絕對(duì)值大小來(lái)決定的,而空間模態(tài)的時(shí)間變化特征則是通過(guò)時(shí)間系數(shù)的正負(fù)來(lái)體現(xiàn)。時(shí)間系數(shù)為正值表示降水時(shí)空變化一致,負(fù)值則相反,絕對(duì)值越大表明該模態(tài)的變化規(guī)律越顯著。對(duì)兩種降水空間模態(tài)的數(shù)據(jù)(圖4)分析可知,斯里蘭卡降水場(chǎng)主要有4種表現(xiàn)類型,分別為:全域多雨、全域少雨、西南部多雨?yáng)|北部少雨和西南部少雨?yáng)|北部多雨。對(duì)2008—2020年4種降水分布類型進(jìn)行統(tǒng)計(jì)分析,得出全年多雨有7年,全年少雨有4年,全年西南部多雨有2年,具體年份如表3所示。近年來(lái)降水場(chǎng)分布主要以第1模態(tài)(全年多雨)為主,時(shí)間占比為84.6%,在第1模態(tài)下,趨勢(shì)線斜率為126.24,R2為0.092,表明近年來(lái)斯里蘭卡降水增長(zhǎng)趨勢(shì)顯著,未來(lái)發(fā)生極端降水的概率進(jìn)一步增加;模態(tài)2的時(shí)間占比為15.4%,趨勢(shì)線斜率為-25.09,R2為0.033,表明斯里蘭卡西南濕潤(rùn)區(qū)降水有減小趨勢(shì),東北干旱區(qū)降水有增加趨勢(shì)。

表3 4種降水場(chǎng)表現(xiàn)類型年份分布Table 3 Year distribution of four precipitation field performance types

圖4 前2個(gè)降水空間模態(tài)的時(shí)間系數(shù)Fig.4 Time coefficients of the first 2 spatial modes of precipitation

3.2 極端降水指數(shù)特征分析

3.2.1 時(shí)間特征分析

依據(jù)斯里蘭卡10個(gè)氣象站的降水系列數(shù)據(jù),選取五日最大降水量RX5d、普通日降水強(qiáng)度SDII、極端降水日數(shù)R20、極端降水量R95p、持續(xù)干燥指數(shù)CDD和持續(xù)濕潤(rùn)指數(shù)CWD6個(gè)極端降水指標(biāo),計(jì)算當(dāng)?shù)貥O端降水量的時(shí)空分布及其變化趨勢(shì),如圖5所示。

圖5 極端降水指數(shù)年際變化Fig.5 Interannual variation of extreme precipitation index

由圖5可以看出,RX5d、SDII、R20、R95p、CWD均呈上升趨勢(shì)。其中,R95p變化率為61.60 mm/10 a,RX5d變化率為11.80 mm/10 a,R95p和RX5d是極端降水量的重要指標(biāo)因子,這兩項(xiàng)指數(shù)顯著上升的變化趨勢(shì)能夠說(shuō)明極端降水總量及頻率將會(huì)進(jìn)一步增加,相應(yīng)的顯著上升的站點(diǎn)比重較高。SDII、R20和CWD指數(shù)變化率分別為1.66、8.05、2.29 mm/10 a,呈小幅上升趨勢(shì),但存在明顯震蕩周期;CDD指數(shù)的變化率為-18.44 mm/10 a,持續(xù)干旱情況呈明顯減少趨勢(shì)。綜合6個(gè)極端降水指數(shù)的年際變化特征,斯里蘭卡整體呈濕潤(rùn)化趨勢(shì),發(fā)生極端降水的概率和等級(jí)變大,未來(lái)面臨暴雨洪澇災(zāi)害風(fēng)險(xiǎn)增強(qiáng),應(yīng)做好防汛風(fēng)險(xiǎn)應(yīng)急預(yù)案。同時(shí)通過(guò)斯里蘭卡國(guó)內(nèi)大型水利樞紐工程,完善流域防洪工程體系,提升當(dāng)?shù)厮禐?zāi)害防御能力。

3.2.2 空間特征分析

斯里蘭卡獨(dú)特的地理位置,以及所屬的熱帶氣候造就了當(dāng)?shù)亟邓诳臻g分布上的不均衡。從圖6中可以看出,在極端降水量和降水貢獻(xiàn)率方面,6個(gè)極端降水指數(shù)在空間格局上分布差異。RX5d指數(shù)大致呈由西南向東北方向逐漸降低趨勢(shì),該指數(shù)的高值區(qū)包含加勒和拉特納普勒地區(qū);SDII指數(shù)由北向南逐漸降低,亭可馬里SDII指數(shù)多年平均值為19.53 mm/d,為區(qū)域最高值,在漢班托塔出現(xiàn)最低值,為12.13 mm/d,表明亭可馬里的極端降水頻率較高。R20、R95p指數(shù)和CWD指數(shù)空間分布情形類似,西南部地區(qū)均為高值中心,西北地區(qū)為低值區(qū)域,R95p在拉特納普勒(高值中心)與漢班托塔(低值區(qū))相差562.42 mm,在兩地間的CWD差值為9.4 d,空間差異顯著,而拉特納普勒及以西地區(qū)為斯里蘭卡降雨事件的多發(fā)區(qū)。

圖6 極端降水指數(shù)空間變化趨勢(shì)Fig.6 Spatial trends of extreme precipitation index

3.2.3 周期變化分析

運(yùn)用Morlet小波變換對(duì)各極端降水指數(shù)進(jìn)行分析,揭示極端指數(shù)的周期性變化規(guī)律[20]。圖7為小波系數(shù)實(shí)部值的高低,能清晰地反映各極端降水指數(shù)隨時(shí)間的交替變化趨勢(shì)?？梢钥闯?各項(xiàng)極端降水指數(shù)隨時(shí)間變化存在著高低振蕩現(xiàn)象,周期變化較為明顯,6個(gè)極端降水指標(biāo)除CDD外,其余指標(biāo)振蕩主周期均為6～7 a。RX5d和SDII存在3個(gè)相同的振蕩主周期,分別在2、4、7 a,表現(xiàn)為多峰譜型,其中第一主周期振蕩強(qiáng)度最大。R20和R95p兩個(gè)指標(biāo)的振蕩周期基本一致,但表現(xiàn)出更加強(qiáng)烈的振蕩強(qiáng)度,周期性趨勢(shì)也更加顯著。結(jié)合周期變化規(guī)律和前文的趨勢(shì)分析,可進(jìn)一步判斷未來(lái)10 a斯里蘭卡極端降水事件總體變化趨勢(shì)將呈現(xiàn)上升趨勢(shì)。

3.2.4 R/S時(shí)間變化趨勢(shì)預(yù)測(cè)

R/S分析無(wú)法檢驗(yàn)時(shí)間序列的變化趨勢(shì),但可根據(jù)Hurst指數(shù)H的大小來(lái)判斷時(shí)間序列未來(lái)的持續(xù)性[21-23]。運(yùn)用R/S分析法對(duì)斯里蘭卡極端降水指數(shù)進(jìn)行預(yù)測(cè),結(jié)果如表4所示。R20的H=0.5,表明歷史變化與未來(lái)變化趨勢(shì)無(wú)相關(guān)。此外除了CWD的H<0.5,其余指標(biāo)的H均大于0.5,說(shuō)明斯里蘭卡持續(xù)濕潤(rùn)指數(shù)變化將呈反持續(xù)性,而其他指標(biāo)表現(xiàn)出了持續(xù)性,未來(lái)較長(zhǎng)時(shí)間內(nèi)斯里蘭卡重點(diǎn)區(qū)域降水將保持現(xiàn)有增加的趨勢(shì)。SDII的H在0.6～0.7,表現(xiàn)出較強(qiáng)的持續(xù)性。

表4 斯里蘭卡各極端降水指數(shù)R/S趨勢(shì)分析Table 4 Analysis of R/S trend of each extreme precipitation index in Sri Lanka

由于基礎(chǔ)數(shù)據(jù)的時(shí)間序列較短,各極端降水指數(shù)的年際變化特征總體盡管表現(xiàn)出降水頻率加大的趨勢(shì)性,但其指導(dǎo)意義和科學(xué)性相對(duì)較弱[24-26]。但另一方面,小波分析結(jié)果表現(xiàn)出來(lái)的數(shù)據(jù)周期性較強(qiáng),尤其RX5d和SDII的第一主振蕩周期僅為2 a,表明在13 a的歷史時(shí)間序列中已展現(xiàn)出較高的重復(fù)性,分析結(jié)果的可靠性也更高。

結(jié)合小波分析及Hurst指數(shù)分析結(jié)果:SDII的主振蕩周期短、R/S分析的R2達(dá)到0.973,預(yù)測(cè)可信度最高;其次為R95p和RX5d兩個(gè)指標(biāo),分別具有較高的R2或較短的主振蕩周期,也具有一定的預(yù)測(cè)可信度;CDD指標(biāo)主振蕩周期不明顯,盡管R2達(dá)到0.910,其預(yù)測(cè)可信度也略低;而R20和CWD指標(biāo)的R2較低、主振蕩周期較長(zhǎng),其預(yù)測(cè)可信度也較低。

基于分析預(yù)測(cè)結(jié)果,未來(lái)斯里蘭卡的普通日降水強(qiáng)度SDII將呈現(xiàn)較強(qiáng)增加趨勢(shì),極端降水量R95p和五日最大降水量RX5d均將呈現(xiàn)微弱增加趨勢(shì),而其他極端降水指數(shù)的變化趨勢(shì)可信度相對(duì)較低。氣候演化是一個(gè)復(fù)雜的過(guò)程,斯里蘭卡未來(lái)極端降水將表現(xiàn)出持續(xù)的增加趨勢(shì),但隨機(jī)性和波動(dòng)性較大,而極端降水指數(shù)6～7 a的主周期特征可作為關(guān)注點(diǎn),基于重現(xiàn)性角度進(jìn)行水資源管理并不斷在工作中進(jìn)行優(yōu)化驗(yàn)證。

4 結(jié)論

以氣象點(diǎn)觀測(cè)資料為基礎(chǔ),分析了斯里蘭卡極端降水時(shí)空分布特點(diǎn),得到以下結(jié)論,可為當(dāng)?shù)厣鷳B(tài)環(huán)境保護(hù)和防災(zāi)減災(zāi)等工作提供參考依據(jù)。

但對(duì)于“一帶一路”合作倡議沿線國(guó)家或地區(qū),氣象數(shù)據(jù)質(zhì)量的參差不齊給研究帶來(lái)了一定難度,未來(lái)需要一些新的方法來(lái)高效處理類似降水序列缺失地區(qū)的數(shù)據(jù),同時(shí)后續(xù)還需進(jìn)一步探究極端氣候的差異性及對(duì)不同災(zāi)害事件的影響。得出如下結(jié)論。

(1)EOF分析表明,近年來(lái)斯里蘭卡降水增長(zhǎng)趨勢(shì)顯著,未來(lái)發(fā)生極端降水的概率進(jìn)一步增加。斯里蘭卡降水場(chǎng)的空間分布類型以EOF1為主,降水在空間變化上呈高度一致性,西南部地區(qū)降水增多趨勢(shì)較東北部更強(qiáng)。

(2)時(shí)間上,近10年斯里蘭卡極端降水事件呈現(xiàn)較明顯增加趨勢(shì),極端降水量指數(shù)(R95p和RX5d)的變化傾向率分別為71.03 mm/10 a和21.10 mm/10 a;在極端降水日指數(shù)(SDII、R20、CDD、CWD)中,除連續(xù)干燥日CDD呈下降趨勢(shì)外,其他指標(biāo)呈弱上升趨勢(shì)?？臻g上,斯里蘭卡極端降水指數(shù)空間差異性明顯,西南部地區(qū)均為高值中心,西北地區(qū)為低值區(qū)域,R95p在兩個(gè)區(qū)域內(nèi)相差552.81 mm,而拉特納普勒地區(qū)為斯里蘭卡降雨事件的多發(fā)區(qū)。

(3)周期變化上,各項(xiàng)極端降水指數(shù)隨時(shí)間變化存在著高低振蕩現(xiàn)象,周期性趨勢(shì)顯著,RX5d和SDII存在3個(gè)相同的振蕩主周期,分別在2、4、7 a,表現(xiàn)為多峰譜型,其中第一主周期振蕩強(qiáng)度最大。

(4)綜合R/S分析及各指標(biāo)的周期變化規(guī)律,未來(lái)斯里蘭卡的普通日降水強(qiáng)度SDII將呈現(xiàn)較強(qiáng)增加趨勢(shì),極端降水量R95p和五日最大降水量RX5d均將呈現(xiàn)微弱增加趨勢(shì),而其他極端降水指數(shù)的變化趨勢(shì)可信度相對(duì)較低。表斯里蘭卡未來(lái)極端降水總體將表現(xiàn)出持續(xù)的增加趨勢(shì),極端降水指數(shù)6～7 a的主周期特征可作為當(dāng)?shù)厮Y源管理的關(guān)注點(diǎn)并不斷在工作中進(jìn)行優(yōu)化驗(yàn)證。