胡 悅，杜靈通，候 靜，劉 可，朱玉果

(1.寧夏大學(xué)西北土地退化與生態(tài)系統(tǒng)恢復(fù)省部共建國家重點實驗室培育基地， 寧夏 銀川 750021；2.寧夏大學(xué)西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點實驗室， 寧夏 銀川 750021)

基于SPI指數(shù)的寧夏中部干旱帶1960—2012年干旱特征研究

胡 悅1,2，杜靈通1,2，候 靜1,2，劉 可1,2，朱玉果1,2

(1.寧夏大學(xué)西北土地退化與生態(tài)系統(tǒng)恢復(fù)省部共建國家重點實驗室培育基地， 寧夏 銀川 750021；2.寧夏大學(xué)西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點實驗室， 寧夏 銀川 750021)

研究寧夏中部干旱帶的干旱特征可為區(qū)域農(nóng)牧業(yè)發(fā)展和防災(zāi)減災(zāi)提供科學(xué)決策依據(jù)。利用1960—2012年的各站點氣象觀測資料，計算標準化降水指數(shù)(SPI)，并借助集合經(jīng)驗?zāi)B(tài)分解、功率譜分析和重新標度極差分析等方法，對近50年來的干旱特征和變化趨勢進行了研究。研究結(jié)果表明，寧夏中部干旱帶在1960—2012年間發(fā)生中旱以上的干旱事件達到了總干旱次數(shù)的30%以上；干旱波動周期的多尺度分解顯示，該區(qū)域存在著0.5 a尺度的輕旱，1 a尺度的中旱，3 a尺度的重旱和10 a尺度的特旱；平均干旱強度的年際波動周期特征中，有平均超過65%的貢獻來自于3.25-4 a左右的干旱周期波動，而平均超過14%的貢獻來自于5.78-6.5 a左右的干旱周期波動；總體來看，寧夏中部干旱帶的干旱趨勢處于增強態(tài)勢，且具有長程依賴性，即未來還將持續(xù)干旱，但不同季節(jié)的干旱趨勢特征又表現(xiàn)出一定的差異。

干旱；標準化降水指數(shù)；集合經(jīng)驗?zāi)B(tài)分解；功率譜分析；寧夏中部干旱帶

全球氣候變暖現(xiàn)已成為一個不爭的事實，氣候變暖使得各類災(zāi)害頻發(fā)，嚴重威脅著人類社會的經(jīng)濟和環(huán)境的可持續(xù)發(fā)展。在各類災(zāi)害中，旱災(zāi)是我國目前主要的農(nóng)業(yè)氣象災(zāi)害，平均每年受災(zāi)面積最高可達2 200萬hm2[1]。寧夏是中國西部典型的生態(tài)脆弱區(qū)和氣候變化敏感區(qū)[2]，干旱半干旱的氣候特征導(dǎo)致了寧夏干旱災(zāi)害頻發(fā)，而中部干旱帶又是寧夏干旱災(zāi)害最為嚴重的區(qū)域，頻發(fā)的氣象干旱事件常常導(dǎo)致該地區(qū)的農(nóng)牧生產(chǎn)受損。因此，研究其干旱變化特征和未來演化趨勢，不僅可為地方政府抗旱行動提供決策支持，而且可為區(qū)域生態(tài)建設(shè)提供理論依據(jù)。

旱災(zāi)因其具有影響范圍廣、持續(xù)時間長、發(fā)生頻率高的特點，引起了國內(nèi)外學(xué)者的廣泛關(guān)注，并開發(fā)了一系列的監(jiān)測指標對區(qū)域干旱災(zāi)害開展研究。目前國內(nèi)外發(fā)展成熟的氣象干旱監(jiān)測和評估指標有標準化降水指數(shù)(Standardized Precipitation Index，SPI)、帕默爾干旱指數(shù)(Palmer Drought Severity Index，PDSI)和地表濕潤指數(shù)等。馬柱國等[3]利用地表濕潤指數(shù)對我國北方極端干旱的分布特征進行了研究，發(fā)現(xiàn)東北和華北兩個地區(qū)在近10年中極端干旱頻率顯著增加。白樺等[4]使用12月標準化降水指數(shù)(SPI12)評估了渭河流域氣象干旱程度及水文旱澇演變規(guī)律。黃妙芬[5]用帕默爾干旱指數(shù)分析了中國黃土高原西北部地區(qū)的干旱特征。現(xiàn)有的研究和實踐應(yīng)用證明，標準化降水指數(shù)是一種簡單易行且監(jiān)測效率較高的指數(shù)，適宜于不同氣候區(qū)的干旱監(jiān)測與變化特征分析。

為了挖掘時間序列氣候數(shù)據(jù)的內(nèi)在特征和變化趨勢[6]，近年來發(fā)展起了許多新的分析方法，集合經(jīng)驗?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition，EEMD)，通過對時間序列數(shù)據(jù)的有效分解，可從不同尺度上探索時間序列數(shù)據(jù)的內(nèi)在波動特征和趨勢變化[7-8]；以傅里葉變換為基礎(chǔ)的功率譜分析能夠有效判斷時序數(shù)據(jù)的波動周期[6]；而重新標度極差(Rescaled Range Analysis，R/S)分析構(gòu)建的赫斯特指數(shù)(Hurst，H)，能用于判斷時間序列數(shù)據(jù)是否具有自相似和長程依賴性，廣泛用于分析預(yù)測氣候要素的未來變化趨勢[9]。這些氣候診斷技術(shù)和方法，為現(xiàn)在的氣候?qū)W研究，特別是干旱變化特征研究提供了強有力的手段。為此，本研究基于標準化降水指數(shù)和這些先進的氣候診斷技術(shù)與分析方法，對寧夏中部干旱帶的氣候干旱內(nèi)在特征和演變趨勢展開研究，以期為區(qū)域農(nóng)牧業(yè)和經(jīng)濟社會發(fā)展提供科學(xué)依據(jù)。

1 研究區(qū)概況、數(shù)據(jù)及方法

1.1 研究區(qū)概況

寧夏位于中國西北地區(qū)的東部，處在東部季風(fēng)區(qū)、西北干旱區(qū)和青藏高寒區(qū)的交匯過渡帶。寧夏的氣候特征南北差異較大，根據(jù)自然地理特征和農(nóng)牧業(yè)生產(chǎn)特征可將全省劃分為北部引黃灌區(qū)、中部干旱帶和南部山區(qū)三個地理單元(圖1)。中部干旱帶在行政區(qū)劃上包括鹽池、同心和紅寺堡3個縣(區(qū))的全境及中寧、靈武、海原、利通區(qū)、沙坡頭區(qū)等縣(市、區(qū))的部分區(qū)域，由于氣候為典型大陸性氣候，年平均降水量在200～400 mm之間，日照強烈、蒸發(fā)量大，故被稱為中部干旱帶。中部干旱帶是農(nóng)牧交錯區(qū)，自然植被主要為干草原和荒漠草原，農(nóng)業(yè)以旱作為主，輔以少量揚黃灌溉和井灌區(qū)。地帶性土壤為棕鈣土和灰鈣土，腐殖質(zhì)層較薄，肥力低。由于降水量少，且年內(nèi)降水分布不均勻，常發(fā)生氣象干旱，干旱災(zāi)害也是中部干旱帶最主要的自然災(zāi)害之一，對區(qū)域農(nóng)牧業(yè)生產(chǎn)影響很大。

圖1 研究區(qū)概況圖

Fig.1 Location map of the study area

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

寧夏中部干旱帶上有鹽池、中寧、同心、韋州和興仁堡5個常規(guī)氣象站，其中韋州和興仁堡的數(shù)據(jù)序列較短，且區(qū)域代表性不強。雖然中寧縣境內(nèi)大部分屬于引黃灌區(qū)，但由于中寧氣象站位于縣境南側(cè)，其降水特點和地理位置更接近于中部干旱帶的鹽池和同心，因此也作為中部干旱帶的典型站點[10]。本文選取了基本均勻分布于研究區(qū)的鹽池、中寧和同心3個典型氣象站開展研究，從中國氣象科學(xué)數(shù)據(jù)共享網(wǎng)(http://www.escience.gov.cn/)獲取了3個氣象站1960年1月到2012年12月的逐月氣象觀測數(shù)據(jù)，觀測要素包括降水量和平均氣溫。

1.3 研究方法

1.3.1 標準化降水指數(shù) 標準化降水指數(shù)(Standardized Precipitation Index，SPI)是一種應(yīng)用廣泛的氣象干旱指數(shù)，他是表征某時段降水量出現(xiàn)的概率多少的指標，適合于月以上時間尺度的干旱監(jiān)測與評估。由于SPI采用Γ函數(shù)的標準化降水累積頻率分布來描述降水量的變化，因此其值在不同地區(qū)和不同時間段之間具有可比性。由于自然降水量是一種偏態(tài)分布，所以在計算標準化降水指數(shù)時，先將偏態(tài)概率分布的降水量進行了正態(tài)標準化處理，再用標準化降水累積頻率分布來劃分干旱等級[11-12]，計算公式如下[13]：

(1)

(2)

式中，γ、β分別為形狀和尺度參數(shù)。

本研究利用鹽池、中寧和同心3個縣1960—2012年的逐月降水資料，計算了逐月標準化降水指數(shù)，并按照SPI數(shù)值的大小劃分出不同的干旱等級，根據(jù)《氣象干旱等級GB/T20481-2006》中的分級標準分為特旱(SPI≤-2.0)、重旱(-2.0

1.3.2 集合經(jīng)驗?zāi)B(tài)分解方法 集合經(jīng)驗?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition，EEMD)是基于經(jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)發(fā)展起來的一種自適應(yīng)時間序列分析技術(shù)，其改進了經(jīng)驗?zāi)B(tài)分解中由于信號間隔不連續(xù)而引起的模態(tài)混疊問題，適用于非線性、非平穩(wěn)的時間序列數(shù)據(jù)分析[14]，如分析時間序列氣候要素，提取干旱災(zāi)害等氣候變化的信息[15-17]。EEMD在原始信號中加入若干次白噪聲，把信號和噪聲的組合作為一個待分解信號，再利用EMD將原始數(shù)據(jù)分解成不同尺度的內(nèi)在模態(tài)分量(Intrinsic Mode Function，IMF)和一個殘余趨勢分量(RES)，步驟如下[18]：

1) 通過給目標信號序列x(t)加上一組有限振幅的白噪聲w(t)來獲得一個總體序列X(t)：

X(t)=x(t)+w(t)

(3)

2) 對X(t)進行EMD分解，得到第一個IMF分量c1，代表原始序列中最高頻的分量。用X(t)減去c1，得到去掉高頻成分的殘余序列r1，同理可得c2，c3，…，cn，直到rn滿足預(yù)先給定的終止準則，分解終止，得到X(t)：

(4)

3) 將不同的白噪聲wi(t)加入到目標信號中，并重復(fù)以上步驟：Xi(t)=x(t)+wi(t)，分解后得到各自的IMF分量組：

(5)

4) 取相應(yīng)IMF的均值作為最終的IMF組：

(6)

其中，N為總體的個數(shù)。

本研究利用集合經(jīng)驗?zāi)B(tài)分解法對1960—2012年鹽池、中寧和同心縣的SPI數(shù)據(jù)進行分解，獲取干旱信號的時間尺度特征，分解過程中將擾動白噪聲與原始信號的信噪比設(shè)置為0.2，集合平均次數(shù)設(shè)置為150。

1.3.3 功率譜分析方法 功率譜是一種提取序列數(shù)據(jù)周期的分析方法，用于研究信號在頻域中的各種特征。功率譜分析是以傅里葉變換為基礎(chǔ)的頻域分析方法，將時間序列數(shù)據(jù)的總能量分解到不同頻率的分量上，根據(jù)不同頻率波的方差貢獻診斷出序列數(shù)據(jù)的主周期，從而確定序列數(shù)據(jù)隱含的各個顯著周期[6]。目前功率譜分析方法已在氣象分析和氣候變化研究中得到廣泛應(yīng)用[19-22]，計算公式如下[23]：

1) 自相關(guān)系數(shù)計算：

(7)

(8)

(9)

2) 粗普估計：

(10)

1.3.4 赫斯特指數(shù) 重新標度極差分析方法是組成分形理論的重要部分，最早由水文學(xué)家Hurst在研究尼羅河水文流量時提出，通過計算赫斯特指數(shù)(H)來分析時間序列數(shù)據(jù)的分形特征和長期記憶過程[24]，計算公式如下[25]：

考慮一個時間序列x(t)，t=1，2，…，N，對于任意正整數(shù)n≥1，定義均值序列：

(11)

1) 計算累積離差：

(12)

2) 計算極差：

R(n)=max1

(13)

3) 計算標準差：

(14)

考慮比值R(n)/S(n)?R/S，若R/S∝nH存在，則說明時間序列x(t)存在Hurst現(xiàn)象，H為Hurst指數(shù)，由n、R/S的對數(shù)序列通過最小二乘法擬合得到，值在0-1之間。0

2 結(jié)果與分析

2.1 干旱強度及頻次變化

在計算1960—2012年鹽池、中寧和同心3縣逐月SPI的基礎(chǔ)上，依據(jù)SPI值大小劃分不同的干旱等級，統(tǒng)計近53年間不同干旱等級發(fā)生的頻次(表1)。結(jié)果表明，鹽池縣共發(fā)生輕旱以上等級干旱215次；中寧縣共發(fā)生輕旱以上等級干旱206次；同心縣共發(fā)生輕旱以上等級干旱189次。從發(fā)生干旱的等級來看，近53年鹽池、中寧和同心3縣發(fā)生的干旱事件中，超過中旱以上的干旱事件達到了總干旱次數(shù)的30%以上；其中中旱平均發(fā)生周期在0.93—1.2a之間，基本上為年年有中旱，特大干旱發(fā)生周期在8.83—17.67a之間，基本上為十年一大旱。盡管不同等級干旱的發(fā)生在近53年間存在不均性，但從總體干旱強度和發(fā)生頻率來看，中部干旱帶各縣是寧夏發(fā)生干旱頻率較高的地區(qū)，這與前人通過空間遙感數(shù)據(jù)研究的結(jié)果一致[26]。

表1 1960—2012年干旱情況統(tǒng)計

2.2 干旱災(zāi)害多尺度變化特征

采用集合經(jīng)驗?zāi)B(tài)分解分別對寧夏中部干旱帶的鹽池、中寧和同心3縣1960—2012年的逐月SPI數(shù)據(jù)進行逐步分解，3縣均可分解獲得8個具有不同波動周期的固有模態(tài)函數(shù)分量(IMF1-8)和一個殘余趨勢分量(RES)(圖2)。其中各IMF分量分別包含了由高到低的不同頻率信息，即鹽池、中寧和同心3縣的SPI序列數(shù)據(jù)包含有多個時間尺度的波動特征，且各個IMF組分分別體現(xiàn)原始序列的不同波動周期的局部化特性，最后分解殘余的趨勢項表示的是原始數(shù)據(jù)序列總體上隨時間變化的趨勢。在8個IMF分量中，由于IMF1代表的是EEMD分解時模型添加的白噪聲信息，所以圖2中剔除高頻分量IMF1。從圖中可以看出，其余各IMF分量均有其相對穩(wěn)定的準周期，在相同的時間段內(nèi)，不同時間尺度的準周期振蕩隨著時間也呈現(xiàn)出或強或弱的非均勻變化。

圖2 1960—2012年逐月SPI的EEMD分解結(jié)果

Fig.2EEMDdecompositionresultsofmonthlySPIfrom1960to2012

通過求取IMF2-8各分量的平均周期，發(fā)現(xiàn)鹽池、中寧和同心3縣的SPI均存在0.5～35a的多種時間尺度波動周期(表2)。其中IMF2分量的波動周期為0.55～0.57a，與發(fā)生輕旱的頻率(0.37～0.48a)基本一致，故可推斷IMF2分量為發(fā)生輕旱的波動；IMF3分量的波動周期為1.12～1.23a，與發(fā)生中旱的頻率(0.93～1.20a)基本一致，因此其可能為發(fā)生中旱的波動；IMF4則是發(fā)生重旱的波動；而約十年一遇的特旱則與IMF6分量的周期基本一致。研究計算了每個IMF分量對原始信號的方差貢獻率，用于定量化表征每種尺度信號波動頻率和振幅對原始數(shù)據(jù)總體特征的影響程度，同時用皮爾遜相關(guān)系數(shù)來衡量每個IMF分量和原始序列的相關(guān)性，結(jié)果如表2所示。從中可以看出，IMF2的方差貢獻率最高，隨著IMF分量周期的增大，其方差貢獻率也相應(yīng)的降低，同時IMF分量與原始序列的相關(guān)系數(shù)也隨著降低。這一結(jié)果表明，寧夏中部干旱帶鹽池、中寧和同心3縣的干旱災(zāi)害存在著0.5a尺度的輕旱，1a尺度的中旱，3a尺度的重旱和10a尺度的特旱。

2.3 基于功率譜的年際干旱周期分析

通過EEMD分解后獲取的IMF各分量的波動周期，是從IMF分量的波譜曲線中主觀判別出來的，由于各IMF分量波譜并不是具有穩(wěn)定周期和振幅的波動函數(shù)，因此計算出的IMF分量的周期存在一定的主觀性。為獲取年際干旱的準確周期，本研究利用功率譜分析方法，從時序信號的頻率域?qū)}池、中寧和同心3縣的干旱發(fā)生周期進行定量分析。由于年平均SPI序列數(shù)據(jù)的IMF1和IMF2分量的方差貢獻率最大(表3)，因此，利用功率譜重點分析了鹽池、中寧和同心3縣的年平均SPI原始序列以及IMF1、IMF2兩個高頻分量序列。

表2 中部干旱帶各縣逐月SPI序列EEMD分解的各IMF分量特征

注(Note)：*P<0.05，**P<0.01。

表3 年平均SPI各分量方差貢獻率

功率譜的分析結(jié)果顯示(圖3)，鹽池縣年平均SPI原始序列主要存在一個4 a的準周期；同心存在5.78、3.25 a和2.36 a的準周期，其中5.78 a的準周期最強；而中寧縣則表現(xiàn)出6個不同的準周期特征，但最強的為13、4.33 a和2.17 a的準周期。對原始年平均SPI序列方差貢獻率在67.35%～74.93%的IMF1分量和方差貢獻率在14.03%～18.64%的IMF2分量分析的結(jié)果顯示，鹽池、中寧和同心3縣年平均SPI的IMF1分量存在3.25～4 a的準周期，而IMF2分量存在5.78～6.5 a的準周期。這一結(jié)果說明，在寧夏中部干旱帶各縣的平均干旱強度年際波動周期特征中，有平均超過65%的貢獻來自于3.25～4 a左右的干旱周期波動，而平均超過14%的貢獻來自于5.78～6.5 a左右的干旱周期波動。

圖3 中部干旱帶各縣年平均SPI序列功率譜分析結(jié)果

Fig.3 Power spectrum analysis results of annual SPI in arid zone of middle Ningxia

2.4 干旱趨勢分析

為了探討寧夏中部干旱帶近50年來干旱的變化趨勢特征，本研究利用一元線性回歸對鹽池、中寧和同心3縣全年SPI平均值的變化趨勢進行了分析，同時也分析了各站點春季(3—5月)、夏季(6—8月)、秋季(9—11月)和冬季(12月到次年2月)SPI平均值的變化趨勢，以此來探討寧夏中部干旱帶不同地區(qū)、不同季節(jié)干旱演變特征，結(jié)果如表4所示。從全年SPI平均值的變化趨勢來看，3縣的SPI值在近50年中均呈現(xiàn)降低趨勢，而SPI值越低則表征該縣氣象干旱越嚴重，由此可知在近50年中，寧夏中部干旱帶的氣象干旱一直處于增強的趨勢，且以同心縣為最強，一元線性回歸斜率為-0.0050。但不同季節(jié)的干旱趨勢特征又表現(xiàn)出差異，其中3縣的春旱、秋旱表現(xiàn)為增強趨勢；鹽池和中寧的冬旱有所減弱，但同心的冬旱在近50年基本不變；鹽池、同心的夏旱在減弱，但中寧的夏旱卻在增強。但近50年的干旱線性趨勢中，只有同心的春旱增強趨勢通過P<0.05的顯著性檢驗，秋旱增強趨勢通過P<0.01的極顯著性檢驗，其他趨勢均未通過顯著性檢驗。

同時，研究計算了近50年鹽池、中寧和同心3縣全年及四個季節(jié)SPI均值的赫斯特指數(shù)，用于分析寧夏中部干旱帶干旱趨勢的長程依賴性和持續(xù)性，并依據(jù)赫斯特值的大小，劃分出較弱持續(xù)(0.50.7)和較弱逆轉(zhuǎn)(0.4≤H<0.5)等不同的持續(xù)特征。從表4可以看出，3縣年干旱強度均具有持續(xù)的增強趨勢，且同心達到了中等持續(xù)；從季節(jié)來看，鹽池和中寧2縣春旱增強的趨勢在未來會出現(xiàn)較弱逆轉(zhuǎn)，而同心春旱增強的趨勢則繼續(xù)保持較弱的持續(xù)性；夏、秋和冬季的赫斯特指數(shù)均大于0.5，即這三個季節(jié)鹽池、中寧和同心3縣的干旱會繼續(xù)保持當前50年的趨勢，并表現(xiàn)出不同的持續(xù)性。

表4 近50年干旱變化趨勢及持續(xù)特征

注：*和**分別表示P<0.05、P<0.01顯著水平。

Note: * and ** indicate significance atP<0.05 andP<0.01, respectively.

3 結(jié) 論

本研究利用寧夏中部干旱帶鹽池、中寧和同心3縣1960—2012年的氣象資料，在計算標準化降水指數(shù)(SPI)的基礎(chǔ)上，利用集合經(jīng)驗?zāi)B(tài)分解、功率譜分析和重新標度極差分析等方法，對中部干旱帶近50年來的干旱變化特征進行了研究，研究取得如下幾點結(jié)論：

1) 不同干旱等級的統(tǒng)計結(jié)果表明，在1960—2012年間，中部干旱帶鹽池、中寧和同心3縣發(fā)生的干旱事件中，超過中旱以上的干旱事件達到了總干旱次數(shù)的30%以上；其中中旱平均0.93—1.2 a發(fā)生一次，為年年有中旱，特大干旱8.83—17.67 a發(fā)生一次，為十年一大旱。

2) 集合經(jīng)驗?zāi)B(tài)分解獲取的干旱災(zāi)害多尺度變化特征表明，寧夏中部干旱帶鹽池、中寧和同心3縣的干旱災(zāi)害存在著0.5 a尺度的輕旱，1 a尺度的中旱，3 a尺度的重旱和10 a尺度的特旱。

3) 功率譜分析得出，在寧夏中部干旱帶各縣的平均干旱強度年際波動周期特征中，有平均超過65%的貢獻來自于3.25-4 a左右的干旱周期波動，而平均超過14%的貢獻來自于5.78-6.5 a左右的干旱周期波動。

4) 寧夏中部干旱帶的干旱趨勢總體處于增強態(tài)勢，但不同季節(jié)的干旱趨勢特征又表現(xiàn)出差異，鹽池、同心的夏旱在減弱，但中寧的夏旱卻在增強；鹽池和中寧的冬旱在減弱，但同心的冬旱基本不變。赫斯特指數(shù)分析顯示，中部干旱帶各縣年平均干旱強度的增強趨勢均具有長程依賴性，即未來還將持續(xù)干旱。

[1] 徐建文，居 輝，劉 勤,等.黃淮海地區(qū)干旱變化特征及其對氣候變化的響應(yīng)[J].生態(tài)學(xué)報,2014,34(2):460-470.

[2] 陳曉光,李劍萍,韓穎娟,等.寧夏近20年來植被覆蓋度及其與氣溫降水的關(guān)系[J].生態(tài)學(xué)雜志,2007,26(9):1375-1383.

[3] 馬柱國,華麗娟,任小波.中國近代北方極端干濕事件的演變規(guī)律[J].地理學(xué)報,2003,58(s):69-74.

[4] 白 樺,穆興民,王 飛,等.渭河流域氣象及水文旱澇演變規(guī)律分析[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(2):237-241.

[5] 黃妙芬.黃土高原西北部地區(qū)的旱度模式[J].氣象,1990,17(1):23-28.

[6] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)[M].北京:氣象出版社,2007.

[7] Wu Z H, Schneider E K, Kirtman B P, et al. The modulated annual cycle:an alternative reference frame for climate anomalies[J]. Climate Dynamics, 2008,31(7-8):823-841.

[8] Breaker L C, Ruzmaikin A. The 154-year record of sea level at San Francisco: extracting the long-term trend, recent changes, and other tidbits[J]. Climate Dynamics, 2011,36(3-4):545-559.

[9] 江田漢,鄧蓮堂.Hurst指數(shù)估計中存在的若干問題——以在氣候變化研究中的應(yīng)用為例[J].地理科學(xué),2004,24(2):177-182.

[10] 信忠保,謝志仁.寧夏氣候變化對ENSO事件的響應(yīng)[J].干旱區(qū)地理,2005,28(2):239-243.

[11] 黃晚華,楊曉光,李茂松,等.基于標準化降水指數(shù)的中國南方季節(jié)性干旱近58 a演變特征[J].農(nóng)業(yè)工程學(xué)報,2010,26(7):50-59.

[12] 杜靈通,田慶久,黃 彥,等.基于TRMM數(shù)據(jù)的山東省干旱監(jiān)測及其可靠性檢驗[J].農(nóng)業(yè)工程學(xué)報,2012,28(2):121-126.

[13] 龔艷冰,張繼國,劉高峰,等.基于SPI指數(shù)與R/S分析的曲靖市干旱特征研究[J].干旱地區(qū)農(nóng)業(yè)研究,2015,33(3):273-277.

[14] Huang N E, Wu Z H. A review on Hilbert-Huang transform: method and its applications to geophysical studies[J]. Reviews of Geophysics, 2008,46(2):2006.

[15] 孫銀鳳,陸寶宏.基于EEMD的南京市降水特征分析[J].中國農(nóng)村水利水電,2013,(3):5-9.

[16] 薛春芳,侯 威,趙俊虎,等.集合經(jīng)驗?zāi)B(tài)分解在區(qū)域降水變化多尺度分析及氣候變化響應(yīng)研究中的應(yīng)用[J].物理學(xué)報,2013,62(10):504-511.

[17] 李艷萍,陳昌春,張余慶,等.基于EEMD的明代北京地區(qū)干旱災(zāi)害特征分析[J].中國沙漠,2014,(3):835-840.

[18] 孫 陽,陳元芳,程 龍,等.基于EEMD的枯季入庫徑流預(yù)報分析[J].中國農(nóng)村水利水電,2012,(2):34-37.

[19] 黃嘉佑.氣象中的譜分析[M].北京:氣象出版社,1984.

[20] Priestley M B. Spectral Analysis and Time Series[M]. New York: Academic Press Inc, 1981.

[21] 黃忠恕.波譜分析方法及其在水文氣象中的應(yīng)用[M].北京:氣象出版社,1983.

[22] 雷兆崇.數(shù)值模式中的譜方法[M].北京:氣象出版社,1991.

[23] 吳月祥,董立紅,王 丹.基于功率譜的扎龍濕地降水周期分析[J].齊齊哈爾大學(xué)學(xué)報,2014,30(3):79-82.

[24] 覃邑龍,應(yīng)益榮.重標極差方法下時變Hurst指數(shù)的構(gòu)建和實證研究[J].系統(tǒng)管理學(xué)報,2011,20(5):620-626.

[25] 張 羽中,任志遠.黃土高原地區(qū)植被覆蓋變化的時空差異及未來趨勢[J].資源科學(xué),2011,33(11):2143-2149.

[26] 杜靈通,候 靜,胡 悅,等.基于遙感溫度植被干旱指數(shù)的寧夏2000—2010年旱情變化特征[J].農(nóng)業(yè)工程學(xué)報,2015,31(14):209-216.

Drought characteristics in arid zone of middle Ningxia from 1960 to 2012 base on SPI index

HU Yue1,2, DU Ling-tong1,2, HOU Jing1,2, LIU Ke1,2, ZHU Yu-guo1,2

(1.BreedingBaseforStateKeyLaboratoryofLandDegradationandEcologicalRestorationinNorthwestChina,NingxiaUniversity,Yinchuan,Ningxia750021,China; 2.KeyLaboratoryforRestorationandReconstructionofDegradedEcosysteminNorthwestChinaofMinistryofEducation,NingxiaUniversity,Yinchuan,Ningxia750021,China)

With the purpose to provide scientific basis for regional agricultural development and disaster prevention and mitigation, by using the meteorological data from three in-situ stations during 1960 to 2012, the drought characteristics and variation trends in arid zone of middle Ningxia were studied by the method of standardized precipitation index (SPI), along with ensemble empirical mode decomposition (EEMD), power spectrum analysis and rescaled range analysis (R/S). The results showed that more than 30% of drought events, which happened in arid zone of middle Ningxia from 1960 to 2012, were moderate or even more serious drought. The multi-scale decomposition of drought period showed that there existed four time scales for the change cycle of different intensities of drought, i.e., 0.5 a for mild drought, 1 a for moderate drought, 3 a for severe drought, and 10 a for extreme drought. The EEMD analysis showed that the inter-annual fluctuation of drought was mainly derived from a 3.25-4 a period which contributed more than 65% to the whole variance, and followed by a 5.78-6.5 a period which contributed more than 14%. On the whole, the drought in arid zone of middle Ningxia is enhancing currently and this trend has a long-range dependence. Therefore, this region will keep a sustained drought in the future, but the characteristic of drought trend in various seasons is different.

drought; standardized precipitation index; ensemble empirical mode decomposition; power spectrum analysis; arid zone of middle Ningxia

1000-7601(2017)02-0255-08

10.7606/j.issn.1000-7601.2017.02.41

2015-12-24基金項目：國家自然科學(xué)基金(41201438)；寧夏高等學(xué)?？茖W(xué)研究重點項目(NGY2014001)；寧夏大學(xué)研究生創(chuàng)新項目(GIP2015041)

胡 悅(1992—)，女，寧夏銀川人，碩士研究生，主要從事遙感災(zāi)害研究。 E-mail: 365493815@qq.com。

杜靈通(1980—)，男，寧夏同心人，博士，副研究員，主要從事資源環(huán)境遙感研究。 E-mail:dult80@qq.com。

S162.1

A