基于WRF模式的暴雨天氣過程的數(shù)值模擬及診斷分析

袁成松，王秋云，包云軒，嚴(yán)明良

(1.江蘇省氣象科學(xué)研究所，江蘇南京210008;2.江蘇省農(nóng)業(yè)氣象重點(diǎn)實(shí)驗(yàn)室，江蘇南京210044;3.南京信息工程大學(xué)應(yīng)用氣象學(xué)院，江蘇南京210044)

0 引言

2009年9月17日江蘇省南部各城市都出現(xiàn)了大暴雨，這是進(jìn)入9月以來最強(qiáng)的一次降雨過程，大雨幾乎持續(xù)了一整天。這場暴雨來得突然，給生產(chǎn)和生活帶來諸多不便。由于這場暴雨發(fā)生在交通高峰期，因此，對江蘇南部地區(qū)的交通產(chǎn)生了巨大影響，使交通一度癱瘓。暴雨是影響人類的嚴(yán)重災(zāi)害性天氣之一，所以對暴雨的準(zhǔn)確預(yù)報(bào)顯得十分重要(Molinari and Dudek，1992;孫建華等，2004)。

目前暴雨的準(zhǔn)確預(yù)報(bào)難度較大，而暴雨的研究主要基于實(shí)測資料和數(shù)值模擬兩種途徑(張曉紅等，2009;張羽等，2010)。近年來，已經(jīng)有很多學(xué)者利用MM5或WRF中尺度氣象數(shù)值預(yù)報(bào)模式來模擬暴雨天氣過程，結(jié)果表明:利用中尺度數(shù)值模式對降雨天氣過程的模擬具有一定的能力(陳業(yè)國和農(nóng)孟松，2010;沈桐立等，2010)。在我國由于計(jì)算機(jī)條件的限制，利用高分辨率數(shù)值模式進(jìn)行暴雨等中尺度天氣現(xiàn)象的研究還比較困難，采用多重嵌套是在現(xiàn)有計(jì)算機(jī)條件下盡可能提高重點(diǎn)預(yù)報(bào)區(qū)域分辨率的有效途徑(Cui et al.，2007;史薈燕等，2009)。王建捷和李澤椿(2002)對1998年6月的一次梅雨鋒暴雨中尺度系統(tǒng)進(jìn)行了模擬和診斷分析;張小玲等(2002)研究表明，中尺度系統(tǒng)強(qiáng)烈發(fā)展引發(fā)了1998年7月20—21日武漢地區(qū)的突發(fā)性暴雨，并確認(rèn)了暴雨發(fā)生時(shí)期的中小尺度特征;王智等(2003)對一次西南渦及其伴隨低空急流的發(fā)展演變進(jìn)行了數(shù)值模擬，揭示了西南渦和低空急流與1999年初夏江淮流域多日持續(xù)性暴雨的關(guān)系。但是以上研究所采用的模式分辨率都較低，且著重分析某一大范圍區(qū)域暴雨的中尺度特征，并沒有針對暴雨對公路交通能見度變化的影響進(jìn)行探討。

本文利用新一代中尺度天氣研究和預(yù)報(bào)模式WRF，對2009年9月17日發(fā)生在江蘇南部地區(qū)、覆蓋滬寧高速公路的一次大暴雨天氣過程進(jìn)行了數(shù)值模擬研究，與過去的研究相比，分辨率大大提高。文中結(jié)合相關(guān)實(shí)測資料與模式輸出的物理量，分析了暴雨形成的環(huán)流背景和水汽、動(dòng)力、熱力條件等，從而為公路交通降雨預(yù)報(bào)提供了一些定量化的參考指標(biāo)。

1 雨情概況及大尺度環(huán)流背景

1.1 雨情概況

根據(jù)江蘇省地面氣象觀測網(wǎng)提供的實(shí)測降水資料，2009年9月17日江蘇大部普降大到暴雨。圖1給出了2009年9月17日08—19時(shí)南京站的逐時(shí)降水量，可知南京站12 h累積降雨量超過了50 mm。降水主要在17日08—09時(shí)、13—15時(shí)集中發(fā)生，下班高峰期(17時(shí)左右)也出現(xiàn)了較強(qiáng)降雨。

圖1 南京站2009年9月17日08—19時(shí)實(shí)測的逐時(shí)降雨量(單位:mm)Fig.1 Observed hourly precipitation at Nanjing station from 08:00 BST to 19:00 BST 17 September 2009(units:mm)

根據(jù)AWMS系統(tǒng)提供的實(shí)時(shí)監(jiān)測數(shù)據(jù)，圖2給出了2009年9月17日02時(shí)—18日02時(shí)24 h滬寧高速公路黃粟墅、河陽、常州北三站的水平能見度變化曲線。17日06時(shí)之前，河陽站水平能見度在3 000～3 500 m范圍內(nèi)波動(dòng)，此時(shí)未降雨;07時(shí)開始，降雨發(fā)生，能見度陡降至500 m以下，之后能見度小幅度緩慢回升，但仍維持在1 000 m以下;11時(shí)左右，雨勢轉(zhuǎn)小，能見度開始上升;13時(shí)左右，雨勢開始增強(qiáng)，能見度隨之降低;隨后能見度維持在1 000 m左右。黃粟墅站在17日09時(shí)、13時(shí)、17時(shí)出現(xiàn)的能見度低谷，正好對應(yīng)著3個(gè)降雨高峰期。常州北站的能見度也在多個(gè)時(shí)段出現(xiàn)了1 000 m以下的低值。由于各站所處的局地地形不同，出現(xiàn)降雨的時(shí)段不同，所以它們的能見度變化存在一些差異，但總體來說，能見度大小與降水強(qiáng)度存在著明顯的反相關(guān)關(guān)系。

1.2 大尺度環(huán)流背景

圖2 2009年9月17日02時(shí)—18日02時(shí)滬寧高速公路自動(dòng)氣象站觀測的水平能見度逐時(shí)變化(單位:m)Fig.2 Hourly change of horizontal visibility by auto-meteorological monitoring stations in the Shanghai-Nanjing expressway from 02:00 BST 17 to 02:00 BST 18 September 2009(units:m)

圖3給出了2009年9月17日02時(shí)500 hPa位勢高度場與850 hPa風(fēng)場?？梢?，暴雨發(fā)生前，高緯度烏拉爾山至西西伯利亞地區(qū)和俄羅斯遠(yuǎn)東至鄂霍次克海地區(qū)分別為強(qiáng)大穩(wěn)定的阻塞高壓所占據(jù)，中、東西伯利亞上空為一大低渦控制，極地冷空氣強(qiáng);中低緯地區(qū)，帶狀副熱帶高壓位于日本以南西北太平洋至我國東南沿海地區(qū)，脊線位于33°N附近。副高西北側(cè)盛行的西南低空急流將大量暖濕空氣不斷地輸送到江淮地區(qū)，與從極渦中分裂南下的干冷空氣相遇，從而為暴雨的產(chǎn)生提供了極為有利的條件。

圖3 2009年9月17日02時(shí)500 hPa位勢高度場(綠色實(shí)線;單位:gpm)與850 hPa風(fēng)矢量場(箭矢;單位:m·s－1)(彩色區(qū)表示風(fēng)速大于12 m·s－1;灰色區(qū)表示地形高于1 500 m)Fig.3 Geopotential height at 500 hPa(green solid lines;units:gpm)and wind vector at 850 hPa(arrows;units:m·s－1)at 02:00 BST 17 September 2009(color shaded areas denote the wind speed greater than 12 m·s－1and gray shaded areas the terrain above 1 500 m)

2 資料來源與模式簡介

2.1 資料來源

實(shí)測資料來源于江蘇省氣候中心提供的常規(guī)觀測資料以及江蘇省氣象科學(xué)研究所自主研發(fā)、并在滬寧高速公路上按每10 km布設(shè)1套的自動(dòng)天氣監(jiān)測系統(tǒng)(AWMS)每分鐘實(shí)時(shí)監(jiān)測的氣象數(shù)據(jù)，包括降水、大氣水平能見度、氣壓、氣溫、地溫、相對濕度、風(fēng)向、風(fēng)速等氣象要素(吳贊平等，2006)。

運(yùn)行WRF模式所用的NCEP資料來自NCAR/NCEP(美國國家大氣研究中心和國家環(huán)境預(yù)測中心)每日4個(gè)時(shí)次(02、08、14、20時(shí))的1°×1°氣象再分析資料。

2.2 WRF模式簡介

由美國國家大氣研究中心(NCAR)和美國環(huán)境預(yù)測中心(NCEP)等單位聯(lián)合開發(fā)的新一代中尺度數(shù)值天氣預(yù)報(bào)模式WRF(Weather Research and Forecast Model)，在目前世界的天氣預(yù)報(bào)業(yè)務(wù)和大氣科學(xué)研究中得到了廣泛應(yīng)用(章國材，2006)。該模式不僅具有獨(dú)特的數(shù)值化動(dòng)力框架，還提供了多種可供耦合的先進(jìn)的物理過程方案。與前一代中尺度數(shù)值預(yù)報(bào)模式MM5相比，WRF模式考慮了比較詳細(xì)的陸面過程，能描述不同下墊面的熱量、水分等要素的傳輸過程，并將其耦合到邊界層參數(shù)化方案中，加強(qiáng)了邊界層物理過程的模擬，進(jìn)而改善了整個(gè)模式的模擬性能(Joseph，2004;Janjic，2004;Paula et al.，2005)。

本文應(yīng)用了中尺度數(shù)值預(yù)報(bào)模式WRF2.2(Weather Research and Forecast Model for version 2.2)。模式采用雙重嵌套方式，水平范圍以(120.0°E，33.0°N)為中心，模擬區(qū)域包括了116～124°E、30～36°N，且覆蓋滬寧高速公路。粗網(wǎng)格格距為9 km，格點(diǎn)數(shù)為75×75;細(xì)網(wǎng)格格距為3 km，格點(diǎn)數(shù)為121×91;積分步長為60 s。在模式微物理過程的參數(shù)選擇上，粗細(xì)網(wǎng)格均采用Ferrier(new Eta)方案;長波輻射均選用RRTM方案;短波輻射均選用Dudhia方案;近地面方案選用Monin-Obukhov方案;陸面過程采用Noah陸面參數(shù)化方案，并耦合了城市冠層模型;積云參數(shù)化方案粗網(wǎng)格選擇Betts-Miller-Janjic方案，細(xì)網(wǎng)格不采用積云參數(shù)化方案;邊界層方案選擇MRF方案。模式模擬個(gè)例的時(shí)間為2009年9月17日02時(shí)(北京時(shí)間，下同)至18日02時(shí)共24 h，粗網(wǎng)格每3 h輸出一次結(jié)果，細(xì)網(wǎng)格每1 h輸出一次結(jié)果。

3 數(shù)值模擬結(jié)果與診斷分析

3.1 降水模擬結(jié)果的驗(yàn)證

圖4 2009年9月17日02時(shí)至18日02時(shí)的24 h累積降水量(單位:mm)a.實(shí)況;b.模擬Fig.4 Distribution of 24-hr accumulated precipitation from 02:00 BST 17 to 02:00 BST 18 September 2009(units:mm)a.observed;b.simulated

評價(jià)一個(gè)中尺度天氣模式的好壞，主要看它對降水的模擬和預(yù)報(bào)能力。圖4a、b分別是實(shí)測和模擬的2009年9月17日02時(shí)至18日02時(shí)的24 h累積降水量?？煽闯瞿Ｊ綄τ陰ё呦?、降雨落區(qū)和雨量大小的模擬比較好，與實(shí)況基本一致。但是在整個(gè)模擬區(qū)域內(nèi)，WRF模擬的降水量比實(shí)況偏低10 mm左右;且對于(120.0°E，32.4°N)附近的降雨落區(qū)，WRF模擬出幾個(gè)虛中心;此外，WRF未模擬出位于(120.3°E，32.0°N)的降水中心;對于(121.3°E，31.7°N)附近的暴雨中心，模擬的位置略偏西北，比實(shí)況范圍偏小;在(119.3°E，31.6°N)、(119.7°E，31.0°N)、(118.9°E，30.9°N)等附近的降水模擬都與實(shí)況存在一些微小差異。這可能是模式地形與實(shí)際情況存在誤差、地面實(shí)測站點(diǎn)有限且分布不均勻等因素有關(guān)?？傮w來說模擬結(jié)果基本體現(xiàn)了此次暴雨過程的主要特征，模擬效果較為理想。

3.2 暴雨的流場特征診斷

分析不同高度的水平流場可發(fā)現(xiàn)，長江中下游地區(qū)存在一個(gè)β中尺度渦旋的發(fā)展、移動(dòng)過程。2009年9月17日12時(shí)有一氣旋性渦旋在119.8～120.4°E、31.2～31.5°N生成(圖5a)，北緣緊逼滬寧高速公路中段，其水平尺度在70 km左右，垂直尺度為10 km左右(結(jié)合模式輸出的多層位勢高度場和水平流場估算得到)，形態(tài)比(即中尺度系統(tǒng)的垂直尺度H與水平尺度L之比H/L;陸漢城和楊國祥，2000)約為0.143，700 hPa高度渦旋區(qū)內(nèi)最大水平氣壓梯度約為2.0 hPa/(10 km)，此渦旋屬于典型的β中尺度渦旋，與此次暴雨過程降水的加強(qiáng)、維持關(guān)系密切;之后的4 h內(nèi)渦旋有明顯增強(qiáng)，并緩慢東移，于16時(shí)達(dá)到最強(qiáng)(圖5b)，此時(shí)，渦旋位于119.6～120.5°E、30.8～31.4°N，水平尺度達(dá)110 km以上，垂直尺度達(dá)11 km左右，形態(tài)比(H/L)約為0.1，700 hPa高度渦旋區(qū)內(nèi)最大水平氣壓梯度約為1.3 hPa/(10 km)，影響區(qū)覆蓋整個(gè)蘇南丘陵地區(qū)，此時(shí)蘇南中西部暴雨強(qiáng)度達(dá)最強(qiáng);17時(shí)前后渦旋開始明顯減弱，渦旋有一短時(shí)加速東移過程，影響范圍也東移，18時(shí)以后移速趨緩，渦旋區(qū)范圍達(dá)最小，降水強(qiáng)度明顯減弱;至20時(shí)(圖5c)，渦旋位于120.8～121.2°E、30.8～31.0°N，水平尺度只有45 km左右，垂直尺度不到6 km，形態(tài)比(H/L)約為0.133，700 hPa高度渦旋區(qū)內(nèi)最大水平氣壓梯度不到1.1 hPa/(10 km)，此時(shí)它對這場暴雨的影響也接近尾聲，蘇南地區(qū)的降水也漸漸停止;21時(shí)以后，它完全東移入海，對此次降水過程的影響徹底終止。由此可見:β中尺度渦旋對此次暴雨過程的加強(qiáng)與維系起著重要的作用，這一結(jié)果與史小康等(2007)的研究發(fā)現(xiàn)一致。

圖5 模擬的9月17日12時(shí)(a)、16時(shí)(b)和20時(shí)(c)850 hPa流場Fig.5 The simulated streamline at 850 hPa at(a)12:00 BST，(b)16:00 BST，and(c)20:00 BST 17 September

3.3 暴雨的水汽條件診斷

圖6為2009年9月17日08時(shí)相對濕度沿120°E的垂直剖面?？梢?，對流層中低層相對濕度高達(dá)80%～95%，暴雨區(qū)上空中低層的潮濕空氣為大暴雨的形成提供了較好的水汽條件。在整個(gè)暴雨過程中，江淮流域一直存在著水汽輸送帶。模擬結(jié)果顯示:2009年9月17日08時(shí)對流低層(850 hPa、925 hPa)有兩支明顯的東西向的水汽輸送通道存在，一支在暴雨區(qū)西南側(cè)，從西南向東北輻合進(jìn)入暴雨區(qū)，另一支從偏東向西呈扇形向暴雨區(qū)輻合(圖略);到17日14時(shí)，850 hPa水汽通量場中江蘇的沿江和蘇南地區(qū)(即暴雨區(qū))成為極顯著的水汽匯聚區(qū)(圖略)，之后一直源源不斷的有水汽向該區(qū)域輸送，為暴雨的增強(qiáng)與維持創(chuàng)造了條件。水汽通量散度場能更好地反映這一點(diǎn)，2009年9月17日08時(shí)850 hPa沿江和蘇南地區(qū)開始出現(xiàn)水汽輻合區(qū)，到了14時(shí)，水汽輻合達(dá)最強(qiáng)(圖7)。圖8為17日14時(shí)的水汽通量散度沿120°E的垂直剖面，可看出水汽輻合帶在500 hPa以下非常顯著，且隨高度增加水

3.4 暴雨的動(dòng)力條件診斷

汽輻合帶向北傾斜，在暴雨區(qū)形成了深厚的高濕環(huán)境，為對流的不穩(wěn)定增長以及暴雨的維持奠定了基礎(chǔ)(易軍等，2009)。

3.4.1 垂直速度

大氣中能量的輸送主要是通過空氣運(yùn)動(dòng)得以實(shí)現(xiàn)，垂直運(yùn)動(dòng)對水汽、熱量、動(dòng)量等物理量的輸送及對天氣系統(tǒng)的發(fā)展都起著極為重要的作用(Hoskins，1974)。由2009年9月17日14時(shí)垂直速度沿120°E的垂直剖面(圖9)可知，暴雨區(qū)低層到高層都出現(xiàn)了較為強(qiáng)烈的垂直運(yùn)動(dòng)。500～700 hPa是對流發(fā)展最旺盛的大氣層，其垂直運(yùn)動(dòng)的速度大小反映了系統(tǒng)對流運(yùn)動(dòng)狀況并影響高層的垂直運(yùn)動(dòng)。對模擬結(jié)果診斷分析發(fā)現(xiàn):17日14時(shí)前后，垂直運(yùn)動(dòng)發(fā)展迅速，500～700 hPa垂直運(yùn)動(dòng)的強(qiáng)度最強(qiáng)、范圍最廣。強(qiáng)烈的垂直運(yùn)動(dòng)使得大氣層結(jié)不穩(wěn)定能量快速釋放，致使對流天氣增強(qiáng)并得以維系。

3.4.2 渦度和散度

圖6 模擬的2009年9月17日08時(shí)相對濕度沿120°E的經(jīng)向垂直剖面(單位:%)Fig.6 Latitude-height cross-section of simulated relative humidity along 120°E at 08:00 BST 17 September 2009(units:%)

圖7 模擬的2009年9月17日14時(shí)850 hPa的水汽通量散度(單位:10－8g·cm－2·hPa－1·s－1)Fig.7 Simulated divergence of moisture flux at 850 hPa at 14:00 BST 17 September 2009(units:10－8g·cm－2·hPa－1·s－1)

圖8 模擬的2009年9月17日14時(shí)的水汽通量散度沿120°E的經(jīng)向垂直剖面(單位:10－8g·cm－2·hPa－1·s－1)Fig.8 Latitude-height cross-section of simulated divergence of moisture flux along 120°E at 14:00 BST 17 September 2009(units:10－8g·cm－2·hPa－1·s－1)

從高低空渦度、散度場的配置來看:暴雨開始前(17日02時(shí))，低層正渦度的強(qiáng)度較弱、范圍較小，中層500 hPa為弱輻散區(qū);08時(shí)，800～900 hPa的正渦度開始增大，500 hPa仍為弱輻散區(qū)，但高層200 hPa上輻散增強(qiáng)，低層900 hPa附近出現(xiàn)小范圍的輻合;之后的幾小時(shí)內(nèi)，低空正渦度隨著降水的增強(qiáng)而增強(qiáng)，輻合層逐漸加厚、輻合中心加強(qiáng);到了14時(shí)，在暴雨區(qū)上空900 hPa散度場出現(xiàn)高達(dá)－40×10－5s－1的強(qiáng)輻合中心(圖10)，200 hPa附近存在20×10－5s－1的強(qiáng)輻散中心，可見暴雨所在地具備中尺度環(huán)境下有利于對流發(fā)展的低層輻合、高層輻散的動(dòng)力配置結(jié)構(gòu)(陳紅專和湯劍平，2009)。同時(shí)，暴雨區(qū)上方850 hPa上正渦度顯著增強(qiáng)并達(dá)到最大值(圖11)，在200 hPa附近的負(fù)渦度中心仍高達(dá)－10×10－5s－1，中心位于120.3°E、31.4°N附近，正渦度區(qū)從900 hPa一直向高層延伸到350 hPa，在暴雨區(qū)形成了一個(gè)從下到上的正渦管，促使強(qiáng)對流發(fā)展。說明渦度場與散度場在垂直結(jié)構(gòu)配置上一致，對于暴雨的加強(qiáng)和維持有重要的促進(jìn)作用。

圖9 模擬的2009年9月17日14時(shí)垂直速度沿120°E的經(jīng)向垂直剖面(單位:cm·s－1)Fig.9 Latitude-height cross-section of simulated vertical velocity along 120°E at 14:00 BST 17 September 2009(units:cm·s－1)

圖10 2009年9月17日14時(shí)900 hPa的散度場(單位:10－5s－1)Fig.1 0Divergence field at 900 hPa at 14:00 BST 17 September 2009(units:10－5s－1)

3.5 暴雨的熱力條件診斷

圖11 2009年9月17日14時(shí)850 hPa渦度場(單位:10－5s－1)Fig.1 1Vorticity field at 850 hPa at 14:00 BST 17 September 2009(units:10－5s－1)

假相當(dāng)位溫θse是在大氣的干、濕絕熱過程中都守恒的一個(gè)重要特征參數(shù)，θse的垂直分布可以反映大氣的對流性不穩(wěn)定。當(dāng)氣層中θse隨高度減小時(shí)，整層空氣抬升后，氣層表現(xiàn)為整層位勢不穩(wěn)定;反之，當(dāng)θse隨高度增大時(shí)，整個(gè)氣層被抬升后，氣層能量將會變得更加穩(wěn)定(朱乾根等，2000)。由17日14時(shí)θse沿120°E的垂直剖面(圖12)可知，對流層中部500～700 hPa近似等于零，大氣為中性層結(jié);而低層＜0，為位勢不穩(wěn)定大氣，所以氣層內(nèi)有對流發(fā)展。注意到31.8～32.0°N附近，500 hPa以下θse幾乎不變，表明該地區(qū)中低層垂直方向上有強(qiáng)烈的上升運(yùn)動(dòng)，使低層水汽、能量向上混合輸送。

3.6 暴雨的配置場特征診斷

圖13模擬了暴雨天氣過程中動(dòng)力和水汽條件的配置情況。在31.8～32.2°N的暴雨區(qū)附近，1 000～400 hPa形成一條狹長的垂直運(yùn)動(dòng)帶，上升運(yùn)動(dòng)劇烈，中心速度達(dá)到60 cm·s－1(圖13a)。同時(shí)，在此區(qū)域內(nèi)，出現(xiàn)大范圍水汽通量輻合(圖13b)，可見強(qiáng)烈上升運(yùn)動(dòng)和充分水汽條件的合理配置對該區(qū)的暴雨發(fā)生至關(guān)重要(王巍巍等，2007)。

4 結(jié)論與討論

應(yīng)用中尺度數(shù)值預(yù)報(bào)模式WRF2.2，結(jié)合常規(guī)氣象觀測資料、AWMS系統(tǒng)監(jiān)測資料、NECP氣象再分析資料，模擬分析了2009年9月17日發(fā)生在江蘇大部(包括滬寧高速公路在內(nèi))的一次暴雨天氣過程。通過對模擬結(jié)果診斷分析，得到如下結(jié)論:

1)分析不同高度的水平流場發(fā)現(xiàn)，長江中下游地區(qū)存在一個(gè)β中尺度渦旋的發(fā)展、移動(dòng)過程，此β中尺度渦旋對暴雨過程降水的加強(qiáng)、維持起著重要的作用。

圖12 模擬的2009年9月17日14時(shí)假相當(dāng)位溫沿120°E的經(jīng)向垂直剖面(單位:K)Fig.1 2Latitude-height cross-section of simulated pseudoequivalent temperature along 120°E at 14:00 BST 17 September 2009(units:K)

2)暴雨區(qū)上空對流中低層空氣潮濕，為大暴雨的形成提供了較好的水汽條件。水汽輻合帶在500 hPa以下非常顯著，且隨高度增加水汽輻合帶向北傾斜，在暴雨區(qū)形成了深厚的高濕環(huán)境，為對流的不穩(wěn)定增長以及暴雨的維持奠定了基礎(chǔ)。

圖13 模擬的2009年9月17日14時(shí)垂直速度(a;cm·s－1)和水汽通量散度(b;10－8g·cm－2·hPa－1·s－1)沿120°E的經(jīng)向垂直剖面Fig.1 3Latitude-height cross-sections of simulated(a)vertical velocity(cm·s－1)and(b)divergence of moisture flux(10－8g·cm－2·hPa－1·s－1)along 120°E at 14:00 BST 17 September 2009

3)雨區(qū)低層到高層都出現(xiàn)了較為強(qiáng)烈的垂直運(yùn)動(dòng)，使得大氣層結(jié)不穩(wěn)定能量快速釋放，致使對流天氣增強(qiáng)并得以維系。暴雨區(qū)低層輻合、高層輻散，低空為正渦度、高空為負(fù)渦度，且前期隨著降水的增強(qiáng)，低空正渦度逐漸增強(qiáng)，輻合層逐漸加厚、輻合中心加強(qiáng)。渦度場與散度場在垂直結(jié)構(gòu)配置上一致，對暴雨的加強(qiáng)和維持有重要的促進(jìn)作用。

4)暴雨區(qū)對流層中低層有強(qiáng)烈的垂直上升運(yùn)動(dòng)，促使低層水汽、能量等向上混合輸送。對流層中層大氣為中性層結(jié)，而低層大氣為位勢不穩(wěn)定，所以整層大氣有對流發(fā)展，有利于暴雨的形成。

通過個(gè)例研究發(fā)現(xiàn):水汽條件、動(dòng)力條件和熱力條件是暴雨形成的重要物理基礎(chǔ)，調(diào)整WRF模式中的相關(guān)參數(shù)可以成功地模擬高速公路上的暴雨天氣過程，結(jié)合實(shí)測數(shù)據(jù)可以較好地預(yù)測高速公路暴雨災(zāi)害產(chǎn)生的時(shí)間、范圍和強(qiáng)度等，為交通氣象部門準(zhǔn)確地預(yù)測道路災(zāi)害性天氣提供科學(xué)依據(jù)。模擬的降雨分布與實(shí)況基本一致，因此，可以認(rèn)為對此次天氣過程的模擬是較成功的。當(dāng)然，模擬也與實(shí)測之間存在一些誤差，這是由于WRF模式是針對中尺度天氣系統(tǒng)設(shè)計(jì)的，再者，NCEP資料的時(shí)空分辨率較粗，也導(dǎo)致了模式輸出與實(shí)測之間有一定的誤差。對于模式數(shù)值試驗(yàn)中存在的誤差，未來將從下列兩個(gè)方面力求改進(jìn):一是對模式的初始場資料進(jìn)行加密、同化和調(diào)整以提高預(yù)報(bào)準(zhǔn)確率;二是對模式物理過程和參數(shù)化方案不斷進(jìn)行試驗(yàn)、選擇，找出最佳配置以便更精確、更科學(xué)、更合理地反映實(shí)際發(fā)生的狀況。

陳紅專，湯劍平.2009.一次突發(fā)性特大暴雨的中尺度分析和診斷［J］.氣象科學(xué)，29(6):797-803.

陳業(yè)國，農(nóng)孟松.2010.2008年6月廣西持續(xù)性暴雨的診斷與數(shù)值模擬［J］.氣象科學(xué)，30(2):250-255.

陸漢城，楊國祥.2000.中尺度天氣原理和預(yù)報(bào)［M］.北京:氣象出版社.

沈桐立，曾瑾瑜，朱偉軍，等.2010.2006年6月6—7日福建特大暴雨數(shù)值模擬和診斷分析［J］.大氣科學(xué)學(xué)報(bào)，33(1):14-24.

史薈燕，潘曉濱，謝海軍，等.2009.一次華南暴雨過程中不同區(qū)域降水特征的模擬分析［J］.氣象科學(xué)，29(1):58-63.

史小康，李耀東，高守亭，等.2007.一次江淮暴雨的數(shù)值模擬及暴雨落區(qū)的診斷分析［J］.氣象科學(xué)，27(1):26-34.

孫建華，張小玲，齊琳琳，等.2004.2002年中國暴雨試驗(yàn)期間一次低渦切變上發(fā)生發(fā)展的中尺度對流系統(tǒng)研究［J］.大氣科學(xué)，28(5):675-690.

王建捷，李澤椿.2002.1998年一次梅雨鋒暴雨中尺度對流系統(tǒng)的模擬與診斷分析［J］.氣象學(xué)報(bào)，60(2):146-151.

王巍巍，張艷玲，壽紹文.2007.一次局地大暴雨過程的數(shù)值模擬和診斷分析［J］.南京氣象學(xué)院學(xué)報(bào)，30(3):411-416.

王智，高坤，翟國慶.2003.一次與西南低渦相聯(lián)系的低空急流的數(shù)值研究［J］.大氣科學(xué)，27(1):75-85.

吳贊平，王宏偉，袁成松，等.2006.滬寧高速公路氣象決策支援系統(tǒng)［J］.現(xiàn)代交通技術(shù)(5):91-94.

易軍，壽紹文，高守亭，等.2009.一次梅雨鋒暴雨過程的中尺度特征分析［J］.氣象科學(xué)，29(6):761-767.

章國材.2006.美國WRF模式的進(jìn)展和應(yīng)用前景［J］.氣象，30(12):27-31.

張曉紅，賈天山，胡雯.2009.2007年淮河流域雨季暴雨多尺度環(huán)流分析［J］.南京氣象學(xué)院學(xué)報(bào)，32(2):321-326.

張小玲，陶詩言，張慶云.2002.1998年7月20—21日武漢地區(qū)梅雨鋒上突發(fā)性中β系統(tǒng)的發(fā)生發(fā)展分析［J］.應(yīng)用氣象學(xué)報(bào)，13(4):385-397.

張羽，牛生杰，于華英，等.2010.雷州半島“07.8”致洪特大暴雨的數(shù)值模擬［J］.大氣科學(xué)學(xué)報(bào)，33(1):47-57.

朱乾根，林錦瑞，壽紹文，等.2000.天氣學(xué)原理和方法［M］.北京:氣象出版社.

Cui X，Gao S，Li X.2007.Diagnostic analysis of mesoscale rainstorms in the Jiang-Huai valley of China with convective vorticity vector［J］.Progress in Natural Science，17(5):559-568.

Hoskins B J.1974.The role of potential vorticity in symmetric stability and instability［J］.Quart J Roy Meteor Soc，100(425):480-482.

Joseph B K.2004.Weather research and forecasting model:A technical overview［C］//84th AMS Annual Meeting，Seattle，U.S.A.，A-merican Meteorological Society，January 10—15，2004.

Molinari J，Dudek M.1992.Parameterization of convective precipitation in mesoscale numerical models:A critical review［J］.Mon Wea Rev，120:326-344.

Paula T M，John R S，Brent S.2005.User’s guide for advanced research WRF(ARW)modeling system version 2.1［M］.Boulder，Colorado:Weather News，Norman Oklahoma.

Janjic Z I.2004.The NCEP WRF core［C］//84th AMS Annual Meeting，Seattle，U.S.A，American Meteorological Society，January 10—15，2004.