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

    WRF模式中不同參數(shù)化方案對(duì)華中地區(qū)一次強(qiáng)降水天氣過(guò)程模擬的影響

    2017-12-01 09:47:25莫欣岳
    關(guān)鍵詞:積云實(shí)況強(qiáng)降水

    杜 韜,莫欣岳,李 歡,張 鐳*

    (1.蘭州大學(xué)大氣科學(xué)學(xué)院,半干旱氣候變化教育部重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000;2.蘭州大學(xué)信息科學(xué)與工程學(xué)院,甘肅 蘭州 730000)

    WRF模式中不同參數(shù)化方案對(duì)華中地區(qū)一次強(qiáng)降水天氣過(guò)程模擬的影響

    杜 韜1,莫欣岳1,李 歡2,張 鐳1*

    (1.蘭州大學(xué)大氣科學(xué)學(xué)院,半干旱氣候變化教育部重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000;2.蘭州大學(xué)信息科學(xué)與工程學(xué)院,甘肅 蘭州 730000)

    【目的】利用中尺度數(shù)值模式WRF(V3.6.1)和NCEP的FNL資料,模擬2013年7月5-6日華中地區(qū)的一次強(qiáng)降水天氣過(guò)程?!痉椒ā客ㄟ^(guò)WRF模式中2種積云對(duì)流參數(shù)化方案和3種云微物理參數(shù)化方案的6種不同組合,考察不同微物理和積云對(duì)流參數(shù)化方案對(duì)此次降水過(guò)程模擬的影響,將模擬結(jié)果和24 h累計(jì)降水實(shí)況以及雷達(dá)回波觀(guān)測(cè)進(jìn)行對(duì)比分析。【結(jié)果】發(fā)現(xiàn)采用不同微物理和積云對(duì)流參數(shù)化方案的組合都能大致模擬出降水的落區(qū)和發(fā)展趨勢(shì),但模擬的強(qiáng)降水中心的位置和強(qiáng)度與實(shí)況相比存在一定偏差?!窘Y(jié)論】總體來(lái)說(shuō),KFN方案與Kessler方案的組合以及G3方案與Lin方案的組合的模擬效果較好。結(jié)合ETS評(píng)分,不同的參數(shù)化方案對(duì)不同等級(jí)降水模擬的效果不同,整體模擬效果比較穩(wěn)定的是WSM6方案和G3方案的組合。

    WRF模式;微物理參數(shù)化方案;強(qiáng)降水天氣過(guò)程;ETS評(píng)分

    【研究意義】降水是多種因子起作用的復(fù)雜物理過(guò)程,它涉及到云的微物理過(guò)程和宏觀(guān)動(dòng)力過(guò)程。利用高分辨率的中尺度數(shù)值天氣模式模擬和預(yù)報(bào)區(qū)域降水過(guò)程,具有不可替代的作用[1-3]。中尺度模式的參數(shù)化方案主要包括對(duì)流參數(shù)化過(guò)程方案、云微物理方案以及邊界層方案等,選擇不同的參數(shù)化方案對(duì)模擬結(jié)果的準(zhǔn)確度影響很大。其中云微物理過(guò)程是中尺度數(shù)值模式中最重要的非絕熱加熱物理過(guò)程之一,在決定大尺度大氣溫度、濕度場(chǎng)的垂直結(jié)構(gòu)中起著關(guān)鍵作用,也是降水預(yù)報(bào)的關(guān)鍵所在[4]。積云對(duì)流也是數(shù)值模式中不可缺少的物理過(guò)程,對(duì)一次強(qiáng)降水天氣過(guò)程的模擬有著重要的作用。WRF(Weather Research Forecast)模式系統(tǒng)是由美國(guó)國(guó)家大氣研究中心(NCAR)、國(guó)家大氣海洋局預(yù)報(bào)系統(tǒng)實(shí)驗(yàn)室、國(guó)家大氣環(huán)境研究中心(FSL,NCEP/NOAA)和俄克拉荷馬大學(xué)暴雨分析中心等多單位聯(lián)合發(fā)展起來(lái)的新一代非靜力平衡、高分辨率、科研和業(yè)務(wù)預(yù)報(bào)統(tǒng)一的中尺度預(yù)報(bào)和資料同化模式[5-6],在天氣預(yù)報(bào)、大氣化學(xué),以及區(qū)域氣候等方面有著廣泛應(yīng)用前景?!厩叭搜芯窟M(jìn)展】國(guó)內(nèi)外許多學(xué)者對(duì)WRF模式運(yùn)用于區(qū)域降水的模擬研究進(jìn)行了不同的研究。何宏讓等[7]用WRF模式的3種云微物理參數(shù)化方案模擬2010年舟曲特大暴雨,選擇不同的云微物理參數(shù)化方案對(duì)強(qiáng)降水強(qiáng)度的模擬有很大差別,kessler方案模擬的降水范圍和走向以及累積降水量與實(shí)況基本符合,而Lin方案和Morrison方案模擬的累積降水量與實(shí)況相差較大。王婷婷等[8]運(yùn)用WRF模式對(duì)三峽庫(kù)區(qū)的一次降水過(guò)程進(jìn)行了中尺度數(shù)值模擬,發(fā)現(xiàn)WSM3簡(jiǎn)單冰方案和Ferrier微物理方案在面平均降水量、降水分布等指標(biāo)上能更好的反映實(shí)況,而WSM5方案的模擬結(jié)果與實(shí)測(cè)值差異較大。朱格力等[9]用WRF模式的8種不同云微物理參數(shù)化方案完成對(duì)華南地區(qū)一次暴雨過(guò)程的數(shù)值模擬和性能分析,發(fā)現(xiàn)WSM3方案對(duì)小到大雨和大暴雨的模擬效果最好,對(duì)暴雨的模擬效果最差,WSM5方案對(duì)暴雨模擬效果最好,結(jié)合TS評(píng)分和誤差分析結(jié)果,整體模擬效果最好的是WSM5方案,最差的是Lin方案。Rajeevan et al.[10]通過(guò)分析WRF模式的4個(gè)不同的云微物理過(guò)程對(duì)印度南部雷雨事件的敏感性發(fā)現(xiàn),Thompson方案模擬的降水與實(shí)況符合最好,其他3種方案模擬結(jié)果均偏高。Jason等[11]用WSSR-88D雷達(dá)反射率資料推斷降水量,與WRF模式模擬的降水量進(jìn)行比較,發(fā)現(xiàn)在許多極端降水事件中,降水的量值、位置及局地降水頻率的最大值、降水頻率的日變化強(qiáng)度和總的空間分布模擬與觀(guān)測(cè)是相似的。黃海波等[12]用WRF模式研究了不同的云微物理參數(shù)化方案以及水平分辨率對(duì)降水預(yù)報(bào)效果的影響,發(fā)現(xiàn)模式水平分辨率的提高存在明顯的閾值?!颈狙芯壳腥朦c(diǎn)】云微物理參數(shù)化方案和積云對(duì)流參數(shù)化方案對(duì)不同類(lèi)型天氣系統(tǒng)的發(fā)展有不同的效果,而且在某些降水過(guò)程中,積云對(duì)流參數(shù)化方案幾乎完全不起作用。如何對(duì)WRF模式的參數(shù)化方案進(jìn)行優(yōu)化和選擇,有待于做進(jìn)一步系統(tǒng)的研究和評(píng)估?!緮M解決的關(guān)鍵問(wèn)題】本文利用WRF模式模擬2013年7月5日14:00至6日14:00華中地區(qū)一次強(qiáng)降水過(guò)程,選擇3種云微物理參數(shù)化方案和2種積云對(duì)流參數(shù)化方案相組合成6種參數(shù)化方案對(duì)強(qiáng)降水天氣過(guò)程進(jìn)行數(shù)值模擬試驗(yàn),對(duì)比分析不同微物理和積云對(duì)流參數(shù)化方案對(duì)強(qiáng)降水模擬的影響,對(duì)華中地區(qū)如何合理選擇微物理和積云對(duì)流參數(shù)化方案進(jìn)行初步探討。

    1 材料與方法

    1.1 天氣過(guò)程概述

    2013年7月5-6日,受高空槽、中低層切變線(xiàn)和西南低渦的共同影響下,我國(guó)華中地區(qū)出現(xiàn)了一次強(qiáng)降水天氣過(guò)程(圖1~2)。在500 hPa位勢(shì)高度場(chǎng)上(圖1a和圖2a),歐亞中高緯地區(qū)為兩槽一脊形勢(shì),在西伯利亞地區(qū)有一個(gè)大槽,隨著時(shí)間的推移,大槽不斷減弱,最終發(fā)展為一個(gè)橫槽,同時(shí)東北地區(qū)也有一個(gè)大槽,大槽不斷東移減弱,長(zhǎng)江地區(qū)存在一個(gè)南支槽,華中地區(qū)位于槽前脊后的位置,有上升氣流,有利于對(duì)流的發(fā)展,6日雖然華中地區(qū)的大槽有所減弱,但仍處于槽前脊后的位置(圖2a),表明之后降水強(qiáng)度有所減弱,但仍有降水發(fā)生;在低層,即在700(或850)hPa位勢(shì)高度場(chǎng)上,存在一個(gè)東北大槽,大槽不斷東移減弱,同時(shí)江淮地區(qū)有切變線(xiàn),偏北氣流和來(lái)自西南地區(qū)的水汽在華中地區(qū)發(fā)生輻合,有利于降水的產(chǎn)生;在海平面氣壓場(chǎng)(圖2d)上,貝加爾湖附近有一個(gè)低壓系統(tǒng),江淮地區(qū)應(yīng)有冷鋒存在,華中地區(qū)處于東北低壓的控制范圍下,從5日到6日,可看出低壓強(qiáng)度不斷增強(qiáng),表明水汽輻合上升增強(qiáng),有利于強(qiáng)降水的產(chǎn)生和發(fā)展。

    (a)500 hPa位勢(shì)高度場(chǎng),(b)700 hPa 位勢(shì)高度場(chǎng),(c)850 hPa位勢(shì)高度場(chǎng),(d)海平面氣壓場(chǎng)(a) 500 hPa geopotential height field, (b) 700 hPa geopotential height field, (c)850 hPa geopotential height field, (d) Sea level pressure field圖1 2013年7月5日8:00天氣形勢(shì)Fig.1 The weather chart at 8:00 on 5 July 2013

    (a)500 hPa位勢(shì)高度場(chǎng),(b)700 hPa 位勢(shì)高度場(chǎng),(c)850 hPa位勢(shì)高度場(chǎng),(d)海平面氣壓場(chǎng)(a) 500 hPa geopotential height field, (b) 700hPa geopotential height field, (c)850 hPa geopotential height field, (d) Sea level pressure field圖2 2013年7月6日8:00天氣形勢(shì)圖Fig.2 The weather chart at 8:00 on 6 July 2013.

    1.2 試驗(yàn)設(shè)計(jì)

    本文采用WRF模式V3.6.1版本,模式選用單向嵌套式網(wǎng)格,試驗(yàn)區(qū)域中心定為(100°E,32.5°N),水平格距為27 km,格點(diǎn)數(shù)為296×226,垂直方向分為30層,模式層頂為50 hPa。模式的初始和側(cè)邊界場(chǎng)資料來(lái)自NCEP的FNL資料,水平分辨率為1°×1°,時(shí)間間隔為6 h,采用Lambert地圖投影[13-16]。為了討論不同微物理和積云對(duì)流參數(shù)化方案在同一分辨率下對(duì)一次強(qiáng)降水的影響,選用RRTM長(zhǎng)波輻射方案、Dudhia短波輻射方案、YSU邊界層方案和Noah陸面過(guò)程方案,同時(shí)采用3種云微物理參數(shù)化方案和兩種對(duì)流參數(shù)化方案的6種組合(表1)。

    2 結(jié)果與分析

    2.1 雷達(dá)回波拼圖的模擬對(duì)比

    對(duì)于一些降水系統(tǒng)(如MCS),從初始對(duì)流系統(tǒng)到可分辨尺度系統(tǒng)轉(zhuǎn)變的正確模擬是再現(xiàn)許多觀(guān)測(cè)特征和動(dòng)力特征的關(guān)鍵所在[17]。從圖3可以看出,在7月5日14:00(圖3a),雷達(dá)回波帶呈東-西走向,江蘇南部和安徽中部以及湖北、四川局部地區(qū)降水強(qiáng)度大,降水范圍較廣;到5日20:00 (圖3b),降雨帶范圍減小,同時(shí)強(qiáng)降水區(qū)域位于江蘇與安徽的南部地區(qū),這表明強(qiáng)降水的中心位置向東轉(zhuǎn)移,降水中心范圍減小;到6日2:00(圖3c),在江蘇南部、安徽中南部以及安徽武漢交界處有強(qiáng)降水,表明強(qiáng)降水中心范圍有所擴(kuò)大,同時(shí)降水分布范圍進(jìn)一步縮??;到 6日8:00 (圖3d),強(qiáng)降水中心位于江蘇南部、安徽中部以及湖南、湖北交界處,強(qiáng)降水范圍進(jìn)一步增大;到6日14:00 (圖3e),強(qiáng)降水中心區(qū)域減少,強(qiáng)降水主要位于安徽和湖南交界地帶。

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

    (a)14:00,(b)20:00,(c) 02:00,(d) 08:00,(e) 14:00 (a)at 14:00 on 5, (b) at 20:00 on 5, (c) at 02:00 on 6,(d) at 08:00 on 6,(e) at 14:00 on 6圖3 2013年7月5-6日的累積降水量實(shí)況全國(guó)雷達(dá)拼圖Fig.3 Cumulative precipitation realistic national radar puzzle from 5 to 6 July 2013

    由圖4可知,將6種不同微物理過(guò)程模擬的20:00全國(guó)雷達(dá)拼圖與實(shí)況圖對(duì)比,6種方案都大致模擬出了雷達(dá)回波的位置及走向,方案1和方案2模擬的回波范圍較實(shí)況偏少,對(duì)回波中心位置的模擬也有較大偏差,總體來(lái)看,方案5和方案6的模擬效果較好,都大致模擬出了雷達(dá)回波帶的分布范圍,實(shí)況中的回波中心位于江蘇以及安徽南部地區(qū),模擬出的回波中心位置位于湖南地區(qū),同時(shí)模擬的回波強(qiáng)度偏低。

    綜合以上分析可知,方案5和方案6的模擬效果最好,基本模擬出了此次強(qiáng)降水的雨帶走向以及降水范圍,對(duì)降水過(guò)程的連續(xù)過(guò)程模擬效果較好,其次是方案4和方案3,這表明在相同的對(duì)流參數(shù)化方案下,云微物理方案中WSM6方案和Lin方案對(duì)強(qiáng)降水的模擬效果較好,Kessler方案的模擬效果一般。

    2.2 24 h累積降水量的模擬對(duì)比

    由圖5可知,降水系統(tǒng)呈東-西走向,降水主要集中在江蘇、安徽、湖北南部地區(qū)以及湖南北部地區(qū),貴州地區(qū)也有少量降水,其中強(qiáng)降水中心位于安徽南部地區(qū),可達(dá)到200 mm左右;次強(qiáng)降水中心位于湖北南部地區(qū),可達(dá)到160 mm左右;在湖南和貴州局部地區(qū),降水強(qiáng)度可達(dá)到100 mm左右。

    對(duì)比圖5~6可以發(fā)現(xiàn),6個(gè)方案基本都能模擬出本次降水的主要雨帶及其走向,但降水區(qū)東西向區(qū)域相對(duì)實(shí)況明顯伸長(zhǎng)。具體來(lái)看,模式結(jié)果都較好的模擬出了位于安徽南部的強(qiáng)降水中心位置,沒(méi)有模擬出湖北南部的次強(qiáng)降水中心,方案1、方案3 和方案4模擬的強(qiáng)降水中心位置與強(qiáng)度與實(shí)況基本吻合;方案2模擬的強(qiáng)降水中心強(qiáng)度偏弱;方案5和方案6模擬出強(qiáng)降水中心位于湖北和安徽的交界處,強(qiáng)降水中心范圍相對(duì)實(shí)況有偏差。

    綜合以上分析可知,對(duì)于華中地區(qū)的一次強(qiáng)降水的模擬,各個(gè)方案均模擬出了降水的落區(qū),方案1和方案4的模擬效果與實(shí)況基本相符合,其次是方案6和方案3的模擬效果較好,這表明KFN與Kessler的組合方案以及G3與Lin的組合方案模擬的降水區(qū)走向以及24 h累積降水量與實(shí)況基本吻合,而Kessler方案與G3方案的組合所模擬的累積降水量與實(shí)況相差較大,不能很好的模擬出強(qiáng)降水中心的位置及強(qiáng)度。

    2.3 ETS評(píng)分檢驗(yàn)

    ETS評(píng)分方法是目前氣象部門(mén)檢驗(yàn)數(shù)值預(yù)報(bào)產(chǎn)品并對(duì)降水預(yù)報(bào)性能進(jìn)行評(píng)估的方法之一,采用ETS評(píng)分來(lái)對(duì)各組參數(shù)化方案數(shù)值模擬得到的降雨進(jìn)行檢驗(yàn)。

    從表2可以看出,不同微物理和積云對(duì)流參數(shù)化方案組合對(duì)不同等級(jí)降水的模擬效果不同。對(duì)于小雨等級(jí)的降水,方案6的模擬效果最好;對(duì)于中雨等級(jí)的降水,也是方案6的模擬效果最好;對(duì)于大雨等級(jí)的降水,方案1和方案4的模擬效果較好,其次是方案6;對(duì)于暴雨量級(jí)的降水,方案4的模擬效果較好,其次是方案1;對(duì)于大暴雨量級(jí)的降水,方案6的模擬效果最好,其次是方案4。總體來(lái)說(shuō),方案6對(duì)各個(gè)等級(jí)降水的模擬效果最好,模擬結(jié)果比較穩(wěn)定,即WSM6方案和G3方案的組合的模擬效果較好,其次是方案4,方案3和方案5的模擬效果也較好,方案1對(duì)大雨和暴雨的模擬效果很好,方案3對(duì)中雨的模擬效果較好。

    (a)方案1,(b)方案2,(c)方案3,(d)方案4,(e)方案5,(f)方案6(a) Scheme 1, (b) Scheme 2,(c) Scheme 3,(d) Scheme 4,(e) Scheme 5,(f) Scheme 6圖4 6種不同微物理過(guò)程模擬的7月5日20:00全國(guó)雷達(dá)拼圖(單位:dBZ)Fig.4 National radar puzzle in six different microphysical process simulations at 20:00 on 5 July (unit: dBZ)

    圖5 2013年7月5-6日的24 h累積降水量實(shí)況圖(單位:mm)Fig.5 The 24-hour cumulative precipitation from 5 to 6 July 2013

    綜合以上分析可知,對(duì)于大雨和暴雨等強(qiáng)降水的模擬,積云對(duì)流參數(shù)化方案對(duì)強(qiáng)降水模擬的影響大于云微物理參數(shù)化方案對(duì)強(qiáng)降水的影響。同時(shí)可知,從2種積云參數(shù)化方案的模擬結(jié)果評(píng)分來(lái)看,并沒(méi)有最優(yōu)的積云參數(shù)化方案,從微物理過(guò)程方案來(lái)比較,WSM6方案和Lin方案的總體ETS評(píng)分較好,Kessler方案的總體ETS評(píng)分較低。

    (a)方案1,(b)方案2,(c)方案3,(d)方案4,(e)方案5,(f)方案6(a)scheme 1, (b) scheme 2,(c) scheme 3,(d) scheme 4,(e) scheme 5,(f) scheme 6圖6 6種不同微物理過(guò)程模擬降水的空間分布(單位:mm)Fig.6 Spatial distribution of simulated precipitation in six different microphysical processes

    降水等級(jí)不同微物理參數(shù)化方案123456小雨0.2630.1910.3070.3300.3300.353中雨0.2610.4010.2070.1800.1800.401大雨0.5880.2170.3480.4490.3480.415暴雨0.4020.2240.3790.4250.3790.334大暴雨0.0410.1250.1670.2530.2100.296

    3 討 論

    為比較WRF模式的不同微物理和積云對(duì)流參數(shù)化方案對(duì)一次強(qiáng)降水天氣過(guò)程模擬結(jié)果的影響,利用6種不同微物理和積云對(duì)流參數(shù)化方案的組合,對(duì)2013年7月5日14:00~6日14:00華中地區(qū)一次強(qiáng)降水天氣過(guò)程進(jìn)行模擬。

    (1)通過(guò)對(duì)比分析模擬雷達(dá)回波拼圖和實(shí)況圖,發(fā)現(xiàn)采用不同微物理和積云對(duì)流參數(shù)化方案的組合都能大致模擬出降水的落區(qū)和發(fā)展趨勢(shì),但模擬的強(qiáng)降水中心的位置和強(qiáng)度與實(shí)況相比存在一定偏差,總的來(lái)說(shuō),WSM6、Lin方案這2種微物理參數(shù)化方案與KFN、G3這2種積云對(duì)流參數(shù)化方案搭配所得到的方案較好的模擬出了此次強(qiáng)降水系統(tǒng)的發(fā)生發(fā)展過(guò)程,與實(shí)況較接近。

    (2)WRF模式中不同的微物理和積云對(duì)流參數(shù)化方案的選擇對(duì)強(qiáng)降水落區(qū)模擬的影響不大,但是對(duì)強(qiáng)降水中心區(qū)域以及強(qiáng)度的模擬有較大差別,模擬結(jié)果表明KFN與Kessler的組合方案以及G3與Lin的組合方案模擬的降水走向和24 h累積降水量與實(shí)況基本吻合,而Kessler方案與G3方案的組合所模擬的累積降水量和強(qiáng)降水中心的位置及強(qiáng)度與實(shí)況存在偏差。

    (3)從試驗(yàn)結(jié)果的ETS評(píng)分可以知道,對(duì)于小雨等級(jí)的降水,WSM6與G3組合方案的模擬效果最好;對(duì)于中雨等級(jí)的降水,WSM6、Lin方案這2種微物理方案搭配G3方案的模擬效果最好;對(duì)于大雨等級(jí)的降水,Kessler方案搭配KFN方案組合的模擬效果較好;對(duì)于暴雨量級(jí)的降水,Lin方案搭配G3方案的模擬效果較好;對(duì)于大暴雨量級(jí)的降水,WSM6方案搭配G3方案的模擬效果最好。

    4 結(jié) 論

    由模擬結(jié)果可知,在相同的積云對(duì)流參數(shù)化方案下,微物理參數(shù)化方案中對(duì)各個(gè)降水等級(jí)的模擬效果較穩(wěn)定的是WSM6方案和Lin方案,Kessler方案對(duì)不同等級(jí)降水模擬的效果相差較大。同時(shí)對(duì)于相同的云微物理參數(shù)化方案,選取不同的積云對(duì)流參數(shù)化方案對(duì)模擬結(jié)果也有影響。

    本文的結(jié)論只是WRF模式對(duì)一次強(qiáng)降水天氣過(guò)程的模擬結(jié)果,由于不同天氣過(guò)程的主要影響因素不同,模擬的準(zhǔn)確性除了受不同微物理過(guò)程影響外,還與模式的水平分辨率的高低、是否選取合適的嵌套方案和其他物理參數(shù)化方案有關(guān),因此,結(jié)論是否適用于其他強(qiáng)降水過(guò)程,還有待于進(jìn)一步檢驗(yàn)。

    [1]CraigK,Bornstein R.MM5 Simulation of Urban Induced Convective Precipitation overAtlanta [A].Preprints, Fourth Conference on the Urban Environment[C].Norfolk, American Meteorology Society,2002:5-6.

    [2]王 歡,倪允琪.2003年淮河汛期一次中尺度強(qiáng)暴雨過(guò)程的診斷分析和數(shù)值模擬研究[J].氣象學(xué)報(bào),2006,64(6): 734-742.

    [3]王 婷,吳池勝,馮瑞權(quán).2005年6月廣東一次暴雨過(guò)程的中尺度對(duì)流系統(tǒng)的數(shù)值研究[J].大氣科學(xué),2008,32(1):184-196.

    [4]陳德輝,胡志晉,徐大海,等.CAMS大氣數(shù)值預(yù)報(bào)模式系統(tǒng)研究[M].氣象科技,2004:35-66.

    [5]Emanuel Kerry A,Zivkovic-Rothman. Development and evaluation of a convection scheme for use in climate models[J].Journal of the Atmospheric Science,1999,56:1 766-1 782.

    [6]Skamarock W, Klemp J, Dudhia J, et al. A description of the advanced research WRF version 2[J].NCAR Technical Note, NCAR/TN-47.5+STR,2005.

    [7]李安泰,何宏讓.不同云微物理參數(shù)化方案對(duì)舟曲“8.8”暴雨過(guò)程的影響[J].氣象與減災(zāi)研究,2011,34(3):9-11.

    [8]王婷婷,周建中,曾小凡,等.WRF模式云微物理過(guò)程對(duì)三峽庫(kù)區(qū)降水預(yù)報(bào)的影響[J].水電能源科學(xué),2013,31(6):17-20.

    [9]朱格利,林萬(wàn)濤,曹艷華.用WRF模式中不同云微物理參數(shù)化方案對(duì)華南一次暴雨過(guò)程的數(shù)值模擬和性能分析[J].大氣科學(xué),2014,38(3):513-523.

    [10]Rajeevan M, Kesarkar A, Thampi S B, et al.Sensitivity of WRF cloud microphysics to simulations of a severe thunder storm event over Southeast India[J].Ann.Geophys,2010,28:603-619.

    [11]Jason C Knievel. The Diurnal Mode of Summer Rainfall Across the Conterminous United States in 10-km Simulation by the WRF Model[C]//84th AMS Annual Meeting , U.S.A.: Seattle ,2004:10-15.

    [12]黃海波,陳春艷,朱雯娜.WRF模式不同云微物理參數(shù)化方案對(duì)降水預(yù)報(bào)效果的影響[J].氣象科技,2011,39(5):529-536.

    [13]Lin Y L, Farley R, Orville H D. Bulk parameterization of the snow field in a cloud model[J]. J. Climate Appl. Meteor.,1983,22,1065-1029.

    [14]閆之輝,鄧蓮堂.WRF模式中的微物理過(guò)程及其預(yù)報(bào)對(duì)比試驗(yàn)[J].沙漠與綠洲氣象,2007,1(6):1-6.

    [15]譚冶冰.不同陸面過(guò)程及邊界層參數(shù)化方案的城市邊界層模擬實(shí)驗(yàn)[J].咸陽(yáng)師范學(xué)院學(xué)報(bào),2010,25(2):56-61.

    [16]胡向軍,陶健紅.WRF模式物理過(guò)程參數(shù)化方案簡(jiǎn)介[J].甘肅科技,2008,20(1):73-75.

    [17]侯瑞欽,程麟生.積云對(duì)流參數(shù)化方案對(duì)梅雨鋒暴雨過(guò)程模擬的影響[J].氣象科技,2007,35(4):458-463.

    (責(zé)任編輯 李 潔)

    InfluenceofSimulatingHeavyRainfallProcessinCentralChinabyUsingDifferentMicrophysicsandCumulusParameterizationSchemesinWRFModel

    DU Tao1, MO Xin-yue1, LI Huan2, ZHANG Lei1*

    (1.College of Atmospheric Sciences, Lanzhou University, Key Laboratory for Semi-Arid Climate Change of the Ministry of Education, Gansu Lanzhou 730000, China;2.College of Information Science and Engineering,Lanzhou University,Gansu Lanzhou 730000,China)

    【Objective】In this paper, a heavy rainfall occurred in central China during June 5-6 in 2013 was investigated by using the mesoscale numerical model WRF (V3.6.1) and FNL NCEP data.【Method】Based on six different combinations of two cumulus convective parameterization schemes and three cloud microphysics parameterization schemes, the influences of different parameterization schemes on the simulation of the precipitation process were studied, by comparing the observed 24-hour accumulated precipitation and radar echo with simulations.【Result】Results showed that the movements and the distribution of precipitation could roughly be simulated by using different parameterization schemes. However, relative to the observations, the simulated location and intensity of heavy precipitation center has some basis.【Conclusion】In general, the simulation results of the combination of KFN scheme and Kessler scheme and the combination of Lin scheme and G3 scheme are pretty well. Combining the equitable threat score, the simulation effects of different parameterization schemes are different for different levels of precipitation and the combination of WSM6 scheme and G3 scheme has stable simulation effects.

    WRF model; Microphysics parameterization schemes; Heavy rainfall process; Equitable threat score

    P445

    A

    1001-4829(2017)11-2583-07

    10.16213/j.cnki.scjas.2017.11.033

    2016-12-13

    國(guó)家自然科學(xué)基金面上項(xiàng)目(41475008);國(guó)家自然科學(xué)基金創(chuàng)新研究群體科學(xué)基金(41521004);國(guó)家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(41605005);蘭州大學(xué)中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金(lzujbky-2016-k06)

    杜 韜(1993-),男,甘肅蘭州人,碩士研究生,研究方向:大氣物理學(xué)與大氣環(huán)境,Tel:15117125544,E-mail:545965475@qq.com,*為通訊作者:張 鐳(1960-),男,教授,博士生導(dǎo)師,研究方向:大氣物理學(xué)與大氣環(huán)境,Tel:13893615553,E-mail:zhanglei@lzu.edu.cn。

    猜你喜歡
    積云實(shí)況強(qiáng)降水
    2020年江淮地區(qū)夏季持續(xù)性強(qiáng)降水過(guò)程分析
    一次東移型西南低渦引發(fā)的強(qiáng)降水診斷分析
    WRF模式中積云對(duì)流參數(shù)化方案對(duì)南海土臺(tái)風(fēng)“Ryan”模擬的影響研究
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實(shí)——宜君鄉(xiāng)村教育實(shí)況
    積云——晴天的使者
    Warning signs of an approaching hurricane
    天舟一號(hào)貨運(yùn)飛船發(fā)射實(shí)況掠影
    可愛(ài)潮咖們的獨(dú)門(mén)彩妝實(shí)況直播
    四川盆地西南部短時(shí)強(qiáng)降水天氣特征分析
    牽掛
    999久久久精品免费观看国产| 99国产精品一区二区三区| 国产亚洲欧美98| 精品无人区乱码1区二区| 别揉我奶头~嗯~啊~动态视频| 黄片大片在线免费观看| a级毛片在线看网站| 男女之事视频高清在线观看| 亚洲av日韩精品久久久久久密| 男女之事视频高清在线观看| 最近在线观看免费完整版| 国产亚洲欧美98| 伦理电影免费视频| 国产av麻豆久久久久久久| 看免费av毛片| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久久电影 | 久99久视频精品免费| 亚洲成av人片免费观看| 国产亚洲精品久久久com| 禁无遮挡网站| 熟女人妻精品中文字幕| 国产视频内射| svipshipincom国产片| 久久热在线av| 麻豆一二三区av精品| 国产伦一二天堂av在线观看| 岛国在线观看网站| 一级黄色大片毛片| av在线蜜桃| 国产av在哪里看| 亚洲片人在线观看| 床上黄色一级片| 日韩av在线大香蕉| 国产aⅴ精品一区二区三区波| 黄色日韩在线| 丰满的人妻完整版| 免费人成视频x8x8入口观看| 中文在线观看免费www的网站| 国产精品1区2区在线观看.| 一a级毛片在线观看| 精品免费久久久久久久清纯| 欧美中文综合在线视频| 女同久久另类99精品国产91| 天堂动漫精品| 午夜激情福利司机影院| 不卡av一区二区三区| 村上凉子中文字幕在线| 99久久无色码亚洲精品果冻| 欧美中文综合在线视频| 国产一区二区三区视频了| 亚洲,欧美精品.| 一本久久中文字幕| 亚洲av中文字字幕乱码综合| 久久久国产欧美日韩av| 免费看美女性在线毛片视频| 欧美极品一区二区三区四区| 国产精品亚洲一级av第二区| 国产伦精品一区二区三区视频9 | 麻豆久久精品国产亚洲av| 一个人看的www免费观看视频| 亚洲av五月六月丁香网| 精品不卡国产一区二区三区| 国内少妇人妻偷人精品xxx网站 | 99久国产av精品| 九色国产91popny在线| 国产人伦9x9x在线观看| 啦啦啦观看免费观看视频高清| 老司机深夜福利视频在线观看| 动漫黄色视频在线观看| 欧美又色又爽又黄视频| 亚洲欧美日韩无卡精品| www.自偷自拍.com| 国产一区在线观看成人免费| 丝袜人妻中文字幕| xxx96com| 神马国产精品三级电影在线观看| 岛国在线免费视频观看| 熟女少妇亚洲综合色aaa.| 成年版毛片免费区| 九九久久精品国产亚洲av麻豆 | netflix在线观看网站| 婷婷精品国产亚洲av在线| 久久久精品欧美日韩精品| 两人在一起打扑克的视频| 欧美乱色亚洲激情| 久久精品国产清高在天天线| 欧美日韩一级在线毛片| 熟妇人妻久久中文字幕3abv| 窝窝影院91人妻| 午夜久久久久精精品| 成人特级黄色片久久久久久久| 精品乱码久久久久久99久播| 99久久国产精品久久久| 国产亚洲精品一区二区www| 亚洲欧美日韩无卡精品| 高清毛片免费观看视频网站| 色综合婷婷激情| 2021天堂中文幕一二区在线观| 欧美日韩乱码在线| 国产精品久久电影中文字幕| 身体一侧抽搐| 制服丝袜大香蕉在线| 国产v大片淫在线免费观看| 曰老女人黄片| 亚洲国产欧洲综合997久久,| 色视频www国产| 免费看光身美女| 久久这里只有精品中国| 久久久国产成人免费| 一区二区三区国产精品乱码| 免费看a级黄色片| 日本a在线网址| 99热这里只有精品一区 | 久久久水蜜桃国产精品网| 无限看片的www在线观看| 97碰自拍视频| 亚洲av第一区精品v没综合| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 亚洲精品粉嫩美女一区| 免费在线观看成人毛片| 国产免费男女视频| 一区福利在线观看| 99在线视频只有这里精品首页| 怎么达到女性高潮| 午夜福利18| 亚洲国产精品成人综合色| 日本三级黄在线观看| av福利片在线观看| 91麻豆精品激情在线观看国产| 午夜免费激情av| 成人18禁在线播放| 国产免费男女视频| 国内少妇人妻偷人精品xxx网站 | 性欧美人与动物交配| 一进一出抽搐动态| 久久久久性生活片| 亚洲精品乱码久久久v下载方式 | 少妇的逼水好多| 又紧又爽又黄一区二区| 美女午夜性视频免费| 色综合站精品国产| 最近最新免费中文字幕在线| av片东京热男人的天堂| 国产 一区 欧美 日韩| 国产伦在线观看视频一区| 久久久久久人人人人人| 国产高清视频在线观看网站| 日本免费a在线| 极品教师在线免费播放| 亚洲av第一区精品v没综合| 一本一本综合久久| 伦理电影免费视频| 亚洲一区二区三区不卡视频| 好看av亚洲va欧美ⅴa在| 99国产精品一区二区三区| 婷婷精品国产亚洲av| 亚洲成人免费电影在线观看| 久久久久久人人人人人| 日韩欧美国产一区二区入口| 久久精品夜夜夜夜夜久久蜜豆| 少妇丰满av| 国产主播在线观看一区二区| 一a级毛片在线观看| 国产主播在线观看一区二区| 一二三四社区在线视频社区8| 琪琪午夜伦伦电影理论片6080| 18禁美女被吸乳视频| 岛国在线免费视频观看| 亚洲av电影在线进入| 国产三级黄色录像| 欧美黑人欧美精品刺激| 久久久久久国产a免费观看| 99国产综合亚洲精品| 国产成人影院久久av| 校园春色视频在线观看| 亚洲av成人一区二区三| 男人舔女人的私密视频| 国产精品av视频在线免费观看| 国产精品精品国产色婷婷| 美女 人体艺术 gogo| 亚洲在线观看片| 国产精品亚洲一级av第二区| 欧美午夜高清在线| 91在线精品国自产拍蜜月 | 国产亚洲av高清不卡| 99国产精品一区二区三区| 国产精品乱码一区二三区的特点| 无限看片的www在线观看| 麻豆成人av在线观看| 999精品在线视频| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久大精品| 母亲3免费完整高清在线观看| 午夜福利成人在线免费观看| 国产亚洲精品一区二区www| 国产亚洲av高清不卡| 97人妻精品一区二区三区麻豆| 亚洲av电影在线进入| 国产久久久一区二区三区| 悠悠久久av| 亚洲欧美日韩卡通动漫| 国产精品野战在线观看| 欧美午夜高清在线| 亚洲欧美精品综合久久99| 亚洲中文字幕一区二区三区有码在线看 | 日韩高清综合在线| 欧美日韩瑟瑟在线播放| 久久久久久大精品| 怎么达到女性高潮| 国产精品一及| 老司机福利观看| 12—13女人毛片做爰片一| 国产高清有码在线观看视频| 国产亚洲精品综合一区在线观看| 中文字幕人成人乱码亚洲影| 最近最新中文字幕大全免费视频| 国产69精品久久久久777片 | 热99在线观看视频| 日韩三级视频一区二区三区| 90打野战视频偷拍视频| 亚洲成av人片免费观看| 精品99又大又爽又粗少妇毛片 | 成人欧美大片| 手机成人av网站| 精品电影一区二区在线| а√天堂www在线а√下载| 亚洲在线自拍视频| 他把我摸到了高潮在线观看| 久久久久国产一级毛片高清牌| 久久久久久久久久黄片| 中出人妻视频一区二区| 又粗又爽又猛毛片免费看| 麻豆av在线久日| 欧美日本视频| 国产午夜精品久久久久久| 蜜桃久久精品国产亚洲av| 亚洲 国产 在线| 欧美激情在线99| 看片在线看免费视频| 黄色女人牲交| 中出人妻视频一区二区| 精品人妻1区二区| 禁无遮挡网站| 在线a可以看的网站| 久久中文字幕一级| 小蜜桃在线观看免费完整版高清| 在线永久观看黄色视频| 99久久综合精品五月天人人| 国产三级在线视频| 国产麻豆成人av免费视频| 99视频精品全部免费 在线 | 老汉色∧v一级毛片| 免费搜索国产男女视频| 两性午夜刺激爽爽歪歪视频在线观看| 噜噜噜噜噜久久久久久91| 久久久久久久久免费视频了| 黄色女人牲交| 亚洲男人的天堂狠狠| 此物有八面人人有两片| 操出白浆在线播放| 亚洲色图av天堂| 亚洲国产欧洲综合997久久,| 免费看光身美女| 色在线成人网| 狠狠狠狠99中文字幕| 亚洲av电影不卡..在线观看| 又黄又爽又免费观看的视频| 丰满的人妻完整版| 一本一本综合久久| 成人一区二区视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲精品中文字幕一二三四区| 黄色丝袜av网址大全| 国内精品美女久久久久久| or卡值多少钱| 99久久综合精品五月天人人| 久久国产精品影院| 五月玫瑰六月丁香| x7x7x7水蜜桃| 欧美不卡视频在线免费观看| 国产精品免费一区二区三区在线| 久99久视频精品免费| 国内精品一区二区在线观看| 美女黄网站色视频| 精品久久久久久久久久久久久| 色在线成人网| 日本黄色片子视频| 午夜福利在线在线| 99热6这里只有精品| 不卡一级毛片| 一本一本综合久久| 亚洲一区二区三区色噜噜| 麻豆av在线久日| 久久久久国产精品人妻aⅴ院| 90打野战视频偷拍视频| 小蜜桃在线观看免费完整版高清| 成年女人永久免费观看视频| 免费在线观看日本一区| 1024香蕉在线观看| 熟女少妇亚洲综合色aaa.| 色在线成人网| 又爽又黄无遮挡网站| 国产精品久久久久久亚洲av鲁大| 国产成人一区二区三区免费视频网站| 日韩人妻高清精品专区| 别揉我奶头~嗯~啊~动态视频| 夜夜看夜夜爽夜夜摸| 老鸭窝网址在线观看| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 亚洲国产欧美网| 变态另类丝袜制服| 日日夜夜操网爽| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 一a级毛片在线观看| 女人被狂操c到高潮| 人人妻,人人澡人人爽秒播| 可以在线观看的亚洲视频| 精品一区二区三区视频在线观看免费| 两性夫妻黄色片| 三级男女做爰猛烈吃奶摸视频| 成人午夜高清在线视频| 手机成人av网站| 91在线精品国自产拍蜜月 | 久久精品91蜜桃| a级毛片在线看网站| 国产精品国产高清国产av| 制服丝袜大香蕉在线| 91九色精品人成在线观看| 欧美性猛交黑人性爽| 成人一区二区视频在线观看| 国产一区二区在线av高清观看| 人妻丰满熟妇av一区二区三区| 美女黄网站色视频| 天堂网av新在线| 丰满人妻一区二区三区视频av | 成人一区二区视频在线观看| 女生性感内裤真人,穿戴方法视频| 天堂网av新在线| 美女高潮喷水抽搐中文字幕| 夜夜躁狠狠躁天天躁| 亚洲av美国av| 99热精品在线国产| 成人特级av手机在线观看| 脱女人内裤的视频| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 亚洲中文av在线| 免费一级毛片在线播放高清视频| 90打野战视频偷拍视频| 国产熟女xx| 俄罗斯特黄特色一大片| 免费在线观看视频国产中文字幕亚洲| 国产成人精品久久二区二区91| 成人高潮视频无遮挡免费网站| 一个人免费在线观看电影 | 国产成人影院久久av| 国产精品亚洲一级av第二区| 两性夫妻黄色片| 精品国产超薄肉色丝袜足j| 午夜福利视频1000在线观看| 1000部很黄的大片| 午夜激情欧美在线| 最好的美女福利视频网| 性欧美人与动物交配| 日韩欧美精品v在线| 国产成人欧美在线观看| avwww免费| 久久精品国产综合久久久| 九色成人免费人妻av| 一级毛片精品| 国产精品精品国产色婷婷| 国产熟女xx| 免费观看人在逋| 日韩国内少妇激情av| 狠狠狠狠99中文字幕| 又黄又爽又免费观看的视频| 欧美xxxx黑人xx丫x性爽| 国产精品亚洲一级av第二区| 好男人在线观看高清免费视频| 免费在线观看影片大全网站| 在线免费观看不下载黄p国产 | 精品久久久久久久久久久久久| 丰满的人妻完整版| 国产成人av激情在线播放| 亚洲18禁久久av| 国产探花在线观看一区二区| 国产一区二区三区在线臀色熟女| 超碰成人久久| 亚洲精华国产精华精| 三级男女做爰猛烈吃奶摸视频| 俺也久久电影网| 国产一级毛片七仙女欲春2| 好男人在线观看高清免费视频| 日韩欧美 国产精品| 亚洲成人久久爱视频| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 亚洲一区二区三区色噜噜| 成人av在线播放网站| 男人舔奶头视频| 日本黄色片子视频| 久久精品夜夜夜夜夜久久蜜豆| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 国内精品久久久久精免费| 一区福利在线观看| 亚洲人成伊人成综合网2020| 精品国内亚洲2022精品成人| 一区二区三区国产精品乱码| 欧美又色又爽又黄视频| 久久精品影院6| e午夜精品久久久久久久| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 9191精品国产免费久久| 国产成人系列免费观看| 久99久视频精品免费| 国产精品1区2区在线观看.| 国产精品精品国产色婷婷| 国产高清三级在线| 亚洲成人中文字幕在线播放| bbb黄色大片| 亚洲中文字幕一区二区三区有码在线看 | 午夜福利在线观看吧| 亚洲av成人一区二区三| 久久精品影院6| 中文在线观看免费www的网站| 午夜福利欧美成人| x7x7x7水蜜桃| 每晚都被弄得嗷嗷叫到高潮| 91久久精品国产一区二区成人 | 最近最新中文字幕大全电影3| 在线观看一区二区三区| 最新在线观看一区二区三区| 国产精品综合久久久久久久免费| 十八禁人妻一区二区| 亚洲国产精品合色在线| a级毛片a级免费在线| 中文字幕久久专区| 亚洲av美国av| svipshipincom国产片| 巨乳人妻的诱惑在线观看| 97超级碰碰碰精品色视频在线观看| 久久久久久人人人人人| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲一级av第二区| 免费看十八禁软件| 精品福利观看| 久久久国产精品麻豆| 丝袜人妻中文字幕| 少妇熟女aⅴ在线视频| 久久中文字幕一级| 一个人看视频在线观看www免费 | 亚洲精品美女久久av网站| 精品久久久久久久久久久久久| 成人三级黄色视频| 午夜免费激情av| 老司机在亚洲福利影院| av在线天堂中文字幕| 国产aⅴ精品一区二区三区波| 国产乱人视频| 国产精华一区二区三区| 成年版毛片免费区| 亚洲九九香蕉| 国产精品电影一区二区三区| 国产精品久久久人人做人人爽| 欧美一级毛片孕妇| 国产av在哪里看| 熟女人妻精品中文字幕| 在线观看66精品国产| 成人午夜高清在线视频| 成人高潮视频无遮挡免费网站| 精品免费久久久久久久清纯| 日本黄色视频三级网站网址| 男插女下体视频免费在线播放| 亚洲色图av天堂| 成人国产综合亚洲| 级片在线观看| 怎么达到女性高潮| 在线免费观看的www视频| 免费搜索国产男女视频| 99国产精品一区二区三区| 国产一区二区在线av高清观看| 一区二区三区国产精品乱码| 日韩欧美免费精品| 国产乱人伦免费视频| 国产三级黄色录像| 欧美一级毛片孕妇| 小蜜桃在线观看免费完整版高清| 制服人妻中文乱码| 精品国产亚洲在线| 欧美日韩乱码在线| 亚洲在线自拍视频| 啦啦啦免费观看视频1| 久9热在线精品视频| 动漫黄色视频在线观看| 精品无人区乱码1区二区| 亚洲国产看品久久| 91老司机精品| 九九热线精品视视频播放| 一本综合久久免费| 久久精品91无色码中文字幕| 身体一侧抽搐| 在线免费观看不下载黄p国产 | 麻豆国产97在线/欧美| 成人永久免费在线观看视频| 亚洲精品久久国产高清桃花| 老熟妇仑乱视频hdxx| 中出人妻视频一区二区| 88av欧美| 久久精品亚洲精品国产色婷小说| 久久精品91无色码中文字幕| 此物有八面人人有两片| 人妻久久中文字幕网| 五月玫瑰六月丁香| 中文资源天堂在线| 久久久久国内视频| 成人国产一区最新在线观看| 日本三级黄在线观看| 19禁男女啪啪无遮挡网站| 一个人看视频在线观看www免费 | 两个人看的免费小视频| 亚洲性夜色夜夜综合| 一级黄色大片毛片| 很黄的视频免费| 黄色女人牲交| 人人妻人人看人人澡| 一个人观看的视频www高清免费观看 | 激情在线观看视频在线高清| 两性夫妻黄色片| 色在线成人网| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 中文亚洲av片在线观看爽| x7x7x7水蜜桃| 香蕉av资源在线| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩卡通动漫| 身体一侧抽搐| 校园春色视频在线观看| 女生性感内裤真人,穿戴方法视频| 婷婷亚洲欧美| 久久久国产成人免费| 搡老岳熟女国产| 高潮久久久久久久久久久不卡| 精品国产乱码久久久久久男人| 午夜两性在线视频| 国产成人一区二区三区免费视频网站| 亚洲片人在线观看| 99国产精品99久久久久| 日日干狠狠操夜夜爽| 精品一区二区三区视频在线观看免费| 首页视频小说图片口味搜索| 亚洲精品色激情综合| 亚洲自偷自拍图片 自拍| 香蕉国产在线看| 国产精品久久久av美女十八| 夜夜躁狠狠躁天天躁| 搡老妇女老女人老熟妇| 午夜视频精品福利| 欧洲精品卡2卡3卡4卡5卡区| 九九热线精品视视频播放| 国产亚洲av高清不卡| 国产精品亚洲av一区麻豆| 97超视频在线观看视频| 欧美日本视频| 嫁个100分男人电影在线观看| 亚洲欧美日韩无卡精品| 中文字幕久久专区| 免费观看人在逋| 夜夜看夜夜爽夜夜摸| 国产成人欧美在线观看| 久久精品国产亚洲av香蕉五月| xxx96com| 亚洲av第一区精品v没综合| 又紧又爽又黄一区二区| 国产精品亚洲av一区麻豆| 天堂√8在线中文| 两性午夜刺激爽爽歪歪视频在线观看| 99在线视频只有这里精品首页| 国产高清视频在线播放一区| 国产成人啪精品午夜网站| 舔av片在线| 精品国产超薄肉色丝袜足j| 首页视频小说图片口味搜索| 国产av不卡久久| 网址你懂的国产日韩在线| 在线观看免费午夜福利视频| 草草在线视频免费看| 在线观看舔阴道视频| 欧美成狂野欧美在线观看| 国产一区二区在线观看日韩 | 欧美不卡视频在线免费观看| bbb黄色大片| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 国产淫片久久久久久久久 | 757午夜福利合集在线观看| 国产真实乱freesex| 国产主播在线观看一区二区| 欧美日本视频| 麻豆成人午夜福利视频| 在线视频色国产色| 日日夜夜操网爽| 欧美成人一区二区免费高清观看 | 久久久久精品国产欧美久久久| 亚洲九九香蕉| 99久久久亚洲精品蜜臀av| 丰满人妻熟妇乱又伦精品不卡| 99国产综合亚洲精品| 免费在线观看影片大全网站|