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

    WRF模式對(duì)西南地區(qū)干旱事件的模擬研究

    2015-03-20 00:50:52王紅麗廖留峰
    中低緯山地氣象 2015年6期
    關(guān)鍵詞:陸面西南地區(qū)降水

    王紅麗,廖留峰

    (1.貴州省氣象服務(wù)中心,貴州 貴陽 550002;2.貴州省山地環(huán)境氣候研究所,貴州 貴陽 550002)

    ?

    WRF模式對(duì)西南地區(qū)干旱事件的模擬研究

    王紅麗1,廖留峰2

    (1.貴州省氣象服務(wù)中心,貴州 貴陽 550002;2.貴州省山地環(huán)境氣候研究所,貴州 貴陽 550002)

    2009—2010年在中國西南地區(qū)發(fā)生了嚴(yán)重的秋冬春三季持續(xù)性干旱,為了評(píng)估WRF模式兩種陸面方案在中國西南地區(qū)干旱事件中的模擬性能,該文采用Noah陸面方案和Pleim_Xiu陸面方案兩種WRF模擬方案對(duì)西南干旱事件進(jìn)行模擬研究。兩種陸面方案都合理的模擬出了西南地區(qū)極端干旱事件中的降水空間分布與時(shí)間序列,但與觀測資料相比仍然存在一定的誤差。兩種陸面方案的相互比較表明,在極端干旱事件模擬中,Noah陸面方案在降水模擬中有著一定的改進(jìn)作用。在對(duì)西南干旱時(shí)期的蒸發(fā)模擬當(dāng)中,Noah和Px方案相比有著更低的誤差,采用Noah陸面方案能更好的模擬蒸發(fā)狀況。

    WRF模式;陸面方案;西南干旱

    1 引言

    陸地表面結(jié)構(gòu)特征多變,植被性質(zhì)復(fù)雜,下墊面分布極不均勻,是地球系統(tǒng)中一個(gè)重要的部分。陸面與大氣及其他圈層之間存在著各種復(fù)雜的相互作用,動(dòng)量、水汽、能量等物理量的交換,以及輻射傳輸對(duì)于大氣環(huán)流及氣候變化產(chǎn)生顯著的影響[1-3],許多研究工作表明陸面的地表特征、物理和生化過程對(duì)于天氣及氣候變化的影響是顯著的[4-6]。中國有著如此廣闊的幅員,因此陸面狀況必將體現(xiàn)出紛呈復(fù)雜的特征。

    大量研究認(rèn)為陸面過程對(duì)于氣候變化具有重要意義,丁一匯等[7]通過陸面過程和區(qū)域氣候模式RegCM的耦合試驗(yàn),對(duì)中國中部和東部夏季降水區(qū)域氣候進(jìn)行了模擬分析,發(fā)現(xiàn)耦合了陸面過程的區(qū)域氣候模式能較好的描述中國夏季強(qiáng)降水區(qū)與氣候特征。胡婭敏等[8]等在用區(qū)域氣候模式對(duì)中國東部夏季區(qū)域氣候進(jìn)行模擬時(shí),發(fā)現(xiàn)通過修改陸面過程土壤濕度處理,加入土壤濕度同化資料,能提高地表氣溫和降水的模擬效果。

    關(guān)于陸面過程在WRF模式模擬暴雨中的作用有很多研究[9-10]。然而目前對(duì)于WRF模式在干旱事件的模擬中表現(xiàn)如何還不甚清楚,對(duì)于何種陸面方案在干旱事件模擬中表現(xiàn)更好也缺乏深入的了解。近年來西南坪區(qū)頻發(fā)干旱災(zāi)害[11,12],有必要采取區(qū)域模式耦合陸面方案來對(duì)此加強(qiáng)研究。眾所周知,在2009年秋季至2010年春季,中國西南地區(qū)遭受了嚴(yán)重的秋冬春三季持續(xù)性干旱,造成了嚴(yán)重的經(jīng)濟(jì)損失,我們選取這個(gè)事件來驗(yàn)證WRF模式不同陸面方案模擬西南地區(qū)極端干旱事件的性能。

    2 模式及方案設(shè)計(jì)

    2.1 WRF模式及物理過程介紹

    WRF模式(Weather Research and Forecasting model)是美國國家大氣研究中心(NCAR)、美國國家海洋大氣局(NOAA)的環(huán)境預(yù)報(bào)中心(NCEP)和預(yù)測系統(tǒng)實(shí)驗(yàn)室(FSL)等多個(gè)研究機(jī)構(gòu)聯(lián)合開發(fā)的新一代中尺度天氣和氣候模式。

    WRF 模式內(nèi)部物理參數(shù)化方案比較豐富,考慮的物理過程也比較全面。主要的物理過程及參數(shù)化方案包含陸面過程方案、輻射方案、微物理過程方案,積云對(duì)流參數(shù)化方案、邊界層方案。

    輻射方案中包含了短波輻射和長波輻射兩種方案。短波方案主要有RRTMG方案、Dudhia方案、Goddard方案、CAM方案、New Goddard方案等,長波輻射包括RRTM方案、RRTMG方案、CAM方案、New Goddard方案、FLG(UCLA)方案、Held-Suarez方案、GFDL等方案。

    微物理過程包括水蒸汽、云及降水等過程,在WRFV3.4以后的版本中,包含了十余種微物理過程,包括Purdue Lin、WSM3、WSM5、WSM6、Eta(Grid-scale Cloud and Precipitation)GCP、Thompson、Goddard和Morrison 2-Moment等。

    WRF模式中的積云對(duì)流參數(shù)化方案也很豐富,主要有Kain-Fritsch方案、Betts-Miller-Janjic方案、Grell-Devenyi方案、Grell-3方案、Zhang-McFarlane方案、New SAS方案和New SAS(HWRF)方案等。

    土地表面和大氣之間的動(dòng)量、熱量和水汽的交換主要發(fā)生在邊界層,所以邊界層方案也是非常重要的,邊界層的描述是否合理直接影響模式對(duì)于對(duì)流層低層的氣象要素的模擬,WRF模式中有MRF方案、YSU方案、MYJ方案、GFS方案、QNSE方案、MYNN2等邊界層方案。

    2.2 陸面方案介紹

    在模式模擬中,缺乏陸面過程的參與會(huì)顯著的影響模擬結(jié)果的準(zhǔn)確性,經(jīng)過氣象工作者多年的研究,陸續(xù)發(fā)展了許多描述陸面與大氣之間動(dòng)量、水汽、熱量等交換過程的陸面過程方案。陸面過程經(jīng)過長期的發(fā)展,通過不斷的修改和完善,LSM在陸—?dú)馕锢?、生化相互作用反饋機(jī)制方面更具有有效性和合理性,能夠模擬多種陸面過程,如冰面,雪蓋,植被冠層,土壤以及地表水文過程。

    基于俄勒岡州立大學(xué)的OSU陸面模式,Noah陸面方案由此發(fā)展而來。Noah方案在初始化處理中使用了一部分季或年平均的植被和土壤資料,考慮了植被影響、冠層水分和積雪等下墊面情況,并描述了與之相關(guān)的能量及水分收支過程,對(duì)于土壤溫度和土壤濕度分為4層進(jìn)行預(yù)報(bào),分別為0~10 cm,10~30 cm,30~60 cm以及60~100 cm。Noah方案還能輸出地面徑流和地下徑流等水文過程。

    WRF模式中的Pleim_Xiu陸面方案包括兩層土壤,1 cm厚的地表上層和1 cm~1 m的次層土壤,采用強(qiáng)迫—恢復(fù)法計(jì)算土壤溫度和濕度,該方案參數(shù)化的一個(gè)特征是,水汽通量過程的描述有3個(gè)來源:蒸散、土壤蒸發(fā)及植被冠層蒸發(fā),考慮了植被的作用,方案中植被和土壤參數(shù),是基于土地利用、土地覆蓋變化及土壤質(zhì)地類型獲取的。

    2.3 模擬方案設(shè)計(jì)

    本文中采用WRF Version 3.5,模擬區(qū)域采用蘭伯特坐標(biāo)投影,水平格距分辨率為35 km,中心點(diǎn)為104°E,27°N,區(qū)域格點(diǎn)數(shù)為68×58。由NCEP Final Operational Global Analysis 全球再分析資料(FNL)提供了模式的初始驅(qū)動(dòng)場,空間分辨率為1.0°×1.0°,時(shí)間分辨率為6 h(UTC00時(shí)、06時(shí)、12時(shí)和18時(shí)),模擬時(shí)間是2009年5月1日—2010年4月30日。

    為了評(píng)估WRF模式陸面方案在中國西南地區(qū)干旱事件中的模擬性能,我們設(shè)計(jì)了兩種WRF模擬方案,分別為Noah陸面方案、Pleim-Xiu陸面方案(以下稱Px),其余的參數(shù)化方案均保持一致。輻射方案中采取的長波輻射為CAM方案,短波輻射為Dudhia方案;微物理過程中采取WSM6方案;積云對(duì)流參數(shù)采取Kain-Fritsch方案;行星邊界層方案采取YSU方案。

    WRF模式模擬的降水和蒸發(fā)將同觀測資料進(jìn)行對(duì)比以評(píng)估兩種陸面方案在極端干旱事件模擬當(dāng)中的能力。觀測降水為1901—2013年Climatic Research Unit(CRU)3.22版本0.5°×0.5°降水資料和中國氣象局0.5°×0.5°降水資料。地面蒸發(fā)資料分辨率為0.25°×0.25°,來自全球陸地?cái)?shù)據(jù)同化系統(tǒng)GLDAS(Global Land Data Assimilation System)。本文分析內(nèi)容集中于2009年9月1日—2010年3月30日這一嚴(yán)重干旱時(shí)段。

    3 結(jié)果分析

    3.1 降水

    圖1 2009年9月—2010年3月西南地區(qū)(23°~32°N,100°~108°E)區(qū)域平均降水相較CRU 1901—2013年平均降水的偏差

    為了定量的說明2009—2010年之間中國西南地區(qū)發(fā)生的干旱情況,本文采用了CRU 3.22版本降水資料,分析2009年9月—2010年3月降水量相較于1901—2013年平均對(duì)應(yīng)月降水的異常狀況。從圖1可以看出,在整個(gè)干旱時(shí)段,月降水量相較歷史同期水平都是偏少的,其中9月和10月均偏少40 mm以上,其余月份降水量的偏差相對(duì)較小,但考慮到冬季(12月,1月,2月)降水本就較少,這樣的偏差仍是可觀的。前期的秋季(9月,10月,11月)降水已持續(xù)偏少,進(jìn)入春季(3月)以后降水仍較歷史同期水平偏少,產(chǎn)生較強(qiáng)的累積降水偏少現(xiàn)象,于是造成這一持續(xù)三季的極端干旱事件。

    圖2 西南地區(qū)2009年9月—2010年3月平均觀測降水:(a)觀測,(b)CRU,(c)Noah,(d)Px(單位:mm)

    圖2顯示的是中國西南地區(qū)2009年9月1日—2010年3月30日的平均觀測和模擬降水。通過兩種觀測資料,可以清楚的看出除成都平原和重慶東部區(qū)域以外,整個(gè)西南地區(qū)在2009年9月—2010年3月這7個(gè)月之間累積降水不足300 mm,干旱最嚴(yán)重的區(qū)域中,云南東部,貴州大部,四川南部和廣西西北部累積降水甚至不足200 mm,CRU降水資料和觀測降水資料的整體空間分布格局大致相同,即西南地區(qū)的西—南部區(qū)域降水較少,東—北部區(qū)域降水較多。通過兩種WRF模擬結(jié)果與觀測降水進(jìn)行對(duì)比分析可知,兩種陸面方案的模擬在云南、貴州以及廣西和四川的西北部的效果較好,然而Noah陸面方案和Px陸面方案都高估了四川南部的降水,另外還過高的預(yù)報(bào)了四川盆地的干旱狀況。當(dāng)然值得指出的是降水的大值誤差帶多出現(xiàn)于地形陡峭的區(qū)域,比如川西高原和橫斷山脈,這些區(qū)域陡峭的地形可能是導(dǎo)致降水誤差的原因之一。Noah和Px兩種陸面方案都在西南地區(qū)東北部存在較大誤差,但相比Px陸面方案在四川西北部也存在較大偏差的情況下,Noah陸面方案在這一地區(qū)的降水模擬結(jié)果有所改進(jìn)。

    為了進(jìn)一步評(píng)估WRF模式耦合的Noah陸面方案和Px陸面方案的降水模擬性能,本文對(duì)兩種陸面方案的時(shí)間序列的模擬結(jié)果與觀測資料進(jìn)行了對(duì)比分析(圖3)。從圖3可以看出,在干旱時(shí)期,觀測與CRU兩種降水資料的區(qū)域平均結(jié)果除10月份有一定差別外,其余月份十分接近。兩種陸面方案總體來說都合理的模擬出了降水的時(shí)間序列,但都存在一定的誤差,Px陸面方案表現(xiàn)為每個(gè)月的模擬結(jié)果都較觀測降水偏多,而相較Px陸面方案,Noah陸面方案除了9月的降水較觀測偏多外,其余各月模擬降水都與觀測降水較為接近,特別是10月,11月,12月,1月這4個(gè)月的模擬降水與觀測降水基本上一致。

    圖3 2009年10月至2010年4月西南地區(qū)(23~32°N,100~108°E)區(qū)域平均觀測和模擬降水

    3.2 地面蒸發(fā)

    除了降水,蒸發(fā)在干旱的形成中起著關(guān)鍵的作用,也是干旱狀況的一個(gè)重要表征量,因此我們也評(píng)估了兩種陸面方案模擬的蒸發(fā)。由GLDAS蒸發(fā)資料可以得知,整體空間分布格局為西南地區(qū)的西南部區(qū)域蒸發(fā)較大,東—北部區(qū)域蒸發(fā)較少。通過兩種WRF模擬結(jié)果與GLDAS蒸發(fā)資料進(jìn)行對(duì)比分析可知,Px陸面方案表現(xiàn)為在云南和貴州以及廣西西北部模擬偏少,而在四川的模擬偏多,而相較Px陸面方案,Noah陸面方案則較為合理的給出了蒸發(fā)的空間分布形態(tài),但也存在一定的偏差。

    相對(duì)于GLDAS蒸發(fā)資料,Noah,Px兩種陸面方案在西南地區(qū)干旱時(shí)期(2009年9月—2010年3月)平均蒸發(fā)的均方根誤差(RMSEs)分別為0.45,0.67(圖4)。Noah陸面方案在蒸發(fā)的模擬效果上比Px陸面方案表現(xiàn)較好。

    圖4 西南地區(qū)2009年10月—2010年4月平均地面蒸發(fā)(a)GLDAS,(b)Noah,(c)Px(單位:mm/day)

    4 結(jié)論與討論

    文中我們評(píng)估了WRF模式耦合的Noah和Px兩種陸面方案在中國西南地區(qū)2009年秋季至2010年春季持續(xù)性干旱事件中的模擬性能。主要的結(jié)論如下:

    兩種陸面方案都較為合理的模擬出了西南地區(qū)降水的空間和時(shí)間形態(tài),但與觀測資料相比仍然存在一定的誤差。比較兩種陸面方案對(duì)于降水模擬,發(fā)現(xiàn)對(duì)于降水的空間分布和時(shí)間序列模擬上,比之Px陸面方案,Noah陸面方案表現(xiàn)較好。這表明在西南地區(qū)干旱事件模擬中,考慮采用Noah陸面方案較之Px陸面方案在降水模擬中有著一定的改進(jìn)作用。兩種陸面方案的對(duì)地面蒸發(fā)的模擬對(duì)比分析表明,由于與Px方案相比,Noah陸面方案有著更低的均方根誤差(RMSEs),在對(duì)西南極端干旱時(shí)期的蒸發(fā)模擬當(dāng)中,采用Noah陸面方案更為合理。

    由于我國陸面特征的顯著差異,在具體應(yīng)用陸面過程方案中,恰當(dāng)?shù)奶幚矸绞綉?yīng)是分別考慮不同的地理地形和下墊面特征。但目前來說,現(xiàn)在應(yīng)用的陸面過程模式,考慮的陸面特征比較單一,這樣的陸面過程配置在模擬應(yīng)對(duì)西南地區(qū)這樣具有復(fù)雜下墊面環(huán)境的區(qū)域,難免會(huì)存在不恰當(dāng)和不細(xì)致之處。未來的工作中應(yīng)注意根據(jù)不同的區(qū)域采用適宜的陸面方案,并嘗試發(fā)展能更加全面細(xì)致描述下墊面特征的陸面方案。

    [1] 包慶,劉屹岷,周天軍,等.LASG/IAP大氣環(huán)流譜模式對(duì)陸面過程的敏感性試驗(yàn)[J].大氣科學(xué),2006,30(6):1 078-1 089.

    [2] 布和朝魯,紀(jì)立人,崔茂常.夏季我國干旱、半干旱區(qū)陸面過程能量平衡及其局地大氣環(huán)流[J].氣候與環(huán)境研究,2002,7(1):61-73.

    [3] 符淙斌,董文杰,溫剛,等.全球變化的區(qū)域響應(yīng)和適應(yīng)[J].氣象學(xué)報(bào),2003,61(2):245-250.

    [4] 符淙斌,溫剛.中國北方干旱化的幾個(gè)問題[J].氣候與環(huán)境研究,2002,7(1):22-29.

    [5] 呂達(dá)仁,陳佐忠,陳家宜,等.內(nèi)蒙古半干旱草原土壤-植被-大氣相互作用綜合研究[J].氣象學(xué)報(bào),2005,3(5):571-593.

    [6] 劉輝志,董文杰,符淙斌,等.半干旱地區(qū)吉林通榆“干旱化和有序人類活動(dòng)”長期觀測研究[J].氣候與環(huán)境研究,2004,9(2):378-389.

    [7] 丁一匯,張晶,趙宗慈.一個(gè)改進(jìn)的陸面過程模式及其模擬試驗(yàn)研究第二部分:陸面過程模式與區(qū)域氣候模式的耦合模擬試驗(yàn)[J].氣象學(xué)報(bào),1998,56(4):385-400.

    [8] 胡婭敏,丁一匯,廖菲.土壤濕度資料同化對(duì)中國東部夏季區(qū)域氣候模擬的改進(jìn)[J].科學(xué)通報(bào),2009,54:2 388-2 394.

    [9] 馬紅云,郭品文,宋潔.耦合不同陸面方案的WRF模式對(duì)2007年7月江淮強(qiáng)降水過程的模擬[J].大氣科學(xué),2009,33,557-567.

    [10] 李安泰,何宏讓,張?jiān)?WRF模式陸面參數(shù)擾動(dòng)對(duì)一次西北暴雨影響的數(shù)值模擬[J].高原氣象,2012,31(1):65-75.

    [11] 崔庭,吳古會(huì),岑啟林,熊方,楊玲.2009—2010年黔西南州特大干旱成因分析[J].貴州氣象,2011,35(2):26-29.

    [12] 黃桂東,佘燕.黔南州2013年和2011年兩次夏季干旱對(duì)比分析[J].貴州氣象,2014(6):20-22.

    WRF model’s simulation study on the drought in Southwest China

    WANG Hongli1,LIAO Liufeng2

    (1.Meteorological service center of Guizhou province,Guiyang 550002;
    2.Guizhou Instftute of Mountainous Climate and Envioronment,Guiyang 550002)

    During the 2009—2010,the Southwest China suffered a persistent drought event.To evaluate the performance of two land surface schemes in the severe drought,two experiments were examined to investigate the drought event using the Noah land surface scheme and Pleim_Xiu land surface scheme.The WRF simulations with two land surface schemes can reasonably simulate the spatial pattern and temporal variation of precipitation in Southwest China,although the simulations still show several errors relative to observations.The comparison of the WRF simulations with two different land surface schemes indicated that the Noah land surface scheme showed a slightly better performance in the precipitation during the severe drought simulations.Noah land surface scheme can better reproduce the evaporation than Pleim_Xiu land surface scheme due to the lower errors.

    WRF model;Land surface scheme;Southwest drought

    2015-06-28

    王紅麗(1963—),女,工程師,主要從事氣象服務(wù)與應(yīng)用氣象工作。

    1003-6598(2015)06-0001-05

    P435

    A

    資助信息:國家科技計(jì)劃課題(西南突發(fā)性災(zāi)害應(yīng)急與防控技術(shù)集成與示范)2012BAD20B06。

    猜你喜歡
    陸面西南地區(qū)降水
    基于MCI指數(shù)的西南地區(qū)近60年夏季干旱特征
    西南地區(qū)干濕演變特征及其對(duì)ENSO事件的響應(yīng)
    黑龍江省玉米生長季自然降水與有效降水對(duì)比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    戴永久:砥礪前行,永遠(yuǎn)奮斗
    戴永久: 砥礪前行,永遠(yuǎn)奮斗
    相似廓線法在星載雷達(dá)陸面衰減訂正中的應(yīng)用
    山西隊(duì)?wèi)蛳蛭髂系貐^(qū)的傳播
    中華戲曲(2019年2期)2019-02-06 06:54:16
    基于區(qū)域氣候模式RegCM4的長江流域降水預(yù)報(bào)效果分析
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    Завершено строительство крупнейшего комплексного транспортного узла на юго
    --западе Китая
    中亞信息(2016年7期)2016-10-20 01:41:30
    一级黄色大片毛片| 国产精品野战在线观看| 97超级碰碰碰精品色视频在线观看| av视频在线观看入口| 18禁国产床啪视频网站| 黄色成人免费大全| 国产在线精品亚洲第一网站| 亚洲欧美激情综合另类| 亚洲一码二码三码区别大吗| 免费高清视频大片| 亚洲美女黄片视频| 免费在线观看视频国产中文字幕亚洲| 国产成人啪精品午夜网站| 无人区码免费观看不卡| videosex国产| 午夜福利18| 久久久久久免费高清国产稀缺| 99热这里只有精品一区 | 久久中文看片网| 久久国产精品男人的天堂亚洲| 精品国产亚洲在线| www日本在线高清视频| aaaaa片日本免费| 18美女黄网站色大片免费观看| АⅤ资源中文在线天堂| 午夜视频精品福利| 黑人操中国人逼视频| 午夜久久久在线观看| 1024手机看黄色片| 香蕉丝袜av| 国产精品,欧美在线| 天天躁狠狠躁夜夜躁狠狠躁| 99国产精品一区二区三区| 两个人看的免费小视频| 国产精品亚洲av一区麻豆| 日韩高清综合在线| 黄色视频不卡| 少妇被粗大的猛进出69影院| www.自偷自拍.com| 亚洲三区欧美一区| 国产亚洲精品综合一区在线观看 | 操出白浆在线播放| 国产精品自产拍在线观看55亚洲| 国产高清激情床上av| www.自偷自拍.com| 老司机在亚洲福利影院| 9191精品国产免费久久| 国产在线精品亚洲第一网站| 91大片在线观看| 成人午夜高清在线视频 | 老司机午夜十八禁免费视频| 18美女黄网站色大片免费观看| 啦啦啦 在线观看视频| 国产蜜桃级精品一区二区三区| 国产成人系列免费观看| 亚洲中文字幕一区二区三区有码在线看 | 欧美成人免费av一区二区三区| 久久久久国内视频| 国产一区二区三区在线臀色熟女| 国产午夜精品久久久久久| avwww免费| 亚洲三区欧美一区| 757午夜福利合集在线观看| 99国产精品99久久久久| 午夜福利成人在线免费观看| 国产成人一区二区三区免费视频网站| 两性夫妻黄色片| 午夜福利在线观看吧| 久久精品国产亚洲av高清一级| 亚洲成av片中文字幕在线观看| 1024香蕉在线观看| 搡老岳熟女国产| 自线自在国产av| 亚洲成人久久性| 黄色成人免费大全| 无限看片的www在线观看| 亚洲午夜理论影院| 国产精品野战在线观看| 国产精品 欧美亚洲| 亚洲精品国产一区二区精华液| 在线看三级毛片| 国产97色在线日韩免费| 国产极品粉嫩免费观看在线| 在线观看午夜福利视频| 免费搜索国产男女视频| 久久亚洲真实| 日本a在线网址| 听说在线观看完整版免费高清| 免费高清视频大片| 久久久久久久久久黄片| 亚洲成av片中文字幕在线观看| 在线免费观看的www视频| 一二三四社区在线视频社区8| 美女高潮到喷水免费观看| 亚洲无线在线观看| 久久久国产精品麻豆| 每晚都被弄得嗷嗷叫到高潮| 久久精品亚洲精品国产色婷小说| 国产成+人综合+亚洲专区| 国产亚洲精品av在线| 午夜成年电影在线免费观看| 亚洲国产精品久久男人天堂| 久久精品国产99精品国产亚洲性色| 免费电影在线观看免费观看| 成人亚洲精品一区在线观看| 色综合站精品国产| 欧美日韩中文字幕国产精品一区二区三区| 中文资源天堂在线| 国内久久婷婷六月综合欲色啪| 婷婷丁香在线五月| 亚洲第一欧美日韩一区二区三区| 国产v大片淫在线免费观看| 无人区码免费观看不卡| 欧美色视频一区免费| 欧美成人性av电影在线观看| 欧美性猛交黑人性爽| 欧美激情 高清一区二区三区| 国产精品久久久久久精品电影 | 可以在线观看的亚洲视频| 一边摸一边抽搐一进一小说| 国产精品爽爽va在线观看网站 | 午夜福利免费观看在线| 久久中文看片网| 亚洲激情在线av| 亚洲欧美日韩无卡精品| 精品久久久久久,| 97碰自拍视频| 神马国产精品三级电影在线观看 | 熟女电影av网| 欧美激情久久久久久爽电影| 曰老女人黄片| 亚洲午夜精品一区,二区,三区| 国产精品一区二区精品视频观看| 亚洲午夜精品一区,二区,三区| 熟妇人妻久久中文字幕3abv| 亚洲专区国产一区二区| 搞女人的毛片| 曰老女人黄片| 在线观看免费视频日本深夜| 精品一区二区三区视频在线观看免费| 国产99白浆流出| 自线自在国产av| 亚洲美女黄片视频| 巨乳人妻的诱惑在线观看| 国产又爽黄色视频| 搡老岳熟女国产| 久久精品国产清高在天天线| 午夜免费观看网址| 亚洲精品色激情综合| 欧美黑人欧美精品刺激| 国产三级黄色录像| 午夜亚洲福利在线播放| 亚洲aⅴ乱码一区二区在线播放 | 变态另类丝袜制服| 亚洲精品粉嫩美女一区| 国产欧美日韩精品亚洲av| 19禁男女啪啪无遮挡网站| 亚洲精品粉嫩美女一区| av有码第一页| 琪琪午夜伦伦电影理论片6080| 变态另类丝袜制服| 国产一区二区激情短视频| 久久久久久久午夜电影| 久久中文字幕一级| 老鸭窝网址在线观看| 人成视频在线观看免费观看| 久久久久久久精品吃奶| 精品久久久久久久末码| 男人的好看免费观看在线视频 | 一本久久中文字幕| 嫁个100分男人电影在线观看| 十八禁人妻一区二区| 可以在线观看的亚洲视频| 一级作爱视频免费观看| 国产成人一区二区三区免费视频网站| 欧美成人免费av一区二区三区| 国产又黄又爽又无遮挡在线| 怎么达到女性高潮| 在线观看舔阴道视频| av电影中文网址| 欧美国产精品va在线观看不卡| 午夜两性在线视频| 韩国精品一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲男人天堂网一区| 亚洲av电影在线进入| 日韩视频一区二区在线观看| 国产精品自产拍在线观看55亚洲| 亚洲七黄色美女视频| 搞女人的毛片| 男女那种视频在线观看| 脱女人内裤的视频| 久久 成人 亚洲| 搡老熟女国产l中国老女人| 激情在线观看视频在线高清| 又黄又粗又硬又大视频| 美女国产高潮福利片在线看| 一区二区三区激情视频| www.自偷自拍.com| 丁香欧美五月| 亚洲第一av免费看| 天堂动漫精品| 一个人免费在线观看的高清视频| 欧美激情高清一区二区三区| 夜夜爽天天搞| 99riav亚洲国产免费| 亚洲五月婷婷丁香| 免费在线观看成人毛片| 黄色丝袜av网址大全| 久久久久久久久中文| 日韩国内少妇激情av| 亚洲 欧美 日韩 在线 免费| ponron亚洲| 欧美av亚洲av综合av国产av| 免费在线观看完整版高清| av有码第一页| 亚洲欧美一区二区三区黑人| 99国产精品一区二区蜜桃av| 久久精品91无色码中文字幕| 亚洲国产毛片av蜜桃av| 国产av一区二区精品久久| 国产av一区二区精品久久| 国产亚洲欧美在线一区二区| 精品国产乱码久久久久久男人| 男女视频在线观看网站免费 | 国产精品av久久久久免费| 一二三四在线观看免费中文在| 热re99久久国产66热| 熟妇人妻久久中文字幕3abv| 老司机在亚洲福利影院| 中文亚洲av片在线观看爽| 99久久久亚洲精品蜜臀av| 国产av在哪里看| 国产黄a三级三级三级人| 18禁美女被吸乳视频| cao死你这个sao货| av福利片在线| 老司机福利观看| 国产人伦9x9x在线观看| 久久 成人 亚洲| 色老头精品视频在线观看| 一区二区三区高清视频在线| 人人妻人人澡欧美一区二区| 久热这里只有精品99| 亚洲熟妇熟女久久| 成人免费观看视频高清| 99久久无色码亚洲精品果冻| 无人区码免费观看不卡| 高清在线国产一区| 色播在线永久视频| 国产一区二区三区视频了| 亚洲精品久久成人aⅴ小说| 久久人妻福利社区极品人妻图片| 黄片大片在线免费观看| 免费看日本二区| 久久热在线av| 精品福利观看| 亚洲欧美日韩无卡精品| 90打野战视频偷拍视频| 一本久久中文字幕| 成人国语在线视频| 中文字幕av电影在线播放| 777久久人妻少妇嫩草av网站| 亚洲精品在线观看二区| 99精品欧美一区二区三区四区| 日韩av在线大香蕉| 国产精品美女特级片免费视频播放器 | 美女扒开内裤让男人捅视频| 成人精品一区二区免费| 精品国产超薄肉色丝袜足j| 亚洲国产精品999在线| 久久久水蜜桃国产精品网| 国产高清有码在线观看视频 | 搡老岳熟女国产| 亚洲精品国产精品久久久不卡| 国产精品久久视频播放| 性欧美人与动物交配| 一级毛片精品| 999久久久精品免费观看国产| 国产高清激情床上av| 免费观看人在逋| 精品欧美国产一区二区三| www.精华液| 91在线观看av| 一区二区三区高清视频在线| 色综合亚洲欧美另类图片| 国产成人精品无人区| 久久这里只有精品19| 国产99久久九九免费精品| 在线观看66精品国产| 在线观看一区二区三区| 天天添夜夜摸| 黄色片一级片一级黄色片| 久久午夜综合久久蜜桃| 久久久国产精品麻豆| 一级片免费观看大全| 成熟少妇高潮喷水视频| 韩国av一区二区三区四区| 亚洲中文字幕一区二区三区有码在线看 | 一进一出好大好爽视频| 国产一区二区在线av高清观看| 18禁裸乳无遮挡免费网站照片 | 欧美性长视频在线观看| 日韩欧美 国产精品| 日韩欧美三级三区| 免费人成视频x8x8入口观看| 国产亚洲欧美98| 日日夜夜操网爽| 欧美日韩黄片免| 精品乱码久久久久久99久播| 亚洲精品久久国产高清桃花| 欧美黄色淫秽网站| 女人爽到高潮嗷嗷叫在线视频| ponron亚洲| 亚洲,欧美精品.| av视频在线观看入口| e午夜精品久久久久久久| 啦啦啦观看免费观看视频高清| 熟女电影av网| 草草在线视频免费看| 欧美在线一区亚洲| 熟女少妇亚洲综合色aaa.| 国产精品免费视频内射| 国产精品久久久人人做人人爽| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲av嫩草精品影院| 国产高清videossex| 国产乱人伦免费视频| av欧美777| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 女性被躁到高潮视频| 欧美在线黄色| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区| 中文资源天堂在线| 午夜福利18| 一个人免费在线观看的高清视频| 精品无人区乱码1区二区| 色综合婷婷激情| 亚洲专区中文字幕在线| 国产三级黄色录像| 国产成人一区二区三区免费视频网站| 中文字幕高清在线视频| www.自偷自拍.com| 午夜福利免费观看在线| 欧美色欧美亚洲另类二区| 日本熟妇午夜| 中文字幕人妻丝袜一区二区| 久久性视频一级片| 国产aⅴ精品一区二区三区波| 欧美成人午夜精品| 哪里可以看免费的av片| 在线观看日韩欧美| 精品卡一卡二卡四卡免费| 2021天堂中文幕一二区在线观 | 精品久久久久久久久久免费视频| 国产99白浆流出| 亚洲中文av在线| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| 美女午夜性视频免费| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 国产精品98久久久久久宅男小说| 欧美午夜高清在线| 黄片大片在线免费观看| 看免费av毛片| 91成人精品电影| 中亚洲国语对白在线视频| 久久热在线av| 97碰自拍视频| 欧美日韩乱码在线| 亚洲第一欧美日韩一区二区三区| 欧美激情久久久久久爽电影| 国产精品九九99| 美女扒开内裤让男人捅视频| 亚洲专区中文字幕在线| 俺也久久电影网| 成人欧美大片| 999久久久国产精品视频| av有码第一页| 国产v大片淫在线免费观看| 在线十欧美十亚洲十日本专区| 精华霜和精华液先用哪个| 麻豆av在线久日| 夜夜躁狠狠躁天天躁| 嫩草影视91久久| 亚洲精品中文字幕一二三四区| 亚洲国产中文字幕在线视频| 成人特级黄色片久久久久久久| 在线天堂中文资源库| 久久这里只有精品19| 日韩国内少妇激情av| 欧美在线黄色| 最近最新免费中文字幕在线| 脱女人内裤的视频| 国产三级在线视频| 久久精品国产亚洲av高清一级| 国产久久久一区二区三区| 午夜两性在线视频| 久久久久久国产a免费观看| 成年免费大片在线观看| 搞女人的毛片| 少妇被粗大的猛进出69影院| a级毛片在线看网站| 一级毛片女人18水好多| 99国产极品粉嫩在线观看| 欧美zozozo另类| 夜夜躁狠狠躁天天躁| 国产三级在线视频| 精品久久久久久久末码| 看黄色毛片网站| 男人舔女人下体高潮全视频| 嫁个100分男人电影在线观看| 国产成人精品久久二区二区免费| 国产精品久久久av美女十八| 色精品久久人妻99蜜桃| 国产精品一区二区免费欧美| 国产午夜福利久久久久久| 91成年电影在线观看| 精品久久久久久成人av| 身体一侧抽搐| 日本免费a在线| 日本一本二区三区精品| 免费高清在线观看日韩| 一级作爱视频免费观看| 国产熟女xx| 久久精品国产清高在天天线| 18禁裸乳无遮挡免费网站照片 | 国产极品粉嫩免费观看在线| 久久精品国产亚洲av高清一级| 91av网站免费观看| 久久久久久久午夜电影| 久久狼人影院| 亚洲av日韩精品久久久久久密| 99久久国产精品久久久| 成人一区二区视频在线观看| 搡老岳熟女国产| 大香蕉久久成人网| 亚洲午夜理论影院| 久久精品国产亚洲av香蕉五月| 99热6这里只有精品| 99久久久亚洲精品蜜臀av| av有码第一页| 成熟少妇高潮喷水视频| 波多野结衣高清无吗| 日本 欧美在线| 精品卡一卡二卡四卡免费| 亚洲激情在线av| 最新在线观看一区二区三区| 欧美成人午夜精品| 91av网站免费观看| 精品一区二区三区av网在线观看| 午夜影院日韩av| 亚洲色图 男人天堂 中文字幕| 久久久久国产精品人妻aⅴ院| 天堂√8在线中文| 日本在线视频免费播放| 国产男靠女视频免费网站| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 成年女人毛片免费观看观看9| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 国内精品久久久久精免费| 亚洲成人精品中文字幕电影| 久久久久九九精品影院| 不卡av一区二区三区| 长腿黑丝高跟| √禁漫天堂资源中文www| 欧美又色又爽又黄视频| a级毛片在线看网站| 欧美日本视频| 国产精品一区二区三区四区久久 | 美国免费a级毛片| 成人亚洲精品一区在线观看| 欧美中文综合在线视频| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 18禁美女被吸乳视频| 波多野结衣巨乳人妻| 久久中文看片网| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 欧美成人一区二区免费高清观看 | 女性生殖器流出的白浆| 国产又黄又爽又无遮挡在线| 亚洲国产精品999在线| 久9热在线精品视频| 亚洲色图av天堂| 91av网站免费观看| 亚洲在线自拍视频| 黄色女人牲交| 精品午夜福利视频在线观看一区| 老汉色∧v一级毛片| 国产精品 国内视频| 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密| 色尼玛亚洲综合影院| 精品久久久久久,| 99久久国产精品久久久| 亚洲 欧美 日韩 在线 免费| 成人午夜高清在线视频 | 国产精品二区激情视频| 国产精品久久电影中文字幕| 国产精品久久久av美女十八| 一进一出抽搐动态| 婷婷精品国产亚洲av| 一夜夜www| 在线观看舔阴道视频| 天堂√8在线中文| 欧美大码av| 在线永久观看黄色视频| 久久草成人影院| 波多野结衣av一区二区av| av福利片在线| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| 搡老熟女国产l中国老女人| 亚洲中文av在线| 国产成人av激情在线播放| 国内揄拍国产精品人妻在线 | 亚洲av成人一区二区三| 国产欧美日韩精品亚洲av| 一本精品99久久精品77| 18禁黄网站禁片免费观看直播| 非洲黑人性xxxx精品又粗又长| 999久久久精品免费观看国产| 母亲3免费完整高清在线观看| 天天一区二区日本电影三级| or卡值多少钱| 丁香六月欧美| 亚洲成av片中文字幕在线观看| 精品久久久久久久久久久久久 | 天天添夜夜摸| 一级毛片高清免费大全| 亚洲精品美女久久av网站| 国内久久婷婷六月综合欲色啪| 一a级毛片在线观看| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 99在线视频只有这里精品首页| 成人特级黄色片久久久久久久| 最近最新免费中文字幕在线| 免费无遮挡裸体视频| 欧美成人一区二区免费高清观看 | 黄片大片在线免费观看| 国产精品1区2区在线观看.| 久久久久久久精品吃奶| 国产亚洲精品一区二区www| 欧美日韩中文字幕国产精品一区二区三区| 我的亚洲天堂| 亚洲精品美女久久久久99蜜臀| 国产精品永久免费网站| 日韩欧美一区视频在线观看| 久久草成人影院| 国产又爽黄色视频| 国产精品亚洲av一区麻豆| av中文乱码字幕在线| bbb黄色大片| 999久久久国产精品视频| 亚洲人成电影免费在线| 亚洲专区中文字幕在线| 免费观看精品视频网站| 中文字幕人妻丝袜一区二区| 我的亚洲天堂| 天堂动漫精品| 国产激情欧美一区二区| 黄色视频不卡| 波多野结衣巨乳人妻| 日韩视频一区二区在线观看| 亚洲国产日韩欧美精品在线观看 | 国产成人精品无人区| 成人国产综合亚洲| 18美女黄网站色大片免费观看| 一级a爱视频在线免费观看| 亚洲五月天丁香| 亚洲国产高清在线一区二区三 | 天堂影院成人在线观看| 丰满的人妻完整版| 级片在线观看| 国产私拍福利视频在线观看| 国产又黄又爽又无遮挡在线| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 欧美黄色片欧美黄色片| 亚洲人成网站高清观看| 久久久久久久久久黄片| 国内精品久久久久精免费| 亚洲精品粉嫩美女一区| 色播在线永久视频| 香蕉丝袜av| 老司机在亚洲福利影院| 18禁裸乳无遮挡免费网站照片 | 国产又色又爽无遮挡免费看| 国内久久婷婷六月综合欲色啪| 香蕉av资源在线| 我的亚洲天堂| 亚洲av日韩精品久久久久久密| 亚洲av成人一区二区三| 久久精品国产亚洲av高清一级| 国产精品1区2区在线观看.| a在线观看视频网站| 在线观看www视频免费| 99精品久久久久人妻精品| 国内少妇人妻偷人精品xxx网站 | 久久久国产欧美日韩av| 一卡2卡三卡四卡精品乱码亚洲| 免费电影在线观看免费观看| 两个人看的免费小视频| 91成年电影在线观看|