基于SPI指數(shù)的云南省干旱時(shí)空特征分析

唐 怡，卯昌書(shū)，蘇建廣

(云南省水利水電勘測(cè)設(shè)計(jì)研究院，云南 昆明 650021)

0 引 言

干旱一般分為氣象干旱、農(nóng)業(yè)干旱、水文干旱以及經(jīng)濟(jì)社會(huì)干旱等，氣象干旱是指某時(shí)段內(nèi)由于蒸散量和降水量的收支不平衡，水分支出大于水分收入而造成地表水分短缺的現(xiàn)象，是其他類(lèi)型干旱的起因和監(jiān)測(cè)評(píng)估的基礎(chǔ)[1]。氣象干旱指數(shù)是表征氣象干旱特征的定量化指標(biāo)，目前國(guó)內(nèi)外已發(fā)展的用以描述干旱事件的干旱指數(shù)近上百種[2]，常用的氣象干旱指數(shù)主要有降水距平百分率(Pa)、標(biāo)準(zhǔn)化降水指數(shù)(SPI)[3]、帕爾默旱度指數(shù)(PDSI)、相對(duì)濕潤(rùn)度指數(shù)(MI)、綜合氣象干旱指數(shù)(CI)和標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)等。其中SPI指數(shù)由于資料獲取相對(duì)簡(jiǎn)單、計(jì)算成果穩(wěn)定、適用性較好，被廣泛應(yīng)用于干旱研究中[4-6]。

干旱一直是影響中國(guó)農(nóng)業(yè)生產(chǎn)與經(jīng)濟(jì)發(fā)展的主要自然災(zāi)害，也是國(guó)內(nèi)學(xué)者的研究熱點(diǎn)。李楊等[7]通過(guò)分析各年代干旱頻率重心的變化得出，中國(guó)的干旱重心由西北部逐步轉(zhuǎn)移到了西南部地區(qū)。賀晉云等[8]研究發(fā)現(xiàn)，21世紀(jì)以來(lái)，西南地區(qū)年極端干旱的發(fā)生次數(shù)有所增加，區(qū)域間差異卻顯著減小。趙蘭蘭等[9]認(rèn)為西南地區(qū)干旱呈明顯加重趨勢(shì)，其中云南西南部氣象干旱呈顯著加重趨勢(shì)；楊蕊等[10]研究了云南省冬春干旱，認(rèn)為滇中是冬春干旱的主要發(fā)生區(qū)，冬春干旱中重旱與中旱的范圍呈緩慢擴(kuò)大趨勢(shì)。

云南省干旱頻發(fā)，近年來(lái)呈現(xiàn)出干旱范圍擴(kuò)大、程度加重的趨勢(shì)，對(duì)社會(huì)經(jīng)濟(jì)造成嚴(yán)重影響。為全面了解云南省氣象干旱的時(shí)空演變特征，本文基于云南省104個(gè)氣象站點(diǎn)1959～2018年的逐月降雨數(shù)據(jù)，采用標(biāo)準(zhǔn)化降水指數(shù)(SPI)對(duì)60 a來(lái)云南省干旱的強(qiáng)度變化、空間分布、周期波動(dòng)等進(jìn)行分析，旨在為云南省干旱監(jiān)測(cè)和防災(zāi)減災(zāi)提供參考。

1 研究區(qū)域概況

云南省地處中國(guó)西南邊陲，簡(jiǎn)稱(chēng)滇或云，位于東經(jīng)97°31′～106°11′，北緯21°08′～29°15′，跨越10個(gè)經(jīng)度帶和9個(gè)緯度帶，北回歸線橫穿南部。云南省總面積39.41萬(wàn)km2，占全國(guó)陸地總面積的4.1%，居全國(guó)第8 位。云南省地處低緯度高原，冬季受干燥大陸季風(fēng)控制，夏季盛行濕潤(rùn)海洋季風(fēng)，屬低緯山原季風(fēng)氣候。全省氣候類(lèi)型多樣，由南至北依次分布有北熱帶、南亞熱帶、中亞熱帶、北亞熱帶、南溫帶、中溫帶和高原寒溫帶7個(gè)氣候類(lèi)型。全省多年平均降水量1 278.8 mm，多年平均水面蒸發(fā)量1 057.0 mm，年均日照時(shí)數(shù)1 000～2 800 h。全省年溫差小，年平均溫度5～24 ℃左右；日溫差大，早晚較涼，中午較熱，日溫差達(dá)12～20 ℃。云南省的氣候特點(diǎn)有利于各類(lèi)農(nóng)作物生長(zhǎng)，但干季和濕季過(guò)于集中，垂直變化突出，干旱、洪澇、低溫冷凍、霜凍、冰雹與風(fēng)災(zāi)等自然災(zāi)害頻發(fā)。

2 研究數(shù)據(jù)及方法

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

從云南省126個(gè)氣象站點(diǎn)中篩選出104個(gè)資料系列滿足要求、代表性較好的站點(diǎn)，對(duì)其1959～2018年共計(jì)60 a的逐月降雨量數(shù)據(jù)進(jìn)行分析，部分站點(diǎn)缺失數(shù)據(jù)進(jìn)行插值處理。研究區(qū)域范圍和站點(diǎn)分布情況見(jiàn)圖1。

圖1 研究區(qū)域及站點(diǎn)分布Fig.1 Study area and site distribution

2.2 研究方法

2.2.1 標(biāo)準(zhǔn)化降水指數(shù)SPI

標(biāo)準(zhǔn)化降水指數(shù)SPI 通過(guò)計(jì)算某時(shí)段內(nèi)的降水量Γ分布概率，并進(jìn)行正態(tài)標(biāo)準(zhǔn)化處理后，用標(biāo)準(zhǔn)化降水累積頻率分布來(lái)劃分干旱等級(jí)。SPI 計(jì)算方法如下[1]：

(1)

F(x)由以下Γ分布函數(shù)概率密度積分公式得到：

(2)

式中：r，β分別為Γ分布函數(shù)的形狀和尺度參數(shù)。

GB/T 20481-2017《氣象干旱等級(jí)標(biāo)準(zhǔn)》中，依據(jù)標(biāo)準(zhǔn)化降水指數(shù)劃分的干旱等級(jí)見(jiàn)表1。

表1 標(biāo)準(zhǔn)化降水指數(shù)干旱等級(jí)劃分

2.2.2 Mann-Kendall(M-K)檢驗(yàn)及小波分析

Mann-Kendall(M-K)檢驗(yàn)法是世界氣象組織推薦的非參數(shù)檢驗(yàn)方法，常用于氣象、水文等領(lǐng)域，可檢驗(yàn)時(shí)間序列的趨勢(shì)變化。本文采用M-K檢驗(yàn)法分析不同尺度SPI系列的突變性與趨勢(shì)性。

小波分析是Morlet 在20 世紀(jì)80 年代初提出的一種具有時(shí)頻多分辨功能的分析工具，能夠揭示時(shí)間序列中的多種變化周期，可以反映系統(tǒng)在不同時(shí)間尺度中的變化趨勢(shì)，近年來(lái)已廣泛應(yīng)用于時(shí)間序列分析[11-12]。本文采用Morlet小波分析不同尺度SPI系列的周期變化規(guī)律。

2.2.3 干旱特征指標(biāo)

(1) 干旱頻率(Pi)。Pi表示干旱發(fā)生的頻繁程度，計(jì)算公式為

Pi=n/N×100%

(3)

式中：Pi為某站點(diǎn)的干旱頻率；n為某站點(diǎn)發(fā)生不同程度干旱的年數(shù)；N為某站點(diǎn)計(jì)算總年數(shù)。

(2) 干旱站次比(Pj)。Pj用于評(píng)價(jià)干旱影響范圍的大小，計(jì)算公式為

Pj=m/M×100%

(4)

式中：Pj為干旱站次比；m為發(fā)生不同程度干旱的站點(diǎn)數(shù)；M為總站點(diǎn)數(shù)。當(dāng)Pj≥50% 時(shí)，為全域性干旱；當(dāng) 50%>Pj≥33%時(shí)，為區(qū)域性干旱；當(dāng) 33%>Pj≥25%時(shí)，為部分區(qū)域性干旱；當(dāng) 25%>Pj≥10%時(shí)，為局域性干旱；當(dāng)Pj<10%時(shí)，為無(wú)明顯干旱[11]。

(3) 干旱強(qiáng)度(Sij)。Sij用于評(píng)價(jià)干旱嚴(yán)重程度，單站的某時(shí)段內(nèi)的干旱強(qiáng)度可由SPI的絕對(duì)值表示，SPI絕對(duì)值越大，表示干旱越嚴(yán)重。Sij計(jì)算公式為

(5)

式中：SPIi為發(fā)生干旱時(shí)SPI的絕對(duì)值。當(dāng)1>Sij≥0.5時(shí)，為輕度干旱；當(dāng)1.5>Sij≥1時(shí)，為中度干旱；當(dāng)2>Sij≥1.5時(shí)，為重度干旱；當(dāng)Sij≥2時(shí)，為特旱。

3 結(jié)果分析

3.1 時(shí)間變化特征

3.1.1 多尺度SPI時(shí)間變化特征

根據(jù)云南省面雨量計(jì)算月尺度、季尺度和年尺度的SPI值(圖2)，即SPI1，SPI3，SPI12，從圖2可以看出SPI1受短期降雨影響大，在0值上下頻繁轉(zhuǎn)換，波動(dòng)幅度最大，SPI3變化幅度逐漸減緩，SPI12綜合前期降雨變化的累計(jì)影響，變化幅度最小，較為穩(wěn)定。

圖2 云南省多尺度SPI時(shí)間變化趨勢(shì)Fig.2 Multi-scale SPI temporal variation trend in Yunnan Province

從M-K檢驗(yàn)結(jié)果來(lái)看，月尺度的突變點(diǎn)在1961，2000年附近，1959～1967年總體呈下降趨勢(shì)，1960年下降趨勢(shì)明顯，1968～1979年呈上升趨勢(shì)，1980～2018年總體呈下降趨勢(shì)，2008年以后下降趨勢(shì)明顯。季尺度的突變點(diǎn)在1961，2002年附近，1959～1966年總體呈下降趨勢(shì)，1960年下降趨勢(shì)明顯，1967～1979年呈上升趨勢(shì)，1980～1996年總體呈下降趨勢(shì)，1997～2009年呈上升趨勢(shì)，2009年以后呈下降趨勢(shì)，且下降趨勢(shì)明顯。年尺度的突變點(diǎn)在1961，2003年附近，1959～1966年總體呈下降趨勢(shì)，1960年下降趨勢(shì)明顯，1967～1979年呈上升趨勢(shì)，1980～2000年總體呈下降趨勢(shì)，2001～2009年呈上升趨勢(shì)，2009年以后呈下降趨勢(shì)，且下降趨勢(shì)明顯。各尺度M-K檢驗(yàn)結(jié)果均反映出1959～1966年、1980～1996年、2009～2018年云南省干旱呈加重趨勢(shì)，特別是2009年以后加重趨勢(shì)明顯。

對(duì)不同尺度SPI值進(jìn)行小波分析計(jì)算，結(jié)果顯示各尺度SPI值均呈現(xiàn)出一定的周期性變化規(guī)律(圖3)。從小波系數(shù)實(shí)部圖(圖4)看出：SPI1值在5～10 a，10～15 a，20～30 a時(shí)間尺度上均有較明顯的周期變化，其中20～30 a尺度表現(xiàn)最明顯，中心時(shí)間尺度為27 a左右，旱澇交替變化；SPI3值在5～10 a，10～15 a，20～30 a時(shí)間尺度上均有較明顯的周期變化，其中20～30 a尺度表現(xiàn)最明顯，中心時(shí)間尺度為26 a左右，旱澇交替變化；SPI12值在10～15 a，15～30 a時(shí)間尺度上均有較明顯的周期變化，其中15～30 a尺度表現(xiàn)最明顯，中心時(shí)間尺度為24 a左右，旱澇交替變化。由圖3可知：SPI1在4，7，14，27 a附近存在震蕩周期，其中27 a為主周期；SPI3在14，26 a附近存在震蕩周期，其中26 a為主周期；SPI12在7，24 a附近存在震蕩周期，其中24 a為主周期。

圖3 不同尺度小波方差Fig.3 Square difference of wavelets of different scales

圖4 不同尺度小波系數(shù)實(shí)部圖Fig.4 Real part diagram of wavelet coefficients at different scales

3.1.2 干旱強(qiáng)度時(shí)間變化特征

根據(jù)各站降雨量，計(jì)算年尺度、季尺度的SPI值。從年尺度來(lái)看，1959～2018年間云南省干旱強(qiáng)度在0.59～1.77之間變化，平均為1.03；基本上每年都發(fā)生干旱，大部分為中、輕度干旱，2009，2011年發(fā)生重度干旱(圖5)。

圖5 年尺度干旱強(qiáng)度趨勢(shì)Fig.5 Trend of drought intensity at annual scale

從季節(jié)來(lái)看，春季干旱強(qiáng)度在0～2.13之間變化，平均為0.90，多為中度以下干旱，干旱強(qiáng)度呈緩慢下降趨勢(shì)。夏季干旱強(qiáng)度在0.77～1.73之間變化，平均為1.04，多為中度以下干旱，干旱強(qiáng)度呈緩慢上升趨勢(shì)。秋季干旱強(qiáng)度在0.54～2.02之間變化，平均為1.04，多為中度以下干旱，干旱強(qiáng)度呈緩慢上升趨勢(shì)。冬季干旱強(qiáng)度在0～1.84之間變化，平均為0.96，多為中度以下干旱，干旱強(qiáng)度呈緩慢下降趨勢(shì)(圖6)。

圖6 季尺度干旱強(qiáng)度趨勢(shì)Fig.6 Trend of seasonal drought intensity

3.1.3 干旱站次比時(shí)間變化特征

從年尺度的干旱站次比來(lái)看，云南省多為部分區(qū)域性干旱，干旱站次比多年平均為30.8%。1959～2000年間干旱站次比大于50%的有7次，全省平均每10 a有2次全域性干旱，2000年以后全域性干旱次數(shù)明顯增加，2000～2018年共發(fā)生6次，2009～2012年連續(xù)4 a發(fā)生了全域性干旱，其中2009年干旱站次比高達(dá)86.5%(圖7)。

從不同干旱強(qiáng)度的站次比來(lái)看，輕旱站次比最高，多年平均為14.4%、其次為中旱，為9.31%，重旱、特旱分別為4.73%，2.39%，各種程度的干旱均呈現(xiàn)緩慢上升趨勢(shì)；干旱站次比越高的年份，重旱、特旱比例越大，以2011年為例，干旱站次比83.65%，特旱、重旱站次比分別為31.7%，19.2%。這說(shuō)明發(fā)生干旱的范圍越大，干旱的強(qiáng)度越大(圖7)。

圖7 年尺度干旱范圍(站次比)趨勢(shì)Fig.7 Trend of annual scale drought extent (station’s drought proportions)

從季節(jié)尺度的干旱站次比來(lái)看，春季干旱站次比在0～91.2%之間變化，平均為29.8%，總體呈下降趨勢(shì)；夏季干旱站次比在2.88%～87.5%之間變化，平均為30.3%，總體呈上升趨勢(shì)；秋季干旱站次比在0～92.3%之間變化，平均為29.7%，總體呈上升趨勢(shì)；冬季干旱站次比在0～94.2%之間變化，平均為28.6%，總體呈下降趨勢(shì)(圖8)。

圖8 不同季節(jié)干旱范圍(站次比)趨勢(shì)Fig.8 Trend of drought range (station’s drought proportions) in different seasons

3.2 干旱空間分布特征

近60 a來(lái)，云南省各站點(diǎn)發(fā)生干旱的頻率為21.7%～40.0%，平均為30.8%。其中發(fā)生輕旱的頻率在3.3%～23.3%之間，平均為14.4%；發(fā)生中旱的頻率在3.3%～20.0%之間，平均為9.3%；發(fā)生重旱的頻率在0～13.3%之間，平均為4.7%；發(fā)生特旱的頻率在0～6.7%之間，平均為2.4%(圖9)。發(fā)生輕旱頻率較高的州市有德宏州、保山市、臨滄市、昭通市，主要集中在長(zhǎng)江、瀾滄江、怒江流域；發(fā)生中旱頻率較高的有楚雄州、昆明市、西雙版納州、文山州，主要集中在長(zhǎng)江、瀾滄江、紅河流域；發(fā)生重旱頻率較高的有玉溪市、普洱市、臨滄市、怒江州，主要集中在紅河、瀾滄江、怒江流域；發(fā)生特旱頻率較高的有昭通市、曲靖市、玉溪市、紅河州，主要集中在長(zhǎng)江、珠江、紅河流域。整體來(lái)看，全省干旱以中、輕度為主，重旱以上干旱發(fā)生的頻率不高，但基本每個(gè)州市都會(huì)發(fā)生。

圖9 云南省干旱頻率(單位：%)Fig.9 Drought frequency map of Yunnan Province

圖10 云南省不同時(shí)段干旱頻率(單位：%)Fig.10 Drought frequency map of Yunnan Province in different years

從不同年代的干旱頻率(圖10)來(lái)看，20世紀(jì)60～70年代云南省干旱發(fā)生頻率較高的地區(qū)主要分布在昆明、楚雄、麗江、臨滄、文山、德宏等州市，干旱頻率為28%～33%；20世紀(jì)80～90年代主要分布在昆明、玉溪、楚雄、西雙版納、昭通等州市，干旱頻率為30%～37%之間；進(jìn)入21世紀(jì)以后，德宏、保山、文山、紅河等州市成為干旱頻率發(fā)生較高的地區(qū)，近20 a干旱頻率為40%～46%。

4 結(jié) 論

本文利用云南省104個(gè)氣象站點(diǎn)1959～2018年共60 a的逐月降雨量數(shù)據(jù)，計(jì)算出云南省面雨量不同尺度的SPI值及各站點(diǎn)不同尺度的SPI值，分析云南省干旱的時(shí)空變化特征，得到如下結(jié)論。

(1) 不同尺度的SPI值隨時(shí)間尺度的增加波動(dòng)逐漸減小，M-K檢驗(yàn)結(jié)果反映出1959～1966年、1980～1996年、2009～2018年云南省干旱呈加重趨勢(shì)，特別是2009年以后加重趨勢(shì)明顯。不同尺度的SPI值呈現(xiàn)出一定的周期性變化特征，其中月尺度以27 a為主周期，季尺度以26 a為主周期；年尺度以24 a為主周期。

(2) 從時(shí)間上來(lái)看，云南省基本每年都會(huì)有地區(qū)發(fā)生干旱，大部分以輕旱、中旱為主，但重旱、特旱發(fā)生的比例也在增加；近60 a來(lái)干旱強(qiáng)度和范圍均呈緩慢上升趨勢(shì)，且干旱的范圍越大，干旱的強(qiáng)度越強(qiáng)。近60 a來(lái)，云南省夏、秋季節(jié)干旱強(qiáng)度、干旱范圍均呈上升趨勢(shì)，春、冬季節(jié)呈下降趨勢(shì)。

(3) 從空間上來(lái)看，云南省重度以上干旱在各州市均有發(fā)生，極端干旱主要集中在長(zhǎng)江、珠江、紅河流域；近年來(lái)，文山、紅河等滇東南區(qū)，德宏、保山等滇西南區(qū)成為云南省干旱高發(fā)區(qū)。