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

    WRF模式對(duì)浙江2011年夏季降水和溫度預(yù)報(bào)評(píng)估及其濕過(guò)程敏感性分析*

    2012-12-07 06:51:16董美瑩冀春曉林惠娟
    浙江氣象 2012年3期
    關(guān)鍵詞:積云實(shí)況降水

    陳 鋒 董美瑩 冀春曉 林惠娟 查 賁 駱 陽(yáng)

    (1.浙江省氣象科學(xué)研究所,浙江 杭州310017;2.蘇州市氣象局,江蘇 蘇州215021;3.杭州市氣象局,浙江杭州310051;4.浙江省氣象信息網(wǎng)絡(luò)中心,浙江杭州310017)

    0 引言

    隨著高性能計(jì)算能力和天氣預(yù)報(bào)技術(shù)的逐步發(fā)展,數(shù)值天氣預(yù)報(bào)已成為現(xiàn)代天氣預(yù)報(bào)業(yè)務(wù)的重要基礎(chǔ)和天氣預(yù)報(bào)業(yè)務(wù)發(fā)展的主流方向。由于受初始場(chǎng)、邊界場(chǎng)、物理過(guò)程及模式本身設(shè)計(jì)等諸多因素影響,模式預(yù)報(bào)產(chǎn)品不可避免存在誤差。因此,檢驗(yàn)評(píng)估工作已成為數(shù)值天氣預(yù)報(bào)的一個(gè)重要組成部分,評(píng)估結(jié)果有助于模式本身的改進(jìn)和預(yù)報(bào)人員更好地使用模式產(chǎn)品。

    王雨(2006)利用國(guó)家基本氣象站資料評(píng)估了中央臺(tái)5個(gè)常用模式的定量降水預(yù)報(bào)結(jié)果;王麗等(2007)利用湖北省77個(gè)氣象觀測(cè)站資料對(duì)日本、德國(guó)和T213 3種數(shù)值模式在強(qiáng)降水過(guò)程中的預(yù)報(bào)進(jìn)行了統(tǒng)計(jì)比較;藥明和晏曉英(2008)利用代表不同區(qū)域地理環(huán)境的6個(gè)站點(diǎn),檢驗(yàn)并分析了吉林省數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)短期預(yù)報(bào)各時(shí)段降水、溫度、氣壓和風(fēng)速等要素預(yù)報(bào)性能;段旭等(2011)在建立了本地化的中尺度WRF數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)后,也對(duì)系統(tǒng)預(yù)報(bào)的降水和氣溫做了評(píng)估??傮w而言,已有的數(shù)值預(yù)報(bào)模式檢驗(yàn)評(píng)估主要側(cè)重于對(duì)定量降水的評(píng)估,且評(píng)估所用觀測(cè)資料多以常規(guī)站資料為主,利用中尺度自動(dòng)站詳細(xì)資料針對(duì)浙江省數(shù)值預(yù)報(bào)的評(píng)估工作尚未見(jiàn)報(bào)告。浙江地處我國(guó)東南沿海,地理地貌獨(dú)特、天氣復(fù)雜,究竟數(shù)值預(yù)報(bào)模式對(duì)浙江省夏季降水和溫度有多大預(yù)報(bào)能力,存在多大的誤差,在何種程度上能為預(yù)報(bào)員提供可信的參考信息,已成為模式開(kāi)發(fā)人員和預(yù)報(bào)員共同關(guān)注的重要問(wèn)題,使用更詳盡觀測(cè)資料評(píng)價(jià)高分辨模式預(yù)報(bào)性能的精細(xì)化評(píng)估技術(shù)值得探索。

    另外,國(guó)內(nèi)外大量研究表明,模式本身設(shè)計(jì)及其參數(shù)化方案對(duì)模式的預(yù)報(bào)性能有著重要影響(黃安寧等,2009;肖玉華等,2010;張麗霞等,2011)。模式的濕過(guò)程選擇對(duì)于濕熱夏季的預(yù)報(bào)效果十分重要(許習(xí)華和丁一匯,1990;Wang and Seaman,1997;劉一鳴,1998;樓小鳳等,2003;伍華平等,2009;Rajeevan et al.,2010)。模式的不同濕過(guò)程方案選擇會(huì)對(duì)預(yù)報(bào)結(jié)果有多大影響?究竟哪些濕過(guò)程的微物理方案和積云參數(shù)化方案更適合浙江的天氣特點(diǎn)?這對(duì)更好地開(kāi)展WRF模式在浙江省的本地化適用研究,提高模式預(yù)報(bào)性能,具有重要意義。

    因此,本文首先利用全省1500多個(gè)包括區(qū)域自動(dòng)站在內(nèi)的站點(diǎn)觀測(cè)資料,應(yīng)用空間分布和時(shí)間演變及統(tǒng)計(jì)檢驗(yàn)3種方法來(lái)評(píng)價(jià)WRF模式對(duì)浙江省2011年夏季降水和溫度的預(yù)報(bào)性能;進(jìn)而通過(guò)不同濕過(guò)程參數(shù)化方案數(shù)值模擬的對(duì)比,探討不同微物理參數(shù)化方案和積云參數(shù)化方案對(duì)模式預(yù)報(bào)性能的影響;以期為高分辨率數(shù)值模式的改進(jìn)和系統(tǒng)用戶做精細(xì)化天氣預(yù)報(bào)提供較為客觀全面的參考信息。

    1 模式簡(jiǎn)介

    浙江省中尺度數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)采用WRF v3.1 模式(Weather Research and Forecasting model version 3.1)。該模式(Skamarock et al.,2008)是由美國(guó)NCAR等多家研究單位聯(lián)合研制的新一代數(shù)值模式,是一個(gè)完全可壓縮、非靜力模式,控制方程組為通量形式,水平網(wǎng)格采用有利于在高分辨率模擬中提高數(shù)值計(jì)算精度的Arakawa C 格式(Kalnay,2001,Janjic,2004)。WRF模式適用于業(yè)務(wù)天氣預(yù)報(bào)和大氣研究等需要,能有效模擬天氣系統(tǒng)尤其是中小天氣系統(tǒng)的發(fā)生、發(fā)展、消亡過(guò)程及其機(jī)理。同時(shí),模式系統(tǒng)具有可移植、易維護(hù)、可擴(kuò)充、高效率、使用方便等諸多特性,被大量應(yīng)用于理想試驗(yàn)、個(gè)例數(shù)值試驗(yàn)等研究和數(shù)值天氣預(yù)報(bào)實(shí)時(shí)業(yè)務(wù)之中(Klemp,2004;Dimego,2004;Isidora et al.,2005;Lo et al.,2008;范水勇等,2008;馬紅云等,2009;吳偉等,2011;林惠娟等,2010),表現(xiàn)出較好的模擬能力和較可靠的預(yù)報(bào)性能。

    根據(jù)浙江省現(xiàn)有計(jì)算條件和地理位置特征,模式選用雙重嵌套,分辨率各為9 km和3 km,垂直層次為35層。9 km分辨率的大區(qū)域(華東區(qū)域,大致范圍為 106°E ~128°E,22°N ~42°N)重點(diǎn)覆蓋影響浙江的大尺度系統(tǒng),3 km分辨率的小區(qū)域(浙江區(qū)域,大致范圍為117°E~124°E,26°N ~32°N)重在識(shí)別影響浙江的中尺度系統(tǒng)(圖1)。模式積云參數(shù)化方案在大區(qū)域中采用Betts-Miller-Janjic方案,在小區(qū)域不采用參數(shù)化方案;微物理方案均采用WSM 6-Class方案;陸面過(guò)程均使用Noah方案;行星邊界層均采用Yonsei University(YSU)參數(shù)化方案;表面層均使用基于Monin-Obukhov的MM5相似理論;長(zhǎng)波、短波輻射均選用RRTM快速輻射傳輸方案和 Dudhia方案。初始場(chǎng)由 NCEP/GFS(Globe Forecast System)的分析場(chǎng)提供,側(cè)邊界場(chǎng)由NCEP/GFS間隔3 h的預(yù)報(bào)場(chǎng)提供。

    圖1 模式研究區(qū)域及地形

    2 2011年夏季降水和溫度預(yù)報(bào)結(jié)果評(píng)估

    浙江地屬亞熱帶季風(fēng)氣候,冬季干冷,夏季濕熱,降水量的季節(jié)分配大體上與溫度一致,表現(xiàn)為“雨熱同步”的氣候特征(祝啟桓,1992)。其中,夏季的主要影響系統(tǒng)為梅雨、臺(tái)風(fēng)和強(qiáng)對(duì)流等高影響、災(zāi)害性天氣,對(duì)這些天氣的預(yù)報(bào)能力直接關(guān)系到防災(zāi)減災(zāi)工作的成敗。因此夏季降水和溫度預(yù)報(bào)水平是衡量數(shù)值預(yù)報(bào)模式性能的一個(gè)重要指標(biāo)。本節(jié)對(duì)2011年夏季(6—8月)模式預(yù)報(bào)的降水和溫度結(jié)果進(jìn)行定量評(píng)估。

    2.1 資料和方法

    實(shí)況資料選用浙江省內(nèi)1500多個(gè)站點(diǎn)觀測(cè)資料(包括3個(gè)基準(zhǔn)站、20個(gè)基本站、47個(gè)一般站和1500多個(gè)區(qū)域自動(dòng)站,如圖2所示),評(píng)估要素為每日08:00—次日08:00的24 h累計(jì)降水和平均2 m氣溫。預(yù)報(bào)結(jié)果為與實(shí)況相對(duì)應(yīng)的每天08:00起報(bào)的預(yù)報(bào)時(shí)效為24 h,48 h及72 h的24 h累計(jì)降水和平均2 m氣溫。

    圖2 觀測(cè)站點(diǎn)分布(Level1——國(guó)家基準(zhǔn)氣候站;Level2——國(guó)家基本氣象站;Level3——國(guó)家一般氣象站;Level4——區(qū)域自動(dòng)站)

    為更好地評(píng)估模式結(jié)果,文中給出了降水和溫度的空間分布和時(shí)間序列演變檢驗(yàn)。對(duì)于降水,同時(shí)給出了TS,ETS(Equitable Threat Score),HSS(Heidke Skill Score),擊中率(PODY),空?qǐng)?bào)率(FAR),漏報(bào)率(FOM)等多個(gè)檢驗(yàn)量(丁金才,1995)共同描述模式的降水預(yù)報(bào)性能。

    2.1 降 水

    2011年夏季降水過(guò)程比較明顯,6月份梅雨期間出現(xiàn)了4次典型的梅雨降水過(guò)程,7—8月主要是以臺(tái)風(fēng)降水(“梅花”和“南瑪都”)和局地強(qiáng)對(duì)流降水為主。由圖3(a)夏季平均降水實(shí)況可知,浙江省2011年夏季降水主要集中于浙西北部山區(qū)、寧波地區(qū)及浙東南沿海地區(qū)。圖3(b-d)分別給出了模式預(yù)報(bào)的預(yù)報(bào)時(shí)效為24 h,48 h及72 h的夏季日平均降水分布。與實(shí)況對(duì)比可知,模式24 h預(yù)報(bào)報(bào)出了浙西北部和東南沿海地區(qū)的降水中心及雨帶的基本走向,預(yù)報(bào)和實(shí)況的空間相關(guān)系數(shù)為0.52;但預(yù)報(bào)降水量比實(shí)況整體偏大,平均偏大為5.62 mm/d。48 h預(yù)報(bào)和72 h預(yù)報(bào)也能在一定程度上預(yù)報(bào)出降水帶的位置,空間相關(guān)系數(shù)分別為0.36和0.42;降水量級(jí)相對(duì)于24 h預(yù)報(bào)結(jié)果更接近于實(shí)況,平均偏差分別為3.85和2.17 mm/d。這說(shuō)明模式在48 h和72 h的預(yù)報(bào)結(jié)果仍然是有一定的參考價(jià)值的。但是,在浙西南地區(qū),不論24 h,48 h還是72 h均存在空?qǐng)?bào)現(xiàn)象,數(shù)值預(yù)報(bào)產(chǎn)品的參考性存在地域差異,需要謹(jǐn)慎分析。

    圖3 2011年夏季日平均累計(jì)降水預(yù)報(bào)及實(shí)況(單位:mm/d):(a)實(shí)況;(b)24 h預(yù)報(bào);(c)48 h預(yù)報(bào);(d)72 h預(yù)報(bào)

    夏季的日平均降水分布反映了模式對(duì)季節(jié)內(nèi)降水預(yù)報(bào)空間分布情況的預(yù)報(bào)能力,而區(qū)域平均逐日降水演變則能較好反映模式對(duì)降水過(guò)程隨時(shí)間發(fā)展趨勢(shì)的預(yù)報(bào)能力。浙江省2011年夏季逐日降水的時(shí)間演變(圖4)表明,模式各預(yù)報(bào)時(shí)效的預(yù)報(bào)均給出了較為準(zhǔn)確的降水發(fā)生、發(fā)展及減弱過(guò)程,尤其是在6月和8月的幾次較大降水過(guò)程,模式預(yù)報(bào)的降水量級(jí)與實(shí)況相當(dāng)。但是,在降水強(qiáng)度為小到中雨時(shí)(如7月期間,31~61 d),模式預(yù)報(bào)降水偏大較明顯,尤其是24 h預(yù)報(bào)結(jié)果。這可能跟模式參數(shù)化方案的閾值偏低有關(guān),需要進(jìn)一步研究改進(jìn)。

    圖4 2011年夏季逐日降水量預(yù)報(bào)及實(shí)況

    利用 TS,ETS,HSS,擊中率、空?qǐng)?bào)率、漏報(bào)率等參數(shù),定量降水累加檢驗(yàn)結(jié)果列于表1。由表1分析可知,模式對(duì)降水的預(yù)報(bào)能力隨著降水量級(jí)的增大而減小,24 h預(yù)報(bào)的TS評(píng)分從小雨以上級(jí)別的0.59逐步下降到大暴雨以上級(jí)別的0.02;模式預(yù)報(bào)結(jié)果較實(shí)況偏大,空?qǐng)?bào)多,漏報(bào)少,空?qǐng)?bào)率高出漏報(bào)率約0.21,尤其是24 h預(yù)報(bào)結(jié)果偏大較多。比較24 h,48 h及72 h預(yù)報(bào)的TS評(píng)分來(lái)看,模式24 h預(yù)報(bào)結(jié)果較好,且48 h和72 h預(yù)報(bào)結(jié)果雖略有下降但仍存在一定的應(yīng)用價(jià)值。

    2.2 溫 度

    溫度預(yù)報(bào)是日常天氣預(yù)報(bào)業(yè)務(wù)中的另一個(gè)重要內(nèi)容。浙江省2011年夏季日平均2 m氣溫預(yù)報(bào)和實(shí)況顯示(圖5):2 m氣溫的空間分布受海拔和緯度影響較大,高分辨率模式能較好地刻畫(huà)浙江省的地形分布,對(duì)浙北地區(qū)、金衢盆地及浙東南沿海平原地區(qū)的高溫區(qū)有很好的預(yù)報(bào)能力,24 h預(yù)報(bào)與實(shí)況的空間相關(guān)系數(shù)達(dá)0.77。尤其是模式對(duì)金衢盆地的高值區(qū)的預(yù)報(bào),不僅在空間分布上與實(shí)況很相似,而且在量級(jí)上也較接近。但是,模式對(duì)其它地區(qū)還存在著1~2℃左右的低估,24 h預(yù)報(bào)的平均偏差為-1.12℃。值得一提的是,模式48 h和72 h預(yù)報(bào)結(jié)果與24 h預(yù)報(bào)結(jié)果相差很小(偏差分別為-1.21℃和 -1.20℃),這說(shuō)明模式對(duì) 2 m 平均氣溫有較長(zhǎng)的預(yù)報(bào)時(shí)效。

    表1 2011年夏季浙江省降水預(yù)報(bào)累加檢驗(yàn)結(jié)果

    圖5 2011年夏季日平均2 m氣溫預(yù)報(bào)及實(shí)況(單位:℃):(a)實(shí)況;(b)24 h預(yù)報(bào);(c)48 h預(yù)報(bào);(d)72 h預(yù)報(bào)

    如圖6所示,浙江省2011年夏季逐日2 m平均氣溫的時(shí)間演變表明,模式能較為準(zhǔn)確的預(yù)報(bào)出氣溫的逐日變化過(guò)程。對(duì)比圖4和圖6發(fā)現(xiàn),氣溫變化與降水變化有較好的反相關(guān)關(guān)系,模式在6月和8月降水較多時(shí)期氣溫較低,而在6月末至8月初降水較少時(shí)期氣溫較高,這種聯(lián)系在物理上是合理的。仔細(xì)對(duì)比還可看到,模式對(duì)降水較大的6月中上旬和8月下旬預(yù)報(bào)的氣溫與實(shí)況較為一致,但在降水量級(jí)為小到中雨的7月,模式預(yù)報(bào)氣溫較實(shí)況偏低1~2℃,且與降水預(yù)報(bào)偏差在時(shí)間上對(duì)應(yīng)較為一致。初步分析認(rèn)為,這可能是由于在弱降水期間,實(shí)況云量較少,而模式預(yù)報(bào)的云量偏多,預(yù)報(bào)誤差較大;但在較大降水期間,實(shí)況云量多,覆蓋了全省大部,此時(shí)預(yù)報(bào)云量的誤差小于弱降水期間的誤差(即弱降水期間云量預(yù)報(bào)偏多誤差大于強(qiáng)降水);而云量預(yù)報(bào)偏多將導(dǎo)致到達(dá)地面的凈輻射偏小,溫度偏低,因此會(huì)出現(xiàn)在弱降水期間溫度預(yù)報(bào)偏低誤差大于強(qiáng)降水期間的現(xiàn)象。

    圖6 2011年夏季逐日2m平均氣溫預(yù)報(bào)及實(shí)況

    3 濕過(guò)程參數(shù)化方案敏感性分析

    前文的評(píng)估給出了對(duì)2011年夏季浙江省中尺度數(shù)值預(yù)報(bào)模式預(yù)報(bào)水平的整體認(rèn)識(shí)。研究表明,模式的濕過(guò)程選擇對(duì)于模式降水的預(yù)報(bào)效果十分重要,不同濕過(guò)程的微物理方案和積云參數(shù)化方案會(huì)對(duì)模式預(yù)報(bào)性能產(chǎn)生影響。因此,本節(jié)選取2011年6月4—7日的浙江省典型梅雨過(guò)程,著重探討不同濕過(guò)程的微物理方案和積云參數(shù)化方案的組合對(duì)模式預(yù)報(bào)性能的影響和由此帶來(lái)的不確定性,進(jìn)而加深對(duì)于數(shù)值模式預(yù)報(bào)能力的認(rèn)識(shí)。敏感性試驗(yàn)設(shè)計(jì)見(jiàn)表2,即選取了適用于高分辨中尺度模式的微物理方案和積云參數(shù)化方案各3種,設(shè)計(jì)了9個(gè)敏感性試驗(yàn),對(duì)6月4—7日的梅雨過(guò)程的降水預(yù)報(bào)結(jié)果進(jìn)行對(duì)比分析。

    表2 敏感性試驗(yàn)設(shè)計(jì)

    圖7 各濕過(guò)程組合方案下模式對(duì)2011年6月4日08:00—5日08:00降水預(yù)報(bào)及實(shí)況:(a)S1-WSM6+KF方案;(b)S2-WSM6+BMJ方案;(c)S3-WSM6+Grell方案;(d)S4-WSM3+KF方案;(e)S5-WSM3+BMJ方案;(f)S6-WSM3+Grell方案;(g)S7-Lin+KF方案;(h)S8-Lin+BMJ方案;(i)S9-WSM6+Grell方案;(j)實(shí)況。

    圖7 是各濕過(guò)程組合方案下模式對(duì)2011年6月4日08:00—5日08:00降水預(yù)報(bào)和實(shí)況對(duì)比。由圖可知,不同濕過(guò)程參數(shù)化方案對(duì)降水預(yù)報(bào)結(jié)果差異較大,各類(lèi)方案對(duì)降水預(yù)報(bào)的偏差為4.1~25.6 mm/d不等,空間相關(guān)系數(shù)為-0.11~0.59不等。相比而言,降水對(duì)于積云參數(shù)化方案的敏感性明顯大于對(duì)微物理方案的敏感性:積云參數(shù)化方案的不同可以導(dǎo)致整體雨帶分布的明顯不同,而微物理方案不同對(duì)于雨帶影響不大,它主要影響雨帶內(nèi)強(qiáng)降水中心的位置和強(qiáng)度。KF方案和Grell方案在浙江北部及浙中地區(qū)報(bào)了兩個(gè)雨帶,尤其是在浙中地區(qū)出現(xiàn)了嚴(yán)重的空?qǐng)?bào),BMJ方案所報(bào)雨帶集中在浙江北部,空?qǐng)?bào)相對(duì)較少。在選定BMJ方案下,微物理方案對(duì)雨帶位置和降水中心強(qiáng)度進(jìn)行了微調(diào)。從總體上看 BMJ方案配合WSM3,WSM6和Lin方案所預(yù)報(bào)的結(jié)果均與實(shí)況較為接近,從誤差角度分析,BMJ方案配合WSM3方案最小,WSM6方案次之,但是比較兩者預(yù)報(bào)和實(shí)況位于浙北中部的東西向的強(qiáng)雨帶可見(jiàn),BMJ方案配合WSM6方案所預(yù)報(bào)的降水雨帶和中心明顯更接近實(shí)況,它的相關(guān)系數(shù)也是最高的。圖8給出了各方案組合下模式降水預(yù)報(bào)的評(píng)分結(jié)果。對(duì)比各個(gè)方案來(lái)看,模式降水的預(yù)報(bào)能力隨著降水量級(jí)的增加而下降,且在中小雨級(jí)別得分均較高,TS評(píng)分均在0.7以上;各方案對(duì)大雨(25 mm)和暴雨(50 mm)級(jí)別的預(yù)報(bào)性能差異較為明顯,方案2(BMJ與WSM6組合方案)相對(duì)其他方案有較好表現(xiàn),其TS評(píng)分、ETS評(píng)分、HSS評(píng)分和擊準(zhǔn)率(PODY)均高于其他方案,而空?qǐng)?bào)率(FAR)和漏報(bào)率(FOM)均低于其他方案??傮w上,采用方案2(BMJ積云對(duì)流參數(shù)化方案和WSM6微物理方案)預(yù)報(bào)的降水落區(qū)和強(qiáng)度較其他方案更接近于實(shí)況。

    圖8 各濕過(guò)程組合方案下模式對(duì)2011年6月4日08:00—5日08:00降水預(yù)報(bào)評(píng)分結(jié)果

    為分析不同參數(shù)化方案對(duì)降水形成原因,圖9給出了各方案下模式預(yù)報(bào)的格點(diǎn)尺度降水與積云對(duì)流降水的對(duì)比。由圖可見(jiàn),格點(diǎn)尺度降水與積云對(duì)流降水的比值受積云參數(shù)化方案影響較大,BMJ方案預(yù)報(bào)的降水以格點(diǎn)尺度降水為主,積云對(duì)流降水僅占總降水量的5%左右;KF方案預(yù)報(bào)的格點(diǎn)尺度降水比積云對(duì)流降水略多30%左右;而Grell方案預(yù)報(bào)的格點(diǎn)尺度降水和積云對(duì)流降水相當(dāng)。

    圖9 各濕過(guò)程組合方案下模式預(yù)報(bào)的格點(diǎn)尺度降水與積云對(duì)流降水對(duì)比

    分析實(shí)況表明,本次降水過(guò)程主要與西南暖濕氣流和北方弱冷空氣在浙北地區(qū)形成的緯向切變線有關(guān)。對(duì)比各方案組合下模式積分24 h后850 hPa風(fēng)場(chǎng)可以看到(圖10),方案2預(yù)報(bào)出了(BMJ與WSM6組合方案)在30°N附近的切變,較其余方案更接近實(shí)況;同時(shí)方案2預(yù)報(bào)的850 hPa濕度場(chǎng)在浙北和浙西較高、浙東南較低的分布形態(tài)也與實(shí)況較符合。這可能是該方案雨帶位置相對(duì)準(zhǔn)確的原因之一。由此可見(jiàn),濕過(guò)程參數(shù)化方案的不同可以導(dǎo)致風(fēng)場(chǎng)產(chǎn)生差異,進(jìn)而影響到切變線和雨帶的分布,從而對(duì)降水的預(yù)報(bào)產(chǎn)生重要影響。

    圖10 各濕過(guò)程組合方案下模式積分24 h(5日08時(shí))500 hPa位勢(shì)高度場(chǎng)及850 hPa風(fēng)場(chǎng)和濕度場(chǎng):(a)S1-WSM6+KF方案;(b)S2-WSM6+BMJ方案;(c)S3-WSM6+Grell方案;(d)S4-WSM3+KF方案;(e)S5-WSM3+BMJ方案;(f)S6-WSM3+Grell方案;(g)S7-Lin+KF方案;(h)S8-Lin+BMJ方案;(i)S9-WSM6+Grell方案;(j)實(shí)況。

    4 結(jié)語(yǔ)

    綜上所述,本文利用全省1500多個(gè)包括區(qū)域自動(dòng)站在內(nèi)的站點(diǎn)觀測(cè)資料,從空間分布和時(shí)間演變及統(tǒng)計(jì)檢驗(yàn)3個(gè)方面對(duì)浙江省2011年夏季W(wǎng)RF模式預(yù)報(bào)的降水和溫度進(jìn)行了定量評(píng)估;并通過(guò)不同濕過(guò)程參數(shù)化方案數(shù)值試驗(yàn)的對(duì)比,探討了不同微物理參數(shù)化方案和積云參數(shù)化方案對(duì)模式預(yù)報(bào)性能的影響,得到以下初步認(rèn)識(shí)。

    (1)WRF模式能基本預(yù)報(bào)出降水的空間分布及時(shí)間演變規(guī)律,對(duì)降水落區(qū)預(yù)報(bào)和降水發(fā)生演變有較好的指示意義,但模式預(yù)報(bào)范圍和量值均偏大,其中在小到中雨量級(jí)期間降水量值偏大較明顯。

    (2)模式能很好地預(yù)報(bào)出2 m平均氣溫的空間分布和時(shí)間演變過(guò)程,但存在著一定的低估現(xiàn)象,尤其是在浙江省弱降水期間,模式預(yù)報(bào)偏低1~2℃。

    (3)就區(qū)域平均而言,模式在實(shí)況降水量較大期間預(yù)報(bào)性能較好,但在小到中雨期間較差,主要表現(xiàn)為降水量的高估和氣溫的低估。這可能是由于模式在小到中雨期間對(duì)云量預(yù)報(bào)偏多引起,還需要進(jìn)一步的研究。

    (4)積云參數(shù)化方案對(duì)于降水影響明顯,它可以導(dǎo)致整體雨帶分布的明顯不同;而微物理方案不同對(duì)于雨帶影響不大,它主要影響雨帶內(nèi)強(qiáng)降水中心的位置和強(qiáng)度;總體上,采用Betts-Miller-Janjic積云對(duì)流參數(shù)化方案和WSM6微物理方案預(yù)報(bào)的降水落區(qū)和強(qiáng)度更接近于實(shí)況。

    (5)不同參數(shù)化方案組合可以引起預(yù)報(bào)流場(chǎng)的差異,從而影響到雨帶位置,對(duì)降水的預(yù)報(bào)性能產(chǎn)生不可忽視的影響。

    總之,浙江省高分辨率中尺度數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)對(duì)降水和溫度的空間分布和時(shí)間演變規(guī)律有較好的預(yù)報(bào)能力,對(duì)預(yù)報(bào)員具有一定的參考價(jià)值。根據(jù)各量級(jí)降水TS,ETS等各項(xiàng)評(píng)估指標(biāo)顯示,模式的預(yù)報(bào)水平和國(guó)家氣象中心常用數(shù)值模式的預(yù)報(bào)能力相當(dāng)(王雨,2006)?;谡憬?500多個(gè)觀測(cè)站的詳細(xì)觀測(cè)資料,對(duì)于水平格距為3 km×3 km的高分辨數(shù)值模式預(yù)報(bào)結(jié)果進(jìn)行檢驗(yàn),以上評(píng)估結(jié)果是令人鼓舞的。誠(chéng)然,模式系統(tǒng)還存在降水預(yù)報(bào)偏大等問(wèn)題,還需要逐步改進(jìn),尤其是對(duì)已有評(píng)估結(jié)果,還需要進(jìn)一步探究原因。

    [1] Dimego G.WRF Development Activities at NCEP[C].The 84thAMS Annual Meeting.Seattle,U S A,Jan 10 -15,2004.

    [2] Isidora J,A G JWilliam,S.Moti,et al.The impact of different WRF model physical parameterizations and their interactions on warm season MCSrainfall[J].Wea.Forecasting,2005,20,1048 -1060.

    [3] Janjic Z I.The NCEPWRF Core[C].The84thAMSAnnual Meeting.Seattle,U SA,Jan 10 -15,2004.

    [4] Kalnay E.Numerical weather forecasting and predictability[M].Scotland:Cambridge University Press,2001:341-343.

    [5] Klemp JB.Weather research and forecasting model:a technical overview[C].The 84thAMSAnnual Meeting.Seattle,U SA,Jan 10-15,2004:32-36.

    [6] Lo JCF,ZL Yang,R A Pielke Sr.Assessment of three dynamical climate downscaling methods using the Weather Research and Forecasting(WRF)model[J].J Geophys Res,2008,113,D09112,doi:10.1029/JD009216.

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

    [8] Skamarock W C,JB Klemp,JDudhia,et al.A description of the Advanced Research WRF version 3.NCAR Tech Notes-475+STR,2008.

    [9] Wang Wand N L Seaman.A comparison study of convective parameterization schemes in a mesoscale model[J].Mon Wea Rev,1997,125:252 - 278,doi:http://dx.doi.org/10.1175/1520 -0493(1997)125 <0252:ACSOCP >2.0.CO;2.

    [10] 丁金才.天氣預(yù)報(bào)評(píng)分方法評(píng)述[J].南京氣象學(xué)院學(xué)報(bào),1995,18(1):143 -150.

    [11] 段旭,王曼,陳新梅,等.中尺度WRF數(shù)值模式系統(tǒng)本地化業(yè)務(wù)試驗(yàn)[J].氣象,2011,37(1):39 -47.

    [12] 范水勇,郭永潤(rùn),陳敏,等.高分辨率WRF三維變分同化在北京地區(qū)降水預(yù)報(bào)中的應(yīng)用[J].高原氣象,2008,27(6):1181-1188.

    [13] 黃安寧,張耀存,朱堅(jiān).中國(guó)夏季不同強(qiáng)度降水模擬對(duì)不同積云對(duì)流參數(shù)化方案的敏感性研究[J].大氣科學(xué),2009,33(6):1212 -1224.

    [14] 林惠娟,冀春曉.WRF參數(shù)化方案對(duì)臺(tái)風(fēng)路徑和強(qiáng)度模擬的影響[C].第七屆長(zhǎng)三角氣象科技論壇文集,2010.

    [15] 劉一鳴.中國(guó)關(guān)于積云參數(shù)化方案的應(yīng)用[J].氣象學(xué)報(bào),1998,56(2):247 -255.

    [16] 樓小鳳,胡志晉,王鵬云,等.中尺度模式云降水物理方案介紹[J].應(yīng)用氣象學(xué)報(bào),2003,14(增刊):49 -59.

    [17] 馬紅云,郭品文,宋潔.耦合不同陸面方案的WRF模式對(duì)2007年7月江淮強(qiáng)降水過(guò)程的模擬[J].大氣科學(xué),2009,33(3):557 -567.

    [18] 王麗,金琪,柯怡明.三種數(shù)值預(yù)報(bào)產(chǎn)品短期強(qiáng)降水預(yù)報(bào)定量誤差評(píng)估[J].暴雨災(zāi)害,2007,26(4):301 -305.

    [19] 王雨.2004年主汛期各數(shù)值預(yù)報(bào)模式定量降水預(yù)報(bào)評(píng)估[J].應(yīng)用氣象學(xué)報(bào),2006,17(3):316 -324.

    [20] 吳偉,鄧蓮堂,王式功.“0911”華北暴雪的數(shù)值模擬及云微物理特征分析[J].氣象,2011,37(8):991 -998.

    [21] 伍華平,束炯,顧瑩,等.暴雨模擬中積云對(duì)流參數(shù)化方案的對(duì)比試驗(yàn)[J].熱帶氣象學(xué)報(bào),2009,25(2):175 -180.

    [22] 肖玉華,何光碧,顧清源,等.邊界層參數(shù)化方案對(duì)不同性質(zhì)降水模擬的影響[J].高原氣象,2010,29(2):331-339.

    [23] 許習(xí)華,丁一匯.一種考慮云相互作用的積云對(duì)流參數(shù)化方案及其在降水預(yù)報(bào)中的應(yīng)用[J].中國(guó)科學(xué)(B輯),1990,9:998 -1008.

    [24] 藥明,晏曉英.吉林省數(shù)值預(yù)報(bào)業(yè)務(wù)系統(tǒng)要素預(yù)報(bào)性能分析[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2008,38(增刊):209-213.

    [25] 張麗霞,周天軍,曾先鋒,等.積云參數(shù)化方案對(duì)熱帶降水年循環(huán)模態(tài)模擬的影響[J].大氣科學(xué),2011,35(4):777-790.

    [26] 祝啟桓,張淑云,顧強(qiáng)民,等.浙江省災(zāi)害性天氣預(yù)報(bào)[M].北京:氣象出版社,1992:1-5.

    猜你喜歡
    積云實(shí)況降水
    黑龍江省玉米生長(zhǎng)季自然降水與有效降水對(duì)比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    WRF模式中積云對(duì)流參數(shù)化方案對(duì)南海土臺(tái)風(fēng)“Ryan”模擬的影響研究
    鄉(xiāng)村小學(xué)的愿望與現(xiàn)實(shí)——宜君鄉(xiāng)村教育實(shí)況
    積云——晴天的使者
    Warning signs of an approaching hurricane
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    天舟一號(hào)貨運(yùn)飛船發(fā)射實(shí)況掠影
    可愛(ài)潮咖們的獨(dú)門(mén)彩妝實(shí)況直播
    降水現(xiàn)象儀模擬軟件設(shè)計(jì)與實(shí)現(xiàn)
    牽掛
    欧美一区二区精品小视频在线| 亚洲人成伊人成综合网2020| 色播在线永久视频| 久热爱精品视频在线9| 99久久人妻综合| av视频免费观看在线观看| 亚洲国产精品一区二区三区在线| 国产97色在线日韩免费| 91九色精品人成在线观看| 18禁观看日本| 老汉色av国产亚洲站长工具| 国产色视频综合| 欧美一级毛片孕妇| 免费日韩欧美在线观看| 精品国产一区二区三区四区第35| 夜夜爽天天搞| 欧美日韩亚洲国产一区二区在线观看| 老司机亚洲免费影院| 999久久久国产精品视频| 亚洲av美国av| 69精品国产乱码久久久| 美女福利国产在线| 久久久久久人人人人人| 琪琪午夜伦伦电影理论片6080| 国产成人影院久久av| 成年女人毛片免费观看观看9| 91精品国产国语对白视频| 国产一区二区三区视频了| 一进一出好大好爽视频| 青草久久国产| 脱女人内裤的视频| 亚洲全国av大片| 一二三四社区在线视频社区8| 久久精品成人免费网站| 人人澡人人妻人| 在线观看免费视频日本深夜| 91麻豆精品激情在线观看国产 | 韩国av一区二区三区四区| 精品电影一区二区在线| 露出奶头的视频| 国产欧美日韩综合在线一区二区| 啦啦啦在线免费观看视频4| 不卡一级毛片| 五月开心婷婷网| 精品熟女少妇八av免费久了| 精品国产国语对白av| 嫩草影视91久久| 欧美激情极品国产一区二区三区| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| 成人三级黄色视频| a级毛片在线看网站| 精品久久久久久成人av| 久久久久国产精品人妻aⅴ院| 国产在线观看jvid| 久久久久久久精品吃奶| 欧美另类亚洲清纯唯美| 久久精品91蜜桃| 国产精品免费一区二区三区在线| 中文字幕人妻熟女乱码| 精品久久久精品久久久| 亚洲美女黄片视频| 欧美成人性av电影在线观看| 午夜成年电影在线免费观看| 亚洲欧美精品综合一区二区三区| 91麻豆精品激情在线观看国产 | 好男人电影高清在线观看| 国产一区二区在线av高清观看| 午夜精品久久久久久毛片777| xxxhd国产人妻xxx| 国产男靠女视频免费网站| 黑丝袜美女国产一区| 97超级碰碰碰精品色视频在线观看| 在线国产一区二区在线| 亚洲第一av免费看| 母亲3免费完整高清在线观看| www.精华液| 亚洲专区国产一区二区| 国产精品二区激情视频| 国产成人av教育| 精品高清国产在线一区| 丁香欧美五月| 日日干狠狠操夜夜爽| 91麻豆av在线| 日韩国内少妇激情av| 欧美日韩亚洲高清精品| 精品国产一区二区三区四区第35| 国产一区二区在线av高清观看| 一级作爱视频免费观看| 欧美乱码精品一区二区三区| avwww免费| 日韩人妻精品一区2区三区| 国产精品永久免费网站| 999精品在线视频| 成年人免费黄色播放视频| 别揉我奶头~嗯~啊~动态视频| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 国产一区二区在线av高清观看| 黄色片一级片一级黄色片| 亚洲午夜理论影院| 波多野结衣av一区二区av| 亚洲七黄色美女视频| 人妻久久中文字幕网| 精品福利观看| 99国产精品免费福利视频| 国产av又大| 国产激情欧美一区二区| 亚洲国产精品999在线| 亚洲精品中文字幕一二三四区| 91麻豆av在线| 国产无遮挡羞羞视频在线观看| 亚洲精品美女久久久久99蜜臀| 99在线视频只有这里精品首页| 久久国产亚洲av麻豆专区| 久久人人97超碰香蕉20202| 成在线人永久免费视频| 亚洲一区高清亚洲精品| 欧美成人午夜精品| 久久久水蜜桃国产精品网| 99热国产这里只有精品6| 9191精品国产免费久久| 法律面前人人平等表现在哪些方面| 老司机在亚洲福利影院| 免费高清视频大片| 国产成人一区二区三区免费视频网站| www日本在线高清视频| 久热这里只有精品99| 十八禁人妻一区二区| 亚洲欧美日韩高清在线视频| 亚洲美女黄片视频| 在线观看免费日韩欧美大片| 真人做人爱边吃奶动态| 国产又色又爽无遮挡免费看| 午夜精品久久久久久毛片777| 亚洲午夜理论影院| 精品国产美女av久久久久小说| 丝袜人妻中文字幕| 国产1区2区3区精品| 人人妻,人人澡人人爽秒播| 免费女性裸体啪啪无遮挡网站| 婷婷丁香在线五月| 一级黄色大片毛片| 色综合站精品国产| 国产精品 欧美亚洲| 不卡一级毛片| 免费在线观看完整版高清| 国产成人精品在线电影| 午夜福利影视在线免费观看| 亚洲精品成人av观看孕妇| 亚洲精品久久成人aⅴ小说| 久久久国产精品麻豆| 99香蕉大伊视频| www.www免费av| 国产极品粉嫩免费观看在线| 欧美成人免费av一区二区三区| 久久国产亚洲av麻豆专区| 久久久久久人人人人人| 999久久久国产精品视频| 国产在线观看jvid| 啦啦啦 在线观看视频| 交换朋友夫妻互换小说| 首页视频小说图片口味搜索| av天堂久久9| 亚洲一区二区三区色噜噜 | 日韩视频一区二区在线观看| 免费在线观看黄色视频的| av国产精品久久久久影院| 国产一区二区三区视频了| 老汉色av国产亚洲站长工具| 午夜福利在线观看吧| 如日韩欧美国产精品一区二区三区| 黄片播放在线免费| 9热在线视频观看99| 久久欧美精品欧美久久欧美| 欧美日韩亚洲综合一区二区三区_| 久久精品国产99精品国产亚洲性色 | 天堂俺去俺来也www色官网| av欧美777| 又大又爽又粗| 免费在线观看日本一区| 黑人巨大精品欧美一区二区蜜桃| 香蕉丝袜av| av片东京热男人的天堂| 两个人看的免费小视频| 亚洲欧美一区二区三区久久| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 18禁黄网站禁片午夜丰满| 黑人操中国人逼视频| 亚洲第一av免费看| 视频在线观看一区二区三区| 久久影院123| 亚洲 欧美 日韩 在线 免费| 巨乳人妻的诱惑在线观看| 天天添夜夜摸| 久久久国产成人精品二区 | 国产亚洲精品久久久久5区| 在线观看免费午夜福利视频| 久久人人爽av亚洲精品天堂| 夜夜看夜夜爽夜夜摸 | 美女 人体艺术 gogo| 亚洲精品av麻豆狂野| 欧美日韩黄片免| 亚洲中文av在线| 夜夜爽天天搞| 免费av中文字幕在线| 日本黄色视频三级网站网址| 两个人免费观看高清视频| 日本免费a在线| 国产亚洲欧美精品永久| 日本黄色视频三级网站网址| 国产欧美日韩精品亚洲av| 高清黄色对白视频在线免费看| 伦理电影免费视频| 热re99久久精品国产66热6| 亚洲成人精品中文字幕电影 | 成年人黄色毛片网站| 欧美大码av| a级毛片黄视频| 日本 av在线| 欧美在线一区亚洲| 男女高潮啪啪啪动态图| 亚洲国产精品合色在线| 波多野结衣av一区二区av| 免费观看人在逋| 首页视频小说图片口味搜索| 欧美乱色亚洲激情| 亚洲自拍偷在线| 黄频高清免费视频| 老司机午夜福利在线观看视频| 香蕉久久夜色| 国产男靠女视频免费网站| 香蕉丝袜av| 国产亚洲欧美精品永久| 在线观看66精品国产| 亚洲精品一二三| 最近最新免费中文字幕在线| 咕卡用的链子| 99国产精品免费福利视频| 欧美成人免费av一区二区三区| 国产高清视频在线播放一区| 麻豆av在线久日| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 淫妇啪啪啪对白视频| 国产xxxxx性猛交| 亚洲色图 男人天堂 中文字幕| 不卡av一区二区三区| 天堂√8在线中文| 精品久久久久久成人av| 国产成人系列免费观看| 亚洲欧美激情综合另类| 少妇 在线观看| 嫩草影院精品99| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 成人手机av| 免费在线观看视频国产中文字幕亚洲| 免费高清在线观看日韩| 看片在线看免费视频| 一区福利在线观看| 亚洲人成伊人成综合网2020| 午夜精品国产一区二区电影| 大型av网站在线播放| 国产成人精品久久二区二区免费| 天堂动漫精品| 亚洲精品国产色婷婷电影| 久久国产乱子伦精品免费另类| 国产欧美日韩一区二区精品| 精品国产亚洲在线| 亚洲黑人精品在线| 高清欧美精品videossex| 他把我摸到了高潮在线观看| 日韩欧美在线二视频| 精品人妻1区二区| 一区二区三区精品91| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| av电影中文网址| 久久久水蜜桃国产精品网| 高清在线国产一区| 日韩成人在线观看一区二区三区| 免费在线观看亚洲国产| 91字幕亚洲| 国产黄色免费在线视频| 欧美精品啪啪一区二区三区| 性色av乱码一区二区三区2| 啦啦啦在线免费观看视频4| 亚洲成人国产一区在线观看| 在线播放国产精品三级| av在线播放免费不卡| 精品久久久久久电影网| 97碰自拍视频| 国产成人av教育| 亚洲一区二区三区色噜噜 | 一本大道久久a久久精品| 亚洲熟妇熟女久久| 免费观看精品视频网站| 99精品欧美一区二区三区四区| 午夜久久久在线观看| 乱人伦中国视频| 久久久久久久久久久久大奶| 国产精品电影一区二区三区| 久久狼人影院| 欧美成狂野欧美在线观看| 91在线观看av| 午夜福利欧美成人| 在线观看午夜福利视频| 真人做人爱边吃奶动态| 亚洲激情在线av| 女警被强在线播放| 超碰成人久久| 80岁老熟妇乱子伦牲交| 久9热在线精品视频| 搡老熟女国产l中国老女人| avwww免费| 欧美人与性动交α欧美精品济南到| 午夜日韩欧美国产| 99re在线观看精品视频| 波多野结衣一区麻豆| 在线国产一区二区在线| 久久精品国产亚洲av香蕉五月| 桃红色精品国产亚洲av| 99久久99久久久精品蜜桃| 在线观看免费视频网站a站| 99在线视频只有这里精品首页| 无限看片的www在线观看| 18禁国产床啪视频网站| 亚洲人成伊人成综合网2020| 中文亚洲av片在线观看爽| videosex国产| 亚洲av片天天在线观看| 90打野战视频偷拍视频| 中文字幕高清在线视频| 国产高清激情床上av| 69精品国产乱码久久久| a级毛片黄视频| 免费在线观看完整版高清| 女人被狂操c到高潮| 欧美精品啪啪一区二区三区| 欧美激情 高清一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 色老头精品视频在线观看| 啦啦啦 在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 精品电影一区二区在线| 电影成人av| 久久久久久久久中文| 波多野结衣高清无吗| 亚洲aⅴ乱码一区二区在线播放 | 无限看片的www在线观看| 人人妻人人爽人人添夜夜欢视频| 日韩一卡2卡3卡4卡2021年| 又黄又爽又免费观看的视频| 欧美日本中文国产一区发布| 大码成人一级视频| 丰满饥渴人妻一区二区三| 国内毛片毛片毛片毛片毛片| 最新美女视频免费是黄的| 日韩免费av在线播放| а√天堂www在线а√下载| 日韩大尺度精品在线看网址 | 在线免费观看的www视频| 在线av久久热| 亚洲成人久久性| 18美女黄网站色大片免费观看| 久久 成人 亚洲| avwww免费| 国产成人av教育| 大型黄色视频在线免费观看| 精品国产一区二区久久| 国产男靠女视频免费网站| aaaaa片日本免费| 嫩草影视91久久| 午夜老司机福利片| 久久人人精品亚洲av| 在线av久久热| 在线观看一区二区三区| 亚洲精品国产色婷婷电影| 精品国产美女av久久久久小说| 黄色丝袜av网址大全| 久久天堂一区二区三区四区| 淫妇啪啪啪对白视频| 亚洲成av片中文字幕在线观看| 夜夜躁狠狠躁天天躁| 国产欧美日韩精品亚洲av| 精品国内亚洲2022精品成人| 国产三级黄色录像| 午夜91福利影院| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看 | 老司机亚洲免费影院| 亚洲av熟女| 国产精品国产高清国产av| 久久久久久久久久久久大奶| 夜夜看夜夜爽夜夜摸 | 9热在线视频观看99| 亚洲国产欧美一区二区综合| 女同久久另类99精品国产91| 午夜精品在线福利| 自线自在国产av| 很黄的视频免费| 国产男靠女视频免费网站| 女人被狂操c到高潮| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 激情在线观看视频在线高清| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 十八禁网站免费在线| 精品久久久久久成人av| 婷婷精品国产亚洲av在线| 青草久久国产| 女性被躁到高潮视频| 久久精品人人爽人人爽视色| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品久久二区二区91| 一级毛片高清免费大全| 1024香蕉在线观看| 99久久综合精品五月天人人| 老司机靠b影院| 视频在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 美女高潮到喷水免费观看| 91大片在线观看| 男女下面进入的视频免费午夜 | 亚洲人成网站在线播放欧美日韩| 俄罗斯特黄特色一大片| 精品国产乱子伦一区二区三区| 80岁老熟妇乱子伦牲交| 在线av久久热| 嫩草影视91久久| 99精品久久久久人妻精品| 天天躁夜夜躁狠狠躁躁| 国产单亲对白刺激| 日韩欧美一区二区三区在线观看| 精品国产乱子伦一区二区三区| 女性生殖器流出的白浆| 99国产精品一区二区三区| 一级片免费观看大全| 黑人巨大精品欧美一区二区mp4| 国产色视频综合| 亚洲专区中文字幕在线| 在线播放国产精品三级| 欧美在线黄色| 久久午夜综合久久蜜桃| 成人免费观看视频高清| 欧美 亚洲 国产 日韩一| 亚洲国产欧美一区二区综合| 亚洲精品在线美女| 日本wwww免费看| 久久久国产欧美日韩av| 高清在线国产一区| 国产黄a三级三级三级人| netflix在线观看网站| 亚洲av第一区精品v没综合| 一级a爱视频在线免费观看| 亚洲五月婷婷丁香| 乱人伦中国视频| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区日韩欧美中文字幕| 看免费av毛片| 俄罗斯特黄特色一大片| 一二三四在线观看免费中文在| tocl精华| 国产野战对白在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品影院6| 欧美乱码精品一区二区三区| a级毛片在线看网站| 国产一区二区三区视频了| 精品一区二区三区四区五区乱码| 日韩精品免费视频一区二区三区| 我的亚洲天堂| 精品国产乱码久久久久久男人| 777久久人妻少妇嫩草av网站| 久久午夜亚洲精品久久| 亚洲精品一区av在线观看| 欧美另类亚洲清纯唯美| 老汉色av国产亚洲站长工具| 又黄又粗又硬又大视频| 麻豆av在线久日| 久久久久久人人人人人| 99久久综合精品五月天人人| 欧美黑人精品巨大| 欧美日韩亚洲国产一区二区在线观看| 午夜福利欧美成人| 亚洲精品国产精品久久久不卡| 成年人黄色毛片网站| 国产精品亚洲一级av第二区| 亚洲五月天丁香| 又黄又爽又免费观看的视频| 欧美日韩亚洲综合一区二区三区_| 国产成+人综合+亚洲专区| 男女之事视频高清在线观看| 窝窝影院91人妻| 欧美老熟妇乱子伦牲交| 婷婷丁香在线五月| av天堂在线播放| 亚洲色图 男人天堂 中文字幕| 757午夜福利合集在线观看| 国产亚洲欧美在线一区二区| 午夜久久久在线观看| 亚洲国产欧美日韩在线播放| 久久人妻熟女aⅴ| 国产av一区在线观看免费| 亚洲国产欧美日韩在线播放| 91av网站免费观看| 成人18禁在线播放| 嫩草影视91久久| 久久久精品欧美日韩精品| 久久久国产一区二区| 国产99白浆流出| 国产亚洲精品久久久久5区| 中文字幕精品免费在线观看视频| 亚洲av第一区精品v没综合| 热99国产精品久久久久久7| 精品免费久久久久久久清纯| 久久国产亚洲av麻豆专区| 亚洲国产精品999在线| 亚洲五月天丁香| 亚洲精品国产色婷婷电影| 少妇粗大呻吟视频| 午夜福利在线观看吧| 国产亚洲精品一区二区www| 一级a爱片免费观看的视频| 久久九九热精品免费| 国产成人精品久久二区二区免费| 精品一区二区三卡| 精品无人区乱码1区二区| 一级a爱视频在线免费观看| 免费观看人在逋| 9色porny在线观看| 国产精品综合久久久久久久免费 | 成年版毛片免费区| avwww免费| 午夜老司机福利片| 亚洲成人久久性| 午夜亚洲福利在线播放| 国内久久婷婷六月综合欲色啪| 亚洲精品成人av观看孕妇| www.999成人在线观看| 丰满迷人的少妇在线观看| 久久人人精品亚洲av| 欧美色视频一区免费| 日韩免费av在线播放| 97超级碰碰碰精品色视频在线观看| 久久国产精品男人的天堂亚洲| 欧美另类亚洲清纯唯美| 两个人免费观看高清视频| 黄色丝袜av网址大全| 长腿黑丝高跟| 在线国产一区二区在线| 久久中文字幕人妻熟女| 又黄又爽又免费观看的视频| 国产精华一区二区三区| 久久婷婷成人综合色麻豆| 99香蕉大伊视频| 12—13女人毛片做爰片一| 免费不卡黄色视频| 国产主播在线观看一区二区| 国产成人av激情在线播放| 精品久久久久久久毛片微露脸| 亚洲国产精品合色在线| 久久精品亚洲av国产电影网| 欧美+亚洲+日韩+国产| 久久久国产成人精品二区 | 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区综合在线观看| 好看av亚洲va欧美ⅴa在| 日韩av在线大香蕉| 亚洲欧美激情在线| 国产成人av教育| 精品国产一区二区三区四区第35| 级片在线观看| 88av欧美| 国产1区2区3区精品| 欧美日韩乱码在线| 精品少妇一区二区三区视频日本电影| 久久久国产成人免费| 成人三级黄色视频| 最新美女视频免费是黄的| 亚洲精品粉嫩美女一区| 91老司机精品| 一个人免费在线观看的高清视频| 亚洲av五月六月丁香网| 99精国产麻豆久久婷婷| 久久国产精品影院| 亚洲片人在线观看| 老司机福利观看| 国产欧美日韩一区二区三| 国产高清国产精品国产三级| 99在线人妻在线中文字幕| 18禁黄网站禁片午夜丰满| 三上悠亚av全集在线观看| 日本免费a在线| 精品一区二区三区av网在线观看| 欧美激情 高清一区二区三区| 黄色怎么调成土黄色| 久久久精品欧美日韩精品| 亚洲va日本ⅴa欧美va伊人久久| 丝袜人妻中文字幕| ponron亚洲| 热99re8久久精品国产| 男女午夜视频在线观看| 老司机福利观看| 国产精品国产av在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产成年人精品一区二区 | 国产精品香港三级国产av潘金莲| 高清av免费在线| 男女做爰动态图高潮gif福利片 | 国产精品 国内视频| 日韩中文字幕欧美一区二区|