• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    1951-2015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析

    2017-05-16 02:28:00胡毅鴻李景保
    農(nóng)業(yè)工程學報 2017年7期
    關(guān)鍵詞:洞庭湖區(qū)旱澇洪澇

    胡毅鴻,李景保

    (1. 湖南師范大學資源與環(huán)境科學學院,長沙 410081; 2. 北京師范大學減災與應急管理研究院,北京 100875)

    1951-2015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析

    胡毅鴻1,2,李景保1※

    (1. 湖南師范大學資源與環(huán)境科學學院,長沙 410081; 2. 北京師范大學減災與應急管理研究院,北京 100875)

    旱澇災害是制約洞庭湖區(qū)經(jīng)濟發(fā)展的重要因素,為預測湖區(qū)旱澇災害變化趨勢,該文選取1951-2015年洞庭湖區(qū)及其周邊地區(qū)35個氣象站點的逐月降水量數(shù)據(jù),利用標準化降水指數(shù)和夏季長周期旱澇急轉(zhuǎn)指數(shù),輔以小波分析方法,分析近65 a來研究區(qū)域的旱澇演變及典型年份旱澇急轉(zhuǎn)特征。結(jié)果表明:該區(qū)旱澇事件發(fā)生的頻率高,不同時間尺度旱澇事件發(fā)生的平均頻率為42.62%。在季節(jié)尺度上,夏季洪澇事件發(fā)生的頻率最高,為33.33%,冬季干旱事件發(fā)生的頻率最高,為41.54%;在年代際尺度上,20世紀70年代干旱事件發(fā)生的頻率最高,為25.83%,其次是2001-2015年,頻率為22.99%,90年代洪澇事件發(fā)生的頻率最高,為30%。該區(qū)域旱澇變化有4個特征時間尺度,分別為4、10、25和31 a,預測湖區(qū)在未來幾年降水偏多。洞庭湖區(qū)3個區(qū)域的旱澇變化大抵一致,但在不同時期仍存在一定的差異。湖區(qū)近65 a來有15 a的LDFAI絕對值大于1,旱澇急轉(zhuǎn)現(xiàn)象較頻繁,不管是旱澇變化還是旱澇急轉(zhuǎn)變化,東、西洞庭湖區(qū)的變化趨于一致,南洞庭湖區(qū)與東洞庭湖區(qū)旱澇急轉(zhuǎn)變化差異較大。

    氣象;干旱;降水;標準化降水指數(shù);旱澇急轉(zhuǎn)指數(shù);洞庭湖區(qū)

    胡毅鴻,李景保. 1951-2015年洞庭湖區(qū)旱澇演變及典型年份旱澇急轉(zhuǎn)特征分析[J]. 農(nóng)業(yè)工程學報,2017,33(7):107-115.doi:10.11975/j.issn.1002-6819.2017.07.014 http://www.tcsae.org

    Hu Yihong, Li Jingbao. Analysis on evolution of drought-flood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(7): 107-115. (in Chinese with English abstract)doi:10.11975/j.issn.1002-6819.2017.07.014 http://www.tcsae.org

    0 引 言

    IPCC第五次評估報告指出,氣候系統(tǒng)變暖是毋容置疑的事實[1]。氣候變暖不僅直接影響溫度極值波動,而且已經(jīng)導致高溫、干旱和暴雨洪澇等極端氣候事件發(fā)生頻率增高、強度加大,尤其是在氣候變化響應的敏感區(qū)和脆弱區(qū)[2-3]。近幾年來許多學者開展了很多具有區(qū)域特色的旱澇事件研究,并取得了一定成果[4-10]。

    素稱“漁米之鄉(xiāng)”的洞庭湖區(qū)受“馬蹄形”盆地格局和不穩(wěn)定的氣候系統(tǒng)(副熱帶高壓的北跳南移,西風環(huán)流的南侵北退,以及東南季風與西南季風的輻合交匯)的綜合作用,年降水量分配不均,年際變率大,溫度季節(jié)變化顯著,年際波動大,造成極端氣候事件頻發(fā),成為中國旱澇災害發(fā)生頻繁且嚴重的地區(qū)之一。

    研究表明,洞庭湖區(qū)旱澇致災因子的危險性大、承災體易損率高,因旱澇造成的直按經(jīng)濟損失巨大。在1951-2000年間由洪澇導致的多年平均成災面積為24.56萬hm2,其中1954、1996、1998年因洪澇造成的直按經(jīng)濟損失分別為121億元、303億元、197億元[11-15]。特別在干旱減少的年際背景下,極端干旱有上升趨勢[16],如2001年干旱成災面積18.3萬hm2,2006年達21.6萬hm2。可見,旱澇災害已成可持續(xù)發(fā)展的嚴重制約因素。

    長江中下游地區(qū)旱澇急轉(zhuǎn)事件頻繁,1960—2011年共有5 a出現(xiàn)了典型的旱澇急轉(zhuǎn)事件[17],其中以2011年旱澇急轉(zhuǎn)事件最為典型。2011年6月之前,該地區(qū)出現(xiàn)近60 a來最嚴重的冬春持續(xù)干旱,6月3日后出現(xiàn)5輪強降水過程,由大旱快速轉(zhuǎn)成大澇。這種極端現(xiàn)象造成多人遇難,經(jīng)濟損失巨大[18],嚴重影響工農(nóng)業(yè)生產(chǎn)和自然生態(tài)系統(tǒng)[19]。旱澇急轉(zhuǎn)是多種因素互相配合、持續(xù)異常的結(jié)果[19]。因此,旱澇急轉(zhuǎn)的規(guī)律、影響機制及預測一直是短期氣候與旱澇關(guān)系研究領(lǐng)域中的一個熱點和中心課題。

    鑒于此,基于近65 a來洞庭湖區(qū)及周邊地區(qū)35個氣象站點逐月降水數(shù)據(jù),利用標準化降水指數(shù)(standardized precipitation index,SPI)和夏季長周期旱澇急轉(zhuǎn)指數(shù)(long-cycle drought-flood abrupt alternation index,LDFAI),從時空尺度上,系統(tǒng)分析該地區(qū)旱澇時空變化特征以及夏季旱澇急轉(zhuǎn)特征。這對于較準確地把握洞庭湖區(qū)旱澇出現(xiàn)規(guī)律、轉(zhuǎn)換機制以及指導眾多水庫科學調(diào)度,緩解發(fā)電、灌溉與防汛間的矛盾具有重要的學術(shù)價值和實際意義。

    1 資料來源與研究方法

    1.1 研究區(qū)域及數(shù)據(jù)來源

    洞庭湖區(qū)(28°03′—30°20′N,110°40′—113°30′E)跨湘、鄂兩省的沖積平原和湖泊水網(wǎng)區(qū),總面積18 780 km2,其中湖南省15 200 km2,管轄岳陽、常德、益陽、長沙、湘潭、株洲6個地級市的36個縣(區(qū)、市)和14個國營農(nóng)場;湖北省3 580 km2,管轄荊州市的松滋縣、公安縣、石首市和江陵縣部分。氣象站點包括沅江、石門、常德、南縣、荊州、岳陽、平江等共35個(圖1),1951—2015年65 a各站點逐月降水量數(shù)據(jù)來自于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.gov.cn)提供的中國地面氣候資料月值數(shù)據(jù)集。

    圖1 氣象站點分布圖Fig.1 Distribution of meteorological statio ns

    1.2 研究方法

    1.2.1 標準化降水指數(shù)(SPI)

    降水量分布一般不是正態(tài)分布,而是一種偏態(tài)分布。所以在進行降水分析和干旱監(jiān)測、評估中,采用 Г分布概率來描述降水量的變化。標準化降水指數(shù) SPI是在計算某時段降水量的 Г分布概率后,再進行正態(tài)標準化處理,最終用標準化降水累積頻率分布來劃分干旱等級[20]。標準化降水指數(shù)具體計算方法參照《氣象干旱等級GB/T20481-2006》國家標準[21],具體計算步驟為:

    假設某一時段的降水量為隨機變量x,則其Г分布的概率密度函數(shù)g(x)為

    式中α>0,β>0,α和β分別為形狀參數(shù)和尺度參數(shù),x為降水量,Γ(α)是Gamma函數(shù)。最佳的α、β估計值可采用極大似然估計方法求得。確定概率密度函數(shù)中的參數(shù)后,對于某一年的降水量x0,對g(x)從0至x0進行數(shù)值積分可以計算出隨機變量x小于x0概率的近似估計值。由于Gamma方程不包含x=0的情況,所以降水量為0時的事件概率G由下式估計:G(x=0)=m/n,式中m為降水量為0的樣本數(shù),n為總樣本數(shù)。

    對Г分布概率進行正態(tài)標準化處理,將求得的概率值代入標準化正態(tài)分布函數(shù)即

    近似求解得到Z值,即為SPI值。

    SPI是McKee等提出來的,是一種功能強大、應用靈活、計算簡單的指數(shù),能很好地反映干旱的強度和持續(xù)時間[22],是一種非常有價值的干旱強度評估指數(shù)[23],優(yōu)于Z指數(shù)[24],目前已廣泛運用于各國的氣象旱澇檢測[25-29]。本文利用SPI分析洞庭湖區(qū)的旱澇變化特征,并分別計算1、3、6、12個月的SPI,其中SPI1表示1個月的SPI值,以反映短時間的降水狀況,比較接近土壤濕度;SPI3表示 3個月的 SPI值,以反映短期氣象干旱的特征,即降水季節(jié)性的變化;SPI6表示6個月的SPI值,反映中長時期的降水狀況;SPI12表示12個月的SPI值,以此指數(shù)來反映長時間降水的演變特征[30-31]。旱澇類型根據(jù)SPI值劃分為9個等級[20](表1)。

    表1 標準化降水指數(shù)的旱澇類型等級劃分Table1 Classification of grades of drought and flood for standardized precipitation index (SPI)

    1.2.2 降水量距平百分率

    降水量距平百分率(precipitation abnormity percentage,Pa)反映了某一時段降水與同期平均狀態(tài)的偏離程度,其計算公式如下[21]

    式中P為某時段降水量,mm,為多年平均同期降水量,mm。

    1.2.3 夏季長周期旱澇急轉(zhuǎn)指數(shù)

    旱澇急轉(zhuǎn)是指在某一段時期內(nèi),一段時間出現(xiàn)干旱,另一段時間出現(xiàn)洪澇,旱澇交替出現(xiàn)的情形,用夏季長周期旱澇急轉(zhuǎn)指數(shù)(long-cycle drought-flood abrupt alternation index,LDFAI)來衡量。

    1)將降水量進行標準化處理

    式中R為標準化降水量;Pi為降水量原值,mm,i=1,2,…,N;是降水量平均值,mm;N為樣本總數(shù)[32]。

    2)計算LDFAI

    式中R78是7—8月標準化降水量;R56是5—6月標準化降水量;(R78?R56)為旱澇急轉(zhuǎn)強度項;(|R56|+|R78|)為旱澇強度項;是權(quán)重系數(shù),其作用是增加長周期旱澇急轉(zhuǎn)事件所占權(quán)重,降低全旱或全澇事件權(quán)重[33-36]。

    1.2.4 LDFAI強度

    LDFAI能反映該區(qū)域是澇轉(zhuǎn)旱還是旱轉(zhuǎn)澇的過程,但是不能充分反映旱澇急轉(zhuǎn)的強度,尤其是在年代際尺度上,這是由于LDFAI有正有負,直接對其求平均會造成正負相抵,不能反映出旱澇急轉(zhuǎn)的強度,所以在計算時先對某年代各個站點的LDFAI取絕對值,再計算其平均值得到某一年代的LDFAI強度,來更好地反映該區(qū)域年代際尺度上旱澇急轉(zhuǎn)的強度變化。

    2 結(jié)果與分析

    2.1 不同時間尺度下 SPI的變化特征

    洞庭湖區(qū)SPI1在0線上下波動頻繁(圖2a),這是由于SPI1反映的是短時間降水狀況,沒有考慮前期降水的影響,降水時間的持續(xù)性弱,受短時間降水影響大,所以SPI1變化大。1951-2015年間有2.19%的月份屬于嚴重或極端旱澇月份,而這些極端和嚴重干旱月份78.57%是出現(xiàn)在秋冬季節(jié)。由此表明,秋冬季節(jié)易出現(xiàn)短時期的干旱事件,這是由于洞庭湖區(qū)地屬亞熱帶季風氣候,秋冬季節(jié)降水量較少。

    SPI3呈震蕩波動狀態(tài),旱澇交替頻繁(圖 2b)。近65 a間,2.60%屬于嚴重或極端旱澇事件,典型的季節(jié)性嚴重和極端干旱事件的年份為2011、1979和1999年,典型的嚴重和極端洪澇事件的年份為1954、1993、1998和2002年,這與洞庭湖區(qū)旱澇發(fā)生的實際情況基本一致[15,18]。由此認為,根據(jù)長時間序列逐月降水量數(shù)據(jù)計算得出的 SPI3能較好地反映洞庭湖區(qū)實際的旱澇情況。

    據(jù)洞庭湖區(qū)SPI6(圖2c),近65 a,1.56%屬于嚴重或極端旱澇事件,相比 SPI3其頻率稍有降低,但旱澇持續(xù)時間明顯延長。由此表明,隨著SPI時間尺度的增長,旱澇的等級會發(fā)生變化,且旱澇的開始時間與終止時間會相應地延后,這意味著時間尺度越長,越能充分反映前期降水對旱澇的累積影響。

    對于SPI12而言(圖2d),洞庭湖區(qū)近65 a的旱澇變化過程大致分為5個階段:1951—1955年相對濕潤,干旱事件出現(xiàn)頻率低;1956—1970年旱澇程度均不大,但干旱程度略高于洪澇;1971—1990年旱澇交替出現(xiàn),旱澇程度較前一階段大;1991—2003年為濕潤期,洪澇事件發(fā)生頻率高;2003年后又由洪澇向干旱轉(zhuǎn)變,且干旱事件發(fā)生頻率較高。長時間連續(xù)的洪澇或干旱事件是制約洞庭湖區(qū)農(nóng)業(yè)發(fā)展的重要因素,連續(xù) 8個月及以上洪澇事件大約每個年代分別出現(xiàn)一次(除90年代出現(xiàn)2次),維持時間最長的是1969年8月—1971年3月,長達20個月。20世紀60和80年代各出現(xiàn)1次連續(xù)8個月及以上干旱,70年代出現(xiàn)2次,而近十幾年來出現(xiàn)4次,這表明干旱事件有明顯增加的趨勢。

    根據(jù)不同時間尺度的SPI,洞庭湖區(qū)旱澇事件出現(xiàn)的頻率較高,不同時間尺度旱澇事件發(fā)生的平均頻率為42.62%,其中極端和嚴重旱澇事件發(fā)生的頻率較低,為1.84%,輕微旱澇事件發(fā)生的頻率最高,為31.85%。

    圖2 1951-2015年洞庭湖區(qū)不同時間尺度旱澇變化特征Fig.2 Characteristics of drought-flood changes with different time scales in Dongting Lake area during 1951-2015

    2.2 旱澇的時間變化特征

    2.2.1 不同等級旱澇頻次的季節(jié)變化

    洞庭湖區(qū)降水量年內(nèi)分配很不均勻,降水主要集中在夏季,近65 a來夏季平均降水量為498.69 mm,其次是春季為456.19 mm,秋季為247.20 mm,冬季降水最少,為149.54 mm,這與每個季節(jié)的干旱事件和洪澇事件發(fā)生的頻率大致對應(圖3)。春季,無極端干旱和洪澇事件,正常的頻率在四季中最高。夏季,降水普遍偏多,但在近10 a來夏季降水有減少的趨勢,極端和嚴重洪澇事件較多,無極端和嚴重干旱事件,洪澇事件發(fā)生的頻率最高,為 33.33%,干旱事件發(fā)生的頻率最少,為7.69%。秋季,無極端旱澇事件,有嚴重旱澇事件,干旱事件發(fā)生的頻次大于洪澇事件。冬季,降水普遍較少,極端和嚴重干旱事件發(fā)生最多,無極端洪澇事件,干旱事件發(fā)生的頻率最高,為 41.54%,洪澇事件發(fā)生的頻率最低,為7.69%,在四季中也是旱澇事件發(fā)生頻率最高的季節(jié)。近十幾年來,四季干旱事件發(fā)生的頻率都有所增高。

    圖3 1951—2015年洞庭湖區(qū)旱澇季節(jié)變化特征Fig.3 Seasonal variation characteristics of drought and flood in Dongting Lake area during 1951-2015

    2.2.2 不同等級旱澇頻次的年代際變化

    該地區(qū)不同年代極端旱澇事件發(fā)生頻率最低(表2),其頻率平均值僅為0.47%,正常的頻率最高,為55.34%,其次是輕微旱澇事件,為33.41%。其中90年代降水最多,平均降水量為1 422.73 mm,該年代發(fā)生了2次特大洪水,分別是1996和1998年,洪澇事件發(fā)生的頻率高達30%,嚴重干旱事件發(fā)生的頻率最高,為4.17%,該年代旱澇事件發(fā)生的頻率最高,為52.50%。其次是60年代,降水也偏多,平均降水量為1 366.13 mm,輕微洪澇事件發(fā)生的頻率最高,為23.33%。2001—2015年,降水最少,平均降水量為1 316.33 mm,極端干旱事件發(fā)生的頻率最高,為1.15%,干旱事件發(fā)生的頻率為22.99%。其次是70年代,平均降水量為1 325.50 mm,干旱事件發(fā)生的頻率最高,為25.83%。不同年代洪澇事件發(fā)生的頻率的變化大于干旱事件。

    表2 1951-2015年洞庭湖區(qū)不同旱澇等級出現(xiàn)頻率的年代際變化Table2 Interdecadal variation of frequency of drought and flood with different grades in Dongting Lake area during 1951-2015 %

    2.2.3 旱澇的周期變化

    將1951-2015年的SPI進行Morlet小波分析,以揭示洞庭湖區(qū)近65 a來旱澇的周期變化規(guī)律。圖4a的紅色(藍色)渲染表示小波系數(shù)實部為正數(shù)(負數(shù)),表明降水偏多(偏少),顏色越深,表示旱澇程度越大。圖 4b的小波系數(shù)模部平方相當于小波能量譜,數(shù)值越大(紅色越深)表示波動能量越強、周期越顯著。圖4c的小波方差能反映 SPI的波動能量隨時間尺度的分布情況,可確定旱澇變化過程中存在的主周期。在圖4b中有4個能量最聚集的中心,它們代表波動能量變化的特性,這 4個中心分別是:①尺度范圍在23~31 a,波動能量在90年代后期以來表現(xiàn)最為強烈。②中心尺度為31 a,其波動能量較強且貫穿整個時期。③中心尺度為10 a,波動能量主要影響的時域是70和80年代。④中心尺度為4 a,波動能量主要影響的時域是60年代后期—80年代前期以及2001年后。

    圖4 1951-2015年洞庭湖區(qū)旱澇變化小波分析Fig.4 Morlet wavelet analysis of drought-flood change in Dongting Lake area during 1951-2015

    這表明旱澇周期變化的局部化特征明顯,即在不同時域,變化周期不完全相同。在圖4c中有4個比較明顯的峰值,分別為4、10、25和31 a,說明這4個尺度的準周期在洞庭湖區(qū)近65 a來的旱澇變化中起主要作用,其中最大峰值對應著31 a特征時間尺度,表明其周期震蕩性最強、周期最顯著,為該地區(qū)旱澇變化的第一主周期,25、10和4 a分別為旱澇變化的第二、三和四主周期,這一結(jié)果與上述圖4b的分析結(jié)果大抵一致。具體表現(xiàn)在:31 a特征時間尺度大概經(jīng)歷了3個旱澇交替期,旱澇變化的平均周期為22 a左右,旱澇中心為1963、1974、1985、1998、2007、2015年;在25 a特征時間尺度上,大概有4個周期的旱澇交替期,平均變化周期為16 a左右;在10 a特征時間尺度上,大約有9個周期的旱澇交替期,平均變化周期為7 a左右;在4 a特征時間尺度上,大概有23個周期的旱澇交替期,平均變化周期為3 a左右。無論是哪個時間尺度,2015年后小波系數(shù)實部均呈現(xiàn)正值的趨勢,這意味著在未來一段時間內(nèi)洞庭湖區(qū)降水偏多。

    2.3 旱澇的空間分異特征

    為便于分析旱澇空間變化特征,按習慣將洞庭湖區(qū)劃分為東洞庭湖區(qū)、南洞庭湖區(qū)和西洞庭湖區(qū)。其中東、西和南洞庭湖區(qū)多年平均降水量分別為 1 364.95、1 373.18、1 531.33 mm,南洞庭湖區(qū)降水最豐富。為揭示三大湖區(qū)旱澇的分異特征,對各湖區(qū)的SPI12做散點圖進行分析(圖5)。東洞庭湖區(qū)與南、西洞庭湖區(qū)的相關(guān)系數(shù)都較高,說明東洞庭湖區(qū)與其他兩大湖區(qū)的旱澇變化有較強的相關(guān)性,但在20世紀80年代以來西洞庭湖區(qū)的干旱程度明顯大于東洞庭湖區(qū),這是因為80年代以來西洞庭湖區(qū)的降水量比東洞庭湖區(qū)偏少,西洞庭湖區(qū)1980—2015年的年均降水量為1 375.64 mm,而東洞庭湖區(qū)的降水量為1 395.57 mm。南洞庭湖區(qū)與東洞庭湖區(qū)的旱澇變化在20世紀90年代以前趨于一致,但在90年代后旱澇事件的差異變大。南洞庭湖區(qū)與西洞庭湖區(qū)的相關(guān)系數(shù)相對較小,旱澇變化的差異較大,甚至存在反相變化,在90年代出現(xiàn)過幾次明顯的反相變化。對相關(guān)系數(shù)進行顯著性檢驗,均通過0.01的信度檢驗。就整個湖區(qū)而言,旱澇的空間變化具有較強的一致性。

    圖5 洞庭湖區(qū)不同區(qū)域SPI相關(guān)分析Fig.5 Correlation of SPI in different areas of Dongting Lake area

    2.4 夏季長周期旱澇急轉(zhuǎn)特征

    2.4.1 典型年份夏季旱澇與旱澇急轉(zhuǎn)

    為揭示夏季旱澇與旱澇急轉(zhuǎn)特征,選取洞庭湖區(qū)近65 a來逐年5-8月的降水量,分別計算出5-8、5-6月以及7-8月的降水量距平百分率,作為反映該地區(qū)夏季降水變化的情況,同時選出5-8月Pa絕對值在15%以上的年份,以及5-6月和7-8月的Pa之差的絕對值在 30%以上的年份,從中找出典型的洪澇、干旱、旱轉(zhuǎn)澇以及澇轉(zhuǎn)旱事件各 3 a,分別得到這些典型年份的 Pa以及LDFAI(表3)。

    通過分析可得出以下認識:

    1)洪澇嚴重的年份,5-8月的Pa大于30%,降水量比長時間序列的平均降水量明顯偏多,整個夏季的降水都偏多。洪澇最嚴重的年份是1954年,5-8月的降水量偏多67.28%。該年洪澇事件屬于梅雨型澇災,由于受到強大且穩(wěn)定的蒙古高壓和西太平洋副熱帶高壓的控制,冷暖氣流長期在長江流域交鋒徘徊[37],使其入汛時間早、雨期長且雨量大、湖區(qū)連續(xù)暴雨,為洞庭湖區(qū)百年以來從未發(fā)生過的特大型洪澇事件,導致直接經(jīng)濟損失巨大。

    表3 洞庭湖區(qū)夏季典型年份的旱澇與旱澇急轉(zhuǎn)特征Table3 Features of drought-flood and its abrupt alternation in Dongting Lake area in summer of typical years

    2)干旱嚴重的年份,5-8月的Pa在?15%以下,7 -8月的Pa低于?40%,降水明顯偏少,嚴重影響該地區(qū)的農(nóng)業(yè)生產(chǎn)。

    3)典型澇轉(zhuǎn)旱年份,LDFAI小于?1,5-8月Pa的絕對值低于20%,5-6月的Pa為正值,7-8月的Pa為負值,且低于?30%,也就是說夏季前一段時間降水偏多后一段時間降水顯著偏少,導致 2個時期的降水量差距大,于是就形成澇轉(zhuǎn)旱事件。

    4)典型旱轉(zhuǎn)澇年份,LDFAI大于1,5-8月的Pa絕對值在10%以下,其降水量接近正常,但5-6月的Pa為負值,7-8月的Pa為正值,這與澇轉(zhuǎn)旱事件相反。綜上,典型旱澇年份夏季降水量普遍偏多或偏少,典型旱轉(zhuǎn)澇(澇轉(zhuǎn)旱)年份,5-6月降水偏少(偏多),7-8月降水偏多(偏少)。

    綜上,夏季偏澇(偏旱),整個夏季降水偏多(偏少),LDFAI的絕對值小于1;夏季出現(xiàn)旱轉(zhuǎn)澇(澇轉(zhuǎn)旱),一個時期降水偏多,另一個時期降水偏少,降水分配不均勻,LDFAI的絕對值大于1,且值越大,旱澇急轉(zhuǎn)的強度越大,旱澇急轉(zhuǎn)現(xiàn)象越顯著。對比洞庭湖區(qū)夏季降水與旱澇災情,LDFAI能較好地反映夏季5—6月和7—8月的旱澇急轉(zhuǎn)特征。

    2.4.2 旱澇急轉(zhuǎn)的時間變化

    據(jù)統(tǒng)計,20世紀60和90年代分別有2 a的LDFAI大于1,2 a的LDFAI小于-1,50和70年代及2001—2015年分別有1 a的LDFAI大于1,1 a的LDFAI小于-1,80年代只有1 a的LDFAI大于1,65 a中有15 a 的LDFAI絕對值大于1,每個年代出現(xiàn)的旱轉(zhuǎn)澇和澇轉(zhuǎn)旱現(xiàn)象頻率大抵一致。該地區(qū)不同年代的 LDFAI平均值均較小,LDFAI強度的平均值均較大(表 4),并且 80年代的LDFAI強度最小,出現(xiàn)旱澇急轉(zhuǎn)現(xiàn)象的頻率最低,90年代的LDFAI強度最大,出現(xiàn)旱澇急轉(zhuǎn)現(xiàn)象的頻率最高,說明不同年代的LDFAI強度能在一定程度上反映旱澇急轉(zhuǎn)現(xiàn)象出現(xiàn)的頻率。90年代的Pa最大,為該區(qū)域夏季最濕潤的一個年代,80年代為夏季最干燥的一個年代,其次是2001—2015年。

    表4 洞庭湖區(qū)夏季降水和長周期旱澇急轉(zhuǎn)指數(shù)的年代際統(tǒng)計特征Table4 Interdecadal statistics features of summer precipitation and long-cycle drought-flood abrupt alternation index (LDFAI) in Dongting Lake area

    夏季LDFAI有微弱上升的趨勢(圖6),洞庭湖區(qū)旱澇急轉(zhuǎn)現(xiàn)象大抵可分為4個階段,第1階段是1951—1970年,這一階段LDFAI呈波動性振蕩,旱轉(zhuǎn)澇和澇轉(zhuǎn)旱交替頻繁,且旱澇急轉(zhuǎn)強度較大,旱澇急轉(zhuǎn)現(xiàn)象顯著;第2階段是1971-1989年,這一階段除了1980年旱轉(zhuǎn)澇強度較大,總體來說LDFAI的波動較小,且強度也較小,旱澇急轉(zhuǎn)現(xiàn)象不顯著;第3階段是1990-2002年,這一階段LDFAI呈波動性振蕩,澇轉(zhuǎn)旱的強度呈逐漸減小的趨勢;第4階段是2003-2015年,除了2008年旱轉(zhuǎn)澇顯著,總體來說旱澇急轉(zhuǎn)現(xiàn)象不顯著。

    圖6 1951-2015年洞庭湖區(qū)夏季LDFAI變化Fig.6 Variation of LDFAI in Dongting Lake area during 1951-2015

    2.4.3 旱澇急轉(zhuǎn)的空間分異特征

    對各區(qū)域LDFAI做散點圖來分析其相關(guān)性,可知:西洞庭湖區(qū)和東洞庭湖區(qū)的相關(guān)系數(shù)較高(r=0.70,P<0.01),旱澇急轉(zhuǎn)指數(shù)的趨勢變化較一致,東洞庭湖區(qū)的澇轉(zhuǎn)旱程度大于西洞庭湖區(qū),西洞庭湖區(qū)的旱轉(zhuǎn)澇程度大于東洞庭湖區(qū),這說明該湖區(qū)不同區(qū)域雖然旱轉(zhuǎn)澇或澇轉(zhuǎn)旱的趨勢趨于一致,但旱澇急轉(zhuǎn)的強度不一致。西洞庭湖區(qū)與南洞庭湖區(qū)的相關(guān)性較弱(r=0.52,P<0.01),東洞庭湖區(qū)與南洞庭湖區(qū)的相關(guān)系數(shù)最小(r=0.35,P<0.01),其旱澇急轉(zhuǎn)指數(shù)趨勢變化的差異較大。不管是旱澇變化還是旱澇急轉(zhuǎn)變化,東、西洞庭湖區(qū)的趨勢變化趨于一致。

    3 結(jié) 論

    1)洞庭湖區(qū) 1951-2015年旱澇變化過程大致可以分為5個階段:1951-1955年相對濕潤;1956-1970年,旱澇程度不高;1971-1990年,旱澇交替出現(xiàn),旱澇程度較前一階段大;1991-2003年為較濕潤;2003年后較干旱,近十幾年來干旱事件有明顯增加的趨勢。

    2)洞庭湖區(qū)夏季洪澇事件發(fā)生的頻率最高,冬季干旱事件發(fā)生的頻率最高,近十幾年來,四季的干旱事件普遍增加。20世紀90年代降水最多,不同等級洪澇事件發(fā)生的頻率都較高,2001-2015年降水最少,極端干旱事件發(fā)生的頻率最高,70年代降水也偏少,干旱事件發(fā)生的頻率最高。

    3)洞庭湖區(qū)旱澇交替變化存在4、10、25和31 a這樣4個時間尺度,無論是哪個時間尺度,2015年后小波系數(shù)實部均呈現(xiàn)正值的趨勢,表明在未來一段時間內(nèi)該地區(qū)降水偏多。

    4)對比洞庭湖區(qū)夏季降水與旱澇災情,LDFAI能較好地反映夏季5—6月和7—8月的旱澇急轉(zhuǎn)特征,65 a中有15 a的LDFAI絕對值大于1。不管是旱澇變化還是旱澇急轉(zhuǎn)變化,東、西洞庭湖區(qū)的變化較為接近。南洞庭湖區(qū)與東洞庭湖區(qū)旱澇急轉(zhuǎn)變化差異較大。

    [1] IPCC. Climate Change 2013: The Physical Science Basis, the Summary for Policymakers of the Working Group I Contribution to the Fifth Assessment Report[M]. Cambridge: Cambridge University Press, 2013.

    [2] 第三次氣候變化國家評估報告編寫委員會. 第三次氣候變化國家評估報告[M]. 北京:科學出版社,2015.

    [3] 黃榮輝,杜振彩. 全球變暖背景下中國旱澇氣候災害的演變特征及趨勢[J]. 自然雜志,2010,32(4):187-195. Huang Ronghui, Du Zhencai. Evolution characteristic and trend of droughts and floods in china under the background of global warming[J]. Chinese Journal of Nature, 2010, 32(4): 187-195. (in Chinese with English abstract)

    [4] 李雙雙,楊賽霓,劉憲鋒. 1960-2013年北京旱澇變化特征及其影響因素分析[J]. 自然資源學報,2015,30(6):951-962. Li Shuangshuang, Yang Saini, Liu Xianfeng. The characteristics of drought-flood variation and its influence factors in Beijing during 1960-2013[J]. Journal of Natural Resources, 2015, 30(6): 951-962. (in Chinese with English abstract)

    [5] 劉琳,徐宗學. 西南地區(qū)旱澇特征及其趨勢預測[J]. 自然資源學報,2014,29(10):1791-1801. Liu Lin, Xu Zongxue. Drought/flood characteristic and trend prediction in Southwest China[J]. Journal of Natural Resources, 2014, 29(10): 1791-1801. (in Chinese with English abstract)

    [6] 楊金龍,張強,王勁松,等. 近60年來西南地區(qū)旱澇變化及極端和持續(xù)性特征認識[J]. 地理科學,2015,35(10):1336-1340. Yang Jinlong, Zhang Qiang, Wang Jinsong, et al. Extreme and persistent feature of drought and flood of Southwest China in past 60 years[J]. Scientia Geographica Sinica, 2015, 35(10): 1336-1340. (in Chinese with English abstract)

    [7] 何慧,廖雪萍,陸虹,等. 華南地區(qū)1961-2014年夏季長周期旱澇急轉(zhuǎn)特征[J]. 地理學報,2016,71(1):130-141. He Hui, Liao Xueping, Lu Hong, et al. Features of long-cycle drought-flood abrupt alternation in South China during summer in 1961-2014[J]. Acta Geographica Sinica, 2016, 71(1): 130-141. (in Chinese with English abstract)

    [8] 陳晶,顧世祥,何大明,等. 近60多年來滇池流域干旱特性及重現(xiàn)期分析[J]. 長江流域資源與環(huán)境,2016,25(9):1457-1465. Chen Jing, Gu Shixiang, He Daming, et al. Return period of drought in Dianchi Basin in last 69 years[J]. Resources and Environment in the Yangtze Basin, 2016, 25(9): 1457-1465. (in Chinese with English abstract)

    [9] 周玉良,周平,金菊良,等. 基于供水水源的干旱指數(shù)及在昆明干旱頻率分析中應用[J]. 水利學報,2014,45(9):1038-1047. Zhou Yuliang, Zhou Ping, Jin Juliang, et al. Establishment of hydrological drought index based on sources of regional water supply and its application to drought frequency analysis for Kunming[J]. Journal of Hydraulic Engineering, 2014, 45(9): 1038-1047. (in Chinese with English abstract)

    [10] 楊曉靜,左德鵬,徐宗學. 基于標準化降水指數(shù)的云南省近55年旱澇演變特征[J]. 資源科學,2014,36(3):473-480. Yang Xiaojing, Zuo Depeng, Xu Zongxue. Characteristics of droughts and flood analyzed using the standardized precipitation index in Yunan Province during the past 55 years[J]. Resource Science, 2014, 36(3): 473-480. (in Chinese with English abstract)

    [11] 劉世奇,段煉中,余元君. 洞庭湖水患災害分析及其治理策略[J]. 人民長江,2009,40(14):76-78. Liu Shiqi, Duan Lianzhong, Yu Yuanjun. Analysis on flood disaster in Dongting Lake and its control strategy[J]. Yangtze River, 2009, 40(14): 76-78. (in Chinese with English abstract)

    [12] 李景保,余果,歐朝敏,等. 洞庭湖區(qū)農(nóng)業(yè)水旱災害演變特征及影響因素:60年來的災情診斷[J]. 自然災害學報,2011,20(2):74-81. Li Jingbao, Yu Guo, Ou Chaomin, et al. Evolution characters and influence factors of agricultural flood and drought in Dongting Lake area diagnosis of disaster situation in recent 60 years[J]. Journal of Natural Disaster, 2011, 20(2): 74-81. (in Chinese with English abstract)

    [13] 黃菊梅,鄒用昌,彭嘉棟,等. 1960-2011年洞庭湖區(qū)年降水量變化特征[J]. 氣象與環(huán)境學報,2013,29(6):81-86. Huang Jumei, Zou Yongchang, Peng Jiadong, et al. Variation characteristic of annual precipitation from 1960 to 2011 in Dongting Lake area[J]. Journal of Meteorology and Environment, 2013, 29(6): 81-86. (in Chinese with English abstract)

    [14] 中華人民共和國水利部. 2014年中國水旱災害公報[R]. 北京:中國水利水電出版社,2011.

    [15] 王婷,章新平,黎祖賢,等. 近52年來洞庭湖流域氣象干旱的時空分布特征[J]. 長江流域資源與環(huán)境,2016,25(3):514-522. Wang Ting, Zhang Xinping, Li Zuxian, et al. Temporal and spatial distribution characteristic of meteorological drought for recent 52 years in Dongting Lake Basin[J]. Resources and Environment in the Yangtze Basin, 2016, 25(3): 514-522. (in Chinese with English abstract)

    [16] 李景保,王克林,楊燕,等. 洞庭湖區(qū)2000年2007年農(nóng)業(yè)干旱災害特點及成因分析[J]. 水資源與水工程學報,2008,19(6):1-5. Li Jingbao, Wang Kelin, Yang Yan, et al. Characteristic and causes of agricultural drought disasters from 2000 to 2007 in the areas of Dongting Lake[J]. Journal of Water Resources and Water Engineering, 2008, 19(6): 1-5. (in Chinese with English abstract)

    [17] 程智,丁小俊,徐敏,等. 長江中下游地區(qū)典型旱澇急轉(zhuǎn)氣候特征研究[J]. 長江流域資源與環(huán)境,2012,21(增刊2):115-120.Cheng Zhi, Ding Xiaojun, Xu Min, et al. Analysis of characteristic of sharp turn from drought to flood in the middle and lower reaches of the Yangtze River[J]. Resources and Environment in the Yangtze Basin, 2012, 21(Supp.2): 115-120. (in Chinese with English abstract)

    [18] 閃麗潔,張利平,陳心池,等. 長江中下游流域旱澇急轉(zhuǎn)時空演變特征分析[J]. 長江流域資源與環(huán)境,2015,24(12):2100-2107. Shan Lijie, Zhang Liping, Chen Xinchi, et al. Spatiotemporal evolution characteristics of drought-flood abrupt alternation in the middle and lower reaches of the Yangtze River Basin[J]. Resources and Environment in the Yangtze Basin, 2015, 24(12): 2100-2107. (in Chinese with English abstract)

    [19] 牛建利,何紫云,張?zhí)煊睿? 旱澇急轉(zhuǎn)對生產(chǎn),生活與生態(tài)的影響及應對措施效果分析[J]. 長江流域資源與環(huán)境,2013,22(增刊1):108-115. Niu Jianli, He Ziyun, Zhang Tianyu, et al. Impacts of sudden drought to flood change on manufacture livelihood and ecology[J]. Resources and Environment in the Yangtze Basin, 2013, 22(Supp.1): 108-115. (in Chinese with English abstract)

    [20] 史建國,嚴昌榮,何文清,等. 氣象干旱指數(shù)計算方法研究概述[J]. 中國農(nóng)業(yè)氣象,2007,28(增刊):191-195. Shi Jianguo, Yan Changrong, He Wenqing, et al. Overview on calculation methods of meteorological drought index[J]. Chinese Journal of Agrometeorology, 2007, 28(Supp.): 191-195. (in Chinese with English abstract)

    [21] 中國氣象局.氣象干旱等級GB/T20481-2006 [M]. 北京:中國標準出版社,2006.

    [22] McKee T B, Doesken N J, Kleist J. The relationship of drought frequency and duration to time scales[C]// Proceedings of the 8th Conference on Applied Climatology. Boston, USA: American Meteorological Society, 1993: 179-183.

    [23] 史建國,嚴昌榮. 氣象干旱指數(shù)計算方法研究概述[J]. 中國農(nóng)業(yè)氣象,2007,28(增刊):191-195. Shi Jianguo, Yan Changrong. Summary on the computing method of meteorological index of drought[J]. Chinese Journal of Agrometeorology, 2007, 28(Supp.): 191-195. (in Chinese with English abstract)

    [24] 袁文平,周廣勝. 標準化降水指標與Z指數(shù)在我國應用的對比分析[J]. 植物生態(tài)學報,2004,28(4):523-529. Yuan Wenping, Zhou Guangsheng. Comparison between standardized precipitation index and Z index in China[J]. Acta Phytoecologica Sinica, 2004, 28(4): 523-529. (in Chinese with English abstract)

    [25] 翟祿新,馮起. 基于 SPI的西北地區(qū)氣候干濕變化[J]. 自然資源學報,2011,26(5):847-857. Zhai Luxin, Feng Qi. Dryness/wetness climate variation based on standardized precipitation index in Northwest China[J]. Journal of Natural Resources, 2011, 26(5): 847-857. (in Chinese with English abstract)

    [26] 劉可晶,王文,朱燁,等. 淮河流域過去60年干旱趨勢特征及其與極端降水的聯(lián)系[J]. 水利學報,2012,43(10): 1179-1187. Liu Kejing, Wang Wen, Zhu Ye, et al. Trend of drought and its relationship with extreme precipitation in Huaihe Basin over the last 60 years[J]. Journal of Hydraulic Engineering, 2012, 43(10): 1179-1187. (in Chinese with English abstract)

    [27] 孫德亮,吳建峰,李威,等. 基于SPI指數(shù)的近50年重慶地區(qū)干旱時空分布特征[J]. 水土保持通報,2016,36(4):197-203. Sun Deliang, Wu Jianfeng, Li Wei, et al. Temporal and spatial patterns of droughts in recent 50 years of Chongqing city based on SPI index[J]. Bulletin of Soil and Water Conservation, 2016, 36(4): 197-203. (in Chinese with English abstract)

    [28] 王鶯,李耀輝,胡田田. 基于SPI指數(shù)的甘肅省河東地區(qū)干旱時空特征分析[J]. 中國沙漠,2014,34(1):244-253. Wang Ying, Li Yaohui, Hu Tiantian. Analysis on spatial and temporal patterns of drought based on standardized precipitation index in Hedong area in Gansu Province[J]. Journal of Desert Research, 2014, 34(1): 244-253. (in Chinese with English abstract)

    [29] 唐寶琪,延軍平,李雙雙,等. 近55年來華東地區(qū)旱澇時空變化特征[J]. 長江流域資源與環(huán)境,2016,25(3):497-505. Tang Baoqi, Yan Junping, Li Shuangshuang, et al. Evolution of spatiotemporal of floods and droughts in east China from 1960 to 2015[J]. Resources and Environment in the Yangtze Basin, 2016, 25(3): 497-505. (in Chinese with English abstract)

    [30] Heim Jr R R. A review of twentieth-century drought indices used in the United States[J]. Bulletin of the American Meteorological Society, 2002, 83(8): 1149.

    [31] Patel N R, Chopra P, Dadhwal V K. Analyzing spatial patterns of meteorological drought using standardized precipitation index[J]. Meteorological Applications, 2007, 14(4): 329-336.

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

    [33] 張水鋒,張金池,閔俊杰,等. 基于徑流分析的淮河流域汛期旱澇急轉(zhuǎn)研究[J]. 湖泊科學,2012,24(5):679-686. Zhang Shuifeng, Zhang Jinchi, Min Junjie, et al. Drought-flood abrupt alternation based on runoff in the Huaihe River Basin during rainy season[J]. Journal of Lake Science, 2012, 24(5): 679-686. (in Chinese with English abstract)

    [34] Wu Z, Li J, He J, et al. Occurrence of droughts and floods during the normal summer monsoons in the mid‐and lower reaches of the Yangtze River[J]. Geophysical Research Letters, 2006, 33, L05813. doi:10.1029/2005GL024487.

    [35] Wu Z, Li J, He J, et al. Large-scale atmospheric singularities and summer long-cycle droughts-floods abrupt alternation in the middle and lower reaches of the Yangtze River[J]. Chinese Science Bulletin, 2006, 51(16): 2027-2034.

    [36] 吳志偉,李建平,何金海,等. 大尺度大氣環(huán)流異常與長江中下游夏季長周期旱澇急轉(zhuǎn)[J]. 科學通報,2006,51(14):1717-1724.Wu Zhiwei, Li Jianping, He Jinhai, et al. The relationship between the circulation anomaly in large scale and long-periodic sharp turns from flood to drought over the middle and lower reaches of the Yangtze River in summer[J]. Chinese Science Bulletin, 2006, 51(14): 1717-1724. (in Chinese with English abstract)

    [37] 彭際作,毛德華. 洞庭湖區(qū)1998年與1954年特大洪澇災害比較研究[J]. 熱帶地理,2000,20(2):116-120. Peng Jizuo, Mao Dehua. A comparative study of the flood and waterlogging disasters occurred in 1998 and 1954 in the Dongting Lake region[J]. Tropical Geography, 2000, 20(2): 116-120. (in Chinese with English abstract)

    Analysis on evolution of drought-flood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area

    Hu Yihong1,2, Li Jingbao1※
    (1.College of Resource and Environmental Science,Hunan Normal University,Changsha410081,China; 2.Academy of Disaster Reduction and Emergency Management,Beijing Normal University,Beijing100875,China)

    Drought-flood disaster is one of vital factors restricting Dongting Lake area’s economic development. In order to predict the future drought-flood change tendency, using the monthly precipitation data of 36 meteorological stations in Dongting Lake and its surrounding areas, this paper attempted to describe the characteristics of evolution of drought-flood and its abrupt alternation in typical year from 1951 to 2015 in Dongting Lake area. Standardized precipitation index (SPI) and long cycle drought-flood abrupt alternation index (LDFAI), Morlet wavelet analysis and other climate diagnosis method were adopted. The results showed that the variation of drought and flood could be roughly divided into 5 periods. The first period was from 1951 to 1955, which was relatively humid. The second period was from 1956 to 1970. In this period, the alternation of flood and drought was frequent and the calamity degree was not severe. The third period was from 1971 to 1990. The alternation was also frequent and the calamity degree was severer than the previous period. The fourth was a humid period from 1991-2003. The last was an arid period after 2003. In this period, drought disaster has showed an obvious increasing trend. In the study area, flood and drought occurred frequently, and its alternation was also frequent during 1951-2015. The average frequency of different time scales of drought and flood events was 42.62%. The distribution of precipitation was uneven and concentrated in summer. The highest frequency of flood events was 33.33%, which occurred in summer. In winter, the frequency of drought events was 41.54%. Flood events showed an increasing trend in spring and summer in the 1990s. But in recent years, the rainfall has kept decreasing in 4 seasons. On an interdecadal scale, the 1990s was the wettest decade and had the highest frequency of flood events which was 30%. The precipitation was the least and the frequency of drought event was 22.99% in 2001-2015. In this period, the frequency of extreme drought event was the highest. The highest frequency of drought events was 25.83% in the 1970s. In this period, the precipitation was also low. Variation of drought-flood had 4 characteristic time scales, respectively, 4, 10, 25 and 31 a, which revealed that Dongting Lake area would be have more precipitation in the next few years. The variations of drought-flood in the 3 areas including East, West and South Dongting Lake area were basically the same, but had some differences in different periods. The variation of drought-flood even showed an obviously reverse phase in South and West Dongting Lake area in 1990s. There were 15 years with the absolute value of LDFAI is greater than 1 in recent 65 years in this area, which meant drought- flood abrupt alternation phenomenon was frequent. For either the variation of drought-flood or its abrupt alternation, West Dongting Lake area and East Dongting Lake area were basically the same. The variation of drought-flood abrupt alternation in South Dongting Lake area was obviously different from East and West Dongting Lake area.

    meteorology; drought; precipitation; standardized precipitation index; long-cycle drought-flood abrupt alternation index; Dongting Lake area

    10.11975/j.issn.1002-6819.2017.07.014

    P426.616

    A

    1002-6819(2017)-07-0107-09

    2016-08-15

    2017-01-25

    國家自然科學基金項目(41571100);湖南省重點學科地理學建設項目

    胡毅鴻,女,湖南沅江人,主要從事氣候變化研究。北京 北京師范大學減災與應急管理研究院,100875。

    Email:18374848655@163.com

    ※通信作者:李景保,男,湖南桂陽人,教授,主要從事水文水資源和自然災害教學與研究。長沙 湖南師范大學資源與環(huán)境科學學院,410081。

    Email:lijingbao1951@126.com

    猜你喜歡
    洞庭湖區(qū)旱澇洪澇
    洪澇造成孟加拉損失25.4萬噸大米
    衡陽暴雨洪澇特征及對后期農(nóng)業(yè)干旱影響
    四川省1960-2017年旱澇時空變化特征分析
    黔東南州洪澇發(fā)生規(guī)律研究
    基于GIS的洞庭湖區(qū)水土流失變化監(jiān)測
    氣候變化背景下西北干旱區(qū)旱澇的變化規(guī)律
    基于FloodArea 模型的龍須河流域暴雨洪澇淹沒模擬研究
    北方農(nóng)牧交錯帶氣候變化與旱澇響應特征
    寧夏旱澇災害對氣候變化的響應
    洞庭湖區(qū)荸薺高產(chǎn)優(yōu)質(zhì)栽培技術(shù)
    作物研究(2014年6期)2014-03-01 03:39:13
    黄色女人牲交| 欧美xxxx性猛交bbbb| 真人做人爱边吃奶动态| 少妇被粗大猛烈的视频| aaaaa片日本免费| 国产精品久久久久久久电影| 女人十人毛片免费观看3o分钟| 色5月婷婷丁香| 乱人视频在线观看| 午夜日韩欧美国产| 丰满乱子伦码专区| 99久久99久久久精品蜜桃| 长腿黑丝高跟| 国产精品一区二区性色av| 国产乱人伦免费视频| 日本三级黄在线观看| 久久久精品欧美日韩精品| 久久久久亚洲av毛片大全| 国产白丝娇喘喷水9色精品| 性色avwww在线观看| 99久久久亚洲精品蜜臀av| 三级毛片av免费| 亚洲美女黄片视频| 女人被狂操c到高潮| 成人三级黄色视频| 内地一区二区视频在线| 中文字幕人成人乱码亚洲影| 国产久久久一区二区三区| 免费黄网站久久成人精品 | 真人做人爱边吃奶动态| 老司机午夜福利在线观看视频| 99riav亚洲国产免费| 国产精品三级大全| 搡女人真爽免费视频火全软件 | 男人狂女人下面高潮的视频| 婷婷色综合大香蕉| www.熟女人妻精品国产| 日韩欧美三级三区| 在线十欧美十亚洲十日本专区| 亚洲自偷自拍三级| 成人欧美大片| 极品教师在线免费播放| 黄色一级大片看看| 国产成人影院久久av| 我的老师免费观看完整版| 免费在线观看亚洲国产| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 国产真实乱freesex| 久久人人精品亚洲av| 精品99又大又爽又粗少妇毛片 | 亚洲国产高清在线一区二区三| 欧美日韩瑟瑟在线播放| 一个人看视频在线观看www免费| 人妻制服诱惑在线中文字幕| 免费大片18禁| 精品不卡国产一区二区三区| 国产欧美日韩一区二区精品| 国产精品久久久久久亚洲av鲁大| 99久国产av精品| 在线免费观看的www视频| 少妇的逼水好多| 日日干狠狠操夜夜爽| 日本黄色视频三级网站网址| 国产麻豆成人av免费视频| av国产免费在线观看| 老女人水多毛片| 99国产综合亚洲精品| 亚洲专区国产一区二区| 三级毛片av免费| 婷婷六月久久综合丁香| 脱女人内裤的视频| 国产高清有码在线观看视频| 亚洲成人免费电影在线观看| 三级毛片av免费| 精品一区二区免费观看| 国产精品人妻久久久久久| 免费人成视频x8x8入口观看| 国产精品一区二区性色av| 蜜桃亚洲精品一区二区三区| 深夜精品福利| 国产成+人综合+亚洲专区| 美女cb高潮喷水在线观看| 日本与韩国留学比较| 观看免费一级毛片| 乱人视频在线观看| 久久精品91蜜桃| 男人狂女人下面高潮的视频| 国产精品一区二区三区四区久久| 成人永久免费在线观看视频| 美女cb高潮喷水在线观看| 亚洲精品日韩av片在线观看| 久久人妻av系列| 国内毛片毛片毛片毛片毛片| 午夜福利在线观看免费完整高清在 | av黄色大香蕉| 国产精品女同一区二区软件 | 国产探花在线观看一区二区| 最近中文字幕高清免费大全6 | 三级毛片av免费| 黄片小视频在线播放| netflix在线观看网站| 性色avwww在线观看| 无遮挡黄片免费观看| 国产成人a区在线观看| av视频在线观看入口| 97碰自拍视频| 久久久久久九九精品二区国产| 国产亚洲精品久久久久久毛片| 日韩欧美在线乱码| 国产高清三级在线| 91麻豆av在线| 久久伊人香网站| 国产高清有码在线观看视频| 国产精品久久久久久亚洲av鲁大| 十八禁国产超污无遮挡网站| 男人的好看免费观看在线视频| 美女大奶头视频| 亚洲中文日韩欧美视频| www.www免费av| 午夜a级毛片| 国产伦在线观看视频一区| 久久性视频一级片| 亚洲第一欧美日韩一区二区三区| 亚洲成av人片在线播放无| 亚洲精品在线美女| 香蕉av资源在线| 女人十人毛片免费观看3o分钟| 国产视频内射| 成人欧美大片| 在线观看免费视频日本深夜| 国产野战对白在线观看| 亚洲熟妇熟女久久| 两个人视频免费观看高清| 69人妻影院| 一边摸一边抽搐一进一小说| 午夜福利成人在线免费观看| 欧美日本视频| 欧美日韩黄片免| 午夜福利欧美成人| 久久久精品欧美日韩精品| 成年人黄色毛片网站| 在线免费观看不下载黄p国产 | 俺也久久电影网| 99热这里只有是精品在线观看 | av欧美777| or卡值多少钱| 99精品久久久久人妻精品| 国产成人啪精品午夜网站| 搡老熟女国产l中国老女人| 动漫黄色视频在线观看| 国产伦精品一区二区三区视频9| 中文字幕熟女人妻在线| 男人的好看免费观看在线视频| 国产真实乱freesex| 国产视频内射| 成人一区二区视频在线观看| 免费无遮挡裸体视频| 亚洲 欧美 日韩 在线 免费| 色哟哟哟哟哟哟| 国产精品一及| 69av精品久久久久久| 国产淫片久久久久久久久 | 欧洲精品卡2卡3卡4卡5卡区| 9191精品国产免费久久| 精品人妻1区二区| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频| 午夜福利18| 欧美成人一区二区免费高清观看| 少妇的逼水好多| 婷婷色综合大香蕉| 国产淫片久久久久久久久 | 国产成人福利小说| 亚洲avbb在线观看| 精品人妻偷拍中文字幕| 日日干狠狠操夜夜爽| 国产精品一区二区免费欧美| 日韩中字成人| 国产美女午夜福利| 变态另类成人亚洲欧美熟女| 免费大片18禁| 欧美潮喷喷水| 国产欧美日韩一区二区三| 欧美日韩福利视频一区二区| 少妇熟女aⅴ在线视频| 我的女老师完整版在线观看| 中出人妻视频一区二区| 大型黄色视频在线免费观看| 精品久久久久久久末码| 中文资源天堂在线| 99在线人妻在线中文字幕| 成年女人毛片免费观看观看9| 在线天堂最新版资源| or卡值多少钱| 怎么达到女性高潮| 热99在线观看视频| 欧美xxxx性猛交bbbb| 国产精品一区二区三区四区免费观看 | 两性午夜刺激爽爽歪歪视频在线观看| 麻豆国产97在线/欧美| 首页视频小说图片口味搜索| 国产免费一级a男人的天堂| 国产黄色小视频在线观看| 又黄又爽又刺激的免费视频.| 在线观看av片永久免费下载| 丝袜美腿在线中文| 国产伦精品一区二区三区四那| 婷婷色综合大香蕉| 久久久色成人| 在线观看66精品国产| a级毛片免费高清观看在线播放| 欧美高清性xxxxhd video| 成人av一区二区三区在线看| 夜夜爽天天搞| 好看av亚洲va欧美ⅴa在| 国产精品,欧美在线| 成年免费大片在线观看| 麻豆久久精品国产亚洲av| а√天堂www在线а√下载| 亚洲欧美日韩高清在线视频| 日韩中字成人| xxxwww97欧美| 免费在线观看影片大全网站| 18美女黄网站色大片免费观看| 在线国产一区二区在线| 国产成年人精品一区二区| 国产精品伦人一区二区| 精品人妻偷拍中文字幕| 青草久久国产| 国内精品一区二区在线观看| 90打野战视频偷拍视频| 天堂影院成人在线观看| 91av网一区二区| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| 波多野结衣高清作品| 91麻豆av在线| 国产三级在线视频| 日韩有码中文字幕| 国产午夜福利久久久久久| 欧美日韩乱码在线| 精品久久久久久久人妻蜜臀av| 国内揄拍国产精品人妻在线| 欧美中文日本在线观看视频| 嫩草影视91久久| 99riav亚洲国产免费| 亚洲一区二区三区不卡视频| 国产一区二区激情短视频| 全区人妻精品视频| 露出奶头的视频| 乱码一卡2卡4卡精品| 久久国产乱子伦精品免费另类| 免费观看的影片在线观看| 少妇熟女aⅴ在线视频| 久久久久久国产a免费观看| 丰满人妻熟妇乱又伦精品不卡| 午夜免费成人在线视频| 18禁黄网站禁片午夜丰满| 欧美日韩综合久久久久久 | 国产精品久久久久久人妻精品电影| 又黄又爽又刺激的免费视频.| 婷婷色综合大香蕉| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久久电影| 久久6这里有精品| 国产精品电影一区二区三区| www.色视频.com| 三级毛片av免费| 日韩中文字幕欧美一区二区| 性色av乱码一区二区三区2| 人妻丰满熟妇av一区二区三区| 亚洲片人在线观看| 老司机福利观看| 午夜福利在线在线| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 亚洲av一区综合| 国产探花在线观看一区二区| 亚洲,欧美精品.| 人妻久久中文字幕网| 中文字幕av成人在线电影| 97热精品久久久久久| 熟妇人妻久久中文字幕3abv| 久久久精品欧美日韩精品| 亚洲三级黄色毛片| 亚洲,欧美精品.| 亚洲男人的天堂狠狠| 自拍偷自拍亚洲精品老妇| 久久久色成人| 永久网站在线| 久久九九热精品免费| h日本视频在线播放| 国产一区二区三区在线臀色熟女| 人妻丰满熟妇av一区二区三区| 久久这里只有精品中国| 欧美中文日本在线观看视频| 久久精品人妻少妇| 国产精品美女特级片免费视频播放器| 亚洲中文字幕日韩| 久久午夜福利片| 亚洲最大成人中文| 超碰av人人做人人爽久久| 免费一级毛片在线播放高清视频| 变态另类丝袜制服| 精品久久久久久久久久免费视频| 最新中文字幕久久久久| 91麻豆av在线| 五月伊人婷婷丁香| 久久草成人影院| 婷婷精品国产亚洲av在线| 久久久色成人| 无人区码免费观看不卡| 麻豆成人av在线观看| 免费看a级黄色片| 欧美成人a在线观看| 午夜久久久久精精品| 久久热精品热| 偷拍熟女少妇极品色| 国产成人aa在线观看| 欧美绝顶高潮抽搐喷水| 欧美黄色片欧美黄色片| 欧美3d第一页| 动漫黄色视频在线观看| 国产精品自产拍在线观看55亚洲| 首页视频小说图片口味搜索| 亚洲色图av天堂| 亚洲熟妇熟女久久| 成人一区二区视频在线观看| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 欧美精品国产亚洲| 亚洲三级黄色毛片| 久久人人爽人人爽人人片va | 亚洲乱码一区二区免费版| 精品一区二区三区视频在线| 欧美日韩综合久久久久久 | 日本一二三区视频观看| 亚洲成人免费电影在线观看| 午夜福利视频1000在线观看| 999久久久精品免费观看国产| 成人无遮挡网站| 99久久九九国产精品国产免费| 观看美女的网站| 久久久久免费精品人妻一区二区| 极品教师在线免费播放| 欧美乱色亚洲激情| 最近中文字幕高清免费大全6 | 内射极品少妇av片p| 久久久久久大精品| 在线观看舔阴道视频| a级毛片免费高清观看在线播放| 日日夜夜操网爽| 久久人人爽人人爽人人片va | 久久亚洲精品不卡| 直男gayav资源| avwww免费| 99久国产av精品| 美女大奶头视频| 久久精品国产99精品国产亚洲性色| 波多野结衣高清无吗| 国产美女午夜福利| 亚洲精品在线观看二区| 无人区码免费观看不卡| av中文乱码字幕在线| 51国产日韩欧美| 久久久久久九九精品二区国产| 久久久国产成人免费| 免费在线观看影片大全网站| 亚洲人与动物交配视频| av中文乱码字幕在线| aaaaa片日本免费| 日本五十路高清| 99热6这里只有精品| 亚洲av成人不卡在线观看播放网| 99久久精品国产亚洲精品| 日韩欧美国产一区二区入口| 麻豆av噜噜一区二区三区| 2021天堂中文幕一二区在线观| 在线看三级毛片| 人人妻人人澡欧美一区二区| 极品教师在线视频| 久久精品国产亚洲av香蕉五月| 亚洲国产精品久久男人天堂| 免费电影在线观看免费观看| 亚洲欧美日韩高清在线视频| 久久精品综合一区二区三区| 国产亚洲精品综合一区在线观看| 欧美性猛交╳xxx乱大交人| 国产亚洲精品综合一区在线观看| 免费观看人在逋| 99国产精品一区二区三区| 天美传媒精品一区二区| 老熟妇乱子伦视频在线观看| 性色avwww在线观看| 日本黄色视频三级网站网址| 简卡轻食公司| 日本五十路高清| 亚洲一区高清亚洲精品| 男人舔女人下体高潮全视频| 久久伊人香网站| 国产精品久久久久久精品电影| 日韩欧美国产在线观看| 久久精品影院6| 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| 国产精品不卡视频一区二区 | 三级男女做爰猛烈吃奶摸视频| 在线观看午夜福利视频| av在线蜜桃| 国产精品美女特级片免费视频播放器| av中文乱码字幕在线| 亚洲一区高清亚洲精品| 日本 av在线| 天堂影院成人在线观看| 亚洲熟妇熟女久久| 日日摸夜夜添夜夜添av毛片 | 麻豆av噜噜一区二区三区| 精品人妻视频免费看| 成熟少妇高潮喷水视频| 色综合亚洲欧美另类图片| 男女那种视频在线观看| 亚洲av美国av| 99精品在免费线老司机午夜| 九九热线精品视视频播放| 欧美高清成人免费视频www| 日韩欧美免费精品| 国产精品98久久久久久宅男小说| 亚洲成av人片免费观看| 久久欧美精品欧美久久欧美| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久精品电影| 男女那种视频在线观看| 丰满乱子伦码专区| 在线十欧美十亚洲十日本专区| 性插视频无遮挡在线免费观看| 中文亚洲av片在线观看爽| 久久性视频一级片| 男人舔奶头视频| 成人鲁丝片一二三区免费| 国产不卡一卡二| 国产av一区在线观看免费| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 日本a在线网址| 国产毛片a区久久久久| 国模一区二区三区四区视频| 一区二区三区高清视频在线| 欧美黄色淫秽网站| 男女做爰动态图高潮gif福利片| 久久久久精品国产欧美久久久| 亚洲国产欧洲综合997久久,| 日韩欧美在线乱码| 91久久精品国产一区二区成人| 69av精品久久久久久| 欧美激情国产日韩精品一区| 精品久久久久久久末码| 每晚都被弄得嗷嗷叫到高潮| 97超级碰碰碰精品色视频在线观看| 精华霜和精华液先用哪个| 直男gayav资源| 久久中文看片网| 99精品在免费线老司机午夜| 亚洲av一区综合| 黄色一级大片看看| 国产亚洲精品久久久com| 亚洲专区中文字幕在线| 国产v大片淫在线免费观看| 九九久久精品国产亚洲av麻豆| 九色国产91popny在线| 国产精品三级大全| 天天一区二区日本电影三级| 国模一区二区三区四区视频| 午夜福利高清视频| 欧美3d第一页| 久久99热这里只有精品18| 又黄又爽又刺激的免费视频.| 国产精品日韩av在线免费观看| 久久久久久久久久成人| 欧美成人一区二区免费高清观看| 国产爱豆传媒在线观看| 性欧美人与动物交配| 国产免费av片在线观看野外av| 淫妇啪啪啪对白视频| 久久久久久久午夜电影| 欧美精品啪啪一区二区三区| 毛片女人毛片| av天堂在线播放| 蜜桃久久精品国产亚洲av| 亚洲自偷自拍三级| 色5月婷婷丁香| 国产精品久久视频播放| 亚洲专区中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 国产欧美日韩一区二区三| 天堂av国产一区二区熟女人妻| 18禁在线播放成人免费| 国产精品98久久久久久宅男小说| 能在线免费观看的黄片| 五月伊人婷婷丁香| 亚洲欧美日韩高清在线视频| 观看美女的网站| 久久香蕉精品热| 少妇人妻精品综合一区二区 | 精品一区二区三区视频在线| 深夜a级毛片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色av中文字幕| 桃红色精品国产亚洲av| 人人妻,人人澡人人爽秒播| 精品久久久久久成人av| 亚洲精品粉嫩美女一区| 成人国产综合亚洲| 成人av在线播放网站| 亚洲国产精品999在线| 久久国产精品人妻蜜桃| or卡值多少钱| 欧美最新免费一区二区三区 | 日本成人三级电影网站| x7x7x7水蜜桃| 九九久久精品国产亚洲av麻豆| 丁香六月欧美| 丝袜美腿在线中文| 老熟妇仑乱视频hdxx| 欧美成人免费av一区二区三区| 精品午夜福利在线看| 熟妇人妻久久中文字幕3abv| 国产亚洲av嫩草精品影院| 亚洲乱码一区二区免费版| 国产aⅴ精品一区二区三区波| 亚洲狠狠婷婷综合久久图片| 欧美精品啪啪一区二区三区| 久9热在线精品视频| 久久国产乱子伦精品免费另类| 乱码一卡2卡4卡精品| 舔av片在线| 日本成人三级电影网站| 99久久精品国产亚洲精品| 国产欧美日韩一区二区三| 免费av观看视频| 性插视频无遮挡在线免费观看| 国产乱人视频| 直男gayav资源| 午夜福利免费观看在线| h日本视频在线播放| 国产麻豆成人av免费视频| 国内揄拍国产精品人妻在线| 精品久久久久久久人妻蜜臀av| 男女视频在线观看网站免费| 亚洲精品亚洲一区二区| 国产精品永久免费网站| a级毛片a级免费在线| 少妇高潮的动态图| 可以在线观看的亚洲视频| 色噜噜av男人的天堂激情| av在线老鸭窝| 超碰av人人做人人爽久久| 精品一区二区三区视频在线观看免费| 国产精品伦人一区二区| 天堂√8在线中文| 黄色一级大片看看| 亚洲人与动物交配视频| 国产精品久久久久久精品电影| 亚洲成av人片在线播放无| www.www免费av| 亚洲乱码一区二区免费版| 精品久久久久久久久久免费视频| 69人妻影院| 一进一出抽搐gif免费好疼| 欧美另类亚洲清纯唯美| 亚洲无线在线观看| 国产视频一区二区在线看| 日本撒尿小便嘘嘘汇集6| 99热只有精品国产| 欧美精品啪啪一区二区三区| 精品人妻熟女av久视频| 久久婷婷人人爽人人干人人爱| 9191精品国产免费久久| 久久性视频一级片| eeuss影院久久| 国产av麻豆久久久久久久| 国产精品久久视频播放| 久久人人精品亚洲av| 亚洲欧美日韩高清在线视频| 欧美在线一区亚洲| 人人妻,人人澡人人爽秒播| av在线老鸭窝| 国产高清有码在线观看视频| 国产乱人视频| 国产大屁股一区二区在线视频| 99国产综合亚洲精品| 动漫黄色视频在线观看| 亚洲av不卡在线观看| 欧美性猛交黑人性爽| 亚洲人成网站在线播放欧美日韩| 夜夜夜夜夜久久久久| 久久精品影院6| 黄色日韩在线| 成人鲁丝片一二三区免费| 99在线人妻在线中文字幕| av福利片在线观看| 性欧美人与动物交配| 精品久久国产蜜桃| 99视频精品全部免费 在线| 老司机午夜福利在线观看视频| 国产一区二区激情短视频| 国产老妇女一区| 如何舔出高潮| 90打野战视频偷拍视频| 国产极品精品免费视频能看的| 中文字幕免费在线视频6| 舔av片在线| 1000部很黄的大片| 在线播放国产精品三级| 国产伦人伦偷精品视频|