基于SPI的旱澇氣候演變及對洪湖水位的影響分析

鄧艷君 周守華 楊超

摘要： 利用洪湖市氣象站逐月降水?dāng)?shù)據(jù)和洪湖水位站水位資料，計(jì)算不同時(shí)間尺度標(biāo)準(zhǔn)化降水指數(shù)（SPI）及其干旱等級，分析洪湖市旱澇氣候演變及不同時(shí)間尺度旱澇氣候事件演變特征，探討降水量、不同時(shí)間尺度SPI與洪湖水位的相關(guān)性。結(jié)果表明，洪湖市近62年來干旱年和洪澇年頻次和等級均相差較小，不同季節(jié)旱澇變化及等級頻次分布存在一定差異，洪湖市呈濕潤化演變趨勢，SPI以0.086/10年的速率弱增長，春、秋季SPI呈下降趨勢，干旱略加劇，夏、冬季SPI呈增長趨勢，洪澇加劇;洪湖水位與降水變化趨勢一致，但時(shí)間上有所滯后，洪湖水位與SPI具有較好的相關(guān)性，相關(guān)程度因SPI的時(shí)間尺度存在一定差異，其中SPI3與水位相關(guān)性最好，利用3個(gè)月時(shí)間尺度的SPI3來監(jiān)測洪湖水位變化，兩者具有很好的一致性，可以采用SPI3指數(shù)監(jiān)測洪湖水位。

關(guān)鍵詞：旱澇; SPI; 洪湖; 水位

中圖分類號：P467;P426.616 ? ?文獻(xiàn)標(biāo)識碼：A

文章編號：0439-8114（2020）10-0032-06

DOI：10.14088/j.cnki.issn0439-8114.2020.10.007 ? ? ? ? ? ? 開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Abstract： Based on the monthly precipitation data from Honghu meteorological station and water level data from Honghu lake， variables time-scale standard precipitation index（SPI） was calculated. Using SPI， the changes of drought-flood and the temporal characteristics of drought-flood events were analyzed in Honghu lake. The correlation among precipitation， various time-scale SPI and water lever were investigated. The results showed that the frequency and grade of the drought and flood events in Honghu lake in the past 62 years were approximate on the annual scale， little difference on seasonal scale. There was a wetting trend in Honghu， and annual SPI decreased with an average rate of 0.086/10a. SPI declined in spring and autumn， and slightly intensified in drought; SPI increased in summer and winter， and increased in flood and waterlogging. The variability between the water level and the precipitation data was consistent， but the water level later than the precipitation data three months. Variables time-scale SPIand water level were significant correlation. The 3 month scale SPI3 was most closely correlated with water level of Honghu lake. The 3 month scale SPI3 agrees well with the variability of the water level. SPI3 can be used to monitor the water level of Honghu lake.

Key words： drought and flood; standard precipitation index（SPI）; Honghu lake; water level

在近百年以全球變暖為主要特征的氣候變化影響下，洪湖濕地氣候與生態(tài)出現(xiàn)了顯著的變化[1]，極端旱澇氣候事件對湖泊濕地造成嚴(yán)重威脅、破壞[2]，2010―2011年旱澇急轉(zhuǎn)、2016年全流域洪澇，使洪湖濕地及其生態(tài)系統(tǒng)遭受極大危害，已經(jīng)引起當(dāng)?shù)卣坝嘘P(guān)學(xué)者的關(guān)注。國內(nèi)外對極端天氣氣候事件進(jìn)行了許多分析研究[3-5]，但主要關(guān)注的是極端氣候事件的時(shí)間空間變化特點(diǎn)。對于極端氣候事件的影響分析中，對農(nóng)業(yè)的影響已有了相關(guān)研究[6-11]，針對極端氣候事件對湖泊濕地等自然生態(tài)系統(tǒng)的影響研究相對較少。胡毅鴻等[12]利用標(biāo)準(zhǔn)化降水指數(shù)和夏季長周期旱澇急轉(zhuǎn)指數(shù)，分析洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征;劉敏等[13]分析不同時(shí)間尺度SPEI/SPI干旱指數(shù)與洪湖水位的相關(guān)性，發(fā)現(xiàn)基于5個(gè)月時(shí)間尺度的SPEI/SPI干旱指數(shù)監(jiān)測歷史時(shí)期洪湖水位變化具有很好的一致性;程玉菲等[14]研究疏勒河流域極端水文事件對極端氣候事件的響應(yīng)關(guān)系，發(fā)現(xiàn)極端洪水事件主要受控于極端降水事件，特別是極端降水總量;黃菊梅等[15]基于標(biāo)準(zhǔn)化降水指數(shù)分析了洞庭湖區(qū)干旱時(shí)空分布特征，計(jì)算了標(biāo)準(zhǔn)化降水指數(shù)、降水量和水位之間的相關(guān)性。本研究以洪湖為例，計(jì)算不同尺度標(biāo)準(zhǔn)化降水指數(shù)（SPI）及其干旱等級，基于SPI分析洪湖市旱澇氣候演變特征及不同時(shí)間尺度旱澇氣候事件變化特征，探討SPI、降水量和水位之間的相關(guān)性，分析降水對洪湖水位的影響。這對把握洪湖市旱澇出現(xiàn)規(guī)律、指導(dǎo)洪湖科學(xué)調(diào)度、防汛抗旱、防災(zāi)減災(zāi)具有參考意義。

1 ?研究區(qū)概況

洪湖市位于湖北省江漢平原腹地，以境內(nèi)最大的湖泊——洪湖而命名，屬亞熱帶濕潤季風(fēng)氣候，其特點(diǎn)是冬夏長，春秋短，四季分明，光照充足，雨量充沛，年均降水量1 174 mm，降水主要集中于春夏，最大月降水量出現(xiàn)在6月，最小月降水量出現(xiàn)在12月（圖1）。洪湖位于四湖流域下游入口，是湖北省第一大淡水湖，年平均水位24.6 m，高水位出現(xiàn)在7―9月，低水位出現(xiàn)在12月至次年4月，滯后于月降水量2～3個(gè)月。

2 ?材料與方法

2.1 ?資料來源

氣候資料采用洪湖氣象站1957年1月至2018年12月逐月平均降水量。洪湖水位資料采用洪湖湖區(qū)水位代表站挖溝咀水位站2000年1月至2018年12月月平均觀測資料。為消除年周期及長期趨勢，對原始水位資料進(jìn)行了標(biāo)準(zhǔn)化預(yù)處理。

2.2 ?研究方法

采用具有多時(shí)間尺度特征的標(biāo)準(zhǔn)化降水指數(shù)作為衡量洪湖市氣候旱澇指標(biāo)。SPI計(jì)算簡單，僅需降水資料即可計(jì)算，具有良好的穩(wěn)定性，適用于不同地區(qū)和不同時(shí)間尺度的旱澇監(jiān)測評估服務(wù)，被廣泛應(yīng)用于旱澇評估中[16-20]。標(biāo)準(zhǔn)化降水指數(shù)假定降水量的變化符合Γ分布，運(yùn)用數(shù)理方法，對其經(jīng)過正態(tài)標(biāo)準(zhǔn)化處理得到指數(shù)。具體計(jì)算方法參考?xì)庀蟾珊档燃塆B/T 20481―2017[21]。本研究利用1957―2018年近62年的降水量資料，計(jì)算了SPI1、SPI3、SPI6、SPI12（分別為1、3、6和12個(gè)月尺度的SPI）。季節(jié)SPI取每季最后1個(gè)月的SPI3，即取5月、8月、11月、次年2月的SPI3，分別代表春、夏、秋、冬的SPI。年SPI取每年12月的SPI12，代表該年的SPI。

根據(jù)上述方法計(jì)算月、季、年尺度的標(biāo)準(zhǔn)化降水指數(shù)，并參考GB/T 20481―2017[21]和相關(guān)文獻(xiàn)[16-18]中分級標(biāo)準(zhǔn)劃分旱澇等級，如表1所示。

Pearson相關(guān)系數(shù)是反映變量之間相關(guān)關(guān)系密切程度的統(tǒng)計(jì)指標(biāo)，通過兩個(gè)離差乘積反映兩個(gè)變量之間的相關(guān)程度，著重研究線性的單相關(guān)系數(shù)。本研究利用統(tǒng)計(jì)學(xué)中常用的Pearson相關(guān)系數(shù)來衡量降水資料與水位資料的相關(guān)性，并采用t檢驗(yàn)進(jìn)行相關(guān)系數(shù)顯著性檢驗(yàn)，其中月尺度相關(guān)分析樣本量為228。

3 ?結(jié)果與分析

3.1 ?不同時(shí)間尺度SPI特征

干旱與洪澇的發(fā)生和發(fā)展往往表現(xiàn)在不同的時(shí)間尺度上，SPI1正負(fù)波動頻繁，能夠反映短期降水不足，SPI3反映季節(jié)旱澇變化特征、SPI6和SPI12旱澇波動較平緩，對長時(shí)間旱澇有較好的反映，這個(gè)特點(diǎn)在SPI的時(shí)間序列上具有明顯的體現(xiàn)。洪湖氣象站1957―2018年1、3、6、12個(gè)月4個(gè)時(shí)間尺度的SPI時(shí)間變化序列見圖2。由圖2可以看出，不同時(shí)間尺度SPI變化趨勢有顯著的一致性，但在發(fā)生頻率和持續(xù)時(shí)間上有明顯差異，SPI1波動頻繁，干旱、洪澇事件發(fā)生頻率高，持續(xù)時(shí)間短，較長時(shí)間尺度的SPI3、SPI6和SPI12干旱、洪澇事件發(fā)生頻率越來越低，持續(xù)時(shí)間越來越長，多種時(shí)間尺度SPI相結(jié)合可以實(shí)現(xiàn)對旱澇的綜合檢測。以2010―2011年旱澇急轉(zhuǎn)的典型年份為例，SPI顯示2010年1―2月為輕旱，3―4月為輕澇到重澇，5月正常，6月輕旱，7月極澇，8月正常，9―10月中澇，11月中旱，12月正常，2011年1―5月輕旱到極旱，6月重澇，7月中旱，8―11月正常，12月輕旱;SPI3顯示2010年1月正常，2月輕旱，3月輕澇，4―5月中澇，6月正常，7―11月中澇到重澇，12月輕澇，2011年1―5月重旱到極旱，6―8月正常，9―12月輕旱到中旱;SPI6顯示2010年1―2月重旱到極旱，3月正常，4―5月輕澇，6月正常，7―12月重澇到極澇，2011年1―2月正常，3―9月輕旱到極旱，10―11月正常，12月重旱;SPI12顯示2010年1―6月正常到輕旱，2010年7月到2011年6月為輕澇到極澇，2011年7―12月為輕旱到中旱?？梢婋S著時(shí)間尺度的增加，旱澇起止時(shí)間相應(yīng)延后，旱澇等級發(fā)生變化，充分反映了前期降水變化的累積影響，這與學(xué)者們基于SPI研究旱澇特征的結(jié)論[16-20]基本一致。

3.2 ?洪湖市旱澇變化特征

3.2.1 旱澇年際變化特征 ? ?年SPI可以監(jiān)測旱澇年份及等級，由圖3可知，洪湖市近62年來年SPI呈弱增長趨勢，速率為0.086/10年，干旱年和洪澇年分別占21年和22年，頻次和等級均相差較?。▓D4）。極旱年發(fā)生在1968、1963和1984年，為重旱年;極澇年發(fā)生在1996、2002和2010年，為重澇年。

3.2.2 旱澇季節(jié)變化特征 ? ?洪湖市不同季節(jié)旱澇變化（圖5）及等級頻次（圖6）的分布存在較大差異，春、夏、秋、冬季速率分別為-0.022/10年、0.095/10年、-0.011/10年和0.093/10年，春、秋季SPI呈降低趨勢，干旱略加劇，夏、冬季SPI呈增長趨勢，洪澇加劇，且夏、冬季對年濕潤化的貢獻(xiàn)明顯，年SPI也呈增長趨勢。中度及以上旱澇頻率春、夏、秋、冬季分別為27.4%、37.1%、33.9%、30.6%。極旱季出現(xiàn)在1997年春、2011年春、1979年秋、2007年秋、1978年冬、1998年冬;極澇季出現(xiàn)在1973年春、2002年春、1996年夏、1972年秋、1958年冬、2002年冬。

由表2可以看出，洪湖市1968年極旱年由春、夏、秋、冬四季持續(xù)中度到重度干旱所致;1963年和1984年重旱年主要由夏季中旱和冬季重旱所致;2002年重澇年由春、冬季極澇導(dǎo)致;2010年春、夏、秋三季連續(xù)中度洪澇，雖冬季轉(zhuǎn)為重旱，年旱澇等級仍為重澇;春季微旱、秋冬季正常、夏季極澇導(dǎo)致1996年洪澇程度最為嚴(yán)重。

3.3 ?洪湖水位變化特征及與降水、SPI的關(guān)系

2010―2018年近9年洪湖水位時(shí)間序列變化如圖7所示，2011年上半年水位逐漸下降，是對2010年冬季重旱和2011年春季極旱的滯后響應(yīng);2010年7―8月持續(xù)高水位是由2010年春、夏、秋三季連續(xù)中澇導(dǎo)致;2016年洪湖春季微澇、夏季中澇，但受超強(qiáng)厄爾尼諾事件和拉尼娜事件的先后影響，長江流域汛期時(shí)間長、降雨過程多，受本地降水和長江流域降水的共同影響，導(dǎo)致7―8月洪湖持續(xù)高水位，7月3日就超過警戒水位的26.2 m，達(dá)26.9 m。

利用2000―2018年逐月228個(gè)樣本資料進(jìn)行降水量、SPI和水位的相關(guān)分析，從表3可以看出，月降水量距平與SPI的相關(guān)性最好，高達(dá)0.93，降水量與距平和SPI1相關(guān)系數(shù)分別為0.69和0.67，距平與SPI3的相關(guān)系數(shù)為0.53，SPI1與SPI3的相關(guān)系數(shù)為0.55，SPI3與SPI6的相關(guān)系數(shù)為0.71，SPI6與SPI12的相關(guān)系數(shù)為0.61，且均通過0.01的顯著性檢驗(yàn)，因此多時(shí)間尺度SPI與降水顯著相關(guān)。SPI3與洪湖水位的相關(guān)系數(shù)為0.51，相關(guān)性最好，其次是SPI6，相關(guān)系數(shù)為0.41，且均通過0.01的顯著性檢驗(yàn)。

依據(jù)上述分析結(jié)果，采用3個(gè)月時(shí)間尺度SPI3對歷史洪湖水位進(jìn)行檢測，結(jié)果見圖8。由圖8可以看出，SPI3指數(shù)與洪湖水位總體波動基本一致，個(gè)別年份和季節(jié)的SPI3和標(biāo)準(zhǔn)化水位距平程度略有差異。2000―2018年19年228個(gè)樣本的歷史序列中，水位距平相比SPI3偏高的時(shí)間段有2000年秋季、2013年冬至2014年春、2018年下半年，標(biāo)準(zhǔn)化水位距平相比SPI3指數(shù)偏低的時(shí)間段有2003年冬至2004年春、2006年秋、2008年冬至2009年春。水位監(jiān)測效果偏差的時(shí)段大部分出現(xiàn)在非汛期，考慮大規(guī)模水利工程的建設(shè)，洪湖水位除受自然降水等氣候因子變化的影響外，人工調(diào)控、人類活動的影響也非常顯著[22-25]，與劉敏等[13]、劉可群等[25]結(jié)論一致。因此，利用3個(gè)月時(shí)間尺度SPI3監(jiān)測洪湖水位，基本可以反映旱澇氣候變化對洪湖水位的影響。

4 ?小結(jié)與討論

具有多時(shí)間尺度的標(biāo)準(zhǔn)化降水指數(shù)使用數(shù)據(jù)少，計(jì)算簡便，普適性強(qiáng)，廣泛應(yīng)用于旱澇氣候事件監(jiān)測。利用多時(shí)間尺度SPI，分析洪湖市近62年旱澇氣候事件變化規(guī)律，并探討降水量、SPI和水位之間的相關(guān)關(guān)系，選取適當(dāng)時(shí)間尺度SPI監(jiān)測洪湖水位，得到了如下結(jié)論。

1）洪湖降水主要集中在春、夏季，最大月降水量出現(xiàn)在6月，最小月降水量出現(xiàn)在12月。洪湖高水位出現(xiàn)在7―9月，低水位出現(xiàn)在12月至次年4月，滯后于降水量2～3個(gè)月。

2）SPI1正負(fù)波動頻繁，能夠反映短期降水不足;SPI3反映季節(jié)旱澇變化特征;SPI6和SPI12旱澇波動較平緩，對長時(shí)間旱澇有較好地反映，且隨著時(shí)間尺度的增加，旱澇起止時(shí)間相應(yīng)延后，旱澇等級發(fā)生變化，充分反映了前期降水變化的累積影響。

3）洪湖市近62年來不同季節(jié)旱澇變化及等級頻次的分布存在較大差異，春、秋季SPI呈降低趨勢，干旱略加劇，夏、冬季SPI呈增長趨勢，洪澇加劇，且夏、冬季對年濕潤化的貢獻(xiàn)明顯，年SPI呈弱增長趨勢，干旱年和洪澇年頻次和等級均相差較小。

4）2010―2018年洪湖水位的時(shí)間序列可以反映出水位對降水的滯后響應(yīng)，特別是2010―2011年旱澇急轉(zhuǎn)、2016年洪澇等旱澇氣候事件對洪湖水位的影響顯著。

5）具有多時(shí)間尺度的SPI與降水量顯著相關(guān)。SPI3與洪湖水位的相關(guān)系數(shù)為0.51，相關(guān)性最好，其次是SPI6，相關(guān)系數(shù)為0.41，且均通過0.01的顯著性檢驗(yàn)。采用3個(gè)月時(shí)間尺度SPI3對歷史洪湖水位進(jìn)行檢測，可以看出SPI3與洪湖水位總體波動基本一致，個(gè)別年份和季節(jié)的SPI3與水位波動程度略有差異，SPI3可以反映洪湖水位變化狀況，可用來監(jiān)測洪湖水位。

參考文獻(xiàn)：

[1] 王慧亮，王學(xué)雷，厲恩華.氣候變化對洪湖濕地的影響[J].長江流域資源與環(huán)境，2010，19（6）：653-658.

[2] 吳秀芹，龍花樓，高吉喜，等.江漢平原濕地功能下降與洪澇災(zāi)害關(guān)系分析[J].生態(tài)環(huán)境，2005（6）：884-889.

[3] IPCC. Special report on global warming of 1.5℃[M]. UK：Cambridge university press，2018.

[4] IPCC. Climate Change 2014： Impacts，adaptation， and vulnerability [M/OL].Cambridge：Cambridge university press，2014.[2014-09-03].

[5] 蔡文香.中國極端氣候事件的趨勢特征與極值分布[D].北京：對外經(jīng)濟(jì)貿(mào)易大學(xué)，2016.

[6] 柏會子，肖登攀，劉劍鋒，等.1965―2014年華北地區(qū)極端氣候事件與農(nóng)業(yè)氣象災(zāi)害時(shí)空格局研究[J].地理與地理信息科學(xué)，2018，34（5）：99-105.

[7] 劉勝利.氣候變化對我國雙季稻區(qū)水稻生產(chǎn)的影響與技術(shù)適應(yīng)研究[D].北京：中國農(nóng)業(yè)大學(xué)，2018.

[8] 白美蘭，郝潤全，高建峰，等.內(nèi)蒙古地區(qū)極端氣候事件分布特征及對農(nóng)業(yè)影響評估[J].干旱地區(qū)農(nóng)業(yè)研究，2009，27（2）：21-27.

[9] 付 ?偉.安徽省極端氣候事件對農(nóng)業(yè)經(jīng)濟(jì)的影響分析[D].安徽蕪湖：安徽師范大學(xué)，2013.

[10] 王明田.氣候變化背景下四川農(nóng)業(yè)季節(jié)性干旱的發(fā)展趨勢及應(yīng)對措施[D].成都：四川農(nóng)業(yè)大學(xué)，2012.

[11] 馮相昭，鄒 ?驥，馬 ?珊，等.極端氣候事件對中國農(nóng)村經(jīng)濟(jì)影響的評價(jià)[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2007（2）：19-25.

[12] 胡毅鴻，李景保.1951—2015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2017，33（7）：107-115.

[13] 劉 ?敏，秦鵬程，劉可群，等.洪湖水位對不同時(shí)間尺度SPEI/SPI干旱指數(shù)的響應(yīng)研究[J].氣象，2013，39（9）：1163-1170.

[14] 程玉菲，程文舉，胡想全，等.疏勒河流域極端水文事件對極端氣候的響應(yīng)[J].高原氣象，2019，38（3）：583-592.

[15] 黃菊梅，徐志偉，蔡海朝，等.基于標(biāo)準(zhǔn)化降水指數(shù)的洞庭湖區(qū)干旱特征[A].S5全球典型干旱半干旱地區(qū)氣候變化及其影響[C].中國氣象學(xué)會，2012.8.

[16] 張苗苗.基于SPI的近51a晉北地區(qū)旱澇變化及干旱事件時(shí)空特征研究[J].干旱區(qū)資源與環(huán)境，2018，32（3）：138-144.

[17] 劉佳旭，李九一，李麗娟，等.基于降水?dāng)?shù)據(jù)的云南省近61年旱澇特性研究[J].熱帶氣象學(xué)報(bào)，2018，34（1）：68-77.

[18] 楊國慶，王佳真，孫萌萌.基于標(biāo)準(zhǔn)化降水指數(shù)的滄州干旱時(shí)空特征[J].干旱氣象，2019，37（2）：218-225，242.

[19] 李虹雨，馬 ?龍，劉廷璽，等.基于標(biāo)準(zhǔn)化降水指數(shù)的內(nèi)蒙古地區(qū)干旱時(shí)空變化特征[J].水文，2018，38（5）：47-51，90.

[20] 車少靜，李春強(qiáng).基于標(biāo)準(zhǔn)化降水指數(shù)的石家莊干旱時(shí)空特征[J].氣象科技，2010，38（1）：66-70.

[21] GB/T 20481―2017，氣象干旱等級[S].

[22] 王學(xué)雷，寧龍梅，肖 ?銳.洪湖濕地恢復(fù)中的生態(tài)水位控制與江湖聯(lián)系研究[J].濕地科學(xué)，2008，6（2）：316-320.

[23] 金衛(wèi)斌，胡秉民.湖北四湖流域景觀結(jié)構(gòu)變化對汛期湖泊水位影響的模擬分析[J].生態(tài)學(xué)報(bào)，2003，23（4）：642-648.

[24] 鄭紹偉，慕長龍，陳祖銘，等.長江上游森林影響流域水文過程模擬分析[J].生態(tài)學(xué)報(bào)，2010，30（11）：3046-3056.

[25] 劉可群，梁益同，周金蓮，等.人類活動與氣候變化對洪湖春旱的影響[J].生態(tài)學(xué)報(bào)，2014，34（5）：1302-1310.