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

    1954—2016年疏勒河上游徑流變化特征及影響因素

    2022-06-14 16:17:18孫美平姚曉軍彭莉斌牛舒婷
    冰川凍土 2022年2期
    關(guān)鍵詞:疏勒河融水徑流量

    孫美平, 張 磊, 姚曉軍, 彭莉斌, 張 浩, 牛舒婷

    (1.西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,甘肅蘭州730070; 2.中國科學(xué)院西北生態(tài)環(huán)境資源研究院,甘肅蘭州730000)

    0 引言

    河川徑流是重要的地表水資源,尤其是在西北干旱區(qū),缺水與生態(tài)問題相互交織,是制約區(qū)域社會(huì)經(jīng)濟(jì)可持續(xù)發(fā)展的主要瓶頸[1-2]。1987 年以來西北干旱區(qū)呈現(xiàn)明顯暖濕化趨勢,致使冰川快速消融、山區(qū)降水增加,水資源時(shí)空分布和水循環(huán)過程發(fā)生了顯著變化[3-5]。IPCC 第六次評(píng)估報(bào)告將氣候變化對(duì)水文過程和情勢的影響作為第一工作組和第二工作組的焦點(diǎn)領(lǐng)域,并重點(diǎn)關(guān)注氣候變化背景下流域尺度上的水資源變化及其影響[6]。

    疏勒河流域位于祁連山西段,其水資源是維持下游綠洲及城鎮(zhèn)穩(wěn)定和繁榮的重要保障。由于特殊的地理位置,疏勒河流域在我國西部生態(tài)安全和經(jīng)濟(jì)發(fā)展中具有十分重要的戰(zhàn)略地位。在氣候變暖背景下,西北干旱區(qū)冰川普遍退縮[7],其融水補(bǔ)給的內(nèi)流水系徑流變化趨勢以及對(duì)氣候變化的響應(yīng)規(guī)律一直備受關(guān)注。早在20 世紀(jì)90 年代初,楊針娘[8]綜合應(yīng)用冰川融水徑流模數(shù)法、流量與氣溫關(guān)系法、對(duì)比觀測實(shí)驗(yàn)法等對(duì)甘肅河西地區(qū)主要水文站點(diǎn)的徑流變化及冰川融水補(bǔ)給進(jìn)行了較系統(tǒng)分析,估算河西內(nèi)流水系冰川融水徑流總量約9.99×108m3,占河川徑流量的3.3%。目前,關(guān)于疏勒河流域徑流變化特征及規(guī)律已形成一定認(rèn)識(shí)。孫棟元等[9]利用線性傾向方法分析了疏勒河干流不同時(shí)間尺度的徑流量變化和突變特征,發(fā)現(xiàn)流域徑流集中于5—9月,季節(jié)變化明顯。李計(jì)生等[10]應(yīng)用坎德爾秩次相關(guān)法分析了1956—2013 年疏勒河流域出山徑流變化規(guī)律,并預(yù)測了其變化趨勢,認(rèn)為2014—2018 年疏勒河干流徑流偏豐。李浩杰[11]基于多源遙感數(shù)據(jù)分析了疏勒河上游冰川面積和高程變化,發(fā)現(xiàn)2000—2015年期間該區(qū)域冰川呈退縮減薄趨 勢,其 中 冰 川 面 積 減 少(57.6±2.68)km2(11.7%),冰川表面高程降低(2.58±0.6)m,冰儲(chǔ)量減 少(1.271±0.307)km3。Zhang 等[12]應(yīng) 用VICCAS 模型模擬了1971—2012 年期間疏勒河冰川融水對(duì)徑流貢獻(xiàn),研究表明,流域冰川面積覆蓋占比僅約4%,但冰川融水對(duì)河流徑流的平均貢獻(xiàn)約為23.6%。Wu 等[13]應(yīng)用SPHY 模型模擬了疏勒河上游1971—2020年水文過程變化,得到冰川徑流貢獻(xiàn)為28%,隨著冰川面積的減少和溫度的升高,冰川和積雪徑流對(duì)總徑流的年貢獻(xiàn)率呈下降趨勢。李洪源等[14]應(yīng)用寒區(qū)水文模型亦定量模擬了疏勒河流域徑流過程,得出冰川融水徑流占流域總徑流量的30.5%,且1971—2015 年徑流增加近70%。盡管學(xué)者們在不同時(shí)期分析和總結(jié)了疏勒河流域徑流與冰川變化特征,但受限于模型方法、時(shí)空尺度以及可獲取的資料,使得對(duì)該流域徑流變化的系統(tǒng)認(rèn)識(shí)仍有待加強(qiáng)。

    本文基于1954—2016 年昌馬堡水文站日徑流資料、疏勒河流域周邊氣象臺(tái)站和探空站觀測資料以及中國兩次冰川編目數(shù)據(jù),以徑流變化特征及對(duì)氣候變化的響應(yīng)為切入點(diǎn),系統(tǒng)分析疏勒河上游徑流長時(shí)間序列變化趨勢、周期特征及關(guān)鍵氣候要素,以期揭示氣候-冰川變化背景下的冰川融水補(bǔ)給型河流徑流變化規(guī)律,從而為西北干旱區(qū)水資源配置優(yōu)化和科學(xué)管理提供參考依據(jù)。

    1 研究區(qū)概況

    疏勒河是河西內(nèi)流水系中僅次于黑河的第二大河流,發(fā)源于祁連山脈西段托來南山與疏勒南山之間的沙果林那穆吉木嶺,在青海省境內(nèi)自東南向西北流經(jīng)高山峽谷地區(qū),進(jìn)入甘肅省境內(nèi)折向北流入昌馬盆地和河西走廊沖積洪積平原(圖1)。昌馬堡水文站以上為疏勒河上游,集水區(qū)面積為10 946 km2,海拔介于1 990~5 826 m,平均海拔3 885 m,最高峰為團(tuán)結(jié)峰(又名崗則吾結(jié),5 826 m)。流域年平均氣溫為5.3°C,年平均降水量為93.1 mm,屬大陸荒漠干旱型氣候[15]。上游祁連山區(qū)降水較豐沛,現(xiàn)代冰川發(fā)育。根據(jù)中國第二次冰川編目,疏勒河流域是河西內(nèi)流水系擁有冰川資源最多的二級(jí)流域,共有冰川660 條,冰川面積和冰儲(chǔ)量分別為509.87 km2和29.66 km3,各占祁連山冰川相應(yīng)總量的24.59%、31.91%和35.11%[16]。需要說明的是,昌馬堡水文站以上集水區(qū)并不包括祁連山面積最大的老虎溝12號(hào)冰川(又名透明夢柯冰川,編碼為5Y448D0012,面積20.42 km2)所在的老虎溝流域。除冰川外,流域上游下墊面多為寒漠和裸地。

    圖1 疏勒河上游及冰川、水文站與氣象站分布Fig. 1 The upper Shule River and distribution of glaciers,hydrological and meteorological stations

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

    2.1 數(shù)據(jù)來源

    本文所用數(shù)據(jù)包括疏勒河上游昌馬堡水文站及流域周邊氣象臺(tái)站、探空站長時(shí)間觀測資料(表1)。其中,昌馬堡水文站1954—2016年日徑流數(shù)據(jù)為甘肅省水文水資源局的實(shí)測整編資料,用于分析疏勒河上游徑流變化特征;流域周邊氣象臺(tái)站包括托勒、酒泉、玉門鎮(zhèn)、安西、敦煌、冷湖和大柴旦7 個(gè)國家標(biāo)準(zhǔn)氣象站點(diǎn),觀測要素為1954—2016年降水量、平均氣溫、潛在蒸發(fā)、相對(duì)濕度、風(fēng)速等日數(shù)據(jù),用于分析疏勒河上游徑流變化特征及與氣候變化間的關(guān)系;疏勒河流域周邊探空站極少,僅有馬鬃山、敦煌和酒泉3 個(gè)探空站,考慮到2010 年后探空觀測儀器更新,本文選用1971—2010 年夏季(6—8月)700 hPa、500 hPa 和400 hPa 等壓面的氣溫、位勢高度等日數(shù)據(jù),用于計(jì)算0°C 層高度及對(duì)汛期疏勒河上游徑流的影響。上述7個(gè)國家標(biāo)準(zhǔn)氣象站點(diǎn)和3 個(gè)探空站數(shù)據(jù)從國家氣象科學(xué)數(shù)據(jù)中心(http://data.cma.cn)獲得。此外,疏勒河上游冰川數(shù)據(jù)來自中國第一次冰川編目[17]和中國第二次冰川編目[18],前者數(shù)據(jù)源為1956—1963 年的1∶5 萬航測圖并經(jīng)過人工修訂,后者數(shù)據(jù)源為2005—2010 年的Landsat TM 遙感影像。

    表1 疏勒河上游水文站、周邊氣象站和探空站信息Table 1 Information of hydrological station in upper Shule River,meteorological and radiosonde stations

    2.2 研究方法

    2.2.1 線性傾向估計(jì)與滑動(dòng)平均

    為揭示1954—2016 年疏勒河上游徑流變化趨勢及年代際波動(dòng)特征,采用線性傾向估計(jì)和10 a 滑動(dòng)平均方法,計(jì)算徑流隨時(shí)間變化相關(guān)系數(shù)并檢驗(yàn)其顯著性水平,從而判斷徑流變化趨勢是否顯著,線性傾向估計(jì)和滑動(dòng)平均計(jì)算方法見文獻(xiàn)[19]。

    2.2.2 Mann-Kendall趨勢檢驗(yàn)

    Mann-Kendall(M-K)趨勢檢驗(yàn)是世界氣象組織(World Meteorological Organization,WMO)推薦并廣泛應(yīng)用于氣象、水文變化趨勢與突變檢驗(yàn)的一種非參數(shù)統(tǒng)計(jì)方法,其優(yōu)點(diǎn)在于樣本不需服從一定的分布,且不受少數(shù)異常值的干擾,定量化程度高。在M-K 趨勢檢驗(yàn)方法中,其趨勢分析和突變檢測通過統(tǒng)計(jì)量UFk和UBk進(jìn)行判別[20-21]。在給定的顯著性水平α下,當(dāng)M-K 檢驗(yàn)曲線中UFk線在臨界線內(nèi)變動(dòng),表明變化趨勢和突變不明顯。當(dāng)UFk值大于0意味著序列呈上升趨勢,反之為下降趨勢;當(dāng)UFk或UBk值超過臨界值時(shí),表明上升或下降趨勢顯著。若UFk和UBk兩條曲線出現(xiàn)交點(diǎn),且交點(diǎn)在臨界線內(nèi),則該交點(diǎn)對(duì)應(yīng)的時(shí)刻即為突變開始年。

    2.2.3 小波分析與經(jīng)驗(yàn)?zāi)B(tài)分解

    小波分析(wavelet)亦稱多分辨分析,可以將信號(hào)和圖像分解成交織在一起的多尺度成分,得到時(shí)間系列相關(guān)參數(shù)(如平均值、頻率、振幅)的變化特征,從而能夠不斷聚焦到所研究對(duì)象的微小細(xì)節(jié)。相較于小波分析,經(jīng)驗(yàn)?zāi)B(tài)分解法(empirical mode decomposition,EMD)能更直觀、直接地得到有物理意義的頻率,其實(shí)質(zhì)是通過特征時(shí)間尺度識(shí)別數(shù)據(jù)資料中所內(nèi)含的所有振動(dòng)模態(tài)[22]。對(duì)于給定的一時(shí)間序列X(t),EMD 方法通過不斷剔除極大值和極小值連接上下包絡(luò)的均值將原數(shù)據(jù)信號(hào)分解為:

    式中:cj(t)為一個(gè)本征模函數(shù)(intrinsic mode function,IMF)分量;Rn(t)為殘余分量,一般為數(shù)據(jù)序列的平均趨勢。EMD方法的關(guān)鍵是經(jīng)驗(yàn)?zāi)J椒纸猓箯?fù)雜資料分解為有限個(gè)IMF,所分解出來的各IMF分量不斷從高頻到低頻進(jìn)行提取,最終得到一個(gè)頻率接近為0的殘余分量項(xiàng)。小波分析和經(jīng)驗(yàn)?zāi)B(tài)分解兩種方法常用于水文序列的時(shí)頻結(jié)構(gòu)分析,以探求徑流階段性和周期性的變化規(guī)律,二者具體算法見文獻(xiàn)[23-24]。

    2.2.4 分層回歸

    分層回歸(hierarchical multiple regression)又稱層次回歸,其本質(zhì)是對(duì)所建立的兩個(gè)及兩個(gè)以上回歸模型進(jìn)行比較,根據(jù)兩個(gè)模型所解釋的變異量差異判定自變量的重要性,且差異量可用統(tǒng)計(jì)顯著性來估計(jì)和檢驗(yàn)。通常情況下,采用容忍度和方差膨脹因子進(jìn)行自變量之間的多重共線性診斷,診斷標(biāo)準(zhǔn)為容忍度小于0.1 且方差膨脹因子大于10 時(shí),則認(rèn)為存在嚴(yán)重的多重共線性。在分層回歸方法中,決定性系數(shù)(R2)是一個(gè)重要指標(biāo),可反映自變量解釋因變量變異的程度,從其改變量可以檢驗(yàn)增加的自變量是否具有統(tǒng)計(jì)學(xué)意義[25]。

    3 疏勒河上游徑流長時(shí)間序列觀測事實(shí)

    3.1 徑流變化趨勢

    圖2為1954—2016 年昌馬堡水文站日徑流序列和年內(nèi)徑流分布及變化趨勢。從日尺度徑流序列[圖2(a)]來看,疏勒河上游徑流量隨著時(shí)間推移呈增大趨勢,尤其是2000 年后徑流量增加明顯,較2000 年前增加50.5%。1954—2016 期間,昌馬堡水文站年徑流最大值與最小值之比為4.18,年徑流變差系數(shù)僅0.29,表明疏勒河上游徑流量年際變化小,即徑流較為穩(wěn)定,這與流域內(nèi)冰雪融水補(bǔ)給調(diào)蓄有關(guān)。然而,疏勒河上游徑流年內(nèi)分配極不均衡[圖2(b)],汛期(5—9月)來水量占全年徑流總量的57%~79%,月徑流量的極值分別達(dá)到0.09×108m3和5.95×108m3。從年內(nèi)徑流變化趨勢[圖2(c)]看,各月徑流自有觀測以來均呈增加趨勢,且主要出現(xiàn)在冰雪消融期,尤其是7—8 月增加幅度最大,達(dá)0.24×108m3·(10a)-1;而1—4 月和10—12 月徑流年內(nèi)變化微弱,其中1—3 月徑流增加幅度僅0.03×108m3·(10a)-1,僅是7—8月徑流增幅的1/8。

    圖2 昌馬堡水文站日徑流序列及年內(nèi)徑流分布與變化Fig. 2 Daily series,intra-annual distribution and variation of runoff at Changmabao hydrological station

    3.2 徑流變化特征檢測

    1954—2016年昌馬堡水文站實(shí)測年徑流量呈明顯的上升趨勢[圖3(a)],變化速率為1.00×108m3·(10a)-1,且通過0.01 顯著性檢驗(yàn)。從10 a 滑動(dòng)平均曲線來看,1970s—1990s中期疏勒河上游來水偏少,年平均徑流量為9.84×108m3;自21 世紀(jì)初以來進(jìn)入豐水期,并維持在高位波動(dòng),年平均徑流量為12.60×108m3,較2000 年 前 增 加 近28%。Mann-Kendall檢驗(yàn)統(tǒng)計(jì)[圖3(b)]同樣表明,疏勒河上游年徑流量自1960年以來呈上升趨勢,且該趨勢在1998年超過顯著性水平0.05臨界線,1999年徑流發(fā)生突變,這與徑流變化曲線一致。無論是汛期還是非汛期,1954—2016 年昌馬堡水文站觀測到的徑流均呈顯著上升趨勢[圖3(c)、3(e)],變化率分別為0.73×108m3·(10a)-1和0.31×108m3·(10a)-1,其中汛期多年平均徑流量為7.39×108m3,遠(yuǎn)高于非汛期多年平均徑流量(2.44×108m3)。與年徑流突變年份相同,汛期徑流突變亦出現(xiàn)在1999年[圖3(d)],而非汛期徑流量不存在明顯的突變年份[圖3(f)]。

    圖3 昌馬堡水文站徑流變化及Mann-Kendall趨勢檢驗(yàn)Fig.3 Runoff variation and Mann-Kendall trend test at Changmabao hydrological station

    3.3 徑流周期分析

    由小波實(shí)部系數(shù)[圖4(a)]和小波方差[圖4(b)]可知,1954—2016 年期間疏勒河上游年徑流量呈現(xiàn)3 個(gè)時(shí)間尺度的周期變化規(guī)律:10~16 a 尺度上表現(xiàn)為枯、豐交替的準(zhǔn)3 次振蕩,且具有全域性;6~10 a尺度上存在準(zhǔn)4次豐、枯振蕩;2~6 a尺度上存在一定交替振蕩,但較為紊亂。從年徑流小波方差判斷,12~15 a尺度上振蕩能量最強(qiáng),7~9 a尺度上振蕩能量較弱,即15 a為第一主周期,7 a為第二主周期。疏勒河上游汛期徑流量的小波實(shí)部系數(shù)表現(xiàn)出與年徑流量相一致的周期情況,分別為準(zhǔn)16 a、8 a 和5 a;非汛期小波方差在小波影響錐有效情況下存在2個(gè)明顯的振蕩,分別為12~13 a和6~7 a。

    對(duì)1954—2016 年昌馬堡水文站年徑流量時(shí)間序列進(jìn)一步分解,可得到7 個(gè)具有不同波動(dòng)周期的IMF 分量和一個(gè)殘余趨勢項(xiàng)Residual[圖4(c)],表明疏勒河上游年徑流量具有多時(shí)間尺度的波動(dòng)特性,且總體表現(xiàn)為隨著階數(shù)增加,IMF分量振幅減小而周期增長。其中,IMF1、IMF2 為高頻分量,IMF1波動(dòng)頻率最高,保持著3 a左右的波動(dòng)周期,而IMF2分量對(duì)應(yīng)準(zhǔn)6 a 周期,二者在1980s 和1990s 波動(dòng)幅度小且穩(wěn)定,1980 年前和2000 年后波動(dòng)幅度較大;IMF3 和IMF4 為中頻分量,對(duì)應(yīng)對(duì)波動(dòng)周期分別為9~10 a、15~20 a;IMF5~I(xiàn)MF7 為低頻分量,三個(gè)分量對(duì)應(yīng)周期分別為準(zhǔn)40 a、50 a 和60 a。殘差項(xiàng)(Residual)表明1954—2016 年疏勒河上游年徑流量亦呈增加趨勢。

    圖4 昌馬堡水文站年徑流Morlet小波分析及經(jīng)驗(yàn)?zāi)B(tài)分解[(a)為小波實(shí)部圖,小波影響椎(圖中弧線)以內(nèi)的小波方差對(duì)應(yīng)周期為有效周期;(b)為小波方差圖,紅色實(shí)線為方差曲線,灰色虛線為傅立葉譜,黑色虛線為95%置信水平;(c)為年徑流量經(jīng)驗(yàn)?zāi)B(tài)分解]Fig.4 Morlet wavelet analysis and empirical modal decomposition of annual runoff at Changmabao hydrological station(subfigure(a)is the real part of Morlet wavelet and the curve corresponding to wavelet variance in centrum represents the effective period;subfigure(b)is wavelet variance,and the red solid line,gray dashed line and black dashed line denote variance,F(xiàn)ourier spectrum and confidence level with 95%,respectively;subfigure(c)is empirical mode decomposition of annual runoff)

    4 疏勒河上游徑流變化影響因素

    4.1 氣候變化對(duì)徑流變化的影響

    昌馬堡水文站以上人類活動(dòng)較少,其徑流變化可真實(shí)反映疏勒河上游水文自然變化過程。圖5反映了疏勒河上游及周邊7 個(gè)氣象站1956—2016 年年均氣溫、正積溫、降水與潛在蒸發(fā)量變化。顯然,1956—2016 年期間疏勒河上游流域年平均氣溫呈顯著升高趨勢,增溫速率為0.29°C·(10a)-1;與之類似,年正積溫也總體呈現(xiàn)顯著增大趨勢,尤其是1990s 末以來,年正積溫超過3 200 °C。同期,年降水量表現(xiàn)為在波動(dòng)中增加趨勢,變化速率為4.9 mm·(10a)-1,且通過0.001 顯著性檢驗(yàn)。與上述三個(gè)要素變化趨勢不同的是,疏勒河上游潛在蒸發(fā)量在研究時(shí)段內(nèi)整體表現(xiàn)為下降趨勢,這可能與這一時(shí)期的降水量增加和相對(duì)濕度增大(0.22%·(10a)-1,P>0.05)導(dǎo)致飽和水汽壓差減小抑制地表蒸發(fā)有關(guān),也與同時(shí)期的風(fēng)速顯著減?。?0.09 m·s-1·(10a)-1,P<0.001)有直接關(guān)系??傮w上,疏勒河上游流域在過去60 a間經(jīng)歷了“暖干→冷濕→暖濕”3個(gè)明顯的氣候變化過程,分別對(duì)應(yīng)1956—1974年、1975—1999年和2000—2016年。

    圖5 疏勒河上游年平均氣溫、正積溫、年降水與潛在蒸發(fā)量變化Fig. 5 Changes in annual mean temperature,cumulative positive temperature,annual precipitation and potential evapotranspiration in upper Shule River

    為辨析影響疏勒河上游徑流的關(guān)鍵氣候要素,采用分層回歸模型檢驗(yàn)上述氣象要素對(duì)昌馬堡水文站徑流的影響程度及重要性(表2)。若僅考慮年降水影響(模型1),決定性系數(shù)R2為0.442,表明降水單一要素僅能解釋徑流貢獻(xiàn)的44%左右;當(dāng)增加正積溫要素(模型2)后,二者對(duì)年徑流均有顯著影響,且通過0.001顯著性水平檢驗(yàn),決定性系數(shù)R2升至0.807,表明疏勒河上游80%以上徑流可由降水和正積溫來解釋。根據(jù)標(biāo)準(zhǔn)化系數(shù),正積溫對(duì)徑流影響強(qiáng)于降水,這與疏勒河上游流域發(fā)育有面積較大的冰川、積雪和凍土密切相關(guān),因正積溫這一熱量指標(biāo)能較好地反映冰凍圈要素的消融強(qiáng)度,即冰川、積雪融水和凍土凍融補(bǔ)給大于大氣降水影響。當(dāng)增加年平均氣溫(模型3)和年蒸發(fā)(模型4)兩個(gè)要素時(shí),二者均未通過顯著性檢驗(yàn),且兩個(gè)模型中的決定性系數(shù)較模型2 變化微小,這意味著年降水和年正積溫是決定疏勒河上游徑流變化的關(guān)鍵氣候因子,尤其是1990s 末以來二者的同步上升導(dǎo)致昌馬堡水文站徑流快速增大。

    表2 疏勒河上游徑流與各氣象要素的分層回歸模型Table 2 Hierarchical multiple regression of runoff with each meteorological element in upper Shule River

    4.2. ℃層高度對(duì)汛期徑流影響

    冰川融水是疏勒河上游徑流的重要補(bǔ)給源。伴隨氣溫上升,冰川消融區(qū)范圍和強(qiáng)度不斷擴(kuò)大和加劇,其融水對(duì)河流徑流補(bǔ)給持續(xù)加強(qiáng)繼而進(jìn)入汛期。冰川消融與0 ℃層高度變化存在密切關(guān)系,后者又稱冰凍層高度、消融層高度,強(qiáng)烈影響著高海拔地區(qū)的凍融過程,對(duì)冰雪消融有明確的物理意義[26]。根據(jù)馬鬃山、酒泉和敦煌三個(gè)探空站700 hPa、500 hPa 和400 hPa 等壓面的氣溫(T)和位勢高度(H)數(shù)據(jù),二者分別存在以下函數(shù)關(guān)系:T=-0.00702H+32.63(馬鬃山),T=-0.00721H+35.15(酒泉),T=-0.00684H+33.06(敦煌)。由圖6(a)可知,1971—2010 年三個(gè)探空站夏季0 ℃層高度均呈顯著升高趨勢,變化速率分別為77.3 m·(10a)-1、52.7 m·(10a)-1和48.6 m·(10a)-1,這與該流域冰川中值面積海拔呈上升趨勢一致[16]。夏季0 ℃層高度的上升意味著暴露在0 ℃層之下的冰川消融區(qū)范圍增大,冰川物質(zhì)虧損使得更多的冰川冰轉(zhuǎn)換為液態(tài)融水。由昌馬堡水文站夏季徑流與0 ℃層高度散點(diǎn)擬合可知,疏勒河上游夏季徑流量與三個(gè)探空站的0 ℃層高度均存在顯著的正相關(guān)關(guān)系,相關(guān)系數(shù)在0.7 以上,且均通過0.01 顯著性檢驗(yàn)[圖7(b)],這與陳忠升等[27]對(duì)西北干旱區(qū)夏季徑流量對(duì)大氣0 ℃層高度的響應(yīng)研究結(jié)果相吻合,表明西北干旱區(qū)夏季0 ℃層高度升降已成為影響冰川融水補(bǔ)給型河流徑流量變化的關(guān)鍵指示因子。

    圖6 疏勒河流域夏季0°C層高度變化及與徑流量變化關(guān)系Fig.6 The change of 0 ℃layer height and relation between runoff in summer and 0 ℃layer height in upper Shule River basin

    4.3 冰川變化對(duì)疏勒河上游徑流影響

    疏勒河上游流域冰川融水主要源自祁連山脈的疏勒南山,1956—2010 年祁連山冰川自西向東呈加速退縮趨勢[16]。由昌馬堡水文站以上疏勒河流域兩次冰川編目數(shù)據(jù)(表3)可知,1966—2006 年期間流域冰川整體呈退縮態(tài)勢,冰川面積減少79.17 km2(-19.06%)。隨著冰川萎縮,間接導(dǎo)致較大規(guī)模冰川轉(zhuǎn)為小規(guī)模冰川,如面積<0.1 km2的冰川數(shù)量增加63 條,面積<1.0 km2的冰川數(shù)量占比由78.02%上升為81.44%,使得該流域冰川數(shù)量以小規(guī)模冰川為主這一特征更加凸顯,而冰川面積減少則以≥1.0 km2的冰川為主,其面積減少(-65.12 km2)遠(yuǎn)大于<1.0 km2的冰川(-14.05 km2)。相比河西內(nèi)流水系其他流域,疏勒河流域冰川平均規(guī)模最大(0.78 km2)[16],其中面積介于5~10 km2、10~20 km2的冰川分別有13 條和4 條,這也是該流域冰川融水對(duì)徑流貢獻(xiàn)相較其他流域大的原因,已有研究得出疏勒河流域冰川融水約占出山徑流比重的23.6%~30.5%[12-14],表明冰川變化對(duì)于疏勒河上游徑流變化具有舉足輕重的作用。

    表3 昌馬堡水文站以上疏勒河流域不同面積等級(jí)冰川變化Table 3 Changes in glaciers of different area sizes in the Shule River basin above Changmabao hydrological station

    5 結(jié)論

    本文基于昌馬堡水文站1954—2016 年徑流觀測資料及疏勒河流域周邊氣象觀測資料,結(jié)合兩次冰川編目數(shù)據(jù),較系統(tǒng)分析了疏勒河上游長序列徑流變化事實(shí)及其影響因素,得出以下主要結(jié)論:

    (1)1954—2016 年疏勒河上游年徑流量總體呈顯著上升趨勢,變化速率為1.00×108m3·(10a)-1,其中1970s—1990s 初徑流偏枯(平均年徑流量為8.86×108m3),2000 年之后徑流處于豐水期(平均年徑流量為13.19×108m3),徑流突變發(fā)生在1999 年;汛期和非汛期徑流亦均呈顯著增加趨勢,汛期徑流突變時(shí)間與年徑流一致,而非汛期徑流不存在明顯的突變年份。疏勒河上游徑流周期振蕩特征顯著,年徑流、汛期和非汛期徑流大體都存在15a 和7a 左右的顯著主周期。

    (2)正積溫和降水是決定疏勒河上游徑流變化的主導(dǎo)氣候因子,二者可解釋80%以上的徑流變化,且正積溫對(duì)徑流的影響強(qiáng)于降水,這與疏勒河上游流域冰川、積雪和凍土普遍發(fā)育有關(guān)。疏勒河上游夏季徑流與0 °C 層高度存在顯著正相關(guān)關(guān)系(r≥0.7),表明夏季0 ℃層高度升降是影響疏勒河上游徑流量變化的關(guān)鍵指示因子,可間接表征疏勒河上游汛期徑流量變化。1966—2006 年疏勒河上游流域冰川冰儲(chǔ)量減少約5.77 km3,冰川變化對(duì)于疏勒河上游徑流變化具有舉足輕重的作用。

    猜你喜歡
    疏勒河融水徑流量
    我在融水唱貝江
    融水糯米柚種植栽培管理與技術(shù)探究
    疏勒河灌區(qū)防洪問題初探
    疏勒河流域水權(quán)交易探索與對(duì)策
    疏勒河流域高效灌溉系統(tǒng)集成技術(shù)
    例舉高中數(shù)學(xué)解題切入點(diǎn)的找尋
    水文比擬法在計(jì)算河川徑流量時(shí)的修正
    疏勒河
    飛天(2016年7期)2016-07-13 07:09:15
    廣西苗鄉(xiāng)“蘆笙鬥馬節(jié)”驚險(xiǎn)刺激
    澳門月刊(2014年11期)2014-11-15 15:53:07
    SCS模型在紅壤土坡地降雨徑流量估算中的應(yīng)用
    非洲黑人性xxxx精品又粗又长| 午夜免费鲁丝| 亚洲色图av天堂| 久久久久久九九精品二区国产 | 久久久久久国产a免费观看| 黄片大片在线免费观看| 免费看日本二区| 久久国产亚洲av麻豆专区| 在线观看66精品国产| 亚洲国产精品sss在线观看| АⅤ资源中文在线天堂| 在线观看免费午夜福利视频| 国产片内射在线| 日本熟妇午夜| 中文字幕人成人乱码亚洲影| 精品国产乱码久久久久久男人| 两性夫妻黄色片| 在线播放国产精品三级| 黑人巨大精品欧美一区二区mp4| 国产不卡一卡二| 免费观看人在逋| 国产精品久久视频播放| 精品久久久久久久久久久久久 | 90打野战视频偷拍视频| 无人区码免费观看不卡| 91九色精品人成在线观看| 伊人久久大香线蕉亚洲五| 久久天堂一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 久久人人精品亚洲av| 国产精品日韩av在线免费观看| 最近在线观看免费完整版| 最近最新中文字幕大全电影3 | 午夜亚洲福利在线播放| 老司机午夜福利在线观看视频| 中文资源天堂在线| 最好的美女福利视频网| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 国产1区2区3区精品| 老汉色∧v一级毛片| 免费在线观看完整版高清| 午夜福利在线观看吧| 国产激情久久老熟女| 免费在线观看视频国产中文字幕亚洲| 亚洲精品在线观看二区| 91av网站免费观看| 成人三级做爰电影| 欧美激情久久久久久爽电影| 91成人精品电影| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 91九色精品人成在线观看| 韩国精品一区二区三区| 性色av乱码一区二区三区2| 欧美黄色淫秽网站| 亚洲av五月六月丁香网| 久久精品aⅴ一区二区三区四区| 日韩av在线大香蕉| 亚洲欧洲精品一区二区精品久久久| 欧美精品亚洲一区二区| 国产91精品成人一区二区三区| 18美女黄网站色大片免费观看| 法律面前人人平等表现在哪些方面| 又黄又爽又免费观看的视频| 亚洲激情在线av| 这个男人来自地球电影免费观看| 日日爽夜夜爽网站| 97碰自拍视频| 久久久久久久精品吃奶| 一本精品99久久精品77| 国产99白浆流出| 中文字幕人妻丝袜一区二区| 午夜激情福利司机影院| 国产伦一二天堂av在线观看| 国产又色又爽无遮挡免费看| 两性夫妻黄色片| 中文字幕久久专区| 亚洲全国av大片| 女性生殖器流出的白浆| e午夜精品久久久久久久| av电影中文网址| 无遮挡黄片免费观看| 性色av乱码一区二区三区2| 色综合欧美亚洲国产小说| 久久久久免费精品人妻一区二区 | 人人妻人人澡欧美一区二区| 亚洲国产毛片av蜜桃av| 啦啦啦 在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 日本精品一区二区三区蜜桃| 欧美绝顶高潮抽搐喷水| 搡老熟女国产l中国老女人| 日韩欧美在线二视频| www国产在线视频色| 久久九九热精品免费| 国产av又大| e午夜精品久久久久久久| 国产精品 欧美亚洲| 精品久久久久久久久久免费视频| 波多野结衣巨乳人妻| 欧美不卡视频在线免费观看 | 一区二区三区国产精品乱码| 女同久久另类99精品国产91| 亚洲精品一区av在线观看| 日韩欧美免费精品| 成人国语在线视频| 欧美精品亚洲一区二区| 变态另类丝袜制服| 久久精品国产99精品国产亚洲性色| 两性午夜刺激爽爽歪歪视频在线观看 | 久久婷婷成人综合色麻豆| 男女床上黄色一级片免费看| tocl精华| 成人一区二区视频在线观看| 少妇粗大呻吟视频| 少妇裸体淫交视频免费看高清 | 亚洲欧洲精品一区二区精品久久久| 欧美三级亚洲精品| a级毛片a级免费在线| 日韩一卡2卡3卡4卡2021年| 黄色毛片三级朝国网站| 18禁黄网站禁片免费观看直播| 亚洲国产精品成人综合色| 国产片内射在线| 黄色 视频免费看| 国产黄a三级三级三级人| 国产成人av教育| 亚洲一区二区三区不卡视频| 在线永久观看黄色视频| 久99久视频精品免费| 9191精品国产免费久久| 男女下面进入的视频免费午夜 | www日本在线高清视频| 亚洲精品久久成人aⅴ小说| 黄片大片在线免费观看| 国产亚洲精品综合一区在线观看 | 狠狠狠狠99中文字幕| 午夜久久久在线观看| 欧美成人午夜精品| 亚洲无线在线观看| 亚洲av电影不卡..在线观看| 亚洲中文字幕一区二区三区有码在线看 | 精品不卡国产一区二区三区| videosex国产| 国产精品久久电影中文字幕| 黄网站色视频无遮挡免费观看| 亚洲熟妇熟女久久| 久久精品国产99精品国产亚洲性色| 一级a爱视频在线免费观看| 99热只有精品国产| 99热6这里只有精品| 亚洲自偷自拍图片 自拍| 男女做爰动态图高潮gif福利片| 真人一进一出gif抽搐免费| 男男h啪啪无遮挡| 国产成人系列免费观看| 在线观看午夜福利视频| 人人妻人人看人人澡| 亚洲国产毛片av蜜桃av| 1024手机看黄色片| 国产成+人综合+亚洲专区| 亚洲精品久久成人aⅴ小说| 精品欧美国产一区二区三| 亚洲av片天天在线观看| 免费av毛片视频| 欧美人与性动交α欧美精品济南到| 亚洲午夜理论影院| 亚洲激情在线av| 国产成人av教育| 老熟妇乱子伦视频在线观看| 亚洲国产欧洲综合997久久, | 国产免费av片在线观看野外av| 亚洲色图 男人天堂 中文字幕| 亚洲国产欧美网| 一a级毛片在线观看| 大香蕉久久成人网| 亚洲最大成人中文| 在线免费观看的www视频| 亚洲人成伊人成综合网2020| 一边摸一边做爽爽视频免费| 叶爱在线成人免费视频播放| 黄色a级毛片大全视频| 午夜福利在线在线| 亚洲成av人片免费观看| 亚洲最大成人中文| 色综合站精品国产| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| 欧美色视频一区免费| 国产熟女xx| 国产蜜桃级精品一区二区三区| av福利片在线| 一级a爱视频在线免费观看| 美女高潮到喷水免费观看| 久久久久国产一级毛片高清牌| 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 亚洲精品av麻豆狂野| 村上凉子中文字幕在线| 国产av又大| 日韩国内少妇激情av| 在线免费观看的www视频| 亚洲免费av在线视频| 亚洲人成网站在线播放欧美日韩| 黑丝袜美女国产一区| 法律面前人人平等表现在哪些方面| 久久久久久九九精品二区国产 | 国产精品 欧美亚洲| 天天躁夜夜躁狠狠躁躁| 欧美黑人巨大hd| 成人一区二区视频在线观看| 欧美午夜高清在线| 国产精品亚洲av一区麻豆| 亚洲精品在线美女| 狠狠狠狠99中文字幕| 亚洲天堂国产精品一区在线| 不卡av一区二区三区| 婷婷丁香在线五月| 国产aⅴ精品一区二区三区波| 欧美日韩瑟瑟在线播放| 午夜两性在线视频| 亚洲午夜精品一区,二区,三区| 免费在线观看日本一区| 精品福利观看| 男人舔女人下体高潮全视频| 亚洲人成电影免费在线| 亚洲国产高清在线一区二区三 | aaaaa片日本免费| 一级黄色大片毛片| 精品免费久久久久久久清纯| 男男h啪啪无遮挡| 免费av毛片视频| 亚洲人成电影免费在线| 久久国产精品影院| av视频在线观看入口| 欧美最黄视频在线播放免费| 亚洲成人国产一区在线观看| 成熟少妇高潮喷水视频| 女人高潮潮喷娇喘18禁视频| 免费在线观看黄色视频的| 午夜激情av网站| 欧美性猛交黑人性爽| 亚洲第一电影网av| 村上凉子中文字幕在线| 无人区码免费观看不卡| 久久婷婷成人综合色麻豆| 69av精品久久久久久| 亚洲av中文字字幕乱码综合 | 波多野结衣巨乳人妻| 精品国产美女av久久久久小说| 在线看三级毛片| 一区二区三区激情视频| 中国美女看黄片| 又黄又粗又硬又大视频| 国产精品爽爽va在线观看网站 | 国产又色又爽无遮挡免费看| 十分钟在线观看高清视频www| 不卡一级毛片| 视频在线观看一区二区三区| 岛国在线观看网站| 性色av乱码一区二区三区2| 黄色毛片三级朝国网站| www.999成人在线观看| 成人亚洲精品av一区二区| 神马国产精品三级电影在线观看 | 久久久久久久久免费视频了| 美女午夜性视频免费| x7x7x7水蜜桃| 成人国产一区最新在线观看| 久9热在线精品视频| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 欧洲精品卡2卡3卡4卡5卡区| 成人三级黄色视频| 禁无遮挡网站| 国产亚洲av高清不卡| 亚洲精品中文字幕一二三四区| 成人一区二区视频在线观看| 亚洲国产欧洲综合997久久, | 国产亚洲欧美在线一区二区| 国内精品久久久久精免费| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 麻豆久久精品国产亚洲av| 视频区欧美日本亚洲| 91大片在线观看| 日韩精品免费视频一区二区三区| 免费看日本二区| 国产精品久久久久久亚洲av鲁大| 亚洲成人久久爱视频| 女同久久另类99精品国产91| 18美女黄网站色大片免费观看| 欧美 亚洲 国产 日韩一| 亚洲精品国产精品久久久不卡| 久久国产精品人妻蜜桃| 国产精品免费一区二区三区在线| 亚洲第一电影网av| 少妇被粗大的猛进出69影院| 午夜久久久在线观看| 欧美av亚洲av综合av国产av| 国产区一区二久久| 亚洲真实伦在线观看| 久久香蕉激情| 悠悠久久av| 一本久久中文字幕| 亚洲成人免费电影在线观看| 视频在线观看一区二区三区| 精品不卡国产一区二区三区| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 亚洲aⅴ乱码一区二区在线播放 | 男女做爰动态图高潮gif福利片| 免费看日本二区| av中文乱码字幕在线| 一进一出抽搐动态| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久人妻精品电影| 亚洲午夜理论影院| 黄色视频,在线免费观看| 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 中文在线观看免费www的网站 | 麻豆av在线久日| 又黄又爽又免费观看的视频| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 久久中文字幕一级| 欧美绝顶高潮抽搐喷水| 国产高清有码在线观看视频 | 少妇熟女aⅴ在线视频| 久久中文字幕一级| 老司机深夜福利视频在线观看| 十八禁网站免费在线| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| 99热6这里只有精品| 亚洲专区中文字幕在线| 可以在线观看的亚洲视频| 999久久久精品免费观看国产| 午夜免费激情av| 亚洲欧美精品综合一区二区三区| a级毛片在线看网站| 成年版毛片免费区| 麻豆av在线久日| 黄色成人免费大全| 国产成人影院久久av| 777久久人妻少妇嫩草av网站| 69av精品久久久久久| 欧美精品啪啪一区二区三区| 色播亚洲综合网| 欧美不卡视频在线免费观看 | 法律面前人人平等表现在哪些方面| 精品国产一区二区三区四区第35| 国产成人啪精品午夜网站| 99re在线观看精品视频| 亚洲avbb在线观看| 欧美黑人精品巨大| 级片在线观看| 欧美黑人欧美精品刺激| 日日夜夜操网爽| 88av欧美| xxx96com| 国产不卡一卡二| av福利片在线| 久久久久亚洲av毛片大全| 欧美日韩福利视频一区二区| 天天躁夜夜躁狠狠躁躁| 日韩一卡2卡3卡4卡2021年| 午夜福利在线在线| 三级毛片av免费| 日本免费一区二区三区高清不卡| 18禁黄网站禁片午夜丰满| 999久久久国产精品视频| 久久香蕉激情| 人人妻,人人澡人人爽秒播| 午夜老司机福利片| 午夜影院日韩av| 性欧美人与动物交配| 美女大奶头视频| 国产精品亚洲av一区麻豆| 欧美av亚洲av综合av国产av| 免费在线观看亚洲国产| 动漫黄色视频在线观看| 久热这里只有精品99| 天天躁夜夜躁狠狠躁躁| 嫩草影院精品99| 熟女少妇亚洲综合色aaa.| 波多野结衣av一区二区av| 1024视频免费在线观看| 精品国产乱码久久久久久男人| 99久久国产精品久久久| 超碰成人久久| 国产不卡一卡二| 欧美成人性av电影在线观看| 级片在线观看| 欧美亚洲日本最大视频资源| 波多野结衣高清作品| 亚洲精品国产一区二区精华液| 亚洲全国av大片| a在线观看视频网站| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区久久 | 中文字幕av电影在线播放| 国产av一区二区精品久久| 97超级碰碰碰精品色视频在线观看| 夜夜看夜夜爽夜夜摸| 大型av网站在线播放| 欧美日本视频| АⅤ资源中文在线天堂| 亚洲国产欧美网| 免费观看人在逋| 亚洲一区中文字幕在线| 一边摸一边抽搐一进一小说| 身体一侧抽搐| 热99re8久久精品国产| 国产三级在线视频| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看| 精品国产乱子伦一区二区三区| АⅤ资源中文在线天堂| 一级毛片高清免费大全| 欧美日韩乱码在线| 欧美日韩精品网址| 亚洲五月婷婷丁香| 欧美一区二区精品小视频在线| 99热6这里只有精品| 久久国产乱子伦精品免费另类| 午夜福利成人在线免费观看| 国产一卡二卡三卡精品| 亚洲 欧美一区二区三区| 国产一区二区三区在线臀色熟女| 国产精品乱码一区二三区的特点| 久久精品成人免费网站| 每晚都被弄得嗷嗷叫到高潮| 非洲黑人性xxxx精品又粗又长| 久久久国产欧美日韩av| 丰满人妻熟妇乱又伦精品不卡| 美女高潮喷水抽搐中文字幕| 欧美另类亚洲清纯唯美| 中文字幕人成人乱码亚洲影| 九色国产91popny在线| 欧美激情极品国产一区二区三区| 一本大道久久a久久精品| 久久久久久久精品吃奶| 妹子高潮喷水视频| 99久久精品国产亚洲精品| 欧美在线黄色| 国产在线精品亚洲第一网站| 一个人免费在线观看的高清视频| 亚洲一区中文字幕在线| 亚洲成av人片免费观看| 此物有八面人人有两片| 特大巨黑吊av在线直播 | 亚洲 欧美一区二区三区| 听说在线观看完整版免费高清| 亚洲一区中文字幕在线| 亚洲成av人片免费观看| 免费看a级黄色片| 99热这里只有精品一区 | 成人特级黄色片久久久久久久| 在线观看免费午夜福利视频| 国产精品久久久av美女十八| 久久国产精品人妻蜜桃| 露出奶头的视频| 欧美激情久久久久久爽电影| 欧美成狂野欧美在线观看| 97碰自拍视频| 国产精品久久电影中文字幕| 哪里可以看免费的av片| 一本综合久久免费| 99riav亚洲国产免费| 精品电影一区二区在线| 免费看日本二区| 亚洲激情在线av| 国产1区2区3区精品| 精品免费久久久久久久清纯| 亚洲国产高清在线一区二区三 | 久久久久久免费高清国产稀缺| 日本一区二区免费在线视频| 国产不卡一卡二| 黄色a级毛片大全视频| 听说在线观看完整版免费高清| 亚洲色图av天堂| 麻豆国产av国片精品| 午夜福利18| 欧美中文综合在线视频| 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 免费看十八禁软件| 国产精品久久久久久亚洲av鲁大| 国产久久久一区二区三区| 国产欧美日韩精品亚洲av| 久久久久免费精品人妻一区二区 | 日本熟妇午夜| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 黑人欧美特级aaaaaa片| 亚洲性夜色夜夜综合| 欧美不卡视频在线免费观看 | 丁香六月欧美| 最近最新中文字幕大全免费视频| 日本黄色视频三级网站网址| e午夜精品久久久久久久| 色综合亚洲欧美另类图片| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av高清一级| 精品国产超薄肉色丝袜足j| av片东京热男人的天堂| 1024香蕉在线观看| 欧美乱码精品一区二区三区| 视频在线观看一区二区三区| 午夜福利在线在线| 国产人伦9x9x在线观看| 亚洲性夜色夜夜综合| 一级黄色大片毛片| 免费在线观看亚洲国产| 欧美午夜高清在线| 亚洲国产日韩欧美精品在线观看 | 日本免费一区二区三区高清不卡| 校园春色视频在线观看| 久久香蕉国产精品| 777久久人妻少妇嫩草av网站| 久99久视频精品免费| 一a级毛片在线观看| 欧美最黄视频在线播放免费| 99国产精品一区二区蜜桃av| 久久久久久九九精品二区国产 | 亚洲欧美精品综合一区二区三区| 99久久国产精品久久久| 国产又爽黄色视频| 啦啦啦韩国在线观看视频| 日本免费a在线| 亚洲熟女毛片儿| aaaaa片日本免费| 亚洲激情在线av| 女人爽到高潮嗷嗷叫在线视频| 99riav亚洲国产免费| 国产野战对白在线观看| 两个人看的免费小视频| 久久中文字幕一级| 欧美日韩黄片免| 婷婷丁香在线五月| 国产精品九九99| 手机成人av网站| 日本一区二区免费在线视频| 午夜福利一区二区在线看| 久久 成人 亚洲| 一个人观看的视频www高清免费观看 | 三级毛片av免费| 丰满的人妻完整版| netflix在线观看网站| 高清毛片免费观看视频网站| 国产男靠女视频免费网站| 最新美女视频免费是黄的| 精品国产乱码久久久久久男人| 久久久久久免费高清国产稀缺| 老司机午夜福利在线观看视频| 亚洲av成人av| 老司机在亚洲福利影院| 很黄的视频免费| 国产精品免费一区二区三区在线| 一本大道久久a久久精品| tocl精华| 亚洲精品久久成人aⅴ小说| 色综合婷婷激情| 高清毛片免费观看视频网站| 欧美av亚洲av综合av国产av| 亚洲 欧美一区二区三区| 母亲3免费完整高清在线观看| 18禁观看日本| 午夜久久久久精精品| 午夜亚洲福利在线播放| 一区二区三区高清视频在线| 色尼玛亚洲综合影院| 久久中文字幕人妻熟女| 亚洲一区二区三区色噜噜| 亚洲avbb在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲国产高清在线一区二区三 | 韩国精品一区二区三区| 国产又黄又爽又无遮挡在线| 国产av不卡久久| 精品国产国语对白av| 成人免费观看视频高清| 香蕉国产在线看| 欧美日韩黄片免| 免费在线观看黄色视频的| 午夜免费成人在线视频| 成人国语在线视频| 日韩欧美一区视频在线观看| 国产精品一区二区精品视频观看| 国产成人精品无人区| 亚洲va日本ⅴa欧美va伊人久久| 国产伦人伦偷精品视频| 亚洲av五月六月丁香网| 久久精品91无色码中文字幕| 一区二区三区激情视频| 欧美色视频一区免费| 神马国产精品三级电影在线观看 | 三级毛片av免费| 成人三级黄色视频| 丰满人妻熟妇乱又伦精品不卡| 日本精品一区二区三区蜜桃| 中文字幕久久专区| 精品国产乱码久久久久久男人| 国产精品免费视频内射| 青草久久国产| 日韩欧美国产一区二区入口| 久久草成人影院| 国产一卡二卡三卡精品| 亚洲天堂国产精品一区在线| 又黄又爽又免费观看的视频|