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

    基于SDSM的撫河流域未來極端氣溫模擬與預估

    2019-08-27 03:02:54朱圣男劉衛(wèi)林萬一帆吳德勝黃一鵬劉麗娜
    人民珠江 2019年8期
    關鍵詞:最低氣溫氣溫尺度

    朱圣男,劉衛(wèi)林,萬一帆,吳德勝,黃一鵬,劉麗娜

    (1.南昌工程學院 江西省水文水資源與水環(huán)境重點實驗室,江西 南昌 330099;2.南昌工程學院 瑤湖學院,江西 南昌 330099)

    氣候變化和人類活動對流域氣候要素的影響已成為熱點問題之一,尤其是對區(qū)域性水文系統(tǒng)的影響更是近年來各國發(fā)展中面臨的重大問題。撫河流域農業(yè)發(fā)達,流域所處的贛撫平原是江西省最大的糧食生產地,水資源開發(fā)利用強度大,意味著該流域水資源安全將達到新的高度。隨著氣候變化、土地利用和人類活動的影響,撫河流域內水資源問題日益嚴重,對區(qū)域內居民生活起居、農業(yè)生產需要、社會經濟發(fā)展乃至整個生態(tài)系統(tǒng)產生深遠影響和一系列社會問題。與此同時,氣候變化導致極端天氣事件頻發(fā),尤其體現(xiàn)在對氣候變化較為敏感和脆弱的地區(qū)[1]。因此,開展撫河流域氣候變化特征研究,對本區(qū)域內的水資源利用、水利工程安全高效運行、保證社會經濟穩(wěn)定高速發(fā)展具有重要意義。

    目前,全球氣候模式(Global Climate Model,GCM)和降尺度方法是研究未來氣候變化的主要工具,其輸出信息只能反映100 km2或以上的大尺度網格氣候變化的平均特征[2],直接使用將會忽略植被、地形等小尺度的差別,輸出結果的空間分辨率較低,對研究精度有影響[3]。因此,使用GCM預測結果作為區(qū)域水文水資源變化評估輸入因子時,需進行降尺度分析,將GCM輸出的大尺度、低分辨率信息轉化為區(qū)域尺度信息。降尺度分為動力降尺度、統(tǒng)計降尺度、動力降尺度與統(tǒng)計降尺度相結合的方法[3]。統(tǒng)計降尺度始于20世紀90年代,它通過使用多元的線性回歸方法推導出預報因子和預報量之間的統(tǒng)計關系,將大尺度氣候模式轉化為可應用的小尺度,提高了區(qū)域分辨率,方法多,計算量小,在國內外得到了廣泛應用。任何與預報量相關的大尺度大氣變量都可以作為預報因子。本文使用的第5次國際耦合模式比較計劃CMIP5比前一階段CMIP3更加精確與復雜[4],是目前較先進的預估未來氣候變化方法[5],該模式對中國區(qū)域模擬,模擬氣溫較降水偏強[6],且在鄱陽湖流域應用較少。

    本文建立撫河流域氣象要素的降尺度模型,選擇CMIP5中的3種情景模式RCP2.6、RCP4.5、RCP8.5,選取A(2006—2035年)、B(2036—2065年)、C(2066—2100年) 3個時段,以1991—2005年為基準期,基于氣象資料對撫河流域RCP2.6、RCP4.5、RCP8.5 3種排放情景下2006—2100年最高氣溫和最低氣溫的可能變化情況進行了模擬和對比分析,為撫河流域氣候變化下的水文響應研究和氣候變化下適應政策的制定提供科學依據(jù)。

    1 研究區(qū)概況及研究方法

    1.1 研究區(qū)域

    圖1 撫河流域

    撫河流域(圖1)位于江西省東部,發(fā)源于武夷山脈西麓廣場新驛前鄉(xiāng)血木嶺[7]。地處東經115°36至117°10′、北緯26°30′至28°50′之間,是鄱陽湖水系的五大河流之一,主河道長348 km,流域面積16 493 km2。南城以上段包括南豐和廣昌為盱江,為撫河上游;南城至廖家灣為中游;撫州以下為撫河下游。流域地形東南高,西北低,南北寬,東西狹。撫河流域屬于中亞熱帶濕潤季風氣候區(qū),氣候溫和、雨量充沛,多年平均年水面蒸發(fā)量894 mm;多年平均年降水量1 732 mm,降水時空分配不均,東南部降水量多于西北部;多年平均氣溫17.8℃,極端最高氣溫42.1℃(1971年)、最低氣溫-12.7℃(1991年)。

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

    本文采用的數(shù)據(jù)主要包括實測氣象站點(南城、廣昌)數(shù)據(jù)、NCEP再分析資料數(shù)據(jù)及GCM(CanESM2)的輸出數(shù)據(jù)。

    a) 實測氣象站點數(shù)據(jù)。研究中用到的實測氣象數(shù)據(jù)來自中國氣象科學數(shù)據(jù)共享服務網(http://cdc.cma.gov.cn/home.do),根據(jù)氣象站的地理位置及氣象資料的完整性,選取了撫河流域2個氣象站1961—2005 年共45 a的氣象資料作為研究對象,經過數(shù)據(jù)處理以后,選取氣溫數(shù)據(jù)序列作為模型的預報量,主要包括日最高氣溫(Tmax)和日最低氣溫(Tmin)。

    b) NCEP再分析數(shù)據(jù)。數(shù)據(jù)為1961—2005年的日數(shù)據(jù),一共有26個因子,主要包括平均海平面氣壓、地表平均氣溫、比濕、相對濕度、地轉風速、風向、緯向風速、經向風速、渦度、散度,500 hPa位勢高度及500 hPa高度場的相對濕度、地轉風速、風向、緯向風速、經向風速、渦度、散度,850 hPa位勢高度及850 hPa高度場的相對濕度、地轉風速、風向、緯向風速、經向風速、渦度、散度。NCEP 網格大小為1.875°×1.875°,與CanESM2不符,因此將NCEP網格數(shù)據(jù)轉換成與CanESM2數(shù)據(jù)格式統(tǒng)一。

    c) GCM數(shù)據(jù)。GCM輸出數(shù)據(jù)采用加拿大環(huán)境與氣候變化網(https://www.canada.ca/en.html)CanESM2數(shù)據(jù),共有128×64個網格,選取6個網格為本研究區(qū)域覆蓋的網格,網格分辨率為2.8125°×2.7906°。選擇3種排放情景RCP2.6(低等溫室氣體排放情景)、RCP4.5(中等溫室氣體排放情景)和RCP8.5(最高溫室氣體排放情景)典型濃度路徑,選擇的預報因子與NCEP相同,對未來情景下最高氣溫和最低氣溫進行模擬。

    2 SDSM統(tǒng)計降尺度模型

    SDSM模型是由英國Wilby等人建立的降尺度工具,融合了天氣發(fā)生器和多元線性回歸技術,是一種轉換函數(shù)與隨機天氣發(fā)生器耦合的降尺度方法,廣泛應用于亞洲、歐洲、美洲的氣象、水文等研究領域。使用SDSM統(tǒng)計降尺度模型主要包括以下6個步驟:①選擇合適的預報量和預報因子;②數(shù)據(jù)檢驗與轉換;③模型參數(shù)率定;④天氣發(fā)生器;⑤生成未來情景;⑥分析數(shù)據(jù)結果。

    SDSM統(tǒng)計降尺度模型的核心思想是確定氣象站點尺度下的預報量(如日最高氣溫、最低氣溫)與大尺度下的預報因子(如平均海平面氣壓mslp、地表平均氣溫temp和比濕shum等)之間的數(shù)學統(tǒng)計關系[8]。SDSM模型中預報因子與預報量間統(tǒng)計關系的基本原理如下:

    (1)

    式中Ti——氣溫變量;Pij——預報因子;n——預報因子數(shù)量;γ——模式參數(shù);ei——模型誤差。

    2.1 預報因子的選擇

    預報因子的選取原則主要是:①所選的預報因子與預報量之間要有強相關性;②所選因子必須能夠表征大尺度氣候場的重要物理過程;③GCM必須能夠準確模擬出所選的預報因子;④所選因子之間應為弱相關[9]。預報因子選擇影響著模擬程度的好壞,選擇相關性更強的預報因子會預報更準確的信息。為了防止出現(xiàn)多個預報因子共線現(xiàn)象,對地面氣溫進行估計時,排除互相關系大的情況[10]。本文主要研究共26個預報因子,對2個站點利用逐步回歸的方法,根據(jù)這4個原則,選取其中最優(yōu)預報因子,見表1。

    本研究預報因子選擇控制在5種,廣昌、南城2個站點所選擇預報因子對于3種情景下氣溫模式大部分相同,表明這些預報因子對氣溫模擬敏感性較高,其中在temp(2 m平均氣溫)的選擇率最高,此外shum、p1_v在2個站點都有被選用。這些因子一方面通過影響大氣環(huán)流的經向和緯向分布,從而影響氣候的經緯地帶性;另一方面海平面氣壓、位勢高度場和地表溫度等氣象環(huán)流因子通過熱量和水分的擴散及輸送過程,進而影響對氣溫產生的變化[11]。

    表1 各站點預報因子選擇

    注:temp表示在2 m平均氣溫;s表示近表面比濕度;p表示海平面氣壓;_v表示經向速度分量;_z表示渦度;_u表示緯向速度分量;p500表示500 hPa位勢高度場;s500表示在500 hPa高度比濕度

    2.2 模型的驗證與率定

    為了評價模型模擬的效果,選取了擬合度R2和標準差相對誤差SE作為評價指標。流域SDSM模型率定結果見表2、3。從表2可知,兩站點氣溫的擬合度均在65%以上,氣溫的標準誤差在1.02~3.72之間。而根據(jù)有關學者研究結果指出擬合度一般高于60%,說明本次氣溫率定結果很好[12]。

    表2 撫河流域SDSM模型率定期擬合度和標準誤差(月) %

    表3 撫河流域SDSM模型率定期擬合度和標準誤差(年) %

    使用獨立檢驗的方法,將觀測序列分為2個時間段:第一時間段為1961—1990年,作為率定期,用于建立統(tǒng)計關系;第二時間段為1991—2005年,作為驗證期,用于驗證模型的可靠性。通過比較,分析兩階段的數(shù)據(jù)信息驗證結果。Hanssen-Bauer 和 Forland[13]兩位學者曾經指出在經驗模型的參數(shù)率定期序列長度應超過5 a,因此,本文應用30、15 a的數(shù)據(jù)來率定和驗證模型是符合要求的[14]。

    圖2為撫河流流域最高氣溫和最低氣溫率定期、驗證期與實測值的對比。根據(jù)圖2可得最高氣溫和最低氣溫的模擬效果較好。率定期與實測值(圖2a、2c)相比較,呈現(xiàn)出震蕩趨勢多年平均相差0.08℃~0.23℃,結果可以接受。在預報量間的對比中可以得出:最低氣溫在2—4、10—11月模擬效果相對其他月份較差,平均差值在0.4℃左右,其中3月最為明顯,在1℃左右。最高氣溫綜合來看擬合效果最佳,僅在率定期11月誤差稍大,其他月份擬合效果良好,擬合度良好月份最多。該模型率定期和驗證期氣溫擬合整體較好,符合預期要求。

    3 未來極端氣溫變化

    3.1 未來最高氣溫與最低氣溫變化

    撫河流域最高、最低氣溫多年平均實測值與預測期2006—2100年的RCP2.6、RCP4.5、RCP8.5模式的模擬值對比見圖3a、3b。折線圖表現(xiàn)的趨勢基本可以完全反映出未來的氣溫變化。對比撫河流域最高氣溫與基準期可知,在RCP2.6、RCP4.5、RCP8.5情景下,未來氣溫相對基準期都將呈現(xiàn)上升趨勢,這與IPCC報告指出的全球氣溫將持續(xù)上升相吻合[15]。各氣溫、各情景下與實測值變化趨勢相近,也進一步表明了本模型是可靠的[16]。

    將2006—2100年的預測數(shù)據(jù)分為3個階段:A(2006—2035年)、B(2036—2065年)、C(2066—2100年)。分別將最高、最低氣溫在3個時段、3種排放情景下比較分析,從圖3可知,在3種情景下,最低氣溫與最高氣溫在1—7月之間增溫相對明顯,約2℃左右;在高排放情景RCP8.5下全年增溫幅度均較大,約2℃~5℃左右,而RCP2.6和RCP4.5排放情景下在1—7月增幅較大,約2℃~4℃,其他月份增幅微??;RCP2.6和RCP8.5情景下在8—12月最低氣溫比多年平均氣溫還低0.5℃~1℃,其他月份增幅微小。

    a) 率定期、實測期最高氣溫

    b) 驗證期、實測期最高氣溫

    c) 率定期、實測期最低氣溫

    d) 驗證期、實測期最低氣溫

    圖2撫河流域率定期、驗證期與實測期氣溫對比

    a) 最高氣溫多年平均實測值與相關模擬值對比

    b) 最低氣溫多年平均實測值與相關模擬值對比

    圖3撫河流域未來情景最低和最高氣溫預測

    3.2 最高氣溫分析

    由圖4可以看出:同基準期一樣,撫河流域未來月最高氣溫的年內變化趨勢基本一致,整體呈先遞增后遞減的變化趨勢,3種情景下最高溫度均出現(xiàn)在7月份,A階段中7—8、12月出現(xiàn)了降溫趨勢,最高下降0.8℃;3—6月增溫在1.6℃~2.9℃。B階段,7—8、11—12月降溫在0.3℃~0.8℃之間;其他月份為升溫趨勢,其中2—6月增溫在2.1℃~4.5℃。C階段,7—8、10—12月降溫幅度在0.1℃~1.7℃;其他月份表現(xiàn)為升溫,最高5.5℃(3月),最低2.8℃(6月),升高幅度大于前兩階段。

    a) RCP2.6圖4 撫河流域情景最高氣溫與基準期相對比

    b) RCP4.5

    c) RCP8.5

    續(xù)圖4撫河流域情景最高氣溫與基準期相對比

    3.3 最低氣溫分析

    由圖5可以看出:3種情景下3個時段未來月最低氣溫的年內變化與最高氣溫變化趨勢一致。在RCP2.6情景下,A階段全年氣溫平均升高約1.1℃,9—10月降溫幅度約1.45℃;1—6、12月升溫在1.5℃~3.3℃,其他月份均低于1℃。B階段,大部分月份呈現(xiàn)升溫趨勢,9—11月表現(xiàn)為降溫,1—6月升溫幅度在1.3℃~4.9℃。RCP4.5情景:C階段全年氣溫平均增加約2.8℃,各月份增幅2℃以上,最高達4.4℃(3月)。B階段幅度較A大,最高增溫6.9℃(3月),且2—6月均增溫5℃以上。C階段全年氣溫平均增加約4.3℃。RCP8.5情景下:A階段全年氣溫平均增加約1.1℃,1—6月增溫幅度較大在1.3℃以上,最高增加2.6℃(3月)。B階段,1—6月份增溫在3℃,最高達4.8℃。C階段全年氣溫平均增加約3.4℃,但在9—11月表現(xiàn)為降溫,1—6月份增溫約4℃,其他月份增溫在2℃以上。

    a) RCP2.6

    b) RCP4.5

    c) RCP8.5

    圖5撫河流域情景最低氣溫與基準期相對比

    利用SDSM模型模擬的最高氣溫與最低氣溫在RCP2.6、RCP4.5和RCP8.5情景下相對于基準期的增量變化見表4。3種情景下的最高氣溫與最低氣溫均呈較為明顯的上升趨勢,最高氣溫在A、B和C 3個時期平均分別升高0.97℃、1.78℃、2.31℃,最低氣溫在3個時期平均分別升高1.70℃、2.40℃、3.12℃,最低氣溫在C上升幅度最大。其中RCP4.5情景下最低氣溫上升4.34℃。3種情景下,撫河流域未來最低氣溫在秋季(9—11月)均低于基準期,其他情景時段下均高于基準期。由此可推測,撫河流域未來存在發(fā)生極端高溫天氣的風險。

    表4 撫河流域未來極端氣溫變化 ℃

    4 結語

    本文以撫河流域為例,利用SDSM統(tǒng)計降尺度模型對南城和廣昌2個站點日最高氣溫和最低氣溫在RCP2.6、RCP4.5、RCP8.5 3種排放情景下流域未來最高及最低氣溫可能變化的情況進行了模擬和對比分析。運用CMIP5模式、降尺度方法,取1961—1990年為率定期,1991—2005年為驗證期,在2個站點篩選出不同的預報因子,模擬撫河流域2006—2100年3個時段A、B、C間RCP2.6、RCP4.5、RCP8.5 3種不同情景下的最高氣溫和最低氣溫變化,經分析后有如下結論。

    a) 撫河流域未來氣溫將持續(xù)升高,并且最低氣溫升高的幅度大于最高氣溫,這一趨勢與全國氣溫變化趨勢總體上保持一致。

    b) 對于撫河流域未來氣溫變化的增量,最高氣溫變化平均增量在3種情景模式下分別為1.43℃、1.42℃、2.21℃;最低氣溫變化增量分別為1.36℃、3.68℃、2.28℃。這說明隨著社會經濟發(fā)展,極端氣候要素的事件也隨之出現(xiàn)并概率增大。

    c) 3種情景下的最高氣溫與最低氣溫均呈較為明顯的上升趨勢,最高氣溫在A、B和C 3個時期平均分別升高0.97℃、1.78℃、2.31℃;最低氣溫在3個時期平均分別升高1.71℃、2.40℃、3.12℃,最低氣溫在C上升幅度最大。

    d) 撫河流域未來氣溫隨著時間段的推移和排放量的增大增溫幅度越來越大,而隨著排放量的增大,升溫月份出現(xiàn)也更加頻繁。最高氣溫下僅有RCP2.6情景下11—12月氣溫低于多年平均氣溫,最低氣溫下RCP2.6情景和RCP8.5情景下9—11月氣溫低于多年平均氣溫,其他情景下均出現(xiàn)升溫趨勢。表明高溫天氣很有可能將提前到來,未來發(fā)生極端高溫天氣的頻率越來越大。

    本文僅使用一種氣候模式數(shù)據(jù)進行氣候變化預估,具有不確定性。在以后的研究中,將增加更多的評價指標、采用不同的降尺度方法、運用多模式集合等方法綜合分析未來氣候變化,以減少模擬結果的不確定性。

    猜你喜歡
    最低氣溫氣溫尺度
    氣溫的神奇“魔法”
    基于FY-3D和FY-4A的氣溫時空融合
    財產的五大尺度和五重應對
    北辰地區(qū)日最低氣溫特征及影響因素研究
    天津科技(2019年3期)2019-03-30 07:17:14
    67年來濟南最高和最低氣溫變化特征
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    數(shù)學選擇題練習
    與氣溫成反比的東西
    9
    室外雕塑的尺度
    雕塑(2000年2期)2000-06-22 16:13:30
    欧美日韩亚洲综合一区二区三区_| 久久精品夜夜夜夜夜久久蜜豆 | √禁漫天堂资源中文www| 两人在一起打扑克的视频| 久久久久性生活片| 久久久久精品国产欧美久久久| 久久热在线av| 亚洲av中文字字幕乱码综合| 麻豆一二三区av精品| 久久久久免费精品人妻一区二区| 国产精品亚洲一级av第二区| 色哟哟哟哟哟哟| 午夜免费成人在线视频| 国产精品,欧美在线| 久久久国产精品麻豆| 精品国产乱码久久久久久男人| 51午夜福利影视在线观看| 十八禁网站免费在线| 黄色视频不卡| 男人舔奶头视频| 亚洲精品久久成人aⅴ小说| 欧美黄色片欧美黄色片| 一个人免费在线观看电影 | 亚洲精品中文字幕一二三四区| 色精品久久人妻99蜜桃| 禁无遮挡网站| 欧美乱妇无乱码| 美女 人体艺术 gogo| 熟女少妇亚洲综合色aaa.| 色av中文字幕| 欧美一级a爱片免费观看看 | 日日摸夜夜添夜夜添小说| av福利片在线观看| 国产三级中文精品| 日韩欧美三级三区| 免费人成视频x8x8入口观看| bbb黄色大片| 久久精品国产99精品国产亚洲性色| 国产av一区在线观看免费| 可以在线观看毛片的网站| 又黄又爽又免费观看的视频| 欧美又色又爽又黄视频| 亚洲18禁久久av| 一个人免费在线观看电影 | 欧美黑人欧美精品刺激| 天天添夜夜摸| 国产精华一区二区三区| 国产亚洲欧美98| 正在播放国产对白刺激| 啦啦啦观看免费观看视频高清| 91九色精品人成在线观看| 久热爱精品视频在线9| 日韩精品青青久久久久久| 国产三级在线视频| 国产私拍福利视频在线观看| 99国产精品一区二区蜜桃av| e午夜精品久久久久久久| 亚洲av成人一区二区三| 国产三级中文精品| 久久久久久国产a免费观看| 1024视频免费在线观看| 亚洲avbb在线观看| 日韩国内少妇激情av| 欧美三级亚洲精品| 最好的美女福利视频网| 欧美午夜高清在线| 90打野战视频偷拍视频| 免费搜索国产男女视频| 国产欧美日韩一区二区精品| 国产精品美女特级片免费视频播放器 | 欧美黑人欧美精品刺激| 中文字幕熟女人妻在线| 十八禁人妻一区二区| 成人av在线播放网站| 国产69精品久久久久777片 | 亚洲天堂国产精品一区在线| 欧美中文综合在线视频| 精品国产乱子伦一区二区三区| 久久热在线av| 女人爽到高潮嗷嗷叫在线视频| 不卡av一区二区三区| 亚洲电影在线观看av| 波多野结衣高清无吗| 丰满的人妻完整版| 黄色视频,在线免费观看| 免费搜索国产男女视频| e午夜精品久久久久久久| 欧美久久黑人一区二区| 特级一级黄色大片| 国产一级毛片七仙女欲春2| 色在线成人网| 久久精品影院6| x7x7x7水蜜桃| 91国产中文字幕| 露出奶头的视频| 欧美一区二区国产精品久久精品 | 国产欧美日韩一区二区三| 天天躁狠狠躁夜夜躁狠狠躁| avwww免费| 91麻豆精品激情在线观看国产| 久久99热这里只有精品18| 免费观看a级毛片全部| 男人舔女人下体高潮全视频| 天堂√8在线中文| 一个人观看的视频www高清免费观看| 亚洲av免费高清在线观看| 成人亚洲精品av一区二区| 免费一级毛片在线播放高清视频| 变态另类成人亚洲欧美熟女| 亚洲欧美精品专区久久| 国产亚洲91精品色在线| 女人十人毛片免费观看3o分钟| 国产美女午夜福利| 亚洲国产精品国产精品| 日日摸夜夜添夜夜添av毛片| 免费观看在线日韩| 久久国产乱子免费精品| 九色成人免费人妻av| 国产精品久久久久久av不卡| 51国产日韩欧美| 中文字幕人妻熟人妻熟丝袜美| 国产伦理片在线播放av一区 | 亚洲一区二区三区色噜噜| 亚洲综合色惰| 91久久精品国产一区二区三区| 男插女下体视频免费在线播放| 国产一级毛片在线| 国产伦精品一区二区三区四那| 欧美最新免费一区二区三区| 熟女电影av网| 亚洲精品日韩av片在线观看| 美女国产视频在线观看| 久久久久久久午夜电影| 91久久精品国产一区二区成人| 在线免费观看不下载黄p国产| 天堂网av新在线| 亚洲最大成人中文| 久久久久久九九精品二区国产| 国产成人91sexporn| 日日撸夜夜添| 日韩成人av中文字幕在线观看| 熟女人妻精品中文字幕| 又爽又黄a免费视频| 日本熟妇午夜| 亚洲精品亚洲一区二区| 一边亲一边摸免费视频| 久久午夜福利片| 久久99热6这里只有精品| 国产成人精品一,二区 | 嫩草影院精品99| 免费人成视频x8x8入口观看| 婷婷色av中文字幕| 成人美女网站在线观看视频| 亚洲av二区三区四区| 久久久午夜欧美精品| 亚洲欧美精品专区久久| 国产一区二区亚洲精品在线观看| 国产黄色视频一区二区在线观看 | 中文亚洲av片在线观看爽| eeuss影院久久| 99九九线精品视频在线观看视频| 蜜桃久久精品国产亚洲av| 中文亚洲av片在线观看爽| 久久久久久国产a免费观看| 国产精品一二三区在线看| 婷婷精品国产亚洲av| 哪里可以看免费的av片| 性欧美人与动物交配| 亚洲经典国产精华液单| 国产成人精品婷婷| 桃色一区二区三区在线观看| 国产精品一区二区三区四区免费观看| 国产精品一及| 最近2019中文字幕mv第一页| 黄色配什么色好看| av免费在线看不卡| 午夜激情福利司机影院| 黄色一级大片看看| 亚洲美女视频黄频| 免费观看的影片在线观看| kizo精华| 国内揄拍国产精品人妻在线| 午夜精品在线福利| kizo精华| 国产蜜桃级精品一区二区三区| 午夜爱爱视频在线播放| 亚洲在久久综合| 亚洲综合色惰| 99热6这里只有精品| 久久久久国产网址| 久久精品夜夜夜夜夜久久蜜豆| 男人的好看免费观看在线视频| 99精品在免费线老司机午夜| 欧美最新免费一区二区三区| ponron亚洲| 午夜福利视频1000在线观看| 亚洲久久久久久中文字幕| 亚洲成人av在线免费| 免费av毛片视频| 国产人妻一区二区三区在| 国产精品久久电影中文字幕| 国产成人freesex在线| ponron亚洲| 国产精品.久久久| 国内精品美女久久久久久| 国产精品国产三级国产av玫瑰| 国产亚洲av嫩草精品影院| 欧美3d第一页| 看黄色毛片网站| 中文字幕av成人在线电影| 在线免费观看的www视频| 一级毛片我不卡| 国产乱人偷精品视频| 一个人看的www免费观看视频| 热99在线观看视频| 精品久久久久久久末码| 青春草国产在线视频 | 国产精品乱码一区二三区的特点| 久久欧美精品欧美久久欧美| 噜噜噜噜噜久久久久久91| 人妻少妇偷人精品九色| 亚洲在线自拍视频| 国产黄片美女视频| 日日撸夜夜添| 久久午夜亚洲精品久久| 成人美女网站在线观看视频| 久久久久久九九精品二区国产| 我要搜黄色片| 欧美成人精品欧美一级黄| 日韩在线高清观看一区二区三区| 亚洲欧美日韩无卡精品| 久久韩国三级中文字幕| 国产在视频线在精品| 一进一出抽搐gif免费好疼| 久久久久久国产a免费观看| a级毛色黄片| 黄色一级大片看看| 美女黄网站色视频| 一区二区三区免费毛片| 欧美3d第一页| 国产成人一区二区在线| 人妻制服诱惑在线中文字幕| 国产高清激情床上av| 精品久久久久久久末码| 国产高清三级在线| 身体一侧抽搐| 狂野欧美激情性xxxx在线观看| 夜夜夜夜夜久久久久| www.av在线官网国产| 日韩欧美精品v在线| 亚洲婷婷狠狠爱综合网| 成年免费大片在线观看| 三级经典国产精品| 国产伦精品一区二区三区视频9| 男插女下体视频免费在线播放| 午夜精品一区二区三区免费看| 日韩av不卡免费在线播放| 亚洲国产欧美人成| 可以在线观看的亚洲视频| 中文在线观看免费www的网站| 免费av不卡在线播放| 五月玫瑰六月丁香| 亚洲av第一区精品v没综合| 蜜桃亚洲精品一区二区三区| 亚洲国产精品合色在线| 此物有八面人人有两片| 久久6这里有精品| 久久久久久久久久久丰满| 人妻夜夜爽99麻豆av| 欧美性感艳星| 日本撒尿小便嘘嘘汇集6| 亚洲国产日韩欧美精品在线观看| 亚洲成人中文字幕在线播放| 男插女下体视频免费在线播放| 午夜视频国产福利| 国产一区二区三区av在线 | 亚洲人成网站在线播| www.av在线官网国产| 国产伦精品一区二区三区视频9| 国产成人a∨麻豆精品| 亚洲国产精品久久男人天堂| 欧美精品一区二区大全| 日韩成人av中文字幕在线观看| 成人亚洲精品av一区二区| 亚洲色图av天堂| 在线观看一区二区三区| 69人妻影院| 久久久国产成人精品二区| 免费一级毛片在线播放高清视频| 亚洲国产精品sss在线观看| 青春草亚洲视频在线观看| 欧美丝袜亚洲另类| 亚洲成av人片在线播放无| 蜜臀久久99精品久久宅男| 日日摸夜夜添夜夜爱| 老司机福利观看| 少妇熟女aⅴ在线视频| 在线观看av片永久免费下载| 久久这里只有精品中国| 少妇猛男粗大的猛烈进出视频 | 91精品一卡2卡3卡4卡| 精品人妻视频免费看| 国产高清有码在线观看视频| 男人狂女人下面高潮的视频| 天天一区二区日本电影三级| 久久99热这里只有精品18| 好男人视频免费观看在线| 黄色欧美视频在线观看| 精品人妻偷拍中文字幕| 久久久久久大精品| 成熟少妇高潮喷水视频| 国产午夜福利久久久久久| 亚洲国产欧洲综合997久久,| 蜜桃亚洲精品一区二区三区| 又粗又爽又猛毛片免费看| 亚洲国产精品成人综合色| 日日摸夜夜添夜夜添av毛片| 99在线人妻在线中文字幕| 日本三级黄在线观看| av国产免费在线观看| 在线免费观看不下载黄p国产| 精品一区二区三区人妻视频| 中国美女看黄片| 国产成人freesex在线| 美女xxoo啪啪120秒动态图| 看片在线看免费视频| 国产高清视频在线观看网站| 色尼玛亚洲综合影院| av免费观看日本| 麻豆av噜噜一区二区三区| 亚洲精品亚洲一区二区| 国产精品一区www在线观看| 欧美日韩综合久久久久久| 男女视频在线观看网站免费| 免费黄网站久久成人精品| 一区二区三区高清视频在线| 嫩草影院精品99| 女人被狂操c到高潮| 看十八女毛片水多多多| 国产精品免费一区二区三区在线| 国产精品久久久久久精品电影小说 | 精品国产三级普通话版| 老司机影院成人| 麻豆成人av视频| 亚洲一级一片aⅴ在线观看| 少妇人妻精品综合一区二区 | 亚洲国产欧洲综合997久久,| 少妇熟女欧美另类| 看十八女毛片水多多多| 国产精品三级大全| 久久精品国产亚洲网站| 少妇熟女aⅴ在线视频| 色综合亚洲欧美另类图片| 一边摸一边抽搐一进一小说| 中文字幕av在线有码专区| 亚洲欧洲国产日韩| 亚洲国产精品成人久久小说 | 一级毛片电影观看 | 国产视频内射| 美女被艹到高潮喷水动态| 国产精品福利在线免费观看| 国产精品久久视频播放| 男女那种视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产av麻豆久久久久久久| 国内精品美女久久久久久| 色噜噜av男人的天堂激情| 夜夜爽天天搞| 精品99又大又爽又粗少妇毛片| av在线亚洲专区| 日韩欧美一区二区三区在线观看| 久99久视频精品免费| 国产午夜精品论理片| 日日撸夜夜添| 成年免费大片在线观看| 插逼视频在线观看| 亚洲最大成人av| 精品久久久久久久久久久久久| 国产爱豆传媒在线观看| 夜夜爽天天搞| 人人妻人人看人人澡| 久久人人爽人人爽人人片va| 日本一本二区三区精品| 国产精品人妻久久久影院| 在线播放国产精品三级| 婷婷六月久久综合丁香| 极品教师在线视频| av国产免费在线观看| 亚洲欧美精品综合久久99| 此物有八面人人有两片| 性色avwww在线观看| 欧美潮喷喷水| 婷婷色综合大香蕉| 国产精品久久视频播放| 亚洲人成网站高清观看| 国产精品一区二区性色av| 少妇裸体淫交视频免费看高清| 观看免费一级毛片| 永久网站在线| 美女xxoo啪啪120秒动态图| 欧美成人免费av一区二区三区| 精品不卡国产一区二区三区| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 三级男女做爰猛烈吃奶摸视频| 哪里可以看免费的av片| 欧洲精品卡2卡3卡4卡5卡区| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 啦啦啦观看免费观看视频高清| 少妇丰满av| 国产成人aa在线观看| 亚洲最大成人中文| 久久久久久久久大av| 99热这里只有是精品在线观看| 熟女电影av网| 国产一区二区激情短视频| 直男gayav资源| 国产午夜精品久久久久久一区二区三区| 毛片一级片免费看久久久久| 中文字幕av在线有码专区| 国产精品,欧美在线| 日韩一区二区三区影片| 久久精品91蜜桃| 亚洲欧美日韩卡通动漫| 国产精品人妻久久久久久| 国产精品av视频在线免费观看| 久久人人爽人人片av| 国产午夜福利久久久久久| 免费人成视频x8x8入口观看| 亚洲无线观看免费| 国产一区二区在线av高清观看| 欧美高清性xxxxhd video| 久久人人精品亚洲av| 亚洲欧美日韩高清专用| 三级男女做爰猛烈吃奶摸视频| 男女边吃奶边做爰视频| 日韩在线高清观看一区二区三区| а√天堂www在线а√下载| 看黄色毛片网站| 九九久久精品国产亚洲av麻豆| 99久久成人亚洲精品观看| 日韩欧美一区二区三区在线观看| 国内精品久久久久精免费| 91久久精品国产一区二区三区| 国产av麻豆久久久久久久| 中文字幕熟女人妻在线| 国产精品国产高清国产av| 国模一区二区三区四区视频| ponron亚洲| 韩国av在线不卡| 美女大奶头视频| 99热这里只有精品一区| 亚洲国产欧美在线一区| 蜜桃亚洲精品一区二区三区| 美女国产视频在线观看| 精品国内亚洲2022精品成人| or卡值多少钱| 国产黄片美女视频| 久久久久网色| av天堂在线播放| 亚洲av不卡在线观看| 国产精品久久久久久精品电影| av在线老鸭窝| 乱人视频在线观看| 欧美激情久久久久久爽电影| 午夜精品在线福利| 亚洲精品乱码久久久久久按摩| 国产v大片淫在线免费观看| 69人妻影院| 日韩成人伦理影院| 长腿黑丝高跟| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆| 可以在线观看的亚洲视频| 精品熟女少妇av免费看| av卡一久久| 91狼人影院| a级毛片免费高清观看在线播放| 国产在视频线在精品| 最新中文字幕久久久久| 欧美色欧美亚洲另类二区| 99久久久亚洲精品蜜臀av| 亚洲五月天丁香| 欧美日韩国产亚洲二区| 麻豆国产97在线/欧美| 老女人水多毛片| 久久久成人免费电影| 人妻久久中文字幕网| av天堂在线播放| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| 美女内射精品一级片tv| 日韩,欧美,国产一区二区三区 | 国产精品一区二区三区四区久久| 老师上课跳d突然被开到最大视频| 五月伊人婷婷丁香| 99久久成人亚洲精品观看| 少妇被粗大猛烈的视频| 热99re8久久精品国产| 日本免费一区二区三区高清不卡| 午夜爱爱视频在线播放| 99视频精品全部免费 在线| 亚洲精品影视一区二区三区av| 两性午夜刺激爽爽歪歪视频在线观看| 长腿黑丝高跟| 国产在线精品亚洲第一网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 男的添女的下面高潮视频| 又爽又黄无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄 | 成年免费大片在线观看| 日韩精品青青久久久久久| 少妇人妻精品综合一区二区 | 国产视频首页在线观看| 亚洲美女搞黄在线观看| 在线观看66精品国产| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 精品一区二区三区视频在线| 国内精品美女久久久久久| 综合色丁香网| 欧美一级a爱片免费观看看| 少妇被粗大猛烈的视频| 51国产日韩欧美| 成年版毛片免费区| 青春草亚洲视频在线观看| 欧美日本视频| 久久综合国产亚洲精品| 91狼人影院| 日本免费a在线| 国产探花极品一区二区| 亚洲av第一区精品v没综合| 少妇高潮的动态图| 一个人观看的视频www高清免费观看| 日本成人三级电影网站| 给我免费播放毛片高清在线观看| 午夜福利高清视频| 日韩一区二区三区影片| 久久国产乱子免费精品| 日本在线视频免费播放| 久久中文看片网| 欧美一级a爱片免费观看看| 麻豆国产97在线/欧美| 亚洲欧美成人精品一区二区| 亚洲av第一区精品v没综合| 国产伦在线观看视频一区| 日本爱情动作片www.在线观看| 变态另类成人亚洲欧美熟女| 午夜免费激情av| 国产精品综合久久久久久久免费| 人妻系列 视频| 国产亚洲5aaaaa淫片| 亚洲一区高清亚洲精品| 日韩中字成人| 又爽又黄无遮挡网站| 成人三级黄色视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品日韩在线中文字幕 | 老熟妇乱子伦视频在线观看| 18禁在线播放成人免费| 国产免费男女视频| 哪里可以看免费的av片| 国产av麻豆久久久久久久| 久久久久久大精品| 日韩在线高清观看一区二区三区| 青春草亚洲视频在线观看| 91精品一卡2卡3卡4卡| 一级毛片久久久久久久久女| 熟妇人妻久久中文字幕3abv| 欧美日韩乱码在线| 久久这里有精品视频免费| 床上黄色一级片| 国产高清三级在线| 久久久久久久久久黄片| 亚洲丝袜综合中文字幕| 变态另类丝袜制服| 久久久久久久亚洲中文字幕| 丰满人妻一区二区三区视频av| 国产大屁股一区二区在线视频| 国产精品久久久久久精品电影| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 内地一区二区视频在线| 免费av毛片视频| 国产成人精品一,二区 | av专区在线播放| 亚洲最大成人av| 搞女人的毛片| 亚洲婷婷狠狠爱综合网| 夜夜看夜夜爽夜夜摸| av女优亚洲男人天堂| 日日啪夜夜撸| 99热全是精品| 欧美又色又爽又黄视频| 国产日本99.免费观看| 国产av在哪里看| 极品教师在线视频| 少妇被粗大猛烈的视频| 老师上课跳d突然被开到最大视频| 我要搜黄色片| 身体一侧抽搐| 直男gayav资源| 超碰av人人做人人爽久久| 寂寞人妻少妇视频99o| 欧美成人免费av一区二区三区| 综合色av麻豆| 国产亚洲精品久久久久久毛片| 国产精品伦人一区二区| 国产一区二区在线av高清观看| 欧美区成人在线视频| 91精品国产九色| 淫秽高清视频在线观看| 久久精品91蜜桃| 欧洲精品卡2卡3卡4卡5卡区|