廣西沿海地形對超強臺風“威馬遜”影響的數值試驗

蘇玉婷，林開平，肖志祥，黃小燕

（1.廣西師范學院，南寧 530022；2.廣西區氣象臺，南寧 530022；3.廣西區氣象減災研究所，南寧 530022）

引 言

臺風是影響我國的主要災害性天氣系統之一，其帶來的大風和強降水給途經地區帶來嚴重的經濟損失和人員傷亡。臺風強度突變（陳聯壽等，2002）和移動路徑突變（倪鐘萍等，2013）是臺風研究和業務預報的難點之一。2014年第9號臺風 “威馬遜”是1973年以來登陸華南的最強臺風，具有移動快、多次登陸、發展迅速等特點，給海南、廣東、廣西造成的經濟損失超過265億元，其登陸減弱為熱帶低壓后帶來的強降水造成了云南37人死亡，9人失蹤的重大災害（徐舒揚等，2016），“威馬遜”在沿海突然加強為超強臺風，給預報帶來了極大的困難。

有關研究表明，90%以上的臺風在近海受到摩擦耗損能量的影響其強度將減弱 （馮錦全和陳多，1995），作為有記錄以來登陸我國的最強臺風“威馬遜”，其近海突然增強的特點得到廣泛關注，大量研究取得了很多有意義的成果。程正泉等（2017）認為臺風下游熱力不穩定在暖海溫、低空急流、越赤道氣流、垂直風切變小等有利條件配合下，“威馬遜”得以在近海快速增強。鄭艷等（2014）指出低層輻合與高層輻散、弱風切變、適宜的海溫及深厚的暖渦是“威馬遜”在近海加強的主要原因。

除了大氣條件、海溫狀況等動力和熱力作用，地形也是影響臺風降水強度和落區的重要因素。數值試驗表明，不同地形高度下臺風路徑及降水空間分布差異很大，地形對暴雨增幅有重要作用 （朱紅芳等，2015）。地形的阻擋作用有利于背風坡降水減弱。臺風“納沙”影響海南期間，五指山地形的阻擋使得海南島東南部降水減弱50-150mm（楊仁勇等，2014）。“威馬遜”在廣西沿海登陸時橫穿十萬大山山脈，在廣西引起大范圍強降水過程，本文嘗試使用中尺度模式WRF開展地形敏感性數值試驗，探討廣西沿海地形對“威馬遜”路徑和降水的影響。

1 試驗設計

超強臺風“威馬遜”7月18日05時在南海加強為超強臺風，19日07時在廣西防城港市再次登陸，隨后強度迅速減弱。“威馬遜”登陸廣西期間，廣西沿海各地出現了暴雨到大暴雨，局部特大暴雨，其中潿洲島19日的最大降水量達到了303.6mm。

廣西沿海地區有著名的十萬大山山脈，“威馬遜”橫穿十萬大山，為考慮十萬大山山脈對“威馬遜”的影響，本文針對十萬大山山脈設計了控制和敏感性試驗兩個方案共三組試驗，控制試驗是從7月18日06時模擬到7月19日18時（世界時），一共積分36小時，使用ERA-Interim再分析資料 （Dee et al.，2011）為初始邊界條件，水平分辨率 0.75°×0.75°，邊界條件每隔6h更新一次。WRF3.4.1使用的主要參數化方案如表1所示。模式為兩層嵌套，外層格點數為360×300，水平分辨率為9公里；內層格點數361×253，水平分辨率為3公里。內層嵌套區域覆蓋了“威馬遜”加強為超強臺風至消亡減弱的整個階段。兩組敏感性試驗方案中，地形被人為移除，其它參數與控制試驗相同。

表1 模式使用的主要參數化方案

具體第一組的控制試驗（圖1a，見彩頁），其結果作為參照，可以看到十萬大山山脈為西南—東北走向，與“威馬遜”移動路徑幾乎垂直，有兩個達700多m的狹長頂峰。地形敏感性的第一組為去掉了狹長的十萬大山山脈試驗（圖1b，見彩頁），第二組則只是移除了南部一些較小山脈的試驗（圖1c，見彩頁）。通過對比分析Exp1、Exp2試驗與控制試驗的結果，可以得到沿海山脈對“威馬遜”移動路徑及降水的影響狀況。

2 試驗結果分析

2.1 地形對臺風移動路徑的影響

圖2為三組數值試驗模擬的臺風移動路徑情況。分析可知“威馬遜”加強為超強臺風后移動路徑較為穩定，一直向西北方向移動，先后在雷州半島和廣西沿海登陸。三組試驗中，模式對臺風初始位置的定位并不準確，且模式模擬的臺風強度遠遠低于觀測（圖略），這可能與初始場的誤差有關，采用WRF模式自帶的 Bogus 技術（Sherrie et al.，2009），對初始場臺風強度和位置進行重定位，能夠有效改善臺風強度和位置模擬。此外，前18小時為模式spin-up階段，臺風路徑擺動較大。控制試驗的臺風登陸點與實況較為接近，總體來說模式能夠較好把握 “威馬遜”的登陸地點及移動方向。

圖2 WRF模擬的臺風“威馬遜”移動路徑觀測：黑色；參考試驗：紅色；地形試驗1：藍色；地形試驗2：綠色

進一步分析圖2可知，去除十萬大山山脈后，模式模擬的臺風登陸地點往東偏移較大，使得登陸后的臺風中心也隨之偏東偏南。而如果只是去除十萬大山南部的一個頂峰（第二組敏感性試驗），臺風登陸前的路徑與控制試驗則較為接近，說明在離臺風較遠時小山脈的影響較小；而在登陸后其移動路徑基本與移除整個十萬大山山脈試驗吻合。

2.2 地形對臺風降水的影響

7月19日“威馬遜”進入北部灣并在防城港市登陸，在廣西沿海及西南地區引發強降水。控制試驗雖然對沿海三市及崇左市南部的暴雨量級模擬偏小，百色南部的暴雨范圍偏大，但是對于南寧市北部的一些單站暴雨卻能夠模擬出來（圖3a，見彩頁），總體而言，模式能夠較好模擬出19日“威馬遜”在廣西造成的暴雨過程。從移除十萬大山地形的兩組敏感性數值試驗可以看到，防城港市和崇左市南部的降水都有所減弱，24小時降水的最大減少量達40mm以上（圖3d-e，見彩頁），兩者對比分析可知，移除整個十萬大山地形的試驗在降水減弱幅度和范圍上均大于只去掉南部山脈的試驗結果。越南范圍的降水在兩組敏感性數值試驗減小的幅度要大于廣西范圍。

從19日水汽輸送及水汽通量散度可以看到，十萬大山附近水汽為往北和往南輸送的輻散區 （圖4a，見彩頁），而Exp2試驗（圖 4b，見彩頁）水汽往北輸送，但是有風向的弱輻散，水汽輻散較Exp1試驗小。而在崇左市的氣旋式異常和南寧市南部的反氣旋式異常及北海市附近的氣旋式異常主要是由于敏感性試驗與控制試驗臺風中心位置（圖2，見彩頁）差異導致。在敏感性試驗中臺風在北海市附近登陸使該處成氣旋式異常，而臺風登陸后移動路徑偏西偏南使得崇左市為氣旋式異常而南寧市南部為反氣旋式異常（圖4，見彩頁）。兩組敏感性數值試驗表明移除沿海山脈后，由于臺風位置的改變，使得環流改變，水汽輻合強度減弱是導致山脈附近降水減少的原因之一。

數值試驗模擬的10m風場表明（圖5，見彩頁），當十萬大山存在時，山脈南部的偏南風較強，有利于低層水汽輸送，在迎風坡有較強的上升運動；而在山脈北部則為較弱的偏北風，也有利于地形抬升；南北風向在十萬大山南北部均受地形抬升作用，有利于強降水的形成（圖5a，見彩頁）。當十萬大山被移除（圖5b，見彩頁）時，輻合區在南風作用下北抬，且臺風中心位置偏西，使得氣旋式環流中心偏北，山脈附近的南北風輻合大大減弱，使得降水減弱。如果只移除十萬大山南部的小山脈，氣旋式環流中心仍然有北抬（圖5c，見彩頁），但是要較Exp1試驗偏南，十萬大山西側仍有較弱的北風受地形抬升的影響，因此其降水減弱幅度較小（圖3e，見彩頁）。

東經108°正好位于十萬大山山脈中部及臺風中心，從垂直剖面圖可以看到在控制試驗中，22°—22.5°N中低層有垂直上升運動，即十萬大山山脈南部和北部受南風及北風的影響，在地形抬升作用下形成的中低層上升運動（圖6a，見彩頁）。而在Exp1試驗中，由于十萬大山山脈被移除，沒有地形阻擋作用偏南風往北推，垂直運動減弱（圖6b，見彩頁），Exp2試驗中仍然可以看到十萬大山山脈低層的弱上升運動，說明十萬大山山脈的存在對低層風的地形抬升作用不可忽視，但是垂直上升運動的厚度明顯減弱（圖 6c，見彩頁）。

3 小結

超強臺風“威馬遜”在近海快速增強后，先后在雷州半島及廣西沿海登陸，隨后以近乎垂直的角度橫穿十萬大山山脈，給廣西沿海及桂西南地區帶來強降水天氣。本文使用區域數值模式WRF3.4.1針對十萬大山山脈開展控制和敏感性數值試驗，結果表明：

（1）WRF模式能夠較好地模擬出“威馬遜”的移動路徑、登陸地點及降水空間分布特征。十萬大山山脈以北的偏北風及以南的偏南風受地形抬升作用，在山脈南北兩側都形成了中低層的垂直上升運動，有利于強降水的出現。

（2）十萬大山山脈移除后，“威馬遜”登陸位置偏東，移動路徑隨之偏南偏西。此外，由于低層臺風中心位置偏北，南風北推，地形抬升作用減弱消失，加上山脈附近水汽輻合減弱，降水強度減小，最大減幅超過40mm。

（3）地形的存在一方面影響臺風移動路徑，從而影響環流的空間分布，另一方面地形阻擋作用的消失使低層風速抬升減弱，兩方面的共同作用減弱了山脈周圍的降水。

在未來的研究工作中，可以考慮采用bogus技術對初始臺風強度及位置進行重定位，從而進一步改善臺風移動路徑及強度的模擬。

參考文獻：

［1］ Dee D P， Uppala S M， Simmons A J， et al. ，The ERA-Interim reanalysis：Configuration and performance of the data assimilation system.Quart J Roy Meteor Soc［J］，2011， 137： 553-597.

［2］陳聯壽，徐祥德，羅哲賢，等.熱帶氣旋動力學引論［M］.北京：氣象出版社，2002，17.Chen Lianshou， Xu Xiangde， Luo Zhexian， etal.Introduction to tropical cyclone dynamics ［M］.Beijing：China Meteorological Press，2002，17. （in Chinese）

［3］倪鐘萍，吳立廣，張玲.2005—2010年臺風突變路徑的預報誤差及其環流背景［J］.氣象，2013，39（6）：719-727.Ni Zhongping，Wu Liguang，Zhang Ling.Analysis on Forecasting Errors and Associated Circulations of Sudden Typhoon Track Changes During 2005-2010 ［J］.Meteorological Monthly，2013，39 （6）：719-727.（in Chinese）

［4］馮錦全，陳多.我國近海熱帶氣旋強度突變的氣候特征分析［J］.熱帶氣象學報，1995，11（1）：35-42.Feng Jinquan，Chen Duo，Analyses of Climatic Characteristics Accompanying Sudden intensity Changes in offshore tropical cyclones in china ［J］.Journal of Tropical Meteorology，1995，11（1）：35-42.（in Chinese）

［5］程正泉，林良勛，楊國杰，等.超強臺風威馬遜快速增強及大尺度環流特征［J］.應用氣象學報，2017，28（3）：318-326.Cheng Zhengquan，Lin Liangxun，YangGuojie，etal.Rapid Intensification and Associated Large-scale Circulation of superTyphoon Rammasun in 2014 ［J］.Journalof Applied Meteorological Science，Beijing：2017，28 （3）：318-326.（in Chinese）

［6］徐舒揚，周德麗，苗紹慧，等.臺風“威馬遜”造成云南強降水災害天氣分析［J］.災害學，2016，31（4）：229-234.Xu Shuyang，Zhou Deli，Miao Shaohui，etal.An Analysis about the Disastrous WeatherofHeavyPrecipitation Resulted from the Typhoon Rammasun in Yunnan Province ［J］.Journal of Catastrophology，2016，31 （4）：229-234.（in Chinese）

［7］朱紅芳，王東勇，婁珊珊，等.地形對臺風“海葵”降水增幅影響的研究［J］.暴雨災害，2015，34（2）：160-167.Zhu Hongfang，Wang Dongyong，Lou Shanshan，etal.Numerical test of topography effect on rainfall amplification associated with typhoon HaiKui ［J］.Torrential Rain and Disasters，2015，34（2）：160-167.（in Chinese）

［8］楊仁勇，閔錦忠，鄭艷.強臺風“納沙”引發的特大暴雨過程數值試驗［J］.高原氣象，2014，33（3）：753-761.Yang Renyong，Min Jinzhong，Zheng Yan.Numerical Simulations of the Extraordinary Rainstorm by Typhoon Nesat［J］.Plateau Meteorology，2014，33（3）：753-761.（in Chinese）

［9］ Sherrie， F， Davis， C， Gill， D.Low-Nam， S.2009.Bogussing of tropical cyclones in WRF Version 3.1.Presented at the 10th Annual WRF Users Workshop［J］， Boulder， Colorado， USA， 10 June 2009.