王碩甫， 鐘蘭頔， 李瑞， 李彩玲

(1．佛山市氣象局，廣東佛山　528000；2.廣東省氣象臺，廣東廣州　510080；3.南昌市氣象局，江西南昌　330000；4.中國民航三亞空管站，海南三亞　572000)

?

臺風“尤特”(1311)暴雨的數(shù)值模擬及螺旋度診斷分析

王碩甫1，2， 鐘蘭頔3， 李瑞4， 李彩玲1

(1．佛山市氣象局，廣東佛山528000；2.廣東省氣象臺，廣東廣州510080；3.南昌市氣象局，江西南昌330000；4.中國民航三亞空管站，海南三亞572000)

摘要：為了研究臺風暴雨的環(huán)流特征，利用NECP/NCAR 1°× 1°再分析資料，應用中尺度WRF模式對“尤特”進行數(shù)值模擬。模擬結(jié)果表明，副高的位置、西南季風急流的水汽輸送和大尺度輻合是造成華南地區(qū)大范圍暴雨的重要原因。利用模擬結(jié)果，對臺風暴雨的垂直螺旋度特征的進一步研究后指出，垂直螺旋度能很好地反映低層大氣輻合、高層大氣輻散的垂直結(jié)構(gòu)，螺旋度大值中心的位置與臺風暴雨落區(qū)有較好的對應關(guān)系，螺旋度迅速增大預示著降水急劇增加，對臺風暴雨的強度和落區(qū)預報有較好的指示作用。

關(guān)鍵詞：天氣學； 臺風暴雨； 數(shù)值模擬； 螺旋度； “尤特”(1311)

近年來，隨著多普勒雷達、氣象衛(wèi)星、中尺度自動站等觀測手段應用于臺風的高分辨率監(jiān)測，對臺風的發(fā)生發(fā)展條件、結(jié)構(gòu)、移動路徑、風雨特征以及臺風暴雨影響因子等有了眾多的研究結(jié)果[1-7]。此后由于數(shù)值模式具有較高的時空分辨率和能較好地反演天氣過程等優(yōu)勢，對臺風過程進行數(shù)值模擬試驗逐漸成為了研究臺風熱力結(jié)構(gòu)、演變過程和臺風渦旋的主要手段[8-9]。周冠博等[10]模擬了臺風“鳳凰”的登陸過程，其診斷結(jié)果表明散度垂直通量能夠描述低層輻合、高層輻散的垂直動力結(jié)構(gòu)，能較好地預報臺風暴雨的落區(qū)；閔穎等[11]對臺風“珍珠”的模擬和診斷結(jié)果表明，臺風暴雨和大風主要集中在螺旋云帶范圍內(nèi)，螺旋云帶內(nèi)的強抽吸作用和毛細效應是上升運動強烈發(fā)展的主要因素。

螺旋度是用來表示流體沿旋轉(zhuǎn)方向運動強弱的物理量，是有利于維持大氣對流運動的重要基本特征之一，被廣泛應用于暴雨、雷暴、臺風、龍卷等強對流天氣的診斷分析[12-14]。許多相關(guān)研究表明，螺旋度中心位置、強度的變化對臺風暴雨具有預報意義[15-17]。本研究在進行大氣環(huán)流分析和數(shù)值模擬的基礎上，利用物理量診斷的方法，分析了2013年第11號臺風“尤特”登陸過程中強降水的分布和螺旋度與強降水的對應關(guān)系，希望能夠?qū)σ院笈_風暴雨的預報提供一定的幫助。

1“尤特”臺風簡介

1311號臺風“尤特”于2013年8月9日18：00(世界時，下同)在菲律賓以東洋面生成后向西北方向移動，強度迅速加強，于11日09：00達到超強臺風級別，11日19：00在菲律賓呂宋島東部沿海登陸，之后繼續(xù)向西北方向移動，以臺風強度進入南海。在南海，“尤特”于13日00：00再次加強為強臺風，并于14日07：50前后登陸廣東省陽西縣，登陸時中心氣壓955 hPa，最大風速42 m/s。

本研究利用NECP/NCAR 1°×1°再分析資料和Micaps降水資料分析大氣環(huán)流背景和降水情況，應用WRF模式對2013年8月13—16日的“尤特”臺風登陸過程開展數(shù)值模擬并進行暴雨過程的螺旋度診斷分析。

2數(shù)值模擬方案設計

本研究利用非靜力中尺度數(shù)值模式WRF3.4.1版本對強臺風“尤特”的登陸過程進行模擬，選取NECP/NCAR 1°×1°再分析資料作為模式初始場，采用雙重嵌套網(wǎng)格，模式中心為22.54°N，115°E，模擬方案設計如表1所示。積分時間為2013年8月13日00：00——16日00：00共72 h，積分步長為180 s。

表1　WRF 模式參數(shù)設定

3模擬結(jié)果及驗證

3.1環(huán)流形勢

由圖1可以看出，在整個積分時間段內(nèi)，模擬的500 hPa副高、臺風閉合環(huán)流的形狀和位置的變化特征與實況一致。850 hPa的模擬流場與實況十分相似，流場中心位置基本吻合，而且也描繪出了由13日00：00流場中心東西兩側(cè)風速不對稱逐漸轉(zhuǎn)為14日00：00的對稱性結(jié)構(gòu)的特征。但不足之處在于，在臺風登陸后，模擬風場陸地上的風速較實況偏大，其原因可能與模式下墊面對風速的影響考慮不足有關(guān)。

圖1　實況(a、b、c)與模擬(d、e、f)的500 hPa高度場(等值線，單位：dagpm)及850 hPa急流(風向桿，單位：m/s)對比

從大尺度環(huán)流特征來看，隨著臺風環(huán)流的西行，在500 hPa上其逐漸轉(zhuǎn)為偏南氣流引導，移動路徑由西北路徑轉(zhuǎn)為偏北路徑。到了15日00：00，中緯度地區(qū)副高西伸，同時菲律賓上空存在塊狀高壓，“尤特”同時受到東風氣流與西風氣流作用，引導氣流不明顯，臺風環(huán)流維持在兩廣交界處。從850 hPa風場的情況可以看到，從13日00：00起，一條寬廣的西南低空急流帶與臺風環(huán)流連接，為臺風環(huán)流的發(fā)展加強提供了豐富的能量和水汽。臺風登陸后，從孟加拉灣經(jīng)中南半島到華南地區(qū)仍維持著急流強度的西南風到偏南風，臺風環(huán)流中心的東側(cè)明顯存在偏南風和偏東風的輻合，華南地區(qū)位于大尺度水汽和風向輻合的中心。

3.2臺風路徑

圖2是“尤特”臺風實況與模擬路徑對比。

圖2　臺風“尤特”的實況(黑線)與模擬(紅線)路徑

可以看出，從模擬開始到15日00：00，模擬與實況的臺風路徑較為一致，模擬試驗中臺風的移速稍慢于實況。實況與模擬路徑都顯示臺風在登陸前后有西折的現(xiàn)象，但是實況路徑的西折發(fā)生在14日06：00到12：00，模擬路徑發(fā)生在14日12：00到18：00。在15日00：00之后臺風路徑開始向東折，但是模擬路徑的東移不是很明顯，速度也明顯慢于實況。

3.3雷達回波與24 h降水

對比2013年8月14日“尤特”登陸期間實況及模擬的雷達回波分布情況可以看到(圖略)，模擬的云團回波與實況的分布形式比較一致，臺風的眼區(qū)、云墻和螺旋雨帶等結(jié)構(gòu)的回波特征十分清晰，也描繪出了臺風中心西南側(cè)的回波明顯比東北側(cè)強度強、組織性高的特征。模擬回波的上岸時間、地點與實況基本吻合，模擬試驗還反映出了“尤特”登陸后回波迅速減弱的特點。不同之處在于，在模擬圖中可以看出00：00和06：00臺風中心東北側(cè)的回波比實況偏強，并且12：00的模擬雷達回波上臺風眼仍清晰可見，但實況雷達回波已難分辨出臺風眼區(qū)。

圖3給出了8月13日00：00—16日00：00的逐24 h實況及模擬降水分布。在“尤特”登陸前，廣東東南部的蓮花山脈南坡及海南東北部地區(qū)已出現(xiàn)暴雨級別降水，而上述地區(qū)的模擬降水量級偏小?！坝忍亍钡顷懏斕?，廣東西南部和廣東東南部沿海地區(qū)出現(xiàn)100 mm以上降水，數(shù)值模擬對廣東西南部的降水落區(qū)吻合較好，但降水中心量級偏大，對于廣東東南部的降水則量級與范圍都明顯偏小。15日，大暴雨區(qū)繼續(xù)向北延伸，珠三角地區(qū)、廣東西南部沿海及廣東西北部、廣西東北部均出現(xiàn)100 mm以上降水，50 mm以上降水范圍也有明顯擴大，模擬降水能很好地反映出當日降水范圍擴大的特征，但存在降水中心位置偏西偏南以及中心量級明顯偏大的問題，對粵東地區(qū)的降水減弱也沒有很好地反映出來。

圖3　13日00：00—14日00：00(a，b)、14日00：00—15日00：00(c，d)、15日00：00—16日00：00(e，f)24 h實況(a，c，e)和模擬(b，d，f)累積降水量(單位：mm)對比

4螺旋度診斷分析

本研究利用模式輸出的高時空分辨率資料，選取垂直螺旋度進行強降水的動力診斷，Z坐標系下垂直螺旋度的表達式為

在對各層垂直螺旋度和散度隨時間的變化進行對比分析后發(fā)現(xiàn)，在13日，高低層大氣的輻合輻散不明顯，此時該區(qū)域內(nèi)的垂直上升運動并不強烈，整層大氣的垂直螺旋度基本為0。14日，大氣垂直方向上形成了低層輻合、高層輻散的散度場配置，在強“抽吸作用”下引起強烈的上升運動，整層大氣的垂直螺旋度迅速增大。隨后低層負散度中心和中高層正散度中心強度有所衰減，正垂直螺旋度明顯減小，說明低層輻合、高層輻散引起的上升運動已明顯減弱。從15日09：00開始，區(qū)域內(nèi)再次形成低層輻合、高層輻散的高低空配置，正垂直螺旋度重新加大。因此可以說明，正垂直螺旋度的變化與大氣的輻合輻散運動有關(guān)，垂直螺旋度能夠反映大氣低層輻合、高層輻散的垂直動力結(jié)構(gòu)，且在臺風暴雨過程中850和300 hPa的垂直螺旋度變化尤為明顯。

850 hPa垂直螺旋度和隨后6 h的降水分布如圖4所示。14日06：00，強降水區(qū)位于廣東西部海面上和廣東陽江的沿海地區(qū)，對應的垂直螺旋度在廣東西部海面存在大值區(qū)，但位置比強降水區(qū)略偏南偏東，范圍也偏小。14日12：00，強降水區(qū)位于廣東西部地區(qū)，廣東西部的垂直螺旋度大值區(qū)的位置與強降水中心基本一致。14日18：00，強降水區(qū)移至粵桂交界地區(qū)，垂直螺旋度在粵桂交界的廣東一側(cè)存在大值區(qū)，位置與強度與廣東的強降水區(qū)對應較好，在廣西境內(nèi)也有大值區(qū)，但位置比廣西的強降水區(qū)明顯偏東。15日00：00，強降水區(qū)繼續(xù)西移至廣西東南部地區(qū)，此時在粵桂交界的廣西一側(cè)有垂直螺旋度大值區(qū)，位置比強降水區(qū)略偏東。300 hPa垂直螺旋度的分布特征與850 hPa垂直螺旋度基本一致(圖略)。總的來說，在強降水區(qū)范圍內(nèi)能找到對應的垂直螺旋度大值區(qū)，垂直螺旋度大值區(qū)內(nèi)的降水比弱螺旋度區(qū)域的降水大，垂直螺旋度大值區(qū)的強度與臺風強降水的強度有較好的對應關(guān)系。

圖4　850 hPa螺旋度(等值線，單位：×10-5m·s-2)(a，c)和滯后6 h降水(填色，單位：mm)(b，d)

圖5顯示了13日00：00—16日00：00暴雨區(qū)300和850 hPa平均垂直螺旋度與3 h累積降水量隨時間的變化?？梢钥吹剑怪甭菪鹊牟ㄐ闻c降水量隨時間的變化十分吻合。13日00：00到14日00：00，高低2層均表現(xiàn)為零或弱垂直螺旋度，此時并未發(fā)生強降水。14日03：00起，300和850 hPa的垂直螺旋度迅速增大，到了15：00 300 hPa正垂直螺旋度超過14×10-5m·s-2，850 hPa正垂直螺旋度超過10×10-5m·s-2，對比降水量發(fā)現(xiàn)，該區(qū)域15：00到18：00的3 h累計降水量達到最大值。隨后垂直螺旋度和降水量均逐漸減小。從15日09：00開始，300 hPa正垂直螺旋度重新加大，最高值為4×10-5m·s-2，850 hPa正垂直螺旋度也于15日18：00加大到3×10-5m·s-2，此時區(qū)域內(nèi)的降水量再次加大，但總體強度小于14日06：00到15日00：00內(nèi)的降水。15日18：00后2層正垂直螺旋度減小的同時，強降水也明顯減弱。這進一步表明了，在時間演變趨勢上，垂直螺旋度與未來3 h的降水強度存在對應關(guān)系，且垂直螺旋度強度的變化提前于降水強度的變化。

圖513—15日暴雨區(qū)300和850 hPa平均垂直螺旋度與3 h降水隨時間變化

5結(jié)論

1)WRF模式能較為成功地模擬出“尤特”臺風登陸過程和降水分布。模擬的500 hPa副高閉合環(huán)流和850 hPa西南急流的位置、形狀、強度等特征與實況較為一致，模擬850 hPa陸上風速較實況偏大。在臺風移向移速方面，模擬的臺風路徑與觀測路徑較為一致，但在移動速度和登陸后的東轉(zhuǎn)上，二者存在一定的差別。在對比模擬與實況的雷達回波圖后發(fā)現(xiàn)，模式能較好反映出回波的結(jié)構(gòu)特征和位置，但在強度上比實況偏強。此外，模擬還反映出了尤特登陸后回波迅速減弱的特點。同時，對臺風中心附近降水分布的模擬與實況降水比較一致，能較好地模擬出14日和15日的強降水落區(qū)，但降水中心量級偏大；另外，模式對于廣東東南部地區(qū)的臺風外圍降水的模擬效果較差。

2)通過研究“尤特”臺風登陸期間的環(huán)流情況，可知臺風登陸后副高的引導氣流不明顯、臺風環(huán)流穩(wěn)定少動、低層西南低空急流帶與臺風環(huán)流持續(xù)連接以及臺風環(huán)流中心東側(cè)存在強風向、水汽輻合區(qū)是造成華南地區(qū)大范圍暴雨的重要原因。

3)利用模式輸出的高時空分辨率資料，對“尤特”的強降水進行了垂直螺旋度的診斷分析。結(jié)果表明，正垂直螺旋度的變化與大氣的輻合輻散運動有關(guān)，且在臺風暴雨過程中850和300 hPa的垂直螺旋度變化特征尤為明顯；在空間分布上，垂直螺旋度中心的位置與未來6 h臺風中心附近的強降水區(qū)有較好的對應關(guān)系；在時間演變趨勢上，850和300 hPa的垂直螺旋度極大值出現(xiàn)的時間與未來3 h短時強降水發(fā)生時間一致，垂直螺旋度大小與降水強度存在較好的對應關(guān)系。這說明垂直螺旋度能反映臺風強降水的發(fā)展狀況和劇烈程度，對預報強降水有較好的指示意義。

參考文獻:

[1]許映龍，矯梅燕，畢寶貴，等.近海臺風雷達定位方法應用研究[J].大氣科學，2006，30(6)：1119-1128.

[2]李霞，林少冰，王捷純，等.利用衛(wèi)星資料分析臺風“浪卡”非對稱結(jié)構(gòu)及降水[J].廣東氣象，2011，33(3)：1-6.

[3]康志明，陳濤，錢傳海，等.0604號強熱帶風暴“碧利斯”特大暴雨的診斷研究[J].高原氣象，2008，27(3)：596-607.

[4]鄧文劍，吳振鵬，吳乃庚，等.臺風“巨爵”近海急劇加強的特征及成因分析[J].廣東氣象，2014，36(3)：9-14.

[5]黃曉瑩，程正泉.臺風“凡亞比”造成馬貴鎮(zhèn)致洪暴雨分析[J].廣東氣象，2011，33(6)：1-3.

[6]程正泉.登陸臺風與環(huán)境因子相互作用對暴雨的影響研究綜述[J].廣東氣象，2008，30(5)：4-7.

[7]程正泉，陳聯(lián)壽，李英.登陸臺風降水的大尺度環(huán)流診斷分析[J].氣象學報，2009，67(5)：840-850.

[8]程銳，葉成志，許愛華，等.臺風“云娜”的熱動力結(jié)構(gòu)模擬試驗研究[J].暴雨災害，2007，26(2)：103-108.

[9]沈新勇，明杰，方珂.臺風渦旋系統(tǒng)的波動性質(zhì)及其數(shù)值模擬[J].氣象科學，2007，27(2)：176-186.

[10]周冠博，崔曉鵬，高守亭.臺風“鳳凰”登陸過程的高分辨率數(shù)值模擬及其降水的診斷分析[J].大氣科學，2012，36(1)：23-34.

[11]閔穎，沈桐立，朱偉軍，等.臺風“珍珠”螺旋雨帶的數(shù)值模擬與診斷分析[J].大氣科學學報，2010，33(2)：227-235.

[12]陸慧娟，高守亭.螺旋度及螺旋度方程的討論[J].氣象學報，2003，61(6)：684-691.

[13]陳華，談哲敏.熱帶氣旋的螺旋度特性[J].熱帶氣象學報，1999，15(1)：81-85.

[14]段朝霞，黃莉莉.“09.5”廣東特大暴雨過程的螺旋度特征[J].廣東氣象，2011，33(3)：16-19.

[15]劉漢華，唐偉民，趙利剛.2008年“鳳凰”臺風暴雨的水汽和螺旋度分析[J].氣象科學，2010，30(3)：344-350.

[16]余貞壽，倪東鴻，閔錦忠.超強臺風“圣帕”(0709)特大暴雨過程的完全螺旋度分析[J].南京氣象學院學報，2009，32(1)：45-53.

[17] 曹楚，王忠東，林念萍，等.2013年“菲特”臺風暴雨的水汽和螺旋度分析[J].廣東氣象，2014，36(5)：41-45.

收稿日期：2015-09-15

基金項目：廣東省氣象局科學技術(shù)研究項目(2014B21)

作者簡介：王碩甫(1990年生)，男，本科，主要從事短期和短時天氣預報。E-Mail：huoyue_anthony@163.com

中圖分類號：P44

文獻標識碼：A

doi:10.3969/j.issn.1007-6190.2016.02.002

Numerical Simulation of Heavy Rain Associated withTyphoon Utor (1311) and Diagnostic Analysis of Its Helicity

WANGShuo-fu1, 2,ZHONGLan-di3,LIRui4,LICai-ling1

(1. Meteorological Bureau of Foshan City, Foshan 528000; 2. Meteorological Observatory of Guangdong Province,Guangzhou 510080; 3. School of Atmospheric Sciences, Chengdu University of Information Technology, ChengduSichuan 610225; 4. Air Traffic Regulator of Sanya, Civil Aviation Administration of China, Sanya Hainan 572000)

Abstract：For the purpose of studying the circulation characteristics of typhoon-related heavy rain, we numerically simulated Typhoon Utor using the NCEP/NCAR 1°×1° reanalysis and the mesoscale WRF model. As shown in the simulation, the location of the subtropical high, the transport of water vapor through the Southwest Monsoon jets and the convergence of large-scale convergence are important contributors to the extensive heavy rain in the south of China. It is discovered, after further analysis of the vertical helicity of typhoon-related heavy rain using the simulation, that the vertical helicity well reflects a vertical structure of low-level convergence versus upper-level divergence; the location of the center of the high helicity values is corresponding well with the area of the typhoon-related heavy rain; and rapid increase of the helicity indicates fast increase of rainfall, making it a good indicator of the intensity and rain area.

Key words：synoptics; typhoon-related heavy rain; numerical simulation; helicity; Utor (1311)

王碩甫， 鐘蘭頔， 李瑞， 等.臺風“尤特”(1311)暴雨的數(shù)值模擬及螺旋度診斷分析[J].廣東氣象，2015,38(2)：6-10.