近45年來瀏陽(yáng)河流域極端降水變化

于澤興, 胡國(guó)華, 陳 肖, 顧慶福, 邵全喜,3, 宋荷花

(1.長(zhǎng)沙理工大學(xué) 水利工程學(xué)院, 長(zhǎng)沙 410114; 2.湖南省水文水資源勘測(cè)局,長(zhǎng)沙 410007; 3.澳大利亞國(guó)家數(shù)據(jù)研究所, 佛意特 WA6014; 4.長(zhǎng)沙水文水資源勘測(cè)局, 長(zhǎng)沙 410001)

近45年來瀏陽(yáng)河流域極端降水變化

于澤興1, 胡國(guó)華1, 陳 肖1, 顧慶福2, 邵全喜1,3, 宋荷花4

(1.長(zhǎng)沙理工大學(xué) 水利工程學(xué)院, 長(zhǎng)沙 410114; 2.湖南省水文水資源勘測(cè)局,長(zhǎng)沙 410007; 3.澳大利亞國(guó)家數(shù)據(jù)研究所, 佛意特 WA6014; 4.長(zhǎng)沙水文水資源勘測(cè)局, 長(zhǎng)沙 410001)

利用1969—2013年瀏陽(yáng)河流域12個(gè)雨量站的降水資料，結(jié)合百分位方法定義極端降水閾值，采用線性回歸、滑動(dòng)平均、Mann-Kendall檢驗(yàn)、累積距平聯(lián)合滑動(dòng)t檢驗(yàn)和Morlet小波分析等方法分析了極端降水指數(shù)的演變規(guī)律。結(jié)果表明：(1) 極端降水量、極端降水強(qiáng)度、極端降水比率表現(xiàn)為上升趨勢(shì)，10年上升率分別為7.3 mm，9.5 mm/d，0.3%，極端降水日數(shù)表現(xiàn)為下降趨勢(shì)，下降率為-0.054 d/10 a。(2) 極端降水量、極端降水強(qiáng)度分別在1991年、1992年存在明顯突變，極端降水日數(shù)、極端降水比率突變不顯著。(3) 極端降水指數(shù)均存在7～8 a,10～11 a和25～27 a的主周期，且豐枯交替變化明顯。研究結(jié)果對(duì)于指導(dǎo)瀏陽(yáng)河流域的防洪減災(zāi)具有重要意義。

極端降水; Mann-Kendall檢驗(yàn); Morlet小波分析; 累積距平聯(lián)合滑動(dòng)t檢驗(yàn); 瀏陽(yáng)河流域

受全球變暖和水循環(huán)加劇的影響，區(qū)域氣候也在發(fā)生著變化，極端降水事件的頻率增加，干旱、洪澇災(zāi)害等極值事件頻發(fā)，對(duì)區(qū)域的社會(huì)經(jīng)濟(jì)和生態(tài)環(huán)境造成了嚴(yán)重影響。分析極端降水事件的演變規(guī)律有利于科學(xué)理解洪旱災(zāi)害的時(shí)空分布，為減少災(zāi)害損失以及水資源管理提供科學(xué)依據(jù)[1-4]。目前針對(duì)不同的地理區(qū)位和氣候環(huán)境，對(duì)中尺度流域極端降水事件的研究取得了許多成果，如對(duì)我國(guó)的淮河流域、東江流域、紅河流域、黑河流域、涇河流域和湘江流域等區(qū)域極端降水事件的時(shí)空變化進(jìn)行了相關(guān)研究[5-10]。極端降水事件在不同的時(shí)空尺度都有其各自的特點(diǎn)，大尺度區(qū)域的氣候變化特點(diǎn)和趨勢(shì)變化具有較高的相似性和一致性，而對(duì)于中、小區(qū)域尺度上的降水變化規(guī)律卻存在著較大的差異和相反的變化特征，由于不同區(qū)域在地形、氣候環(huán)境方面存在著差異，具體的區(qū)域仍急需開展相關(guān)的研究和評(píng)估工作。

瀏陽(yáng)河流域位于湖南省東北部，途徑瀏陽(yáng)市、長(zhǎng)沙縣及長(zhǎng)沙市區(qū)。該區(qū)域經(jīng)濟(jì)發(fā)達(dá)，2015年瀏陽(yáng)市縣域經(jīng)濟(jì)與縣域基本競(jìng)爭(zhēng)力排名躍升至全國(guó)百?gòu)?qiáng)28位。瀏陽(yáng)河流域上游的株樹橋水庫(kù)是長(zhǎng)沙市的飲用水源地，影響著瀏陽(yáng)和長(zhǎng)沙數(shù)百萬人口和工農(nóng)業(yè)生產(chǎn)。但受氣候和地理因素的制約，旱澇頻率高，已成為影響瀏陽(yáng)市農(nóng)業(yè)生產(chǎn)的第一大氣象災(zāi)害，是制約該地區(qū)農(nóng)業(yè)經(jīng)濟(jì)持續(xù)發(fā)展的主要因素[11]。本文使用線性回歸、滑動(dòng)平均、Mann-Kendall檢驗(yàn)、累積距平聯(lián)合滑動(dòng)t聯(lián)合檢驗(yàn)和Morlet小波分析等統(tǒng)計(jì)方法，研究近45 a來瀏陽(yáng)河流域極端降水事件的演變規(guī)律，以期為當(dāng)?shù)卣畱?yīng)對(duì)極端氣候和旱澇災(zāi)害提供科學(xué)依據(jù)。

1 流域概況

瀏陽(yáng)河是湘江的一級(jí)支流，發(fā)源于大圍山北部的橫山坳，有大溪河和小溪河兩個(gè)源流，全長(zhǎng)235 km，流域面積4 665 km2。瀏陽(yáng)河流域地勢(shì)自東北向西南傾斜，以山地丘陵為主。流域?qū)賮啛釒駶?rùn)季風(fēng)氣候，氣候溫和，降水充沛，光照充足，全年四季差異明顯，年平均氣溫17.6℃，年平均降水量1 568 mm。瀏

陽(yáng)河流域水系及雨量站分布見圖1。

圖1瀏陽(yáng)河流域水系與雨量站點(diǎn)圖

2 資料來源與研究方法

2.1 資料來源

本研究采用的原始數(shù)據(jù)為瀏陽(yáng)河流域內(nèi)12個(gè)雨量站1969—2013年實(shí)測(cè)逐日降水?dāng)?shù)據(jù)，數(shù)據(jù)來源于湖南省水文水資源勘測(cè)局。所有數(shù)據(jù)經(jīng)過嚴(yán)格的質(zhì)量控制和錯(cuò)誤修正。為便于分析，依據(jù)世界氣象組織的規(guī)定，采用國(guó)際通用的百分位計(jì)算法來確定極端降水事件的閾值[12]。本文選取極端降水量、極端降水日數(shù)、極端降水強(qiáng)度和極端降水比率4個(gè)降水指數(shù)來反映瀏陽(yáng)河流域極端降水事件的基本特征(表1)。瀏陽(yáng)河流域各雨量站的極端降水指數(shù)通過RClimDex軟件計(jì)算得到。瀏陽(yáng)河全流域的極端降水指數(shù)則通過12個(gè)雨量站的算術(shù)平均值求得，并通過建立極端降水指數(shù)的時(shí)間序列來分析極端降水指數(shù)的年際變化特征。

表1 采用的極端降水指標(biāo)

2.2 研究方法

2.2.1 趨勢(shì)分析 采用線性回歸法和5 a滑動(dòng)平均法進(jìn)行極端降水時(shí)間序列趨勢(shì)分析，該方法能夠有效消除序列波動(dòng)的影響以及簡(jiǎn)單判斷時(shí)間序列的增減趨勢(shì)。

Mann-Kendall 非參數(shù)檢驗(yàn)法[13-14](以下簡(jiǎn)稱MK檢驗(yàn))是由世界氣象組織(WMO)推薦的一種非參數(shù)檢驗(yàn)法，常用于檢驗(yàn)降水、徑流和溫度等水文氣象要素的趨勢(shì)顯著變化，其優(yōu)點(diǎn)在于樣本可以不滿足一定的分布，同樣不受異常值的干擾，且計(jì)算簡(jiǎn)便。該法適用于水文氣象序列服從非正態(tài)分布并且已得到廣泛應(yīng)用。

2.2.2 突變分析 累積距平曲線聯(lián)合滑動(dòng)t聯(lián)合檢驗(yàn)[15]的原理是先通過找出累積距平曲線的最大值和最小值來初步確定為變點(diǎn),然后初步識(shí)別的變點(diǎn)利用滑動(dòng)t檢驗(yàn)進(jìn)行精確識(shí)別。

2.2.3 周期分析 小波變換具有時(shí)域和頻域上同時(shí)局部化的特點(diǎn)，水文氣象序列時(shí)域的變換具有局部性和多尺度性。小波分析給出了水文氣象要素時(shí)間序列不同層次下的變化尺度，通過小波系數(shù)分析可以識(shí)別降雨時(shí)間序列的多時(shí)間尺度的周期性和突變性[16]。因此，本文選用Morlet小波對(duì)瀏陽(yáng)河流域45 a的極端水文氣象序列進(jìn)行分析。

3 結(jié)果與分析

3.1 極端降水變化趨勢(shì)

用線性回歸和5 a滑動(dòng)平均法分析得到的瀏陽(yáng)河流域各極端降水指數(shù)的時(shí)間序列趨勢(shì)變化見圖2。分析可知，極端降水量呈上升趨勢(shì)，上升率為7.3 mm/10 a。極端降水量的多年平均值為468.41 mm，其中20個(gè)年份的極端降水量高于平均值，25個(gè)年份的極端降水量低于平均值，極端降水量最多的年份為1998年的1 102.5 mm，極端降水量最少的年份為2007年的119.5 mm。極端降水日數(shù)呈微弱下降趨勢(shì)，下降率為-0.054 d/10 a，極端降水日數(shù)的多年平均值為9 d，其中有15個(gè)年份的極端降水日數(shù)高于平均值，21個(gè)年份的極端降水日數(shù)低于平均值，4個(gè)年份極端降水日數(shù)等于平均值。極端降水日數(shù)最多的年份為1998年的15 d，極端降水日數(shù)最少的年份為2007年的3 d。極端降水強(qiáng)度總體呈上升趨勢(shì)，上升率為9.5 mm/(d·10 a)，極端降水強(qiáng)度的多年平均值為57.03 mm/d，其中有20個(gè)年份的極端降水強(qiáng)度高于平均值，25個(gè)年份的極端降水強(qiáng)度低于平均值。極端降水強(qiáng)度最大值為1997年的74.9 mm/d，極端降水強(qiáng)度最小值為2007年的44.95 mm/d。極端降水比率呈微弱上升趨勢(shì)，上升率為0.3%/10 a，表明極端降水對(duì)瀏陽(yáng)河流域降水量的影響越來越大。極端降水率的多年平均值為28%，其中21個(gè)年份的極端降水比率高于平均值，23個(gè)年份的極端降水比率低于平均值。極端降水比率最大值為1969年的50%，極端降水比率最小值為2007年的11%。

圖2極端降水指數(shù)的線性趨勢(shì)

計(jì)算瀏陽(yáng)河流域12個(gè)雨量站點(diǎn)的逐年極端降水指數(shù)的算術(shù)平均值，并建立各極端降水指數(shù)的時(shí)間序列。通過MK檢驗(yàn)獲得瀏陽(yáng)河流域各極端降水指數(shù)的統(tǒng)計(jì)量，在95%的置信區(qū)間下，極端降水指數(shù)中極端降水量、極端降水強(qiáng)度和極端降水比率均呈上升趨勢(shì)，極端降水日數(shù)呈下降趨勢(shì)，變化趨勢(shì)均不顯著(表2)，這與線性回歸趨勢(shì)分析的結(jié)果一致。

表2 極端降水指數(shù)MK統(tǒng)計(jì)量Z值

3.2 極端降水突變分析

從圖3可以看出，極端降水量累積距平曲線的最高點(diǎn)年份為1970年，最低點(diǎn)年份為1991年，由于1970年前后資料長(zhǎng)度的限制，無法判定其是否為變異點(diǎn)，因此初步選取1991年為極端降水量的變異點(diǎn)。同理選取極端降水日數(shù)、極端降水強(qiáng)度和極端降水比率的初步變點(diǎn)分別為1991年、1992年和1991年。

對(duì)累積距平曲線法初步確定的變異點(diǎn)，運(yùn)用滑動(dòng)t檢驗(yàn)法進(jìn)行精確識(shí)別。分別在變異點(diǎn)前取10 a和后取10 a，計(jì)算出變異點(diǎn)對(duì)應(yīng)的t值；通過表3可知，極端降水量和極端降水強(qiáng)度分別在1991年、1992年超過了α=0.05的顯著性水平，均發(fā)生顯著性均值突變。極端降水日數(shù)和極端降水比率在初始變點(diǎn)1991年均無顯著性均值突變發(fā)生。

3.3 極端降水周期分析

采用Morlet小波對(duì)近45 a來瀏陽(yáng)河流域極端降水指數(shù)進(jìn)行周期分析，得到極端降水指數(shù)距平序列小波變換系數(shù)的實(shí)部二維等值線圖，見圖4。實(shí)部圖中的實(shí)線表示小波變換系數(shù)為正值，表示極端降水指數(shù)處于偏大期；虛線表示為負(fù)值，表示極端降水指數(shù)處于偏小期。

圖3極端降水指數(shù)累積距平曲線

表3 極端降水指數(shù)滑動(dòng)t檢驗(yàn)

極端降水量存在8 a、11 a和25 a的主周期，以25 a左右周期的相位變化最為顯著，存在2個(gè)偏多期和3個(gè)偏少期，經(jīng)歷了“偏少—偏多—偏少—偏多—偏少”的5個(gè)交替變化。1978—1985年和1997—2002年處于正相位，極端降水為偏多期，其余時(shí)段處于負(fù)相位，極端降水為偏少期，且負(fù)相位的等值線未完全閉合，可以推斷2013年以后的一段時(shí)間極端降水仍然可能處于偏少期。以8 a為中心尺度的相位變化主要發(fā)生在1998—2013年，交替變化十分劇烈。以11 a為中心尺度的相位變化主要發(fā)生在1969—1985年。以上3個(gè)主周期波動(dòng)控制著瀏陽(yáng)河流域極端降水量在整個(gè)時(shí)間域內(nèi)的變化。極端降水日數(shù)存在8 a、11 a和27 a的主周期，以8 a左右周期的相位變化最為顯著。在1990年之前表現(xiàn)為11 a左右的周期變化，在1990—2000年表現(xiàn)為8 a左右的周期振蕩。以27 a左右為中心尺度的相位變化存在2個(gè)偏多期和1個(gè)偏少期。

極端降水強(qiáng)度存在10 a、18 a和27 a的主周期，以27 a左右周期的相位變化最為顯著，經(jīng)歷了“偏少—偏多—偏少—偏多—偏少”的5個(gè)交替變化。以18 a左右為中心尺度的相位變化存在4個(gè)偏多期和3個(gè)偏少期。以10 a左右為中心尺度的相位變化基本貫穿于整個(gè)時(shí)間序列，且交替變化十分劇烈。極端降水比率存在8 a、11 a和27 a的主周期，以27 a左右周期的相位變化最為顯著。在1969—1995年表現(xiàn)為11 a左右的周期變化，在1995年之后表現(xiàn)為8 a左右的周期變化。總之，極端降水指數(shù)的主周期和相位變化大體一致，都體現(xiàn)為小尺度的相位變化嵌套在大尺度的相位變化中。

4 結(jié)論與討論

(1) 通過線性回歸和5 a滑動(dòng)平均法分析得到的瀏陽(yáng)河流域極端降水量、極端降水強(qiáng)度和極端降水比率呈上升趨勢(shì)，極端降水日數(shù)呈下降趨勢(shì)，極端降水對(duì)瀏陽(yáng)河流域降水量的影響越來越大。在95%的置信區(qū)間下，采用MK檢驗(yàn)得到結(jié)果與線性回歸的結(jié)果一致，且變化趨勢(shì)在不同時(shí)期都不顯著。

(2) 通過累積距平聯(lián)合滑動(dòng)t聯(lián)合檢驗(yàn)發(fā)現(xiàn)瀏陽(yáng)河流域極端降水量和極端降水強(qiáng)度分別在1991年、1992年發(fā)生顯著性均值突變。極端降水日數(shù)和極端降水比率的初始變點(diǎn)在1991年，均未發(fā)生顯著性均值突變。

(3) 通過Morlet小波進(jìn)行周期分析可以發(fā)現(xiàn)，各極端降水指數(shù)均存在7～8 a、10～11 a和25～27 a左右的主周期，且豐枯交替變化明顯。各極端降水指數(shù)的相位時(shí)間變化尺度大體一致，均表現(xiàn)為小尺度嵌套于大尺度相位變化中。

圖4極端降水指數(shù)的小波變換系數(shù)二維等值線

(4) 引起極端降水變化的原因是復(fù)雜的，可能與人類活動(dòng)、大氣環(huán)流、太陽(yáng)黑子及地形等因素有關(guān)，各極端降水指數(shù)在空間上是如何變化的，具體的影響機(jī)理是什么，有待進(jìn)一步研究，以期為瀏陽(yáng)市的旱澇管理和防洪決策調(diào)度提供決策依據(jù)。

[1] 任正果,張明軍,王圣杰,等.1961—2011年中國(guó)南方地區(qū)極端降水事件變化[J].地理學(xué)報(bào),2014,69(5):640-649.

[2] 張強(qiáng),李劍鋒,陳曉宏,等.基于Copula函數(shù)的新疆極端降水概率時(shí)空變化特征[J].地理學(xué)報(bào),2011,66(1):3-12.

[3] 汪寶龍,張明軍,魏軍林,等.1960—2009年青海省極端降水事件的變化特征[J].水土保持通報(bào),2012,32(4):92-95.

[4] Wang S, Lei X H, Zhang F. Analysis on variation characteristics of precipitation time series in Weihe River Basin during the past 55 years[J]. Advanced Materials Research, 2013,684:246-252.

[5] 佘敦先,夏軍,張永勇,等.近50年來淮河流域極端降水的時(shí)空變化及統(tǒng)計(jì)特征[J].地理學(xué)報(bào),2011,66(9):1200-1210.

[6] 董滿宇,王炳欽,廖劍宇,等.近50年?yáng)|江流域極端降水事件變化特征[J].資源科學(xué),2013,35(3):521-529.

[7] 李運(yùn)剛,何大明,胡金明,等.紅河流域1960—2007年極端降水事件的時(shí)空變化特征[J].自然資源學(xué)報(bào),2012,27(11):1908-1917.

[8] 李占玲,王武,李占杰.基于GPD分布的黑河流域極端降水頻率特征分析[J].地理研究,2014,33(11):2169-2179.

[9] 李志,劉文兆,鄭粉莉.1965—2005年涇河流域極端降水事件的變化趨勢(shì)分析[J].資源科學(xué),2010,32(8):1527-1532.

[10] 肖艷,黎祖賢,章新平,等.近48年來湘江流域極端降水事件特征分析[J].長(zhǎng)江流域資源與環(huán)境,2010,19(11):1356-1362.

[11] 李春梅,賓素芝,凡安清,等.瀏陽(yáng)市干旱原因及對(duì)策[J].現(xiàn)代農(nóng)業(yè)科技,2013,(8):230-232.

[12] 安華,延軍平,張濤濤,等.增暖背景下華北平原極端降水事件時(shí)空變化特征[J].水土保持通報(bào),2013,33(3):144-148.

[13] Feng X, Zhang G, Yin X. Hydrological responses to climate change in Nenjiang river basin, northeastern China[J]. Water Resources Management, 2011,25(2):677-689.

[14] Chen H, Guo S, Xu C, et al. Historical temporal trends of hydro-climatic variables and runoff response to climate variability and their relevance in water resource management in the Hanjiang basin[J]. Journal of Hydrology, 2007,344(3):171-184.

[15] Shao Q, Li Z, Xu Z. Trend detection in hydrological time series by segment regression with application to Shiyang River Basin[J]. Stochastic Environmental Research and Risk Assessment, 2010,24(2):221-233.

[16] 霍正文,陳文,凡炳文.近54年定西市降水趨勢(shì)及突變分析[J].水文,2012,32(3):88-92.

ChangeinExtremePrecipitationinLiuyangRiverBasinDuringthePast45Years

YU Zexing1, HU Guohua1, CHEN Xiao1, GU Qingfu2, SHAO Quanxi1,3, SONG Hehua4

(1.SchoolofHydraulicEngineering,ChangshaUniversityofScience&Technology,Changsha410114,China; 2.Hydrology&WaterResourcesBureauofHunanProvince,Changsha410007,China; 3.CSIROData,FloreatWA6014,Australia; 4.Hydrology&WaterResourcesBureauofChangsha,Changsha410014,China)

Based on daily precipitation data from 12 weather stations in Liuyang River Basin during the period of 1969—2013 and the percentile-based threshold method, the extreme precipitation indices were analyzed by using linear regression, moving average, Mann-Kendall test, accumulative anomaly combined with movingttest and Morlet wavelet. The results indicated that: (1) extreme precipitation, precipitation intensity and precipitation ratio displayed increasing trends at the rates of 7.3 mm/decade, 9.5 mm/d and 0.3%, respectively, while the number of extreme precipitation days represented a decreasing trend at the rate of -0.054 d/decade; (2) abrupt changes in the extreme precipitation and precipitation intensity were detected in 1991 and 1992, respectively, while the number of extreme precipitation days and the ratio of extreme precipitation ratio to the total precipitation amounts were not statistically significant; (3) the extreme precipitation indices had the main oscillation cycles of 7～8 years, 10～11 years and 25～27 years, respectively, together with clear alternation between dry and wet years. The results are of great importance in flood control and disaster reduction in Liuyang River Basin.

extreme precipitation; Mann-Kendall test; morlet wavelet analysis; accumulative anomaly combined movingttest; Liuyang River Basin

2016-07-23

：2016-10-27

湖南省重大水利科技計(jì)劃項(xiàng)目(2013-243-2);湖南省(長(zhǎng)沙理工大學(xué))研究生科研創(chuàng)新項(xiàng)目(CX2015B359)

于澤興(1991—),男,湖南岳陽(yáng)人,碩士研究生,研究方向?yàn)樗念A(yù)報(bào)研究。E-mail:423552622@qq.com

胡國(guó)華(1965—),男,湖南新田人,教授,碩士研究生導(dǎo)師,研究方向?yàn)樗Y源演變規(guī)律研究。E-mail:g.h.hu@163.com

TV122+.1

：A

：1005-3409(2017)05-0139-05