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

    基于SPEI和SPI指數(shù)的太原多尺度干旱特征與氣候指數(shù)的關(guān)系

    2014-02-09 07:58:31張?jiān)儡?/span>郝智文李明明陳二萍張永瑞

    張?jiān)儡?,郝智文,?雁,李明明,陳二萍,李 芬,張永瑞

    1. 山西省氣象科學(xué)研究所,山西 太原 030002;2. 山西省氣象局影視中心,山西 太原 030002;3. 山西省雷電防護(hù)監(jiān)測中心,山西 太原 030002

    基于SPEI和SPI指數(shù)的太原多尺度干旱特征與氣候指數(shù)的關(guān)系

    張?jiān)儡?,郝智文1,王 雁1,李明明1,陳二萍1,李 芬2,張永瑞3

    1. 山西省氣象科學(xué)研究所,山西 太原 030002;2. 山西省氣象局影視中心,山西 太原 030002;3. 山西省雷電防護(hù)監(jiān)測中心,山西 太原 030002

    干旱是太原地區(qū)發(fā)生最頻繁的自然災(zāi)害之一,出現(xiàn)的次數(shù)多、持續(xù)的時間長,對國民經(jīng)濟(jì)特別是農(nóng)業(yè)生產(chǎn)造成了嚴(yán)重的影響。隨著全球氣候變暖、極端天氣出現(xiàn)得越發(fā)頻繁,太原干旱發(fā)生的頻率和危害程度均呈上升趨勢。基于1951—2012年月平均降水和氣溫資料,引入新的干旱指標(biāo):標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI),應(yīng)用SPEI和SPI(標(biāo)準(zhǔn)化降水指數(shù),簡稱:SPI)定量描述太原地區(qū)的62年來的干濕狀況;基于月尺度SPEI和SPI指數(shù),對太原月季尺度干旱變化做了分析,并利用交叉小波變換(XWT)探討了干旱與大尺度氣候因子之間的關(guān)系。結(jié)果表明:基于降水和蒸散的SPEI可以更靈活地反映月季干旱變化特征;交叉小波變換分析表明,太原地區(qū)的干旱與4個大尺度因子都具有6~12 a年代際主共振周期,在1980s都存在較好的相關(guān)性。SPEI和SPI與NAO通過顯著性檢驗(yàn)的6~12 a共振周期主要表在1985—2000年,序列在此頻段上表現(xiàn)出一定的正位相共振關(guān)系;SPEI和SPI與WP在1955—1960年和1990—2000年分別表現(xiàn)出2~3 a和4~8 a顯著的共振周期,存在明顯的滯后相關(guān),在1970—1990年具有10~16 a正位相顯著共振關(guān)系。SPEI與PDO在1986—2000年之間存在4~6、8~14 a兩個顯著的共振周期,各自表現(xiàn)出負(fù)、正位相共振關(guān)系,在1955—1960年存在2~3 a共振周期,在此頻段上SPEI顯著滯后于PDO。與SPEI相比,SPI與PDO的相關(guān)性較弱,僅在1955—1960、1986—2000年出現(xiàn)了弱的共振關(guān)系。SPEI和SPI與PNA在1983—1995年表現(xiàn)出4~7 a的顯著共振關(guān)系,反映了SPEI和SPI顯著滯后于PNA。

    SPEI;SPI;多尺度;干旱;氣候指數(shù)

    太原(37°27′~38°25′ N、111°30′~113°09′ E)位于山西省中部,太原盆地的北端,于華北地區(qū)黃河流域中部,西、北、東三面環(huán)山,中、南部為河谷平原。太原屬溫帶季風(fēng)性氣候,年平均降雨量在456 mm,年均蒸發(fā)量在743 mm,降水主要集中在5—9月份的汛期,降水的時空分布極不均勻,干旱發(fā)生頻繁,對國民經(jīng)濟(jì)特別是農(nóng)業(yè)影響極大。干旱定義和計(jì)量指標(biāo)的研究歷史久遠(yuǎn)且已碩果累累,但大多數(shù)干旱指標(biāo)都是針對具體的研究目的而設(shè)定的,具有明顯的地域性和特定的時間尺度,目前常用于農(nóng)業(yè)干旱監(jiān)測、預(yù)報(bào)的干旱指數(shù)有帕爾默干旱強(qiáng)度指數(shù)(PDSI)(Palmer,1965)、標(biāo)準(zhǔn)化降水指數(shù)(SPI)(McKee等,1993)、綜合氣象干旱指數(shù)(CI)(李樹巖等,2009;張調(diào)風(fēng)等,2012;曹永強(qiáng)等,2012)。PDSI是基于地表水分平衡原理的干濕指標(biāo),其計(jì)算過程復(fù)雜、參數(shù)區(qū)域性較強(qiáng),且具有固定的時間尺度,這些都使它的應(yīng)用受到一定的限制(姚玉璧等,2011);而SPI只需要降水作為輸入量,計(jì)算簡單易行且具有多時間尺度的特性,使得用同一個干旱指標(biāo)反映不同時間尺度和不同方面的水資源狀況成為可能(王媛媛等,2012)。但是,在全球變化背景下,氣溫的升高已成為加劇干旱過程的重要因子之一(程國棟等,2006),對干旱狀況的客觀表征需要綜合降水和氣溫變化的共同效應(yīng)。為此,Vicente-Serranotal等(2010)在SPI基礎(chǔ)上,提出了一個新的氣候干旱指數(shù):標(biāo)準(zhǔn)化降水蒸散指數(shù)(Standardized Precipitation Evapotranspiration Index,SPEI)。該指數(shù)基于降水和溫度兩大因子,有機(jī)地集成了PDSI對蒸發(fā)需求變化(溫度波動所造成)的靈敏性和SPI的計(jì)算簡單、多時間尺度的屬性,是監(jiān)測干旱化及研究增溫影響干旱化過程較為理想的工具。

    在干旱形成機(jī)理方面,過去的研究主要集中在局地陸氣相互作用對干旱形成的相對貢獻(xiàn)等問題,但現(xiàn)在人們已意識到全球海洋溫度的變化和大尺度氣候因子[如受厄爾尼諾-南方濤動(El Nino-Southern Oscillation,ENSO)]、北大西洋濤動(North Atlantic Oscillation,NAO)、太平洋濤動(Pacific Decadal Oscillation,PDO)等在年代際尺度上對干旱形成的作用可能至關(guān)重要(黃榮輝等,2011)。

    本文應(yīng)用新的干旱指數(shù)SPEI來綜合表征太原地區(qū)1951年以來的干旱狀況;采用交叉小波變換(cross wavelet transform,XWT)探討太原干旱狀況與四大氣候因子(NAO、WP、PDO、PNA)之間的多時間尺度相關(guān)關(guān)系,為未來的干旱狀況客觀評估和準(zhǔn)確預(yù)測提供參考依據(jù)。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)

    本文所用的月平均降水和氣溫資料(1951—2012年)來自于山西省觀象臺(112.5° E,37.73° N)觀測數(shù)據(jù)。NAO、WP、PNA等氣候指數(shù)數(shù)據(jù)來自NOAA(美國國家海洋與大氣管理局)http://www.esrl.noaa.gov/psd/data。PDO采用Mantua (Mantua等,1997)提出的指數(shù), 1900—2010年逐月數(shù)據(jù)取自于http://jisao.washington.edu/pdo/PDO.latest。

    1.2 SPEI和SPI干旱指數(shù)

    采用具有多時間尺度特征的標(biāo)準(zhǔn)化降水指數(shù)(SPI)與在SPI基礎(chǔ)上考慮蒸散影響的標(biāo)準(zhǔn)化降水蒸散指數(shù)(SPEI)作為氣候干濕狀況的表征。不同時間尺度干旱指數(shù)的第一個月對應(yīng)的是當(dāng)前月,對n個月的時間尺度,從時間上自當(dāng)月向前延續(xù)n-1個月。SPI指數(shù)是表征某時段降水量出現(xiàn)概率多少的指標(biāo)之一,該指標(biāo)適合月尺度以上相對于當(dāng)?shù)貧夂驙顩r的干旱監(jiān)測與評估。SPI指數(shù)的構(gòu)建過程如下:

    假設(shè)某一時段的降水量為x,則其Γ分布的概率密度函數(shù)為:

    式中:β為尺度參數(shù);γ為形狀參數(shù);x為降水量;β和γ可用極大似然估計(jì)方法求得:

    式中:xi為降水資料樣本,為降水量氣候平均值,n為計(jì)算序列的長度。確定概率密度函數(shù)中的參數(shù)后,對于某一年的降水量x0,可求出隨機(jī)變量x小于x0事件的概率為:

    利用數(shù)值積分可以計(jì)算用(1)和(5)式后的事件概率近似估計(jì)值。降水量為0時的事件概率為:

    式中:m為降水量為0的樣本數(shù);n為總樣本數(shù)。

    對Γ分布概率進(jìn)行正太標(biāo)準(zhǔn)化處理,即將式(5)、(6)求得的概率值代入標(biāo)準(zhǔn)化正態(tài)分布函數(shù),即:

    對式(7)進(jìn)行近似求解可得:

    SPEI指數(shù)的計(jì)算與SPI指數(shù)類似,計(jì)算時首先以某一時段降水量與潛在蒸散量的差值替代降水量,然后進(jìn)行正態(tài)標(biāo)準(zhǔn)化,潛在蒸散采用Penman-Monteith公式(Allen等,1998)計(jì)算。

    根據(jù)國家氣象干旱等級標(biāo)準(zhǔn),劃分為5個等級,并確定相應(yīng)的SPEI/SPI界限值,各級干旱指標(biāo)見表1。

    2 結(jié)果分析

    2.1 SPEI和SPI驗(yàn)證:表征干旱的多尺度特征

    從圖1可以看出,不同時間尺度的SPEI(SPEI-1、SPEI-12、SPEI-24和SPEI-48分別表示時間尺度為1、12、24和48個月時間尺度的SPEI序列)值隨時間變化的敏感性明顯不同,時間尺度越小,干濕變化越顯著,其值會發(fā)生較大變化,甚至是正負(fù)波動。相反,時間尺度越大則干濕交替轉(zhuǎn)換越平緩,只有多次持續(xù)的降水或無雨、高溫等情況下才會使之發(fā)生轉(zhuǎn)變,這對于監(jiān)測長時期的干旱狀況是合理的(劉敏等,2013)。SPEI表征干旱多時間尺度特性的能力為其在太原地區(qū)干旱監(jiān)測與預(yù)報(bào)預(yù)警的推廣應(yīng)用奠定了良好基礎(chǔ)。

    1個月尺度的SPEI因?yàn)閷τ诙唐诮邓蜏囟茸兓急容^敏感,其數(shù)值波動頻率較高,波動幅度較大(圖1a),充分體現(xiàn)了太原月季尺度旱澇頻繁交替的特性。太原干旱12個月尺度(圖1b)的干旱發(fā)生頻率也比較高,在1958、1965、1972、1986、1997年尤為嚴(yán)重。另外,太原地區(qū)也存在24、48個月尺度的干旱特點(diǎn)(圖1c、d),如1964—1970、1980—1990和1994—2001年,據(jù)文獻(xiàn)記載1986年汾河年均流量僅占多年平均徑流量的39%,1980s—1990s,則是建國以來最嚴(yán)重的連續(xù)干旱年,水庫水位急速下降,河道大多干枯,地下水位急劇下降(溫克剛等,2005)。而SPEI-24和SPEI-48在這3個時期內(nèi)持續(xù)保持較大的負(fù)值,與歷史記錄的較大規(guī)模的連續(xù)干旱有較好的對應(yīng)關(guān)系。由此可見,多時間尺度的SPEI能夠清楚地反映不同方面的干濕演變與可獲得水資源狀況;SPEI-1主要受每月的水分和氣溫變化的影響,可以準(zhǔn)確地反映土壤水分狀況,可用于適時排水灌溉和保證農(nóng)業(yè)生產(chǎn);SPEI-12對于下層土壤水分和河流徑流量等有較好地反映;SPEI-24和SPEI-48則更深入地涉及到地下水位、水庫蓄水量及其民生問題。綜上所述,多時間尺度的SPEI值可以有效地反映太原地區(qū)旱澇程度及其持續(xù)時間。

    表1 SPEI和SPI干旱等級劃分Table 1 SPEI and SPI grade standard of drought

    圖1 1951—2012年時間尺度分別為1、12、24和48個月的太原歷年SPEI序列值對比Fig.1 Comparison of SPEI series in the time scales of 1, 12, 24 and 48 months in Taiyuan from 1951 to 2012

    2.2 SPEI和SPI干旱指數(shù)月季變化特征分析

    基于月尺度SPEI和SPI指數(shù),對太原市干旱變化做了分析。SPEI中位數(shù)表現(xiàn)出了明顯的月季變化(圖2a),冬季月份SPEI干旱指數(shù)收斂性較強(qiáng),變化幅度較小,中位數(shù)總體上位于-1.0±0.1之間,春季和秋季SPEI干旱指數(shù)較為離散,變化幅度較大,5和9月份SPEI中位數(shù)處于0.8±0.1之間,而3和11月份SPEI中位數(shù)則處于-0.8±0.1之間;夏季月份SPEI中位數(shù)在1.0附近,SPEI變化幅度在1.5左右。與SPEI指數(shù)相比,SPI中位數(shù)全年都集中于0±0.2之間(圖2b),基本上沒有反映出月(季)變化,SPI指數(shù)振幅較大,春季、夏季和秋季表現(xiàn)均比較明顯,冬季SPI變化幅度較小。上述分析表明,基于降水和蒸散的SPEI可以更靈活的反映月季干旱變化特征,而SPI指數(shù)忽略了溫度變化對蒸散的影響。

    從SPEI/SPI概率密度分布來看(圖2c),SPEI和SPI的概率分布較為吻合,出現(xiàn)重旱等級為10%~15%,出現(xiàn)中旱等級的頻率15%~20%,出現(xiàn)輕旱等級的頻率在25%~30%,無旱等級出現(xiàn)頻率為 30%~35%;出現(xiàn)特旱等級概率有所區(qū)別,SPEI為5%~10%,而SPI為2%~10%。

    圖2 不同月份SPEI和SPI干旱指數(shù)箱線圖比較及其概率密度分布Fig.2 Drought index box plot comparison in different months and probability density

    3 太原干旱與氣候因子的關(guān)系

    基于月尺度的SPEI和SPI分析了不同月份太原干旱與相應(yīng)時間序列氣候指數(shù)的相關(guān)性(圖3)。從圖3a中可以看出,SPEI和SPI對NAO指數(shù)的響應(yīng)在1—4月份比較顯著,這可能是因?yàn)镹AO在冬半年影響較大,而在夏半年影響較弱(左金清等,2012)。SPEI和SPI與WP指數(shù)相關(guān)性最為密切的時間在2、4、7和11月份(圖3b),當(dāng)冬季W(wǎng)P型出現(xiàn)時,對應(yīng)的500 hPa高度場表現(xiàn)出以45°N為界線的南北偶極子式異常型,同時,東亞-西太平洋地區(qū)對流層中低層的大氣溫度場和高層西風(fēng)急流,以及西北太平洋海表溫度異常也表現(xiàn)出一定程度的南北偶極子型分布,而在海平面氣壓場則表現(xiàn)為北太平洋濤動,對應(yīng)于冬季W(wǎng)P型的高、低指數(shù)年,亞洲大陸上空的西伯利亞高壓、中緯度西風(fēng),以及東亞冬季風(fēng)呈現(xiàn)出一定程度的反相變化特征(李勇等,2007)。SPEI和SPI與PNA指數(shù)相關(guān)性最為密切的時間在1和2月份(圖3c),厄爾尼諾、拉尼娜與WP、PNA遙相關(guān)環(huán)流型的強(qiáng)度關(guān)系密切。厄爾尼諾發(fā)展到盛期冬季大氣環(huán)流具有明顯的特征是:弱WP,強(qiáng)PNA,大氣環(huán)流及遙相關(guān)型強(qiáng)度的這種年代際變化是中國冬季氣候變化的一個重要原因(張?jiān)儡姷龋?013)。SPEI和SPI與PDO在9和10月份呈現(xiàn)較為顯著的負(fù)相關(guān)關(guān)系(圖3d),在PDO冷(暖)位相期,當(dāng)ENSO事件處于發(fā)展階段,華北地區(qū)降水偏多(少)(朱益民等,2003)。

    SPEI和SPI與各氣候指數(shù)的相關(guān),總體上保持了較好的一致性,但是在個別月份SPI和氣候指數(shù)的相關(guān)略弱與SPEI和氣候指數(shù)的相關(guān),這種情況可能反應(yīng)了溫度變化對干旱的影響作用。

    太原主雨季(6—9月份)SPEI和SPI干旱指數(shù)與同期大尺度氣候指數(shù)的交叉小波分析,重點(diǎn)突出太原干旱變異與各大尺度氣候因子在時頻域中高能量區(qū)的相互關(guān)系(圖4)。從圖4a、b中看出,SPEI和SPI與NAO通過顯著性檢驗(yàn)的6~12 a共振周期主要表現(xiàn)在1985—2000年,平均位相角右下接近45°,SPEI和SPI略超前于NAO,序列在此頻段上表現(xiàn)出一定的正位相共振關(guān)系,即存在正相關(guān)關(guān)系。SPEI和SPI與WP在1955—1960和1990—2000年,分別表現(xiàn)出2~3和4~8 a顯著的共振周期,位相角垂直向上,SPEI和SPI顯著滯后于WP,1970—1990年具有10~16 a正位相顯著共振關(guān)系,即存在正相關(guān)關(guān)系(圖4c、d)。SPEI與PDO在1986—2000年之間存在4~6、8~14 a兩個顯著的共振周期,各自表現(xiàn)出負(fù)、正位相共振關(guān)系,在1955—1960年存在2~3 a共振周期,平均位相角垂直向上90°,說明SPEI顯著滯后于PDO(圖4e),與SPEI相比SPI和PDO的相關(guān)性較弱,僅在1955—1960、1986—2000年出現(xiàn)了弱的共振關(guān)系(圖4f)。SPEI和SPI與PNA在1983—1995年表現(xiàn)出4~7 a的顯著共振關(guān)系,平均位相角垂直向上90°,反映了SPEI和SPI均滯后于PNA(圖4g、h)。

    從上述分析可以看出,1980s—1990s太原干旱與四種氣候指數(shù)都存在顯著的關(guān)系。大氣遙相關(guān)(NAO、WP、PNA)現(xiàn)象可以在一定程度上解釋全球相關(guān)聯(lián)的大范圍天氣氣候異常,同時,對于某個區(qū)域氣候變化來說,遙相關(guān)型的氣候指數(shù)也是有用的指標(biāo)(施能,1996)。

    圖3 SPEI和SPI與氣候指數(shù)的相關(guān)系數(shù)Fig.3 The correlation coefficient SPEI/ SPI and the climate index

    4 結(jié)論

    本文引入了綜合降水和氣溫變化共同效應(yīng)的一個新氣候干旱指數(shù)SPEI,它集合了PDSI對溫度波動的靈敏性和SPI的計(jì)算簡單、多時間尺度的自然屬性。經(jīng)過與太原歷史數(shù)據(jù)的對比驗(yàn)證表明,SPEI 可以在多時間尺度上有效地反映太原地區(qū)旱澇程度及其持續(xù)時間,是研究該地區(qū)干旱化過程較為理想的工具。

    1)基于降水和蒸散的SPEI可以更靈活的反映月季干旱變化特征,而SPI指數(shù)沒有溫度變化對蒸散的影響,不利于表征月季干旱變化特征。

    2)SPEI和SPI與各氣候指數(shù)的相關(guān),總體上保持了較好的一致性,但是在個別月份SPI和氣候指數(shù)的相關(guān)略弱與SPEI和氣候指數(shù)的相關(guān)。

    3)交叉小波變換分析表明,太原地區(qū)的干旱與4個大尺度因子都具有6~12 a年代際主共振周期,在1980s都存在較好的相關(guān)性。SPEI和SPI與NAO通過顯著性檢驗(yàn)的6~12 a共振周期主要表現(xiàn)在1985—2000年,序列在此頻段上表現(xiàn)出一定的正位相共振關(guān)系;SPEI和SPI與WP在1955—1960和1990—2000年,分別表現(xiàn)出2~3和4~8 a顯著的共振周期,存在明顯的滯后相關(guān),在1970—1990年具有10~16 a正位相顯著共振關(guān)系。SPEI與PDO在1986—2000年之間存在4~6、8~14 a兩個顯著的共振周期,各自表現(xiàn)出負(fù)、正位相共振關(guān)系,在1955—1960年存在2~3 a共振周期,在此頻段上SPEI顯著滯后于PDO。與SPEI相比SPI與PDO的相關(guān)性較弱,僅在1955—1960、1986—2000年出現(xiàn)了弱的共振關(guān)系。SPEI和SPI與PNA在1983—1995年表現(xiàn)出4~7 a的顯著共振關(guān)系,反映了SPEI和SPI均顯著滯后于PNA。

    圖4 SPEI和SPI與氣候指數(shù)的交叉小波功率譜Fig. 4 Cross wavelet power spectrum of SPEI/ SPI and climate index

    總之,太原地區(qū)干旱的年際、年代際變化主要受大尺度的遙相關(guān)型影響,但干旱事件的發(fā)生及其持續(xù)是大尺度氣候因子與局地陸氣相互作用等多因子綜合作用的結(jié)果。因此,關(guān)于太原地區(qū)干旱成因尚需深入研究。

    PALMER W C. 1965. Meteorological drought US [J]. Weather Bureau Paper, 10(2):45-48.

    MCKEE T B, DOES KEN N J, KLEIST J. 1993. The relationship of drought frequency and duration to time scales [C]. Anaheim, California: Eighth Conference on Applied Climatology, January 17-22, 1993.

    VICENTE-SERRANO S M, SANTIAGO BEGUERíA, JUAN I. LOPEZ-MORENO, 2010. A Multiscale drought index sensitive to global warming: The Standardized Precipitation Evapotranspiration Index-SPEI [J]. Journal of Climate 23: 1696-1718.

    MANTUA N J, HARE S R, ZHANG Y, et al. 1997. A Pacific decadal climate oscillation with impacts on salmon [J]. Bulletin of the American Meteorological Society, 78: 1069-1079.

    RICHARD G. ALLEN, LUIS S. PEREIRA, DIRK RAES & MARTIN SMITH, 1998. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements. FAO Irrigation and drainage paper 56. Rome, Italy: Food and Agriculture Organization of the United Nations. ISBN 92-5-104219-5.

    李樹巖, 劉榮花, 師麗魁, 等. 2009. 基于CI指數(shù)的河南省近40 a干旱特征分析[J]. 干旱氣象, 27(2): 97-102.

    張調(diào)風(fēng), 張勃, 王小敏, 等. 2012. 基于綜合氣象干旱指數(shù)(CI)的干旱時空動態(tài)格局分析——以甘肅省黃土高原區(qū)為例[J]. 生態(tài)環(huán)境學(xué)報(bào), 21(1):13-20.

    曹永強(qiáng),張?zhí)m霞,張?jiān)儡?等. 2012. 基于CI指數(shù)的遼寧省氣象干旱特征分析[J]. 資源科學(xué), 34(2):265-272.

    姚玉璧,王潤元,楊金虎,等. 2011. 黃土高原陸地表層最大可能蒸散量的變化特征[J]. 生態(tài)環(huán)境學(xué)報(bào), 20(8): 1189-1195.

    王媛媛,張勃. 2012. 基于標(biāo)準(zhǔn)化降水指數(shù)的近40年隴東地區(qū)旱澇時空特征[J]. 自然資源學(xué)報(bào), 27(12): 2135-2144.

    程國棟,王根緒. 2006. 中國西北地區(qū)的干旱與旱災(zāi)——變化趨勢與對策[J]. 地學(xué)前緣, 13(1):3-14.

    黃榮輝,陳際龍,劉永. 2011. 我國東部夏季降水異常主模態(tài)的年代際變化及其與東亞水汽輸送的關(guān)系[J]. 大氣科學(xué), 35(4): 589-606.

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

    溫克剛,謝璞. 2005. 中國氣象災(zāi)害大典: 山西卷[M]. 氣象出版社.

    左金清,李維京,任宏利,等. 2012. 春季北大西洋濤動與東亞夏季風(fēng)年際關(guān)系的轉(zhuǎn)變及其可能成因分析[J]. 地球物理學(xué)報(bào), 55(2): 384-395.

    張?jiān)儡?余錦華,劉征宇,等. 2013.2009/2010年冬季中國氣溫異常及其對海表溫度的遙響應(yīng)[J]. 氣候與環(huán)境研究, 18(5): 626-638.

    李勇,何金海,姜愛軍,等. 2007. 冬季西太平洋遙相關(guān)型的環(huán)流結(jié)構(gòu)特征及其與我國冬季氣溫和降水的關(guān)系[J]. 氣象科學(xué), 27(2):119-125.

    朱益民,楊修群. 2003. 太平洋年代際振蕩與中國氣候變率的聯(lián)系[J]. 氣象學(xué)報(bào), 61(6):641-654.

    施能,朱乾根,吳彬貴. 1996. 近40年東亞夏季風(fēng)及我國夏季大尺度天氣氣候異常[J]. 大氣科學(xué), 20(5):575-583.

    Multiscale Characteristics of Drought Based on SPEI and SPI in Association with Climate Index in Taiyuan

    ZHANG Yuejun1, HAO Zhiwen1, WANG Yan, LI Mingming1, CHEN Erping, LI Fen2, ZHANG Yongrui3
    1. Shanxi Province Meteorological Science Research Institute, Taiyuan 030002, China; 2. Media center of Shanxi Province Meteorological Administration, Taiyuan 030002, China; 3. Lightning Protection and Monitoring Center of Shanxi Province, Taiyuan 030002, China

    Drought is one of the most frequent natural disasters occurred in Taiyuan on long duration, especially to the national economy caused serious influence of agricultural production. With the global warming, more frequent extreme weather occurs, Drought frequency and damage are on the rise in Taiyuan. Based on the Station data of average precipitation and temperature from 1951 to 2012, SPEI, SPI index is used to quantitatively describe the Taiyuan drought in 62 years. Based on monthly scale SPEI, SPI index to analyze the changes of Taiyuan drought, and discuss the relationship between drought and large-scale climate factors with the cross-wavelet transform (XWT). The results show that the SPEI based on precipitation and evapotranspiration can be more flexible to reflect the characteristics of the drought changes; Cross wavelet transform analysis manifested a primary resonance cycle in 6-12 a between drought in Taiyuan and four large scale factors. A 6-12 a resonance cycle between SPEI, SPI and NAO was through the test of significance in the 1985—2000, and with certain positive correlation; The relationship between SPEI, SPI and WP in 1955—1960, 1990—2000 and 1970—1990 respectively has been presented in 2-3, 4-8 and 10-16 a significant resonance cycle, and exist lag correlation in the former; The pertinence between SPEI and PDO in1986—2000 had been displayed by two significant resonance cycles of 4-6 and 8-14 a, which respectively showed that the relationship negative and positive phase. In 1955—1960, a resonance spectrum on SPEI significantly lagged behind the PDO. Compared with SPEI, the correlation between SPI and PDO is weak, the relationship of the weak resonance was only in 1955—1960 and 1960—2000 showed. Correlation between SPEI, SPI and PNA in 1983—1995 was emerged by significant resonance cycles of 4-7 a, which reflected the SPEI, SPI significantly lagged behind the PNA.

    SPEI; SPI; multiscale; drought; climate index

    P462.6

    A

    1674-5906(2014)09-1418-07

    張?jiān)儡姡轮俏?,王雁,李明明,陳二萍,李芬,張永? 基于SPEI和SPI指數(shù)的太原多尺度干旱特征與氣候

    指數(shù)的關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 2014, 23(9): 1418-1424.

    ZHANG Yuejun, HAO Zhiwen, WANG Yan, LI Mingming, CHEN Erping, LI Fen, ZHANG Yongrui. Multiscale Characteristics of Drought Based on SPEI and SPI in Association with Climate Index in Taiyuan [J]. Ecology and Environmental Sciences, 2014, 23(9): 1418-1424.

    山西省農(nóng)業(yè)科技攻關(guān)項(xiàng)目(20130312012);公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201306038);山西省氣象局青年基金(SXKQNQH20138765)

    張?jiān)儡姡?985出生),男,工程師,碩士研究生,主要從事氣候變化研究。E-mail:ccolos@163.com

    2014-06-10

    国产av精品麻豆| 亚洲一区中文字幕在线| 亚洲成a人片在线一区二区| 99国产精品一区二区蜜桃av | 伊人久久大香线蕉亚洲五| 岛国在线观看网站| 国产一区二区 视频在线| 热99re8久久精品国产| 桃红色精品国产亚洲av| 国产一区二区激情短视频| 一个人免费看片子| 12—13女人毛片做爰片一| 欧美性长视频在线观看| 窝窝影院91人妻| 午夜成年电影在线免费观看| 少妇的丰满在线观看| 亚洲人成伊人成综合网2020| 成人18禁在线播放| 人妻久久中文字幕网| 久久青草综合色| 国产aⅴ精品一区二区三区波| 精品少妇内射三级| 人妻久久中文字幕网| 亚洲色图综合在线观看| 亚洲欧美日韩高清在线视频 | 麻豆成人av在线观看| 不卡av一区二区三区| 黄色怎么调成土黄色| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人一区二区三区免费视频网站| 免费看十八禁软件| 国产精品av久久久久免费| 露出奶头的视频| 国产欧美亚洲国产| 黄片小视频在线播放| 精品一品国产午夜福利视频| 成人特级黄色片久久久久久久 | 欧美+亚洲+日韩+国产| 丰满饥渴人妻一区二区三| 国产亚洲精品一区二区www | 99热网站在线观看| 精品一区二区三区av网在线观看 | 无人区码免费观看不卡 | 久热爱精品视频在线9| 菩萨蛮人人尽说江南好唐韦庄| 一个人免费在线观看的高清视频| 久久国产亚洲av麻豆专区| 午夜福利视频在线观看免费| 久久香蕉激情| 黄色视频,在线免费观看| 男女之事视频高清在线观看| 啪啪无遮挡十八禁网站| 国产欧美日韩一区二区精品| 欧美日韩成人在线一区二区| 麻豆国产av国片精品| 国产成人免费观看mmmm| 欧美大码av| 青青草视频在线视频观看| 王馨瑶露胸无遮挡在线观看| 老司机亚洲免费影院| av免费在线观看网站| 韩国精品一区二区三区| 美女午夜性视频免费| 涩涩av久久男人的天堂| 亚洲熟女毛片儿| 国产精品亚洲一级av第二区| 国产精品一区二区精品视频观看| 99久久99久久久精品蜜桃| 亚洲国产欧美网| 久久精品亚洲精品国产色婷小说| 一区二区三区国产精品乱码| 看免费av毛片| 黄色a级毛片大全视频| 国产高清国产精品国产三级| 色综合欧美亚洲国产小说| 精品午夜福利视频在线观看一区 | 一本—道久久a久久精品蜜桃钙片| 精品午夜福利视频在线观看一区 | 夜夜骑夜夜射夜夜干| 美女福利国产在线| 国精品久久久久久国模美| 免费一级毛片在线播放高清视频 | 国产不卡一卡二| 亚洲专区国产一区二区| 老司机在亚洲福利影院| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色成人免费大全| 精品国产超薄肉色丝袜足j| 欧美黄色淫秽网站| 亚洲va日本ⅴa欧美va伊人久久| 一区二区三区激情视频| 十八禁高潮呻吟视频| 欧美日本中文国产一区发布| 久久性视频一级片| 国产精品一区二区免费欧美| 午夜福利乱码中文字幕| 久久久久久久久久久久大奶| 亚洲欧美日韩高清在线视频 | 飞空精品影院首页| 90打野战视频偷拍视频| 亚洲伊人久久精品综合| 欧美精品啪啪一区二区三区| 久久中文看片网| 大香蕉久久成人网| 欧美大码av| 国产精品1区2区在线观看. | 国产成人精品在线电影| 久久精品国产99精品国产亚洲性色 | 欧美乱码精品一区二区三区| 欧美在线黄色| 一级片免费观看大全| 欧美日韩成人在线一区二区| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久| 久久久久久免费高清国产稀缺| 亚洲成人国产一区在线观看| 操出白浆在线播放| netflix在线观看网站| 午夜精品国产一区二区电影| 丝袜美腿诱惑在线| 国产精品 欧美亚洲| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区 | 久久精品国产a三级三级三级| 成人免费观看视频高清| 桃花免费在线播放| 成年版毛片免费区| 人妻一区二区av| 亚洲色图 男人天堂 中文字幕| 精品国产乱子伦一区二区三区| 久久婷婷成人综合色麻豆| 中文字幕人妻丝袜一区二区| 99久久国产精品久久久| 久久久久网色| 日韩欧美三级三区| 99热国产这里只有精品6| 真人做人爱边吃奶动态| 精品亚洲成a人片在线观看| 99国产综合亚洲精品| 首页视频小说图片口味搜索| 国产精品一区二区精品视频观看| 天堂动漫精品| 极品人妻少妇av视频| 午夜福利影视在线免费观看| 考比视频在线观看| 香蕉国产在线看| xxxhd国产人妻xxx| 十八禁高潮呻吟视频| 国产伦人伦偷精品视频| 欧美国产精品一级二级三级| 国产成人欧美| 国产精品自产拍在线观看55亚洲 | 无限看片的www在线观看| 精品少妇一区二区三区视频日本电影| 婷婷成人精品国产| 亚洲熟妇熟女久久| 久久久久国产一级毛片高清牌| 欧美性长视频在线观看| 国产欧美日韩一区二区三| 亚洲精品自拍成人| 亚洲成人手机| 亚洲精品久久成人aⅴ小说| 麻豆av在线久日| 天堂8中文在线网| 18禁观看日本| 天天操日日干夜夜撸| 久久ye,这里只有精品| 999久久久精品免费观看国产| 中文字幕人妻丝袜一区二区| 老司机午夜福利在线观看视频 | 日韩欧美一区二区三区在线观看 | 久久久久久亚洲精品国产蜜桃av| 又紧又爽又黄一区二区| 一个人免费在线观看的高清视频| 2018国产大陆天天弄谢| 老司机靠b影院| 90打野战视频偷拍视频| 国产片内射在线| 在线观看www视频免费| 国产精品九九99| 成人国语在线视频| 久热爱精品视频在线9| 十八禁网站网址无遮挡| 91麻豆av在线| 美女视频免费永久观看网站| 精品国产一区二区三区久久久樱花| 色综合欧美亚洲国产小说| 亚洲精品成人av观看孕妇| 久久久久精品人妻al黑| 亚洲综合色网址| 中文亚洲av片在线观看爽 | 99热国产这里只有精品6| 老司机深夜福利视频在线观看| 久热这里只有精品99| 久久九九热精品免费| 真人做人爱边吃奶动态| 欧美午夜高清在线| 女人精品久久久久毛片| 亚洲精品在线美女| 18禁美女被吸乳视频| 欧美精品av麻豆av| 欧美性长视频在线观看| 夫妻午夜视频| 午夜精品国产一区二区电影| 国产成人精品无人区| 青青草视频在线视频观看| av线在线观看网站| 国产精品 国内视频| 国产精品电影一区二区三区 | 亚洲,欧美精品.| 午夜老司机福利片| 在线观看免费视频日本深夜| 男女之事视频高清在线观看| 久久青草综合色| 丝袜喷水一区| 美女主播在线视频| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 国产aⅴ精品一区二区三区波| 大型黄色视频在线免费观看| 精品国产超薄肉色丝袜足j| 亚洲精品av麻豆狂野| 免费看a级黄色片| av超薄肉色丝袜交足视频| 免费在线观看日本一区| 久热这里只有精品99| 国产欧美日韩综合在线一区二区| 手机成人av网站| 欧美黑人欧美精品刺激| 老司机福利观看| 婷婷成人精品国产| 麻豆av在线久日| 精品卡一卡二卡四卡免费| 人妻 亚洲 视频| 可以免费在线观看a视频的电影网站| 九色亚洲精品在线播放| 国内毛片毛片毛片毛片毛片| 免费在线观看日本一区| 久久天堂一区二区三区四区| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 纵有疾风起免费观看全集完整版| 欧美日本中文国产一区发布| 国产人伦9x9x在线观看| 99re在线观看精品视频| 国产成人av教育| 亚洲午夜精品一区,二区,三区| 欧美老熟妇乱子伦牲交| 2018国产大陆天天弄谢| 国产日韩欧美亚洲二区| 欧美人与性动交α欧美精品济南到| av不卡在线播放| 91精品三级在线观看| 一本色道久久久久久精品综合| 激情视频va一区二区三区| 精品国产乱码久久久久久男人| 欧美大码av| 欧美av亚洲av综合av国产av| 久久这里只有精品19| 一二三四在线观看免费中文在| 国产三级黄色录像| 亚洲情色 制服丝袜| 丁香欧美五月| 精品国产一区二区久久| 999久久久精品免费观看国产| 精品免费久久久久久久清纯 | 成人永久免费在线观看视频 | 天天躁日日躁夜夜躁夜夜| 欧美老熟妇乱子伦牲交| 亚洲欧洲日产国产| 亚洲第一av免费看| 久久这里只有精品19| 黄网站色视频无遮挡免费观看| 亚洲一区二区三区欧美精品| 王馨瑶露胸无遮挡在线观看| 精品人妻熟女毛片av久久网站| 成人三级做爰电影| 99国产极品粉嫩在线观看| 国产成人影院久久av| 香蕉久久夜色| 成人免费观看视频高清| 午夜免费鲁丝| 欧美老熟妇乱子伦牲交| 免费av中文字幕在线| 国产精品影院久久| 美女主播在线视频| 999久久久精品免费观看国产| 高清黄色对白视频在线免费看| 少妇猛男粗大的猛烈进出视频| 久久精品熟女亚洲av麻豆精品| 亚洲国产av新网站| 国产欧美日韩一区二区三区在线| 国产色视频综合| 岛国在线观看网站| 久久精品人人爽人人爽视色| 一级毛片女人18水好多| 少妇被粗大的猛进出69影院| 一级片'在线观看视频| 丝瓜视频免费看黄片| 91国产中文字幕| 老司机亚洲免费影院| 亚洲av日韩精品久久久久久密| 久久久久精品国产欧美久久久| 午夜福利视频在线观看免费| 亚洲人成伊人成综合网2020| 中文字幕制服av| 一边摸一边抽搐一进一小说 | 日韩成人在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 国产视频一区二区在线看| 男女之事视频高清在线观看| 日韩免费高清中文字幕av| 欧美日本中文国产一区发布| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 国产区一区二久久| 乱人伦中国视频| 欧美日韩成人在线一区二区| 国产视频一区二区在线看| 亚洲七黄色美女视频| 男男h啪啪无遮挡| 一二三四社区在线视频社区8| 黑人猛操日本美女一级片| 新久久久久国产一级毛片| 亚洲人成电影免费在线| 黄色 视频免费看| 视频在线观看一区二区三区| 亚洲专区国产一区二区| 一级,二级,三级黄色视频| 国产精品熟女久久久久浪| av有码第一页| 真人做人爱边吃奶动态| 亚洲中文av在线| 色视频在线一区二区三区| 欧美变态另类bdsm刘玥| 国产精品欧美亚洲77777| 少妇精品久久久久久久| 成人国产av品久久久| 国产成人影院久久av| 手机成人av网站| 黄频高清免费视频| 一本大道久久a久久精品| 国产亚洲欧美在线一区二区| 在线观看免费高清a一片| 国产精品二区激情视频| 最近最新中文字幕大全电影3 | 国产欧美日韩一区二区三| 亚洲免费av在线视频| 伊人久久大香线蕉亚洲五| 一区在线观看完整版| 亚洲中文字幕日韩| 两性夫妻黄色片| 中文亚洲av片在线观看爽 | 黑人猛操日本美女一级片| 精品一区二区三卡| 久久久久国产一级毛片高清牌| 免费黄频网站在线观看国产| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区综合在线观看| 久久99热这里只频精品6学生| 国产伦理片在线播放av一区| 人人澡人人妻人| 人妻一区二区av| 精品第一国产精品| 国产伦理片在线播放av一区| 日本av手机在线免费观看| 国产精品二区激情视频| 国产精品电影一区二区三区 | 操出白浆在线播放| 国产深夜福利视频在线观看| 亚洲成a人片在线一区二区| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 18禁观看日本| 色婷婷av一区二区三区视频| 亚洲美女黄片视频| 婷婷丁香在线五月| 国产淫语在线视频| 黄色丝袜av网址大全| 一区二区日韩欧美中文字幕| 丰满人妻熟妇乱又伦精品不卡| 久久精品aⅴ一区二区三区四区| 国产精品免费大片| 91麻豆精品激情在线观看国产 | 精品午夜福利视频在线观看一区 | 亚洲国产av新网站| 久久久国产欧美日韩av| 人妻久久中文字幕网| 99re在线观看精品视频| www.熟女人妻精品国产| 久久av网站| 国产精品一区二区免费欧美| 欧美激情极品国产一区二区三区| 亚洲成人手机| 免费观看av网站的网址| 12—13女人毛片做爰片一| 制服诱惑二区| 亚洲国产成人一精品久久久| 日本撒尿小便嘘嘘汇集6| 国产精品影院久久| 中文字幕色久视频| avwww免费| av网站在线播放免费| tube8黄色片| 日韩中文字幕视频在线看片| 两个人免费观看高清视频| 色尼玛亚洲综合影院| 国产不卡av网站在线观看| 制服人妻中文乱码| 十八禁网站免费在线| 啦啦啦中文免费视频观看日本| 美女高潮到喷水免费观看| 一区二区三区国产精品乱码| 少妇裸体淫交视频免费看高清 | 中国美女看黄片| 女性被躁到高潮视频| 男女之事视频高清在线观看| 久久久久网色| 丝袜人妻中文字幕| 国产视频一区二区在线看| 亚洲人成77777在线视频| 伦理电影免费视频| 男女床上黄色一级片免费看| 久久99热这里只频精品6学生| 久久久精品免费免费高清| 极品人妻少妇av视频| 久久久国产欧美日韩av| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| 中文字幕最新亚洲高清| 在线天堂中文资源库| 考比视频在线观看| 久久久久视频综合| 两个人免费观看高清视频| 欧美大码av| 精品人妻熟女毛片av久久网站| av免费在线观看网站| 国产伦人伦偷精品视频| 久久久久网色| 一区二区三区精品91| 日韩一区二区三区影片| av福利片在线| 丁香六月天网| 欧美一级毛片孕妇| 久久毛片免费看一区二区三区| 丰满迷人的少妇在线观看| 99在线人妻在线中文字幕 | 手机成人av网站| 丝袜人妻中文字幕| 三上悠亚av全集在线观看| 亚洲国产av影院在线观看| 免费黄频网站在线观看国产| 成人18禁在线播放| 91精品三级在线观看| 午夜日韩欧美国产| 纵有疾风起免费观看全集完整版| 动漫黄色视频在线观看| 亚洲伊人久久精品综合| 精品一品国产午夜福利视频| 国产一区二区 视频在线| 国产三级黄色录像| 亚洲va日本ⅴa欧美va伊人久久| 男女边摸边吃奶| 妹子高潮喷水视频| 亚洲成a人片在线一区二区| 亚洲国产看品久久| 国产精品一区二区免费欧美| 国产男女超爽视频在线观看| 久久久久网色| 香蕉国产在线看| av线在线观看网站| 欧美中文综合在线视频| 十八禁网站免费在线| 亚洲av成人一区二区三| 69av精品久久久久久 | 人人妻人人澡人人看| 丁香欧美五月| 欧美久久黑人一区二区| 国产精品一区二区在线观看99| 成人永久免费在线观看视频 | 亚洲专区中文字幕在线| 自拍欧美九色日韩亚洲蝌蚪91| 热99久久久久精品小说推荐| 真人做人爱边吃奶动态| cao死你这个sao货| 亚洲国产欧美日韩在线播放| 夜夜夜夜夜久久久久| 成年动漫av网址| 黑人欧美特级aaaaaa片| 久久中文字幕人妻熟女| 亚洲欧美日韩另类电影网站| 亚洲人成电影免费在线| 99re在线观看精品视频| 日本av手机在线免费观看| 久久中文字幕一级| 日韩一区二区三区影片| aaaaa片日本免费| 一区福利在线观看| 精品国产乱子伦一区二区三区| 精品久久久久久久毛片微露脸| 老熟妇乱子伦视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品麻豆人妻色哟哟久久| 飞空精品影院首页| 在线观看免费视频网站a站| 久久ye,这里只有精品| 欧美日韩中文字幕国产精品一区二区三区 | 无限看片的www在线观看| 亚洲av国产av综合av卡| 狠狠精品人妻久久久久久综合| 老司机午夜福利在线观看视频 | 欧美另类亚洲清纯唯美| 亚洲精品美女久久av网站| 国产成+人综合+亚洲专区| 丝袜人妻中文字幕| 午夜免费成人在线视频| 丰满人妻熟妇乱又伦精品不卡| 美女高潮到喷水免费观看| 欧美激情久久久久久爽电影 | 无遮挡黄片免费观看| 国产片内射在线| 亚洲av日韩精品久久久久久密| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 91成人精品电影| 成人影院久久| 俄罗斯特黄特色一大片| 在线观看免费视频日本深夜| 日本撒尿小便嘘嘘汇集6| 他把我摸到了高潮在线观看 | 一级片免费观看大全| 日韩免费av在线播放| 人妻久久中文字幕网| 老汉色av国产亚洲站长工具| www日本在线高清视频| 成人黄色视频免费在线看| 波多野结衣av一区二区av| av天堂久久9| 人人妻人人澡人人爽人人夜夜| 黄色视频在线播放观看不卡| 免费在线观看黄色视频的| avwww免费| 国产1区2区3区精品| 久久久久久久精品吃奶| 欧美 日韩 精品 国产| 国产亚洲欧美精品永久| 国产在线免费精品| 97人妻天天添夜夜摸| 无人区码免费观看不卡 | 久久99热这里只频精品6学生| 亚洲精品中文字幕一二三四区 | 可以免费在线观看a视频的电影网站| 免费在线观看完整版高清| 女人爽到高潮嗷嗷叫在线视频| a在线观看视频网站| 国产日韩欧美亚洲二区| 久久人人爽av亚洲精品天堂| 狠狠婷婷综合久久久久久88av| 动漫黄色视频在线观看| 大码成人一级视频| 久久久水蜜桃国产精品网| 在线 av 中文字幕| 日本一区二区免费在线视频| 日韩欧美一区视频在线观看| 99国产精品99久久久久| 一级,二级,三级黄色视频| 精品人妻1区二区| 国产日韩一区二区三区精品不卡| 国产伦人伦偷精品视频| 欧美 日韩 精品 国产| 免费久久久久久久精品成人欧美视频| 国产在线免费精品| 少妇 在线观看| 最新在线观看一区二区三区| 99国产精品免费福利视频| 亚洲中文av在线| 我要看黄色一级片免费的| 欧美成狂野欧美在线观看| 日韩视频在线欧美| 真人做人爱边吃奶动态| 国产黄色免费在线视频| 国产日韩欧美亚洲二区| 午夜免费鲁丝| 久久久久久久久久久久大奶| av天堂在线播放| 精品一区二区三区av网在线观看 | av国产精品久久久久影院| 欧美在线一区亚洲| 国产深夜福利视频在线观看| 国产aⅴ精品一区二区三区波| 久久人妻熟女aⅴ| 日韩免费av在线播放| 国产高清videossex| 精品国产一区二区三区久久久樱花| 热re99久久精品国产66热6| 精品久久蜜臀av无| 在线天堂中文资源库| 久久国产精品影院| e午夜精品久久久久久久| 99精品欧美一区二区三区四区| 在线十欧美十亚洲十日本专区| 女警被强在线播放| a级片在线免费高清观看视频| 老司机深夜福利视频在线观看| 亚洲人成电影免费在线| 国产精品国产av在线观看| 一本久久精品| 精品熟女少妇八av免费久了| 日韩欧美一区视频在线观看| 一区二区三区乱码不卡18| 欧美人与性动交α欧美精品济南到| 80岁老熟妇乱子伦牲交| 高清欧美精品videossex| 亚洲一卡2卡3卡4卡5卡精品中文| 在线永久观看黄色视频|