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

    氣候變化和人類活動對瑪納斯河流域徑流及干旱的影響

    2022-07-31 11:04:44肖森元仇澤豪羅雅潔黃晶晶
    人民珠江 2022年7期
    關(guān)鍵詞:人類影響活動

    肖森元,蘇 軍,楊 廣*,劉 兵,張 迪,仇澤豪,羅雅潔,黃晶晶

    (1.石河子大學(xué)水利建筑工程學(xué)院/寒旱區(qū)生態(tài)水利工程兵團重點實驗室,新疆 石河子 832000;2.第八師石河子市水利工程管理服務(wù)中心,新疆 石河子 832000)

    氣候變化和人類活動改變了陸地水文循環(huán)的過程,影響了水文和水資源系統(tǒng)的結(jié)構(gòu)和功能,并對水資源的開發(fā)和利用構(gòu)成重大挑戰(zhàn)[1-2]。近幾十年來,氣溫變暖加速了全球水文循環(huán),從而影響了水循環(huán)過程并改變了盆地降水量與水循環(huán)之間的時空格局[3-4]。人類活動則進一步促進了降水和徑流的變化,導(dǎo)致流域水資源的時空分布發(fā)生了變化[5-6],水資源供需平衡遭到破壞,干旱問題日益突出[7]。因此揭示氣候變化和人類活動影響下流域徑流和干旱演變特征成為當(dāng)前亟待解決的重要問題[8]。

    目前,國內(nèi)外學(xué)者在氣候變化和人類活動對徑流和干旱的影響方面做了大量研究,并取得豐碩的成果。目前應(yīng)用較為廣泛的主要有數(shù)理統(tǒng)計法、彈性系數(shù)法及水文模型法。數(shù)理統(tǒng)計法以徑流和氣候因子間的線性關(guān)系為基礎(chǔ),分離氣候變化和人類活動對徑流影響。劉可群等[9]對洪湖春季干旱分析,指出洪湖春旱增多加重主要受氣候變化和人類活動影響,而人類活動產(chǎn)生的影響更顯著。Alifujiang Y[10]等采用雙累積曲線對伊塞克湖盆地徑流量變化進行等量分析,發(fā)現(xiàn)氣候變化和人類活動對徑流減少的貢獻率分別為64.7%和35.3%,氣候變化是導(dǎo)致伊塞克湖盆地徑流減少的主要原因。彈性系數(shù)法主要基于水熱耦合平衡方程,其中基于Budyko假設(shè)的分析方法應(yīng)用最為廣泛。薛聯(lián)青等[11]基于Budyko 假設(shè),對塔里木河流域徑流變化對氣候變化和人類活動的響應(yīng)進行了敏感性分析。Zheng等[12]基于Budyko假設(shè)的互補關(guān)系方法探討氣候變化和人類活動對黃河流域徑流變化的貢獻率,結(jié)果表明黃河流域徑流減少主要是由于氣候變化,其貢獻為60%~70%。水文模型是分析流域水文過程的工具,可實現(xiàn)不同氣候情景下的徑流模擬。 Seguis等[13]應(yīng)用分布式水文模型對尼日爾的薩赫勒小流域徑流進行模擬研究發(fā)現(xiàn),人類活動對徑流的影響比氣候因素對徑流的影響重要。Lehner等[14]通過采用WaterGAP模型模擬未來水流不同流態(tài)認為,未來歐洲北部到東北部在全球變化的影響下易發(fā)生洪水,而歐洲南部和東南部則易發(fā)生干旱。

    瑪納斯河流域是典型內(nèi)陸河流域,地處中國西北干旱區(qū),水文生態(tài)平衡更加脆弱敏感。近年來,隨著水資源的日益消耗,水問題變得非常尖銳,出現(xiàn)了嚴重的生態(tài)環(huán)境問題,如地下水位下降、農(nóng)田次生鹽漬化、植被大面積退化等。目前,在瑪納斯河流域展開的氣候變化和人類活動研究主要集中在對水文水資源產(chǎn)生的影響進行定量區(qū)分方面[15],并未涉及水文干旱演變,而如何加強人類活動對瑪納斯河流域水文變化的正向引導(dǎo)作用,實現(xiàn)區(qū)域水資源可持續(xù)發(fā)展,亟待對流域水文及干旱演變規(guī)律進行研究?;诖?,本研究通過Mann-Kendall非參數(shù)檢驗法、雙積累曲線法和SRI指數(shù)法對氣候變化和人類活動影響下瑪納斯河流域徑流和水文干旱進行分析,并對比分析了氣候變化和人類活動影響下流域干旱狀況,為流域水土資源合理利用提供基礎(chǔ)依據(jù)。

    1 數(shù)據(jù)來源與研究方法

    1.1 研究區(qū)域與數(shù)據(jù)來源

    瑪納斯河流域位于新疆天山北麓、準噶爾盆地南緣,地理位置為43°27′~45°21′N,85°01′~86°32′E(圖1)。流域地勢東南高西北低,自東向西分布有塔西河、瑪納斯河、寧家河、金溝河以及巴音溝河,其中瑪納斯河是流域內(nèi)流量最大的河流,年均徑流量為12.9×108m3,流域年均降雨量為115~200 mm,屬于冰雪融水和降雨混合補給型河流。由于流域地處新疆腹地,形成氣溫差異大,蒸發(fā)量大的溫帶荒漠氣候,流域年均氣溫為4.7~5.7℃,年均蒸發(fā)量為1 500~2 100 mm。

    本研究數(shù)據(jù)來源于瑪納斯河流域肯斯瓦特水文站與石河子和烏蘭烏蘇氣象站提供的1956—2020年實測的逐月徑流和逐月降雨數(shù)據(jù)。

    圖1 研究區(qū)位置

    1.2 研究方法

    a)徑流變化趨勢及變異分析。Mann-Kendall是基于秩的非參數(shù)檢驗,它對異常值不敏感,并且不要求數(shù)據(jù)呈正態(tài)分布,因此廣泛應(yīng)用于降水、溫度及徑流等氣象水文序列的變化趨勢檢驗[16]。在分析時間序列的趨勢和突變點時,首先計算樣本統(tǒng)計量UF與UB,確定時間序列的變化趨勢。當(dāng)UF超過顯著性置信度臨界線(±1.96)時,該序列在統(tǒng)計上存在顯著的增加/減少趨勢,反之則不表現(xiàn)出顯著的趨勢。若UF和UB 2條曲線的交點在臨界線內(nèi),則交點對應(yīng)的時刻被視為時間序列的突變點。具體步驟見文獻[17-18]。

    b)徑流變化歸因分析。雙累積曲線法是研究2個序列間變化關(guān)系和判別其是否具有一致性的方法,由于該方法簡單、直觀、有效而在水文氣象要素分析中得到廣泛應(yīng)用。首先,將2個相同長度的序列進行累加,并在直角坐標系中繪制關(guān)系線,依據(jù)關(guān)系線的斜率是否發(fā)生變化出現(xiàn)拐點,判斷其是否發(fā)生異常,而出現(xiàn)的拐點對應(yīng)的時間則為發(fā)生異常的年份。以異常年份為時間點,對長度序列進行比較,若相對于變化前,變化后序列出現(xiàn)偏離,則表示變化后序列可能受外界影響出現(xiàn)顯著性變化,否則表示兩序列一致性較好,未出現(xiàn)顯著性變化。而后,將第一個突變點之前有實測資料的時段作為自然期,之后作為變化期,根據(jù)自然期內(nèi)降水-徑流累積關(guān)系,反推變化期受氣候變化影響徑流量,實現(xiàn)氣候變化和人類活動對徑流影響的分離,從而探究人類活動和氣候變化對突變期相對于基準期徑流量變化的貢獻率,對徑流量變化進行歸因分析。具體步驟見文獻[19-21]。

    c)水文干旱分析。標準化徑流指數(shù)(SRI)是持續(xù)時間內(nèi)累積的水文徑流百分比的單位標準正態(tài)偏差,流域內(nèi)不同時間尺度下的徑流情況和干旱程度可以通過SRI的大小反映[22-23]。其計算方法如下。

    首先獲取徑流的時間序列,并對時間序列所代表的樣本進行概率分布擬合:

    (1)

    然后對其進行累加,得到累積概率計算公式見式(2):

    (2)

    最后將累積概率分布函數(shù)轉(zhuǎn)換為標準正態(tài)分布函數(shù):

    (3)

    由于標準化徑流指數(shù)(SRI)的計算方法與標準化降水指數(shù)(SPI)類似,故可參照國家氣候中心頒布的氣象干旱等級分類標準[24]確定SRI指數(shù)的干旱等級。具體等級劃分見表1。

    表1 標準化徑流指數(shù)干旱等級劃分

    2 結(jié)果與分析

    2.1 徑流變化趨勢和變異分析

    根據(jù)1956—2020年實測徑流和降水資料對瑪納斯河流域年徑流變化進行分析可知,從整體上看,研究時段內(nèi)瑪納斯河流域年徑流深呈不規(guī)則周期波動(圖2a),其中1999年徑流深為歷年最大值,達到了435.46 mm,而1992年僅有201.78 mm,為歷年最小值。年降水量在1955—1965 年波動幅度較大且有先增加后減少的趨勢,1966—1977 年呈先波動增加后減少,在1977 年達到最小值,從1978 年起降水量呈波動增加的變化趨勢。從年際上看,65 a來瑪納斯河流域年徑流深呈現(xiàn)出比較明顯的線性增加趨勢,尤其是1995年以后遞增幅度最大。從年代上來看(圖2b),流域年徑流深呈現(xiàn)出“先增后減再增”的變化趨勢。除20世紀60年代到70年代徑流深呈現(xiàn)出減少趨勢外,其他年代徑流深均呈現(xiàn)出不同程度遞增狀態(tài),尤其是20世紀80年代到90年代增加最為明顯,變化幅度高達17.27%。進入21世紀后,流域年徑流變化趨勢有所減緩,10 a變化幅度約為2%。在這一時期,降雨量則呈現(xiàn)“先減后增再平穩(wěn)”的變化趨勢。

    a)流域徑流

    b)降雨量

    利用Mann-Kendall非參數(shù)統(tǒng)計檢驗法對瑪納斯河流域年徑流分析可知(圖3),在0.05顯著性水平下,流域年徑流在1956—2004年均未出現(xiàn)顯著性變化,在2005—2020年統(tǒng)計檢驗指標超過臨界線,表明在此時間段內(nèi)年徑流呈顯著性增加趨勢。UF和UB 2條曲線相交的時間節(jié)點是1995年,且交點位于顯著性置信度臨界線內(nèi),表明流域年徑流發(fā)生突變的時間是1995年。1995年恰是瑪納斯河流域推行節(jié)水灌溉技術(shù)的時期,流域規(guī)模化節(jié)水措施作為重要的人類活動,在迅速推動綠洲化進程的同時,加劇了流域徑流的年際變化趨勢。

    圖3 年徑流量Mann-Kendall檢驗

    2.2 徑流變化影響因素分析

    2.2.1氣候變化和人類活動對徑流的影響

    氣候變化和人類活動是影響地表徑流的2個重要因素,為識別不同驅(qū)動因素的貢獻,本研究利用降水-徑流雙累積曲線來分析量化氣候變化和人類活動對瑪納斯河流域徑流的貢獻比例(圖4)。

    圖4 瑪納斯河流域年降水-徑流雙累積曲線

    由圖4可知,降水-徑流雙累積曲線在1995年出現(xiàn)拐點,說明1995年為其確定的突變年份,這與利用Mann-Kendall非參數(shù)統(tǒng)計檢驗結(jié)論相一致。因此,依據(jù)判定的突變年份,以1995年為時間節(jié)點,將年徑流序列劃分為1956—1995年和1996—2014年2個子序列。選取1956—1995年作為自然期,認為在該時期徑流受自然環(huán)境變化影響較大,受到人類活動影響較小,選取1996—2020年作為變化期,認為在該時期徑流受氣候變化和人類活動的雙重影響。經(jīng)計算基準期和變化期2個子序列的年徑流均值可知,基準期流域年徑流均值為250.93 mm,變化期流域年徑流均值較基準期增加62.43 mm,約增加24.88%,可見,流域徑流在變化期增加較為明顯。

    2.2.2氣候變化和人類活動對徑流的貢獻

    為定量分析氣候變化和人類活動對流域徑流的影響程度,以自然期月累積降水-徑流為原始數(shù)據(jù),構(gòu)建月降水-徑流雙累積關(guān)系模型(∑R=0.77633·∑p-264.62,R2=0.99),反推變化期受氣候變化影響徑流量,實現(xiàn)氣候變化和人類活動對徑流影響的分離,而后通過計算實測值和擬合值差值占變化值的百分比,確定其貢獻率,結(jié)果見表2。

    由表2可知,氣候變化對月徑流量變化的影響量為2.95 mm,貢獻率為56.73%,人類活動對月徑流量變化的影響量為2.25 mm,貢獻率為43.37%,表明氣候變化是導(dǎo)致瑪納斯河流域徑流量變化的主要因素。

    表2 氣候變化和人類活動對徑流變化的貢獻率

    2.3 水文干旱演變特征分析

    在進行水文干旱分析時,采用與SPI指數(shù)相類似的方法,所以SRI指數(shù)和SPI指數(shù)一樣具有多時間尺度特征,可進行不同時間尺度下的干旱特征分析。考慮到研究區(qū)主要發(fā)展綠洲農(nóng)業(yè),干旱事件的發(fā)生又多以季節(jié)性干旱為主,因此,在上述徑流分析的基礎(chǔ)上,利用3個月尺度SRI指數(shù)對氣候變化和人類活動影響下流域水文干旱演變進行分析。

    為了便于進行不同影響條件下季節(jié)性干旱變化分析,根據(jù)研究區(qū)地理和氣候特征將12個月依次劃分為春季(3—5月)、夏季(6—8月)、秋季(9—11月)和冬季(12月至次年2月)[25]。并繪制不同影響條件下不同季節(jié)SRI變化過程,見圖5。

    a)春季 b)夏季

    c)夏季 d)夏季

    春季,實測SRI和受氣候變化影響下的SRI在春季的變化趨勢存在較大差異,實測SRI主要呈現(xiàn)出輕度干旱化趨勢,受氣候變化影響下的SRI則呈現(xiàn)出無旱化趨勢。比較春季SRI值可知,前者共出現(xiàn)輕旱事件23次,在1996—2007年和2009—2020年出現(xiàn)連續(xù)春旱,后者僅2012年出現(xiàn)1次輕旱。說明在人類活動影響下,流域春季出現(xiàn)持續(xù)干旱化趨勢。這主要是由于瑪納斯河屬季節(jié)性河流,春季上游來水量少,而農(nóng)田播種又需要消耗大量的水資源,故使得流域春季呈現(xiàn)出干旱化趨勢。

    夏季,實測SRI和受氣候變化影響下的SRI在夏季的變化趨勢大致相同,總體上均呈下降趨勢。比較夏季SRI值可以發(fā)現(xiàn),在所取的研究時段內(nèi),兩者均無干旱事件發(fā)生,但后者SRI值波動幅度明顯大于前者。分析原因可知,膜下滴灌技術(shù)的實施降低了作物的棵間蒸發(fā),增加了土壤的濕度,使得流域夏季干旱化趨勢得到有效控制。

    秋季,實測SRI和受氣候變化影響下的SRI變化趨勢大致相似,主要呈現(xiàn)出無旱化趨勢。比較秋季SRI值可以發(fā)現(xiàn),前者在所取的研究時段內(nèi),無干旱事件發(fā)生;后者處于輕旱以上共出現(xiàn)8次,分別在1996—1998年、2005—2006年、2008年、2013年以及2019年出現(xiàn)秋旱,說明人類活動減少了流域干旱事件的發(fā)生。這主要得益于流域綠洲區(qū)人工綠洲面積的增長,改善了流域的生態(tài)環(huán)境,從而使得秋季無旱化趨勢更加穩(wěn)定。

    冬季,實測SRI和受氣候變化影響下的SRI值均在-0.5以下,干旱程度明顯。進一步分析比較冬季SRI值可知,前者處于中旱以上共出現(xiàn)3次,是在1996—1997年和2014年;后者處于中旱以上共出現(xiàn)18次,其中在1996—1997年、2001年、2011年和2019年出現(xiàn)了重度干旱。經(jīng)分析可知,冰川融雪水對瑪納斯河徑流有著重要的補給作用,而人類活動加速了冰川消融,提高了瑪納斯河冬季的徑流量,從而在短期內(nèi)改善了流域冬季的干旱程度。

    對發(fā)生的季節(jié)性干旱事件進行頻次和頻率統(tǒng)計可知(圖6),春季實測SRI以輕旱為主,無旱以及輕旱發(fā)生的頻次和頻率分別為2次、23次和8%和92%;春季氣候變化影響下SRI僅發(fā)生1次輕旱,其余均為無旱,發(fā)生輕旱和無旱的頻率分別為4%和96%。夏季實測SRI和氣候變化影響下SRI均無干旱事件發(fā)生,統(tǒng)計的25次均處于無旱。秋季實測SRI無干旱事件發(fā)生;秋季氣候變化影響下SRI無重旱事件發(fā)生,無旱、輕旱和中旱發(fā)生的頻次和頻率分別為17次、7次、1次和68%、28%、4%。冬季實測SRI發(fā)生3次中旱和21次輕旱,發(fā)生中旱和輕旱的頻率為12.5%和87.5%;冬季氣候變化影響下SRI發(fā)生6次輕旱、13次中旱和5次重旱,發(fā)生的頻率為25%、54.2%和20.8%。

    對比分析可知,實測SRI易在春季發(fā)生輕度干旱事件,而氣候變化影響下SRI則更易在冬季發(fā)生中旱及其以上干旱事件。可見人類活動一方面使流域春季出現(xiàn)干旱化趨勢,需要人類警惕;另一方面緩解了流域冬季的干旱程度,對于改善季節(jié)性干旱起到推動作用。

    a)實測

    b)受氣候變化影響

    3 討論

    流域徑流的形成和變化不僅與降水的變化有關(guān),而且還與人類活動引起的地表條件的變化以及人類活動的直接影響有關(guān)[26]。本研究從氣候變化和人類活動的角度出發(fā),定量分析了這2種影響因素對流域徑流的影響程度,并通過對比分析探討了在這2種影響因素下流域水文干旱演變情況。流域徑流量的增加受氣候變化和人類活動的雙重影響,相對于人類活動所產(chǎn)生的影響,氣候變化對流域流量變化的影響量為2.95 mm,貢獻率為56.73%,而人類活動對徑流量變化的影響量為2.24 mm,貢獻率為43.37%,可見,氣候變化對流域流量變化的影響更加顯著。這與陳伏龍等[15]、高培等[27]的研究結(jié)果保持一致。同時需要指出的是,雖然降水-徑流雙累積關(guān)系模型引入了累積量以減少年際波動對測量數(shù)據(jù)的影響且年與累積之間的線性擬合具有較高的相關(guān)系數(shù)(R2=0.99),但本研究僅將降水視為氣候因素,實際上應(yīng)包括日照,溫度,風(fēng)速和其他因素[28],因此對氣候因素和人類活動的分解缺乏徹底性,定量研究仍存在不確定性。

    Mann-Kendall非參數(shù)統(tǒng)計檢驗和降水-徑流雙累積曲線表明瑪納斯河流域的徑流量在1995年發(fā)生突變,這與流域?qū)嵭泄?jié)水灌溉的時間相符,流域內(nèi)人口和耕地面積大幅增加(圖7),人類活動對徑流的影響加劇。另外瑪納斯河屬冰川融雪及降雨混合補給型河流,由于人類對自然資源不合理的開發(fā)利用,導(dǎo)致全球氣溫持續(xù)升高,山區(qū)冰川消融加劇。研究[29]發(fā)現(xiàn)瑪納斯河流域已有 32.5%的冰川面積消失,雖然短期內(nèi)冰川融雪水對徑流補給量的增加能夠降低流域干旱事件的發(fā)生,但從長遠來看冰川融水量先增后減的拐點一旦出現(xiàn),以降水和冰川融水為主要補給的瑪納斯河徑流量會出現(xiàn)陡降,在枯水季節(jié)或年份可能出現(xiàn)區(qū)域性缺水危機。同時由于人類活動的影響,現(xiàn)階段瑪納斯河流域春季呈現(xiàn)出持續(xù)干旱化趨勢,針對這一現(xiàn)象,瑪納斯河流域應(yīng)合理控制開墾面積,結(jié)合研究區(qū)實際情況建議調(diào)整農(nóng)作物種植結(jié)構(gòu),減少高耗水農(nóng)作物種植面積,發(fā)展低耗水農(nóng)作物。更為重要的是,應(yīng)大力培養(yǎng)高效農(nóng)業(yè)技術(shù)管理人員,盡管研究區(qū)高效節(jié)水技術(shù)已經(jīng)發(fā)展的相當(dāng)成熟,但相應(yīng)的技術(shù)管理水平還有待進一步提高,因此,培養(yǎng)高效農(nóng)業(yè)技術(shù)管理人員是接下來發(fā)展的首要任務(wù)。

    注:橙色實虛線分別代表1975—1995、1996—2020年的平均值。

    4 結(jié)論

    a)近65 a來瑪納斯河流域徑流量總體上呈增加趨勢,Mann-Kendall突變檢驗法和降水-徑流雙累積曲線法識別出流域年徑流在1995年左右發(fā)生突變,2005—2020年徑流量呈顯著性增加趨勢。

    b)氣候變化是引起瑪納斯河流域徑流增加的主要原因,對流域月徑流量增加的貢獻率為56.73%;人類活動主要通過流域開墾及農(nóng)業(yè)灌溉等方式間接影響徑流變化,對流域月徑流量增加的貢獻率為43.37%。

    c)通過利用SRI對變化期徑流進行季尺度分析發(fā)現(xiàn),相對于受單一氣候變化的影響,人類活動降低了流域干旱事件發(fā)生的頻率,短期內(nèi)改善了流域的干旱狀況;但人類活動同樣打破了原有的水文生態(tài)平衡,使得流域春季出現(xiàn)了干旱化趨勢。

    猜你喜歡
    人類影響活動
    “六小”活動
    少先隊活動(2022年5期)2022-06-06 03:45:04
    “活動隨手拍”
    是什么影響了滑動摩擦力的大小
    人類能否一覺到未來?
    行動不便者,也要多活動
    中老年保健(2021年2期)2021-08-22 07:31:10
    哪些顧慮影響擔(dān)當(dāng)?
    人類第一殺手
    好孩子畫報(2020年5期)2020-06-27 14:08:05
    1100億個人類的清明
    三八節(jié),省婦聯(lián)推出十大系列活動
    海峽姐妹(2018年3期)2018-05-09 08:20:40
    擴鏈劑聯(lián)用對PETG擴鏈反應(yīng)與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    精品一区二区三区视频在线| 中文在线观看免费www的网站| 观看免费一级毛片| 国产免费一级a男人的天堂| 免费人成在线观看视频色| 丰满少妇做爰视频| 午夜免费男女啪啪视频观看| 99热这里只有是精品在线观看| 人人澡人人妻人| 99视频精品全部免费 在线| 成人午夜精彩视频在线观看| 三级国产精品片| 91精品国产国语对白视频| 久久久国产一区二区| 亚洲精品,欧美精品| 看十八女毛片水多多多| 国产乱人偷精品视频| 亚洲av免费高清在线观看| 丰满乱子伦码专区| 看非洲黑人一级黄片| 国产av精品麻豆| 国产色爽女视频免费观看| a级毛色黄片| 卡戴珊不雅视频在线播放| 夫妻午夜视频| 成人黄色视频免费在线看| 欧美xxⅹ黑人| 亚洲久久久国产精品| 99热这里只有是精品50| 亚洲av综合色区一区| 国产视频内射| 熟妇人妻不卡中文字幕| 久久精品国产a三级三级三级| av线在线观看网站| 精品人妻一区二区三区麻豆| 少妇 在线观看| 免费在线观看成人毛片| 性色av一级| 久久6这里有精品| 国产黄频视频在线观看| 观看免费一级毛片| 亚洲欧美日韩另类电影网站| 亚洲,一卡二卡三卡| 99久久精品热视频| 人妻夜夜爽99麻豆av| 大香蕉久久网| 国产成人精品一,二区| 欧美3d第一页| 色婷婷久久久亚洲欧美| 黄色毛片三级朝国网站 | 国产免费一区二区三区四区乱码| 最近手机中文字幕大全| 国语对白做爰xxxⅹ性视频网站| 丰满迷人的少妇在线观看| 各种免费的搞黄视频| 少妇被粗大猛烈的视频| 久久久亚洲精品成人影院| 黄色日韩在线| 国产伦理片在线播放av一区| 久久国产精品男人的天堂亚洲 | 免费播放大片免费观看视频在线观看| 男女边摸边吃奶| 人人妻人人澡人人爽人人夜夜| 久久亚洲国产成人精品v| 精品亚洲成国产av| 亚洲av不卡在线观看| 久久精品夜色国产| 我要看黄色一级片免费的| 一区二区三区四区激情视频| 亚洲欧洲日产国产| av免费在线看不卡| 国产极品粉嫩免费观看在线 | 中文资源天堂在线| 国产成人精品一,二区| 久久精品国产a三级三级三级| 久久精品国产亚洲网站| 国产老妇伦熟女老妇高清| 一级毛片我不卡| 美女国产视频在线观看| 日本与韩国留学比较| 成人毛片a级毛片在线播放| av天堂中文字幕网| 精品少妇黑人巨大在线播放| 另类亚洲欧美激情| 国产91av在线免费观看| 韩国高清视频一区二区三区| 男人爽女人下面视频在线观看| 王馨瑶露胸无遮挡在线观看| 午夜免费鲁丝| 久久久久人妻精品一区果冻| 黑人高潮一二区| 日本wwww免费看| 欧美日韩精品成人综合77777| 亚洲精品日韩在线中文字幕| 大话2 男鬼变身卡| 欧美区成人在线视频| 99国产精品免费福利视频| av视频免费观看在线观看| 久久国产精品男人的天堂亚洲 | 国产精品国产三级国产av玫瑰| 亚洲自偷自拍三级| 又粗又硬又长又爽又黄的视频| 草草在线视频免费看| 69精品国产乱码久久久| 一区二区三区精品91| 在线观看av片永久免费下载| 极品教师在线视频| 国产精品一二三区在线看| 简卡轻食公司| 国产精品嫩草影院av在线观看| 插阴视频在线观看视频| 丰满少妇做爰视频| 精华霜和精华液先用哪个| 成人毛片60女人毛片免费| 肉色欧美久久久久久久蜜桃| .国产精品久久| 女人久久www免费人成看片| 国产精品99久久久久久久久| 91精品一卡2卡3卡4卡| 成人午夜精彩视频在线观看| 亚洲丝袜综合中文字幕| 久久狼人影院| 中文字幕久久专区| 亚洲第一av免费看| 欧美精品高潮呻吟av久久| 国产乱人偷精品视频| 精品亚洲乱码少妇综合久久| 国产乱人偷精品视频| 国产乱人偷精品视频| 高清在线视频一区二区三区| 久久精品国产a三级三级三级| 老熟女久久久| 老熟女久久久| 国产乱人偷精品视频| av在线app专区| 久久久国产欧美日韩av| 少妇精品久久久久久久| 国产日韩欧美亚洲二区| 91成人精品电影| 久久久国产欧美日韩av| 久久久久久久亚洲中文字幕| 国产成人精品一,二区| www.av在线官网国产| 纵有疾风起免费观看全集完整版| 天天操日日干夜夜撸| 中文天堂在线官网| 一级,二级,三级黄色视频| 91午夜精品亚洲一区二区三区| 亚洲天堂av无毛| 国产探花极品一区二区| 精品视频人人做人人爽| 亚洲丝袜综合中文字幕| 成人午夜精彩视频在线观看| 国产成人午夜福利电影在线观看| 熟女av电影| 国产欧美日韩精品一区二区| 亚洲第一av免费看| 在线观看人妻少妇| 亚洲av成人精品一区久久| 精品久久久噜噜| 日韩精品免费视频一区二区三区 | 国产淫语在线视频| av在线老鸭窝| 又爽又黄a免费视频| 国产91av在线免费观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲av中文av极速乱| 天天躁夜夜躁狠狠久久av| 亚洲精品乱码久久久久久按摩| 你懂的网址亚洲精品在线观看| 蜜桃在线观看..| 精品亚洲乱码少妇综合久久| 欧美日韩亚洲高清精品| 99久久人妻综合| 老熟女久久久| av黄色大香蕉| 伊人久久国产一区二区| 一个人免费看片子| 精品一区二区三卡| 亚洲精华国产精华液的使用体验| 国产精品国产av在线观看| 亚洲精品一二三| 麻豆成人av视频| 一个人看视频在线观看www免费| 日本黄色日本黄色录像| www.av在线官网国产| 夫妻性生交免费视频一级片| 久久久久久久久久人人人人人人| 精品国产露脸久久av麻豆| 国产有黄有色有爽视频| 欧美xxⅹ黑人| 午夜日本视频在线| 亚洲真实伦在线观看| 亚洲熟女精品中文字幕| 久久国产乱子免费精品| 精华霜和精华液先用哪个| 欧美bdsm另类| 亚洲av成人精品一二三区| 久久精品久久久久久噜噜老黄| 人人妻人人澡人人看| 桃花免费在线播放| 最近最新中文字幕免费大全7| 大香蕉久久网| 亚洲性久久影院| 一个人看视频在线观看www免费| 欧美激情极品国产一区二区三区 | 成人特级av手机在线观看| 欧美少妇被猛烈插入视频| 国产日韩一区二区三区精品不卡 | 亚洲国产精品一区二区三区在线| 男的添女的下面高潮视频| 人妻夜夜爽99麻豆av| 日日啪夜夜撸| 精品国产乱码久久久久久小说| 99热国产这里只有精品6| 大码成人一级视频| 国产伦精品一区二区三区视频9| 亚洲精品,欧美精品| 欧美bdsm另类| 国产极品粉嫩免费观看在线 | 夜夜爽夜夜爽视频| 老女人水多毛片| 高清毛片免费看| 欧美一级a爱片免费观看看| 国产欧美日韩精品一区二区| 人人妻人人爽人人添夜夜欢视频 | av播播在线观看一区| 久久人人爽人人爽人人片va| 亚洲av电影在线观看一区二区三区| 能在线免费看毛片的网站| 久久久久久久国产电影| 晚上一个人看的免费电影| 国产91av在线免费观看| 亚洲av在线观看美女高潮| 国产乱来视频区| 另类精品久久| 亚洲精品国产av蜜桃| 国产精品麻豆人妻色哟哟久久| 一个人看视频在线观看www免费| 大话2 男鬼变身卡| 观看av在线不卡| 国产亚洲91精品色在线| 视频中文字幕在线观看| 免费观看的影片在线观看| 国产亚洲一区二区精品| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av涩爱| 亚洲第一区二区三区不卡| 国产精品偷伦视频观看了| 99久久中文字幕三级久久日本| a级毛片免费高清观看在线播放| av国产精品久久久久影院| 免费黄色在线免费观看| 99热这里只有精品一区| 男女无遮挡免费网站观看| 欧美日本中文国产一区发布| 男女免费视频国产| 99久久精品热视频| 国产精品女同一区二区软件| 我要看黄色一级片免费的| 国产视频内射| 99热这里只有是精品50| 内射极品少妇av片p| 天堂中文最新版在线下载| 2022亚洲国产成人精品| 精品国产一区二区久久| a级毛片免费高清观看在线播放| 欧美日韩国产mv在线观看视频| 妹子高潮喷水视频| 国产淫语在线视频| 久久国产精品大桥未久av | 美女内射精品一级片tv| 97在线视频观看| 3wmmmm亚洲av在线观看| 国产一区二区在线观看日韩| 大片电影免费在线观看免费| 一二三四中文在线观看免费高清| 国产精品嫩草影院av在线观看| 午夜精品国产一区二区电影| 国产精品不卡视频一区二区| 亚洲欧美一区二区三区黑人 | 欧美日韩在线观看h| 亚洲国产欧美在线一区| 亚洲婷婷狠狠爱综合网| 国产日韩欧美视频二区| 青春草视频在线免费观看| 免费av不卡在线播放| 一级,二级,三级黄色视频| 草草在线视频免费看| 日本av手机在线免费观看| 久久精品夜色国产| 精品亚洲成a人片在线观看| 欧美变态另类bdsm刘玥| 国产白丝娇喘喷水9色精品| 一级,二级,三级黄色视频| 日韩中文字幕视频在线看片| 国产日韩一区二区三区精品不卡 | 男的添女的下面高潮视频| 亚洲av.av天堂| 黄色毛片三级朝国网站 | 制服丝袜香蕉在线| 久久国产乱子免费精品| 激情五月婷婷亚洲| 51国产日韩欧美| 国产精品麻豆人妻色哟哟久久| 久久女婷五月综合色啪小说| 深夜a级毛片| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡 | 大片电影免费在线观看免费| 人妻人人澡人人爽人人| 99热这里只有精品一区| 欧美+日韩+精品| 久久久午夜欧美精品| 夫妻午夜视频| 老司机影院成人| 又黄又爽又刺激的免费视频.| 夫妻性生交免费视频一级片| 成人国产麻豆网| av卡一久久| 熟妇人妻不卡中文字幕| 欧美激情极品国产一区二区三区 | 国产精品伦人一区二区| 欧美老熟妇乱子伦牲交| 亚洲av国产av综合av卡| 亚洲欧洲日产国产| 国产 精品1| 久久这里有精品视频免费| 国产女主播在线喷水免费视频网站| 18禁裸乳无遮挡动漫免费视频| 春色校园在线视频观看| 欧美精品高潮呻吟av久久| 欧美日韩一区二区视频在线观看视频在线| 少妇人妻 视频| 六月丁香七月| 久久精品国产亚洲av天美| 国产在线免费精品| 久久久久精品久久久久真实原创| 日日爽夜夜爽网站| 中文字幕亚洲精品专区| 黄色一级大片看看| 69精品国产乱码久久久| 日日爽夜夜爽网站| 蜜桃久久精品国产亚洲av| 国产在线视频一区二区| 美女国产视频在线观看| 一区二区三区免费毛片| 又爽又黄a免费视频| 九九在线视频观看精品| 大码成人一级视频| 在线观看www视频免费| 夜夜看夜夜爽夜夜摸| 国产精品女同一区二区软件| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 热re99久久国产66热| 成年女人在线观看亚洲视频| 国产精品偷伦视频观看了| 最后的刺客免费高清国语| 久久久欧美国产精品| 一区二区av电影网| 亚洲欧美一区二区三区国产| 伦理电影大哥的女人| 一区二区av电影网| 简卡轻食公司| 国产精品.久久久| 丰满人妻一区二区三区视频av| 久久人人爽人人片av| 三级经典国产精品| 亚洲内射少妇av| 美女视频免费永久观看网站| 精品熟女少妇av免费看| 亚洲情色 制服丝袜| 国产乱人偷精品视频| 国产黄色免费在线视频| 亚洲婷婷狠狠爱综合网| 久热这里只有精品99| 国产极品粉嫩免费观看在线 | 一级毛片 在线播放| 日韩亚洲欧美综合| 国产精品麻豆人妻色哟哟久久| 亚洲经典国产精华液单| 国产黄色免费在线视频| 久久人人爽人人爽人人片va| 国产亚洲午夜精品一区二区久久| 美女xxoo啪啪120秒动态图| 久久99热这里只频精品6学生| 最后的刺客免费高清国语| 国产精品久久久久久精品古装| 伦理电影大哥的女人| 国产精品人妻久久久影院| 国产精品久久久久成人av| a级一级毛片免费在线观看| 国产伦精品一区二区三区视频9| 久久久久国产精品人妻一区二区| 永久网站在线| 亚洲精品久久午夜乱码| 免费人成在线观看视频色| 极品人妻少妇av视频| 国产精品99久久99久久久不卡 | 久久久久国产精品人妻一区二区| 国产在线免费精品| 亚洲欧洲日产国产| 一个人免费看片子| 国产一区二区在线观看日韩| 一级黄片播放器| 你懂的网址亚洲精品在线观看| 亚洲精品久久午夜乱码| 久久婷婷青草| 亚洲欧美成人综合另类久久久| 免费黄色在线免费观看| 成人亚洲欧美一区二区av| 9色porny在线观看| 成年女人在线观看亚洲视频| 精品国产一区二区久久| 在线观看www视频免费| 纵有疾风起免费观看全集完整版| 91久久精品国产一区二区成人| 欧美高清成人免费视频www| 国产一区亚洲一区在线观看| 18禁动态无遮挡网站| 男女边摸边吃奶| 中文字幕精品免费在线观看视频 | 99九九在线精品视频 | 美女xxoo啪啪120秒动态图| 亚洲国产毛片av蜜桃av| 成人国产麻豆网| 青春草国产在线视频| 亚洲av日韩在线播放| 99久久综合免费| 国产精品不卡视频一区二区| 亚洲av福利一区| videos熟女内射| 少妇被粗大猛烈的视频| 国产精品99久久99久久久不卡 | 亚洲精品一区蜜桃| 观看av在线不卡| 性色av一级| 男人舔奶头视频| 午夜老司机福利剧场| 成人免费观看视频高清| 搡老乐熟女国产| 欧美97在线视频| 最新中文字幕久久久久| 久久人人爽人人片av| 我的老师免费观看完整版| 高清av免费在线| 久久久久国产精品人妻一区二区| 成人免费观看视频高清| 天堂8中文在线网| 日韩精品免费视频一区二区三区 | 麻豆成人午夜福利视频| 色网站视频免费| 高清不卡的av网站| 久久亚洲国产成人精品v| 一级a做视频免费观看| av在线老鸭窝| 日韩不卡一区二区三区视频在线| a级毛色黄片| 一级,二级,三级黄色视频| 色婷婷久久久亚洲欧美| 黑丝袜美女国产一区| 青青草视频在线视频观看| 少妇精品久久久久久久| 极品人妻少妇av视频| 欧美日韩亚洲高清精品| 国产黄频视频在线观看| 乱系列少妇在线播放| 全区人妻精品视频| 国精品久久久久久国模美| 国产精品不卡视频一区二区| 多毛熟女@视频| 国产在线一区二区三区精| av播播在线观看一区| 日韩 亚洲 欧美在线| 日本黄大片高清| 久久6这里有精品| 另类亚洲欧美激情| 亚洲国产欧美在线一区| 久久精品久久久久久久性| 日韩一本色道免费dvd| 春色校园在线视频观看| av天堂中文字幕网| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看| 国产亚洲午夜精品一区二区久久| 少妇精品久久久久久久| 日日撸夜夜添| 久久久久久久久久久丰满| 欧美bdsm另类| 香蕉精品网在线| 日本与韩国留学比较| 精品人妻熟女av久视频| 伊人亚洲综合成人网| 成年人午夜在线观看视频| 亚洲欧美一区二区三区黑人 | 精品一区二区三区视频在线| 国产av码专区亚洲av| 国产免费一区二区三区四区乱码| 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 97超碰精品成人国产| 国产亚洲午夜精品一区二区久久| 精品久久久噜噜| 精品午夜福利在线看| 亚洲第一av免费看| 啦啦啦在线观看免费高清www| 日本欧美视频一区| xxx大片免费视频| 一级爰片在线观看| 丰满人妻一区二区三区视频av| 涩涩av久久男人的天堂| 日韩不卡一区二区三区视频在线| 中文字幕制服av| 免费观看a级毛片全部| 亚洲成人一二三区av| 亚洲欧美精品专区久久| 另类精品久久| 国产一区二区三区av在线| 性色av一级| 日本与韩国留学比较| 久久精品国产亚洲网站| 国产精品国产三级专区第一集| 亚洲欧美一区二区三区国产| 男女啪啪激烈高潮av片| 老司机影院成人| 夫妻性生交免费视频一级片| 能在线免费看毛片的网站| 亚洲经典国产精华液单| 王馨瑶露胸无遮挡在线观看| 老司机影院毛片| 国模一区二区三区四区视频| 国产在线视频一区二区| 国产综合精华液| 精品一品国产午夜福利视频| 日本wwww免费看| 乱人伦中国视频| 精品一区二区免费观看| 丁香六月天网| 爱豆传媒免费全集在线观看| 国产一区有黄有色的免费视频| 777米奇影视久久| 免费少妇av软件| 色哟哟·www| 国产又色又爽无遮挡免| 亚洲国产精品一区三区| 99国产精品免费福利视频| 亚洲色图综合在线观看| 性高湖久久久久久久久免费观看| 美女内射精品一级片tv| 最近中文字幕2019免费版| 六月丁香七月| 尾随美女入室| 嫩草影院新地址| 久久久久国产网址| 久久国产亚洲av麻豆专区| 看十八女毛片水多多多| av卡一久久| 国产高清不卡午夜福利| 性色avwww在线观看| 一级毛片黄色毛片免费观看视频| 日韩欧美精品免费久久| 国产精品一区二区在线观看99| videossex国产| 九九在线视频观看精品| 一本一本综合久久| a级毛片在线看网站| 国产视频内射| av福利片在线| 国产日韩欧美亚洲二区| 大码成人一级视频| 国产av一区二区精品久久| 狂野欧美激情性xxxx在线观看| 内射极品少妇av片p| 91精品伊人久久大香线蕉| 全区人妻精品视频| 黄色视频在线播放观看不卡| 欧美亚洲 丝袜 人妻 在线| 美女视频免费永久观看网站| 亚洲av成人精品一二三区| av在线老鸭窝| 欧美国产精品一级二级三级 | 色吧在线观看| 男女无遮挡免费网站观看| 国产视频首页在线观看| 国产探花极品一区二区| 91精品伊人久久大香线蕉| 成人亚洲欧美一区二区av| 国产精品久久久久成人av| 国产女主播在线喷水免费视频网站| 熟女人妻精品中文字幕| 国产极品天堂在线| 久久精品夜色国产| 国产永久视频网站| 国产乱来视频区| 欧美bdsm另类| av在线app专区| 亚洲内射少妇av| 蜜桃在线观看..| 亚洲国产精品一区三区| 国产在线视频一区二区| 国产成人a∨麻豆精品| 老司机影院毛片| 国产成人a∨麻豆精品| 国产高清不卡午夜福利| 春色校园在线视频观看| 午夜福利在线观看免费完整高清在| 男人和女人高潮做爰伦理| 日韩av不卡免费在线播放| 亚洲精品中文字幕在线视频 | 国产成人免费无遮挡视频| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩另类电影网站|