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

    克拉瑪依大風(fēng)數(shù)值預(yù)報(bào)中的重力波拖曳方案應(yīng)用研究

    2022-04-02 01:01:31張涵斌史永強(qiáng)崔麗娜楊程志
    大氣科學(xué)學(xué)報(bào) 2022年1期

    張涵斌 史永強(qiáng) 崔麗娜 楊程志

    摘要 新疆克拉瑪依地區(qū)位于背風(fēng)坡,長(zhǎng)期受翻山氣流形成的局地大風(fēng)影響。該地區(qū)日常大風(fēng)預(yù)報(bào)業(yè)務(wù)主要依托克拉瑪依氣象局的精細(xì)化數(shù)值預(yù)報(bào)系統(tǒng),為進(jìn)一步完善該系統(tǒng)的預(yù)報(bào)效果,本文開展了重力波拖曳參數(shù)化方案應(yīng)用試驗(yàn),對(duì)比分析了三組試驗(yàn)方案:不開啟重力波拖曳、只開啟外層區(qū)域重力波拖曳以及兩重嵌套設(shè)置下均開啟重力波拖曳的結(jié)果。試驗(yàn)結(jié)果表明:在克拉瑪依大風(fēng)預(yù)報(bào)個(gè)例中,對(duì)于山區(qū)地形背風(fēng)坡處的克拉瑪依站,不啟動(dòng)重力波拖曳會(huì)導(dǎo)致風(fēng)速偏大和起風(fēng)時(shí)間較早;開啟外層重力波拖曳方案可以延遲克拉瑪依站起風(fēng)時(shí)間,兩重區(qū)域均開啟可以進(jìn)一步延遲該站的起風(fēng)時(shí)間,使起風(fēng)時(shí)間更接近實(shí)況;對(duì)于非山區(qū)地形背風(fēng)坡處的站點(diǎn),在內(nèi)層開啟重力波拖曳會(huì)導(dǎo)致一定的負(fù)效果;對(duì)內(nèi)層區(qū)域的統(tǒng)計(jì)檢驗(yàn)結(jié)果表明不開啟重力波拖曳,會(huì)使地形復(fù)雜區(qū)域風(fēng)場(chǎng)預(yù)報(bào)存在一定的正偏差,在外層開啟重力波拖曳會(huì)使正偏差有一定改善,而兩重嵌套均開啟重力波拖曳方案會(huì)導(dǎo)致非地形陡峭區(qū)域的負(fù)偏差增加。綜合全區(qū)域內(nèi)大風(fēng)預(yù)報(bào)檢驗(yàn)結(jié)果表明,只在外層開啟重力波拖曳會(huì)獲得最小的誤差和較好的預(yù)報(bào)效果。

    關(guān)鍵詞重力波拖曳;數(shù)值模式;參數(shù)化方案;大風(fēng)預(yù)報(bào);克拉瑪依

    隨著現(xiàn)代高分辨數(shù)值預(yù)報(bào)技術(shù)的發(fā)展,模式中開始越來(lái)越多地考慮地形作用(王宇和黃思訓(xùn),2018),該作用主要通過(guò)次網(wǎng)格參數(shù)化方案來(lái)實(shí)現(xiàn),主要有“包絡(luò)地形”方案(Wallace et al.,1983;Tibaldi,1986)以及“重力波拖曳”方案(Gravity Wave Drag,GWD;Phillips,1984;Palmer et al.,1986)。對(duì)于“重力波拖曳”參數(shù)化方案,其理論是當(dāng)氣流經(jīng)過(guò)地形陡峭區(qū)時(shí)可能激發(fā)出地形重力波,從而產(chǎn)生水平動(dòng)量耗散,數(shù)值模式中需要通過(guò)次網(wǎng)格參數(shù)化方案來(lái)描述這種地形引起的動(dòng)量耗散過(guò)程,以使風(fēng)場(chǎng)的模擬更加準(zhǔn)確(McFarlane,1987;Baines and Manins,1989;Miller et al.,1989;程胡華,2017)。

    Palmer et al.(1986)首次提出了歐洲中期預(yù)報(bào)中心全球模式的重力波拖曳方案,該方案主要描述了地形陡峭處的重力波激發(fā)機(jī)制,隨后許多學(xué)者對(duì)重力波拖曳方案在模式中的應(yīng)用進(jìn)行了研究。Baines and Palmer(1990)利用非流線形動(dòng)力學(xué)理論對(duì)重力波拖曳方案進(jìn)行了進(jìn)一步完善;Boer et al.(1984)首次實(shí)現(xiàn)了將重力波拖曳參數(shù)化引入到大尺度模式中;Kim and Arakawa(1995)對(duì)重力波拖曳方案進(jìn)一步發(fā)展(KA95方案),并引入到WRF中尺度模式中,該方案主要描述了因地形方差引起的大氣高層和低層重力波破碎將波能傳播到下層并引起拖曳力增強(qiáng)的過(guò)程;Lott and Miller(1997)對(duì)原KA95方案進(jìn)行了進(jìn)一步改進(jìn),在重力波破碎模型的基礎(chǔ)上引入了阻塞拖曳作用,并進(jìn)行了一系列數(shù)值試驗(yàn),結(jié)果表明新方案的預(yù)報(bào)技巧優(yōu)于原方案。此后許多研究也表明,模式通過(guò)引入重力波拖曳方案可以更好地描述大尺度山脈對(duì)對(duì)流層動(dòng)量收支的影響,從而改進(jìn)預(yù)報(bào)效果(Kim and Doyle,2005;Hong et al.,2008;Sandu et al.,2013;Choi and Hong,2015)。

    國(guó)內(nèi)也開展了許多地形重力波拖曳參數(shù)化相關(guān)的研究,并闡述了重力波參數(shù)化方案的作用。周天軍和錢永甫(1996)指出模式地形對(duì)于預(yù)報(bào)質(zhì)量十分關(guān)鍵,需要在數(shù)值預(yù)報(bào)中考慮地形效應(yīng);錢永甫(2000)、劉華強(qiáng)和錢永甫(2001)在氣候模式中針對(duì)重力波拖曳開展了試驗(yàn),研究結(jié)果表明區(qū)域氣候模擬中考慮地形重力波拖曳可以有效改善模擬結(jié)果;徐國(guó)強(qiáng)等(2010)通過(guò)全球天氣模式研究了一次降水個(gè)例,指出引入地形重力波拖曳過(guò)程可以有效改善大氣風(fēng)場(chǎng),對(duì)降水預(yù)報(bào)的落區(qū)有所改善。綜合而言,國(guó)內(nèi)針對(duì)重力波拖曳方案在中尺度數(shù)值模式中的研究較少,尚有待進(jìn)一步開展。

    新疆克拉瑪依地區(qū)地處準(zhǔn)噶爾盆地西北邊,西有加依爾山,中部、東部為平坦戈壁,受特殊地理位置影響,大風(fēng)盛行(辛渝等,2015),其中中部及南部主要受背風(fēng)坡的強(qiáng)下坡風(fēng)影響,春季大風(fēng)較多,最大瞬時(shí)風(fēng)速可達(dá)42.2 m/s。許多學(xué)者針對(duì)大地形影響下得氣流過(guò)山及伴隨的強(qiáng)下坡風(fēng)進(jìn)行了研究,指出由于背風(fēng)波破碎而出現(xiàn)的湍流混合使得來(lái)自底層的重力內(nèi)波在波破碎區(qū)產(chǎn)生反射進(jìn)而導(dǎo)致下坡風(fēng)暴的形成(Klemp and Lilly,1975;Lilly,1978;Clark and Peltier,1984;Saito,1993);孫東霞等(2008)指出克拉瑪依的背風(fēng)坡地形是導(dǎo)致該地區(qū)強(qiáng)大風(fēng)盛行的關(guān)鍵因素;盧冰等(2014)利用高分辨率WRF模式,重點(diǎn)闡述了克拉瑪依地區(qū)背風(fēng)坡風(fēng)暴產(chǎn)生的機(jī)理,對(duì)冬季穩(wěn)定層結(jié)及冷空氣背景下重力波破碎導(dǎo)致的背風(fēng)坡大風(fēng)進(jìn)行了細(xì)致分析,結(jié)果表明重力波效應(yīng)是克拉瑪依地區(qū)大風(fēng)預(yù)報(bào)的關(guān)鍵。

    目前克拉瑪依市氣象局業(yè)務(wù)運(yùn)行數(shù)值天氣預(yù)報(bào)系統(tǒng),為克拉瑪依地區(qū)大風(fēng)天氣的預(yù)報(bào)提供有效的業(yè)務(wù)支撐,但目前對(duì)克拉瑪依地區(qū)大風(fēng)模擬的強(qiáng)度和發(fā)生時(shí)間有待改進(jìn),本研究基于克拉瑪依業(yè)務(wù)數(shù)值預(yù)報(bào)系統(tǒng),探索了重力波拖曳方案在克拉瑪依大風(fēng)預(yù)報(bào)的應(yīng)用效果。本研究不僅對(duì)克拉瑪依當(dāng)?shù)卮箫L(fēng)預(yù)報(bào)的發(fā)展具有重要意義,也可為業(yè)務(wù)數(shù)值預(yù)報(bào)系統(tǒng)中大風(fēng)預(yù)報(bào)改進(jìn)提供新方法新思路,具有良好的應(yīng)用前景。

    1 資料和方法

    1.1 克拉瑪依精細(xì)化區(qū)域數(shù)值預(yù)報(bào)系統(tǒng)簡(jiǎn)介

    克拉瑪依市氣象局發(fā)展的精細(xì)化區(qū)域數(shù)值預(yù)報(bào)系統(tǒng)基于美國(guó)環(huán)境預(yù)測(cè)中心(NCEP)和美國(guó)國(guó)家大氣研究中心(NCAR)共同開發(fā)的非靜力中尺度模式WRF V4.2版構(gòu)建,模式區(qū)域設(shè)置為兩重嵌套,其中外層D01水平分辨率15 km,模擬區(qū)域范圍為中亞及新疆地區(qū),共238×179個(gè)格點(diǎn),垂直方向57層,覆蓋新疆及中亞部分區(qū)域,兩層區(qū)域單向嵌套無(wú)反饋;內(nèi)層D02水平分辨率3 km,共295×200個(gè)格點(diǎn),主要覆蓋新疆北部區(qū)域(圖1)。該系統(tǒng)每天從00時(shí)、12時(shí)(世界時(shí),下同)起報(bào)兩次,時(shí)間步長(zhǎng)90 s,預(yù)報(bào)時(shí)效36 h,輸出間隔1 h。系統(tǒng)采用NCEP全球預(yù)報(bào)場(chǎng)GFS資料(水平分辨率0.5°×0.5°,間隔6 h)作為初值驅(qū)動(dòng)場(chǎng)和側(cè)邊界條件,物理過(guò)程參數(shù)化方案選取為:微物理過(guò)程采用Thompson方案(Thompson et al.,2004),邊界層采用YSU方案(Hong et al.,2006),短波輻射和長(zhǎng)波輻射采用RRTMG方案(Bae et al.,2016),D01積云方案采用New Tiedtke方案(Zhang et al.,2011),D02積云方案關(guān)閉。

    1.2 重力波拖曳方案介紹

    WRF次網(wǎng)格地形重力波拖曳參數(shù)化方案中,包含重力波破碎和阻塞拖曳兩種作用(Kim and Doyle,2005),即在“gwd_opt=1”設(shè)置下,可開啟兩種參數(shù)化方案。

    1.2.1 重力波破碎

    WRF模式中的KA95方案不僅考慮了高層波破碎對(duì)重力波拖曳的影響,也考慮了因低層波破碎對(duì)下游重力波拖曳的影響:翻越地形陡峭區(qū)域的氣流在參考層會(huì)產(chǎn)生相應(yīng)的應(yīng)力τ0,并激發(fā)垂直傳播的重力波,重力波在不穩(wěn)定環(huán)境中波破碎,此時(shí)拖曳力會(huì)受到一定損失,剩余的拖曳應(yīng)力τ將在波破碎后繼續(xù)擴(kuò)散,使得下游也存在一定的拖曳應(yīng)力。依據(jù)Kim and Arakawa(1995),參考層上的GWD為

    其中:ρ為密度;N是Brunt-Vaisala頻率;U是水平風(fēng)速;Δx為模式的水平尺度;下標(biāo)0代表是參考層,也就是重力波產(chǎn)生層;E為參考層上拖曳力的增強(qiáng)因子,用來(lái)增加受到低層波破碎或是抑制作用導(dǎo)致的拖曳力,主要受地形的非對(duì)稱以及臨界Froude數(shù)所影響;m代表網(wǎng)格內(nèi)山的數(shù)量,G是漸進(jìn)函數(shù)用來(lái)斷定氣流阻塞與否,因?yàn)榻?jīng)過(guò)修正,所以以m′、G′表示(鐘水新等,2014)。

    當(dāng)重力波不滿足飽和假設(shè)時(shí),波能在上傳過(guò)程中不損失,即上下層的GWD不變(τi=τi+1),據(jù)此獲得上層波振幅(Palmer et al.,1986),并基于計(jì)算的上層波振幅求取局地修正后的Richardson數(shù)Ri。依據(jù)Lindzen(1981)的飽和假設(shè),當(dāng)滿足Ri<Ric(飽和臨界點(diǎn)的Richardson數(shù))時(shí),重力波在該層產(chǎn)生波破碎,部分能量損失使得τ減小,直至恢復(fù)到飽和臨界點(diǎn)。當(dāng)未達(dá)飽和條件時(shí),重力波在該層不破碎,該層τ維持不變。

    臨界振幅hd的計(jì)算通過(guò)將Ri代入下式獲得:

    求得臨界振幅后,計(jì)算該層剩余GWD:

    獲得GWD的垂直方向上所有層次(往上直到τi=0或模式層頂)的分布后,代入式(4)、(5)完成地形重力波拖曳對(duì)風(fēng)場(chǎng)傾向的計(jì)算(鐘水新等,2014)。

    1.2.2 阻塞拖曳

    拖曳參數(shù)化中的阻塞拖曳(Blocked-Layer Drag,BLD)是對(duì)KA95方案進(jìn)行進(jìn)一步擴(kuò)展,公式基本上遵循Lott and Miller(1997)以及Scinocca and McFarlane(2000)的參數(shù)化研究,基于尺度分析獲得,即如下拖曳形式:

    其中:Δ2x為網(wǎng)格面積;Cd為一階體積阻力系數(shù);Δ⊥x是垂直于風(fēng)向的網(wǎng)格長(zhǎng)度;L⊥x是垂直于風(fēng)向的次網(wǎng)格地形的寬度,近似為臨界地形高度測(cè)得的地形寬度;hB為阻塞層高度,定義為:

    地形的各向異性通過(guò)hB考慮,hB的計(jì)算包括Froude數(shù)Fr0,F(xiàn)rc(≈0.8)為臨界Froude數(shù),當(dāng)Fr0>Frc時(shí),該式應(yīng)用到最低模式面以及地面層,并往上呈線性遞減。

    圖2給出了克拉瑪依數(shù)值模式D02范圍內(nèi)的地形分布以及WRF4.2中重力波拖曳參數(shù)化方案引起的地形拖曳應(yīng)力作用,由圖2a可以看出加依爾山-天山一帶為地形陡峭區(qū),其中克拉瑪依位于加依爾山和準(zhǔn)噶爾盆地交界處(46°N一線),由圖2b可以看出次網(wǎng)格地形方差較大區(qū)域即為山區(qū)陡峭地形處,而地形重力波參數(shù)化方案引起的拖曳應(yīng)力與地形大值區(qū)對(duì)應(yīng)較好。

    1.3 試驗(yàn)方案設(shè)計(jì)

    為了驗(yàn)證重力波拖曳方案在克拉瑪依精細(xì)化模式系統(tǒng)中的應(yīng)用效果,基于重力波拖曳參數(shù)化方案開展了大風(fēng)預(yù)報(bào)試驗(yàn),重點(diǎn)針對(duì)重力波拖曳參數(shù)化方案在3 km區(qū)域內(nèi)的表現(xiàn)進(jìn)行分析。系統(tǒng)模式分辨率及范圍設(shè)置均與1.1中介紹的一致。

    設(shè)計(jì)了三種試驗(yàn)方案,方案一在D01和D02區(qū)域內(nèi)關(guān)閉了GWD參數(shù)化方案,其余設(shè)置與1.1中介紹的方案一致,稱為GWD00方案;方案二在D01區(qū)域開啟GWD參數(shù)化方案,D02區(qū)域關(guān)閉GWD方案,稱為GWD10方案,方案三在D01和D02范圍內(nèi)均開啟重力波拖曳,稱為GWD11方案,如表1所示。開展了2021年2月1個(gè)月的連續(xù)試驗(yàn),其中以2月24日—25日發(fā)生在克拉瑪依的大風(fēng)天氣過(guò)程為重點(diǎn)研究個(gè)例。

    2 試驗(yàn)結(jié)果分析

    2.1 大風(fēng)個(gè)例分析

    2021年2月24日,受極渦分裂短波東移影響,克拉瑪依城區(qū)發(fā)生了一次大風(fēng)天氣,25日06時(shí),在冷空氣南下影響下,克拉瑪依城區(qū)起風(fēng),逐小時(shí)最大風(fēng)速可達(dá)21.3 m/s,已達(dá)到大風(fēng)級(jí)別(近地面層平均風(fēng)速17.2~20.7 m/s)。圖3給出了三種重力波拖曳方案配置下,克拉瑪依站(51243)10 m風(fēng)速預(yù)報(bào)和實(shí)況,預(yù)報(bào)時(shí)段為2021年2月24日15—25日12時(shí)。可以看出對(duì)于實(shí)況,25日04時(shí)之前風(fēng)速僅為5 m/s,25日06時(shí)整點(diǎn)風(fēng)速增加到21.3 m/s;對(duì)于GWD00方案,其起風(fēng)時(shí)間明顯過(guò)早,25日01時(shí)已達(dá)20 m/s左右,起風(fēng)時(shí)間相對(duì)于實(shí)況提前了5 h,此后數(shù)小時(shí)基本上處于20 m/s以上,個(gè)別時(shí)間可達(dá)24 m/s,在實(shí)況起風(fēng)時(shí)刻06時(shí)之后的風(fēng)速大小也明顯強(qiáng)于實(shí)況;對(duì)于GWD10,大風(fēng)起風(fēng)時(shí)間相對(duì)于實(shí)況依然偏早,在03時(shí)其風(fēng)速也達(dá)18 m/s,比實(shí)況提前3 h左右,但在06時(shí)之后與實(shí)況風(fēng)速十分接近;對(duì)于GWD11設(shè)置,由于D02加入了重力波拖曳,起風(fēng)時(shí)間相對(duì)于GWD10進(jìn)一步延后,風(fēng)速增加到20 m/s的時(shí)間為25日05時(shí),但06時(shí)之后風(fēng)速相對(duì)于實(shí)況也偏高。06時(shí)之后,不管是實(shí)況還是三種方案預(yù)報(bào),均呈現(xiàn)明顯減小的趨勢(shì)。綜合以上結(jié)果,D01增加GWD可以有效改善背風(fēng)坡山腳處克拉瑪依站的起風(fēng)時(shí)間和大風(fēng)持續(xù)時(shí)段的風(fēng)速預(yù)報(bào),而D02增加GWD可以進(jìn)一步改善大風(fēng)的起風(fēng)時(shí)間。

    圖4給出了三種方案不同預(yù)報(bào)時(shí)效(19~30 h)對(duì)克拉瑪依市區(qū)風(fēng)速預(yù)報(bào)的分布,可以看出模式預(yù)報(bào)的大風(fēng)不斷越過(guò)加依爾山背風(fēng)坡侵入克拉瑪依城區(qū),但是三種方案風(fēng)速大值區(qū)進(jìn)入的時(shí)間明顯不同,GWD00方案在21~24 h預(yù)報(bào)時(shí)效的風(fēng)速大值即位于加依爾山下坡處,侵入克拉瑪依西部的克拉瑪依主城區(qū),在27~30 h時(shí)段風(fēng)速也明顯較強(qiáng),在27~30 h時(shí)段可以發(fā)現(xiàn)在克拉瑪依西北部存在明顯的大值中心,量級(jí)可達(dá)24 m/s以上;GWD10方案在21~24 h風(fēng)速大值區(qū)也影響到了克拉瑪依城區(qū),在27~30 h時(shí)段風(fēng)速也明顯較強(qiáng),風(fēng)速大值中心可達(dá)24 m/s以上,只是大值范圍略小于GWD00方案;GWD11方案相對(duì)于前兩種方案,風(fēng)速明顯減弱,27 h預(yù)報(bào)時(shí)效之前,風(fēng)速大值區(qū)基本沒有影響到克拉瑪依城區(qū),且量級(jí)較小,直到28 h之后,克拉瑪依城區(qū)預(yù)報(bào)才可達(dá)到18 m/s以上,說(shuō)明兩重嵌套均開啟重力波拖曳,可以顯著降低風(fēng)速預(yù)報(bào),從而影響大風(fēng)進(jìn)入克拉瑪依城區(qū)的時(shí)間和量級(jí)。

    盧冰等(2014)對(duì)克拉瑪依強(qiáng)下坡風(fēng)形成機(jī)制的分析中指出,克拉瑪依上游地區(qū)中高層到低層存在強(qiáng)西風(fēng)氣流并在低層伴有冷平流的情況下,隨著風(fēng)速不斷增大,氣流能夠翻越加依爾山在背風(fēng)坡形成重力波,重力波能夠往上傳播影響到高空氣流并形成波破碎,破碎區(qū)湍流運(yùn)動(dòng)可導(dǎo)致上層能量向下傳播。為進(jìn)一步研究本次大風(fēng)個(gè)例中重力波拖曳方案對(duì)克拉瑪依翻山氣流大風(fēng)降溫的影響,給出了27 h預(yù)報(bào)時(shí)效的2021年2月25日03時(shí)GWD00和GWD10方案大風(fēng)進(jìn)入城區(qū)時(shí)刻沿加依爾山至克拉瑪依市區(qū)一線46°N的風(fēng)向和位溫預(yù)報(bào)垂直剖面,如圖5所示。從三個(gè)方案的剖面可以看出,首先對(duì)于高空風(fēng)場(chǎng)GWD00方案(圖5a)加依爾山上空6 km處基本處于較強(qiáng)的西風(fēng)氣流之中,垂直風(fēng)切變較小,說(shuō)明重力波傳播基本已經(jīng)結(jié)束;對(duì)于GWD10方案(圖5b),受外層重力波拖曳方案引起的風(fēng)場(chǎng)變化影響,地形背風(fēng)坡上空存在較明顯的垂直運(yùn)動(dòng),5 km以上存在上升氣流;GWD11方案(圖5c)在5 km以上垂直運(yùn)動(dòng)更為劇烈,具有明顯的垂直向上分量。5 km以下,所有方案中翻山氣流均向山麓推進(jìn),其中GWD00氣流翻山后在近地面以西風(fēng)分量為主,冷空氣快速進(jìn)入,如284 K等位溫線明顯向山腳延伸;GWD10受重力波破碎導(dǎo)致能量向下傳播影響,在5 km高度以下具有明顯的下沉氣流特征,山腳下克拉瑪依站風(fēng)速水平分量有所減弱,而垂直分量有所增加,冷空氣也明顯向山腳下延申;而GWD11方案翻山氣流水平風(fēng)速明顯小于GWD00和GWD10,風(fēng)的垂直分量更為顯著。以上分析說(shuō)明當(dāng)模式內(nèi)外層均未加入重力波拖曳作用,背風(fēng)坡上空的次網(wǎng)格湍流垂直擴(kuò)散未考慮,背風(fēng)坡處的水平風(fēng)速偏強(qiáng);外層加入重力波拖曳后,重力波垂直湍流擴(kuò)散顯著,削弱水平風(fēng)速項(xiàng),導(dǎo)致下坡風(fēng)的延遲和減弱;在內(nèi)層區(qū)域開啟重力波拖曳的方案,進(jìn)一步延緩了翻山氣流的水平風(fēng)速,并延遲了冷平流帶來(lái)的局地降溫。

    利用克拉瑪依附近14個(gè)區(qū)域自動(dòng)站實(shí)況資料對(duì)克拉瑪依市及其周邊的站點(diǎn)預(yù)報(bào)進(jìn)行檢驗(yàn),計(jì)算27 h預(yù)報(bào)時(shí)效(對(duì)應(yīng)時(shí)刻為25日03時(shí))的區(qū)域內(nèi)各站點(diǎn)上的預(yù)報(bào)偏差分布以及區(qū)域平均的平均誤差(Mean Error,ME,或稱之為偏差)以及平均絕對(duì)誤差 (Absolute Mean Error,AME),如圖6所示。對(duì)于GWD00(圖6a),加依爾山山腳下的克拉瑪依站具有明顯的風(fēng)速正偏差,可達(dá)10 m/s以上,其他站的預(yù)報(bào)偏差有正有負(fù),范圍內(nèi)總體的ME為2.19 m/s,AME為3.45 m/s,表現(xiàn)出明顯的正偏差;對(duì)于GWD10方案(圖6b),可以看出克拉瑪依站風(fēng)速正偏差依然為10 m/s以上范圍,而總體上的偏差有一定減少,ME達(dá)到1.57 m/s;GWD11方案,克拉瑪依站預(yù)報(bào)效果進(jìn)一步改善,正偏差降到2~5 m/s范圍內(nèi),總體偏差特征ME為-0.43 m/s,負(fù)偏差站點(diǎn)有所增多。

    總體上看GWD10整體上改善了GWD00的正偏差,但對(duì)克拉瑪依站的起風(fēng)時(shí)間預(yù)報(bào)依然過(guò)早;而GWD11方案進(jìn)一步改善了克拉瑪依站及其附近站點(diǎn)的正偏差,但帶來(lái)了一定的負(fù)偏差。

    2.2 檢驗(yàn)結(jié)果

    對(duì)批量試驗(yàn)時(shí)段內(nèi)的預(yù)報(bào)結(jié)果進(jìn)行評(píng)分統(tǒng)計(jì),進(jìn)一步定性和定量分析不同重力波拖曳設(shè)置下的模式預(yù)報(bào)效果。

    首先,定性分析重力波拖曳方案對(duì)預(yù)報(bào)效果的影響,本節(jié)的實(shí)況分析場(chǎng)采用檢驗(yàn)時(shí)刻對(duì)應(yīng)的WRF初始分析場(chǎng)(由GFS初始場(chǎng)驅(qū)動(dòng)并同化了相應(yīng)的探空和地面站觀測(cè)資料)。圖7給出了批量試驗(yàn)時(shí)段內(nèi)統(tǒng)計(jì)平均的24 h預(yù)報(bào)時(shí)效的D02區(qū)域內(nèi)緯向U風(fēng)場(chǎng)預(yù)報(bào)與預(yù)報(bào)時(shí)刻對(duì)應(yīng)的分析場(chǎng)之差沿85°E剖面,該剖面位置為加依爾山背風(fēng)坡。

    從圖7a可知GWD00方案的U風(fēng)預(yù)報(bào)在加依爾山背風(fēng)坡相對(duì)于分析場(chǎng)具有明顯的正偏差,這種正偏差從高空200 hPa一直延續(xù)到地面背風(fēng)坡上,說(shuō)明無(wú)GWD作用下,山區(qū)地形上空緯向風(fēng)速具有顯著偏大特征;在D01方案中開啟GWD,即GWD10方案,其對(duì)D02緯向風(fēng)速也有顯著影響,如200 hPa左右正偏差得到了顯著削弱,甚至出現(xiàn)了一定的負(fù)偏差,對(duì)于350 hPa以下,46°~47.4°N范圍內(nèi),部分區(qū)域偏差已經(jīng)降低到-2~2 m/s;當(dāng)D02開啟GWD方案后,即GWD11,可以看出,對(duì)于高空風(fēng)正偏差區(qū)域范圍進(jìn)一步縮小,而負(fù)偏差特征進(jìn)一步增大,需要注意的是,GWD11設(shè)置下對(duì)于背風(fēng)坡近地面處的風(fēng)速偏差確實(shí)有明顯改善,如750 hPa高度的山坡區(qū)域風(fēng)速正偏差僅為2 m/s左右,顯著小于GWD00及GWD10方案,而對(duì)于無(wú)地形區(qū)域的近地面及高空風(fēng)場(chǎng),GWD11相對(duì)于前兩種方案顯著增加了風(fēng)速負(fù)偏差,這也說(shuō)明開啟D02的重力波拖曳作用對(duì)全場(chǎng)預(yù)報(bào)效果具有較大的不一致性。

    為了檢驗(yàn)三種方案統(tǒng)計(jì)預(yù)報(bào)效果,計(jì)算了三種方案試驗(yàn)時(shí)段內(nèi)D02區(qū)域不同要素統(tǒng)計(jì)平均的均方根誤差(RMSE)以及偏差(BIAS),包括2 m高度氣溫(T2m)、10 m高度緯向風(fēng)(U10m)以及10 m高度徑向風(fēng)(V10m),實(shí)況采用D02區(qū)域內(nèi)所有地面自動(dòng)站,間隔6 h。圖8給出了三種方案36 h預(yù)報(bào)時(shí)效內(nèi)D02區(qū)域地面要素預(yù)報(bào)的RMSE和BIAS的演變特征,可以看出,對(duì)于RMSE,不管是溫度預(yù)報(bào)還是風(fēng)場(chǎng)預(yù)報(bào),GWD10在各個(gè)預(yù)報(bào)時(shí)效均能夠獲得最低的RMSE,而GWD11方案對(duì)不同要素預(yù)報(bào)的RMSE均較大,GWD00方案居中,如對(duì)于U風(fēng)場(chǎng)的24 h預(yù)報(bào),GWD00誤差為3.25 m/s,GWD10為2.86 m/s,GWD11為3.35 m/s。從BIAS可以看出,對(duì)于風(fēng)場(chǎng)預(yù)報(bào)(圖8b和圖8c),GWD00表現(xiàn)出一定的正偏差,尤其是12~30 h預(yù)報(bào)時(shí)效較為明顯,如18 h預(yù)報(bào)時(shí)效的BIAS為0.23 m/s;GWD10相對(duì)于GWD00風(fēng)場(chǎng)預(yù)報(bào)的正偏差有一定減小,18 h預(yù)報(bào)時(shí)效的BIAS為-0.16 m/s;而GWD11方案相對(duì)于GWD10方案的風(fēng)場(chǎng)有了明顯的負(fù)偏差,18 h預(yù)報(bào)時(shí)效的BIAS為-0.47 m/s,這也與上文中的定性分析結(jié)論相一致。以上分析說(shuō)明整體上GWD10設(shè)置下的預(yù)報(bào)效果較好,而GWD11設(shè)置下會(huì)對(duì)風(fēng)場(chǎng)存在“過(guò)度矯正”問(wèn)題,引起偏差增大。

    圖9給出了試驗(yàn)時(shí)段內(nèi)探空觀測(cè)檢驗(yàn)的12、24、36 h三個(gè)預(yù)報(bào)時(shí)效平均的RMSE垂直分布,其他非整點(diǎn)時(shí)次由于探空資料較少,因此不計(jì)入統(tǒng)計(jì)。首先可以看出對(duì)于U風(fēng)預(yù)報(bào)(圖9a),三種方案各個(gè)層次誤差量級(jí)不同,其中在400 hPa誤差最大,GWD00方案誤差達(dá)到6.6 m/s,而GWD10和GWD11方案誤差為5.5 m/s左右,顯著小于GWD00方案,對(duì)于模式低層,GWD10的誤差相對(duì)也較小,如850 hPa誤差為1.63 m/s,小于GWD00的2.17 m/s和GWD11的1.87 m/s,對(duì)于另外兩個(gè)要素,兩種重力波方案相對(duì)于不加入重力波方案也顯示出一定的優(yōu)勢(shì),綜合而言GWD10在不同等壓面均可獲得更好的要素預(yù)報(bào),尤其是對(duì)于邊界層而言,GWD10的預(yù)報(bào)可以獲得更好的熱動(dòng)力狀況,從而有效改善克拉瑪依城區(qū)翻山下坡風(fēng)的預(yù)報(bào)。

    綜上所述,在外層區(qū)域D01開啟重力波拖曳,內(nèi)層區(qū)域D02關(guān)閉重力波拖曳,可以獲得較好的綜合效果。雖然對(duì)于克拉瑪依單站預(yù)報(bào),D02也開啟重力波拖曳有可能會(huì)獲得較好的效果,但是會(huì)影響其他非地形陡峭區(qū)域的風(fēng)場(chǎng)預(yù)報(bào)。而D01和D02均不開啟重力波拖曳,會(huì)使風(fēng)場(chǎng)存在固有的系統(tǒng)性正偏差,因此從日常需求來(lái)看,在D01開啟重力波拖曳是比較合理的設(shè)置。

    3 結(jié)論和討論

    基于克拉瑪依數(shù)值系統(tǒng)開展了冬季大風(fēng)預(yù)報(bào)試驗(yàn),對(duì)比了不同重力波拖曳方案設(shè)置即兩重嵌套均關(guān)閉重力波拖曳(GWD00)、僅開啟D01重力波拖曳(GWD10)和兩重嵌套均關(guān)閉重力波拖曳(GWD11)三組方案對(duì)克拉瑪依大風(fēng)預(yù)報(bào)的效果,重點(diǎn)分析了2月25日大風(fēng)預(yù)報(bào)個(gè)例,并通過(guò)地面要素檢驗(yàn)指標(biāo)對(duì)三種方案的預(yù)報(bào)效果進(jìn)行綜合評(píng)估,得出以下結(jié)論:

    1)對(duì)于克拉瑪依冬季穩(wěn)定層結(jié)下的一次大風(fēng)個(gè)例,三種方案的起風(fēng)時(shí)間各不相同,其中GWD00方案起風(fēng)過(guò)早,大風(fēng)延續(xù)時(shí)間最長(zhǎng),大風(fēng)量級(jí)相比實(shí)況有正偏差;GWD10方案起風(fēng)時(shí)間有所延遲,且大風(fēng)量級(jí)與實(shí)況最接近;GWD11方案起風(fēng)時(shí)間最接近實(shí)況。說(shuō)明本次個(gè)例,重力波拖曳可以有效減弱地形陡峭區(qū)域大風(fēng)預(yù)報(bào)正偏差。

    2)對(duì)克拉瑪依周邊單站預(yù)報(bào)進(jìn)行分析,表明雖然山區(qū)背風(fēng)坡的克拉瑪依站在兩重嵌套中引入重力波拖曳可以有效減小近地面風(fēng)場(chǎng)誤差,但是對(duì)其他區(qū)域站點(diǎn)預(yù)報(bào)則帶來(lái)一定的負(fù)效果,說(shuō)明兩重嵌套均開啟重力波拖曳并不能全面改善預(yù)報(bào)區(qū)域內(nèi)的效果。

    3)進(jìn)行了批量試驗(yàn)檢驗(yàn),定性分析表明:未開啟重力波拖曳的情況下會(huì)在山區(qū)地形處存在明顯正偏差,在外層區(qū)域開啟重力波拖曳則會(huì)在平流層至地面均有效減少風(fēng)場(chǎng)預(yù)報(bào)正偏差,而內(nèi)外層均開啟重力波拖曳,雖然會(huì)有效減小背風(fēng)坡近地面風(fēng)速正偏差,但會(huì)在非山地區(qū)域增加負(fù)偏差;通過(guò)統(tǒng)計(jì)檢驗(yàn)發(fā)現(xiàn)相對(duì)于未開啟重力波拖曳的設(shè)置,在外層開啟重力波拖曳可以有效降低不同預(yù)報(bào)時(shí)效地面要素均方根誤差,但是進(jìn)一步開啟內(nèi)層區(qū)域的重力波拖曳,會(huì)導(dǎo)致均方根誤差有所增加,主要是由于開啟外層重力波拖曳增加了一定的負(fù)偏差。

    以上結(jié)論表明WRF中重力波拖曳方案,由于包含重力波破碎和阻塞拖曳的模型均針對(duì)大尺度山脈對(duì)對(duì)流層動(dòng)量收支的影響,因此在外層粗分辨率區(qū)域開啟重力波拖曳,內(nèi)層關(guān)閉重力波拖曳時(shí),區(qū)域整體預(yù)報(bào)具有較好的效果,但對(duì)于克拉瑪依這種處于山區(qū)背風(fēng)坡位置的站點(diǎn),開啟內(nèi)層重力波拖曳可獲得更好的預(yù)報(bào)效果(而其他站點(diǎn)預(yù)報(bào)效果會(huì)變差)。如果既提高背風(fēng)坡站點(diǎn)預(yù)報(bào),又保持其他非地形區(qū)域預(yù)報(bào)技巧,則需要引入針對(duì)高分辨率地形網(wǎng)格的重力波拖曳效應(yīng),充分描述穩(wěn)定邊界層內(nèi)的缺失拖曳作用。Steeneveld et al.(2008)提出相對(duì)較淺的穩(wěn)定邊界層內(nèi)的拖曳作用可以通過(guò)小尺度地形來(lái)獲取,并形成了穩(wěn)定邊界層重力波拖曳方案GWDSBL,該方案已初步耦合到WRF V4.3中(Tsiringakis et al.,2017)。此種針對(duì)穩(wěn)定邊界層的重力波拖曳方案在克拉瑪依高分辨率數(shù)值預(yù)報(bào)系統(tǒng)中的表現(xiàn)如何,未來(lái)需要進(jìn)一步開展試驗(yàn)。

    參考文獻(xiàn)(References)

    Bae S Y,Hong S Y,Lim K S S,2016.Coupling WRF double-moment 6-class microphysics schemes to RRTMG radiation scheme in weather research forecasting model[J].Adv Meteorol,2016:1-11.doi:10.1155/2016/5070154.

    Baines P G,Manins P C,1989.The principles of laboratory modeling of stratified atmospheric flows over complex terrain[J].J Appl Meteor,28(11):1213-1225.doi:10.1175/1520-0450(1989)028<1213:tpolmo>2.0.co;2.

    Baines P G,Palmer,T N,1990.Rationale for a new physically based parametrization of sub-grid scale orographic effects[J].Ecmwf Tech.doi:10.21957/h4h36b3u.

    Boer G J,McFarlane N A,Laprise R,et al.,1984.The Canadian Climate Centre spectral atmospheric general circulation model[J].Atmosphere-Ocean,22(4):397-429.doi:10.1080/07055900.1984.9649208.

    程胡華,2017.不同擾動(dòng)場(chǎng)對(duì)大氣重力波參數(shù)結(jié)果影響的初步探討[J].大氣科學(xué)學(xué)報(bào),40(3):401-411. Cheng H H,2017.A preliminary investigation on the atmospheric gravity wave parameters obtained from different disturbance fields[J].Trans Atmos Sci,40(3):401-411.doi:10.13878/j.cnki.dqkxxb.20160504002.(in Chinese).

    Clark T L,Peltier W R,1984.Critical level reflection and the resonant growth of nonlinear mountain waves[J].J Atmos Sci,41(21):3122-3134.doi:10.1175/1520-0469(1984)041<3122:clratr>2.0.co;2.

    Choi H J,Hong S Y,2015.An updated subgrid orographic parameterization for global atmospheric forecast models[J].J Geophys Res:Atmos,120(24):12445-12457.doi:10.1002/2015JD024230.

    Hong S Y,Noh Y,Dudhia J,2006.A new vertical diffusion package with an explicit treatment of entrainment processes[J].Mon Wea Rev,134(9):2318-2341.doi:10.1175/mwr3199.1.

    Hong S Y,Choi J,Chang E C,et al.,2008.Lower-tropospheric enhancement of gravity wave drag in a global spectral atmospheric forecast model[J].Wea Forecasting,23(3):523-531.doi:10.1175/2007waf2007030.1.

    Kim Y J,Arakawa A,1995.Improvement of orographic gravity wave parameterization using a mesoscale gravity wave model[J].J Atmos Sci,52(11):1875-1902.doi:10.1175/1520-0469(1995)052<1875:ioogwp>2.0.co;2.

    Kim Y J,Doyle J D,2005.Extension of an orographic-drag parametrization scheme to incorporate orographic anisotropy and flow blocking[J].Quart J Roy Meteor Soc,131(609):1893-1921.doi:10.1256/qj.04.160.

    Klemp J B,Lilly D R,1975.The dynamics of wave-induced downslope winds[J].J Atmos Sci,32(2):320-339.doi:10.1175/1520-0469(1975)032<0320:tdowid>2.0.co;2.

    Lott F,Miller M J,1997.A new subgrid-scale orographic drag parametrization:Its formulation and testing[J].Quart J Roy Meteor Soc,123(537):101-127.doi:10.1002/qj.49712353704.

    Lilly D K,1978.A severe downslope windstorm and aircraft turbulence event induced by a mountain wave[J].J Atmos Sci,35(1):59-77.doi:10.1175/1520-0469(1978)035<0059:asdwaa>2.0.co;2.

    Lindzen R S,1981.Turbulence and stress owing to gravity wave and tidal breakdown[J].J Geophys Res:Oceans,86(C10):9707-9714.doi:10.1029/JC086iC10p09707.

    劉華強(qiáng),錢永甫,2001.包絡(luò)地形和重力波拖曳對(duì)區(qū)域氣候模擬效果的影響[J].大氣科學(xué),25(2):209-220. Liu H Q,Qian Y F,2001.Effects of envelope orography and gravity wave drag parameterization on regional climate simulations[J].Chin J Atmos Sci,25(2):209-220.(in Chinese).

    盧冰,史永強(qiáng),王光輝,等,2014.新疆克拉瑪依強(qiáng)下坡風(fēng)暴的機(jī)理研究[J].氣象學(xué)報(bào),72(6):1218-1230. Lu B,Shi Y Q,Wang G H,et al.,2014.Numerical study of severe downslope winds at Kelamayi,Xinjiang[J].Acta Meteorol Sin,72(6):1218-1230.(in Chinese).

    Miller M J,Palmer T N,Swinbank R,1989.Parametrization and influence of subgridscale orography in general circulation and numerical weather prediction models[J].Meteor Atmos Phys,40(1/2/3):84-109.doi:10.1007/BF01027469.

    McFarlane N A,1987.The effect of orographically excited gravity wave drag on the general circulation of the lower stratosphere and troposphere[J].J Atmos Sci,44(14):1775-1800.doi:10.1175/1520-0469(1987)044<1775:teooeg>2.0.co;2.

    Palmer T N,Shutts G J,Swinbank R,1986.Alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization[J].Quart J Roy Meteor Soc,112(474):1001-1039.doi:10.1002/qj.49711247406.

    Phillips D S,1984.Analytical surface pressure and drag for linear hydrostatic flow over three-dimensional elliptical mountains[J].J Atmos Sci,41(6):1073-1084.doi:10.1175/1520-0469(1984)041<1073:aspadf>2.0.co;2.

    錢永甫,2000.包絡(luò)地形和重力波拖曳對(duì)氣候模擬效果的影響[J].應(yīng)用氣象學(xué)報(bào),11(1):13-20. Qian Y F,2000.Effects of envelope orography and gravity wave drag on performance of climate modeling[J].J Appl Meteor Sci,11(1):13-20.doi:10.3969/j.issn.1001-7313.2000.01.002.(in Chinese).

    Saito K,1993.A numerical study of the local downslope wind “yamaji-kaze” in Japan[J].Journal Meteorological Society Japan,71(2):247-272.doi:10.2151/jmsj1965.71.2_247.

    Sandu I,Beljaars A,Bechtold P,et al.,2013.Why is it so difficult to represent stably stratified conditions in numerical weather prediction (NWP) models?[J].J Adv Modeling Earth Syst,5(2):117-133.doi:10.1002/jame.20013.

    Scinocca J F,McFarlane N A,2000.The parametrization of drag induced by stratified flow over anisotropic orography[J].Quart J Roy Meteor Soc,126(568):2353-2393.doi:10.1002/qj.49712656802.

    Steeneveld G J,Holtslag A A M,Nappo C J,et al.,2008.Exploring the possible role of small-scale terrain drag on stable boundary layers over land[J].J Appl Meteorol Climatol,47(10):2518-2530.doi:10.1175/2008jamc1816.1.

    Tsiringakis A,Steeneveld G J,Holtslag A A M,2017.Small-scale orographic gravity wave drag in stable boundary layers and its impact on synoptic systems and near-surface meteorology[J].Quart J Roy Meteor Soc,143(704):1504-1516.doi:10.1002/qj.3021.

    孫東霞,謝小紅,郭曉靜,2008.克拉瑪依特強(qiáng)大風(fēng)的氣候特征及天氣分析與預(yù)報(bào)[J].沙漠與綠洲氣象,2(4):18-21. Sun D X,Xie X H,Guo X J,2008.The climatic characteristics of the extreme strong wind in Kelamayi and its weather analysis and forecast[J].Desert Oasis Meteorol,2(4):18-21.doi:10.3969/j.issn.1002-0799.2008.04.003.(in Chinese).

    Tibaldi S,1986.Envelope orography and maintenance of the quasi-stationary circulation in the ECMWF global models[J].Adv Geophys,29:339-374.doi:10.1016/S0065-2687(08)60045-X.

    Thompson G,Rasmussen R M,Manning K,2004.Explicit forecasts of winter precipitation using an improved bulk microphysics scheme.part Ⅰ:description and sensitivity analysis[J].Mon Wea Rev,132(2):519-542.doi:10.1175/1520-0493(2004)132<0519:efowpu>2.0.co;2.

    Wallace J M,Tibaldi S,Simmons A J,1983.Reduction of systematic forecast errors in the ECMWF model through the introduction of an envelope orography[J].Quart J Roy Meteor Soc,109(462):683-717.doi:10.1002/qj.49710946202.

    王宇,黃思訓(xùn),2018.幾種特殊地形對(duì)氣流運(yùn)動(dòng)影響的解析研究[J].大氣科學(xué)學(xué)報(bào),41(6):743-752. Wang Y,Huang S X,2018.Analytical study of the influence of several special topographies on airflow movement[J].Trans Atmos Sci,41(6):743-752.doi:10.13878/j.cnki.dqkxxb.20180114001.(in Chinese).

    辛渝,周穎真,宋麗莉,等,2015.新疆風(fēng)區(qū)下墊面動(dòng)力學(xué)粗糙度估測(cè)值與模式定義值的對(duì)比[J].大氣科學(xué)學(xué)報(bào),38(6):811-818. Xin Y,Zhou Y Z,Song L L,et al.,2015.Aerodynamic roughness comparison between the estimated based on masts and the defined in model over Xinjiang wind areas[J].Trans Atmos Sci,38(6):811-818.doi:10.13878/j.cnki.dqkxxb.20141125002.(in Chinese).

    徐國(guó)強(qiáng),楊學(xué)勝,黃麗萍,等,2010.GRAPES中地形重力波拖曳物理過(guò)程的引進(jìn)和應(yīng)用試驗(yàn)[J].氣象學(xué)報(bào),68(5):631-639. Xu G Q,Yang X S,Huang L P,et al.,2010.Introducing and application testing of the orographic gravity wave drag parameterization physics in the GRAPES[J].Acta Meteorol Sin,68(5):631-639.(in Chinese).

    Zhang C X,Wang Y Q,Hamilton K,2011.Improved representation of boundary layer clouds over the southeast Pacific in ARW-WRF using a modified tiedtke cumulus parameterization scheme[J].Mon Wea Rev,139(11):3489-3513.doi:10.1175/mwr-d-10-05091.1.

    鐘水新,陳子通,戴光豐,等,2014.地形重力波拖曳參數(shù)化對(duì)熱帶氣旋強(qiáng)度和路徑預(yù)報(bào)影響的研究[J].大氣科學(xué),38(2):273-284. Zhong S X,Chen Z T,Dai G F,et al.,2014.Impacts of orographic gravity wave drag parameterization on typhoon intensity and path forecasting[J].Chin J Atmos Sci,38(2):273-284.doi:10.3878/j.issn.1006-9895.2013.13131.(in Chinese).

    周天軍,錢永甫,1996.地形效應(yīng)影響數(shù)值預(yù)報(bào)結(jié)果的試驗(yàn)研究[J].大氣科學(xué),20(4):452-462. Zhou T J,Qian Y F,1996.An experimental study on the effects of topography on numerical prediction[J].Sci Atmos Sin,20(4):452-462.(in Chinese).

    (責(zé)任編輯:袁東敏)

    波多野结衣巨乳人妻| 搡老岳熟女国产| 99久久精品热视频| 一边摸一边抽搐一进一小说| 性欧美人与动物交配| 成人欧美大片| 色哟哟哟哟哟哟| 美女午夜性视频免费| 欧美性猛交黑人性爽| 麻豆国产av国片精品| www.精华液| 亚洲av片天天在线观看| 国产精品综合久久久久久久免费| 亚洲熟妇中文字幕五十中出| a级毛片在线看网站| 特大巨黑吊av在线直播| 国产精品99久久久久久久久| 国产精品久久久久久亚洲av鲁大| 亚洲av免费在线观看| 一边摸一边抽搐一进一小说| 久久午夜综合久久蜜桃| 精品国产超薄肉色丝袜足j| av国产免费在线观看| 免费观看人在逋| 可以在线观看的亚洲视频| 校园春色视频在线观看| 日本一本二区三区精品| 亚洲午夜精品一区,二区,三区| 亚洲欧美日韩东京热| 国产三级在线视频| 天天添夜夜摸| 99久久精品国产亚洲精品| 脱女人内裤的视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久人人做人人爽| 免费看美女性在线毛片视频| 99久国产av精品| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| 久久久久久久午夜电影| 操出白浆在线播放| 国产1区2区3区精品| 久久精品91蜜桃| 久久中文看片网| 欧美成人一区二区免费高清观看 | 亚洲激情在线av| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 亚洲无线观看免费| 国产野战对白在线观看| 噜噜噜噜噜久久久久久91| 亚洲人与动物交配视频| 国产伦一二天堂av在线观看| 男女下面进入的视频免费午夜| 国产精品av久久久久免费| 操出白浆在线播放| 午夜日韩欧美国产| 久久精品人妻少妇| 俄罗斯特黄特色一大片| www.精华液| www日本在线高清视频| 日本五十路高清| 久久精品综合一区二区三区| 最新在线观看一区二区三区| 欧美乱码精品一区二区三区| 午夜福利18| 国内毛片毛片毛片毛片毛片| 亚洲精品乱码久久久v下载方式 | 18禁国产床啪视频网站| 99riav亚洲国产免费| 国产人伦9x9x在线观看| 国产99白浆流出| 国产成年人精品一区二区| 欧美成狂野欧美在线观看| 一本精品99久久精品77| 国产高潮美女av| 床上黄色一级片| 狂野欧美白嫩少妇大欣赏| 小蜜桃在线观看免费完整版高清| 两个人的视频大全免费| 久久久久久国产a免费观看| 亚洲国产色片| 国产激情偷乱视频一区二区| 午夜福利在线在线| 亚洲七黄色美女视频| 国产单亲对白刺激| 欧美精品啪啪一区二区三区| avwww免费| av视频在线观看入口| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产看品久久| www.www免费av| 最近最新中文字幕大全电影3| 一二三四在线观看免费中文在| 国内揄拍国产精品人妻在线| 好男人在线观看高清免费视频| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| 久久久久久大精品| www.精华液| 一个人看的www免费观看视频| 国产精品香港三级国产av潘金莲| www.www免费av| 淫妇啪啪啪对白视频| 12—13女人毛片做爰片一| 久久香蕉国产精品| 精品无人区乱码1区二区| 国产伦一二天堂av在线观看| 亚洲激情在线av| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| 成人国产一区最新在线观看| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩高清在线视频| 精品久久久久久,| 久久久久久久精品吃奶| 欧美最黄视频在线播放免费| 丁香欧美五月| 久久午夜亚洲精品久久| 少妇的逼水好多| 欧美日本亚洲视频在线播放| 欧美在线一区亚洲| 日本一本二区三区精品| 亚洲成人精品中文字幕电影| 欧美在线黄色| 久久亚洲真实| 国产探花在线观看一区二区| av国产免费在线观看| 精品久久久久久久末码| 啦啦啦观看免费观看视频高清| 999精品在线视频| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 无限看片的www在线观看| 成人三级黄色视频| 真人一进一出gif抽搐免费| 中文字幕高清在线视频| 美女免费视频网站| 香蕉国产在线看| 国产又黄又爽又无遮挡在线| 日本五十路高清| 欧美日韩乱码在线| 男女那种视频在线观看| 波多野结衣高清无吗| 日本免费a在线| www.精华液| 国产精品亚洲一级av第二区| 叶爱在线成人免费视频播放| 日本三级黄在线观看| 欧美日韩黄片免| 久久香蕉国产精品| 国产久久久一区二区三区| 夜夜看夜夜爽夜夜摸| 国产一区二区在线av高清观看| 熟女人妻精品中文字幕| 99国产精品一区二区三区| 国产高清激情床上av| 午夜福利视频1000在线观看| 国产精品99久久久久久久久| 国产av不卡久久| 夜夜夜夜夜久久久久| 欧美不卡视频在线免费观看| 欧美性猛交黑人性爽| 色综合欧美亚洲国产小说| 白带黄色成豆腐渣| 欧美性猛交黑人性爽| 成人亚洲精品av一区二区| 动漫黄色视频在线观看| 免费观看人在逋| 亚洲精品乱码久久久v下载方式 | 精品久久久久久久毛片微露脸| 欧美丝袜亚洲另类 | 给我免费播放毛片高清在线观看| 亚洲午夜理论影院| 久久亚洲精品不卡| 欧美日韩一级在线毛片| www日本在线高清视频| 午夜免费观看网址| 欧美3d第一页| 黑人巨大精品欧美一区二区mp4| 五月玫瑰六月丁香| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 我要搜黄色片| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 夜夜爽天天搞| 国产欧美日韩精品一区二区| 99久久精品热视频| 亚洲欧美激情综合另类| 麻豆成人av在线观看| 亚洲中文av在线| 成人亚洲精品av一区二区| 欧美日韩黄片免| 精品久久久久久久毛片微露脸| 成人三级做爰电影| 99视频精品全部免费 在线 | 国产视频一区二区在线看| 99在线人妻在线中文字幕| 精品久久蜜臀av无| 999精品在线视频| 51午夜福利影视在线观看| 日韩欧美国产在线观看| 啦啦啦免费观看视频1| 久久久色成人| 欧美不卡视频在线免费观看| 亚洲av日韩精品久久久久久密| 久久久久国产一级毛片高清牌| 麻豆成人av在线观看| 国产精品爽爽va在线观看网站| 91麻豆av在线| 亚洲专区中文字幕在线| 最近最新免费中文字幕在线| 国产伦在线观看视频一区| 黄片小视频在线播放| 美女 人体艺术 gogo| 国产乱人伦免费视频| 女生性感内裤真人,穿戴方法视频| 国产单亲对白刺激| 国产高清视频在线播放一区| 久久久久国内视频| 国内精品久久久久精免费| 日本黄大片高清| 91av网站免费观看| 欧美极品一区二区三区四区| 久久香蕉精品热| 国产精品精品国产色婷婷| 国产综合懂色| 男人舔女人下体高潮全视频| 欧美黄色淫秽网站| 午夜视频精品福利| 免费人成视频x8x8入口观看| 在线观看舔阴道视频| 国产91精品成人一区二区三区| 国产精品美女特级片免费视频播放器 | 美女大奶头视频| 99在线视频只有这里精品首页| 香蕉久久夜色| 亚洲精品国产精品久久久不卡| 成年女人毛片免费观看观看9| 全区人妻精品视频| 两个人视频免费观看高清| 国产伦精品一区二区三区四那| 91av网一区二区| 久9热在线精品视频| 在线十欧美十亚洲十日本专区| 亚洲精品色激情综合| 99久国产av精品| 日韩人妻高清精品专区| 国产亚洲av嫩草精品影院| 欧美绝顶高潮抽搐喷水| 日韩三级视频一区二区三区| 国产 一区 欧美 日韩| 久久久国产欧美日韩av| 亚洲五月婷婷丁香| 欧美激情久久久久久爽电影| 国内毛片毛片毛片毛片毛片| 亚洲中文日韩欧美视频| 国产精品免费一区二区三区在线| 久久久国产成人精品二区| 亚洲成人久久性| 一个人观看的视频www高清免费观看 | 又粗又爽又猛毛片免费看| 亚洲精华国产精华精| 国产亚洲精品一区二区www| АⅤ资源中文在线天堂| 黑人欧美特级aaaaaa片| 久久99热这里只有精品18| 美女高潮的动态| 悠悠久久av| 啪啪无遮挡十八禁网站| 亚洲国产精品合色在线| 国产精品女同一区二区软件 | 国产一区二区激情短视频| 久久久久精品国产欧美久久久| 91老司机精品| 午夜免费成人在线视频| 午夜免费激情av| www日本在线高清视频| 国产一区二区在线av高清观看| 99热这里只有精品一区 | 18美女黄网站色大片免费观看| 亚洲国产欧美网| 免费av毛片视频| 中出人妻视频一区二区| 亚洲自拍偷在线| bbb黄色大片| 三级男女做爰猛烈吃奶摸视频| 淫秽高清视频在线观看| 天堂网av新在线| 一区福利在线观看| 久久久久性生活片| 神马国产精品三级电影在线观看| 久久午夜亚洲精品久久| 动漫黄色视频在线观看| 午夜福利视频1000在线观看| av福利片在线观看| 狠狠狠狠99中文字幕| 韩国av一区二区三区四区| 国产一区二区在线av高清观看| 久久伊人香网站| 少妇人妻一区二区三区视频| 成人鲁丝片一二三区免费| ponron亚洲| 香蕉国产在线看| 免费无遮挡裸体视频| 亚洲九九香蕉| 99久久99久久久精品蜜桃| 国产私拍福利视频在线观看| 精品久久久久久久久久久久久| 久久人人精品亚洲av| 一级作爱视频免费观看| 午夜视频精品福利| 一本久久中文字幕| 国产精品久久久久久久电影 | 婷婷丁香在线五月| 村上凉子中文字幕在线| 国产精品亚洲av一区麻豆| 国产高潮美女av| 国产亚洲av嫩草精品影院| 精品不卡国产一区二区三区| 一边摸一边抽搐一进一小说| 国产亚洲av高清不卡| 一夜夜www| 9191精品国产免费久久| АⅤ资源中文在线天堂| 国产激情欧美一区二区| 国产三级黄色录像| 天堂网av新在线| 精品乱码久久久久久99久播| 真实男女啪啪啪动态图| 亚洲av免费在线观看| 中文字幕最新亚洲高清| 香蕉av资源在线| svipshipincom国产片| 午夜影院日韩av| 在线看三级毛片| 97碰自拍视频| e午夜精品久久久久久久| 免费在线观看影片大全网站| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 国产精品久久视频播放| 国产亚洲av高清不卡| 91久久精品国产一区二区成人 | 黄色日韩在线| 一边摸一边抽搐一进一小说| 国产精品久久久久久久电影 | 熟女少妇亚洲综合色aaa.| 人人妻,人人澡人人爽秒播| bbb黄色大片| 日本黄大片高清| 禁无遮挡网站| 两个人的视频大全免费| 女警被强在线播放| 999久久久国产精品视频| 欧美成狂野欧美在线观看| 变态另类成人亚洲欧美熟女| 中国美女看黄片| 国产毛片a区久久久久| 国产黄色小视频在线观看| 熟妇人妻久久中文字幕3abv| 一夜夜www| 亚洲欧美日韩东京热| 一夜夜www| 在线十欧美十亚洲十日本专区| av在线天堂中文字幕| 91九色精品人成在线观看| 欧美性猛交黑人性爽| 一区二区三区国产精品乱码| av在线天堂中文字幕| 国产亚洲精品一区二区www| 亚洲九九香蕉| 欧美乱色亚洲激情| 最近最新中文字幕大全免费视频| 亚洲国产精品合色在线| a级毛片在线看网站| 三级男女做爰猛烈吃奶摸视频| 免费人成视频x8x8入口观看| 欧美日韩综合久久久久久 | 中文字幕av在线有码专区| 国产一区二区激情短视频| 一区二区三区国产精品乱码| 国产精品女同一区二区软件 | 欧美在线一区亚洲| 国语自产精品视频在线第100页| 18禁美女被吸乳视频| 一个人看视频在线观看www免费 | 国产v大片淫在线免费观看| 美女免费视频网站| 欧美国产日韩亚洲一区| 日本精品一区二区三区蜜桃| 亚洲av中文字字幕乱码综合| 成在线人永久免费视频| 一本综合久久免费| 午夜亚洲福利在线播放| h日本视频在线播放| 又大又爽又粗| 久久精品91蜜桃| 麻豆国产av国片精品| 国内精品美女久久久久久| 麻豆国产av国片精品| 色哟哟哟哟哟哟| 亚洲成人久久爱视频| 日韩欧美国产一区二区入口| 老司机午夜福利在线观看视频| 村上凉子中文字幕在线| 精品久久久久久成人av| 丰满人妻一区二区三区视频av | 国产亚洲精品av在线| 久久中文字幕一级| 天天一区二区日本电影三级| 免费人成视频x8x8入口观看| 日韩中文字幕欧美一区二区| 日韩有码中文字幕| 久久人人精品亚洲av| 特大巨黑吊av在线直播| 国产精品一及| 国产欧美日韩精品一区二区| 久久中文字幕一级| www日本在线高清视频| av视频在线观看入口| av福利片在线观看| 久9热在线精品视频| 午夜视频精品福利| 久久这里只有精品19| www.自偷自拍.com| 欧美黄色淫秽网站| 国产不卡一卡二| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久久久精品电影| 黄色丝袜av网址大全| 深夜精品福利| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 狠狠狠狠99中文字幕| 变态另类成人亚洲欧美熟女| 麻豆成人午夜福利视频| 人人妻人人澡欧美一区二区| 欧美中文综合在线视频| 丰满人妻一区二区三区视频av | 国产精品亚洲美女久久久| 国产成人精品久久二区二区91| 亚洲色图av天堂| 日本一二三区视频观看| 此物有八面人人有两片| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久人人做人人爽| 午夜久久久久精精品| 中文字幕人妻丝袜一区二区| 波多野结衣巨乳人妻| 国产黄色小视频在线观看| 欧美不卡视频在线免费观看| 一边摸一边抽搐一进一小说| 午夜福利高清视频| 久久久久国产一级毛片高清牌| 国产不卡一卡二| 国产又黄又爽又无遮挡在线| 午夜福利免费观看在线| 国产成+人综合+亚洲专区| 亚洲黑人精品在线| 国产视频一区二区在线看| 欧美乱妇无乱码| 欧美又色又爽又黄视频| 亚洲男人的天堂狠狠| 在线免费观看的www视频| 精品久久久久久久人妻蜜臀av| 噜噜噜噜噜久久久久久91| 又黄又爽又免费观看的视频| 男插女下体视频免费在线播放| 1000部很黄的大片| www国产在线视频色| 免费av毛片视频| 窝窝影院91人妻| 日本免费a在线| 小说图片视频综合网站| 午夜激情欧美在线| 精品国产三级普通话版| 国产激情欧美一区二区| 日日干狠狠操夜夜爽| 午夜两性在线视频| 99在线视频只有这里精品首页| 成人国产一区最新在线观看| 美女扒开内裤让男人捅视频| 99热只有精品国产| 好看av亚洲va欧美ⅴa在| 色精品久久人妻99蜜桃| 麻豆久久精品国产亚洲av| 国产视频内射| 亚洲激情在线av| 成人精品一区二区免费| 亚洲国产精品成人综合色| 动漫黄色视频在线观看| 最新在线观看一区二区三区| 宅男免费午夜| 亚洲中文日韩欧美视频| 亚洲av熟女| 精品一区二区三区视频在线观看免费| 国产精品美女特级片免费视频播放器 | 久久性视频一级片| 一区二区三区高清视频在线| 久久久国产成人精品二区| 国产91精品成人一区二区三区| 日日夜夜操网爽| 国产精品av久久久久免费| 白带黄色成豆腐渣| 在线观看免费午夜福利视频| 亚洲国产精品999在线| xxx96com| 欧美日韩乱码在线| 国产极品精品免费视频能看的| 中文字幕精品亚洲无线码一区| 丁香欧美五月| 精品国产乱子伦一区二区三区| 国产精品爽爽va在线观看网站| 欧美精品啪啪一区二区三区| 亚洲色图 男人天堂 中文字幕| 国产精品影院久久| 亚洲第一电影网av| 午夜福利免费观看在线| 成人国产综合亚洲| 久99久视频精品免费| 男人舔女人下体高潮全视频| 99久久精品一区二区三区| 亚洲国产精品成人综合色| 午夜视频精品福利| 欧美一级毛片孕妇| 天堂动漫精品| 又黄又爽又免费观看的视频| 久99久视频精品免费| 久久亚洲精品不卡| 可以在线观看毛片的网站| cao死你这个sao货| 麻豆久久精品国产亚洲av| 婷婷丁香在线五月| 天堂动漫精品| 亚洲国产精品合色在线| 亚洲人成电影免费在线| 美女免费视频网站| 啦啦啦韩国在线观看视频| 69av精品久久久久久| 日本免费一区二区三区高清不卡| 特大巨黑吊av在线直播| 国产一区二区在线av高清观看| 天天一区二区日本电影三级| 麻豆成人午夜福利视频| 日本 av在线| 看黄色毛片网站| www.www免费av| 午夜影院日韩av| 久久久久久久久免费视频了| av在线天堂中文字幕| 亚洲电影在线观看av| 久久草成人影院| 亚洲熟妇熟女久久| 国语自产精品视频在线第100页| 成人性生交大片免费视频hd| 久久精品国产亚洲av香蕉五月| 丁香六月欧美| а√天堂www在线а√下载| 两个人的视频大全免费| 亚洲男人的天堂狠狠| 90打野战视频偷拍视频| 免费av不卡在线播放| av片东京热男人的天堂| 18禁国产床啪视频网站| 欧美日本视频| 久久性视频一级片| 日韩欧美精品v在线| 怎么达到女性高潮| 久久亚洲真实| 国产男靠女视频免费网站| 精品久久久久久久久久免费视频| 91av网站免费观看| 国产高潮美女av| 国产精品永久免费网站| 亚洲av电影在线进入| 小说图片视频综合网站| 又爽又黄无遮挡网站| 长腿黑丝高跟| 日韩精品中文字幕看吧| 亚洲人成伊人成综合网2020| 久久久久久久久久黄片| 国产毛片a区久久久久| 非洲黑人性xxxx精品又粗又长| 国内揄拍国产精品人妻在线| 黑人操中国人逼视频| 午夜精品在线福利| 国产高清激情床上av| 欧美乱码精品一区二区三区| 男插女下体视频免费在线播放| 久久久国产欧美日韩av| 国产激情欧美一区二区| 免费看十八禁软件| 男女午夜视频在线观看| 日本一本二区三区精品| 亚洲国产精品成人综合色| 婷婷精品国产亚洲av在线| 午夜精品一区二区三区免费看| 一夜夜www| 欧美+亚洲+日韩+国产| 国产av一区在线观看免费| 制服丝袜大香蕉在线| 99久久久亚洲精品蜜臀av| 高清在线国产一区| 男女之事视频高清在线观看| 日韩人妻高清精品专区| 久久热在线av| 51午夜福利影视在线观看| 一区二区三区高清视频在线| 亚洲欧洲精品一区二区精品久久久| 亚洲九九香蕉| 999精品在线视频|