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

    2種邊界層參數(shù)化方案對1013號超強臺風(fēng)“鲇魚”數(shù)值模擬的敏感性研究

    2017-07-13 05:25胡昊
    安徽農(nóng)業(yè)科學(xué) 2017年24期
    關(guān)鍵詞:鲇魚數(shù)值模擬敏感性

    胡昊

    摘要 利用中國臺灣地區(qū)“侵臺臺風(fēng)之飛機偵察及投落送觀測實驗(DOTSTAR)”探空數(shù)據(jù)與WRF模式(3.8版本),針對1013號超強臺風(fēng)“鲇魚”開展了邊界層參數(shù)化方案的數(shù)值模擬敏感性研究。結(jié)果表明,非局地YSU方案在臺風(fēng)強度、路徑、螺旋雨帶、內(nèi)核尺度以及臺風(fēng)外圍位溫、水汽混合比、風(fēng)向這些方面的模擬效果優(yōu)于局地MYNN2方案;但在風(fēng)速的模擬上,YSU方案與MYNN2方案沒有明顯的優(yōu)劣之分。

    關(guān)鍵詞 超強臺風(fēng)“鲇魚”;數(shù)值模擬;邊界層參數(shù)化方案;敏感性

    Sensitivity Research of Two Planetary Boundary Layer Parameterization Schemes in Numerical Simulation of Super Typhoon 1013 “Megi”

    HU Hao

    (College of Meteorology and Oceanography, PLA University of Science and Technology, Nanjing, Jiangsu 211101)

    Abstract Using the dropsonde data from Dropwindsonde Observations for Typhoon Surveillance near the Taiwan Region (DOTSTAR) and the WRF model (version 3.8), sensitivity research of planetary boundary layer parameterization schemes in numerical simulation of super typhoon “Megi” was carried out. The result showed that the non-local YSU scheme was more capable of simulating typhoon intensity, track, spiral rain band, inner-core size and typhoon periphery structure of potential temperature, water vapor mixing ratio, wind direction than the local MYNN2 scheme. However, there was no obvious difference between YSU and MYNN2 in the simulation of wind speed.

    Key words Super typhoon“Megi”;Numerical simulation;Planetary boundary layer parameterization scheme;Sensitivity

    由于大氣運動學(xué)方程中存在著湍流通量項,使得方程組無法閉合,為了解決方程組的閉合問題,就需要引入邊界層參數(shù)化方案[1]。WRF模式作為一種主流的中尺度數(shù)值模式,其中一共引入了13種邊界層參數(shù)化方案。前人對于邊界層參數(shù)化方案在降水[2-4]、邊界層結(jié)構(gòu)[5-8]、污染擴散[9-10]、臺風(fēng)[11-12]等方面的適用開展了大量的模擬研究。筆者選取非局地YSU方案與局地MYNN2方案這2種常用的邊界層參數(shù)化方案對1013號超強臺風(fēng)“鲇魚”進行模擬與分析,以期為臺風(fēng)周邊區(qū)域邊界層熱力與動力結(jié)構(gòu)的數(shù)值模擬研究提供參考。

    1 資料與方法

    1.1 DOTSTAR計劃觀測情況

    中國臺灣地區(qū)“侵臺臺風(fēng)之飛機偵察及投落送觀測實驗(DOTSTAR,Dropwindsonde Observations for Typhoon Surveillance near the TAiwan Region)”計劃于2010年10月17日05∶40—09∶46共釋放了16個探空儀對臺風(fēng)“鲇魚”進行觀測(圖1)。根據(jù)探空位置和臺風(fēng)中心之間的關(guān)系DOTSTAR觀測航程大致可以分為3段:首先,第1個至第3個探空儀自西向東在臺風(fēng)西側(cè)外圍釋放進行觀測;接著,環(huán)繞臺風(fēng)中心外圍釋放了第4個至第13個探空儀;最后,在臺風(fēng)北側(cè)釋放了第14個至第16個探空儀。于是大致上臺風(fēng)西側(cè)有3個探空儀進行觀測,環(huán)繞臺風(fēng)中心外圍有10個探空儀進行觀測,在臺風(fēng)北側(cè)有3個探空儀進行觀測。

    在探空觀測進行的時候,臺風(fēng)“鲇魚”位于大陸高壓與海上高壓組成的高壓帶南側(cè),在高壓帶南側(cè)偏東風(fēng)的引導(dǎo)下逐漸向西移動,強度上從強臺風(fēng)等級逐漸加強為超強臺風(fēng)。

    1.2 模擬設(shè)置

    為了研究參數(shù)化方案對臺風(fēng)周邊邊界層結(jié)構(gòu)的影響,選取了YSU方案和MYNN2方案這2種較為常用的邊界層參數(shù)化方案對于臺風(fēng)“鲇魚”進行模擬,并將二者分別記為試驗YSU與試驗MYNN2。試驗采用三重移動嵌套網(wǎng)格,由外至內(nèi)將水平分辨率設(shè)置為18、6、2 km,并將格點數(shù)分別設(shè)置為250×250、379×337與268×274。具體模擬區(qū)域如圖2所示。每層嵌套在垂直方向上均設(shè)置為45層。為了獲得對邊界層結(jié)構(gòu)更加詳細的模擬結(jié)果,在 2 000 m高度以下對垂直分層進行加密設(shè)置,共設(shè)置30層。2組試驗除了選擇不同的邊界層參數(shù)化方案之外,其余的參數(shù)化方案均設(shè)置了MM5近地層參數(shù)化方案、RRTM長波輻射方案、Dudhia短波輻射方案、Noah陸面模式。同時在最外層粗網(wǎng)格采用Kain-Fritsch積云參數(shù)化方案和WSM6微物理方案。第二層與最內(nèi)層網(wǎng)格分辨率較高,故未設(shè)置積云參數(shù)化方案,采用了WSM3微物理方案。模擬時間設(shè)置為2010年10月16日20∶00—17日14∶00。16日20∶00—17日05∶00預(yù)留作為模式的spin up時間,而之后的模式輸出數(shù)據(jù)用于比對驗證。

    2 結(jié)果與分析

    2.1 臺風(fēng)模擬情況 由圖3可知,模擬期間臺風(fēng)自西向東移動,并逐漸從向西北方向移動轉(zhuǎn)為向西南方向移動。2組試驗均模擬出臺風(fēng)的西行趨勢,但模擬路徑均比實際路徑偏北,移動速度稍慢。模擬期間JTWC觀測數(shù)據(jù)顯示,臺風(fēng)強度逐漸從16日20∶00的941 hPa加強至17日14∶00的914 hPa。在臺風(fēng)強度的模擬上,2組試驗均比實際偏弱(圖4);試驗YSU成功模擬出臺風(fēng)逐漸加強的過程,模擬得到的臺風(fēng)強度從16日20∶00的989.97 hPa加強至17日14∶00的959.21 hPa;試驗MYNN2卻未能模擬出臺風(fēng)的加強趨勢,模擬期間臺風(fēng)強度從16日20∶00的989.97 hPa減弱至17日14∶00的1 004.22 hPa,與實際觀測嚴重不符。

    利用模式第二層6 km水平分辨率細網(wǎng)格輸出的1 000 m高度處雷達反射率識別臺風(fēng)的螺旋雨帶(圖5),并與F-13衛(wèi)星37V和MTSAT衛(wèi)星紅外云圖合成圖(圖6)進行對比,結(jié)果發(fā)現(xiàn),試驗YSU成功模擬出臺風(fēng)“鲇魚”的螺旋雨帶和眼墻結(jié)構(gòu),與實際觀測最為接近;而試驗MYNN2無法模擬出臺風(fēng)完整的眼墻結(jié)構(gòu),且模擬得到的螺旋雨帶強度偏弱,結(jié)構(gòu)松散,與實際觀測相距甚遠。參考Xu等[13]的方法,將2組試驗輸出的10 m高度處水平風(fēng)速25.7 m/s的等值線包圍的區(qū)域定義為臺風(fēng)的內(nèi)核尺度(圖5中的黑色等值線圈出的區(qū)域)。顯然,只有試驗YSU成功模擬出完整的臺風(fēng)內(nèi)核,大致呈邊界規(guī)則的圓形;而試驗MYNN2無法模擬出完整的臺風(fēng)內(nèi)核。螺旋雨帶和內(nèi)核尺度模擬效果的差異主要是2組試驗?zāi)M得到的臺風(fēng)強度的明顯差異造成的。試驗MYNN2模擬得到的臺風(fēng)強度過弱,于是無法模擬得到結(jié)構(gòu)完整、清晰的螺旋雨帶和內(nèi)核尺度。

    綜上所述,在模擬臺風(fēng)路徑上,2組試驗均模擬出臺風(fēng)西移的過程;但在臺風(fēng)強度、螺旋雨帶以及內(nèi)核尺度上,試驗YSU的模擬效果明顯優(yōu)于試驗MYNN2。

    2.2 邊界層結(jié)構(gòu)模擬效果對比

    試驗YSU、試驗MYNN2在對臺風(fēng)自身結(jié)構(gòu)、強度的模擬上表現(xiàn)出明顯的差異,在此繼續(xù)基于模式第二層6 km水平分辨率細網(wǎng)格輸出數(shù)據(jù)分析這2種邊界層參數(shù)化方案對于臺風(fēng)外圍邊界層內(nèi)部結(jié)構(gòu)的影響。

    2.2.1 邊界層頂高。

    由圖7可知,整體而言,2個參數(shù)化方案模擬得到的海洋上的邊界層頂高高于陸地上的邊界層頂高。DOTSTAR在進行探空觀測時對應(yīng)的局地時間為17日08∶00。陸地上有明顯的邊界層周日循環(huán),而海上卻沒有[1]。探測時大陸、島嶼上空正處在從夜間穩(wěn)定邊界層向白天對流型邊界層過渡發(fā)展的階段,地表不斷升溫,同時低空大氣也逐漸因為地表的長波輻射作用而逐漸升溫。于是夜間的穩(wěn)定型邊界層結(jié)構(gòu)逐漸被破壞,白天的對流型邊界層結(jié)構(gòu)逐漸建立起來。但由于距離日出并沒有過去太長的時間,對流型邊界層并沒有完全形成,所以2組試驗?zāi)M得到的陸上邊界層頂高依然較低。試驗YSU模擬得到的大陸、島嶼上空邊界層頂高大致在500 m以下。試驗MYNN2模擬得到的陸地上空的邊界層頂高稍高于試驗YSU的模擬結(jié)果。在海上,由于海洋巨大的熱容量,所以海表溫度的周日循環(huán)變化相比較地表溫度而言小很多,海洋作為一個持續(xù)的熱源持續(xù)對低空大氣進行加熱,始終維持著海洋上空的對流型邊界層結(jié)構(gòu),所以海洋上空沒有明顯的邊界層結(jié)構(gòu)的周日循環(huán)。該研究由于模擬時間較短(僅18 h),所以在模式設(shè)置時沒有選擇海溫更新選項。海洋持續(xù)加熱低空大氣,所以2組試驗?zāi)M得到的海洋上空的邊界層頂高一致明顯高于陸地。由于試驗YSU成功模擬出臺風(fēng)的螺旋雨帶和眼墻結(jié)構(gòu),而在螺旋雨帶之間以及臺風(fēng)眼中盛行明顯的下沉運動,所以在海上螺旋雨帶和臺風(fēng)眼內(nèi),試驗YSU模擬得到的邊界層頂高比海上其他區(qū)域明顯偏低。但對于試驗MYNN2而言,由于模擬得到的臺風(fēng)強度嚴重偏弱,模擬得到的螺旋雨帶結(jié)構(gòu)松散,沒有模擬出明顯的眼墻結(jié)構(gòu),所以海上邊界層頂高的模擬無法看出像試驗YSU中模擬得到的螺旋雨帶之間條帶狀的邊界層頂高低值區(qū)和臺風(fēng)眼內(nèi)的邊界層頂高低值區(qū)。

    2.2.2 位溫廓線。

    上文提及,DOTSTAR計劃在對臺風(fēng)“鲇魚”進行觀測時,航程大致可以分為3段,即接近臺風(fēng)向東飛行、環(huán)繞臺風(fēng)飛行、離開臺風(fēng)向西北方向飛行,相對應(yīng)釋放的探空儀位置分別位于臺風(fēng)西側(cè)外圍、環(huán)繞臺風(fēng)外圍、臺風(fēng)北側(cè)外圍。對于這3段航程,臺風(fēng)西側(cè)外圍選取第2個探空儀進行研究,環(huán)繞臺風(fēng)外圍選取第8個和第13個探空儀進行研究,即臺風(fēng)中心外圍南北方向上各一個探空儀;臺風(fēng)外圍北側(cè)選取第16個探空儀進行研究。為了能夠有一個對于邊界層內(nèi)整體模擬情況的研究,在此并沒有嚴格將各個物理量廓線的研究高度限制在海表至邊界層頂高位置上,而是統(tǒng)一研究 2 000 m以下高度各個物理量的廓線。

    從圖8可以看出,臺風(fēng)西側(cè)與南側(cè)的實際探空廓線在2 000 m以下大致表現(xiàn)為均一的中性層結(jié),位溫隨高度變化不大;而在臺風(fēng)北側(cè)釋放的2個探空儀探測到的位溫廓線在1 000 m以下表現(xiàn)為中性層結(jié),在1 000~1 300 m附近表現(xiàn)為強逆溫層結(jié)。2組試驗中試驗YSU對于位溫廓線的模擬比試驗MYNN2更加與實際探空接近,尤其是模擬出了2個在臺風(fēng)北側(cè)釋放的探空儀觀測到的強逆溫層結(jié),同時也準確模擬出強逆溫層結(jié)所在的垂直位置。整體而言,試驗MYNN2模擬得到的位溫偏低,與探空存在一定差距。并且在臺風(fēng)南側(cè)試驗MYNN2模擬出虛假的強逆溫層結(jié)。所以在位溫模擬上,YSU方案優(yōu)于MYNN2方案。

    2.2.3 水汽混合比廓線。

    由圖9可知,在臺風(fēng)西側(cè)和南側(cè),探空觀測顯示水汽混合比隨著高度的上升而逐漸降低。在臺風(fēng)北側(cè)釋放的2個探空儀觀測到1 000 m以下水汽混合比基本維持不變,只是小幅下降,而在1 000 m附近存在明顯的水汽混合比隨高度銳減層。在對臺風(fēng)西側(cè)和南側(cè)水汽混合比廓線的模擬上,1 300 m高度以下試驗YSU模擬結(jié)果與實際探空較為接近;但在1 300 m高度之上試驗YSU與試驗MYNN2模擬結(jié)果接近,均比實際探空偏干。試驗YSU成功模擬出臺風(fēng)北側(cè)2個探空儀觀測到的水汽混合比隨高度銳減層,且模擬得到的層結(jié)位置也與實際探空吻合得很好;但試驗MYNN2僅能模擬出其中的一個水汽混合比隨高度銳減層,且位置偏高,與實際觀測偏差較大,而且還模擬出臺風(fēng)南側(cè)虛假的水汽混合比隨高度增加層。所以,在水汽混合比的模擬上,YSU方案也同樣優(yōu)于MYNN2方案。

    2.2.4 風(fēng)向廓線。

    從實際探空觀測得到的風(fēng)向廓線(圖10)來看,臺風(fēng)西側(cè)和北側(cè)自下而上為一致的東北風(fēng);而在臺風(fēng)南側(cè),則以自下而上一致的南風(fēng)為主。在臺風(fēng)西側(cè),試驗YSU模擬效果與實際探空觀測較為接近,試驗MYNN2模擬結(jié)果比實際觀測風(fēng)向偏東;在臺風(fēng)南側(cè),1 300 m以下低空試驗MYNN2模擬結(jié)果與實際觀測較為接近,而1 300 m以上2組試驗?zāi)M結(jié)果接近,均比實際風(fēng)向偏東;在臺風(fēng)北側(cè),試驗MYNN2對于第13個探空儀觀測到風(fēng)向的模擬存在自下而上整體偏東的誤差,但對第16個探空儀1 500 m以上高度風(fēng)向的模擬效果卻優(yōu)于試驗YSU。

    2.2.5 風(fēng)速廓線。

    由圖11可見,大多數(shù)情況下試驗YSU模擬得到的風(fēng)速比試驗MYNN2模擬得到的風(fēng)速偏大。除對于第11個探空儀500 m以下和第16個探空儀1 000 m以上風(fēng)速的模擬外,試驗YSU模擬得到的風(fēng)速均比實際探空觀測偏大。2組試驗對于風(fēng)速的模擬效果沒有明顯優(yōu)劣之分。

    3 結(jié)論

    利用DOTSTAR探空數(shù)據(jù)與WRF模式(3.8版本),針對1013號超強臺風(fēng)“鲇魚”開展了邊界層參數(shù)化方案的數(shù)值模擬敏感性研究。結(jié)果表明,非局地YSU方案對于臺風(fēng)強度、路徑、螺旋雨帶、內(nèi)核尺度的模擬明顯優(yōu)于局地MYNN2方案,與實際觀測更加接近。對于臺風(fēng)外圍區(qū)域邊界層內(nèi)部結(jié)構(gòu)的模擬而言,YSU方案模擬位溫、水汽混合比、風(fēng)向效果較好;YSU方案與MYNN2方案對于風(fēng)速的模擬沒有明顯優(yōu)劣區(qū)別。

    參考文獻

    [1] STULL R B.An Introduction to Boundary Layer Meteorology[M].Berlin:Springer Science & Business Media,1988.

    [2] 蔡薌寧,壽紹文,鐘青.邊界層參數(shù)化方案對暴雨數(shù)值模擬的影響[J].南京氣象學(xué)院學(xué)報,2006,29(3):364-370.

    [3] 陳炯,王建捷.邊界層參數(shù)化方案對降水預(yù)報的影響[J].應(yīng)用氣象學(xué)報,2006,17(S1):11-17.

    [4] 徐慧燕,朱業(yè),劉瑞,等.長江下游地區(qū)不同邊界層參數(shù)化方案的試驗研究[J].大氣科學(xué),2013,37(1):149-159.

    [5] 黃文彥,沈新勇,王衛(wèi)國,等.邊界層參數(shù)化方案對邊界層熱力和動力結(jié)構(gòu)特征影響的比較[J].地球物理學(xué)報,2014,57(5):1399-1414.

    [6] 劉夢娟,陳敏.BJ-RUC 系統(tǒng)對北京夏季邊界層的預(yù)報性能評估[J].應(yīng)用氣象學(xué)報,2014,25(2):212-221.

    [7] 黃文彥,沈新勇,王衛(wèi)國,等.邊界層參數(shù)化方案對邊界層熱力和動力結(jié)構(gòu)特征影響的比較[J].地球物理學(xué)報,2014,57(5):1399-1414.

    [8] 黃文彥,吳建秋,吳晶璐,等.多模式對美國中部草原大氣邊界層模擬的對比分析[J].干旱氣象,2016,34(4):710-717.

    [9] 楊貴成.WRF-Chem 中沙塵天氣過程對模式分辨率及邊界層方案的敏感性試驗[J].安徽農(nóng)業(yè)科學(xué),2012,40(6):3462-3466.

    [10] 王穎,隆霄,余曄,等.復(fù)雜地形上氣象場對空氣質(zhì)量數(shù)值模擬結(jié)果影響的研究[J].大氣科學(xué),2013,37(1):14-22.

    [11] 張建海,張立波,龐盛榮.臺風(fēng) “卡努”(0515)加強過程對邊界層參數(shù)化方案的敏感性試驗[J].臺灣海峽,2007,26(1):26-35.

    [12] 王晨稀.邊界層參數(shù)化影響 “梅花” 臺風(fēng)的敏感性試驗[J].地球科學(xué)進展,2013,28(2):197-208.

    [13] XU J,WANG Y.Sensitivity of tropical cyclone inner-core size and intensity to the radial distribution of surface entropy flux[J].Journal of the atmospheric sciences,2010,67(6):1831-1852.

    猜你喜歡
    鲇魚數(shù)值模擬敏感性
    智多寶智斗鲇魚精
    釔對Mg-Zn-Y-Zr合金熱裂敏感性影響
    AH70DB鋼焊接熱影響區(qū)組織及其冷裂敏感性
    如何培養(yǎng)和提高新聞敏感性
    微小RNA與食管癌放射敏感性的相關(guān)研究
    為什么“改革的鲇魚”不生猛?
    精品久久久久久成人av| 成人三级做爰电影| 欧美日韩一级在线毛片| 成人手机av| 久久久久久大精品| 久久亚洲真实| www.精华液| 黄色毛片三级朝国网站| 亚洲熟女毛片儿| 亚洲美女黄片视频| 亚洲中文av在线| 中国美女看黄片| 最新在线观看一区二区三区| 国产精品电影一区二区三区| 成人av一区二区三区在线看| 91麻豆精品激情在线观看国产 | 丰满迷人的少妇在线观看| 欧美人与性动交α欧美软件| 亚洲自拍偷在线| 丝袜在线中文字幕| 一级a爱视频在线免费观看| av免费在线观看网站| 国产亚洲av高清不卡| 久久精品亚洲熟妇少妇任你| 中文字幕高清在线视频| 国产精品亚洲av一区麻豆| 91国产中文字幕| 91av网站免费观看| 嫩草影院精品99| 欧美激情高清一区二区三区| 亚洲一区高清亚洲精品| 一本大道久久a久久精品| 亚洲熟女毛片儿| 丰满饥渴人妻一区二区三| 99久久人妻综合| 精品一区二区三卡| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 国产99白浆流出| 老熟妇仑乱视频hdxx| 最好的美女福利视频网| 日韩欧美国产一区二区入口| 欧美成人午夜精品| 午夜两性在线视频| 18禁国产床啪视频网站| 欧美精品一区二区免费开放| svipshipincom国产片| 欧美最黄视频在线播放免费 | av天堂久久9| 亚洲男人的天堂狠狠| 久久青草综合色| 在线观看www视频免费| 另类亚洲欧美激情| 91国产中文字幕| 亚洲欧洲精品一区二区精品久久久| 亚洲va日本ⅴa欧美va伊人久久| 99国产极品粉嫩在线观看| 99在线人妻在线中文字幕| 嫁个100分男人电影在线观看| 少妇的丰满在线观看| 在线观看舔阴道视频| 亚洲全国av大片| 亚洲精品美女久久av网站| 啪啪无遮挡十八禁网站| 色精品久久人妻99蜜桃| 两性夫妻黄色片| 国产欧美日韩综合在线一区二区| 久久久久久久久中文| 少妇粗大呻吟视频| 露出奶头的视频| 别揉我奶头~嗯~啊~动态视频| 精品第一国产精品| 欧美日韩亚洲高清精品| 淫秽高清视频在线观看| 国产成人av教育| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 国产精品国产av在线观看| 最近最新免费中文字幕在线| 国产高清国产精品国产三级| 热re99久久精品国产66热6| 亚洲五月天丁香| 97超级碰碰碰精品色视频在线观看| 亚洲中文日韩欧美视频| 村上凉子中文字幕在线| 高清黄色对白视频在线免费看| 黄色片一级片一级黄色片| 50天的宝宝边吃奶边哭怎么回事| 好男人电影高清在线观看| 19禁男女啪啪无遮挡网站| 99在线视频只有这里精品首页| 50天的宝宝边吃奶边哭怎么回事| av天堂在线播放| 精品久久久久久久毛片微露脸| 久99久视频精品免费| av国产精品久久久久影院| 成人av一区二区三区在线看| 国产一区二区激情短视频| 中文字幕人妻熟女乱码| 黄色丝袜av网址大全| 韩国av一区二区三区四区| 亚洲熟女毛片儿| 亚洲精品在线观看二区| 亚洲国产精品一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 亚洲五月色婷婷综合| 国产精品久久久人人做人人爽| 免费看十八禁软件| 亚洲精品av麻豆狂野| 成人国产一区最新在线观看| 国产免费现黄频在线看| 精品午夜福利视频在线观看一区| 一边摸一边做爽爽视频免费| 男女做爰动态图高潮gif福利片 | 亚洲中文字幕日韩| 三上悠亚av全集在线观看| 婷婷六月久久综合丁香| 一边摸一边抽搐一进一小说| 精品国产乱子伦一区二区三区| 国产99白浆流出| 成年人免费黄色播放视频| 欧美人与性动交α欧美精品济南到| 9色porny在线观看| 神马国产精品三级电影在线观看 | 在线播放国产精品三级| 三上悠亚av全集在线观看| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 久久久久国产精品人妻aⅴ院| 亚洲国产毛片av蜜桃av| 久久久久国内视频| 看黄色毛片网站| 午夜免费成人在线视频| 精品福利永久在线观看| xxxhd国产人妻xxx| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三区在线| 久久久精品欧美日韩精品| 亚洲免费av在线视频| 国产精品乱码一区二三区的特点 | 人成视频在线观看免费观看| 日日夜夜操网爽| 日韩精品青青久久久久久| av在线天堂中文字幕 | 99国产极品粉嫩在线观看| 少妇裸体淫交视频免费看高清 | 黑人操中国人逼视频| 国产乱人伦免费视频| 成人av一区二区三区在线看| 精品久久久久久成人av| 人人澡人人妻人| 老熟妇乱子伦视频在线观看| 涩涩av久久男人的天堂| 18美女黄网站色大片免费观看| 久久人妻av系列| 国产成人欧美| 男女之事视频高清在线观看| 国产精品永久免费网站| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲综合一区二区三区_| 老司机午夜福利在线观看视频| 97人妻天天添夜夜摸| 国产91精品成人一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影 | 中亚洲国语对白在线视频| 黑人欧美特级aaaaaa片| 午夜精品国产一区二区电影| 美女高潮喷水抽搐中文字幕| 这个男人来自地球电影免费观看| 精品久久蜜臀av无| www.熟女人妻精品国产| 无限看片的www在线观看| 999久久久国产精品视频| 亚洲,欧美精品.| 亚洲五月天丁香| 中出人妻视频一区二区| 免费av中文字幕在线| 一级片免费观看大全| 精品久久久久久成人av| 窝窝影院91人妻| 久久精品91蜜桃| 亚洲片人在线观看| av超薄肉色丝袜交足视频| 久久久国产成人精品二区 | 日本免费一区二区三区高清不卡 | 精品一区二区三区四区五区乱码| 757午夜福利合集在线观看| 久久狼人影院| 天天添夜夜摸| 一区二区三区精品91| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 亚洲av美国av| 亚洲精华国产精华精| 亚洲男人的天堂狠狠| 日韩中文字幕欧美一区二区| 久久久精品欧美日韩精品| 欧美中文综合在线视频| 成人av一区二区三区在线看| 国产精品 国内视频| 51午夜福利影视在线观看| 国产精品永久免费网站| 男女午夜视频在线观看| 久久性视频一级片| a级毛片黄视频| 亚洲激情在线av| 在线观看免费午夜福利视频| 欧美激情久久久久久爽电影 | 国产伦人伦偷精品视频| 国产精品永久免费网站| 国产野战对白在线观看| 十八禁人妻一区二区| 色播在线永久视频| 国产av一区二区精品久久| 久久香蕉激情| 老司机在亚洲福利影院| 成在线人永久免费视频| 亚洲欧美日韩高清在线视频| 久久亚洲精品不卡| 午夜老司机福利片| www日本在线高清视频| 免费在线观看完整版高清| 日本vs欧美在线观看视频| av免费在线观看网站| 99在线人妻在线中文字幕| 欧美中文日本在线观看视频| 一进一出抽搐gif免费好疼 | 精品福利观看| 99精品欧美一区二区三区四区| 桃红色精品国产亚洲av| 亚洲av五月六月丁香网| 国内毛片毛片毛片毛片毛片| 亚洲 欧美一区二区三区| 欧美在线黄色| 国产欧美日韩一区二区精品| 免费在线观看日本一区| 日韩中文字幕欧美一区二区| 欧美日本亚洲视频在线播放| 波多野结衣av一区二区av| 欧美激情高清一区二区三区| 桃色一区二区三区在线观看| 国产1区2区3区精品| 精品国产超薄肉色丝袜足j| 亚洲激情在线av| 久久性视频一级片| 亚洲色图综合在线观看| 丰满的人妻完整版| 中文字幕色久视频| 亚洲 国产 在线| 久久精品成人免费网站| www日本在线高清视频| 亚洲熟妇熟女久久| 久久人妻av系列| 日韩成人在线观看一区二区三区| 精品熟女少妇八av免费久了| 亚洲人成网站在线播放欧美日韩| 51午夜福利影视在线观看| 脱女人内裤的视频| 1024香蕉在线观看| 亚洲国产中文字幕在线视频| 国产精品久久久av美女十八| aaaaa片日本免费| 精品国产亚洲在线| 中国美女看黄片| 亚洲一区二区三区不卡视频| 美女高潮到喷水免费观看| 国产av一区在线观看免费| 亚洲美女黄片视频| 国产精品秋霞免费鲁丝片| 亚洲熟女毛片儿| а√天堂www在线а√下载| 人妻丰满熟妇av一区二区三区| 精品一区二区三卡| e午夜精品久久久久久久| 精品国产亚洲在线| 丝袜人妻中文字幕| 激情在线观看视频在线高清| 欧美日韩乱码在线| 久久精品国产综合久久久| 亚洲av美国av| 久久午夜综合久久蜜桃| 国产av一区二区精品久久| 成人国语在线视频| 波多野结衣高清无吗| 中文欧美无线码| 又黄又爽又免费观看的视频| 久久人妻熟女aⅴ| 麻豆一二三区av精品| 宅男免费午夜| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 色婷婷av一区二区三区视频| 韩国精品一区二区三区| 国产精品美女特级片免费视频播放器 | 成人18禁在线播放| 亚洲精品久久成人aⅴ小说| av片东京热男人的天堂| 亚洲伊人色综图| 男人操女人黄网站| 国产无遮挡羞羞视频在线观看| cao死你这个sao货| 日韩国内少妇激情av| 美女高潮到喷水免费观看| 黑人欧美特级aaaaaa片| 欧美老熟妇乱子伦牲交| 正在播放国产对白刺激| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 他把我摸到了高潮在线观看| 精品久久久久久成人av| 国产不卡一卡二| 夜夜夜夜夜久久久久| 亚洲三区欧美一区| 欧美激情极品国产一区二区三区| 两人在一起打扑克的视频| 丁香六月欧美| 黄色 视频免费看| 99在线人妻在线中文字幕| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产看品久久| 黄片大片在线免费观看| 在线观看免费高清a一片| 成人三级黄色视频| 一边摸一边抽搐一进一小说| 欧美日韩瑟瑟在线播放| 亚洲七黄色美女视频| 9191精品国产免费久久| 午夜福利免费观看在线| 黄片播放在线免费| 美女高潮喷水抽搐中文字幕| 黄片播放在线免费| 日本精品一区二区三区蜜桃| av天堂久久9| 久久亚洲精品不卡| 欧美另类亚洲清纯唯美| 精品久久久久久,| 免费在线观看亚洲国产| 神马国产精品三级电影在线观看 | 色哟哟哟哟哟哟| 69av精品久久久久久| 9191精品国产免费久久| 99精国产麻豆久久婷婷| 欧美乱码精品一区二区三区| 嫩草影院精品99| 在线观看免费高清a一片| 久久精品影院6| 最新美女视频免费是黄的| 啦啦啦免费观看视频1| 97超级碰碰碰精品色视频在线观看| 999久久久精品免费观看国产| 视频在线观看一区二区三区| 欧美精品亚洲一区二区| 麻豆av在线久日| 男女高潮啪啪啪动态图| 欧美av亚洲av综合av国产av| 欧美一级毛片孕妇| 亚洲全国av大片| 麻豆国产av国片精品| 精品国产超薄肉色丝袜足j| 桃红色精品国产亚洲av| 妹子高潮喷水视频| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久国产精品视频| 欧美色视频一区免费| 国产无遮挡羞羞视频在线观看| 嫁个100分男人电影在线观看| 一级黄色大片毛片| 亚洲aⅴ乱码一区二区在线播放 | 免费看十八禁软件| 夜夜看夜夜爽夜夜摸 | 免费日韩欧美在线观看| 欧美成人性av电影在线观看| 久久久久国产精品人妻aⅴ院| 亚洲七黄色美女视频| 日韩国内少妇激情av| 久久人妻熟女aⅴ| 老司机深夜福利视频在线观看| 欧美日韩亚洲综合一区二区三区_| 91在线观看av| 男女之事视频高清在线观看| 免费日韩欧美在线观看| 老鸭窝网址在线观看| 欧美性长视频在线观看| cao死你这个sao货| 久久国产精品影院| 中文欧美无线码| 国产99白浆流出| 制服诱惑二区| 91字幕亚洲| 午夜老司机福利片| 久久久国产欧美日韩av| 成人亚洲精品av一区二区 | 黄片播放在线免费| 国产黄a三级三级三级人| 国产99久久九九免费精品| 99在线人妻在线中文字幕| av国产精品久久久久影院| 日韩三级视频一区二区三区| 在线国产一区二区在线| 日本免费a在线| 中文字幕精品免费在线观看视频| 亚洲七黄色美女视频| 色老头精品视频在线观看| 性色av乱码一区二区三区2| 欧美性长视频在线观看| 亚洲精华国产精华精| 水蜜桃什么品种好| 在线观看免费日韩欧美大片| 久久久久久人人人人人| 性少妇av在线| 欧美在线黄色| 亚洲欧洲精品一区二区精品久久久| 国产欧美日韩一区二区三| 国产日韩一区二区三区精品不卡| 99在线人妻在线中文字幕| 国产精品野战在线观看 | 亚洲午夜精品一区,二区,三区| av网站免费在线观看视频| 国产乱人伦免费视频| 亚洲一区二区三区欧美精品| 亚洲精品中文字幕在线视频| 国产成人一区二区三区免费视频网站| 久久人人精品亚洲av| 大码成人一级视频| 一区二区日韩欧美中文字幕| 国产aⅴ精品一区二区三区波| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 90打野战视频偷拍视频| 欧美人与性动交α欧美软件| 国产精品 国内视频| 中文亚洲av片在线观看爽| 久久九九热精品免费| 色哟哟哟哟哟哟| 99久久国产精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品 欧美亚洲| 国产精品免费一区二区三区在线| 亚洲精品一卡2卡三卡4卡5卡| 啪啪无遮挡十八禁网站| 天堂俺去俺来也www色官网| 男人操女人黄网站| av中文乱码字幕在线| www国产在线视频色| 亚洲精品久久成人aⅴ小说| av天堂久久9| 长腿黑丝高跟| 韩国av一区二区三区四区| 亚洲精品中文字幕一二三四区| 国产精品99久久99久久久不卡| 桃红色精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 狂野欧美激情性xxxx| av欧美777| 国产高清videossex| 88av欧美| 亚洲av成人一区二区三| 久久人妻熟女aⅴ| 亚洲欧美精品综合一区二区三区| 一区在线观看完整版| 久久欧美精品欧美久久欧美| 婷婷丁香在线五月| 18禁国产床啪视频网站| 国产1区2区3区精品| 久久中文字幕人妻熟女| 曰老女人黄片| 国产色视频综合| 91字幕亚洲| 国产精品亚洲av一区麻豆| 男女下面进入的视频免费午夜 | 久久亚洲真实| 国产无遮挡羞羞视频在线观看| 日韩有码中文字幕| 涩涩av久久男人的天堂| 多毛熟女@视频| 岛国视频午夜一区免费看| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 欧美乱码精品一区二区三区| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区mp4| 又大又爽又粗| 久久久久久免费高清国产稀缺| 免费在线观看黄色视频的| 午夜影院日韩av| 亚洲国产中文字幕在线视频| 国产精品一区二区三区四区久久 | 久久午夜亚洲精品久久| 波多野结衣高清无吗| 无人区码免费观看不卡| 国产一区二区三区视频了| 亚洲五月天丁香| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 中文字幕最新亚洲高清| 男人舔女人下体高潮全视频| 国产亚洲欧美98| 久久人人97超碰香蕉20202| tocl精华| 黄网站色视频无遮挡免费观看| 久久久久久久久久久久大奶| 久久久久国产精品人妻aⅴ院| 51午夜福利影视在线观看| 老司机深夜福利视频在线观看| 黑丝袜美女国产一区| 一区二区三区国产精品乱码| 国产欧美日韩精品亚洲av| 欧美久久黑人一区二区| 国产精品久久久人人做人人爽| 久久伊人香网站| 欧美av亚洲av综合av国产av| 国产成人一区二区三区免费视频网站| 国产无遮挡羞羞视频在线观看| 欧美在线一区亚洲| 欧美 亚洲 国产 日韩一| 亚洲人成77777在线视频| 日日摸夜夜添夜夜添小说| 亚洲男人天堂网一区| av天堂久久9| 亚洲激情在线av| 国产高清videossex| 久久精品国产亚洲av香蕉五月| 校园春色视频在线观看| 精品一品国产午夜福利视频| 成人18禁在线播放| 国产av一区二区精品久久| 国产人伦9x9x在线观看| a级片在线免费高清观看视频| 国产高清国产精品国产三级| 亚洲视频免费观看视频| 成人18禁高潮啪啪吃奶动态图| 宅男免费午夜| 国产成人av激情在线播放| 日本黄色日本黄色录像| 看免费av毛片| 久久久久久久午夜电影 | 午夜福利免费观看在线| 91精品国产国语对白视频| 一二三四社区在线视频社区8| 午夜成年电影在线免费观看| 亚洲精品久久成人aⅴ小说| 在线播放国产精品三级| 在线观看一区二区三区激情| 久久草成人影院| 国产又色又爽无遮挡免费看| svipshipincom国产片| 精品一品国产午夜福利视频| 多毛熟女@视频| 亚洲精品一二三| 国产成人av激情在线播放| 亚洲熟妇中文字幕五十中出 | 99国产精品免费福利视频| xxxhd国产人妻xxx| 久久久国产精品麻豆| av电影中文网址| 69av精品久久久久久| 97人妻天天添夜夜摸| 国产精品成人在线| 久久婷婷成人综合色麻豆| 欧美精品一区二区免费开放| 麻豆久久精品国产亚洲av | 人人妻,人人澡人人爽秒播| e午夜精品久久久久久久| 99久久人妻综合| 午夜精品久久久久久毛片777| 老司机亚洲免费影院| 日韩欧美免费精品| 真人做人爱边吃奶动态| 国产精品久久久人人做人人爽| 国产精品 欧美亚洲| 亚洲精品一二三| 99精品在免费线老司机午夜| 国产成人精品在线电影| 成人永久免费在线观看视频| 满18在线观看网站| 一级a爱视频在线免费观看| 嫁个100分男人电影在线观看| 不卡一级毛片| 欧美日韩视频精品一区| 亚洲专区中文字幕在线| 日韩一卡2卡3卡4卡2021年| 久久久久亚洲av毛片大全| 性少妇av在线| 18禁观看日本| 夜夜夜夜夜久久久久| 国产亚洲精品第一综合不卡| 一级作爱视频免费观看| 伊人久久大香线蕉亚洲五| 国产精品亚洲av一区麻豆| 亚洲成人免费av在线播放| 中文字幕色久视频| 黄色女人牲交| 国产在线精品亚洲第一网站| 久久精品影院6| 咕卡用的链子| cao死你这个sao货| 成人特级黄色片久久久久久久| 一级作爱视频免费观看| 日本五十路高清| 国产亚洲欧美精品永久| 亚洲精品一二三| 一级毛片高清免费大全| 视频区欧美日本亚洲| 又黄又粗又硬又大视频| 免费人成视频x8x8入口观看| 久久精品91蜜桃| x7x7x7水蜜桃| 99re在线观看精品视频| 日韩欧美在线二视频| www.999成人在线观看|