• <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)
    牽掛
    国产精品野战在线观看| 69av精品久久久久久| 欧美精品一区二区大全| 中文字幕精品亚洲无线码一区| 91久久精品国产一区二区三区| 欧美一区二区精品小视频在线| 欧美+日韩+精品| 最近最新中文字幕大全电影3| 久久国内精品自在自线图片| 青青草视频在线视频观看| 婷婷六月久久综合丁香| 日本与韩国留学比较| ponron亚洲| 国产av码专区亚洲av| 中文字幕制服av| 在现免费观看毛片| 两个人视频免费观看高清| 亚洲自偷自拍三级| 一区二区三区乱码不卡18| 精品国产三级普通话版| 久久精品91蜜桃| 成人高潮视频无遮挡免费网站| 特级一级黄色大片| 午夜视频国产福利| av黄色大香蕉| 亚洲国产最新在线播放| 亚洲av一区综合| 国产av不卡久久| 午夜a级毛片| 最新中文字幕久久久久| 色综合色国产| 免费av观看视频| 日韩在线高清观看一区二区三区| 国产淫语在线视频| 欧美成人一区二区免费高清观看| 久久韩国三级中文字幕| 国内精品宾馆在线| 男人的好看免费观看在线视频| 亚洲av福利一区| 日日撸夜夜添| 久久久久精品久久久久真实原创| 2022亚洲国产成人精品| 成人毛片a级毛片在线播放| 国产伦一二天堂av在线观看| 欧美变态另类bdsm刘玥| 99久久成人亚洲精品观看| 男人狂女人下面高潮的视频| 国产精品国产三级国产专区5o | 日韩亚洲欧美综合| 嫩草影院精品99| 中文精品一卡2卡3卡4更新| 久久精品夜夜夜夜夜久久蜜豆| 日本与韩国留学比较| 免费av毛片视频| 七月丁香在线播放| 男人舔女人下体高潮全视频| 国产亚洲5aaaaa淫片| 三级国产精品欧美在线观看| 精品国产一区二区三区久久久樱花 | 九九久久精品国产亚洲av麻豆| 只有这里有精品99| 一区二区三区高清视频在线| 插逼视频在线观看| 可以在线观看毛片的网站| 久久久久久久国产电影| 中文字幕制服av| 日日干狠狠操夜夜爽| 变态另类丝袜制服| 日韩欧美三级三区| 我的女老师完整版在线观看| 成人国产麻豆网| 久久精品夜夜夜夜夜久久蜜豆| 日产精品乱码卡一卡2卡三| 日韩人妻高清精品专区| 国产一级毛片七仙女欲春2| 欧美三级亚洲精品| 免费观看在线日韩| 午夜福利在线观看免费完整高清在| 亚洲无线观看免费| 少妇的逼水好多| 午夜精品一区二区三区免费看| 大又大粗又爽又黄少妇毛片口| 久久久久久久久久久丰满| 非洲黑人性xxxx精品又粗又长| 国产免费视频播放在线视频 | 亚洲av成人av| 国产精品av视频在线免费观看| eeuss影院久久| 午夜福利在线在线| 国产av在哪里看| 又粗又爽又猛毛片免费看| 午夜福利高清视频| 22中文网久久字幕| 一区二区三区高清视频在线| 欧美日韩综合久久久久久| 亚洲最大成人中文| 久久久久久久久中文| 国产免费又黄又爽又色| 国产高清三级在线| 国产精品麻豆人妻色哟哟久久 | 国产精品伦人一区二区| 久久精品国产亚洲av天美| videossex国产| 欧美成人a在线观看| 精品久久久久久久末码| 亚洲五月天丁香| 久久欧美精品欧美久久欧美| 亚洲成人av在线免费| 国产伦在线观看视频一区| 免费看光身美女| 精品酒店卫生间| .国产精品久久| 九九在线视频观看精品| 在线观看66精品国产| 国内精品一区二区在线观看| www.av在线官网国产| 岛国毛片在线播放| 国产午夜精品一二区理论片| 精品少妇黑人巨大在线播放 | av又黄又爽大尺度在线免费看 | 毛片女人毛片| 青春草亚洲视频在线观看| av在线天堂中文字幕| 久久久久九九精品影院| 午夜精品在线福利| 国产精品熟女久久久久浪| 亚洲国产色片| 国产乱来视频区| 好男人在线观看高清免费视频| 久久久色成人| 午夜爱爱视频在线播放| 精品人妻熟女av久视频| 一级黄片播放器| 免费看av在线观看网站| 久久人人爽人人片av| 日韩一区二区三区影片| 久久99热6这里只有精品| 三级国产精品片| 欧美丝袜亚洲另类| 少妇熟女欧美另类| 日韩一区二区三区影片| 女的被弄到高潮叫床怎么办| 亚洲国产成人一精品久久久| 国产精品蜜桃在线观看| 99久国产av精品| 看黄色毛片网站| 午夜福利在线观看免费完整高清在| 国内揄拍国产精品人妻在线| 日韩中字成人| 26uuu在线亚洲综合色| 国产伦精品一区二区三区视频9| 国产欧美另类精品又又久久亚洲欧美| 欧美精品国产亚洲| 免费大片18禁| 一二三四中文在线观看免费高清| 亚洲国产最新在线播放| 国产乱来视频区| 少妇人妻精品综合一区二区| 亚洲精品456在线播放app| 亚洲av电影不卡..在线观看| 成人亚洲精品av一区二区| 中文字幕人妻熟人妻熟丝袜美| 免费观看精品视频网站| 久久精品熟女亚洲av麻豆精品 | 国产精华一区二区三区| 秋霞在线观看毛片| 婷婷色av中文字幕| 夜夜爽夜夜爽视频| 尾随美女入室| av天堂中文字幕网| 男人舔奶头视频| 亚洲欧美日韩无卡精品| 久久精品综合一区二区三区| 一级毛片电影观看 | 性插视频无遮挡在线免费观看| 国产亚洲精品久久久com| 国产中年淑女户外野战色| 国产精品蜜桃在线观看| 国产 一区 欧美 日韩| 精品人妻熟女av久视频| 国产美女午夜福利| 一区二区三区乱码不卡18| 国产人妻一区二区三区在| 国产v大片淫在线免费观看| 国产精品永久免费网站| 色5月婷婷丁香| 黄色一级大片看看| 亚洲电影在线观看av| 国产三级中文精品| 哪个播放器可以免费观看大片| 免费av毛片视频| 长腿黑丝高跟| 偷拍熟女少妇极品色| 99久久无色码亚洲精品果冻| 毛片女人毛片| 国产乱来视频区| 国产综合懂色| 内射极品少妇av片p| 22中文网久久字幕| 黄片wwwwww| 国产精品人妻久久久影院| 免费大片18禁| av女优亚洲男人天堂| 欧美成人精品欧美一级黄| 中文字幕久久专区| 国产黄片美女视频| 亚洲成人精品中文字幕电影| 一边亲一边摸免费视频| 一个人观看的视频www高清免费观看| 人妻制服诱惑在线中文字幕| 国产精品福利在线免费观看| 亚洲av一区综合| 午夜爱爱视频在线播放| 成年av动漫网址| 成人美女网站在线观看视频| 久久鲁丝午夜福利片| 午夜福利在线在线| 国产精品,欧美在线| 99久久中文字幕三级久久日本| 欧美性猛交黑人性爽| 一边摸一边抽搐一进一小说| a级毛色黄片| av在线老鸭窝| 日本熟妇午夜| 纵有疾风起免费观看全集完整版 | 国产精品永久免费网站| 中文字幕av在线有码专区| 一级黄色大片毛片| 国产乱人视频| 七月丁香在线播放| 免费观看性生交大片5| 国产在视频线精品| 久久久久久久国产电影| 赤兔流量卡办理| 亚洲中文字幕一区二区三区有码在线看| 国产精品国产三级专区第一集| 精品国产露脸久久av麻豆 | 十八禁国产超污无遮挡网站| 欧美激情国产日韩精品一区| 国产爱豆传媒在线观看| 欧美bdsm另类| 国产视频首页在线观看| 18禁裸乳无遮挡免费网站照片| 免费看av在线观看网站| 青春草视频在线免费观看| 亚洲色图av天堂| 哪个播放器可以免费观看大片| 日韩欧美三级三区| 日韩成人av中文字幕在线观看| 干丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 中文欧美无线码| 国产精品人妻久久久影院| 免费看美女性在线毛片视频| 免费观看的影片在线观看| 精品久久久久久成人av| 久久精品熟女亚洲av麻豆精品 | 在现免费观看毛片| 久久精品国产亚洲av天美| 日韩在线高清观看一区二区三区| 波野结衣二区三区在线| 国产精品av视频在线免费观看| 国产美女午夜福利| 久久精品熟女亚洲av麻豆精品 | 日韩一区二区三区影片| 啦啦啦韩国在线观看视频| 看十八女毛片水多多多| 久久这里只有精品中国| 成人无遮挡网站| 国产在线一区二区三区精 | 精品人妻视频免费看| 色综合站精品国产| 亚洲精品国产av成人精品| 亚洲欧美精品综合久久99| 天堂av国产一区二区熟女人妻| 亚洲国产高清在线一区二区三| www日本黄色视频网| 最近中文字幕高清免费大全6| 亚洲欧美精品综合久久99| 白带黄色成豆腐渣| av国产免费在线观看| 亚洲最大成人av| 好男人在线观看高清免费视频| 精品99又大又爽又粗少妇毛片| 秋霞伦理黄片| 精品不卡国产一区二区三区| 成人三级黄色视频| 国产一区亚洲一区在线观看| av在线播放精品| 一级黄片播放器| 国产高清三级在线| 人人妻人人澡人人爽人人夜夜 | 欧美潮喷喷水| 美女脱内裤让男人舔精品视频| 亚洲成人中文字幕在线播放| 日本猛色少妇xxxxx猛交久久| 欧美xxxx黑人xx丫x性爽| 国产单亲对白刺激| 波多野结衣巨乳人妻| 永久网站在线| 免费av毛片视频| 久久精品久久精品一区二区三区| 最近的中文字幕免费完整| 精品久久久噜噜| 99久久九九国产精品国产免费| 2022亚洲国产成人精品| 精品久久久久久久久av| 久久久久九九精品影院| 精品99又大又爽又粗少妇毛片| 伦精品一区二区三区| 午夜福利在线在线| 亚洲内射少妇av| 麻豆成人午夜福利视频| 国产日韩欧美在线精品| 你懂的网址亚洲精品在线观看 | 超碰97精品在线观看| 久久精品国产亚洲av涩爱| 国产精品一区二区三区四区久久| 免费大片18禁| 中文字幕精品亚洲无线码一区| 联通29元200g的流量卡| 亚洲中文字幕日韩| 精品一区二区免费观看| 女人十人毛片免费观看3o分钟| 99热这里只有是精品在线观看| 国产在线男女| 又粗又爽又猛毛片免费看| 久久精品影院6| 在线免费十八禁| 精品国内亚洲2022精品成人| 男的添女的下面高潮视频| 久久久久久久午夜电影| 麻豆成人午夜福利视频| 日韩视频在线欧美| 国产老妇女一区| 99热网站在线观看| 99国产精品一区二区蜜桃av| 老司机福利观看| 不卡视频在线观看欧美| 日韩欧美精品v在线| 少妇裸体淫交视频免费看高清| 久久精品久久久久久久性| 日韩在线高清观看一区二区三区| 午夜日本视频在线| 亚洲国产欧美人成| 简卡轻食公司| 亚洲丝袜综合中文字幕| 身体一侧抽搐| 国产成人精品婷婷| 国产黄片视频在线免费观看| 又爽又黄无遮挡网站| 久久久久久久久中文| 男女啪啪激烈高潮av片| 床上黄色一级片| 亚洲欧美一区二区三区国产| 看片在线看免费视频| 麻豆久久精品国产亚洲av| 变态另类丝袜制服| 麻豆成人午夜福利视频| 青春草视频在线免费观看| 纵有疾风起免费观看全集完整版 | 国语自产精品视频在线第100页| 日本wwww免费看| 亚洲国产最新在线播放| 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 爱豆传媒免费全集在线观看| 青春草国产在线视频| 欧美不卡视频在线免费观看| 日本色播在线视频| 一夜夜www| 欧美成人午夜免费资源| .国产精品久久| 高清日韩中文字幕在线| 18禁裸乳无遮挡免费网站照片| 国产综合懂色| 永久网站在线| 亚洲国产欧洲综合997久久,| 中文乱码字字幕精品一区二区三区 | 亚洲国产欧美在线一区| 久久久国产成人免费| a级一级毛片免费在线观看| 久久精品人妻少妇| 国产黄片视频在线免费观看| 精品久久久久久成人av| .国产精品久久| 美女黄网站色视频| 亚洲最大成人手机在线| 国产色婷婷99| 黄片无遮挡物在线观看| 亚洲国产精品成人综合色| 老师上课跳d突然被开到最大视频| 大话2 男鬼变身卡| a级毛片免费高清观看在线播放| 久久精品久久久久久久性| 人妻制服诱惑在线中文字幕| 国内揄拍国产精品人妻在线| 女人十人毛片免费观看3o分钟| 丝袜喷水一区| 午夜日本视频在线| 欧美zozozo另类| 国产高清三级在线| 日韩欧美 国产精品| 乱人视频在线观看| 亚洲国产成人一精品久久久| .国产精品久久| av又黄又爽大尺度在线免费看 | 天堂影院成人在线观看| 永久免费av网站大全| 美女xxoo啪啪120秒动态图| 内地一区二区视频在线| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 一卡2卡三卡四卡精品乱码亚洲| 欧美三级亚洲精品| 国产精品精品国产色婷婷| 欧美另类亚洲清纯唯美| 欧美三级亚洲精品| 亚洲国产欧美人成| 一级av片app| 最近中文字幕2019免费版| 中文字幕熟女人妻在线| 日本一本二区三区精品| 国产三级在线视频| 最后的刺客免费高清国语| 亚洲欧美成人精品一区二区| 精品国产三级普通话版| 1024手机看黄色片| 特级一级黄色大片| 国产成人免费观看mmmm| 少妇熟女欧美另类| 国语自产精品视频在线第100页| 久久久久久国产a免费观看| 亚洲国产精品久久男人天堂| 日韩av不卡免费在线播放| 毛片一级片免费看久久久久| 色5月婷婷丁香| 天堂√8在线中文| 一级爰片在线观看| 欧美bdsm另类| 秋霞在线观看毛片| 中文亚洲av片在线观看爽| 91久久精品国产一区二区三区| 欧美高清性xxxxhd video| 一级黄色大片毛片| 91久久精品国产一区二区三区| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 亚洲av免费高清在线观看| 少妇熟女aⅴ在线视频| 精品久久久久久久久av| 九九爱精品视频在线观看| 男人舔奶头视频| 嫩草影院入口| 两性午夜刺激爽爽歪歪视频在线观看| 国产毛片a区久久久久| 日韩三级伦理在线观看| 一二三四中文在线观看免费高清| 亚洲熟妇中文字幕五十中出| 免费看日本二区| 亚洲精品aⅴ在线观看| 色播亚洲综合网| 中文字幕av在线有码专区| 日韩三级伦理在线观看| av在线老鸭窝| 日韩大片免费观看网站 | 色5月婷婷丁香| 国产亚洲最大av| 国产精品野战在线观看| 国产高清视频在线观看网站| 级片在线观看| 日本午夜av视频| 韩国高清视频一区二区三区| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 18+在线观看网站| 亚洲av.av天堂| 一级av片app| 日产精品乱码卡一卡2卡三| 午夜免费激情av| 免费看光身美女| 久久精品综合一区二区三区| 色尼玛亚洲综合影院| 亚洲乱码一区二区免费版| 亚洲va在线va天堂va国产| 日本熟妇午夜| 99九九线精品视频在线观看视频| 精品久久久噜噜| 日韩大片免费观看网站 | 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| 午夜日本视频在线| 国产亚洲av嫩草精品影院| 三级经典国产精品| 国产高清有码在线观看视频| 日日摸夜夜添夜夜添av毛片| 久久久久久国产a免费观看| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱久久久久久| 亚洲不卡免费看| 日韩欧美精品v在线| 99热全是精品| 午夜福利网站1000一区二区三区| 亚洲欧美日韩卡通动漫| 尾随美女入室| 男女边吃奶边做爰视频| 一级爰片在线观看| 国产精品久久久久久精品电影| 久99久视频精品免费| 免费观看在线日韩| 亚洲国产最新在线播放| 欧美激情久久久久久爽电影| 激情 狠狠 欧美| 日韩一区二区视频免费看| 日韩欧美三级三区| 国产精品蜜桃在线观看| 欧美一区二区精品小视频在线| 男人的好看免费观看在线视频| 波野结衣二区三区在线| 狠狠狠狠99中文字幕| 亚洲欧美成人精品一区二区| 亚洲精品国产成人久久av| 成人综合一区亚洲| 亚洲av成人精品一二三区| 啦啦啦啦在线视频资源| 国内少妇人妻偷人精品xxx网站| 禁无遮挡网站| 成人亚洲欧美一区二区av| 校园人妻丝袜中文字幕| 亚洲人成网站在线观看播放| 精华霜和精华液先用哪个| 最近手机中文字幕大全| 女人久久www免费人成看片 | 精品欧美国产一区二区三| av福利片在线观看| 免费观看性生交大片5| 国产伦在线观看视频一区| 综合色av麻豆| 国产精品野战在线观看| 一边亲一边摸免费视频| 国产精品麻豆人妻色哟哟久久 | 精品国产三级普通话版| 女的被弄到高潮叫床怎么办| 日本午夜av视频| 麻豆av噜噜一区二区三区| ponron亚洲| 欧美色视频一区免费| 寂寞人妻少妇视频99o| 国产精华一区二区三区| 少妇高潮的动态图| 日本爱情动作片www.在线观看| 高清视频免费观看一区二区 | 床上黄色一级片| 亚洲不卡免费看| 国产不卡一卡二| 国产精品久久久久久久久免| 精品人妻视频免费看| 麻豆av噜噜一区二区三区| 欧美日韩在线观看h| 亚洲精品456在线播放app| 国产精品福利在线免费观看| 国产精品国产三级国产专区5o | 亚洲综合精品二区| 精品国产一区二区三区久久久樱花 | 国产伦理片在线播放av一区| 亚洲欧洲日产国产| 免费观看精品视频网站| 国产亚洲91精品色在线| 亚洲av成人精品一区久久| 日本wwww免费看| 色尼玛亚洲综合影院| 高清日韩中文字幕在线| 久久人妻av系列| 在线免费观看的www视频| 欧美97在线视频| АⅤ资源中文在线天堂| 国产成人精品一,二区| 亚洲精品色激情综合| 成人综合一区亚洲| 长腿黑丝高跟| 99久久精品一区二区三区| 亚洲精品国产成人久久av| 国产精品不卡视频一区二区| 听说在线观看完整版免费高清| 中文资源天堂在线| 我要看日韩黄色一级片| 久久久久久久久中文| 欧美一区二区国产精品久久精品| 午夜免费男女啪啪视频观看| 色哟哟·www| 赤兔流量卡办理| 九九热线精品视视频播放| 国产午夜精品一二区理论片| 日韩强制内射视频| 在线播放无遮挡| 激情 狠狠 欧美| 欧美+日韩+精品| 亚州av有码| videos熟女内射| 长腿黑丝高跟| 美女被艹到高潮喷水动态| 国产精品日韩av在线免费观看| a级一级毛片免费在线观看| 嫩草影院精品99| 美女被艹到高潮喷水动态| 男女国产视频网站| 91精品一卡2卡3卡4卡| 久久99热这里只频精品6学生 | 精华霜和精华液先用哪个| 亚洲美女搞黄在线观看| 日韩中字成人| 99久久人妻综合| 亚洲一级一片aⅴ在线观看| 欧美潮喷喷水|