• <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
    免费在线观看成人毛片| 在线观看人妻少妇| a级一级毛片免费在线观看| 国产有黄有色有爽视频| 国产伦精品一区二区三区视频9| 综合色丁香网| 妹子高潮喷水视频| 亚洲无线观看免费| 日韩一本色道免费dvd| 黑人高潮一二区| 毛片一级片免费看久久久久| av线在线观看网站| av天堂中文字幕网| 亚洲美女视频黄频| 能在线免费看毛片的网站| 亚洲国产色片| 国产成人a区在线观看| 国产高清不卡午夜福利| 少妇人妻 视频| 国产av国产精品国产| 欧美xxxx性猛交bbbb| 欧美97在线视频| 九九在线视频观看精品| 亚州av有码| 精品久久久精品久久久| 亚洲人成网站高清观看| 高清黄色对白视频在线免费看 | 在线观看国产h片| av女优亚洲男人天堂| 丝袜喷水一区| av.在线天堂| av免费在线看不卡| 99久久精品国产国产毛片| 极品少妇高潮喷水抽搐| 男女免费视频国产| 国国产精品蜜臀av免费| 日韩一区二区三区影片| 小蜜桃在线观看免费完整版高清| 日日啪夜夜撸| 国产黄片视频在线免费观看| 高清毛片免费看| kizo精华| 久久久久久久久久人人人人人人| 成人免费观看视频高清| 中文字幕亚洲精品专区| 国产大屁股一区二区在线视频| 又爽又黄a免费视频| 精品视频人人做人人爽| 免费av不卡在线播放| 国产一区二区三区av在线| 丰满人妻一区二区三区视频av| 国产视频首页在线观看| 午夜老司机福利剧场| 亚洲av不卡在线观看| 在线观看免费视频网站a站| 日本免费在线观看一区| 中文欧美无线码| 色哟哟·www| 亚洲婷婷狠狠爱综合网| 啦啦啦啦在线视频资源| 成人午夜精彩视频在线观看| 国产精品一区二区在线观看99| 亚洲电影在线观看av| 亚洲精品久久久久久婷婷小说| 3wmmmm亚洲av在线观看| av免费在线看不卡| 综合色丁香网| 国产精品国产三级国产专区5o| 午夜免费观看性视频| 欧美最新免费一区二区三区| 日本色播在线视频| 免费久久久久久久精品成人欧美视频 | 人人妻人人添人人爽欧美一区卜 | 在线免费观看不下载黄p国产| 男人爽女人下面视频在线观看| 美女高潮的动态| av黄色大香蕉| 日韩在线高清观看一区二区三区| 大陆偷拍与自拍| 亚洲精品自拍成人| 十分钟在线观看高清视频www | 少妇人妻久久综合中文| 久久人人爽av亚洲精品天堂 | 制服丝袜香蕉在线| 精品国产三级普通话版| 成人漫画全彩无遮挡| av国产久精品久网站免费入址| 99久久人妻综合| 老司机影院毛片| 十分钟在线观看高清视频www | 超碰97精品在线观看| 一个人看视频在线观看www免费| 欧美日韩视频精品一区| 久久青草综合色| 国产91av在线免费观看| 午夜福利在线观看免费完整高清在| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 尾随美女入室| 亚洲怡红院男人天堂| 能在线免费看毛片的网站| 国产高清三级在线| 日韩欧美 国产精品| 免费av中文字幕在线| 久久人人爽人人片av| 色5月婷婷丁香| 五月伊人婷婷丁香| av不卡在线播放| 大又大粗又爽又黄少妇毛片口| 久久久久久九九精品二区国产| 97在线视频观看| 少妇熟女欧美另类| 亚洲av二区三区四区| 一级毛片 在线播放| 大香蕉97超碰在线| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 亚洲精品国产av成人精品| 男人添女人高潮全过程视频| 中文欧美无线码| 久久 成人 亚洲| 精品人妻一区二区三区麻豆| 欧美精品人与动牲交sv欧美| 成人漫画全彩无遮挡| 在线亚洲精品国产二区图片欧美 | 久久久a久久爽久久v久久| 国产一区二区三区综合在线观看 | 亚洲一级一片aⅴ在线观看| 免费av不卡在线播放| 国产在线男女| 亚洲成色77777| 久久国产乱子免费精品| 又粗又硬又长又爽又黄的视频| 日韩大片免费观看网站| 五月玫瑰六月丁香| 国产在视频线精品| 国产日韩欧美亚洲二区| 成人亚洲欧美一区二区av| 亚洲伊人久久精品综合| 亚洲高清免费不卡视频| 亚洲不卡免费看| 国产在线免费精品| 亚洲av日韩在线播放| 99久久中文字幕三级久久日本| 高清不卡的av网站| 国产精品久久久久久久久免| 久久av网站| 91aial.com中文字幕在线观看| 夫妻午夜视频| 免费观看无遮挡的男女| 欧美极品一区二区三区四区| 看免费成人av毛片| 26uuu在线亚洲综合色| 91在线精品国自产拍蜜月| 亚洲精品中文字幕在线视频 | 亚洲精品国产av成人精品| 在线亚洲精品国产二区图片欧美 | 久久6这里有精品| 免费在线观看成人毛片| 久久久久久久亚洲中文字幕| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 美女福利国产在线 | 久久久久久久久久久免费av| 日韩制服骚丝袜av| 男人和女人高潮做爰伦理| 乱码一卡2卡4卡精品| 国产熟女欧美一区二区| 精品午夜福利在线看| 色婷婷久久久亚洲欧美| 久热这里只有精品99| 亚洲欧美精品自产自拍| 亚洲怡红院男人天堂| 免费观看的影片在线观看| 成人综合一区亚洲| 国产永久视频网站| 国产视频首页在线观看| 亚洲国产色片| 日韩 亚洲 欧美在线| 亚洲三级黄色毛片| 1000部很黄的大片| 在线观看三级黄色| 国产在视频线精品| 韩国av在线不卡| 国产在线男女| 亚洲成色77777| 欧美bdsm另类| 高清不卡的av网站| 亚洲av不卡在线观看| 久久亚洲国产成人精品v| 精品人妻偷拍中文字幕| 男女边吃奶边做爰视频| 精品人妻熟女av久视频| 久久这里有精品视频免费| 午夜福利在线观看免费完整高清在| 国产亚洲91精品色在线| videossex国产| 久久精品国产亚洲av天美| 一个人看视频在线观看www免费| 欧美日韩视频精品一区| 亚洲av福利一区| 18禁裸乳无遮挡免费网站照片| 纯流量卡能插随身wifi吗| 欧美老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 黑人猛操日本美女一级片| 精品人妻熟女av久视频| 日韩电影二区| 交换朋友夫妻互换小说| 精品久久国产蜜桃| 国产精品一区www在线观看| 中文字幕久久专区| 欧美xxⅹ黑人| 亚洲国产欧美在线一区| 久久久色成人| 人妻 亚洲 视频| 久久久亚洲精品成人影院| 丰满乱子伦码专区| 国产在线一区二区三区精| 日韩av在线免费看完整版不卡| 久久国产亚洲av麻豆专区| 久久这里有精品视频免费| 日本wwww免费看| 欧美xxxx性猛交bbbb| av天堂中文字幕网| 久久ye,这里只有精品| 99久久人妻综合| 国产亚洲5aaaaa淫片| 夜夜骑夜夜射夜夜干| 久久久久久久久久成人| 欧美国产精品一级二级三级 | 久久精品国产亚洲av涩爱| 97超视频在线观看视频| 九九爱精品视频在线观看| 91午夜精品亚洲一区二区三区| 久久韩国三级中文字幕| 国产成人一区二区在线| 黑人猛操日本美女一级片| 精品人妻一区二区三区麻豆| 26uuu在线亚洲综合色| 最近2019中文字幕mv第一页| 亚洲欧美一区二区三区黑人 | 亚洲图色成人| 深夜a级毛片| 国产精品熟女久久久久浪| 欧美激情极品国产一区二区三区 | 亚洲中文av在线| 久久人人爽人人爽人人片va| 国产爽快片一区二区三区| 天天躁夜夜躁狠狠久久av| 天天躁日日操中文字幕| 久久久久国产网址| 欧美日韩视频高清一区二区三区二| 欧美一级a爱片免费观看看| 能在线免费看毛片的网站| av国产精品久久久久影院| av专区在线播放| 久久久久久久久大av| 久久99热6这里只有精品| 国产免费福利视频在线观看| 成人无遮挡网站| 亚洲无线观看免费| 国产亚洲最大av| 国产高清三级在线| 日日撸夜夜添| 久久久久久久久久久丰满| 亚洲美女黄色视频免费看| 色婷婷久久久亚洲欧美| 欧美三级亚洲精品| 卡戴珊不雅视频在线播放| 91精品一卡2卡3卡4卡| 各种免费的搞黄视频| 国产精品蜜桃在线观看| 深夜a级毛片| 狂野欧美激情性bbbbbb| 成人毛片a级毛片在线播放| 18禁在线无遮挡免费观看视频| 国产毛片在线视频| 搡老乐熟女国产| 老女人水多毛片| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 1000部很黄的大片| 在线 av 中文字幕| 国产免费福利视频在线观看| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 99九九线精品视频在线观看视频| 久久久欧美国产精品| 自拍欧美九色日韩亚洲蝌蚪91 | 在线精品无人区一区二区三 | 女人久久www免费人成看片| 狂野欧美激情性xxxx在线观看| 日韩视频在线欧美| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 视频区图区小说| 亚洲av综合色区一区| 亚洲一区二区三区欧美精品| 下体分泌物呈黄色| 身体一侧抽搐| 国产精品熟女久久久久浪| 成人亚洲欧美一区二区av| 亚洲,一卡二卡三卡| 亚洲无线观看免费| 深爱激情五月婷婷| 在线播放无遮挡| 亚洲国产欧美人成| 成人一区二区视频在线观看| 九色成人免费人妻av| av在线app专区| 视频区图区小说| 免费看日本二区| av在线观看视频网站免费| 国产成人aa在线观看| 日本黄大片高清| 丰满少妇做爰视频| 久久女婷五月综合色啪小说| 91精品国产国语对白视频| 亚洲无线观看免费| 欧美日韩视频高清一区二区三区二| 国产午夜精品一二区理论片| 欧美国产精品一级二级三级 | 国产永久视频网站| 狂野欧美激情性bbbbbb| 午夜日本视频在线| 亚洲欧美日韩另类电影网站 | 美女cb高潮喷水在线观看| av天堂中文字幕网| 97热精品久久久久久| 三级经典国产精品| 欧美成人a在线观看| 黄色一级大片看看| 国产精品久久久久久av不卡| h日本视频在线播放| 国产高清有码在线观看视频| 九草在线视频观看| 天美传媒精品一区二区| 亚洲在久久综合| 一级毛片久久久久久久久女| 国产精品国产av在线观看| 一级片'在线观看视频| 一级爰片在线观看| 亚洲国产欧美在线一区| 精品视频人人做人人爽| 一级毛片电影观看| 日韩在线高清观看一区二区三区| 美女xxoo啪啪120秒动态图| 国产精品不卡视频一区二区| 舔av片在线| 亚洲精华国产精华液的使用体验| 在线观看国产h片| 97超碰精品成人国产| 国产亚洲最大av| 欧美变态另类bdsm刘玥| 伦理电影免费视频| 久热这里只有精品99| 日韩伦理黄色片| 欧美xxxx性猛交bbbb| .国产精品久久| 国产精品爽爽va在线观看网站| 极品教师在线视频| 欧美区成人在线视频| 国产视频首页在线观看| 99热这里只有是精品50| 美女主播在线视频| 中文字幕精品免费在线观看视频 | 欧美日韩一区二区视频在线观看视频在线| 亚洲国产欧美在线一区| 干丝袜人妻中文字幕| 人人妻人人澡人人爽人人夜夜| 欧美精品国产亚洲| 亚洲精品第二区| 大香蕉久久网| 日韩三级伦理在线观看| 日本黄色片子视频| av.在线天堂| 色视频在线一区二区三区| 国产免费又黄又爽又色| 亚洲精品乱码久久久久久按摩| 大片免费播放器 马上看| 久久久精品免费免费高清| 亚洲欧洲国产日韩| 在线免费观看不下载黄p国产| 99精国产麻豆久久婷婷| 中国国产av一级| 亚洲国产毛片av蜜桃av| 欧美日韩在线观看h| 精品人妻熟女av久视频| 国产精品一区二区性色av| 毛片一级片免费看久久久久| 成人午夜精彩视频在线观看| 99久久中文字幕三级久久日本| 亚洲精品456在线播放app| 亚洲欧美一区二区三区黑人 | 午夜日本视频在线| 最近手机中文字幕大全| 99国产精品免费福利视频| 大片电影免费在线观看免费| 久热久热在线精品观看| 精品人妻偷拍中文字幕| 夜夜骑夜夜射夜夜干| 人妻制服诱惑在线中文字幕| 一级毛片电影观看| 亚洲人成网站高清观看| 国产伦精品一区二区三区视频9| 久久久久久人妻| 亚洲综合色惰| 日韩电影二区| 18禁在线播放成人免费| 久久久久国产网址| 精品久久久精品久久久| 交换朋友夫妻互换小说| 色视频www国产| 精品久久久久久久久亚洲| 国产爽快片一区二区三区| 色婷婷av一区二区三区视频| 在线观看一区二区三区激情| 婷婷色麻豆天堂久久| 秋霞在线观看毛片| 欧美日韩精品成人综合77777| 亚洲美女搞黄在线观看| 久久久久视频综合| 久热这里只有精品99| 久久影院123| 亚洲国产日韩一区二区| 国产成人免费无遮挡视频| 99精国产麻豆久久婷婷| 国产精品一区www在线观看| 亚洲精品久久久久久婷婷小说| 一级av片app| 少妇人妻 视频| 国产高清三级在线| 91久久精品国产一区二区三区| 国产毛片在线视频| 大片免费播放器 马上看| 欧美三级亚洲精品| 午夜激情福利司机影院| 亚洲精品久久久久久婷婷小说| 免费看日本二区| 一区二区av电影网| 99热这里只有是精品50| 国产精品一区二区在线不卡| 日本黄色日本黄色录像| 丰满少妇做爰视频| 中文字幕av成人在线电影| 啦啦啦啦在线视频资源| www.色视频.com| 亚洲三级黄色毛片| 超碰av人人做人人爽久久| 免费大片黄手机在线观看| 毛片一级片免费看久久久久| 国产精品一区二区性色av| 免费在线观看成人毛片| 日韩不卡一区二区三区视频在线| 夜夜爽夜夜爽视频| 搡老乐熟女国产| 欧美激情极品国产一区二区三区 | 日产精品乱码卡一卡2卡三| 国产精品伦人一区二区| 99久久中文字幕三级久久日本| 在线观看免费视频网站a站| 深爱激情五月婷婷| 老司机影院毛片| 亚洲精品日本国产第一区| 日韩欧美 国产精品| 久久精品久久久久久久性| 搡女人真爽免费视频火全软件| videossex国产| 成年美女黄网站色视频大全免费 | 夜夜骑夜夜射夜夜干| 中文字幕av成人在线电影| 国产精品伦人一区二区| 丰满人妻一区二区三区视频av| a级一级毛片免费在线观看| 亚洲精品久久午夜乱码| 插逼视频在线观看| 最近的中文字幕免费完整| 观看av在线不卡| 观看免费一级毛片| 亚洲精品国产成人久久av| 岛国毛片在线播放| 一本一本综合久久| 自拍偷自拍亚洲精品老妇| 九九爱精品视频在线观看| 午夜激情久久久久久久| 人妻一区二区av| 欧美成人一区二区免费高清观看| 日韩av不卡免费在线播放| 晚上一个人看的免费电影| 丰满人妻一区二区三区视频av| 国产成人午夜福利电影在线观看| 毛片一级片免费看久久久久| 欧美精品人与动牲交sv欧美| 国产精品蜜桃在线观看| 中文精品一卡2卡3卡4更新| 国产精品不卡视频一区二区| 欧美xxxx黑人xx丫x性爽| 1024视频免费在线观看| 丝瓜视频免费看黄片| 亚洲国产中文字幕在线视频| 日本黄色日本黄色录像| 老司机亚洲免费影院| 香蕉国产在线看| 在线观看人妻少妇| 黑人欧美特级aaaaaa片| 欧美激情高清一区二区三区| 一区二区三区激情视频| 十分钟在线观看高清视频www| 另类亚洲欧美激情| 国产免费一区二区三区四区乱码| 免费看十八禁软件| 午夜免费成人在线视频| 国产视频首页在线观看| 黑丝袜美女国产一区| 又紧又爽又黄一区二区| 精品久久久久久电影网| 成人午夜精彩视频在线观看| 国产成人精品久久二区二区免费| 久久免费观看电影| 亚洲图色成人| 国产人伦9x9x在线观看| 大陆偷拍与自拍| 国产视频一区二区在线看| 精品久久久久久电影网| 亚洲视频免费观看视频| 99九九在线精品视频| 又大又爽又粗| 国产亚洲精品第一综合不卡| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 国产视频一区二区在线看| 色婷婷av一区二区三区视频| 久久亚洲国产成人精品v| 91九色精品人成在线观看| 1024视频免费在线观看| 女人高潮潮喷娇喘18禁视频| 久久午夜综合久久蜜桃| 一级毛片黄色毛片免费观看视频| 黄网站色视频无遮挡免费观看| 亚洲免费av在线视频| 欧美黄色片欧美黄色片| 免费在线观看视频国产中文字幕亚洲 | 日本av免费视频播放| 亚洲国产毛片av蜜桃av| 免费在线观看日本一区| 纵有疾风起免费观看全集完整版| 中文字幕亚洲精品专区| svipshipincom国产片| 久久女婷五月综合色啪小说| 日韩一区二区三区影片| 各种免费的搞黄视频| 久久久久精品人妻al黑| 老司机在亚洲福利影院| 午夜老司机福利片| 亚洲精品一二三| 婷婷色综合www| 香蕉国产在线看| 精品一区二区三卡| 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| av线在线观看网站| 久久热在线av| 黄频高清免费视频| 亚洲国产成人一精品久久久| 国产精品免费大片| 免费高清在线观看日韩| 激情视频va一区二区三区| xxx大片免费视频| 国产精品九九99| 久久国产精品影院| 国产97色在线日韩免费| 久久久久精品人妻al黑| 精品福利观看| 亚洲av男天堂| 91精品三级在线观看| 国产精品久久久久久精品电影小说| 国产成人av教育| 欧美激情 高清一区二区三区| 免费在线观看视频国产中文字幕亚洲 | 操美女的视频在线观看| 在线精品无人区一区二区三| a 毛片基地| 色婷婷av一区二区三区视频| 在线精品无人区一区二区三| a 毛片基地| 亚洲av男天堂| 国产在线视频一区二区| 亚洲欧美精品综合一区二区三区| 亚洲av男天堂| 亚洲国产日韩一区二区| 97精品久久久久久久久久精品| 七月丁香在线播放| 十八禁网站网址无遮挡| 亚洲欧美中文字幕日韩二区| 色婷婷av一区二区三区视频| 一二三四社区在线视频社区8| 777久久人妻少妇嫩草av网站| e午夜精品久久久久久久| 十八禁网站网址无遮挡| 国产精品免费视频内射| 亚洲成国产人片在线观看| 国产在线免费精品| 国产日韩欧美视频二区| 啦啦啦在线观看免费高清www| 女性生殖器流出的白浆| 韩国高清视频一区二区三区| 一边亲一边摸免费视频| 99国产综合亚洲精品| 久久久久久人人人人人| 制服人妻中文乱码| 满18在线观看网站| 十八禁高潮呻吟视频| av国产精品久久久久影院| 国产黄色免费在线视频|