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

    C波段雷達資料同化對GRAPES中尺度數值模式短臨降水預報的影響

    2022-07-26 00:39:32齊大鵬李彥霖周明飛朱文達
    中低緯山地氣象 2022年3期
    關鍵詞:實況強降水反演

    齊大鵬,楊 靜,李彥霖,周明飛,朱文達

    (貴州省氣象臺,貴州 貴陽 550002)

    0 引言

    短臨預報一直是天氣預報業(yè)務關注的重點和難題。中尺度數值模式由于時空分辨率高,同時又實現(xiàn)了快速循環(huán)同化更新,有效地提高了短臨天氣預報的質量。但是,由于模式初始場中沒有包含非絕熱參數,業(yè)務預報系統(tǒng)面臨的一個首要問題是不可避免的“Spin-up”現(xiàn)象,導致模式前幾個小時的預報特別是降水預報的質量普遍不高[1]。為解決這一問題,目前已研究出一些利用多普勒天氣雷達資料反演大氣三維風場、溫度場和云微物理場的方法,使得模式在初始時刻就包含了云的微物理特征,改進了模式積分初始場,提高了模式的短臨天氣預報質量。如王瑾等[2]用ARPS模式的資料三維同化系統(tǒng)ARPS3DVAR和復雜云分析模塊同化雷達資料后,發(fā)現(xiàn)模式在初始時刻就能夠調整出一個合理的云微物理量場,減少了模式的熱啟動時間。張蕾等[3]的研究表明S波段雷達資料同化對一次暴雨過程的降水預報效果有明顯提高。李平等[4]利用ARPS模式進行了C波段雷達反射率資料同化實驗,結果表明雷達資料引入模式初始場后,模式的回波和降水預報效果都有一定提高。上述研究多基于WRF和ARPS模式,而在GRAPES模式方面,張艷霞等[5]用多普勒雷達資料調整GRAPES模式云參數,發(fā)現(xiàn)模式初始場中水汽場的修訂,尤其是飽和水汽的修訂,對模式預報的水汽場有明顯影響,從而影響到了環(huán)流場,改善了模式的雷達回波和降水的預報。張誠忠等[6]在熱啟動的基礎上,考慮同化雷達反演的水汽后,不僅較好地模擬了前6 h的降水,也可以預報出6~12 h的降雨帶變化趨勢。高郁東等[7]在模式中同化雷達降水率資料后不僅改進了降水分布,而且還能清晰地表現(xiàn)出中尺度對流系統(tǒng)的發(fā)展和消亡。張誠忠等[8]在雷達反演潛熱的基礎上,評估了潛熱對模式短臨降水預報的影響,發(fā)現(xiàn)通過同化雷達反演的潛熱顯著提高了模式短臨降水預報效果。

    上述研究表明通過在模式中同化雷達資料,能夠提高模式的短臨預報能力,然而針對貴州復雜地形條件下雷達資料同化對模式的影響研究較少,單純的數值模擬預報偏差較大[9-10]。本文利用貴州區(qū)域C波段雷達資料,通過1個強降水天氣過程的模擬實驗,分析和討論C波段雷達資料同化對GRAPES中尺度數值模式在貴州的短臨降水預報的影響。

    1 資料和方法

    1.1 資料來源

    本文所用資料有雷達數據、雨量觀測數據和NCEP的GFS全球模式預報數據。雷達資料來自2017年7月9日02時貴州區(qū)域內7部CD型號多普勒天氣雷達(貴陽、遵義、都勻、黔東南、銅仁、畢節(jié)、興義)觀測的基本反射率因子,雷達體掃模式為VCP21(6 min掃描9個仰角),空間分辨率250 m,最大探測距離250 km,鑒于雷達觀測存在靜錐區(qū)以及電磁波衰減,只選擇距離雷達站50~150 km內的數據進行反演同化。降水觀測資料用于模式檢驗,選取2017年7月9日02—08時貴州區(qū)域內1562個嚴格考核有質量保障的地面自動站質控后的降水觀測數據,時間分辨率為1 h。NCEP數據為GFS全球模式2017年7月8日20時—9日08時的預報結果,時間間隔6 h,空間分辨率0.25°,該數據用于為模擬實驗提供初始場和邊界條件。本文中的時間如無說明均指北京時(BJT)。

    1.2 雷達資料前處理方法

    雷達資料前處理包括雷達資料質控、坐標轉換、雷達組網拼圖,多部雷達組網拼圖時采用最鄰近最大值原則進行(方法略)。

    1.2.1 雷達反射率質控 貴州地處高原東側,地形復雜,雷達回波地物遮擋嚴重,且C波段雷達雨區(qū)衰減較大,因此雷達反射率資料在同化前需要做必要的質量控制。李豐等[11]在現(xiàn)有S波段雷達地物識別方法的基礎上,對C波段雷達地物回波特征進行分析,改進了識別參數的隸屬函數,建立適合C波段多普勒天氣雷達的地物識別方法(MCC方法)。檢驗表明該方法可顯著提高C波段雷達地物回波的識別正確率,并可減少層狀云降水回波的誤判。

    該方法依據地物回波主要出現(xiàn)在較低的仰角、回波徑向速度較小、沒有明顯移動、回波垂直變化大等特征,計算了4個能明顯區(qū)分C波段雷達地物回波和降水回波的因子(回波紋理變率、回波垂直變率、回波庫間變率、中值濾波徑向速度),利用隸屬函數對各個識別因子進行模糊化處理,并對模糊權重求和得到識別指數IAPI 。當IAPI≥0.6時,識別為地物回波,給予消除。在地物回波識別剔除之前,還進行了孤立回波消除,質控算法詳細介紹參見文獻[9]。

    1.2.2 雷達資料坐標轉換 雷達觀測資料為球坐標,需要轉換成與模式相匹配的三維格點坐標。坐標轉換方法采用三線性插值方法,分別沿方位角θ、徑向r、仰角φ3個方向進行插值得到格點上的反射率因子。

    假設笛卡爾坐標中位于(x,y,z)處的格點,經坐標轉換后位于球坐標(r,θ,φ)處,兩個坐標系中心點都位于雷達所在位置,r為徑向距離,θ為方位角,φ為仰角。則(r,θ,φ)點處的反射率因子Z(r,θ,φ)可以用與其相鄰的8個點的觀測值通過三線性插值獲得。

    (1)

    通過式(1)可得點(r,θ,φ)在它的2個相鄰仰角的高仰角錐面上垂直投影點(r,θ,φj)處的反射率Z(r,θ,φj)。同理可得低仰角上的反射率Z(r,θ,φj-1)。插值過程如圖1所示。

    圖1 雙線性插值示意圖

    再用這2個回波值按仰角線性插值(式2),即可得到笛卡爾坐標中網格點上的值Z(x,y,z)。

    (2)

    1.3 同化方法

    目前有多種雷達資料同化方法,張誠忠等[12],徐道生等[13]比較了幾種不同的雷達回波反演技術,發(fā)現(xiàn)利用LAPS系統(tǒng)云分析方案反演云物理量并用Nudging技術調整模式水汽場是一套比較適合華南區(qū)域高分辨率模式的技術路線。

    1.3.1 云微物理量反演方法 LAPS云分析模塊主要分析產品包括三維的云量場、云水及云冰混合比、云及降水類型、結冰嚴重指數、雹混合比、云底及云頂高度等[2]。本文云微物理量反演采用張?zhí)m等[14]提出的方法,利用LAPS云分析系統(tǒng)對云水(qc) 和雨水(qr) 2個微物理量進行反演。首先利用模式背景場的氣壓和溫度計算空氣密度,然后根據雷達回波反演qr,公式如下:

    (3)

    其中,ρ為空氣密度,dbz為雷達反射率因子。

    云水(qc)由雨水(qr)及其下落末速度計算得到,下落末速度計算公式如下:

    (4)

    其中,Vt為下落末速度,P為各層氣壓,Ps為地面氣壓,qr為雨水。

    為了獲得qc,還需要計算1個與水成物降落有關的量ir:

    (5)

    其中,下標1和2表示相鄰2層變量,其余項意義同上。

    最后,根據ir和qr計算出qc,公式如下:

    (6)

    1.3.2 Nudging方法 Nudging 是在模式積分的同化時段δt內,在預報方程中增加一個線性強迫項,該項與模式預報值和實況值之差成正比,其作用是使模式預報逐漸向觀測逼近[15-16]。

    (7)

    W為模式預報變量,∑Fi代表模式中的所有物理過程變率,α>0為張弛逼近系數,W0為對應時刻的觀測值。將方程(1)對時間進行積分,并采用準隱式分步計算方法寫成離散形式如下:

    (8)

    2 模式設計和個例實驗

    2.1 個例簡介

    2017年7月8—9日,受梅雨鋒西段的影響,貴州中南部出現(xiàn)了1次暴雨天氣過程,部分地區(qū)還出現(xiàn)了大暴雨。此次過程以穩(wěn)定的層狀云降水為主,雷達回波以層狀云回波為主。過程主要影響時段:2017年7月8日20時—9日16時。9日02—08時,在鋒面上有明顯的帶狀回波自西北向東南影響貴州,在貴州中南部造成了暴雨到大暴雨天氣,以下對該時段進行數值模擬實驗。

    2.2 模式和實驗設計

    模式采用華南精細預報模式,該模式是基于GRAPES的高分辨率數值預報模式[17]。模式范圍:96~123.36°E,16~31.36°N,水平分辨率3 km,模式其他參數設定同文獻[6]。

    本文擬通過設計2種實驗方案,對2017年7月9日早晨發(fā)生在貴州中南部的暴雨天氣過程進行模擬,分析雷達資料同化對模式前6 h降水預報的影響。設計方案1(Exp1):當天12 UTC NCEP資料積分6 h,并將第6 h預報產品作為模式初始場;設計方案2(Exp2):Exp1背景場,同化18 UTC雷達反射率因子后的分析場。

    2.3 結果分析

    圖2展示了2017年7月9日02—08時累計降水觀測實況(圖2a)及Exp1(圖2b)和Exp2(圖2c)2組模擬實驗預測的降水空間分布情況。從圖2a可以看到,此次過程強降水主要集中在貴州中南部,強降水區(qū)域呈東北—西南向帶狀分布,降水量級以中到大雨為主,局地出現(xiàn)暴雨和大暴雨,暴雨主要分布在黔南州中部和安順市南部。對比實況和控制實驗(Exp1,圖2b)發(fā)現(xiàn),Exp1成功模擬出了這次強降水過程中東北—西南向的強降水雨帶,但是強降水位置較實況偏北且面積比實況更小。例如Exp1暴雨落區(qū)主要位于安順市西北部,而實況的暴雨面積更大且位置更加偏南。對比控制實驗Exp1和對照實驗(Exp2,圖2c)可以發(fā)現(xiàn),同化雷達資料后的實驗Exp2對強降水空間分布做了調整??梢悦黠@看出,Exp1實驗在黔南州中部、安順市南部沒有預測暴雨,存在漏報現(xiàn)象。Exp2實驗在這些區(qū)域預測了多個暴雨點,雖然這些暴雨點面積較實況偏小,但是較Exp1實驗而言更加接近實況。

    圖2 2017年7月9日02—08時實況降水(a) 和實驗模擬6 h降水(b.Exp1;c.Exp2)

    為定量描述2組實驗模擬的降水在雨帶位置偏差和降水強度方面的差異,利用2017年7月9日02—08時的累計降水實況和對應時段的預報降水,分別計算了5~50 mm共6個等級的降水偏差幅度和站數比。降水實況選取貴州區(qū)域1562個質控后的雨量觀測站數據,預報降水為利用線性插值方法將模式輸出的格點降水插值到1562個觀測站上的數據,結果如圖3所示。

    降水偏差幅度用于描述預報降水與實況之間空間位置的偏差程度,計算公式如下:

    (9)

    其中,NA為預報正確的站(次)數、NB為空報站(次)數、NC為漏報站(次)數。NA、NB、NC由表1給出,表中T為降水量級。B為介于0~1之間的數,數值越小說明預測的降水空間分布越接近實況,預報效果越好。

    表1 降水檢驗評定表

    站數比通過計算實況或預測降水中大于某個降水量級T的站數和總降水觀測站數的比乘以100得到。

    圖3a為各降水等級對應的偏差幅度,橫坐標為降水等級,縱坐標為偏差幅度得分。從圖3a中可以看出,未同化雷達的Exp1實驗對于5 mm和10 mm等級降水預報較好,偏差幅度小于0.1。但是當降水量級大于25 mm后,偏差幅度快速增大雨帶位置偏離實況越來越多,在50 mm處,最大約等于0.5。同化雷達后Exp2實驗在大于25 mm的強降水量級中的偏差幅度得分較Exp1大幅減少,全部評分都低于0.1。雨帶位置得到調整,相比Exp1更加接近實況分布。但是,Exp2實驗在較弱的降水量級(5、10 mm)中的偏差幅度得分較Exp1反而升高了,尤其是5 mm對應的評分超過0.2,預報效果較Exp1差。

    圖3b為各降水等級對應的站數比率。從圖中可以看到,Exp1實驗預測的5 mm量級降水站數比較實況偏大。當降水量級大于10 mm以后,Exp1降水站數比低于實況,說明對于較強的降水量級控制實驗降水強度比實況偏弱。對比Exp1和Exp2發(fā)現(xiàn),同化雷達資料后,所有量級的降水站數比都增加了。結合實況的站數比分析發(fā)現(xiàn),當降水大于25 mm后,Exp2的站數比與實況幾乎一致,比Exp1的預測效果更好。而對于5 mm和10 mm的降水,Exp2站數比與實況的偏差更大,存在過度預報。

    圖3 2017年7月9日02—08時分等級降水預報相對實況偏差幅度(a);各等級降水站數相對總雨量站數比率(b)

    為考察同化雷達資料前后模式逐小時降水預報效果差異,計算了7月9日02—08時逐小時降水站數比、偏差幅度、分等級TS評分和相關系數,結果如圖4所示。

    TS評分描述預報的命中情況,計算公式如下:

    (10)

    式中,NA、NB、NC意義與公式9相同,計算方法見表1。

    相關系數反映了預報與實況之間相關關系的密切程度,是一個常用的統(tǒng)計學量,限于篇幅,計算方法此處省略。

    圖4a顯示了5mm降水量級逐小時的站數比,橫坐標為預報時間,縱坐標為站數比。從圖中可以看到,隨著時間推移,實際的降水站數先升高后下降,峰值出現(xiàn)在預報的第2 h,降水在前1~3 h較強,其后逐漸減弱。從Exp1實驗來看,5 mm以上的站數比也呈現(xiàn)先升高后下降的趨勢,峰值出現(xiàn)在第5 h,前3 h Exp1站數較實況小,第4 h后站數高于實況,6 h后逐漸降低趨近實況??傮w而言,由于初始場中沒有足夠的水汽,導致控制實驗(Exp1)在開始的前3 h降水預報偏弱,降水預報較實況存在滯后性,也就是存在spin-up問題。對比Exp2發(fā)現(xiàn),同化雷達后站數比也呈現(xiàn)先上升后下降得趨勢,但是前3 h Exp2比Exp1預報了更多的降水,與實況更接近。而且站數比峰值出現(xiàn)在第4 h,較Exp1提前了1 h。通過吸收雷達資料,調整初始場中的水汽,Exp2縮短了模式spin-up時間。

    相應的,從前6 h內逐小時5 mm量級降水預報相對實況的偏差幅度(圖4b)可以看出,同化雷達后在預報前3 h Exp2對5 mm降水量級的預報偏差幅度小于Exp1,同樣說明吸收雷達資料對緩解模式spin-up問題有積極作用,但這種作用主要集中在前3 h。

    圖4 7月9日02—08時逐小時5 mm量級降水站數比(a);逐小時5 mm量級降水偏差幅度(b);分等級逐小時降水預報TS評分差(c);逐小時預報降水與實況的相關系數(d)

    3 h后,雷達資料同化產生了副作用,例如在預報的第4~6 h,Exp2比Exp1預報了更多的降水(圖4a),比Exp1偏離實況更多,圖4b中4~6 h的偏差幅度對比也能印證這一點。結合圖3b的分析可以看出,偏多的降水主要由10 mm以下的弱降水貢獻,且主要來源于預報的4~6 h。同化雷達能夠彌補模式前6 h強降水預報偏弱的不足,但是產生了更多的弱降水,導致模式對10 mm以下的弱降水預報較實況偏差幅度更大(圖3a)。

    圖4c為利用實況和預報的降水分別檢驗了2組模擬實驗7月9日02—08時之間逐小時累計降水預報分等級TS評分,然后對TS評分進行相減,用來定量考察雷達資料同化對模式短臨降水預報效果的改進,預報降水仍然采用線性插值方法插值到1562個觀測站上與實況比較。圖中虛線為同化了雷達資料的對照實驗Exp2分等級降水TS評分減去未同化雷達資料的控制實驗Exp1對應時段的TS評分,正值表示同化雷達資料后降水預報質量得到提升,負值相反,值越大表示提升或下降越大。

    從圖4c中可以看到,同化雷達資料后模式前6 h逐小時的降水預報在所有降水等級中的評分都有提高。晴雨預報質量開始時刻提升最大隨后迅速下降,中雨以上量級降水預報質量的提升主要集中在4~6 h,呈先上升后下降的趨勢,峰值出現(xiàn)在4~5 h。從TS評分的變化可以看出,雷達資料的同化明顯改善了模式前6 h的降水預報質量,改進主要由強降水貢獻,且主要集中在第4~6 h。

    從逐小時的相關系數(圖4d)來看,吸收雷達資料后Exp2比Exp1與實況的相關性更高,相關系數都隨時間先升高后下降,2組相關系數的差值在第3~6 h更大。總體而言,雷達資料同化改進了模式前6 h的降水預報質量。

    3 結論與討論

    本文在雷達反射率因子質控的基礎上,利用Laps模式云分析系統(tǒng)對雷達資料做了反演,然后利用Nudging方法將雷達反演的云微物理量引入中尺度模式,通過2017年貴州1次強降水天氣過程模擬實驗,得到如下結論:

    ①雷達資料同化能夠改進中尺度模式的降水預報,TS評分檢驗顯示同化雷達資料后模式前6h的降水預報在所有降水等級中的評分都有提高,同化雷達使得降水預報相關系數更高。降水峰值提前1 h,有助于緩解模式spin-up問題。

    ②模式降水預報的改進效果主要由強降水貢獻,最大改進效果集中在第4~6 h。

    ③同化雷達資料后25 mm以上強降水預報站數比更接近實況,落區(qū)偏差幅度更小,降水落區(qū)和強度向實況方向得到調整。

    ④ 10 mm以下弱降水在吸收雷達資料后,站數比較控制預報比實況增加更多,落區(qū)偏差幅度增大,存在預報過量的問題。弱降水預報過量主要集中在4~6 h,而前3 h對降水預報的改進有積極作用。

    本文的研究表明,C波段雷達資料同化,對GRAPES中尺度模式前6 h降水預報有改善作用,但是存在弱降水預報過量的問題,下一步將深入分析其原因并研究改進方法。另外,本次實驗僅僅基于1個強降水個例,今后還需要選取多種類型的降水過程做更廣泛的驗證。

    猜你喜歡
    實況強降水反演
    反演對稱變換在解決平面幾何問題中的應用
    中等數學(2022年5期)2022-08-29 06:07:38
    2020年江淮地區(qū)夏季持續(xù)性強降水過程分析
    一次東移型西南低渦引發(fā)的強降水診斷分析
    鄉(xiāng)村小學的愿望與現(xiàn)實——宜君鄉(xiāng)村教育實況
    當代陜西(2019年10期)2019-06-03 10:12:24
    基于低頻軟約束的疊前AVA稀疏層反演
    基于自適應遺傳算法的CSAMT一維反演
    天舟一號貨運飛船發(fā)射實況掠影
    可愛潮咖們的獨門彩妝實況直播
    四川盆地西南部短時強降水天氣特征分析
    2014 年5 月末西藏強降水天氣分析
    西藏科技(2015年3期)2015-09-26 12:11:10
    无限看片的www在线观看| 成年版毛片免费区| www日本在线高清视频| 日韩国内少妇激情av| 成人高潮视频无遮挡免费网站| 欧美激情久久久久久爽电影| 老熟妇乱子伦视频在线观看| 国产精品久久久久久久电影 | 特级一级黄色大片| 国产成+人综合+亚洲专区| 免费av不卡在线播放| 国产精品98久久久久久宅男小说| 国产毛片a区久久久久| 在线观看免费视频日本深夜| 亚洲电影在线观看av| 午夜福利高清视频| 午夜免费成人在线视频| 人人妻人人澡欧美一区二区| 色噜噜av男人的天堂激情| 久久这里只有精品中国| 国产精品,欧美在线| АⅤ资源中文在线天堂| 偷拍熟女少妇极品色| 90打野战视频偷拍视频| 51午夜福利影视在线观看| 欧美xxxx黑人xx丫x性爽| 精品熟女少妇八av免费久了| e午夜精品久久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 午夜精品一区二区三区免费看| 亚洲午夜理论影院| 超碰成人久久| 久久国产精品影院| ponron亚洲| 日韩人妻高清精品专区| av在线天堂中文字幕| 久久精品综合一区二区三区| 在线观看美女被高潮喷水网站 | 国产亚洲精品av在线| 亚洲欧美日韩东京热| 午夜两性在线视频| 国产乱人视频| 熟妇人妻久久中文字幕3abv| 人人妻人人看人人澡| 亚洲人成网站高清观看| 99久久精品一区二区三区| 可以在线观看毛片的网站| 中文字幕高清在线视频| 在线国产一区二区在线| 亚洲va日本ⅴa欧美va伊人久久| 老司机深夜福利视频在线观看| 国产精品一区二区精品视频观看| 国产精品,欧美在线| 日本一二三区视频观看| 午夜精品一区二区三区免费看| 国产午夜福利久久久久久| 精品国产三级普通话版| 变态另类丝袜制服| 色综合站精品国产| 日日干狠狠操夜夜爽| 青草久久国产| 久久这里只有精品19| 老司机午夜福利在线观看视频| 久久久久久久久免费视频了| 好看av亚洲va欧美ⅴa在| 一级黄色大片毛片| 日本免费一区二区三区高清不卡| 啦啦啦免费观看视频1| 十八禁网站免费在线| 99国产综合亚洲精品| 男插女下体视频免费在线播放| 在线看三级毛片| 国产人伦9x9x在线观看| 欧美日韩国产亚洲二区| 999久久久精品免费观看国产| 亚洲国产精品sss在线观看| 美女高潮喷水抽搐中文字幕| 国产高清有码在线观看视频| 精品熟女少妇八av免费久了| 丰满人妻熟妇乱又伦精品不卡| 在线观看舔阴道视频| 午夜两性在线视频| 国产精品综合久久久久久久免费| 身体一侧抽搐| 啦啦啦免费观看视频1| 日韩欧美三级三区| 精品免费久久久久久久清纯| 午夜精品久久久久久毛片777| 久久久久国产一级毛片高清牌| 欧美大码av| 国产熟女xx| 全区人妻精品视频| 人妻夜夜爽99麻豆av| 国产成人欧美在线观看| 色视频www国产| 亚洲av成人精品一区久久| 亚洲最大成人中文| 黄色 视频免费看| 日韩欧美国产一区二区入口| avwww免费| 国产欧美日韩精品亚洲av| 91在线精品国自产拍蜜月 | 久久精品91无色码中文字幕| 午夜激情欧美在线| 国产人伦9x9x在线观看| 国产成+人综合+亚洲专区| 动漫黄色视频在线观看| www.自偷自拍.com| 91麻豆精品激情在线观看国产| 天天添夜夜摸| 青草久久国产| 国产激情久久老熟女| 看片在线看免费视频| av天堂中文字幕网| 亚洲av成人不卡在线观看播放网| 99久久精品热视频| 色综合站精品国产| 最近最新免费中文字幕在线| 久久这里只有精品19| 欧美大码av| 91av网一区二区| 国产又色又爽无遮挡免费看| 九色成人免费人妻av| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 婷婷丁香在线五月| 午夜福利成人在线免费观看| 脱女人内裤的视频| 757午夜福利合集在线观看| 国产爱豆传媒在线观看| 国产蜜桃级精品一区二区三区| 嫩草影院精品99| 免费搜索国产男女视频| 别揉我奶头~嗯~啊~动态视频| 久久久国产成人精品二区| 日本一本二区三区精品| 美女午夜性视频免费| 成年女人永久免费观看视频| 精品久久久久久久末码| 日本撒尿小便嘘嘘汇集6| 日本一本二区三区精品| 色老头精品视频在线观看| 国产精品久久久久久亚洲av鲁大| 国产私拍福利视频在线观看| 亚洲av第一区精品v没综合| 国产麻豆成人av免费视频| 级片在线观看| 欧美绝顶高潮抽搐喷水| 99久久久亚洲精品蜜臀av| 99久久精品热视频| 国产伦人伦偷精品视频| 熟女人妻精品中文字幕| 两个人看的免费小视频| 老司机深夜福利视频在线观看| 老熟妇仑乱视频hdxx| 一进一出好大好爽视频| 狠狠狠狠99中文字幕| 狂野欧美白嫩少妇大欣赏| 久久热在线av| 老司机在亚洲福利影院| 久久久国产欧美日韩av| 日韩三级视频一区二区三区| 久久久久九九精品影院| 亚洲欧美精品综合久久99| 欧美高清成人免费视频www| 色哟哟哟哟哟哟| 国产欧美日韩一区二区三| а√天堂www在线а√下载| 亚洲人成伊人成综合网2020| 欧美最黄视频在线播放免费| 国产精品一区二区精品视频观看| 午夜影院日韩av| 欧美性猛交黑人性爽| 1024手机看黄色片| 午夜久久久久精精品| 又黄又爽又免费观看的视频| 国产男靠女视频免费网站| 女生性感内裤真人,穿戴方法视频| 黄片大片在线免费观看| 亚洲男人的天堂狠狠| 国产精品99久久久久久久久| 伊人久久大香线蕉亚洲五| 中文字幕人妻丝袜一区二区| 午夜福利在线观看吧| 校园春色视频在线观看| 香蕉久久夜色| 人人妻,人人澡人人爽秒播| АⅤ资源中文在线天堂| 黄片大片在线免费观看| 伦理电影免费视频| 网址你懂的国产日韩在线| 午夜影院日韩av| 一个人免费在线观看的高清视频| 最近在线观看免费完整版| 久久人妻av系列| 又紧又爽又黄一区二区| 国产精品香港三级国产av潘金莲| 男女之事视频高清在线观看| 免费在线观看亚洲国产| 免费一级毛片在线播放高清视频| 国产一区二区三区视频了| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| 色av中文字幕| 老司机深夜福利视频在线观看| 亚洲真实伦在线观看| 90打野战视频偷拍视频| 成人av在线播放网站| 夜夜躁狠狠躁天天躁| 久久精品91蜜桃| 国产精品免费一区二区三区在线| 中国美女看黄片| 国产av一区在线观看免费| 久久久久国产精品人妻aⅴ院| 天天躁狠狠躁夜夜躁狠狠躁| 国内毛片毛片毛片毛片毛片| 热99re8久久精品国产| 丁香六月欧美| 亚洲专区中文字幕在线| 国产一区二区激情短视频| 麻豆国产av国片精品| 久久热在线av| 成人特级黄色片久久久久久久| 午夜免费观看网址| www.www免费av| 99久久成人亚洲精品观看| 国产乱人伦免费视频| 国产精品久久久av美女十八| 一个人免费在线观看电影 | 神马国产精品三级电影在线观看| 美女黄网站色视频| 日本一二三区视频观看| 中文在线观看免费www的网站| 91九色精品人成在线观看| 国产三级中文精品| 久久香蕉国产精品| 精品欧美国产一区二区三| 国产在线精品亚洲第一网站| 小蜜桃在线观看免费完整版高清| 国产精品亚洲一级av第二区| 男人舔女人的私密视频| 国产精品久久久人人做人人爽| 一个人观看的视频www高清免费观看 | 久久天躁狠狠躁夜夜2o2o| 99热这里只有精品一区 | 午夜亚洲福利在线播放| 性欧美人与动物交配| 精品国内亚洲2022精品成人| 12—13女人毛片做爰片一| 亚洲熟妇中文字幕五十中出| 1024手机看黄色片| 国产1区2区3区精品| 伊人久久大香线蕉亚洲五| 久久久精品欧美日韩精品| 欧美黄色淫秽网站| 天天一区二区日本电影三级| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 中文字幕最新亚洲高清| 亚洲精品在线美女| 19禁男女啪啪无遮挡网站| 热99re8久久精品国产| 亚洲无线观看免费| av国产免费在线观看| 亚洲精华国产精华精| 黄色丝袜av网址大全| 成年女人看的毛片在线观看| 天堂动漫精品| 美女 人体艺术 gogo| 国产免费男女视频| 少妇丰满av| 精品久久久久久久末码| 日韩大尺度精品在线看网址| 日韩免费av在线播放| 久久中文字幕一级| 国产激情久久老熟女| 亚洲熟女毛片儿| 哪里可以看免费的av片| 国内精品久久久久久久电影| 超碰成人久久| 久久亚洲精品不卡| 国产亚洲精品久久久com| 99精品在免费线老司机午夜| 男人舔奶头视频| 日本a在线网址| 级片在线观看| 亚洲专区中文字幕在线| 性色av乱码一区二区三区2| 日本在线视频免费播放| 老汉色∧v一级毛片| or卡值多少钱| 色av中文字幕| 国产精品久久久久久久电影 | 俄罗斯特黄特色一大片| 亚洲男人的天堂狠狠| 成年女人看的毛片在线观看| 久久久久久久午夜电影| www国产在线视频色| 在线观看66精品国产| 十八禁网站免费在线| e午夜精品久久久久久久| 亚洲色图av天堂| 波多野结衣高清作品| 免费在线观看成人毛片| 日本熟妇午夜| 丰满人妻一区二区三区视频av | 亚洲五月天丁香| 欧美最黄视频在线播放免费| 久久人人精品亚洲av| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| 神马国产精品三级电影在线观看| 日本a在线网址| 制服人妻中文乱码| 日韩欧美 国产精品| 极品教师在线免费播放| 亚洲国产精品久久男人天堂| 亚洲成人久久性| 少妇裸体淫交视频免费看高清| 久久伊人香网站| 久久久久久人人人人人| 精品国产乱码久久久久久男人| 色尼玛亚洲综合影院| 18禁观看日本| 国产成人精品无人区| 亚洲人成网站高清观看| 国产三级在线视频| 麻豆av在线久日| 中文字幕av在线有码专区| 国产真人三级小视频在线观看| 午夜福利视频1000在线观看| 亚洲国产日韩欧美精品在线观看 | 五月玫瑰六月丁香| 两人在一起打扑克的视频| 免费av不卡在线播放| 亚洲国产欧美人成| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人免费电影在线观看| 午夜福利免费观看在线| 国产精品久久电影中文字幕| 三级国产精品欧美在线观看 | 国产伦一二天堂av在线观看| 男女做爰动态图高潮gif福利片| 久久久国产成人免费| 日本熟妇午夜| 一二三四社区在线视频社区8| 性欧美人与动物交配| 一区二区三区激情视频| 日韩欧美免费精品| 最近最新中文字幕大全电影3| 成人国产一区最新在线观看| 母亲3免费完整高清在线观看| 欧美最黄视频在线播放免费| 中文字幕人成人乱码亚洲影| tocl精华| 一个人看视频在线观看www免费 | 亚洲国产欧美人成| 熟女少妇亚洲综合色aaa.| 黄色 视频免费看| 日日干狠狠操夜夜爽| 精品国产乱子伦一区二区三区| 亚洲国产精品sss在线观看| 毛片女人毛片| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 性欧美人与动物交配| 国产一区二区三区视频了| 国产黄色小视频在线观看| 91老司机精品| 搡老岳熟女国产| 99在线视频只有这里精品首页| 岛国在线观看网站| 国产一区二区激情短视频| 国产精品久久电影中文字幕| 男人舔女人的私密视频| a级毛片在线看网站| 国产成人精品久久二区二区免费| 久久精品国产亚洲av香蕉五月| 国产成人aa在线观看| av福利片在线观看| 亚洲av电影在线进入| 日韩欧美国产在线观看| av中文乱码字幕在线| 黄片大片在线免费观看| 亚洲中文日韩欧美视频| 国产高潮美女av| 日本免费a在线| 国产精品 欧美亚洲| 麻豆成人av在线观看| 老鸭窝网址在线观看| 麻豆一二三区av精品| 国产欧美日韩一区二区三| 欧美色欧美亚洲另类二区| 夜夜看夜夜爽夜夜摸| 性色av乱码一区二区三区2| 欧美最黄视频在线播放免费| 久久精品夜夜夜夜夜久久蜜豆| 狂野欧美白嫩少妇大欣赏| 国产亚洲av高清不卡| 日本免费a在线| 美女被艹到高潮喷水动态| 国产主播在线观看一区二区| 国内精品久久久久久久电影| 欧美三级亚洲精品| 亚洲欧美日韩卡通动漫| 少妇的逼水好多| 美女黄网站色视频| 精品日产1卡2卡| 又粗又爽又猛毛片免费看| 又大又爽又粗| 成人高潮视频无遮挡免费网站| 99在线人妻在线中文字幕| 国内久久婷婷六月综合欲色啪| 亚洲美女黄片视频| 黄色片一级片一级黄色片| 变态另类成人亚洲欧美熟女| or卡值多少钱| 高潮久久久久久久久久久不卡| 少妇丰满av| 99久久精品国产亚洲精品| 亚洲一区二区三区色噜噜| 色播亚洲综合网| 色综合婷婷激情| 国产精品一区二区三区四区久久| 欧美中文综合在线视频| 欧美日本亚洲视频在线播放| 后天国语完整版免费观看| 欧美不卡视频在线免费观看| 亚洲性夜色夜夜综合| 天堂网av新在线| 在线永久观看黄色视频| 欧美中文日本在线观看视频| 91字幕亚洲| tocl精华| 在线观看免费视频日本深夜| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 99在线人妻在线中文字幕| 久久午夜综合久久蜜桃| 少妇丰满av| 又爽又黄无遮挡网站| 日韩精品青青久久久久久| 久久午夜亚洲精品久久| 一个人看的www免费观看视频| 亚洲九九香蕉| 亚洲av中文字字幕乱码综合| 美女高潮喷水抽搐中文字幕| 日本三级黄在线观看| 亚洲成人中文字幕在线播放| 中文字幕精品亚洲无线码一区| 老汉色∧v一级毛片| 观看美女的网站| 99国产精品99久久久久| 国产精品av久久久久免费| 成人国产综合亚洲| 两性夫妻黄色片| 男人舔女人的私密视频| 婷婷丁香在线五月| 欧美黄色片欧美黄色片| 真人做人爱边吃奶动态| 狠狠狠狠99中文字幕| 小说图片视频综合网站| 老鸭窝网址在线观看| 国产黄片美女视频| 99久久无色码亚洲精品果冻| 国产精品乱码一区二三区的特点| 天天添夜夜摸| 在线国产一区二区在线| 色综合站精品国产| 成人欧美大片| 久久午夜亚洲精品久久| 亚洲成av人片在线播放无| 亚洲中文字幕日韩| 免费观看人在逋| 午夜精品一区二区三区免费看| 97超视频在线观看视频| 欧美中文综合在线视频| 日韩欧美一区二区三区在线观看| 级片在线观看| 黄片小视频在线播放| 精品一区二区三区av网在线观看| 亚洲av美国av| 观看免费一级毛片| 国产蜜桃级精品一区二区三区| 伦理电影免费视频| 黄色女人牲交| svipshipincom国产片| 黄色女人牲交| 日本一本二区三区精品| 亚洲熟妇熟女久久| 国产欧美日韩精品一区二区| 天天添夜夜摸| 免费看光身美女| 久久亚洲真实| 午夜两性在线视频| 毛片女人毛片| 日韩人妻高清精品专区| 免费在线观看成人毛片| 最好的美女福利视频网| 国产v大片淫在线免费观看| www日本黄色视频网| 亚洲国产精品久久男人天堂| 国产欧美日韩一区二区三| 最新美女视频免费是黄的| 国产精品一区二区三区四区久久| 91老司机精品| 国产欧美日韩精品亚洲av| 国产伦精品一区二区三区四那| 欧美性猛交黑人性爽| 国产亚洲精品一区二区www| 少妇裸体淫交视频免费看高清| 日韩国内少妇激情av| 亚洲欧美日韩卡通动漫| 又黄又爽又免费观看的视频| 成人午夜高清在线视频| 一个人看的www免费观看视频| www.www免费av| 不卡一级毛片| 热99re8久久精品国产| 2021天堂中文幕一二区在线观| 亚洲欧美精品综合久久99| 九九久久精品国产亚洲av麻豆 | 国产成人精品久久二区二区免费| 国产黄a三级三级三级人| 亚洲在线自拍视频| 一二三四社区在线视频社区8| 人人妻人人看人人澡| 一个人免费在线观看电影 | 黄频高清免费视频| 亚洲午夜精品一区,二区,三区| 亚洲国产精品成人综合色| 国产精品乱码一区二三区的特点| 俄罗斯特黄特色一大片| 91九色精品人成在线观看| 欧美日韩乱码在线| 757午夜福利合集在线观看| 后天国语完整版免费观看| av国产免费在线观看| 国产探花在线观看一区二区| 岛国在线免费视频观看| 一区二区三区高清视频在线| 欧美中文综合在线视频| 亚洲无线在线观看| 国产精品亚洲一级av第二区| 日本黄色片子视频| 国产精品一及| 国产精品99久久久久久久久| 精品欧美国产一区二区三| 啦啦啦免费观看视频1| 一级毛片女人18水好多| 久久久久久久精品吃奶| 亚洲 国产 在线| 国产av一区在线观看免费| 亚洲性夜色夜夜综合| 亚洲电影在线观看av| 老司机午夜福利在线观看视频| 老汉色av国产亚洲站长工具| 亚洲av五月六月丁香网| 精品无人区乱码1区二区| 国产精品久久久人人做人人爽| 9191精品国产免费久久| 午夜福利视频1000在线观看| 国产成人系列免费观看| 蜜桃久久精品国产亚洲av| 午夜免费观看网址| 精品国内亚洲2022精品成人| 亚洲国产欧美人成| 国产成+人综合+亚洲专区| 男女视频在线观看网站免费| 午夜免费激情av| 1000部很黄的大片| 精品无人区乱码1区二区| 91av网一区二区| 深夜精品福利| 国产人伦9x9x在线观看| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 久久精品影院6| 9191精品国产免费久久| 亚洲精品中文字幕一二三四区| 1000部很黄的大片| 日本三级黄在线观看| av在线天堂中文字幕| 国产一级毛片七仙女欲春2| 欧美成人免费av一区二区三区| 久久久久久久午夜电影| 亚洲成人免费电影在线观看| 亚洲欧美精品综合久久99| 琪琪午夜伦伦电影理论片6080| 99视频精品全部免费 在线 | 又粗又爽又猛毛片免费看| 制服人妻中文乱码| 日日干狠狠操夜夜爽| 欧美+亚洲+日韩+国产| 久久热在线av| 免费人成视频x8x8入口观看| av中文乱码字幕在线| 啦啦啦观看免费观看视频高清| 99久久无色码亚洲精品果冻| 人人妻人人澡欧美一区二区| 免费观看精品视频网站| 亚洲五月天丁香| 欧美av亚洲av综合av国产av| 成人一区二区视频在线观看| 午夜福利在线观看吧| 欧美av亚洲av综合av国产av| 天堂av国产一区二区熟女人妻| 桃红色精品国产亚洲av| 欧美日本亚洲视频在线播放| 亚洲 欧美一区二区三区| 又大又爽又粗| 国产真人三级小视频在线观看| 欧美高清成人免费视频www| 亚洲色图 男人天堂 中文字幕| 日本 av在线|