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

    黃河流域干旱時(shí)空特征及其與ENSO的關(guān)聯(lián)性分析

    2021-12-08 01:24:56黃婷婷林青霞吳志勇汪瑛琪
    人民黃河 2021年11期
    關(guān)鍵詞:干旱黃河流域

    黃婷婷 林青霞 吳志勇 汪瑛琪

    摘 要:為探究黃河流域歷史氣象干旱時(shí)空演變特征,并結(jié)合ENSO事件探究其關(guān)聯(lián)性,選取黃河流域99個(gè)氣象站1961—2020年的逐月降水與氣溫資料,結(jié)合Mann-Kendall突變檢驗(yàn)、Morlet小波分析等方法揭示近60 a黃河流域年、季尺度的干旱時(shí)空特征,并結(jié)合平均海平面溫度距平指數(shù)SSTA,探究黃河流域氣象干旱特征與ENSO事件的相關(guān)性。研究結(jié)果表明:①年尺度上,黃河流域近60 a的標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI以-0.114/10 a的速度變化,干旱現(xiàn)象較嚴(yán)重,空間上,流域西部呈輕微濕潤化趨勢(shì),流域東部呈干旱化趨勢(shì);②季尺度上,全流域春季干旱化趨勢(shì)最為顯著,夏秋季其次,冬季干旱化趨勢(shì)微弱;③Morlet小波分析顯示,黃河流域存在13 a主周期的旱澇交替,2020年的小波系數(shù)實(shí)部圖顯示,黃河流域?qū)⒊掷m(xù)濕潤;④ENSO事件發(fā)生強(qiáng)度與黃河流域SPEI影響因子的相關(guān)性表明,ENSO事件的發(fā)生與降水及溫度的相關(guān)性顯著,且相關(guān)性存在1~6個(gè)月的滯后。

    關(guān)鍵詞:干旱;標(biāo)準(zhǔn)化降水蒸散指數(shù);ENSO事件;黃河流域

    中圖分類號(hào):P426.616; TV882.1

    文獻(xiàn)標(biāo)志碼:A

    doi:10.3969/j.issn.1000-1379.2021.11.010

    引用格式:黃婷婷,林青霞,吳志勇,等.黃河流域干旱時(shí)空特征及其與ENSO的關(guān)聯(lián)性分析[J].人民黃河,2021,43(11):52-58.

    Spatial and Temporal Characteristics of Drought in the Yellow River Basin and Their Correlation with ENSO

    HUANG Tingting1, LIN Qingxia2, WU Zhiyong ?WANG Yingqi3

    (1.Yellow River Institute of Hydraulic Research, Zhengzhou 45000 ?China;

    2.College of Water Conservancy and Environment, China Three Gorges University, Yichang 443002, China;

    3.College of Hydrology and Water Resources, Hohai University, Nanjing 210098, China)

    Abstract: In order to explore the spatiotemporal evolution characteristics of historical meteorological drought in theYellow River Basin and explore its correlation with ENSO events, it selected the monthly precipitation and temperature data of 99 meteorological stations in the Yellow River Basin from 1961 to 2020. Based on Mann-Kendall mutation test and Morlet wavelet analysis, it revealed the temporal and spatial characteristics of annual and seasonal drought in the Yellow River Basin over the past 60 years, and the correlation between meteorological drought characteristics in the Yellow River Basin and ENSO events were explored by combining the mean sea surface temperature anomaly index (SSTA). The results show that a) on the annual scale, the SPEI index of the Yellow River Basin during recent 60 years is decreased at a rate of -0.114/10 a, the drought phenomenon is serious. Spatially, the western part of the basin shows a slight trend of wetness, while the eastern part of the basin shows a trend of aridity. b) At the seasonal scale, the most significant drying trend is spring, followed by summer and autumn, and the drying trend is weak in winter. c) Morlet wavelet analysis shows that there is a 13-year main cycle of drought and flood alternation in the Yellow River Basin. By observing the real part of wavelet coefficients in 2020, it shows that the Yellow River Basin will remain wet. d) The polynomial fitting relationship and correlation between the occurrence intensity of ENSO events and the SPEI influencing factors in the Yellow River Basin indicate that the occurrence of ENSO events is significantly correlated with precipitation and temperature and the correlation has a lag of 1-6 months.

    Key words: drought; SPEI; ENSO events; Yellow River Basin

    干旱是普遍發(fā)生的一種氣象災(zāi)害,是指區(qū)域水分收支或供求不平衡而形成的水分短缺現(xiàn)象。IPCC評(píng)估報(bào)告顯示,近百年來全球地表溫度升高明顯,全球變暖使得全球各地的旱災(zāi)事件呈顯著增加趨勢(shì)。干旱災(zāi)害的頻發(fā)使得我國多地水資源、水生態(tài)環(huán)境以及農(nóng)業(yè)生產(chǎn)等受到嚴(yán)重影響[1-2]。描述干旱特征的干旱指標(biāo)眾多,包括標(biāo)準(zhǔn)化降水指數(shù)SPI[3]、帕默爾干旱指數(shù)PDSI[4]、Z指數(shù)[5]、標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI等。SPI與Z以降水量為單因素表征干旱特征,未考慮溫度、蒸散發(fā)等因素對(duì)區(qū)域氣候變化特征的影響。PDSI的計(jì)算過程較復(fù)雜,對(duì)研究數(shù)據(jù)的要求較高,且干旱等級(jí)劃分的主觀性較大。SPEI的計(jì)算考慮了溫度與水量平衡因子,在干旱評(píng)估中具有優(yōu)勢(shì),被廣泛應(yīng)用與研究[6]。

    黃河流域是我國重要的生態(tài)屏障和經(jīng)濟(jì)帶,在我國生態(tài)建設(shè)與社會(huì)經(jīng)濟(jì)發(fā)展等方面具有不可替代的地位[7]。21世紀(jì)以來,黃河流域干旱事件頻發(fā)[8],干旱導(dǎo)致的河流干枯和斷流現(xiàn)象日益嚴(yán)重,加速了黃河流域土地荒漠化與植被退化[9]。因此,研究黃河流域干旱特征與演變規(guī)律有利于更加準(zhǔn)確地預(yù)測未來黃河流域氣候變化以及水資源時(shí)空變化規(guī)律[10]。已有學(xué)者對(duì)黃河流域的干旱特征進(jìn)行了研究與分析,如:楊肖麗等[11]利用標(biāo)準(zhǔn)化降水指數(shù)SPI對(duì)黃河流域氣象干旱進(jìn)行了評(píng)估與分析;周帥等[12]以黃河流域?yàn)檠芯繀^(qū)域,運(yùn)用游程理論、Morans I指數(shù)、SPI指數(shù)等方法綜合分析了黃河流域干旱時(shí)空變化規(guī)律,并識(shí)別出干旱事件發(fā)生的位置及空間分布格局。但是,上述研究多注重于使用標(biāo)準(zhǔn)化降水指數(shù)SPI分析區(qū)域干旱特征,較少涉及標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI的應(yīng)用。筆者采用SPEI作為研究干旱的指標(biāo),分析近60 a來黃河流域年代際與四季的干旱時(shí)空變化特征,并探究ENSO事件[13]與干旱特征的相關(guān)性,旨在加快推動(dòng)黃河流域高質(zhì)量發(fā)展[14-15]。

    1 研究區(qū)域概況

    黃河發(fā)源于青海省巴顏喀拉山北麓,流經(jīng)青海、四川、甘肅、寧夏、內(nèi)蒙古、山西、陜西、河南、山東9個(gè)省(區(qū)),干流全長5 464 km,流域面積79.5萬km2,地理位置為北緯32°10′— 41°50′、東經(jīng)95°53′—119°5′,黃河流域海拔由西向東逐級(jí)下降,流域降水量東南多、西北少,多年平均降水量約為460 mm,多年平均水面蒸散發(fā)量約為1 000 mm,是我國典型的氣候敏感區(qū)之一。黃河流域水資源相對(duì)匱乏,生態(tài)環(huán)境脆弱,在全球性氣候變暖與人類活動(dòng)影響下,黃河流域的干旱災(zāi)害頻繁發(fā)生,根據(jù)黃河流域歷史重大干旱事件記載[16-17],21世紀(jì)初期黃河流域干旱事件影響范圍較廣且影響歷時(shí)較長。

    2 研究數(shù)據(jù)與方法

    2.1 研究數(shù)據(jù)

    本文選取了黃河流域99個(gè)氣象站1961—2020年的氣象數(shù)據(jù),缺測數(shù)據(jù)采用相鄰站點(diǎn)數(shù)據(jù)插值的方法獲取,氣象數(shù)據(jù)來源于中國氣象數(shù)據(jù)共享服務(wù)中心(http://cdc.nmic.cn/),包括逐月平均溫度和降水量觀測數(shù)據(jù)。選用熱帶太平洋平面溫度距平指數(shù)SSTA[18]與南方濤動(dòng)指數(shù)SOI[19]作為表示厄爾尼諾-南方濤動(dòng)的指標(biāo),SSTA指數(shù)來源于美國地球系統(tǒng)研究實(shí)驗(yàn)室(NOAA-ESRL)。

    2.2 研究方法

    2.2.1 標(biāo)準(zhǔn)化降水蒸散指數(shù)

    基于標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI分析黃河流域近60 a的干旱變化規(guī)律,應(yīng)用Thomthwaite方法計(jì)算逐月潛在蒸散量PETi。計(jì)算SPEI涉及的主要指標(biāo)及公式如下[20]:

    PETi=16×10TiHm(1)

    Di=Pi-PETi(2)

    F(x)=1+αx-ββ-1(3)

    α=ω0-2ω1βΓ1+1βΓ1-1β(4)

    β=2ω1-ω06ω1-ω0-6ω2(5)

    γ=ω0-αΓ1+1βΓ1-1β(6)

    SPEI=ω-c0+c1ω+c2ω21+d1ω+d2ω2+d3ω3(7)

    ω=-2ln P(8)

    式中:PETi為月潛在蒸散量(i為月份);Pi為月降水量;Ti為月平均溫度;H為年熱量指數(shù);m為常數(shù),m=0.49+0.179H-0.000 077 1H2+0.000 000 675H3;α為尺度函數(shù);β為形狀函數(shù);γ為origin參數(shù);P為累計(jì)概率密度;F(x)為log-logistic概率分布的累計(jì)函數(shù);Γ為階乘函數(shù);ω0、ω1、ω2為數(shù)據(jù)序列Di的概率加權(quán)矩;c0、c1、c2、d1、d2、d3均為常數(shù);ω為距離加權(quán)矩。若P≤0.5,則P=F(x);若P>0.5,則P=1-F(x)。

    本文通過計(jì)算研究區(qū)域99個(gè)氣象站點(diǎn)1961年1月至2020年12月不同時(shí)間尺度的SPEI值,分析研究區(qū)域的年度與季節(jié)的干旱時(shí)空演變特征,SPEI干旱等級(jí)劃分標(biāo)準(zhǔn)[21]見表1。

    SPEI具有多時(shí)間尺度,本文主要分析3個(gè)月與12個(gè)月時(shí)間尺度的SPEI值。文中的春旱、夏旱、秋旱和冬旱分別采用5月、8月、11月與次年2月的SPEI-3表征,年尺度干旱特征采用SPEI-12表征。

    2.2.2 Mann-Kendall檢驗(yàn)

    Mann-Kendall(M-K)趨勢(shì)檢驗(yàn)是目前氣象水文領(lǐng)域進(jìn)行趨勢(shì)判斷的主要方法之一[22],假設(shè)時(shí)間序列數(shù)據(jù){Xi}是n個(gè)獨(dú)立的、隨機(jī)變量同分布的樣本,則:

    S=∑ni=2∑i-1j=1sgnXi-Xj(9)

    式中:S為統(tǒng)計(jì)量。

    依據(jù)檢驗(yàn)統(tǒng)計(jì)量Z判斷序列的趨勢(shì)顯著性,其計(jì)算公式為

    Z=S-1 nn-12n+5/18(S>0)

    0(S=0)

    S+1 nn-12n+5/18(S<0)(10)

    若檢驗(yàn)統(tǒng)計(jì)量Z的絕對(duì)值Z>2.58,表明序列在0.01的置信水平上具有極顯著的變化趨勢(shì);若1.96

    2.2.3 小波分析

    在Matlab軟件中采用復(fù)Morlet小波分析,實(shí)現(xiàn)黃河流域年尺度與季尺度的SPEI指數(shù)周期變化規(guī)律研究。小波分析可反映黃河流域SPEI指數(shù)的多尺度周期特征,通過小波變換系數(shù)和小波方差圖可以確定區(qū)域干旱的主要周期以及各時(shí)間尺度的波動(dòng)能量隨該尺度的分布情況[14]。

    3 結(jié)果分析

    3.1 黃河流域SPEI的年際變化趨勢(shì)分析

    黃河流域1961—2020年年度SPEI變化及其M-K檢驗(yàn)曲線如圖1(a)所示,可以看出年度SPEI呈下降趨勢(shì),線性傾向率為-0.114/10 a,表明近60 a來黃河流域呈干旱化趨勢(shì)。1964年和1997年分別為黃河流域近60 a來最濕潤和最干旱的兩個(gè)年份,其SPEI值分別為1.47和-1.09,分別達(dá)到了中度濕潤和中度干旱的水平。在置信區(qū)間內(nèi),UF和UB曲線相交于1985年,表明1985年黃河流域干旱情勢(shì)發(fā)生突變,1985年前流域干旱化趨勢(shì)不顯著,1985年后流域顯著干旱化,且1999年后干旱化趨勢(shì)達(dá)到了顯著性水平。

    圖1(b)~(e)為黃河流域1961—2020年四季SPEI變化及其M-K檢驗(yàn)曲線。春季SPEI的線性傾向率為-0.117/10 a,1964年和2000年分別是最濕潤和最干旱的兩個(gè)年份,其SPEI值分別為1.22和-1.40,分別達(dá)到中度濕潤和中度干旱的水平。夏季SPEI的線性傾向率為-0.047/10 a,1976年和1997年夏季SPEI分別為1.01和-1.04,分別達(dá)到中度濕潤和中度干旱的水平。秋季SPEI的線性傾向率為-0.101/10 a,1961年和1998年秋季SPEI分別為1.38和-1.80,分別達(dá)到中度濕潤和重度干旱的水平。冬季SPEI的線性傾向率為0.068/10 a,1990年和1999年的冬季SPEI分別為1.34和-1.59,分別達(dá)到中度濕潤和重度干旱的水平。綜上可以看出,近60 a來黃河流域的春季、夏季、秋季均表現(xiàn)出由濕潤向干旱轉(zhuǎn)變的趨勢(shì),而冬季則表現(xiàn)出濕潤化的趨勢(shì)。

    3.2 黃河流域SPEI的空間變化趨勢(shì)分析

    圖2為黃河流域1961—2020年年度與四季SPEI空間分布狀況。從年度SPEI空間分布看,以景泰—臨夏—若爾蓋—馬爾康為分界線把流域分為東部與西部地區(qū),西部地區(qū)呈現(xiàn)濕潤化趨勢(shì),東部地區(qū)呈現(xiàn)中度至極端干旱水平的干旱化趨勢(shì)??傮w而言,1961—2020年流域大部分區(qū)域呈現(xiàn)干旱化加重的趨勢(shì),小部分地區(qū)呈現(xiàn)濕潤化加重的趨勢(shì),流域上中游地區(qū)干旱化程度加重,流域上游地區(qū)濕潤化程度加重。

    從春季SPEI空間分布情況看,以景泰—西吉—天水與太原—臨汾—孟津?yàn)閮蓷l分界線把流域分為西部、中部、東部三部分,西部地區(qū)呈濕潤化趨勢(shì),中部地區(qū)呈干旱化趨勢(shì),東部地區(qū)呈輕微干旱化趨勢(shì);從夏季SPEI空間分布看,大部分地區(qū)呈中度干旱化趨勢(shì),只有西北與東南部分地區(qū)呈濕潤化趨勢(shì);從秋季SPEI空間分布看,大部分地區(qū)呈極端干旱化狀態(tài),以景泰—合作—若爾蓋與銀川—橫山—長治為兩條分界線將流域分為三部分,除西南與東北地區(qū)呈輕微濕潤化趨勢(shì)外,其余地區(qū)均具有干旱化趨勢(shì)嚴(yán)重的特征;從冬季SPEI空間分布看,以包頭—橫山—環(huán)縣—長武—天水為分界線把流域分為東西兩部分,中游與下游具有干旱化的趨勢(shì)。整體而言,流域中游地區(qū)容易發(fā)生干旱事件,上游地區(qū)相對(duì)不易發(fā)生干旱事件。春季干旱化最為顯著,夏季與秋季次之,冬季干旱化趨勢(shì)最不顯著。

    3.3 Morlet小波分析

    選用水文氣象分析中常用的Morlet連續(xù)復(fù)小波變換分析法分析黃河流域1961—2020年的年尺度與季節(jié)尺度SPEI的周期演變規(guī)律,見圖3。

    年度SPEI既存在4~8 a與10~15 a的短周期變化規(guī)律,也存在18~27 a的長周期變化規(guī)律。其中,10~15 a的短周期存在3次較明顯的旱濕交替情況,18~27 a的長周期存在2次較明顯的旱濕交替情況。13、5、21 a尺度上震蕩明顯,分別為第一、第二、第三周期;春季SPEI存在10、36、18 a的震蕩周期,分別為第一、第二、第三周期;夏季SPEI存在21、7、13 a的震蕩周期,分別為第一、第二、第三周期,13 a與21 a尺度上流域均處于持續(xù)偏濕的狀態(tài);秋季13 a為第一周期,5 a為第二周期,13 a尺度上未來幾年流域處于持續(xù)干旱的狀態(tài),5 a尺度上流域處于持續(xù)濕潤的狀態(tài);冬季23 a為第一周期,7 a為第二周期,23 a尺度上流域處于持續(xù)濕潤的狀態(tài),7 a尺度上流域處于濕潤結(jié)束期。

    3.4 ENSO事件對(duì)黃河流域干旱的影響

    3.4.1 ENSO事件強(qiáng)度與SPEI指數(shù)影響因子關(guān)系

    ENSO事件主要分為暖事件(厄爾尼諾事件)和冷事件(拉尼娜事件)。對(duì)1961—2020年的ENSO冷暖事件及強(qiáng)度等級(jí)的變化特征進(jìn)行統(tǒng)計(jì),見表2。李曉燕等[23]將ENSO暖冷事件強(qiáng)度的等級(jí)劃分為強(qiáng)、中等、弱3個(gè)等級(jí),其中ENSO暖事件強(qiáng)度分別用3、2、1表示,ENSO冷事件強(qiáng)度分別用-3、-2、-1表示,未受ENSO事件影響的年份強(qiáng)度用0表示。黃河流域1961—2020年的SPEI值、年降水量、年均氣溫與ENSO事件強(qiáng)度的關(guān)系如圖4~圖6所示。

    由圖4可知,ENSO事件的強(qiáng)度與黃河流域SPEI指數(shù)存在顯著負(fù)相關(guān)關(guān)系,經(jīng)相關(guān)分析得出,相關(guān)系數(shù)為-0.33 達(dá)到極顯著性水平。由圖5可知,ENSO事件對(duì)黃河流域降水的影響較強(qiáng),兩者成明顯負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)為-0.270,達(dá)到顯著性水平,ENSO強(qiáng)度增大時(shí)降水量下降,減弱時(shí)降水量上升。由圖6可知,ENSO強(qiáng)度與氣溫之間整體成較明顯的正相關(guān)關(guān)系,相關(guān)系數(shù)為0.121,沒有達(dá)到顯著性水平??梢钥闯觯邓c氣溫受ENSO事件的影響存在一定的時(shí)間滯后。

    圖7為1961—2020年的SPEI-12指數(shù)、SSTA指數(shù)、SOI指數(shù)演變示意圖。可以看出,由厄爾尼諾事件引起的干旱事件,其SSTA指數(shù)與SOI指數(shù)分別對(duì)應(yīng)著相同時(shí)段的SPEI-12指數(shù),且Pearson相關(guān)性系數(shù)分別為-0.306與0.188,均達(dá)到了極顯著性水平,由此可見厄爾尼諾事件的發(fā)生對(duì)黃河流域的干旱特征具有非常明顯的影響。

    3.4.2 ENSO事件對(duì)降水量與氣溫影響的滯后性

    為探討ENSO事件對(duì)區(qū)域降水量與氣溫影響的滯后性,利用與SSTA同期和滯后1~6個(gè)月的月降水量數(shù)據(jù)與月均氣溫?cái)?shù)據(jù)分別與SSTA指數(shù)進(jìn)行相關(guān)性分析,結(jié)果見表3和表4。由表3可知,ENSO暖事件發(fā)生時(shí),同期SSTA指數(shù)與滯后1個(gè)月降水量的相關(guān)性逐漸增強(qiáng),滯后2個(gè)月后相關(guān)性減弱,滯后4~6個(gè)月負(fù)相關(guān)性不顯著。ENSO冷事件發(fā)生時(shí),同期SSTA指數(shù)與同期降水量正相關(guān)性不顯著,滯后1~2個(gè)月正相關(guān)性逐漸增強(qiáng),滯后4~6個(gè)月負(fù)相關(guān)性不顯著。由此可見,ENSO暖事件的發(fā)生對(duì)同期至滯后3個(gè)月的降水量有明顯的減少效應(yīng),ENSO冷事件的發(fā)生對(duì)同期至滯后3個(gè)月的降水量有明顯的增加效應(yīng)。由表4可知,ENSO暖事件發(fā)生時(shí),同期SSTA與滯后1個(gè)月氣溫的負(fù)相關(guān)性逐漸增強(qiáng),與滯后2~3個(gè)月負(fù)相關(guān)性減弱,與滯后4~6個(gè)月正相關(guān)性逐漸增強(qiáng)。ENSO冷事件發(fā)生時(shí),同期SSTA與滯后1個(gè)月氣溫的正相關(guān)性逐漸增強(qiáng),與滯后2~3個(gè)月正相關(guān)性減弱,與滯后4~6個(gè)月負(fù)相關(guān)性逐漸增強(qiáng)。

    4 結(jié) 論

    以黃河流域?yàn)檠芯繉?duì)象,利用黃河流域99個(gè)代表性氣象站1961—2020年氣象數(shù)據(jù),以標(biāo)準(zhǔn)化降水蒸散指數(shù)SPEI為干旱指標(biāo),運(yùn)用Mann-Kendall趨勢(shì)檢驗(yàn)法、復(fù)Morlet小波分析法等對(duì)年度及季節(jié)SPEI的時(shí)空變化規(guī)律和黃河流域干旱時(shí)空格局進(jìn)行了分析。最后,利用Pearson相關(guān)性探究了ENSO事件強(qiáng)度及海平面溫度距平指數(shù)SSTA分別與SPEI及SPEI影響因子的聯(lián)系,得出以下結(jié)論。

    (1)近60 a來黃河流域呈干旱化加劇的趨勢(shì),2000年后干旱趨勢(shì)達(dá)到顯著性水平。整體而言,黃河流域春季干旱化趨勢(shì)最為明顯,秋季與夏季次之,冬季最弱。由于黃河流域春季與秋季降水量較少,年際變化大,夏季農(nóng)作物需水量較大且伴隨氣溫升高,因此流域春、夏、秋季干旱情勢(shì)最為嚴(yán)峻。流域冬季易受ENSO冷事件的影響呈現(xiàn)濕潤化趨勢(shì)。

    (2)Moelet小波分析顯示黃河流域在13 a尺度上呈現(xiàn)出明顯的旱濕交替周期性,且2020年后有偏濕的趨勢(shì)。不同季節(jié)的周期性和趨勢(shì)具有差異性,流域春、夏、秋、冬季分別在10、21、14、23 a尺度上呈現(xiàn)旱濕交替的可能性較大,流域春季與夏季旱澇交替的頻率較大,冬季旱澇交替頻率較小。

    (3)ENSO事件通過改變大氣環(huán)流來影響區(qū)域降水量與氣溫的變化,且其響應(yīng)具有一定的時(shí)間滯后。ENSO暖事件強(qiáng)度越大越有可能發(fā)生干旱事件,ENSO冷事件強(qiáng)度越大越有可能發(fā)生洪澇事件。ENSO暖事件的發(fā)生對(duì)流域降水量的減少存在1~3月的滯后,對(duì)流域氣溫的升高存在4~6月的滯后。

    本文采用SPEI計(jì)算出的黃河流域近60 a干旱特征值與王飛等[24]的研究結(jié)果保持一致,表明SPEI能夠很好地反映出黃河流域不同時(shí)間尺度的時(shí)空變化特征。黃河流域干旱對(duì)ENSO事件的響應(yīng)關(guān)系與周丹等[15]的研究結(jié)果相吻合。因此在厄爾尼諾事件發(fā)生后1~6月需加強(qiáng)防范,避免大面積、高強(qiáng)度的干旱現(xiàn)象發(fā)生。影響干旱的因素眾多,本文僅探究了ENSO事件對(duì)黃河流域干旱特征演變的影響,還缺乏對(duì)其他因素的考慮,因此今后應(yīng)從機(jī)理性角度分析其他因素對(duì)流域干旱特征的影響。

    參考文獻(xiàn):

    [1] 朱新玉.基于SPEI的豫東地區(qū)近50年干旱演變特征[J].自然災(zāi)害學(xué)報(bào),2015,24(4):128-137.

    [2] 屈艷萍,呂娟,蘇志誠,等.抗旱減災(zāi)研究綜述及展望[J].水利學(xué)報(bào),2018,49(1):115-125.

    [3] 袁夢(mèng),暢建霞,黎云云.基于綜合干旱指數(shù)的渭河流域干旱時(shí)空分析[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2018,51(5):401-408.

    [4] 陶然,張珂.基于PDSI的1982—2015年我國氣象干旱特征及時(shí)空變化分析[J].水資源保護(hù),2020,36(5):50-56.

    [5] 吳志勇,程丹丹,何海,等.綜合干旱指數(shù)研究進(jìn)展[J].水資源保護(hù),2021,37(1):36-45.

    [6] 劉勤,嚴(yán)昌榮,張燕卿,等.近50年黃河流域氣溫和降水量變化特征分析[J].中國農(nóng)業(yè)氣象,2012,33(4):475-480.

    [7] 劉昌明,劉小莽,田巍,等.黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展亟待解決缺水問題[J].人民黃河,2020,42(9):6-9.

    [8] 曹闖,任立良,劉懿,等.基于聯(lián)合干旱指數(shù)的黃河流域干旱時(shí)空特征[J].人民黃河,2019,41(5):51-56.

    [9] 張永永,李福軍,彭少明,等.應(yīng)對(duì)干旱的黃河干流骨干水庫群聯(lián)合調(diào)度研究[J].人民黃河,2019,41(9):31-35.

    [10] 何福力,胡彩虹,王紀(jì)軍,等.基于標(biāo)準(zhǔn)化降水、徑流指數(shù)的黃河流域近50年氣象水文干旱演變分析[J].地理與地理信息科學(xué),2015,31(3):69-75.

    [11] 楊肖麗,鄭巍斐,林長清,等.基于統(tǒng)計(jì)降尺度和SPI的黃河流域干旱預(yù)測[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,45(5):377-383.

    [12] 周帥,王義民,暢建霞,等.黃河流域干旱時(shí)空演變的空間格局研究[J].水利學(xué)報(bào),2019,50(10):1231-1241.

    [13] 李志,王健,劉文兆.涇河流域氣候變化及其與ENSO的關(guān)系[J].地理科學(xué)進(jìn)展,2010,29(7):833-839.

    [14] 張欽,唐海萍,崔鳳琪,等.基于標(biāo)準(zhǔn)化降水蒸散指數(shù)的呼倫貝爾草原干旱變化特征及趨勢(shì)分析[J].生態(tài)學(xué)報(bào),2019,39(19):7110-7123.

    [15] 周丹,張勃,安美玲,等.黃河流域不同時(shí)間尺度干旱對(duì)ENSO事件的響應(yīng)[J].中國沙漠,2015,35(3):753-762.

    [16] 溫克剛,丁一匯.中國氣象災(zāi)害大典(綜合卷)[M].北京:氣象出版社,2008:159-229.

    [17] 黃河流域及西北片水旱災(zāi)害編委會(huì).黃河流域水旱災(zāi)害[M].鄭州:黃河水利出版社,1996:275-282,372.

    [18] 劉衛(wèi)林,劉麗娜.修河流域洪水變化特征及其對(duì)氣候變化的響應(yīng)[J].水土保持研究,2018,25(5):306-312.

    [19] 張萌萌,管兆勇,張奔奔.北半球夏季海洋性大陸區(qū)域降水與南方濤動(dòng)聯(lián)系的年代際變化[J].氣象科學(xué),2019,39(6):731-738.

    [20] 劉衛(wèi)林,朱圣男,劉麗娜,等.基于SPEI的1958—2018年鄱陽湖流域干旱時(shí)空特征及其與ENSO的關(guān)系[J].中國農(nóng)村水利水電,2020(4):116-12 128.

    [21] WANG W G, SHAO Q X, PENG S Z, et al. Spatial and Temporal Characteristics of Changes in Precipitation During 1957-2007 in the Haihe River Basin, China[J].Stochastic Environmental Research and Risk Assessmeent,201 5(7):881-895.

    [22] 夏偉,周維博,李文溢,等.氣候變化和人類活動(dòng)對(duì)灃河流域徑流量影響的定量評(píng)估[J].水資源與水工程學(xué)報(bào),2018,29(6):47-52.

    [23] 李曉燕,翟盤茂,任福民.氣候標(biāo)準(zhǔn)值改變對(duì)ENSO事件劃分的影響[J].熱帶氣象學(xué)報(bào),2005,21(1):72-78.

    [24] 王飛,王宗敏,楊海波,等.基于SPEI的黃河流域干旱時(shí)空格局研究[J].中國科學(xué):地球科學(xué),2018,48(9):1169-1183.

    【責(zé)任編輯 張 帥】

    猜你喜歡
    干旱黃河流域
    生態(tài)環(huán)境部啟動(dòng)新一年度黃河流域“清廢行動(dòng)”
    “黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展”專題征稿
    人民黃河(2022年7期)2022-07-07 06:52:24
    在黃河流域生態(tài)保護(hù)和高質(zhì)量發(fā)展中展現(xiàn)陜西擔(dān)當(dāng)
    黃河流域燦爛的齊家文化——陶器
    收藏界(2019年2期)2019-10-12 08:26:10
    增強(qiáng)大局意識(shí) 提升黃河流域生態(tài)保護(hù)發(fā)展水平
    基于距平的白城地區(qū)干旱時(shí)間分布特征分析
    臨夏地區(qū)干旱特征及干濕氣候區(qū)劃
    三十六計(jì)之順手牽羊
    夏季高溫干旱時(shí)節(jié)高山蔬菜種植管理策略
    小說里的世界 虛擬里的真實(shí)
    丝袜人妻中文字幕| 欧美 日韩 精品 国产| 国产区一区二久久| 欧美日韩福利视频一区二区| 亚洲国产欧美一区二区综合| 另类亚洲欧美激情| 好男人电影高清在线观看| 天天躁日日躁夜夜躁夜夜| svipshipincom国产片| 两性夫妻黄色片| 亚洲精品粉嫩美女一区| 人人妻,人人澡人人爽秒播| 男人爽女人下面视频在线观看| 丰满迷人的少妇在线观看| 99热国产这里只有精品6| 大陆偷拍与自拍| 一本一本久久a久久精品综合妖精| 国产成人欧美| 老司机深夜福利视频在线观看 | 久久久久国产精品人妻一区二区| 我的亚洲天堂| 色综合欧美亚洲国产小说| 国产一区二区在线观看av| 日韩人妻精品一区2区三区| 老司机影院毛片| 老汉色av国产亚洲站长工具| h视频一区二区三区| 极品少妇高潮喷水抽搐| 少妇粗大呻吟视频| 性色av一级| 曰老女人黄片| 高清欧美精品videossex| 久久综合国产亚洲精品| av片东京热男人的天堂| 中亚洲国语对白在线视频| 一级,二级,三级黄色视频| 大型av网站在线播放| tube8黄色片| 久久国产精品影院| 在线观看舔阴道视频| 黄色怎么调成土黄色| 国产成人啪精品午夜网站| 成人国语在线视频| 黄色片一级片一级黄色片| 久久天躁狠狠躁夜夜2o2o| 色94色欧美一区二区| 视频区欧美日本亚洲| 久久天堂一区二区三区四区| 亚洲视频免费观看视频| 国产在线视频一区二区| 国产欧美亚洲国产| 精品一品国产午夜福利视频| 久久久欧美国产精品| 国产在线一区二区三区精| 亚洲中文字幕日韩| 亚洲国产精品成人久久小说| 亚洲欧美日韩另类电影网站| 一级a爱视频在线免费观看| 久久久久久久精品精品| 亚洲精品国产一区二区精华液| 国产精品成人在线| 久久人妻熟女aⅴ| 免费观看av网站的网址| 99国产精品一区二区蜜桃av | 亚洲免费av在线视频| 日韩有码中文字幕| 午夜精品久久久久久毛片777| 久久精品亚洲av国产电影网| 亚洲av日韩在线播放| 性高湖久久久久久久久免费观看| 久久久欧美国产精品| 国产又爽黄色视频| 成年av动漫网址| 秋霞在线观看毛片| 下体分泌物呈黄色| 丰满少妇做爰视频| 日日夜夜操网爽| 久久ye,这里只有精品| 爱豆传媒免费全集在线观看| 亚洲国产av影院在线观看| 久久性视频一级片| 亚洲国产欧美一区二区综合| 欧美xxⅹ黑人| 满18在线观看网站| 最近中文字幕2019免费版| 亚洲视频免费观看视频| 美女午夜性视频免费| 国产亚洲精品第一综合不卡| 国产在线一区二区三区精| 久久人人97超碰香蕉20202| 午夜两性在线视频| 精品少妇一区二区三区视频日本电影| 夜夜骑夜夜射夜夜干| 99九九在线精品视频| 亚洲专区字幕在线| 久久精品亚洲熟妇少妇任你| 一二三四在线观看免费中文在| 久久国产精品影院| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品一卡2卡三卡4卡5卡 | 欧美日韩精品网址| 精品少妇久久久久久888优播| www.熟女人妻精品国产| 性少妇av在线| 久久精品久久久久久噜噜老黄| 久久精品久久久久久噜噜老黄| 欧美黄色淫秽网站| 少妇的丰满在线观看| 国产97色在线日韩免费| 亚洲avbb在线观看| av在线播放精品| 中文字幕av电影在线播放| 亚洲欧美日韩高清在线视频 | 国产精品一区二区在线观看99| 啦啦啦免费观看视频1| 十八禁高潮呻吟视频| 免费在线观看影片大全网站| 亚洲精品成人av观看孕妇| 欧美日韩精品网址| 久久国产精品人妻蜜桃| 亚洲国产中文字幕在线视频| 国产成人啪精品午夜网站| 亚洲欧美成人综合另类久久久| 国产精品 国内视频| 十八禁网站免费在线| 免费在线观看完整版高清| 欧美日本中文国产一区发布| 亚洲av电影在线进入| 国产欧美日韩精品亚洲av| 国产有黄有色有爽视频| 亚洲国产欧美一区二区综合| 国产成人影院久久av| 国产黄色免费在线视频| 麻豆av在线久日| 狠狠婷婷综合久久久久久88av| 大香蕉久久网| 19禁男女啪啪无遮挡网站| 91精品三级在线观看| 男女床上黄色一级片免费看| 欧美激情 高清一区二区三区| 亚洲国产精品999| 国产伦理片在线播放av一区| 欧美黑人精品巨大| 日韩欧美一区二区三区在线观看 | 精品人妻在线不人妻| 精品人妻在线不人妻| 久久毛片免费看一区二区三区| 亚洲精品美女久久久久99蜜臀| 免费在线观看完整版高清| 午夜福利,免费看| 最近中文字幕2019免费版| 男女下面插进去视频免费观看| 欧美日韩成人在线一区二区| 超色免费av| 青草久久国产| 一进一出抽搐动态| 三级毛片av免费| 国产一区二区在线观看av| 欧美日韩亚洲国产一区二区在线观看 | 在线看a的网站| 久久精品国产亚洲av高清一级| 性高湖久久久久久久久免费观看| 国产免费现黄频在线看| 亚洲九九香蕉| 日韩三级视频一区二区三区| 欧美黑人欧美精品刺激| 永久免费av网站大全| 色94色欧美一区二区| 国产男女超爽视频在线观看| 国产成人精品在线电影| 亚洲欧美精品自产自拍| 在线看a的网站| 精品福利永久在线观看| tocl精华| 国产亚洲午夜精品一区二区久久| 黄色怎么调成土黄色| 一本综合久久免费| 国产精品一二三区在线看| 日本五十路高清| 老汉色av国产亚洲站长工具| 这个男人来自地球电影免费观看| 欧美一级毛片孕妇| av国产精品久久久久影院| 黄色怎么调成土黄色| 无限看片的www在线观看| 女性被躁到高潮视频| a 毛片基地| 久久久久久免费高清国产稀缺| 蜜桃在线观看..| 日本a在线网址| 热99re8久久精品国产| 亚洲综合色网址| av天堂久久9| 看免费av毛片| 久久精品国产综合久久久| 日韩免费高清中文字幕av| 欧美激情久久久久久爽电影 | 精品久久久久久久毛片微露脸 | 国产成人一区二区三区免费视频网站| 国产区一区二久久| 国产亚洲午夜精品一区二区久久| 制服人妻中文乱码| 日日夜夜操网爽| 国产亚洲精品一区二区www | 国产伦人伦偷精品视频| 国产精品 国内视频| 久久精品亚洲av国产电影网| 菩萨蛮人人尽说江南好唐韦庄| 欧美黄色淫秽网站| 久久久精品94久久精品| 欧美人与性动交α欧美软件| 久热爱精品视频在线9| 不卡一级毛片| 波多野结衣av一区二区av| 国产又色又爽无遮挡免| 日韩熟女老妇一区二区性免费视频| 91av网站免费观看| 亚洲精品日韩在线中文字幕| 一进一出抽搐动态| 午夜成年电影在线免费观看| 国产99久久九九免费精品| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品99久久久久| 不卡av一区二区三区| 亚洲欧美日韩另类电影网站| 国产欧美日韩一区二区三区在线| 国产又爽黄色视频| 91精品三级在线观看| 99国产精品一区二区蜜桃av | 伊人久久大香线蕉亚洲五| 亚洲男人天堂网一区| 91av网站免费观看| 欧美日韩视频精品一区| 国产成人欧美在线观看 | 99国产精品免费福利视频| 69av精品久久久久久 | 美女主播在线视频| 久久综合国产亚洲精品| 国产成+人综合+亚洲专区| 动漫黄色视频在线观看| 国产精品久久久av美女十八| 国产欧美日韩精品亚洲av| 中文欧美无线码| 欧美另类亚洲清纯唯美| 色视频在线一区二区三区| 下体分泌物呈黄色| 国产精品影院久久| 国产免费现黄频在线看| 欧美亚洲 丝袜 人妻 在线| 性少妇av在线| 国产成人av激情在线播放| 欧美国产精品va在线观看不卡| 人人妻,人人澡人人爽秒播| 美女高潮喷水抽搐中文字幕| 国产成人精品久久二区二区91| 精品国产乱码久久久久久小说| 午夜久久久在线观看| 午夜激情av网站| 欧美97在线视频| 国产男女超爽视频在线观看| 免费在线观看黄色视频的| 欧美少妇被猛烈插入视频| 老司机深夜福利视频在线观看 | 国产黄频视频在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲精品国产av成人精品| 三上悠亚av全集在线观看| 一二三四社区在线视频社区8| 亚洲一区二区三区欧美精品| 大片电影免费在线观看免费| 欧美乱码精品一区二区三区| 成人av一区二区三区在线看 | 国产精品秋霞免费鲁丝片| 91大片在线观看| 老熟妇乱子伦视频在线观看 | 视频区欧美日本亚洲| 国产色视频综合| svipshipincom国产片| 超碰成人久久| 在线观看人妻少妇| 欧美激情 高清一区二区三区| 久久99一区二区三区| 欧美一级毛片孕妇| 一本综合久久免费| 国产主播在线观看一区二区| 欧美精品人与动牲交sv欧美| 国产激情久久老熟女| 欧美人与性动交α欧美软件| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| 夜夜夜夜夜久久久久| 精品第一国产精品| 777久久人妻少妇嫩草av网站| 欧美精品高潮呻吟av久久| 一个人免费在线观看的高清视频 | 亚洲av男天堂| 久久 成人 亚洲| av线在线观看网站| 人人妻人人爽人人添夜夜欢视频| 免费看十八禁软件| 国产男人的电影天堂91| 国精品久久久久久国模美| 建设人人有责人人尽责人人享有的| tube8黄色片| 国产成人精品在线电影| 一进一出抽搐动态| 91精品伊人久久大香线蕉| 国产无遮挡羞羞视频在线观看| 美女中出高潮动态图| 久久人妻福利社区极品人妻图片| 后天国语完整版免费观看| 久久久久精品国产欧美久久久 | 美女视频免费永久观看网站| 人妻 亚洲 视频| 欧美精品人与动牲交sv欧美| 2018国产大陆天天弄谢| 99九九在线精品视频| 国产成人欧美| 久久天堂一区二区三区四区| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看 | 日本猛色少妇xxxxx猛交久久| 中文字幕精品免费在线观看视频| av在线老鸭窝| 久久av网站| 亚洲av电影在线进入| 交换朋友夫妻互换小说| 嫁个100分男人电影在线观看| 国产成人啪精品午夜网站| 欧美国产精品va在线观看不卡| 国产日韩欧美亚洲二区| 精品亚洲乱码少妇综合久久| 精品少妇黑人巨大在线播放| 欧美日韩成人在线一区二区| 国产精品免费大片| 国产精品久久久久成人av| 天天影视国产精品| 少妇裸体淫交视频免费看高清 | 亚洲成人免费电影在线观看| 欧美97在线视频| 国产男女超爽视频在线观看| 麻豆国产av国片精品| 国产一区二区 视频在线| 麻豆乱淫一区二区| 国产欧美日韩一区二区三 | 一区二区av电影网| 久久 成人 亚洲| 999久久久精品免费观看国产| 一区二区日韩欧美中文字幕| 男女免费视频国产| 777米奇影视久久| av电影中文网址| 国产精品一二三区在线看| 欧美激情 高清一区二区三区| av福利片在线| 人人妻人人添人人爽欧美一区卜| 黄色视频在线播放观看不卡| 青春草视频在线免费观看| 亚洲国产欧美一区二区综合| 50天的宝宝边吃奶边哭怎么回事| 欧美97在线视频| 国产成人精品无人区| 人妻久久中文字幕网| 下体分泌物呈黄色| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 国产高清videossex| 国产精品熟女久久久久浪| 亚洲国产看品久久| 91av网站免费观看| 国产极品粉嫩免费观看在线| 免费一级毛片在线播放高清视频 | 亚洲人成电影观看| 亚洲欧洲日产国产| 一二三四在线观看免费中文在| 久热这里只有精品99| 亚洲精品日韩在线中文字幕| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播| 日韩欧美免费精品| av网站在线播放免费| 国产精品一区二区精品视频观看| 热99久久久久精品小说推荐| 一区二区av电影网| 国产免费一区二区三区四区乱码| 免费黄频网站在线观看国产| 免费观看人在逋| 亚洲精品国产色婷婷电影| 一区二区三区激情视频| 久久av网站| 国产成人av教育| 99久久99久久久精品蜜桃| 日本91视频免费播放| 亚洲人成电影观看| 国产成人免费无遮挡视频| 99国产极品粉嫩在线观看| 伊人亚洲综合成人网| 日本av手机在线免费观看| 亚洲专区中文字幕在线| 精品免费久久久久久久清纯 | 少妇被粗大的猛进出69影院| 国产熟女午夜一区二区三区| 超碰成人久久| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 久久精品国产综合久久久| 欧美精品高潮呻吟av久久| 人成视频在线观看免费观看| 爱豆传媒免费全集在线观看| 老司机在亚洲福利影院| 欧美黄色片欧美黄色片| 日韩一卡2卡3卡4卡2021年| 老司机亚洲免费影院| 美女中出高潮动态图| 日韩大码丰满熟妇| 久久 成人 亚洲| 啦啦啦免费观看视频1| 又大又爽又粗| 少妇猛男粗大的猛烈进出视频| 无遮挡黄片免费观看| 亚洲精华国产精华精| 精品少妇内射三级| 青春草视频在线免费观看| 如日韩欧美国产精品一区二区三区| 亚洲色图综合在线观看| 国产精品影院久久| cao死你这个sao货| 两个人看的免费小视频| av天堂在线播放| 欧美精品人与动牲交sv欧美| 欧美另类一区| 亚洲熟女毛片儿| 国产精品二区激情视频| 日本av免费视频播放| 国产1区2区3区精品| 黄网站色视频无遮挡免费观看| 日韩欧美一区二区三区在线观看 | 性高湖久久久久久久久免费观看| 亚洲av欧美aⅴ国产| 一级毛片女人18水好多| 91成人精品电影| 亚洲精品日韩在线中文字幕| 久久国产亚洲av麻豆专区| 啦啦啦免费观看视频1| 欧美亚洲 丝袜 人妻 在线| 亚洲中文日韩欧美视频| 国产精品 欧美亚洲| 妹子高潮喷水视频| 丝袜美腿诱惑在线| 中国国产av一级| av在线app专区| 男女边摸边吃奶| 91精品国产国语对白视频| 男女无遮挡免费网站观看| 欧美成狂野欧美在线观看| av电影中文网址| 国产欧美日韩一区二区三 | 欧美少妇被猛烈插入视频| 中文字幕制服av| 亚洲国产精品一区二区三区在线| 91大片在线观看| 新久久久久国产一级毛片| 亚洲av国产av综合av卡| 色播在线永久视频| 久久精品国产亚洲av香蕉五月 | 亚洲成人免费电影在线观看| 成人亚洲精品一区在线观看| 亚洲综合色网址| 国产成人欧美| 国产成人精品久久二区二区91| 久久青草综合色| av视频免费观看在线观看| 免费不卡黄色视频| 少妇粗大呻吟视频| 男女边摸边吃奶| 欧美乱码精品一区二区三区| 国产伦理片在线播放av一区| 亚洲精品乱久久久久久| 另类精品久久| 高清黄色对白视频在线免费看| 中文字幕高清在线视频| 1024香蕉在线观看| 狂野欧美激情性xxxx| 久久久久久久国产电影| 电影成人av| 一区二区三区乱码不卡18| 一本一本久久a久久精品综合妖精| 欧美另类亚洲清纯唯美| 午夜精品国产一区二区电影| 汤姆久久久久久久影院中文字幕| 午夜福利视频精品| 亚洲男人天堂网一区| 色婷婷av一区二区三区视频| 另类精品久久| 狠狠狠狠99中文字幕| 在线观看免费高清a一片| 18禁国产床啪视频网站| 久久久久国产一级毛片高清牌| 老熟妇仑乱视频hdxx| 这个男人来自地球电影免费观看| 国产精品久久久久成人av| 欧美黄色片欧美黄色片| av天堂久久9| 久久狼人影院| 爱豆传媒免费全集在线观看| 91字幕亚洲| 欧美日韩视频精品一区| 手机成人av网站| 成在线人永久免费视频| 日日夜夜操网爽| 国产一卡二卡三卡精品| 亚洲成国产人片在线观看| 美女视频免费永久观看网站| 天天操日日干夜夜撸| 亚洲欧美一区二区三区久久| a级片在线免费高清观看视频| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放| 免费看十八禁软件| 女性被躁到高潮视频| 免费高清在线观看日韩| 亚洲色图综合在线观看| 80岁老熟妇乱子伦牲交| 欧美一级毛片孕妇| 日本91视频免费播放| 夜夜夜夜夜久久久久| 999久久久精品免费观看国产| 久久99一区二区三区| 亚洲av成人一区二区三| 午夜福利在线免费观看网站| 一级毛片女人18水好多| 国产欧美日韩一区二区三 | 国产免费一区二区三区四区乱码| 老汉色av国产亚洲站长工具| 搡老乐熟女国产| 无遮挡黄片免费观看| 欧美日韩一级在线毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜老司机福利片| 超碰97精品在线观看| 波多野结衣一区麻豆| 另类精品久久| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩高清在线视频 | 午夜福利乱码中文字幕| 午夜久久久在线观看| 热99久久久久精品小说推荐| 亚洲综合色网址| 国产亚洲精品第一综合不卡| 99热国产这里只有精品6| 亚洲三区欧美一区| 日韩人妻精品一区2区三区| 美女国产高潮福利片在线看| 丝袜美腿诱惑在线| 免费一级毛片在线播放高清视频 | 国产精品一区二区在线不卡| 美女中出高潮动态图| 亚洲精品乱久久久久久| 嫩草影视91久久| 青草久久国产| 9热在线视频观看99| 国产亚洲精品久久久久5区| 午夜日韩欧美国产| 啦啦啦 在线观看视频| 国产亚洲午夜精品一区二区久久| 少妇精品久久久久久久| 少妇人妻久久综合中文| 国产在线一区二区三区精| 一区二区av电影网| www.自偷自拍.com| 亚洲人成电影免费在线| 999精品在线视频| 久久性视频一级片| 麻豆乱淫一区二区| av超薄肉色丝袜交足视频| a级毛片黄视频| 国产成人精品久久二区二区91| 夜夜骑夜夜射夜夜干| 日韩欧美免费精品| 欧美人与性动交α欧美软件| 欧美激情高清一区二区三区| 色播在线永久视频| 午夜福利在线免费观看网站| 午夜日韩欧美国产| 老司机影院毛片| 精品卡一卡二卡四卡免费| 日韩三级视频一区二区三区| 日本av手机在线免费观看| 日本精品一区二区三区蜜桃| 人人澡人人妻人| 久久av网站| 欧美久久黑人一区二区| 热99久久久久精品小说推荐| 黑人猛操日本美女一级片| 他把我摸到了高潮在线观看 | 欧美中文综合在线视频| 他把我摸到了高潮在线观看 | 久久久久视频综合| 在线观看人妻少妇| 国产成人欧美在线观看 | 91麻豆av在线| 国产淫语在线视频| 一级毛片电影观看| 黄片小视频在线播放| 色视频在线一区二区三区| 成年人免费黄色播放视频| 日日摸夜夜添夜夜添小说| 国产一区二区激情短视频 | 美女高潮到喷水免费观看| 精品卡一卡二卡四卡免费| 高清黄色对白视频在线免费看|