• <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í)
    性色av一级| 又黄又粗又硬又大视频| 9191精品国产免费久久| 日日爽夜夜爽网站| 精品熟女少妇八av免费久了| 国产精品九九99| 91麻豆精品激情在线观看国产 | 纯流量卡能插随身wifi吗| 在线观看www视频免费| 人人妻,人人澡人人爽秒播| 另类亚洲欧美激情| 黄色视频,在线免费观看| 国产成人精品在线电影| 91国产中文字幕| 欧美国产精品一级二级三级| 久久久久国产精品人妻一区二区| 国产精品国产av在线观看| 狂野欧美激情性xxxx| 在线看a的网站| 大陆偷拍与自拍| 天堂中文最新版在线下载| 最近最新免费中文字幕在线| 亚洲精品久久久久久婷婷小说| 各种免费的搞黄视频| 9191精品国产免费久久| 搡老岳熟女国产| 精品国产乱码久久久久久男人| 黄色片一级片一级黄色片| 亚洲第一欧美日韩一区二区三区 | 亚洲精品中文字幕一二三四区 | 女人高潮潮喷娇喘18禁视频| 搡老岳熟女国产| 一区二区日韩欧美中文字幕| 久久亚洲国产成人精品v| 99国产精品免费福利视频| 这个男人来自地球电影免费观看| 黑人欧美特级aaaaaa片| 欧美成人午夜精品| 新久久久久国产一级毛片| 国产精品一区二区精品视频观看| 大香蕉久久成人网| 亚洲七黄色美女视频| av网站免费在线观看视频| 国产精品成人在线| 在线av久久热| 男女国产视频网站| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影 | 操出白浆在线播放| e午夜精品久久久久久久| 欧美精品一区二区免费开放| 99香蕉大伊视频| 久久久国产成人免费| 国产男女内射视频| av在线播放精品| 夜夜骑夜夜射夜夜干| www日本在线高清视频| 女人精品久久久久毛片| 操美女的视频在线观看| 操美女的视频在线观看| 日韩人妻精品一区2区三区| 久久中文字幕一级| 亚洲人成77777在线视频| 香蕉国产在线看| 久久久久国产一级毛片高清牌| 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美激情在线| 久久中文看片网| 亚洲人成电影观看| 亚洲精品国产色婷婷电影| 亚洲欧美精品综合一区二区三区| 91麻豆精品激情在线观看国产 | 少妇粗大呻吟视频| 亚洲久久久国产精品| 午夜福利免费观看在线| 母亲3免费完整高清在线观看| 最近最新中文字幕大全免费视频| 肉色欧美久久久久久久蜜桃| 成年动漫av网址| 日韩视频在线欧美| 80岁老熟妇乱子伦牲交| 在线精品无人区一区二区三| 国产福利在线免费观看视频| 国产欧美日韩综合在线一区二区| 美女扒开内裤让男人捅视频| 黑人欧美特级aaaaaa片| 老司机午夜福利在线观看视频 | 久久人妻福利社区极品人妻图片| 亚洲激情五月婷婷啪啪| 国产国语露脸激情在线看| 国产精品免费大片| 老司机午夜福利在线观看视频 | 两性午夜刺激爽爽歪歪视频在线观看 | 精品亚洲成a人片在线观看| 国产av精品麻豆| e午夜精品久久久久久久| 亚洲欧美清纯卡通| 在线永久观看黄色视频| 桃红色精品国产亚洲av| 男男h啪啪无遮挡| 黑人巨大精品欧美一区二区蜜桃| 久久精品国产亚洲av高清一级| 精品欧美一区二区三区在线| 侵犯人妻中文字幕一二三四区| 国产精品久久久人人做人人爽| 最新的欧美精品一区二区| 中文字幕最新亚洲高清| 日韩欧美国产一区二区入口| 国产精品av久久久久免费| 中文字幕色久视频| a 毛片基地| 一区福利在线观看| 狠狠狠狠99中文字幕| 国产精品偷伦视频观看了| 亚洲国产av影院在线观看| 欧美黄色淫秽网站| 日日爽夜夜爽网站| 一区二区三区激情视频| 日本五十路高清| av视频免费观看在线观看| 国产亚洲欧美在线一区二区| 日本a在线网址| 久久香蕉激情| 亚洲伊人久久精品综合| 丰满人妻熟妇乱又伦精品不卡| 我要看黄色一级片免费的| 三上悠亚av全集在线观看| 欧美黑人欧美精品刺激| 99九九在线精品视频| 十八禁网站网址无遮挡| 高清视频免费观看一区二区| 免费高清在线观看日韩| 欧美亚洲 丝袜 人妻 在线| 欧美精品av麻豆av| 日韩 亚洲 欧美在线| 一个人免费看片子| 叶爱在线成人免费视频播放| 成人黄色视频免费在线看| 制服人妻中文乱码| 一边摸一边抽搐一进一出视频| 欧美 日韩 精品 国产| 国产主播在线观看一区二区| 蜜桃在线观看..| 永久免费av网站大全| 国产97色在线日韩免费| 中文字幕高清在线视频| 亚洲三区欧美一区| 久久久欧美国产精品| 国产成人精品无人区| 黄色片一级片一级黄色片| av又黄又爽大尺度在线免费看| 欧美少妇被猛烈插入视频| 日日爽夜夜爽网站| 亚洲自偷自拍图片 自拍| 午夜91福利影院| 国产欧美日韩一区二区三区在线| 热99久久久久精品小说推荐| 90打野战视频偷拍视频| 久久久久国内视频| 成人三级做爰电影| 国产男人的电影天堂91| 亚洲av电影在线观看一区二区三区| 在线十欧美十亚洲十日本专区| 黄色 视频免费看| 18禁国产床啪视频网站| 一级片免费观看大全| 一级片'在线观看视频| 国产在视频线精品| 国产在线一区二区三区精| 成年人免费黄色播放视频| 成年人黄色毛片网站| 多毛熟女@视频| 超碰成人久久| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 叶爱在线成人免费视频播放| 亚洲欧美成人综合另类久久久| 各种免费的搞黄视频| 中文字幕高清在线视频| 国产精品久久久av美女十八| 叶爱在线成人免费视频播放| 十分钟在线观看高清视频www| 99国产精品99久久久久| 国产精品一二三区在线看| 日韩精品免费视频一区二区三区| 美女脱内裤让男人舔精品视频| 亚洲美女黄色视频免费看| 国产激情久久老熟女| 国产成+人综合+亚洲专区| 国产无遮挡羞羞视频在线观看| 亚洲av美国av| 亚洲色图 男人天堂 中文字幕| 人人澡人人妻人| 丝袜脚勾引网站| 另类精品久久| 亚洲伊人色综图| www.av在线官网国产| 汤姆久久久久久久影院中文字幕| svipshipincom国产片| 亚洲中文日韩欧美视频| 一边摸一边抽搐一进一出视频| 啦啦啦啦在线视频资源| av欧美777| 狠狠婷婷综合久久久久久88av| 国产成人影院久久av| www.熟女人妻精品国产| 岛国在线观看网站| 欧美av亚洲av综合av国产av| 一个人免费在线观看的高清视频 | 国产黄色免费在线视频| 日韩中文字幕视频在线看片| 少妇猛男粗大的猛烈进出视频| 热re99久久精品国产66热6| 日韩欧美国产一区二区入口| 狠狠婷婷综合久久久久久88av| 国产精品免费大片| 日韩 亚洲 欧美在线| 久久毛片免费看一区二区三区| av超薄肉色丝袜交足视频| 国产精品一区二区免费欧美 | 美女大奶头黄色视频| 久久久精品94久久精品| 久久99热这里只频精品6学生| 久久人人97超碰香蕉20202| 精品人妻熟女毛片av久久网站| 制服诱惑二区| 一区二区三区四区激情视频| 免费在线观看视频国产中文字幕亚洲 | 精品国内亚洲2022精品成人 | 91av网站免费观看| 肉色欧美久久久久久久蜜桃| 亚洲熟女精品中文字幕| av又黄又爽大尺度在线免费看| 淫妇啪啪啪对白视频 | 99久久人妻综合| 女人久久www免费人成看片| 国产熟女午夜一区二区三区| 亚洲国产中文字幕在线视频| 丝袜脚勾引网站| 国产xxxxx性猛交| 成年美女黄网站色视频大全免费| 国产男人的电影天堂91| 免费日韩欧美在线观看| 五月开心婷婷网| 最近中文字幕2019免费版| 亚洲国产欧美一区二区综合| 欧美黄色淫秽网站| 中文字幕精品免费在线观看视频| 免费黄频网站在线观看国产| 国产一区二区三区在线臀色熟女 | xxxhd国产人妻xxx| 亚洲伊人色综图| 91av网站免费观看| 欧美激情极品国产一区二区三区| 成人亚洲精品一区在线观看| 中文字幕人妻丝袜一区二区| 日本撒尿小便嘘嘘汇集6| 色老头精品视频在线观看| 99re6热这里在线精品视频| 国产极品粉嫩免费观看在线| 久久精品国产综合久久久| 男人舔女人的私密视频| 国产精品av久久久久免费| 一本色道久久久久久精品综合| 我要看黄色一级片免费的| 91麻豆精品激情在线观看国产 | 亚洲中文日韩欧美视频| 成人国产一区最新在线观看| 久久中文字幕一级| 99精品久久久久人妻精品| 亚洲欧美日韩高清在线视频 | 又紧又爽又黄一区二区| 国产一卡二卡三卡精品| 免费人妻精品一区二区三区视频| 久久久国产一区二区| 超碰成人久久| 日韩欧美一区视频在线观看| 永久免费av网站大全| 9色porny在线观看| 热re99久久精品国产66热6| 久久女婷五月综合色啪小说| 可以免费在线观看a视频的电影网站| 色婷婷久久久亚洲欧美| 国产伦理片在线播放av一区| 女警被强在线播放| 99国产精品一区二区三区| 男女无遮挡免费网站观看| 91国产中文字幕| 午夜福利视频精品| 精品国产一区二区三区四区第35| 中文字幕人妻丝袜制服| 亚洲第一欧美日韩一区二区三区 | xxxhd国产人妻xxx| 精品亚洲乱码少妇综合久久| 啪啪无遮挡十八禁网站| 50天的宝宝边吃奶边哭怎么回事| 人人妻,人人澡人人爽秒播| 最新的欧美精品一区二区| 国产成人a∨麻豆精品| 男人添女人高潮全过程视频| 美女脱内裤让男人舔精品视频| 丰满人妻熟妇乱又伦精品不卡| 飞空精品影院首页| 叶爱在线成人免费视频播放| 视频在线观看一区二区三区| 日本av手机在线免费观看| 啦啦啦啦在线视频资源| 19禁男女啪啪无遮挡网站| 99精品久久久久人妻精品| 午夜成年电影在线免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品日韩在线中文字幕| 丰满迷人的少妇在线观看| 美女视频免费永久观看网站| 欧美激情极品国产一区二区三区| 天堂中文最新版在线下载| 国产区一区二久久| 天堂8中文在线网| 国产精品国产三级国产专区5o| 人妻一区二区av| 免费高清在线观看视频在线观看| 正在播放国产对白刺激| 韩国高清视频一区二区三区| 亚洲av美国av| bbb黄色大片| 老司机亚洲免费影院| 国产精品香港三级国产av潘金莲| 男女下面插进去视频免费观看| 久久午夜综合久久蜜桃| 亚洲伊人色综图| 午夜精品久久久久久毛片777| 亚洲熟女毛片儿| 视频在线观看一区二区三区| 国产日韩欧美视频二区| 久久性视频一级片| 丝袜在线中文字幕| 久9热在线精品视频| 国产亚洲一区二区精品| 国产色视频综合| 少妇 在线观看| 69av精品久久久久久 | 涩涩av久久男人的天堂| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 可以免费在线观看a视频的电影网站| 新久久久久国产一级毛片| 精品久久久久久电影网| 国产日韩欧美亚洲二区| 性少妇av在线| 欧美精品高潮呻吟av久久| av片东京热男人的天堂| av在线老鸭窝| 婷婷丁香在线五月| 久久精品aⅴ一区二区三区四区| 亚洲综合色网址| 丝袜喷水一区| 亚洲国产欧美在线一区| 国产淫语在线视频| 亚洲精品国产一区二区精华液| 男人添女人高潮全过程视频| 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 免费av中文字幕在线| 久热这里只有精品99| 欧美另类一区| 美女高潮喷水抽搐中文字幕| 高潮久久久久久久久久久不卡| 中文字幕人妻丝袜一区二区| 日日摸夜夜添夜夜添小说| 亚洲国产日韩一区二区| 大片免费播放器 马上看| 国产成人精品在线电影| 日韩免费高清中文字幕av| 99精国产麻豆久久婷婷| 欧美精品啪啪一区二区三区 | 在线观看免费高清a一片| 亚洲精品国产精品久久久不卡| 啦啦啦视频在线资源免费观看| 国产精品九九99| 国产成人免费无遮挡视频| 欧美老熟妇乱子伦牲交| 亚洲 欧美一区二区三区| 日韩三级视频一区二区三区| 国产成人a∨麻豆精品| 亚洲av电影在线进入| 久久女婷五月综合色啪小说| 国产av国产精品国产| 国产一区二区三区综合在线观看| 少妇猛男粗大的猛烈进出视频| 成人亚洲精品一区在线观看| 无限看片的www在线观看| 黑丝袜美女国产一区| 9191精品国产免费久久| 亚洲欧美精品自产自拍| 亚洲专区国产一区二区| 香蕉国产在线看| 国产高清视频在线播放一区 | 欧美精品高潮呻吟av久久| av又黄又爽大尺度在线免费看| 亚洲国产看品久久| 国产精品一区二区在线不卡| 91精品三级在线观看| 亚洲精品av麻豆狂野| 国产亚洲欧美在线一区二区| 丁香六月欧美| 国产有黄有色有爽视频| a 毛片基地| 亚洲精品日韩在线中文字幕| 丝袜在线中文字幕| 日韩一卡2卡3卡4卡2021年| 丝袜喷水一区| √禁漫天堂资源中文www| 91麻豆精品激情在线观看国产 | 肉色欧美久久久久久久蜜桃| 亚洲欧美一区二区三区久久| 男女之事视频高清在线观看| 亚洲自偷自拍图片 自拍| 大型av网站在线播放| 午夜激情av网站| 亚洲人成电影观看| 12—13女人毛片做爰片一| 99精品久久久久人妻精品| 日韩欧美国产一区二区入口| 91麻豆av在线| 亚洲成人国产一区在线观看| 一区二区三区四区激情视频| www.av在线官网国产| 大片免费播放器 马上看| 免费观看人在逋| 亚洲午夜精品一区,二区,三区| 国产精品 国内视频| 秋霞在线观看毛片| 免费在线观看黄色视频的| 精品亚洲成a人片在线观看| 777米奇影视久久| 午夜福利乱码中文字幕| 久久久久精品国产欧美久久久 | 欧美少妇被猛烈插入视频| 黄片小视频在线播放| 窝窝影院91人妻| 亚洲性夜色夜夜综合| 亚洲自偷自拍图片 自拍| 电影成人av| 999久久久精品免费观看国产| 国产国语露脸激情在线看| 丝袜喷水一区| 欧美精品av麻豆av| 美女高潮到喷水免费观看| 日韩免费高清中文字幕av| av网站免费在线观看视频| 国产亚洲一区二区精品| 免费在线观看黄色视频的| 涩涩av久久男人的天堂| 欧美少妇被猛烈插入视频| 十八禁网站免费在线| 王馨瑶露胸无遮挡在线观看| 亚洲黑人精品在线| 久久久精品区二区三区| 色婷婷av一区二区三区视频| 男女无遮挡免费网站观看| 亚洲国产精品一区三区| 国产精品一区二区在线观看99| 国产精品熟女久久久久浪| 18在线观看网站| 久久久久久久国产电影| 成年女人毛片免费观看观看9 | 国产精品国产三级国产专区5o| av福利片在线| 天天躁狠狠躁夜夜躁狠狠躁| 制服人妻中文乱码| 一个人免费在线观看的高清视频 | 俄罗斯特黄特色一大片| 国产精品99久久99久久久不卡| 国产一区二区三区av在线| 午夜免费成人在线视频| 美女扒开内裤让男人捅视频| 一本一本久久a久久精品综合妖精| 国产成人a∨麻豆精品| 亚洲av日韩在线播放| 一区二区三区激情视频| 日韩制服骚丝袜av| 两个人免费观看高清视频| 人人澡人人妻人| 最近中文字幕2019免费版| av免费在线观看网站| 99香蕉大伊视频| 国产精品麻豆人妻色哟哟久久| 久久ye,这里只有精品| 19禁男女啪啪无遮挡网站| 超碰成人久久| 日本av免费视频播放| av在线app专区| 男女高潮啪啪啪动态图| 亚洲精品国产精品久久久不卡| 超碰成人久久| 日本av免费视频播放| 在线观看人妻少妇| 伦理电影免费视频| 国产免费现黄频在线看| 97人妻天天添夜夜摸| 中文字幕人妻熟女乱码| 老熟女久久久| 两人在一起打扑克的视频| 国产在线一区二区三区精| 男人操女人黄网站| 一边摸一边抽搐一进一出视频| 亚洲av日韩精品久久久久久密| 成人影院久久| 久久精品aⅴ一区二区三区四区| 亚洲国产av新网站| av天堂在线播放| 菩萨蛮人人尽说江南好唐韦庄| 国产亚洲av片在线观看秒播厂| 欧美精品av麻豆av| 免费高清在线观看日韩| 亚洲精品成人av观看孕妇| 精品国产一区二区久久| 国产免费av片在线观看野外av| 午夜福利视频精品| a级毛片在线看网站| 天天躁夜夜躁狠狠躁躁| 夫妻午夜视频| 大香蕉久久网| 九色亚洲精品在线播放| 日本wwww免费看| h视频一区二区三区| 黄色毛片三级朝国网站| 欧美在线一区亚洲| 欧美 日韩 精品 国产| 国产男人的电影天堂91| 成年人免费黄色播放视频| 女性被躁到高潮视频| 欧美人与性动交α欧美精品济南到| 国产日韩欧美在线精品| 高潮久久久久久久久久久不卡| 天堂中文最新版在线下载| 国产无遮挡羞羞视频在线观看| 亚洲精品乱久久久久久| 亚洲精品美女久久久久99蜜臀| 一本综合久久免费| 欧美成狂野欧美在线观看| av在线老鸭窝| 一级毛片女人18水好多| 男女免费视频国产| 免费在线观看影片大全网站| 亚洲avbb在线观看| 19禁男女啪啪无遮挡网站| 亚洲国产精品成人久久小说| 日本一区二区免费在线视频| 欧美变态另类bdsm刘玥| 国产免费一区二区三区四区乱码| 丝袜人妻中文字幕| 日本91视频免费播放| 我的亚洲天堂| cao死你这个sao货| 999久久久国产精品视频| 一级毛片电影观看| 最近中文字幕2019免费版| 国产亚洲一区二区精品| 欧美一级毛片孕妇| 亚洲国产日韩一区二区| 777久久人妻少妇嫩草av网站| 日本wwww免费看| 18禁国产床啪视频网站| av不卡在线播放| 视频在线观看一区二区三区| 91成人精品电影| 高潮久久久久久久久久久不卡| 亚洲一码二码三码区别大吗| 伊人久久大香线蕉亚洲五| 欧美少妇被猛烈插入视频| 国产av精品麻豆| 亚洲精品国产区一区二| 国产激情久久老熟女| 午夜福利在线免费观看网站| 亚洲三区欧美一区| 大码成人一级视频| 啦啦啦 在线观看视频| 亚洲午夜精品一区,二区,三区| 国产成人一区二区三区免费视频网站| 亚洲精品一二三| 国产在视频线精品| 别揉我奶头~嗯~啊~动态视频 | 女警被强在线播放| 精品人妻一区二区三区麻豆| 国产亚洲欧美精品永久| 国产一卡二卡三卡精品| 国产1区2区3区精品| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 欧美精品一区二区免费开放| 黑人巨大精品欧美一区二区蜜桃| av又黄又爽大尺度在线免费看| kizo精华| 男男h啪啪无遮挡| 男女下面插进去视频免费观看| 一区二区三区激情视频| 精品一区二区三区四区五区乱码| 精品国产一区二区三区四区第35| 我的亚洲天堂| 精品国产乱码久久久久久男人| 亚洲精品第二区| 久久九九热精品免费| 97在线人人人人妻| 男女高潮啪啪啪动态图| 国产一区二区三区综合在线观看| 老司机在亚洲福利影院| 日韩一卡2卡3卡4卡2021年| 99re6热这里在线精品视频| 777久久人妻少妇嫩草av网站| 最新在线观看一区二区三区| 啦啦啦 在线观看视频| 麻豆乱淫一区二区| 色综合欧美亚洲国产小说| 免费高清在线观看日韩|