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

    拉薩河流域徑流對(duì)土地利用和氣候變化的響應(yīng)分析

    2021-10-21 13:28:02霍軍軍伊明啟江玉吉
    關(guān)鍵詞:拉薩河徑流氣候變化

    霍軍軍,伊明啟,2,王 靜,江玉吉

    (1.長(zhǎng)江科學(xué)院 水資源綜合利用研究所,武漢 430010; 2.河海大學(xué) 水文水資源學(xué)院,南京 210098;3.西藏自治區(qū)水文水資源勘測(cè)局,拉薩 850000)

    1 研究背景

    全球氣候在過去一個(gè)多世紀(jì)發(fā)生了顯著變化,并且將在未來一段時(shí)間內(nèi)持續(xù)這種變化[1]。青藏高原平均海拔在4 000 m以上,地形獨(dú)特復(fù)雜,對(duì)氣候變化響應(yīng)迅速,是全球氣候最敏感的區(qū)域之一[2],且對(duì)全球氣候變化趨勢(shì)具有一定的超前性[3],例如其氣候變暖速率自1960年以來為同時(shí)期全球變暖速率的2倍左右[4]。在氣候變化對(duì)水文要素的影響中,徑流的變化最為直觀,但徑流同時(shí)也受土地利用變化等因素的影響。因此,如何量化分析徑流對(duì)氣候變化的響應(yīng)是目前的熱點(diǎn)研究?jī)?nèi)容之一。

    目前相關(guān)研究表明,土地利用和氣候變化是影響流域徑流的主要驅(qū)動(dòng)因素。Wang等[5]利用SWAT(Soil and Water Assessment Tool)模型對(duì)媯水河流域徑流進(jìn)行模擬分析,認(rèn)為氣候變化對(duì)該流域徑流量的變化貢獻(xiàn)率達(dá)到60%;Zhang等[6]也基于SWAT模型在格蘭德河和泰晤士河對(duì)未來氣候情景下的徑流響應(yīng)進(jìn)行了研究,認(rèn)為氣候變化導(dǎo)致降水時(shí)空變異,進(jìn)而加劇徑流的不均勻性。土地利用變化通過改變下墊面條件,也會(huì)對(duì)徑流過程產(chǎn)生較大影響[7]。Chen等[8]運(yùn)用SWAT和HSPF模型開展了不同時(shí)間尺度上土地利用變化對(duì)徑流影響的不確定性分析;Tankpa等[9]在利用CA-Markov模型預(yù)測(cè)未來土地利用變化的基礎(chǔ)上,基于SWAT模型對(duì)Ashi河流域土地利用變化對(duì)徑流的影響進(jìn)行了綜合分析。在不同的區(qū)域,氣候與土地利用變化的程度和速度存在差異,因而對(duì)徑流變化的影響程度也各有不同。在一些區(qū)域,相較于土地利用的變化,氣候變化對(duì)徑流的影響程度更大[10];而在另一些區(qū)域,土地利用對(duì)徑流的影響程度更大[11]。

    拉薩河流域地處青藏高原中南部,水量較為豐沛,人口和耕地較為密集,是西藏自治區(qū)的政治和文化中心。目前對(duì)拉薩河流域的相關(guān)研究側(cè)重在氣溫、降水、冰川凍土等氣候要素的變化分析上[12-13],在土地利用、氣候變化等對(duì)徑流的影響分析方面較少。本文在利用SWAT模型對(duì)拉薩河流域徑流進(jìn)行模擬的基礎(chǔ)上,通過比較不同情景模式下的徑流量變化情況,研究拉薩河徑流變化的敏感因素,并量化分析土地利用和氣候變化對(duì)徑流的影響。

    2 研究區(qū)域與數(shù)據(jù)

    2.1 研究區(qū)概況

    拉薩河(90°05′E—93°20′E,29°20′N—31°15′N)發(fā)源于念青唐古拉山中段南麓,流域面積32 610 km2,河長(zhǎng)551 km,平均海拔4 900 m左右。拉薩河流域?qū)俑咴瓬貛Ц珊蛋敫珊禋夂騾^(qū),植被類型主要有山地灌叢草原和高山草原,徑流來源主要為降水和積雪融水。受西南季風(fēng)影響,流域干濕季節(jié)明顯,5—9月份為雨季,10月—次年4月份則為旱季,多年平均降水量約為650 mm。全年太陽輻射強(qiáng)度較大,日間氣溫變化幅度較大。流域內(nèi)共設(shè)有4個(gè)水文站,分別為羊八井站、唐加站、旁多站與拉薩站,流域內(nèi)及周邊共分布有氣象站5個(gè),分別為尼木站、當(dāng)雄站、拉薩站、澤當(dāng)站和嘉黎站,具體信息詳見表1和圖1。

    表1 拉薩河流域水文站和氣象站基本信息Table 1 Basic information of hydrological andmeteorological stations in Lhasa River Basin

    圖1 拉薩河流域地理位置及流域內(nèi)水文氣象站點(diǎn)分布Fig.1 Map of the geographical location of the LhasaRiver Basin and hydro-meteorological stations inthe basin

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

    為建立拉薩河流域徑流模型,本文采用的數(shù)據(jù)包括數(shù)字高程模型(Digital Elevation Model,DEM)數(shù)據(jù)、土地利用數(shù)據(jù)、土壤數(shù)據(jù)、氣象數(shù)據(jù)和水文數(shù)據(jù)等。DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/),空間分辨率為90 m。土地利用數(shù)據(jù)來源于中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(https:∥www.resdc.cn/),本文用到了2000年和2015年分辨率為1 km的土地利用數(shù)據(jù),其土地利用類型包括林地、耕地、草地、水域、居民地和未利用土地等6個(gè)一級(jí)大類和25個(gè)二級(jí)細(xì)類組成;在此基礎(chǔ)上按SWAT模型要求重分類為耕地、林地、草地、水域、城鎮(zhèn)居民區(qū)、農(nóng)村居民區(qū)、工業(yè)用地、荒地裸地與濕地共9個(gè)類型。土壤類型數(shù)據(jù)來源于聯(lián)合國(guó)糧農(nóng)組織(FAO)發(fā)布的世界土壤數(shù)據(jù)庫(Harmonized World Soil Database version1.1,HWSD)(http:∥westdc.westgis.ac.cn/),空間分辨率為1 km。以上數(shù)據(jù)在本文應(yīng)用中都統(tǒng)一到1 km的空間分辨率。

    氣象數(shù)據(jù)來源為中國(guó)氣象數(shù)據(jù)網(wǎng)(http:∥data.cma.cn/),包括尼木站、當(dāng)雄站、拉薩站、澤當(dāng)站與嘉黎站1976—2016年日降水、氣溫、風(fēng)速、相對(duì)濕度與太陽輻射數(shù)據(jù)。水文數(shù)據(jù)為羊八井站、旁多站、唐加站與拉薩站1976—2016年日徑流數(shù)據(jù)。

    3 研究方法

    3.1 SWAT 模型與徑流模擬

    利用以上地理空間數(shù)據(jù)和氣象水文數(shù)據(jù),構(gòu)建拉薩河流域徑流SWAT模型。該模型將拉薩河流域劃分為27個(gè)子流域,繼而細(xì)分為107個(gè)水量平衡單元。SWAT模型參數(shù)較多,本文根據(jù)相關(guān)研究及區(qū)域特點(diǎn)共選取了11個(gè)參數(shù)進(jìn)行率定和驗(yàn)證,分別為SOL_AWC、CN2、SOL_Z、GW_DELAY、ALPHA_BF、GWQMN、SOL_K、ESCO、CANMX、GW_REVAP和EPCO。模型率定期為1990—1996年,驗(yàn)證期為1997—2003年,并選取決定系數(shù)(R2)、納什效率系數(shù)(Nash-Sutcliffe efficiency,NSE)、誤差百分?jǐn)?shù)(Percent Bias,PBIAS)、均方根誤差與標(biāo)準(zhǔn)誤差比率(Ratio of the root mean square error to the standard deviation of measured data,RSR)用于評(píng)價(jià)模擬精度,評(píng)價(jià)指標(biāo)計(jì)算公式如式(1)—式(4)所示。

    (2)

    (3)

    (4)

    式中:Q1為實(shí)測(cè)徑流量(m3/s);Q2為模擬徑流量(m3/s);Q1m為實(shí)測(cè)平均徑流量(m3/s);Q2m為模擬平均徑流量(m3/s);n為數(shù)據(jù)總個(gè)數(shù)。

    R2和NSE在[0,1]內(nèi)取值,越接近于1,表明擬合程度越好。一般情況下,當(dāng)R2>0.7,NSE>0.6時(shí),模擬結(jié)果與實(shí)測(cè)數(shù)據(jù)擬合情況較好。PBIAS和RSR可對(duì)模型模擬值與實(shí)測(cè)值偏差大小進(jìn)行評(píng)價(jià),一般認(rèn)為PBIAS絕對(duì)值<15%,RSR<0.7時(shí),模擬效果較好。

    我把原來存在舊手機(jī)電池里的電放光,以測(cè)試能否完全靠太陽能工作。這個(gè)應(yīng)急燈就放在客廳的茶幾上,電池板朝向南方,沒有放在陽光直射的窗臺(tái)上。開始的幾天,我每天晚上測(cè)試使用兩個(gè)小時(shí),可以正常工作很多天,后來嘗試了一次整晚點(diǎn)亮應(yīng)急燈,到第二天早晨都能正常發(fā)光。如果把它放到陽臺(tái)上的陽光下,應(yīng)該可以充更多的電。測(cè)試結(jié)果表明,使用這個(gè)家用應(yīng)急燈時(shí),只要放在光線較好的地方,停電時(shí)就可以使用很長(zhǎng)時(shí)間,平時(shí)完全不用惦記充電的事。

    本文選用SWAT-CUP中SUFI-2程序分析模型參數(shù)的全局靈敏度,并采用p值和t值對(duì)靈敏度進(jìn)行評(píng)價(jià)。某一參數(shù)的p值越小,t值絕對(duì)值越大,代表該參數(shù)越靈敏[14]。

    3.2 情景設(shè)置

    本文通過對(duì)土地利用條件與氣候因子設(shè)置不同的情景并進(jìn)行交叉組合,來分析拉薩河流域徑流在不同因素變化及協(xié)同作用下的響應(yīng)情況。

    3.2.1 土地利用變化與氣候因素共同作用情景設(shè)置

    從當(dāng)雄站、嘉黎站、尼木站、拉薩站和澤當(dāng)站1985—2016年間的年均降水、年均氣溫?cái)?shù)據(jù)可以看出,2000—2016年期間年均降水與氣溫較1985—2000年都顯著增加,其中嘉黎站年均降水增加達(dá)22.04 mm,拉薩站年均氣溫升高達(dá)1.17 ℃。因此,為探求拉薩河流域土地利用變化與氣候因素對(duì)徑流的綜合影響,本文擬以2000年為界將1985—2000年、2001—2016年氣象數(shù)據(jù)與2000年、2015年土地利用情況進(jìn)行交叉組合,共設(shè)置S1—S4四種情景,如表2所示。通過比較S1與S3、S2與S4情景結(jié)果以分析土地利用變化情況對(duì)徑流的影響,比較S1與S2、S3與S4情況結(jié)果以分析氣候變化對(duì)徑流的影響,比較S1與S4情景結(jié)果以分析土地利用和氣候變化協(xié)同作用下徑流的變化情況。

    表2 土地利用和氣候變化情景設(shè)置Table 2 Scenarios of land use and climate change

    3.2.2 氣候變化情景設(shè)置

    參考拉薩河流域近30 a氣溫和降水的變化趨勢(shì),氣候情景設(shè)置以1984—2016年實(shí)測(cè)氣溫和降水?dāng)?shù)據(jù)為基礎(chǔ),降水增減10%、20%(P1—P4),氣溫增減1、2 ℃(T1—T4),并交叉組合為情景PT1—PT8,如表3所示。

    表3 氣溫和降水變化情景設(shè)置Table 3 Scenarios for the combination of temperatureand precipitation changes

    4 結(jié)果與分析

    4.1 模擬效果評(píng)價(jià)

    凍土變化通過影響地表徑流、土壤含水量、地下水儲(chǔ)量等因素而使流域水量平衡發(fā)生變化[15]。青藏高原凍土分布較為廣泛,目前對(duì)青藏高原不同區(qū)域的研究結(jié)果表明,凍土退化是徑流變化的驅(qū)動(dòng)因素之一[16-17]。對(duì)于拉薩河流域,Lin等[18]相關(guān)研究表明,凍土變化不是影響徑流的控制性因素。根據(jù)國(guó)家青藏高原科學(xué)數(shù)據(jù)中心(http:∥data.tpdc.ac.cn)發(fā)布的多年凍土分布概率(Permafrost Zonation Index,PZI)結(jié)果(圖2(a)),拉薩河流域無凍土分布(PZI≤0.01)的區(qū)域面積約占41.5%,常年分布凍土(PZI≥0.5)(圖2(b))的區(qū)域面積約占27.5%,凍土年際變化的區(qū)域占比較小。凍土對(duì)拉薩河流域徑流的影響體現(xiàn)在冬季[16],而拉薩河冬季徑流量約占年徑流量的14%,比例較小,因此凍土變化對(duì)拉薩河流域徑流的總體影響較小,故在SWAT模型中未著重考慮凍土因素。

    圖2 拉薩河流域凍土分布情況Fig.2 Distribution of permafrost in the Lhasa Riverbasin

    圖3 拉薩河流域徑流模擬效果評(píng)價(jià)Fig.3 Evaluation results of runoff simulation effect in Lhasa River Basin

    表4 拉薩河流域SWAT模型參數(shù)敏感性分析Table 4 Sensitivity analysis of SWAT model parametersin Lhasa River Basin

    4.2 徑流對(duì)土地利用和氣候因素變化的響應(yīng)分析

    4.2.1 土地利用變化情況

    拉薩河流域2000年和2015年土地利用類型分布情況如圖4所示,2000—2015年的土地利用類型變化轉(zhuǎn)移矩陣如表5所示。2000—2015年拉薩河流域主要土地利用類型為牧場(chǎng)、荒地裸地和林地,且所占面積變化均不大。牧場(chǎng)、荒地裸地和林地面積占比2000年分別為67.99%、23.69%、3.83%,2015年變?yōu)?7.84%、23.69%、3.83%。由表5可以看出,拉薩河流域2000—2015年,牧場(chǎng)變化面積較大,主要向荒地裸地和耕地轉(zhuǎn)移,面積分別為1 017 km2和196 km2;其次荒地裸地變化面積也較大,主要向牧場(chǎng)和水域轉(zhuǎn)移,面積分別為708 km2和41 km2。

    圖4 拉薩河流域2000年和2015年土地利用類型Fig.4 Land use map of Lhasa River Basin in 2000 and 2015

    表5 拉薩河流域2000—2015年土地利用類型變化轉(zhuǎn)移矩陣Table 5 Transfer matrix of land use change in Lhasa River Basin from 2000 to 2015 km2

    4.2.2 徑流對(duì)土地利用變化和氣候因素響應(yīng)分析

    本文利用率定后的SWAT模型分別計(jì)算了表2中各情景下的拉薩河流域徑流過程,并比較了流域總徑流深變化情況,結(jié)果如表6所示??梢钥闯?,在相同的土地利用情況下(情景S1與S2、情景S3與S4),2001—2016年氣候條件較1985—2000年的年均徑流深增加約22.91 mm;在相同氣候條件下(情景S1與S3、情景S2與S4),2015年土地利用狀況較2000年的年均徑流深增加約4.71 mm。氣候變化和土地利用對(duì)年均徑流影響占比分別約為82.9%和17.1%,土地利用變化對(duì)拉薩河流域徑流的影響較小,原因可能在于研究期1985—2016年內(nèi),拉薩河流域土地利用變化較小,而氣溫與降水顯著增加。土地利用狀況與氣候條件都發(fā)生變化時(shí)(情景S1與S4),年均徑流深平均增加約27.61 mm??傮w上看,目前拉薩河流域徑流隨土地利用與氣候條件變化呈同步增加趨勢(shì),且氣候條件變化是影響徑流的主要因素。

    表6 土地利用與氣候變化情景下的拉薩河流域徑流深變化統(tǒng)計(jì)Table 6 Statistics of runoff depth in Lhasa River basinunder land use and climate changes

    拉薩河各子流域表2情景下的徑流深變化情況如圖5所示。圖5(a)為各子流域在相同氣候條件下(情景S1與S3、情景S2與S4)2015年土地利用狀況較2000年的年徑流深變化情況,可以看出:拉薩河上游區(qū)域徑流深呈增加趨勢(shì);中游區(qū)域的支流呈增加趨勢(shì)而干流呈現(xiàn)減少趨勢(shì);下游區(qū)域支流呈減少趨勢(shì)而干流呈增加趨勢(shì)。圖5(b)為各子流域在相同土地利用狀況下(情景S1與S2、情景S3與S4)2001—2016年氣候條件較1985—2000年的年徑流深變化情況,可以看出:上中下游徑流深均呈增加趨勢(shì),下游支流區(qū)域較中上游增加趨勢(shì)更為明顯,原因可能在于相對(duì)中上游區(qū)域,下游區(qū)域2001—2016年的氣溫與降水較1985—2000年增加幅度更大。圖5(c)為各子流域年徑流深對(duì)土地利用與氣候變化的協(xié)同響應(yīng)結(jié)果,總體上徑流深均呈增加趨勢(shì),且下游支流區(qū)域的增加趨勢(shì)最為顯著,這也表明了目前拉薩河流域徑流變化的主導(dǎo)因素是氣候條件。

    圖5 拉薩河流域徑流深變化空間分布Fig.5 Spatial distribution of runoff depth changein Lhasa River basin

    4.3 徑流對(duì)氣候條件變化響應(yīng)分析

    降水(P1—P4)與氣溫(T1—T4)條件下的拉薩河流域徑流深變化結(jié)果如圖6所示。可以看出拉薩河流域徑流深變化與降水、氣溫呈較為顯著的正相關(guān),近似線性擬合下,降水每增加10%致徑流深增加約11.81%,氣溫每升高1 ℃致徑流深增加2.58%。

    圖6 拉薩河流域徑流深變化相關(guān)性分析Fig.6 Correlation analysis of runoff depth changes inLhasa River basin

    各組合情景下(PT1—PT8)的拉薩河流域徑流模擬計(jì)算結(jié)果也表明,流域徑流變化與降水、氣溫存在較為明顯的正相關(guān)關(guān)系,如表7所示。

    表7 降水與氣溫變化情景下的拉薩河流域徑流深變化統(tǒng)計(jì)Table 7 Statistics of runoff depth changes in Lhasa Riverbasin under precipitation and temperature changes

    拉薩河各子流域降水(P1—P4)與氣溫(T1—T4)條件下的徑流深變化結(jié)果如圖7所示。由圖7可以看出,各子流域徑流深均隨降水同步變化,且變化幅度差異較小,降水增大10%時(shí)各子流域徑流深增加約8%~12%。氣溫發(fā)生變化時(shí),中上游徑流深隨氣溫增高而減少,下游徑流深則隨氣溫增高而增加;氣溫每增加1℃,中上游區(qū)域徑流深約減小0.7%,下游區(qū)域徑流深約增加3.6%。

    圖7 降水與氣溫變化情景下的拉薩河流域徑流深變化空間分布Fig.7 Spatial distribution of runoff depth changes in Lhasa River basin under precipitation and temperature changes

    5 結(jié) 論

    本文在建立SWAT模型對(duì)拉薩河流域徑流進(jìn)行模擬的基礎(chǔ)上,通過比較不同情景模式下的徑流變化情況,量化分析土地利用、氣候變化對(duì)拉薩河流域徑流的影響,結(jié)果表明:

    (1)氣候變化與土地利用對(duì)徑流影響占比分別約為82.9%和17.1%,氣候變化對(duì)拉薩河流域徑流的影響較大,主要原因在于近30 a拉薩河流域土地利用情況變化不大,而氣溫和降水顯著增加。

    (2)徑流隨降水增大而增加,降水每增大10%,拉薩河各子流域徑流增加約8%~12%,總徑流增加約11.8%。

    (3)徑流隨氣溫變化的空間差異性較大,氣溫每升高1 ℃,拉薩河中上游徑流減小約0.7%,下游徑流增加約3.6%,但總徑流增加約2.5%。

    猜你喜歡
    拉薩河徑流氣候變化
    《應(yīng)對(duì)氣候變化報(bào)告(2022)》發(fā)布
    拉薩河流域浮游植物生態(tài)學(xué)研究現(xiàn)狀
    拉薩河,從古遠(yuǎn)流來
    布達(dá)拉(2018年4期)2018-05-14 08:53:30
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    拉薩市拉薩河城區(qū)段綜合整治工程4~#閘設(shè)計(jì)
    基于概率統(tǒng)計(jì)和模糊綜合評(píng)價(jià)法的水質(zhì)評(píng)價(jià)模型及其應(yīng)用——以拉薩河水質(zhì)評(píng)價(jià)為例
    Topmodel在布哈河流域徑流模擬中的應(yīng)用
    應(yīng)對(duì)氣候變化需要打通“網(wǎng)關(guān)”
    太陽能(2015年7期)2015-04-12 06:49:50
    探秘“大徑流”
    亚洲综合色惰| 成人毛片a级毛片在线播放| 久久久a久久爽久久v久久| 国产 一区 欧美 日韩| 精品99又大又爽又粗少妇毛片| 精品无人区乱码1区二区| 少妇人妻精品综合一区二区 | 亚洲第一区二区三区不卡| 欧美成人精品欧美一级黄| 插阴视频在线观看视频| 成人漫画全彩无遮挡| 国产亚洲精品久久久久久毛片| 亚洲av中文字字幕乱码综合| 真实男女啪啪啪动态图| 人妻久久中文字幕网| 色视频www国产| 春色校园在线视频观看| 日日干狠狠操夜夜爽| 国产亚洲91精品色在线| 草草在线视频免费看| 成人av一区二区三区在线看| 日韩国内少妇激情av| 成人欧美大片| 精品人妻一区二区三区麻豆 | 欧美另类亚洲清纯唯美| 99热这里只有精品一区| 99热这里只有是精品在线观看| 99在线视频只有这里精品首页| 日韩强制内射视频| 男女边吃奶边做爰视频| 男女边吃奶边做爰视频| 成人漫画全彩无遮挡| 综合色丁香网| 免费看美女性在线毛片视频| 欧美性感艳星| 日韩成人av中文字幕在线观看 | 中国美白少妇内射xxxbb| 亚洲精华国产精华液的使用体验 | 日韩成人伦理影院| 中文资源天堂在线| 可以在线观看毛片的网站| 99热6这里只有精品| 午夜免费激情av| a级毛色黄片| 欧洲精品卡2卡3卡4卡5卡区| 国产aⅴ精品一区二区三区波| 久久久久久久亚洲中文字幕| 赤兔流量卡办理| 久久久久国产精品人妻aⅴ院| 亚洲人成网站在线播放欧美日韩| 一区二区三区高清视频在线| 亚洲av中文av极速乱| 日韩一本色道免费dvd| 天天躁夜夜躁狠狠久久av| 国产高清激情床上av| 国产成年人精品一区二区| 精品久久久久久久久亚洲| 国产成年人精品一区二区| 欧美性猛交╳xxx乱大交人| 午夜精品在线福利| 菩萨蛮人人尽说江南好唐韦庄 | 欧美xxxx黑人xx丫x性爽| 国产精华一区二区三区| 一边摸一边抽搐一进一小说| 亚洲五月天丁香| 69av精品久久久久久| a级毛片a级免费在线| 嫩草影院精品99| 少妇裸体淫交视频免费看高清| 免费看美女性在线毛片视频| 校园人妻丝袜中文字幕| 毛片一级片免费看久久久久| 国模一区二区三区四区视频| videossex国产| 最近手机中文字幕大全| 简卡轻食公司| 一级av片app| 欧美高清性xxxxhd video| 国产又黄又爽又无遮挡在线| 免费观看精品视频网站| 看片在线看免费视频| 亚洲国产日韩欧美精品在线观看| 日韩国内少妇激情av| 在线观看66精品国产| 三级国产精品欧美在线观看| 日韩人妻高清精品专区| 国产精品福利在线免费观看| 少妇被粗大猛烈的视频| 亚洲av.av天堂| 国产精品三级大全| or卡值多少钱| 精品免费久久久久久久清纯| 一级毛片我不卡| 最近2019中文字幕mv第一页| 精品99又大又爽又粗少妇毛片| 国产高清三级在线| 成人精品一区二区免费| 亚洲丝袜综合中文字幕| 欧美人与善性xxx| 99久久九九国产精品国产免费| 国产亚洲精品av在线| 国产毛片a区久久久久| 亚洲国产精品合色在线| 亚洲成人av在线免费| 国内精品美女久久久久久| 日韩人妻高清精品专区| 精品久久久久久久久亚洲| 久久久久久久久久黄片| 国产欧美日韩精品亚洲av| 69av精品久久久久久| 精品人妻熟女av久视频| 亚洲人成网站在线播| 午夜免费男女啪啪视频观看 | 亚洲精品粉嫩美女一区| 日韩强制内射视频| 97热精品久久久久久| 欧美色欧美亚洲另类二区| aaaaa片日本免费| 久久久久精品国产欧美久久久| 国产色婷婷99| 国产精品99久久久久久久久| 日本黄色视频三级网站网址| 大型黄色视频在线免费观看| 99久久精品热视频| a级毛片a级免费在线| 一个人观看的视频www高清免费观看| 给我免费播放毛片高清在线观看| 一进一出抽搐动态| 校园春色视频在线观看| 国内精品久久久久精免费| 97超碰精品成人国产| 精品福利观看| 免费观看在线日韩| 人人妻人人澡欧美一区二区| 久久久久久久久大av| 女生性感内裤真人,穿戴方法视频| 乱人视频在线观看| 最近在线观看免费完整版| 国产精品久久久久久久久免| 亚洲不卡免费看| 中文字幕av在线有码专区| 国产精品1区2区在线观看.| 三级经典国产精品| 精品一区二区三区视频在线观看免费| 中国美女看黄片| 超碰av人人做人人爽久久| 日韩高清综合在线| 少妇人妻精品综合一区二区 | 一边摸一边抽搐一进一小说| 热99re8久久精品国产| 亚洲国产精品国产精品| 日韩一区二区视频免费看| 国产精品乱码一区二三区的特点| 国产精品久久久久久久电影| 日韩人妻高清精品专区| 国产不卡一卡二| 欧美精品国产亚洲| 美女xxoo啪啪120秒动态图| 卡戴珊不雅视频在线播放| 亚洲成人久久爱视频| 亚洲久久久久久中文字幕| 日韩成人伦理影院| 亚洲中文字幕日韩| 精品久久久久久久久亚洲| 国产一区二区亚洲精品在线观看| 免费搜索国产男女视频| 人妻少妇偷人精品九色| 在线免费十八禁| 精品久久久久久久久久久久久| 国产精品亚洲美女久久久| .国产精品久久| 亚洲精品亚洲一区二区| 亚洲18禁久久av| 一夜夜www| 亚洲成a人片在线一区二区| av黄色大香蕉| 黄色配什么色好看| 男人舔奶头视频| 久久中文看片网| 亚洲aⅴ乱码一区二区在线播放| 欧美色欧美亚洲另类二区| 国产在线精品亚洲第一网站| 极品教师在线视频| 丝袜美腿在线中文| 少妇丰满av| 精品熟女少妇av免费看| 欧美高清成人免费视频www| 久久韩国三级中文字幕| 51国产日韩欧美| 日韩中字成人| 日日摸夜夜添夜夜添av毛片| 国产伦精品一区二区三区四那| 黄色一级大片看看| 一进一出抽搐gif免费好疼| 五月伊人婷婷丁香| 少妇的逼好多水| 日韩av在线大香蕉| 麻豆av噜噜一区二区三区| 亚洲中文日韩欧美视频| av天堂中文字幕网| 亚洲自拍偷在线| 少妇的逼好多水| 国产精品久久视频播放| 日本黄色片子视频| 久久久a久久爽久久v久久| 国产午夜精品论理片| 国产女主播在线喷水免费视频网站 | 好男人在线观看高清免费视频| 韩国av在线不卡| 免费观看精品视频网站| 欧美不卡视频在线免费观看| 国产熟女欧美一区二区| 国产高清激情床上av| 免费高清视频大片| 欧美性猛交╳xxx乱大交人| 国产黄色视频一区二区在线观看 | 亚洲av免费在线观看| 日韩欧美免费精品| 午夜福利在线观看吧| 欧美中文日本在线观看视频| 精品国内亚洲2022精品成人| 精品福利观看| 国产一区二区三区av在线 | 亚洲国产欧美人成| 在线观看午夜福利视频| 老司机午夜福利在线观看视频| 国内精品久久久久精免费| 国产精品一区二区三区四区免费观看 | 99久久精品热视频| 嫩草影院新地址| 午夜爱爱视频在线播放| 成人美女网站在线观看视频| 亚洲色图av天堂| 91午夜精品亚洲一区二区三区| 十八禁网站免费在线| 欧美成人一区二区免费高清观看| 非洲黑人性xxxx精品又粗又长| 日本一本二区三区精品| 国内精品一区二区在线观看| 亚洲乱码一区二区免费版| 嫩草影院精品99| 亚洲欧美日韩高清专用| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣巨乳人妻| 一进一出抽搐gif免费好疼| 99热网站在线观看| 97在线视频观看| 男人和女人高潮做爰伦理| 成人亚洲精品av一区二区| 亚洲在线观看片| 国产午夜福利久久久久久| 国产亚洲91精品色在线| 春色校园在线视频观看| 欧美极品一区二区三区四区| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 少妇的逼好多水| 一本一本综合久久| 直男gayav资源| 白带黄色成豆腐渣| 少妇被粗大猛烈的视频| 激情 狠狠 欧美| aaaaa片日本免费| 午夜视频国产福利| 蜜臀久久99精品久久宅男| a级毛片免费高清观看在线播放| 男人的好看免费观看在线视频| 特大巨黑吊av在线直播| 搡老熟女国产l中国老女人| 亚洲精品一区av在线观看| 九色成人免费人妻av| 男人的好看免费观看在线视频| 成熟少妇高潮喷水视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品国产自在天天线| 校园人妻丝袜中文字幕| 热99re8久久精品国产| 99久久精品一区二区三区| 两个人的视频大全免费| 亚洲熟妇熟女久久| 久久精品国产亚洲网站| 看黄色毛片网站| 日本-黄色视频高清免费观看| 亚洲av成人av| 99在线视频只有这里精品首页| 亚洲在线观看片| 日本免费一区二区三区高清不卡| АⅤ资源中文在线天堂| 熟女人妻精品中文字幕| 成人精品一区二区免费| 国产真实伦视频高清在线观看| 在线观看美女被高潮喷水网站| 精品99又大又爽又粗少妇毛片| 亚洲av二区三区四区| eeuss影院久久| 少妇猛男粗大的猛烈进出视频 | 国产男人的电影天堂91| 国产欧美日韩一区二区精品| 国产在线男女| 亚洲精品日韩在线中文字幕 | 18禁在线无遮挡免费观看视频 | 成人国产麻豆网| 男人的好看免费观看在线视频| 亚洲av成人精品一区久久| 国内少妇人妻偷人精品xxx网站| av在线天堂中文字幕| 国产v大片淫在线免费观看| 欧美不卡视频在线免费观看| 91麻豆精品激情在线观看国产| 春色校园在线视频观看| 日日摸夜夜添夜夜添小说| 国产三级中文精品| 搡老妇女老女人老熟妇| 十八禁国产超污无遮挡网站| 变态另类成人亚洲欧美熟女| 日日撸夜夜添| 欧美色欧美亚洲另类二区| 特级一级黄色大片| 国产高清激情床上av| 国产亚洲av嫩草精品影院| 99热精品在线国产| 国产精品久久久久久亚洲av鲁大| 搡老妇女老女人老熟妇| 真实男女啪啪啪动态图| 日本与韩国留学比较| 久久6这里有精品| 99热这里只有是精品50| 又粗又爽又猛毛片免费看| 美女xxoo啪啪120秒动态图| 禁无遮挡网站| 亚洲人成网站在线观看播放| 亚洲av中文av极速乱| 欧美性猛交╳xxx乱大交人| 国产中年淑女户外野战色| 国产精品亚洲美女久久久| 久久精品国产自在天天线| 俺也久久电影网| 青春草视频在线免费观看| 一级黄色大片毛片| 日本免费一区二区三区高清不卡| 激情 狠狠 欧美| 久久草成人影院| 99视频精品全部免费 在线| 国产精品免费一区二区三区在线| 日本一二三区视频观看| 啦啦啦韩国在线观看视频| 91久久精品电影网| 国产精品爽爽va在线观看网站| 成人鲁丝片一二三区免费| 伦理电影大哥的女人| 久久久久久久久久黄片| 51国产日韩欧美| 久久鲁丝午夜福利片| 香蕉av资源在线| 在线免费十八禁| 日韩在线高清观看一区二区三区| 欧美一区二区精品小视频在线| 中国国产av一级| 色av中文字幕| 国产熟女欧美一区二区| 久久人人爽人人爽人人片va| 色视频www国产| 欧美激情久久久久久爽电影| 波野结衣二区三区在线| 亚洲内射少妇av| 国产成人aa在线观看| 日韩在线高清观看一区二区三区| 久久国产乱子免费精品| 亚洲va在线va天堂va国产| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 俺也久久电影网| 国产av麻豆久久久久久久| 精华霜和精华液先用哪个| 99久久精品国产国产毛片| 日韩在线高清观看一区二区三区| 久久亚洲精品不卡| 成人鲁丝片一二三区免费| 亚洲中文字幕一区二区三区有码在线看| 国国产精品蜜臀av免费| 国产v大片淫在线免费观看| 看十八女毛片水多多多| 精品日产1卡2卡| 在线a可以看的网站| 国产精品三级大全| 最新在线观看一区二区三区| 三级男女做爰猛烈吃奶摸视频| 哪里可以看免费的av片| 最新在线观看一区二区三区| 国产高清不卡午夜福利| 久久人人爽人人爽人人片va| 91精品国产九色| 国产极品精品免费视频能看的| 中文字幕熟女人妻在线| 色尼玛亚洲综合影院| 天天躁夜夜躁狠狠久久av| 两个人视频免费观看高清| 真实男女啪啪啪动态图| 91av网一区二区| 老司机影院成人| 日韩精品中文字幕看吧| 亚洲熟妇熟女久久| 成年女人毛片免费观看观看9| 99视频精品全部免费 在线| 一级毛片久久久久久久久女| 日韩欧美国产在线观看| 欧美又色又爽又黄视频| 一区二区三区免费毛片| 黄色欧美视频在线观看| 五月伊人婷婷丁香| 久久午夜亚洲精品久久| 搞女人的毛片| 长腿黑丝高跟| 精品久久久久久成人av| ponron亚洲| 久久久久国内视频| 国产精品人妻久久久久久| 日韩欧美一区二区三区在线观看| a级毛色黄片| 精品一区二区三区人妻视频| 亚洲av第一区精品v没综合| 又爽又黄a免费视频| 男人舔奶头视频| 俺也久久电影网| 一个人看视频在线观看www免费| 国产精品国产三级国产av玫瑰| 国产单亲对白刺激| 1000部很黄的大片| 亚洲va在线va天堂va国产| 久久鲁丝午夜福利片| 日本精品一区二区三区蜜桃| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 狠狠狠狠99中文字幕| 亚洲丝袜综合中文字幕| 一区福利在线观看| 午夜激情福利司机影院| 亚洲一区二区三区色噜噜| 少妇猛男粗大的猛烈进出视频 | 性插视频无遮挡在线免费观看| 成人美女网站在线观看视频| 国产一区二区在线观看日韩| 国产黄色视频一区二区在线观看 | 丰满乱子伦码专区| 搡女人真爽免费视频火全软件 | 淫秽高清视频在线观看| 成人亚洲精品av一区二区| 九九在线视频观看精品| 国产精品精品国产色婷婷| 亚洲图色成人| 午夜福利在线在线| 亚洲精品国产成人久久av| 午夜老司机福利剧场| 精品少妇黑人巨大在线播放 | 国产人妻一区二区三区在| 可以在线观看毛片的网站| 草草在线视频免费看| 精品久久久久久久末码| 97超碰精品成人国产| 一a级毛片在线观看| 国产精品国产三级国产av玫瑰| 久久久色成人| 又黄又爽又刺激的免费视频.| 免费在线观看影片大全网站| 国产精品,欧美在线| 精品人妻熟女av久视频| 欧美人与善性xxx| 天堂动漫精品| 三级毛片av免费| 高清午夜精品一区二区三区 | av在线蜜桃| 欧美一区二区国产精品久久精品| 老女人水多毛片| 我的老师免费观看完整版| 一个人观看的视频www高清免费观看| 蜜桃久久精品国产亚洲av| 又黄又爽又免费观看的视频| 国产精品三级大全| 国产视频内射| 精品少妇黑人巨大在线播放 | 中文字幕久久专区| 国产爱豆传媒在线观看| 久久天躁狠狠躁夜夜2o2o| 免费人成在线观看视频色| 国产不卡一卡二| 春色校园在线视频观看| a级一级毛片免费在线观看| 色在线成人网| 在线观看午夜福利视频| 老女人水多毛片| 热99在线观看视频| 国产欧美日韩一区二区精品| 网址你懂的国产日韩在线| av免费在线看不卡| 狂野欧美白嫩少妇大欣赏| 伊人久久精品亚洲午夜| 亚洲精品日韩在线中文字幕 | 丰满的人妻完整版| 九九热线精品视视频播放| 精品久久久久久久久久久久久| 大型黄色视频在线免费观看| 99久久精品国产国产毛片| 国产精品野战在线观看| 亚洲自偷自拍三级| 51国产日韩欧美| 日本黄大片高清| 少妇的逼好多水| 别揉我奶头 嗯啊视频| 91午夜精品亚洲一区二区三区| 俺也久久电影网| 亚洲av成人精品一区久久| 亚洲成人久久性| 午夜福利视频1000在线观看| 国产精品电影一区二区三区| 日本 av在线| 国产美女午夜福利| 九九在线视频观看精品| 黄色视频,在线免费观看| 丝袜喷水一区| 亚洲av中文av极速乱| 露出奶头的视频| 内地一区二区视频在线| 精品不卡国产一区二区三区| 高清日韩中文字幕在线| 给我免费播放毛片高清在线观看| 欧美日韩国产亚洲二区| 菩萨蛮人人尽说江南好唐韦庄 | 97超级碰碰碰精品色视频在线观看| 人人妻人人看人人澡| 国产黄a三级三级三级人| 成年女人毛片免费观看观看9| 亚洲精品一卡2卡三卡4卡5卡| 少妇裸体淫交视频免费看高清| 99热这里只有是精品50| 国产成人91sexporn| 成人性生交大片免费视频hd| .国产精品久久| 国产精品爽爽va在线观看网站| 美女被艹到高潮喷水动态| 美女免费视频网站| 精品久久国产蜜桃| 午夜视频国产福利| 五月玫瑰六月丁香| 啦啦啦观看免费观看视频高清| 男女下面进入的视频免费午夜| 特大巨黑吊av在线直播| 午夜福利高清视频| 18禁裸乳无遮挡免费网站照片| 国产精品亚洲一级av第二区| 久久精品影院6| 精品午夜福利视频在线观看一区| av中文乱码字幕在线| 韩国av在线不卡| 亚洲一区二区三区色噜噜| 精品一区二区三区视频在线| 久久久久国产网址| 卡戴珊不雅视频在线播放| 国产精品野战在线观看| 有码 亚洲区| 看黄色毛片网站| 亚洲精品日韩在线中文字幕 | 人人妻,人人澡人人爽秒播| 国产黄色视频一区二区在线观看 | 亚洲国产欧洲综合997久久,| 热99re8久久精品国产| 大又大粗又爽又黄少妇毛片口| 日本与韩国留学比较| 我要搜黄色片| 观看免费一级毛片| 国产三级在线视频| 日本欧美国产在线视频| 麻豆久久精品国产亚洲av| 天美传媒精品一区二区| 伊人久久精品亚洲午夜| 亚洲精品色激情综合| 国产探花在线观看一区二区| 免费av毛片视频| 又黄又爽又刺激的免费视频.| 精品国产三级普通话版| 精品人妻一区二区三区麻豆 | 亚洲四区av| 黄色欧美视频在线观看| 美女大奶头视频| 国产淫片久久久久久久久| 免费看a级黄色片| 中出人妻视频一区二区| 少妇高潮的动态图| 卡戴珊不雅视频在线播放| 亚洲在线观看片| 欧美一区二区国产精品久久精品| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 亚洲欧美精品综合久久99| 五月伊人婷婷丁香| 精品人妻熟女av久视频| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放| 日本成人三级电影网站| 日韩欧美免费精品| 国产毛片a区久久久久| 黄片wwwwww| 别揉我奶头 嗯啊视频| 亚洲熟妇中文字幕五十中出| 国产精品永久免费网站| 美女大奶头视频| 日本欧美国产在线视频| 婷婷精品国产亚洲av在线| 在线看三级毛片| 国产探花在线观看一区二区| 日韩欧美一区二区三区在线观看| av在线老鸭窝| 九九爱精品视频在线观看| 99热只有精品国产| 亚洲国产精品国产精品| 九色成人免费人妻av|