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

    開都河流域極端水文事件變化及其對氣候的響應(yīng)

    2014-09-21 08:51:21李佳秀徐長春高沈瞳
    水土保持研究 2014年5期
    關(guān)鍵詞:開都河枯水年際

    李佳秀, 徐長春, 王 曉, 高沈瞳, 趙 杰

    (1.新疆大學(xué) 綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830046;2.新疆大學(xué) 資源與環(huán)境科學(xué)學(xué)院, 烏魯木齊 830046; 3.民族學(xué)院人事處, 四川 康定 616001)

    開都河流域極端水文事件變化及其對氣候的響應(yīng)

    李佳秀1,2, 徐長春1,2, 王 曉3, 高沈瞳1,2, 趙 杰1,2

    (1.新疆大學(xué) 綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室, 烏魯木齊 830046;2.新疆大學(xué) 資源與環(huán)境科學(xué)學(xué)院, 烏魯木齊 830046; 3.民族學(xué)院人事處, 四川 康定 616001)

    基于開都河流域1972—2012年大山口水文站及巴音布魯克氣象站資料,采用累積頻率法定義極端水文事件的閾值,分析了極端水文事件頻次的年內(nèi)年際變化、場次分布特征及其對區(qū)域氣候變化的響應(yīng)。結(jié)果表明:極端水文事件具有極不穩(wěn)定的特性。年內(nèi)分配不均且集中度較高,極端洪水事件主要集中在7月,極端枯水事件主要集中在2月;年際變化極不穩(wěn)定,且偏離正態(tài)分布的程度較高。極端洪水事件總體呈增長趨勢,并在1986年發(fā)生突變;極端枯水事件則呈減少趨勢,突變點(diǎn)為1993年。極端洪水頻次與降水量相關(guān)性較高,極端枯水頻次與氣溫的負(fù)相關(guān)性較高。

    極端洪水事件; 極端枯水事件; 累積頻率; 集中度; 輪次分析

    IPCC(International Panel on Climate Change)第四次評估綜合報(bào)告指出,近百年來,受全球變暖的影響,氣候異?,F(xiàn)象時(shí)有發(fā)生,下墊面的能量和水分循環(huán)特征發(fā)生了很大變化,將改變?nèi)蛩难h(huán)的現(xiàn)狀[1]。導(dǎo)致高溫干旱和暴雨洪澇等極端事件的頻次和強(qiáng)度加強(qiáng),增加水災(zāi)害風(fēng)險(xiǎn)[2]。大量研究表明,世界各地因極端水文事件所造成的直接經(jīng)濟(jì)損失呈指數(shù)上升趨勢[3-5]。極端水文事件比一般水文事件成因上更為復(fù)雜,預(yù)測更困難,因此,愈來愈受到人們的廣泛關(guān)注[6]。在國際上,將極端水文事件定義為:在一定時(shí)期內(nèi)發(fā)生頻率低、突發(fā)性強(qiáng),且對人類社會有重要影響的災(zāi)害性水文事件[7]。我國學(xué)者夏軍等[8]利用Copula函數(shù)構(gòu)建氣候要素同極端水文事件之間的多元統(tǒng)計(jì)模型,分析了不同降水條件下不同等級干旱和洪水的發(fā)生概率變化。孫桂麗等[9-10]研究了新疆極端水文事件的時(shí)空分布特征、年際變化及對氣候變化的響應(yīng),研究結(jié)果表明新疆極端水文事件年內(nèi)年際分配不均,且集中度高,極端水文事件增長趨勢與氣溫、降水變化一致。

    對于流域而言,極端水文事件主要包括極端洪水和極端枯水事件。其實(shí)質(zhì)就是洪水和枯水事件趨于極值化的研究。目前,在水文學(xué)上洪水事件多是從年極大值的角度來進(jìn)行頻率曲線的分析,用頻率的大小來推求洪水的重現(xiàn)期,謝平等基于希爾伯特—黃變換的非一致性洪水頻率計(jì)算方法推求西江中下游大湟江口站過去、現(xiàn)狀、未來3個時(shí)期的洪水頻率分布[11];枯水事件是從年或月極小值的角度來進(jìn)行頻率分析的,方彬等采用混合分布的方法對枯水事件年最小流量的頻率分布進(jìn)行了探究[12]。專門針對極端洪水和極端枯水事件的研究還較少。開都河流域位于天山南坡,以農(nóng)業(yè)經(jīng)濟(jì)發(fā)展為主,其水資源主要來自于開都河上游山區(qū),占到流域水資源總量的84.6%[13]。因此,開都河來水的穩(wěn)定與否決定著整個流域乃至尾閭博斯騰湖和下游孔雀河流域的社會經(jīng)濟(jì)的正常發(fā)展。對開都河來水過程特別是極端水文過程進(jìn)行分析,掌握和了解其發(fā)生發(fā)展規(guī)律,對區(qū)域經(jīng)濟(jì)社會發(fā)展和生態(tài)安全維護(hù)具有重要指導(dǎo)意義。本文基于開都河流域1972—2012年大山口站逐日流量資料以及巴音布魯克站氣象資料,研究分析極端水文事件的時(shí)間分布特征、長期變化趨勢及其對氣候變化的響應(yīng),以期為流域可持續(xù)發(fā)展戰(zhàn)略規(guī)劃和防災(zāi)減災(zāi)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    開都河流域地處歐亞大陸腹地,屬干旱、半干旱性氣候。開都河位于新疆天山南坡焉耆盆地北緣,位于82°58′—86°55′E、41°47′—43°21′N。發(fā)源于天山中部海拔5 000 m的薩爾明山的哈爾尕山和扎克斯臺溝,最高山峰為克勒代烏拉,海拔4 697 m,河源高山區(qū)終年積雪,有現(xiàn)代冰川840條[14],流域高程3 600 m以上終年積雪。出山口以上流域集水面積19 012 km2,冰川覆蓋面積474.98 km2,占出山口以上流域面積的2.5%。大山口水文站以上山區(qū)為開都河的上游區(qū),亦即產(chǎn)流區(qū),該地區(qū)除放牧外無其他人類活動干擾。開都河屬于冰雪融水和雨水混合補(bǔ)給為主的河流,春季隨氣溫升高,以季節(jié)性積雪融化補(bǔ)給為主;夏季以冰川融水和降雨補(bǔ)給為主[15]。多年徑流量達(dá)34.4億m3,其中冰川融水約占年出山徑流量的15.2%,融雪與降水補(bǔ)給合計(jì)約占80%以上[16]。流域總地勢是北高南低,高山、峽谷和盆地交錯分布,地形復(fù)雜多樣。開都河全長560 km,主流自東向西經(jīng)小尤爾都斯盆地至巴音布魯克水文站,而后折轉(zhuǎn)東南,經(jīng)大尤爾都斯盆地至呼斯臺西里,再經(jīng)峽谷段至流域控制站大山口水文站后流出山口,而后進(jìn)入焉耆盆地后流入博斯騰湖,是惟一能常年補(bǔ)給博斯騰湖的河流。開都河作為塔里木河的源流,是實(shí)施北水南調(diào)工程的關(guān)鍵河流,一直承擔(dān)著向塔河下游生態(tài)應(yīng)急輸水的重要任務(wù),因此開都河不僅是南疆綠洲的生態(tài)源,還是重要的生命源。開都河區(qū)域內(nèi)的生態(tài)植被保持完好,天然荒漠植被郁閉度和覆蓋度不斷得到提高,野生動植物種群逐年增加,并形成了人與自然的和諧相處良好局面。開都河與從博斯騰湖流出的孔雀河構(gòu)成面積達(dá)6萬km2的開孔河流域,流域內(nèi)人口和國民生產(chǎn)總值均占到全州的80%以上。

    2 數(shù)據(jù)資料與方法

    2.1 數(shù)據(jù)資料

    本文分析用數(shù)據(jù)主要包括水文數(shù)據(jù)和氣象數(shù)據(jù)。其中,水文數(shù)據(jù)為開都河流域大山口水文站(85°44′E,42°13′N)1972—2012年的日流數(shù)據(jù);氣象數(shù)據(jù)為開都河流域巴音布魯克氣象站(84°09′E,43°02′N)1972—2012年的月均氣溫和降水?dāng)?shù)據(jù)。數(shù)據(jù)資料完整、一致,達(dá)到分析用標(biāo)準(zhǔn)。

    2.2 分析方法

    2.2.1 極端水文事件閾值的確定 國際上對極端水文事件的定義通常根據(jù)流量、水位或造成的損失大小劃分為不同的標(biāo)準(zhǔn)。對極端事件閾值的確定通常采用百分位閾值法,即從概率統(tǒng)計(jì)的角度來定義極端事件。累積頻率法與百分位法相類似,本文運(yùn)用累積頻率法確定極端水文事件的閾值。

    累積頻率(CFD,cumulative frequency distribution):即變量小于某一上限值的出現(xiàn)次數(shù)與總次數(shù)之比[17]??梢员硎緸椋?/p>

    (1)

    (2)

    式中:n——取值范圍內(nèi)的數(shù)值等級數(shù);fi——極端事件在第i個數(shù)值等級內(nèi)變量發(fā)生的頻數(shù);Fi——變量小于等于某上限值的發(fā)生頻數(shù)。若變量為日流量時(shí),則將日流量累積頻率達(dá)到的概率分布所對應(yīng)的流量臨界值定義為極端水文事件的閾值,當(dāng)日流量大于或小于極端水文事件的閾值時(shí),認(rèn)為該日發(fā)生極端水文事件。

    極端水文事件確定的具體步驟如下:首先,確定洪水事件和枯水事件。根據(jù)《中國水資源評價(jià)》所定的三級標(biāo)準(zhǔn),將日流量資料按從小到大的順序排列,正常年相應(yīng)頻率為37.5%≤P≤62.5%,豐水年相應(yīng)頻率為P≥62.5%,枯水年相應(yīng)頻率為P≤37.5%[18]。因此,累積頻率為40%~60%的流量是河道常年發(fā)生頻次最多的一個流量范圍,故將洪水事件和枯水事件的閾值分別定為累積頻率為60%和40%時(shí)所對應(yīng)的流量。若日流量大于洪水事件的閾值時(shí),就記為一次洪水事件;若日流量小于枯水事件的閾值時(shí),就記為一次枯水事件。然后,確定極端洪水事件和極端枯水事件。極端事件即為統(tǒng)計(jì)意義上的小概率事件,其極端程度一般對應(yīng)于90%或10%的累積頻率[2]。故將洪水事件的日流量按從小到大排列,當(dāng)累積頻率達(dá)到90%時(shí),其所對應(yīng)的流量值記為極端洪水事件的閾值;將枯水事件的日流量按從小到大排列,當(dāng)累積頻率達(dá)到10%時(shí),其所對應(yīng)的流量值記為極端枯水事件的閾值。根據(jù)上述定義,計(jì)算得出開都河大山口站洪水事件的閾值為106.26m3/s,極端洪水事件的閾值則為276m3/s;枯水事件的閾值為74.5m3/s,極端枯水事件的閾值則為38.3m3/s。

    2.2.2 集中度與集中期 集中度與集中期的計(jì)算原理和方法為,將一年內(nèi)365天折算為0°~360°圓周角,并按特定的角度將各月極端水文事件頻次分解為水平(A)和垂直(B)兩個方向上的分量,然后分別進(jìn)行水平或垂直方向上的累加,合成矢量向量(C)的方位角(θ)定義為集中期(D),其所對應(yīng)的日期定義為集中日期,合成向量的模與年極端水文事件總頻次的比值定義為極端水文事件集中度(Cn)[10,19]。集中期能夠反映極端水文事件具體集中時(shí)間,集中度則能更好地反映出極端水文事件在年內(nèi)的集中程度。向量夾角及全年每個月代表的角度值如下(圖2)所示。計(jì)算公式見式(3)至式(5)。

    (3)

    D=tan-1(Rx/Ry)

    (4)

    (5)

    式中:Cn——集中度;D——集中期;R——極端水文事件總頻次的向量;Rx——x方向上的合成向量;Ry——y方向的合成向量;W——極端水文事件的總頻次;θi——各月中代表的角度值(其中,i=1,2,3,…,12);ri——第i月極端水文事件頻次的向量。極端水文事件向量夾角及所對應(yīng)的特定角度值θi的具體分布情況如圖1所示。

    圖1 向量夾角及全年每個月代表的角度值(θi)

    2.2.3 輪次分析 輪次分析在水文中的應(yīng)用即將一個連續(xù)的徑流量序列Qt(t=1,2,…,n),被所定義的截取水平Q0所截取,在一個或多個時(shí)間段內(nèi),當(dāng)徑流量Qt依次小于截取水平Q0時(shí),則出現(xiàn)負(fù)輪;反之則出現(xiàn)正輪[20]。負(fù)輪對應(yīng)于枯水事件的研究。采用輪次分析法表示枯水事件的3個指標(biāo)的表達(dá)式見式(6)至式(8)。

    T=n

    (6)

    (7)

    L=max(Q0-Qt) (t=1,2,…,n)

    (8)

    式中:T——枯水歷時(shí),指一個負(fù)輪的長度,表示枯水期持續(xù)的時(shí)間長度;S——枯水總量,指在一個負(fù)輪內(nèi),日流量與閾值相比之差的總和,表示枯水期缺乏水量的總體狀況;L——枯水極值,指在一個負(fù)輪內(nèi),日流量與閾值相比之差的最大值,表示枯水期內(nèi)最嚴(yán)重缺水日的缺水量大小。

    3 極端水文事件特征變化分析

    3.1 年內(nèi)分布特征

    開都河來水年內(nèi)變化很大,4—9月為豐水季節(jié),10月至翌年3月為枯水季節(jié)[21]。通過對極端水文事件頻率的年內(nèi)分析(圖2),可以看出:極端洪水發(fā)生頻率最多的日數(shù)主要集中在5—8月,約占全年的95.6%,最大月出現(xiàn)在7月,約占全年的36.8%。春季(3—5月)約占10.4%,夏季(6—8月)約占88.4%,秋季(9—11月)約占1.3%,冬季(12—2月)為0;極端枯水發(fā)生頻率最多的日數(shù)主要集中在12—3月,約占全年的93.5%,最大月出現(xiàn)在2月,約占全年的31.2%。春季(3—5月)約占20.4%,夏季(6—8月)約占0.2%,秋季(9—11月)約占4.3%,冬季(12—2月)約占75.1%。

    圖2 極端水文事件頻率年內(nèi)分布

    為近一步量化極端洪水與極端枯水事件的年內(nèi)分布情況,借鑒集中度與集中期方法[22]和不均勻系數(shù)[23],進(jìn)一步分析年內(nèi)極端洪水和極端枯水日數(shù)的變化特征。根據(jù)公式計(jì)算集中度、集中期及不均勻系數(shù),計(jì)算結(jié)果見表1。

    表1 開都河流域極端水文事件年內(nèi)變化統(tǒng)計(jì)特征

    從表1可以看出,極端洪水與極端枯水事件的不

    均勻系數(shù)分別為63.35%和60.19%,不均勻程度較高,集中度分別為86.28%和82.03%,集中程度也較高,比較不均勻系數(shù)與集中度可以看出,兩者具有較好的同步性變化規(guī)律,凡是集中度數(shù)值高的不均勻系數(shù)也高。極端洪水的集中期在7月10日,極端枯水的集中期在2月1日,這與圖1的分析結(jié)果一致。可以看出,極端洪水和極端枯水屬于季節(jié)性極端水文事件,年內(nèi)分布不均,且集中程度較高。

    3.2 年際變化分析

    3.2.1 極端洪水事件年際變化 對極端洪水年際頻次距平分析可以看出(圖3a),極端洪水事件在1972—2002年呈上升趨勢,2002—2007年呈急劇下降趨勢,2007年后又有略微上升的趨勢。經(jīng)過統(tǒng)計(jì),多年平均值為14 d,發(fā)生頻次最多的是2002年,正距平67 d,其次為2000年,正距平49 d。1994年之前,負(fù)距平較多,1994—2012年,正距平明顯增多,說明1994—2012年期間,極端洪水發(fā)生率顯著增加。極端洪水頻次的變差系數(shù)為1.30,表明極端洪水年際變化較大;偏態(tài)系數(shù)為2.04,表明極端洪水年際變化偏離正態(tài)分布的程度較高,具有極端不穩(wěn)定的特性(表2)。對極端洪水的年代際分析,可以得出:20世紀(jì)80年代年均頻次較70年代增長了64%,90年代較80年代增長了115.9%,2000—2009年較90年代增長了22.6%,可見90年代增長顯著,極端洪水事件年代際間表現(xiàn)出增長的趨勢。

    圖3 極端洪水與極端枯水頻次的年際距平及5年滑動平均

    對極端洪水頻次做M-K突變檢驗(yàn)(圖4a),可以看出,正向序列UF和反向序列UB在0.01顯著性水平下在信度線內(nèi)交于1點(diǎn),且UF超過信度線,可以確定該點(diǎn)所對應(yīng)的1986年為突變年份。這是因?yàn)闅鉁卦?0年代以前低于距平值,80年代以來明顯高于距平值[10],氣溫升高導(dǎo)致冰雪融水增多,降水量增多,又因?yàn)閺搅髟龆嘤幸欢ǖ臏笞饔?,因此?986年發(fā)生突變。并且,UF在1997—2003年間及2006年之后逐漸超過了信度線,說明極端洪水事件在這些時(shí)段內(nèi)呈顯著增加趨勢。

    3.2.2 極端枯水事件年際變化 對極端枯水年際頻次距平分析可以看出(圖3b),極端枯水事件在1972—1978年呈上升趨勢,1978—1983年呈下降趨勢,1983—1986年呈上升趨勢,1986—2002年呈下降趨勢,2002—2012年呈上升趨勢。多年平均值為13.85 d,發(fā)生頻次最多的是1975年,正距平65 d,其次為1986年,正距平52 d。從圖中明顯可以看出在1975年前后和1985—1989年間發(fā)生極端枯水的頻次較為頻繁。1989年之后基本上都是負(fù)距平,說明1989年之后發(fā)生極端枯水事件的頻次較少。極端枯水頻次的變差系數(shù)為1.30,表明極端枯水年際變化較大;偏態(tài)系數(shù)2.06,表明極端枯水年際變化偏離正態(tài)分布的程度較高,具有極端不穩(wěn)定的特性(表2)。對極端枯水的年代際分析,可以得出,20世紀(jì)80年代年均頻次較70年代減少了6.8%,90年代較80年代減少了63.9%,2000—2009年較90年代減少了18.9%,可見90年代極端枯水事件減少顯著,極端枯水事件年代際間表現(xiàn)出減少的趨勢。

    圖4 極端洪水與極端枯水頻次突變檢驗(yàn)

    對極端枯水頻次做M-K突變檢驗(yàn)(圖4b),可以看出,正向序列UF和反向序列UB在0.01顯著性水平下在信度線內(nèi)交于2點(diǎn),分別為1993年和2007年。由于UF或UB在2007年后并未超過信度線,因此可以排除2007年為突變年份,確定1993年為極端枯水事件的突變點(diǎn)。這與氣候變化下氣溫升高、降水量增多導(dǎo)致徑流量變化有關(guān)。從圖中可以注意到,在1998—2001年UB超過了信度線,說明極端枯水事件在此期間減少顯著。

    表2 開都河流域極端水文事件年際變化統(tǒng)計(jì)特征

    3.3 極端洪水場次分析

    一場洪水要有三要素才能將這場洪水表示完整,洪水三要素包括:洪峰流量、洪量和洪水歷時(shí)[20]。將年內(nèi)徑流量連續(xù)超過極端洪水閾值的時(shí)段記為一場,經(jīng)過統(tǒng)計(jì),得出1972—2012年共發(fā)生極端洪水事件110起,發(fā)生次數(shù)最多的是2008年,發(fā)生了8起;1974年、1977年、1981年、1986年和2003年未發(fā)生極端洪水事件,說明極端洪水事件具有極端不穩(wěn)定性的特征。1972—1999年最大洪峰流量呈增長趨勢,1999—2003年呈減少趨勢,2003年之后又呈增長趨勢。最大洪峰流量793 m3/s,發(fā)生在1999年7月,年最大洪峰流量超過500 m3/s的極端洪水事件有6 a,主要集中在1999—2002年間。對最大洪峰流量的M-K突變檢驗(yàn)結(jié)果表明,洪峰流量在1991年發(fā)生突變。年最大洪水歷時(shí)與年最大洪峰流量基本呈對應(yīng)關(guān)系,年最大洪水歷時(shí)也主要集中在1999—2002年間。最大洪水歷時(shí)80 d,發(fā)生在2002年6—8月,年最大洪水歷時(shí)超過15 d的有8 a。其年際分布如圖5所示。

    圖5 極端洪水事件年最大洪峰流量和最大洪水歷時(shí)

    3.4極端枯水場次分析

    對枯水事件的研究常采用輪次分析法,本文從枯水總量、枯水極值和枯水歷時(shí)三個方面來對枯水事件進(jìn)行表示。對極端枯水事件的統(tǒng)計(jì)得出:1972—2012年間,極端枯水事件共發(fā)生149起,發(fā)生次數(shù)最多的是2008年,發(fā)生了12起,1972年、1980年、1981年、1983年、1988年、1999年、2001年和2002年這8年未發(fā)生極端枯水事件,說明極端枯水事件具有極端不穩(wěn)定性的特征。

    輪次分析結(jié)果表明:2007年11月5日極端枯水極值最大,缺水326.6萬m3。最大枯水歷時(shí)36 d,發(fā)生在1975年2—3月,年最大歷時(shí)超過15 d的有9 a。最大枯水總量為293.24萬m3,發(fā)生在2007年,其次為1985年,枯水總量為1 485.2萬m3,最大枯水總量主要集中在1974—1979年和1984—1986年這兩個時(shí)段內(nèi),分別約占枯水總量的20.2%和21.3%。其年際變化如圖6所示。

    圖6 極端枯水事件年最大枯水總量和最大枯水歷時(shí)

    4 極端水文事件對氣候變化的響應(yīng)

    4.1極端洪水事件對氣候變化的響應(yīng)

    巴音布魯克站位于開都河上游中心位置,其資料可反映流域上游氣候狀況。對巴音布魯克站1972—2012年年均氣溫和年均降水量與極端洪水頻次的關(guān)系(圖7)研究表明:年均氣溫與極端洪水頻次的相關(guān)系數(shù)為0.05,兩者相關(guān)性不大;年降水量年際變化趨

    勢與極端洪水頻次變化趨勢一致,兩者的相關(guān)系數(shù)為0.64,并通過了0.01的顯著性檢驗(yàn),表明降水量對極端洪水事件的貢獻(xiàn)較大。降水量的增加使得地表面水來不及下滲便匯流成河,是導(dǎo)致發(fā)生極端洪水事件的直接因素,而氣溫導(dǎo)致冰雪融水補(bǔ)給河流的比例少,且使河流水量緩慢增加。

    對年內(nèi)極端洪水頻次與氣溫平均值和降水平均值的關(guān)系研究(圖8)表明:極端洪水頻次與氣溫平均值和降水平均值的年內(nèi)變化趨勢一致,表現(xiàn)為1—7月的增長趨勢,7—12月的下降趨勢,均在7月份達(dá)到氣溫、降水和極端洪水頻次的最大值。極端洪水頻次與氣溫平均值的相關(guān)系數(shù)為0.70,通過了0.05的顯著性水平檢驗(yàn);與降水平均值的相關(guān)系數(shù)為0.96,通過了0.01的顯著性水平檢驗(yàn)。因?yàn)殚_都河流域降水量補(bǔ)給比例大,夏季降水多,極端洪水頻次多,冬季降水以固體形式存在,補(bǔ)給量少,發(fā)生極端洪水的頻次少。因此,關(guān)注6—8月的氣溫、降水變化對于預(yù)防極端洪水災(zāi)害至關(guān)重要。

    圖7 極端洪水頻次與年氣溫和年降水量變化關(guān)系

    圖8 年內(nèi)極端洪水頻次與氣溫平均值和降水平均值的變化關(guān)系

    4.2極端枯水事件對氣候變化的響應(yīng)

    對年際極端枯水頻次與年均氣溫和年均降水量的關(guān)系研究(圖9)表明:年際極端枯水頻次與年均氣溫和年均降水總體上呈反相關(guān)。極端枯水頻次與年均氣溫的相關(guān)系數(shù)-0.34,并通過了0.05的顯著性檢驗(yàn);極端枯水頻次與年均降水的相關(guān)系數(shù)為-0.30,相關(guān)性不顯著??梢姡瑲鉁貙O端枯水事件的影響較大。氣溫越高,降水與融雪徑流越多,發(fā)生極端枯水事件的可能性越小。

    對年內(nèi)極端枯水頻次與氣溫平均值和降水平均值的關(guān)系研究(圖10)表明:年內(nèi)極端枯水頻次與氣溫平均值和降水平均值總體上亦呈反相關(guān),6—8月氣溫和降水達(dá)到最大時(shí),發(fā)生極端枯水事件的頻次為0,11月到次年3月氣溫和降水量減少時(shí),發(fā)生極端枯水事件的頻次增多。極端枯水頻次與氣溫平均值的相關(guān)系數(shù)為-0.89,并通過了0.01的顯著性檢驗(yàn);與降水平均值的相關(guān)系數(shù)為-0.58,并通過了0.05的顯著性檢驗(yàn)。表明氣溫對極端枯水事件的影響最大,其次為降水量的影響。這是因?yàn)闅鉁卦谀陜?nèi)各月均增高,而降水主要是夏季增加的比較多,冬季主要以雪等固體形式存在,補(bǔ)給量少,因此極端枯水與氣溫的反相關(guān)比較顯著。

    圖9 極端枯水頻次與年氣溫和年降水量變化關(guān)系

    圖10 年內(nèi)極端枯水頻次與氣溫平均值和降水平均值的變化關(guān)系

    5 結(jié) 論

    本文主要研究了開都河流域過去40 a極端水文事件頻次的年內(nèi)年際變化、場次分布特征以及對區(qū)域氣候變化的響應(yīng)。研究結(jié)果表明該流域極端水文事件具有極端不穩(wěn)定的特性。極端事件頻次和場次的分析結(jié)果一致,表明極端洪水事件主要發(fā)生在7月,極端枯水事件主要發(fā)生在2月,且年內(nèi)分配不均勻。1972—2012年,極端洪水事件總體呈現(xiàn)增長趨勢,且在1986年發(fā)生突變;極端枯水事件總體呈現(xiàn)減少趨勢,且在1993年發(fā)生突變。隨著氣溫的升高,降水和融雪徑流增加,極端洪水事件增加,極端枯水事件減少。極端水文事件與氣溫、降水的關(guān)系研究表明流域近40 a的氣溫是增加的,降水量的增加是導(dǎo)致極端洪水發(fā)生的主要因子。極端水文事件引起的自然災(zāi)害已經(jīng)威脅到開都河流域農(nóng)牧業(yè)、生產(chǎn)、交通安全和水資源供給[22-23],應(yīng)加強(qiáng)氣候變化背景下的水文水資源安全及應(yīng)對措施研究,提高水資源安全保障,減緩氣候變化帶來的可能危害。

    [1] IPCC.Climate change 2007: The physical science basis[M]∥Contribution of Working Group 1 to the Fourth Assessment Report of the Intergovermental Panel on Climate Change Cambridge,United Kingdom and New York,USA: Cambridge University Press,2007.

    [2] 張利平,杜鴻,夏軍,等.氣候變化下極端水文事件的研究進(jìn)展[J].地理科學(xué)進(jìn)展,2011,30(11):1370-1379.

    [3] Akhtar M, Ahmad N, Booij M J. The impact of climate change on the water resources of Hindukush-Karakorum-Himalaya region under different glacier coverage scenarios[J]. Journal of Hydrology,2008,335(1):148-163.

    [4] 趙慶云,張武,王式功,等.西北地區(qū)東部干旱半干旱區(qū)極端降水事件的變化[J].中國沙漠,2005,25(6):904-909.

    [5] 張明軍,汪寶龍,魏軍林,等.近50年寧夏極端氣溫事件的變化研究[J].自然災(zāi)害學(xué)報(bào),2012,21(4):152-160.

    [6] Hong Y, Adler R F. Predicting global landslide spatio temporal distribution: Integrating landslide susceptibility zoning techniques and real-time satellite rainfall estimates[J]. International Journal of Sediment Research,2008,23(3):249-257.

    [7] 李志,劉文兆,鄭粉莉.1965~2005年涇河流域極端溫度事件變化特征[J].地理科學(xué),2010,30(3):469 -474.

    [8] 夏軍,佘敦先,杜鴻.氣候變化影響下極端水文事件的多變量統(tǒng)計(jì)模型研究[J].氣候變化研究進(jìn)展,2012,8(6):397-401.

    [9] 孫桂麗,陳亞寧,李衛(wèi)紅,等.新疆極端水文事件的時(shí)空分布特征[J].自然災(zāi)害學(xué)報(bào),2012,21(3):119-124.

    [10] 孫桂麗,陳亞寧,李衛(wèi)紅,等.新疆極端水文事件年際變化及對氣候變化的響應(yīng)[J].地理科學(xué),2011,31(11):1389-1394.

    [11] 謝平,李析男,許斌,等.基于希爾伯特—黃變換的非一致性洪水頻率計(jì)算方法:以西江大湟江口站為例[J].自然災(zāi)害學(xué)報(bào),2013,22(1):86-91.

    [12] 方彬,武見,肖義,等.混合分布在枯水流量頻率分析中的應(yīng)用研究[J].石河子大學(xué)學(xué)報(bào):自然科學(xué)版,2006,24(2):231-233.

    [13] 鐘瑞森,董新光.新疆博斯騰湖水鹽平衡及水環(huán)境預(yù)測[J].湖泊科學(xué),2008,20(1):58-64.

    [14] 施雅風(fēng).簡明中國冰川目錄[M].上海:上??茖W(xué)技術(shù)普及出版社,2005.

    [15] 藍(lán)永超,鐘英君,吳素芬,等.天山南—北坡典型河流出山徑流對氣候變化響應(yīng)的分析對比[J].中國沙漠,2011,31(1):254-259.

    [16] 鄧銘江,王世江,董新光,等.新疆水資源及可持續(xù)利用[M].北京:中國水利水電出版社,2005:150-167.

    [17] 黃俊利,魏文壽,楊青,等.塔克拉瑪干沙漠南緣極端降水變化趨勢與突變特征[J].沙漠與綠洲氣象,2012,6(2):30-34.

    [18] 陶輝,宋郁東,鄒世平,等.開都河天山出山徑流量年際變化特征與洪水頻率分析[J].干旱區(qū)地理,2007,30(1):43-48.

    [19] 劉賢趙,李嘉竹,宿慶,等.基于集中度與集中期的徑流年內(nèi)分配研究[J].地理科學(xué),2007,27(6):791-695.

    [20] 楊帆.區(qū)域氣候背景下流域極端氣候與水文事件的關(guān)系研究[D].鄭州:鄭州大學(xué),2011.

    [21] 鐘瑞森,吳彬.開都河徑流年內(nèi)分配特征研究[J].灌溉排水學(xué)報(bào),2011,30(2):123-125.

    [22] 賀斌,王國亞,蘇宏超,等.新疆阿爾泰山地區(qū)極端水文事件對氣候變化的響應(yīng)[J].冰川凍土,2012,34(4):927-932.

    [23] 鄭祚芳,張秀麗.北京極端天氣事件及其與區(qū)域氣候變化的聯(lián)系[J].自然災(zāi)害學(xué)報(bào),2007,16(3):55-593.

    ExtremeHydrologicalEventsandItsResponsetoClimateChangeintheKaiduRiverBasin

    LI Jia-xiu1,2, XU Chang-chun1,2, WANG Xiao3, GAO Shen-tong1,2, ZHAO Jie1,2

    (1.KeyLaboratoryofOasisEcologyofMinistryofEducation,XinjiangUniversity,Urumqi830046,China; 2.SchoolofResourcesandEnvironmentalScience,XinjiangUniversity,Urumqi830046,China; 3.DepartmentofPersonnel,SichuanUniversityforNationalities.Kangding,Sichuan626001,China)

    Based on the daily discharge and meteorological date from 1972 to 2012 at the Dashankou hydrological station and Bayinbuluke weather station in the Kaidu River Basin, the characteristics of extreme hydrological events (EHE) were investigated, the threshold value was determined by the cumulative frequency method, the inner-annual and inter-annual changes of EHE, the session distribution of EHE and their responses to climate change were analyzed in the study area. The results showed that the EHE in the region was extremely unstable. The inner-annual distribution was uneven and the concentration degree was higher. The extreme flood events mainly occurred in July, and the extreme low water events mainly appeared in February. The inter-annual variability was quite unstable with a high deviation from the normal distribution. The extreme flood events totally presented an increasing trend and had an abrupt change in 1986; while the extreme low water events totally presented a decreasing trend and had an abrupt change in 1993.The extreme flood frequency had a higher correlation with precipitation, while the extreme low water frequency had a higher negative correlation with temperature.

    extreme flood events; extreme low water events; cumulative frequency method; concentration ratio; runs arithmetics

    2013-11-13

    :2013-12-08

    國家自然科學(xué)基金“基于大氣環(huán)流分型的西北干旱區(qū)極端氣候事件研究”(41305125);國家自然科學(xué)基金“過去500年新疆北部水文變化及其氣候驅(qū)動”(U1203281)

    李佳秀(1988—),女,陜西渭南人,碩士研究生,主要從事極端水文和極端氣候事件研究。E-mail:1428207461@qq.com

    徐長春(1977—),女,江蘇贛榆人,副教授,博士,主要從事干旱區(qū)氣候、水文與水資源研究。E-mail:xcc0110@163.com

    P467;P333

    :A

    :1005-3409(2014)05-0285-07

    猜你喜歡
    開都河枯水年際
    土地利用結(jié)構(gòu)與景觀格局對開都河水質(zhì)的影響
    哲理漫畫
    丹棱縣晚熟桔橙枯水現(xiàn)狀及防控技術(shù)分析
    耕作與栽培(2021年3期)2021-12-01 02:38:13
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    開都河流域新疆裸重唇魚胚胎發(fā)育觀察
    柑橘枯水研究進(jìn)展
    基于逐步聚類分析法的開都河徑流預(yù)測模型
    人民黃河(2019年2期)2019-09-10 07:22:44
    開都河諸小河流水資源管理存在的問題及建議
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    妹子高潮喷水视频| av有码第一页| 黄色视频在线播放观看不卡| 亚洲av中文av极速乱| 欧美一级a爱片免费观看看| 大码成人一级视频| 免费看日本二区| freevideosex欧美| 精品亚洲成a人片在线观看| 在线 av 中文字幕| 大又大粗又爽又黄少妇毛片口| 一级毛片久久久久久久久女| 婷婷色综合大香蕉| 亚洲国产精品国产精品| 日本欧美国产在线视频| 日韩不卡一区二区三区视频在线| 欧美激情国产日韩精品一区| 女性被躁到高潮视频| 亚洲精品第二区| 天堂8中文在线网| 一区二区av电影网| 久久久久久久久久久免费av| 国语对白做爰xxxⅹ性视频网站| 亚洲第一区二区三区不卡| 特大巨黑吊av在线直播| 欧美 亚洲 国产 日韩一| av免费观看日本| 日韩,欧美,国产一区二区三区| 少妇熟女欧美另类| 国产精品女同一区二区软件| av福利片在线| 波野结衣二区三区在线| 国产深夜福利视频在线观看| 啦啦啦视频在线资源免费观看| 久久99蜜桃精品久久| 久久久国产欧美日韩av| 亚洲婷婷狠狠爱综合网| 久久精品久久久久久噜噜老黄| 蜜桃在线观看..| 最近的中文字幕免费完整| 在线观看国产h片| av专区在线播放| 插阴视频在线观看视频| 九九久久精品国产亚洲av麻豆| 内地一区二区视频在线| 看非洲黑人一级黄片| 人人妻人人看人人澡| 国精品久久久久久国模美| 亚洲av欧美aⅴ国产| 日韩成人伦理影院| 少妇猛男粗大的猛烈进出视频| 国产伦理片在线播放av一区| 麻豆成人午夜福利视频| 人人妻人人添人人爽欧美一区卜| 国产精品偷伦视频观看了| 国产爽快片一区二区三区| 视频中文字幕在线观看| 黄色一级大片看看| 一本一本综合久久| 免费观看性生交大片5| 成人国产av品久久久| av国产久精品久网站免费入址| 精品99又大又爽又粗少妇毛片| 多毛熟女@视频| 亚洲av欧美aⅴ国产| 欧美成人精品欧美一级黄| 日本av手机在线免费观看| 制服丝袜香蕉在线| 丝瓜视频免费看黄片| 美女xxoo啪啪120秒动态图| 亚洲天堂av无毛| 久久久久久久亚洲中文字幕| 多毛熟女@视频| 亚洲国产精品一区二区三区在线| 亚洲国产欧美日韩在线播放 | 99久久精品国产国产毛片| 久久久久久久久大av| 欧美97在线视频| av国产久精品久网站免费入址| 成人黄色视频免费在线看| 国产成人一区二区在线| 女性被躁到高潮视频| av播播在线观看一区| 黄色配什么色好看| 午夜福利影视在线免费观看| 国产黄片美女视频| 在线观看人妻少妇| 中国国产av一级| 中国美白少妇内射xxxbb| 国产欧美另类精品又又久久亚洲欧美| 不卡视频在线观看欧美| 成年人午夜在线观看视频| 亚洲美女视频黄频| 伊人久久精品亚洲午夜| 亚洲精品中文字幕在线视频 | 午夜影院在线不卡| 高清欧美精品videossex| 久久久久久久精品精品| 久久99精品国语久久久| 国产又色又爽无遮挡免| 久热久热在线精品观看| 性色avwww在线观看| 少妇丰满av| 免费高清在线观看视频在线观看| 你懂的网址亚洲精品在线观看| 99久久人妻综合| 韩国av在线不卡| 久久国产亚洲av麻豆专区| 国产男人的电影天堂91| 国产精品偷伦视频观看了| 全区人妻精品视频| 国产色爽女视频免费观看| 色视频在线一区二区三区| 亚洲熟女精品中文字幕| 边亲边吃奶的免费视频| 又黄又爽又刺激的免费视频.| 美女cb高潮喷水在线观看| 精品一区在线观看国产| 少妇人妻久久综合中文| 男女边摸边吃奶| 亚洲国产精品999| 成年av动漫网址| av卡一久久| 日本av免费视频播放| 人人妻人人看人人澡| 狠狠精品人妻久久久久久综合| 日韩欧美精品免费久久| 男女国产视频网站| 69精品国产乱码久久久| 在线观看av片永久免费下载| 国产伦理片在线播放av一区| 国产成人精品一,二区| 中文字幕亚洲精品专区| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品自拍成人| 亚洲精品aⅴ在线观看| 免费看不卡的av| 秋霞伦理黄片| 日韩欧美一区视频在线观看 | 国产 精品1| 亚洲国产av新网站| 中文字幕人妻丝袜制服| 观看免费一级毛片| 少妇人妻精品综合一区二区| 男女边摸边吃奶| 国产白丝娇喘喷水9色精品| 亚洲精华国产精华液的使用体验| 新久久久久国产一级毛片| 国产精品久久久久久精品电影小说| 夜夜骑夜夜射夜夜干| 国产av国产精品国产| 精品久久久久久久久av| 最新中文字幕久久久久| 大陆偷拍与自拍| 日韩电影二区| 午夜视频国产福利| 国产成人精品一,二区| 人妻系列 视频| 午夜久久久在线观看| 国产免费福利视频在线观看| 在线精品无人区一区二区三| av黄色大香蕉| 丰满饥渴人妻一区二区三| 大片免费播放器 马上看| 国产精品女同一区二区软件| 一级a做视频免费观看| 搡女人真爽免费视频火全软件| 欧美亚洲 丝袜 人妻 在线| .国产精品久久| 欧美3d第一页| 久久久久精品久久久久真实原创| 中文字幕久久专区| 亚洲国产精品国产精品| 国产精品国产三级国产专区5o| av国产久精品久网站免费入址| 欧美日韩视频高清一区二区三区二| 大码成人一级视频| 久久久久视频综合| 69精品国产乱码久久久| 人妻夜夜爽99麻豆av| 黄色日韩在线| av线在线观看网站| 女的被弄到高潮叫床怎么办| 如何舔出高潮| 狂野欧美激情性bbbbbb| 成人漫画全彩无遮挡| 日本wwww免费看| 久久久午夜欧美精品| 永久网站在线| 黄色日韩在线| 国产精品一区二区在线观看99| 国产欧美亚洲国产| 国产白丝娇喘喷水9色精品| 亚洲av成人精品一区久久| 午夜免费鲁丝| 777米奇影视久久| 99热网站在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美bdsm另类| 久久久午夜欧美精品| 寂寞人妻少妇视频99o| 久久久久久久久久人人人人人人| 午夜视频国产福利| 日韩欧美一区视频在线观看 | 日本黄色日本黄色录像| 男女无遮挡免费网站观看| 国产91av在线免费观看| 少妇猛男粗大的猛烈进出视频| 亚洲精品久久午夜乱码| 国产在视频线精品| 一级毛片电影观看| 久久99精品国语久久久| 日本黄大片高清| av.在线天堂| 精品酒店卫生间| 国产高清三级在线| 97在线人人人人妻| 韩国高清视频一区二区三区| 2021少妇久久久久久久久久久| 亚洲欧洲日产国产| 久久久国产一区二区| 成人国产麻豆网| 国产一区二区在线观看日韩| 免费av不卡在线播放| 在线天堂最新版资源| 在线观看免费视频网站a站| 又黄又爽又刺激的免费视频.| 亚洲av不卡在线观看| 91精品国产九色| 日韩亚洲欧美综合| 少妇猛男粗大的猛烈进出视频| 日本猛色少妇xxxxx猛交久久| 色婷婷av一区二区三区视频| 国产成人91sexporn| 99久久精品热视频| 久久久久精品久久久久真实原创| 国产熟女午夜一区二区三区 | 少妇的逼好多水| 大码成人一级视频| 爱豆传媒免费全集在线观看| 久久国产乱子免费精品| 夜夜看夜夜爽夜夜摸| 欧美激情极品国产一区二区三区 | 久久国内精品自在自线图片| 中文字幕av电影在线播放| 秋霞在线观看毛片| 黄色欧美视频在线观看| 特大巨黑吊av在线直播| 晚上一个人看的免费电影| 晚上一个人看的免费电影| 国产一级毛片在线| 久久99热6这里只有精品| 色视频www国产| 国产69精品久久久久777片| 建设人人有责人人尽责人人享有的| 下体分泌物呈黄色| 国产老妇伦熟女老妇高清| 亚洲精品色激情综合| 亚洲av福利一区| 成人18禁高潮啪啪吃奶动态图 | 国产亚洲午夜精品一区二区久久| 亚洲国产精品一区三区| 日本av手机在线免费观看| 国产亚洲欧美精品永久| 亚州av有码| 精品国产一区二区久久| 自拍欧美九色日韩亚洲蝌蚪91 | 一级二级三级毛片免费看| 偷拍熟女少妇极品色| 啦啦啦视频在线资源免费观看| 国产精品熟女久久久久浪| 欧美3d第一页| 丰满少妇做爰视频| 亚洲精品乱久久久久久| 少妇被粗大猛烈的视频| 国产伦精品一区二区三区视频9| 精品久久久久久久久亚洲| 少妇丰满av| 色网站视频免费| 一级毛片黄色毛片免费观看视频| 制服丝袜香蕉在线| 国产一区二区三区av在线| 亚洲一区二区三区欧美精品| 国产极品粉嫩免费观看在线 | 丰满少妇做爰视频| 日韩免费高清中文字幕av| 日日摸夜夜添夜夜爱| 少妇猛男粗大的猛烈进出视频| 美女cb高潮喷水在线观看| 亚洲av综合色区一区| 亚洲第一av免费看| 欧美97在线视频| 国产精品伦人一区二区| 国产精品一二三区在线看| 亚洲精品456在线播放app| 欧美成人午夜免费资源| 国产成人aa在线观看| 国产精品99久久99久久久不卡 | 在线观看三级黄色| 伦理电影免费视频| 天堂中文最新版在线下载| 老熟女久久久| 成年人午夜在线观看视频| av女优亚洲男人天堂| 国产无遮挡羞羞视频在线观看| 丰满饥渴人妻一区二区三| 内地一区二区视频在线| 午夜影院在线不卡| 欧美老熟妇乱子伦牲交| 国产黄片美女视频| 免费观看av网站的网址| 伊人亚洲综合成人网| 欧美精品一区二区免费开放| 精品久久久久久久久亚洲| 自拍欧美九色日韩亚洲蝌蚪91 | 一级二级三级毛片免费看| 久久6这里有精品| 久久精品夜色国产| 在线免费观看不下载黄p国产| 韩国av在线不卡| 一级黄片播放器| a级毛色黄片| 亚洲美女搞黄在线观看| 精品一区二区三卡| 免费黄频网站在线观看国产| 久久久久久人妻| 欧美老熟妇乱子伦牲交| 日韩欧美 国产精品| 天天躁夜夜躁狠狠久久av| 国产精品一区二区在线观看99| 十八禁网站网址无遮挡 | 成人国产麻豆网| 久久久久久伊人网av| 国产亚洲最大av| 国产日韩一区二区三区精品不卡 | av国产久精品久网站免费入址| 简卡轻食公司| 99热这里只有精品一区| 秋霞伦理黄片| 亚洲av二区三区四区| 亚洲欧美精品自产自拍| 精品国产一区二区三区久久久樱花| 精品人妻偷拍中文字幕| 国产成人91sexporn| 人人妻人人添人人爽欧美一区卜| 亚洲国产精品成人久久小说| 涩涩av久久男人的天堂| 国产毛片在线视频| 亚洲电影在线观看av| 亚洲av欧美aⅴ国产| 26uuu在线亚洲综合色| 一级二级三级毛片免费看| 秋霞在线观看毛片| 中文字幕制服av| 日本91视频免费播放| 人妻制服诱惑在线中文字幕| 制服丝袜香蕉在线| 欧美bdsm另类| 欧美性感艳星| 妹子高潮喷水视频| 一级毛片aaaaaa免费看小| 啦啦啦视频在线资源免费观看| av国产久精品久网站免费入址| 一本一本综合久久| 国产永久视频网站| 男女边摸边吃奶| 69精品国产乱码久久久| 成人亚洲精品一区在线观看| 亚洲一级一片aⅴ在线观看| 亚洲精品中文字幕在线视频 | 亚洲精品,欧美精品| 人人妻人人澡人人看| 日韩大片免费观看网站| 国产精品人妻久久久影院| 精品一区在线观看国产| 少妇猛男粗大的猛烈进出视频| 亚洲欧美一区二区三区国产| 亚洲精品成人av观看孕妇| 日韩成人伦理影院| 久久久久久久久久久免费av| 日本av免费视频播放| 欧美日韩综合久久久久久| 婷婷色麻豆天堂久久| 26uuu在线亚洲综合色| 日韩强制内射视频| 最近中文字幕高清免费大全6| 国产精品嫩草影院av在线观看| a 毛片基地| 蜜桃久久精品国产亚洲av| 国产欧美日韩精品一区二区| 日韩欧美精品免费久久| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久av不卡| 麻豆乱淫一区二区| 欧美 日韩 精品 国产| 中文字幕久久专区| 黄色一级大片看看| 久久久午夜欧美精品| 国产欧美另类精品又又久久亚洲欧美| 女人精品久久久久毛片| 免费大片18禁| 久久久久网色| 亚洲精品国产av蜜桃| 99热这里只有是精品在线观看| 内地一区二区视频在线| 成年女人在线观看亚洲视频| 中国美白少妇内射xxxbb| 国产女主播在线喷水免费视频网站| 99热网站在线观看| 噜噜噜噜噜久久久久久91| 国产高清有码在线观看视频| 成人特级av手机在线观看| 久久青草综合色| 亚洲精品久久午夜乱码| 性高湖久久久久久久久免费观看| 国产免费福利视频在线观看| 久久久精品免费免费高清| 一级毛片黄色毛片免费观看视频| 婷婷色麻豆天堂久久| 欧美日韩综合久久久久久| 日韩av在线免费看完整版不卡| 99热这里只有精品一区| 国产毛片在线视频| 大陆偷拍与自拍| 久久狼人影院| 久久av网站| 亚洲精品亚洲一区二区| 久久久久视频综合| 九九在线视频观看精品| 日韩欧美 国产精品| 亚洲国产成人一精品久久久| 久久久久人妻精品一区果冻| 偷拍熟女少妇极品色| 国产精品熟女久久久久浪| 欧美精品亚洲一区二区| av黄色大香蕉| 国产探花极品一区二区| 亚洲国产成人一精品久久久| 亚洲四区av| 狠狠精品人妻久久久久久综合| 黄色一级大片看看| 欧美三级亚洲精品| 免费播放大片免费观看视频在线观看| a级毛片在线看网站| 欧美精品一区二区免费开放| 国产精品.久久久| 久久久国产一区二区| 欧美性感艳星| 亚洲高清免费不卡视频| 久久韩国三级中文字幕| 国产精品一二三区在线看| 国产免费又黄又爽又色| 高清午夜精品一区二区三区| 十分钟在线观看高清视频www | av在线老鸭窝| 国产日韩一区二区三区精品不卡 | 丝袜在线中文字幕| 一级av片app| 国产成人午夜福利电影在线观看| 成人毛片60女人毛片免费| 国产成人91sexporn| 亚洲精品国产色婷婷电影| 在线观看美女被高潮喷水网站| 亚洲av电影在线观看一区二区三区| 麻豆成人午夜福利视频| 天堂俺去俺来也www色官网| xxx大片免费视频| 日韩中文字幕视频在线看片| 国产精品久久久久久精品电影小说| 久久精品夜色国产| 国产精品国产三级国产专区5o| 国产欧美日韩精品一区二区| 亚洲av日韩在线播放| 九色成人免费人妻av| 一级二级三级毛片免费看| 精品久久国产蜜桃| av国产久精品久网站免费入址| 日韩不卡一区二区三区视频在线| av视频免费观看在线观看| 六月丁香七月| 精品国产一区二区久久| 亚洲一区二区三区欧美精品| 中国美白少妇内射xxxbb| 两个人免费观看高清视频 | freevideosex欧美| 97在线视频观看| 97精品久久久久久久久久精品| 精品亚洲成a人片在线观看| 国产伦在线观看视频一区| 日韩,欧美,国产一区二区三区| 国产伦在线观看视频一区| kizo精华| 美女中出高潮动态图| 久久久久久久久久久久大奶| 男女国产视频网站| 99久国产av精品国产电影| 国内揄拍国产精品人妻在线| 我的老师免费观看完整版| 欧美精品人与动牲交sv欧美| 国产一区二区在线观看av| 久久国产亚洲av麻豆专区| 国产黄片视频在线免费观看| 亚洲欧美成人综合另类久久久| 久久精品久久久久久噜噜老黄| 亚洲精品久久久久久婷婷小说| 亚洲精品aⅴ在线观看| 免费av中文字幕在线| 国产免费一区二区三区四区乱码| 国产日韩欧美在线精品| 永久免费av网站大全| 少妇人妻 视频| 国产精品一区二区在线不卡| 人妻人人澡人人爽人人| 观看美女的网站| 老司机影院成人| 国产精品99久久99久久久不卡 | 九九爱精品视频在线观看| 国产视频内射| 18禁动态无遮挡网站| 久久久久精品性色| 91精品国产九色| 少妇高潮的动态图| 一个人免费看片子| 亚洲第一av免费看| 成人毛片a级毛片在线播放| 建设人人有责人人尽责人人享有的| 五月伊人婷婷丁香| 亚洲av欧美aⅴ国产| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美成人综合另类久久久| 中文天堂在线官网| 午夜激情久久久久久久| 国产高清有码在线观看视频| 久久 成人 亚洲| 又黄又爽又刺激的免费视频.| 久久久国产精品麻豆| 搡女人真爽免费视频火全软件| 91在线精品国自产拍蜜月| 乱系列少妇在线播放| 在线观看人妻少妇| 99热这里只有精品一区| 啦啦啦啦在线视频资源| 麻豆精品久久久久久蜜桃| 日韩免费高清中文字幕av| 2021少妇久久久久久久久久久| 嫩草影院新地址| 美女福利国产在线| 建设人人有责人人尽责人人享有的| 好男人视频免费观看在线| 麻豆乱淫一区二区| 亚洲欧美成人精品一区二区| 亚洲激情五月婷婷啪啪| 久久久国产一区二区| 国产精品成人在线| 日本午夜av视频| 欧美少妇被猛烈插入视频| 久久久久久久大尺度免费视频| 欧美区成人在线视频| 欧美另类一区| 综合色丁香网| 精品国产国语对白av| 久久精品夜色国产| 久久久久久人妻| 日韩熟女老妇一区二区性免费视频| 三级国产精品欧美在线观看| 精品久久久精品久久久| 成人国产av品久久久| 美女脱内裤让男人舔精品视频| 成人亚洲精品一区在线观看| 亚洲av男天堂| 女性被躁到高潮视频| 亚洲三级黄色毛片| 交换朋友夫妻互换小说| 五月开心婷婷网| 国产精品99久久久久久久久| 欧美少妇被猛烈插入视频| 最后的刺客免费高清国语| 免费少妇av软件| 久久av网站| 久久毛片免费看一区二区三区| 美女大奶头黄色视频| 国产精品一区二区在线不卡| 亚洲精品456在线播放app| 国产日韩欧美在线精品| 22中文网久久字幕| 欧美性感艳星| 蜜桃久久精品国产亚洲av| 在线看a的网站| 久久久亚洲精品成人影院| 精品一区二区免费观看| 国产老妇伦熟女老妇高清| 日韩成人伦理影院| 精品国产一区二区久久| 久久青草综合色| 新久久久久国产一级毛片| 成年美女黄网站色视频大全免费 | 免费黄频网站在线观看国产| 伦理电影大哥的女人| 日韩三级伦理在线观看| 麻豆乱淫一区二区| 午夜福利,免费看| 日韩中字成人| 热99国产精品久久久久久7| 日韩视频在线欧美| 精华霜和精华液先用哪个| 中文字幕精品免费在线观看视频 | 99久久中文字幕三级久久日本| 久久久久久久久大av| 精品人妻熟女av久视频| 国产成人freesex在线| 高清av免费在线| 国产精品女同一区二区软件| 我要看黄色一级片免费的| 亚洲欧美一区二区三区黑人 | 在线 av 中文字幕|