袁炳，費(fèi)建芳，王云峰，韓月琪，盧強(qiáng)

（解放軍理工大學(xué) 氣象學(xué)院，江蘇 南京 211101）

一種非對(duì)稱臺(tái)風(fēng)Bogus方法的數(shù)值模擬應(yīng)用

袁炳，費(fèi)建芳，王云峰，韓月琪，盧強(qiáng)

（解放軍理工大學(xué) 氣象學(xué)院，江蘇 南京 211101）

前人研究中BDA方法采用的軸對(duì)稱Bogus臺(tái)風(fēng)不能反映個(gè)別臺(tái)風(fēng)具體特征，也棄掉了背景場(chǎng)的合理成分，有必要改用更精細(xì)的非對(duì)稱臺(tái)風(fēng)Bogus模型對(duì)BDA方案做出改進(jìn)。因此，提出一種充分融合分析場(chǎng)信息和實(shí)際觀測(cè)信息并考慮副高影響的非對(duì)稱臺(tái)風(fēng)Bogus方法，針對(duì)0613、0704臺(tái)風(fēng)個(gè)例，通過(guò)4DVAR技術(shù)，利用MM5模式及其伴隨模式對(duì)此Bogus資料進(jìn)行BDA同化試驗(yàn)和模擬預(yù)報(bào)，數(shù)值結(jié)果表明：BDA方法同化非對(duì)稱臺(tái)風(fēng)Bogus模型，其路徑預(yù)報(bào)效果優(yōu)于同化軸對(duì)稱臺(tái)風(fēng)模型，同化兩種Bogus資料都使臺(tái)風(fēng)強(qiáng)度模擬效果得到極大改善，但改善程度相當(dāng)，故臺(tái)風(fēng)的非對(duì)稱結(jié)構(gòu)及副高等外圍環(huán)境場(chǎng)及背景場(chǎng)弱臺(tái)風(fēng)中合理信息對(duì)臺(tái)風(fēng)的準(zhǔn)確模擬起到關(guān)鍵作用，同時(shí)，Bogus資料同化方法用于遠(yuǎn)洋上的臺(tái)風(fēng)時(shí)，需注意臺(tái)風(fēng)模型與副高的配置關(guān)系。

非對(duì)稱臺(tái)風(fēng)；四維變分資料同化；BDA方案；路徑和強(qiáng)度預(yù)報(bào)

臺(tái)風(fēng)路徑、強(qiáng)度和風(fēng)雨的精確預(yù)報(bào)對(duì)于數(shù)值預(yù)報(bào)而言是個(gè)難題，除了模式的物理過(guò)程尤其是對(duì)海氣相互作用和對(duì)流全面描述的不準(zhǔn)確之外，大量洋面觀測(cè)資料的缺乏使熱帶氣旋內(nèi)部動(dòng)力和熱力結(jié)構(gòu)及大尺度環(huán)流場(chǎng)得不到詳盡描述[1]。因而業(yè)務(wù)上大尺度分析場(chǎng)的臺(tái)風(fēng)強(qiáng)度弱且位置不準(zhǔn)確，國(guó)內(nèi)外的解決辦法基本上是消除弱且不準(zhǔn)確的擾動(dòng)渦旋場(chǎng)，然后根據(jù)觀測(cè)信息和經(jīng)驗(yàn)公式構(gòu)造出一個(gè)三維渦旋環(huán)流場(chǎng)并嵌入模式初始分析場(chǎng)，其Bogus方法宗旨基本上是設(shè)計(jì)一個(gè)軸對(duì)稱臺(tái)風(fēng)渦旋并加入某種能表征臺(tái)風(fēng)及其環(huán)境氣流運(yùn)動(dòng)的非對(duì)稱分量。很多工作表明，這些臺(tái)風(fēng)模型初值化方案使臺(tái)風(fēng)路徑的預(yù)報(bào)水平有了較大的改進(jìn)[2-3]。

各種Bogus方案都遇到臺(tái)風(fēng)模型及有關(guān)參數(shù)尤其環(huán)境氣壓和最大風(fēng)速半徑等的合理確定問題。章家琳等(1989)[4]認(rèn)為由于摩擦作用使風(fēng)向穿越軸對(duì)稱等壓線而產(chǎn)生一個(gè)風(fēng)向內(nèi)偏角，并采用不同氣壓模式求取此內(nèi)偏角進(jìn)而得出風(fēng)場(chǎng)模型，在此基礎(chǔ)上，胡邦輝等(1999)[5]采用藤田公式進(jìn)一步導(dǎo)出了最大風(fēng)速半徑的計(jì)算式，這樣的計(jì)算方法是依賴特定氣壓場(chǎng)模型得出的，氣壓模型的有效性直接影響計(jì)算結(jié)果的正確與否。房文鸞等(1987)[6]依據(jù)天氣圖中的觀測(cè)信息擬合出通用氣壓模型中表征最大風(fēng)速半徑特性的臺(tái)風(fēng)常數(shù)，袁金南等(2007)[7]也引入十級(jí)和七級(jí)風(fēng)圈平均半徑來(lái)改變軸對(duì)稱模型臺(tái)風(fēng)的切向風(fēng)廓線形態(tài)， Lajoie等(2008)[8]利用衛(wèi)星云圖特征能夠較準(zhǔn)確的定位出臺(tái)風(fēng)最大風(fēng)速半徑。然而一方面由于風(fēng)速半徑的觀測(cè)和計(jì)算都存在誤差且較難實(shí)施，另一方面由于模式并不完美，采用與實(shí)際觀測(cè)完全逼近的臺(tái)風(fēng)模型不一定與模式本身相協(xié)調(diào)而使預(yù)報(bào)效果不一定得以改善。所以有學(xué)者試圖尋求避免人為給定臺(tái)風(fēng)常數(shù)的方法，比如陳孔沫等(1989)[9]、王國(guó)民等(1996)[10]、楊支中等(2005)[11]引入臺(tái)風(fēng)外圍閉合特征等壓線來(lái)對(duì)藤田氣壓模型進(jìn)行改進(jìn)，推導(dǎo)出不含臺(tái)風(fēng)常數(shù)的橢圓、對(duì)稱或非對(duì)稱的臺(tái)風(fēng)海面氣壓場(chǎng)。對(duì)于環(huán)境氣壓的給定，前人的多數(shù)研究都是假定其為時(shí)空上的常值，而胡邦輝等(1999)[5]則考慮了臺(tái)風(fēng)受副高影響產(chǎn)生環(huán)境氣壓的方向性不均勻特征。

隨著4DVAR技術(shù)的發(fā)展[12-14]，Zou等(1999)[15]提出的BDA (Bogus data assimilation) 方法使臺(tái)風(fēng)路徑預(yù)報(bào)得到了改善。之后不少學(xué)者對(duì)BDA方案作了進(jìn)一步研究并取得了較好效果[16-18]。然而此方法采用的軸對(duì)稱Bogus臺(tái)風(fēng)不能反映某一個(gè)臺(tái)風(fēng)具體特征，并假設(shè)臺(tái)風(fēng)同化窗內(nèi)靜止不動(dòng)，也拋棄掉了原背景場(chǎng)弱臺(tái)風(fēng)中的合理成分。所以，有必要改用更為精細(xì)的非對(duì)稱臺(tái)風(fēng)Bogus模型來(lái)對(duì)BDA方案做出改進(jìn)。

1 非對(duì)稱臺(tái)風(fēng)Bogus方案

根據(jù)文獻(xiàn)[5]，在以臺(tái)風(fēng)中心為原點(diǎn)的極坐標(biāo)中，考慮海面摩擦作用，假定成熟臺(tái)風(fēng)呈穩(wěn)定狀態(tài)，并在慣性項(xiàng)中考慮臺(tái)風(fēng)移動(dòng)對(duì)曲率半徑的影響，同時(shí)做密度變形即d=1/ ，則臺(tái)風(fēng)域內(nèi)海面任意空氣質(zhì)點(diǎn)(r,)的水平運(yùn)動(dòng)方程為：

圖3 臺(tái)風(fēng)模型中R/r0（最大風(fēng)速半徑與臺(tái)風(fēng)常數(shù)比值）及vmax（最大風(fēng)速）隨a（圖A）和b（圖B）的變化曲線Fig. 3 Variation curve of R/r0 (ratio of MWR and r0) and vmax (max wind) due to a (Fig. A) and b (Fig. B) in the typhoon bogus model

圖4是構(gòu)造出的0613號(hào)臺(tái)風(fēng)海面非對(duì)稱模型。臺(tái)風(fēng)氣壓場(chǎng)非對(duì)稱結(jié)構(gòu)體現(xiàn)了環(huán)境氣壓的影響，而最大風(fēng)速區(qū)則出現(xiàn)在靠近副高一側(cè)的具有較大氣壓梯度的區(qū)域。

圖4 構(gòu)造出的世界時(shí)06年9月14日12時(shí)0613臺(tái)風(fēng)的非對(duì)稱Bogus模型。A： 海平面氣壓（中心氣壓960.3 hPa）; B： 近海面風(fēng)速（最大風(fēng)速41.2 m/s）Fig. 4 Bogus model of 0613 typhoon in UTC 2006-09-14_12：00 A： Sea level pressure (center pressure 960.3 hPa); B： Sea level wind (max wind 41.2 m/s)

2 臺(tái)風(fēng)個(gè)例數(shù)值模擬

為了驗(yàn)證改進(jìn)得到的臺(tái)風(fēng)非對(duì)稱氣壓場(chǎng)和風(fēng)場(chǎng)模型的模擬效果，本文選取0613號(hào)臺(tái)風(fēng)9月14日12時(shí)至17日12時(shí)以及0704號(hào)臺(tái)風(fēng)7月11日18時(shí)至14日18時(shí)的時(shí)段作72小時(shí)模擬。數(shù)值試驗(yàn)中采用的數(shù)值模式為非靜力版本的 MM5模式及其伴隨模式，同化窗口為32分鐘。試驗(yàn)區(qū)域以(124.5E, 25.5N)為中心，兩重固定雙向嵌套。粗網(wǎng)格為75 × 91，格距54 km，細(xì)網(wǎng)格為136 × 130，格距18 km，模式垂直方向?yàn)榉植疾痪鶆虻?3層,資料為6小時(shí)間隔的1o× 1o的Ncep再分析資料，積云參數(shù)化為GRELL方案，邊界層參數(shù)化為MRF方案。試驗(yàn)設(shè)計(jì)了如表1所示的3種初值化方案，非對(duì)稱Bogus氣壓場(chǎng)資料由(14b)式得到，軸對(duì)稱臺(tái)風(fēng)Bogus資料與文獻(xiàn)[15]相同，Bogus臺(tái)風(fēng)模型只包括海平面風(fēng)場(chǎng)和氣壓場(chǎng)。臺(tái)風(fēng)實(shí)測(cè)資料來(lái)自中國(guó)臺(tái)風(fēng)網(wǎng)(www.typhoon.gov.cn) “CMA-STI熱帶氣旋最佳路徑數(shù)據(jù)集”。

表1 數(shù)值試驗(yàn)方案設(shè)計(jì)Tab. 1 Experiment schemes

圖5為路徑預(yù)報(bào)距離誤差。對(duì)于0613臺(tái)風(fēng)，前48小時(shí)內(nèi)經(jīng)兩種BDA方法初值化后預(yù)報(bào)效果都優(yōu)于控制實(shí)驗(yàn)（控制實(shí)驗(yàn)誤差較大且呈無(wú)規(guī)則振蕩，故其72時(shí)誤差優(yōu)于兩個(gè)同化方案可能為偶然性），方案3效果總體又優(yōu)于方案2；而在48小時(shí)之后，采用軸對(duì)稱Bogus模型的方案2開始出現(xiàn)較大誤差，而采用非對(duì)稱Bogus資料的方案3卻仍取得了較好改進(jìn)效果。對(duì)于0704臺(tái)風(fēng)，控制實(shí)驗(yàn)中路徑預(yù)報(bào)出現(xiàn)較大誤差且隨時(shí)間呈無(wú)規(guī)則變化；方案2在18時(shí)至30時(shí)路徑有所改進(jìn)，其余時(shí)段都與實(shí)況發(fā)生嚴(yán)重偏差，誤差超過(guò)了控制實(shí)驗(yàn)；方案3也只在12時(shí)至36時(shí)路徑有所改善，其他時(shí)段也都出現(xiàn)大于控制實(shí)驗(yàn)的誤差。0704臺(tái)風(fēng)路徑模擬出現(xiàn)這種情況的原因可從三種方案的副高演變（圖略）來(lái)分析，初始時(shí)刻大尺度背景場(chǎng)的副高強(qiáng)度偏弱，主體偏東，這與 0704臺(tái)風(fēng)北上路徑的對(duì)應(yīng)關(guān)系與文獻(xiàn)[26]的分析是一致的，0704臺(tái)風(fēng)此時(shí)位于離陸地較遠(yuǎn)的海面上副高引導(dǎo)氣流內(nèi)，其移動(dòng)路徑與副高密切相關(guān)?？刂茖?shí)驗(yàn)（方案1）中弱臺(tái)風(fēng)與副高協(xié)調(diào)發(fā)展，臺(tái)風(fēng)轉(zhuǎn)向前副高維持不變，轉(zhuǎn)向后副高稍有西伸；方案2中初始時(shí)刻副高邊緣動(dòng)力場(chǎng)分布因受到軸對(duì)稱Bogus資料的調(diào)整而發(fā)生很大改變，模擬初期由于同化過(guò)程的協(xié)調(diào)作用能取得好的預(yù)報(bào)效果，但24時(shí)副高受Bogus臺(tái)風(fēng)影響減弱東退（并在模擬后期破裂），臺(tái)風(fēng)也開始發(fā)生轉(zhuǎn)向（實(shí)況中臺(tái)風(fēng)于42時(shí)轉(zhuǎn)向），之后誤差不斷增大；方案 3中由于非對(duì)稱臺(tái)風(fēng)模型考慮了基本場(chǎng)信息及副高影響，副高及臺(tái)風(fēng)外圍環(huán)境場(chǎng)受到Bogus資料的沖擊比方案2小，但由于Bogus臺(tái)風(fēng)強(qiáng)度遠(yuǎn)大于原來(lái)的弱臺(tái)風(fēng)，使得臺(tái)風(fēng)外圍環(huán)境尤其副高的演變也較控制實(shí)驗(yàn)有所不同，介于方案1和方案3之間，臺(tái)風(fēng)模擬36時(shí)后開始轉(zhuǎn)向，之后誤差也開始走高。單獨(dú)比較方案2、3，整個(gè)模擬過(guò)程方案3的路徑誤差幾乎都小于方案2，同樣能體現(xiàn)非對(duì)稱Bogus資料的優(yōu)越性，同時(shí)也說(shuō)明BDA方法尚不能適用于任何海上臺(tái)風(fēng)，需注意臺(tái)風(fēng)與副高的配置關(guān)系。

圖6為預(yù)報(bào)過(guò)程中各個(gè)臺(tái)風(fēng)的中心氣壓變化曲線。兩個(gè)BDA試驗(yàn)?zāi)M的臺(tái)風(fēng)強(qiáng)度都比控制實(shí)驗(yàn)更接近實(shí)況。這種模擬出的強(qiáng)度遠(yuǎn)比控制試驗(yàn)強(qiáng)的情況體現(xiàn)了BDA方法在臺(tái)風(fēng)強(qiáng)度調(diào)整上的特點(diǎn)，改善了臺(tái)風(fēng)成熟期的強(qiáng)度預(yù)報(bào)，但 60小時(shí)之后當(dāng)實(shí)況中臺(tái)風(fēng)強(qiáng)度已大為減弱時(shí)，試驗(yàn)?zāi)M的臺(tái)風(fēng)仍然保持很強(qiáng)的強(qiáng)度，造成氣旋過(guò)度發(fā)展。0613號(hào)臺(tái)風(fēng)前36小時(shí)方案3模擬的強(qiáng)度比方案2偏弱，其余時(shí)段中心氣壓與方案2相當(dāng)，0704臺(tái)風(fēng)中方案3模擬出的中心氣壓都比方案2稍低。兩種同化方案的強(qiáng)度模擬差異較小，然而路徑模擬誤差卻相差巨大，說(shuō)明臺(tái)風(fēng)的非對(duì)稱結(jié)構(gòu)及副高等外圍環(huán)境場(chǎng)以及初始背景場(chǎng)中的合理大尺度信息對(duì)臺(tái)風(fēng)的準(zhǔn)確模擬起到很關(guān)鍵的作用。

圖5 0613，0704臺(tái)風(fēng)路徑預(yù)報(bào)距離誤差，橫坐標(biāo)為時(shí)間 (h)；縱坐標(biāo)為誤差 (km)Fig. 5 Track forecast distance error of 0613, 0704 typhoons. Abscissa axis is time (h); Vertical axis is distance error (km)

圖6 0613，0704臺(tái)風(fēng)中心氣壓變化曲線，橫坐標(biāo)為時(shí)間(h)，縱坐標(biāo)為中心氣壓 (hPa)Fig. 6 Variation curve of center pressure of 0613, 0704 typhoons. Abscissa axis is time (h); Vertical axis is center pressure (hPa)

3 結(jié) 論

a) 提出的非對(duì)稱Bogus方法包括：從考慮摩擦的水平運(yùn)動(dòng)方程導(dǎo)出卡當(dāng)方程形式的風(fēng)速公式并求解；對(duì)臺(tái)風(fēng)外圍閉合等壓線進(jìn)行傅立葉級(jí)數(shù)插值擬合；計(jì)算副高影響下各個(gè)方向的環(huán)境氣壓；對(duì)Myers公式作非對(duì)稱推廣和風(fēng)廓線約束；Bogus風(fēng)場(chǎng)吸收濾除弱渦旋并扣除臺(tái)風(fēng)移速后的基本場(chǎng)成分。該方法全程以非對(duì)稱形式計(jì)算，避免了假定臺(tái)風(fēng)軸對(duì)稱時(shí)引入非對(duì)稱量的難題，也避免了假定臺(tái)風(fēng)呈橢圓時(shí)所帶來(lái)的難以確定風(fēng)場(chǎng)長(zhǎng)短軸方向的難題，還避免了給定最大風(fēng)速半徑時(shí)的不確定性。該方法考慮非對(duì)稱性質(zhì)并吸收基本場(chǎng)信息及估算副高影響，只針對(duì)臺(tái)風(fēng)初值化設(shè)計(jì)，其他場(chǎng)合未必適用。

b) 對(duì)0613、0704臺(tái)風(fēng)的數(shù)值試驗(yàn)對(duì)比發(fā)現(xiàn)， BDA方法同化非對(duì)稱臺(tái)風(fēng)Bogus模型，其路徑預(yù)報(bào)效果優(yōu)于同化軸對(duì)稱臺(tái)風(fēng)模型，體現(xiàn)出新Bogus資料其數(shù)值同化及模擬上的優(yōu)越性，也說(shuō)明了臺(tái)風(fēng)的非對(duì)稱結(jié)構(gòu)及副高等外圍環(huán)境場(chǎng)以及初始背景場(chǎng)中的合理大尺度信息對(duì)臺(tái)風(fēng)的準(zhǔn)確模擬起到關(guān)鍵作用。同化兩種Bogus資料都使臺(tái)風(fēng)強(qiáng)度模擬效果得到極大改善，但渦旋過(guò)度發(fā)展現(xiàn)象需進(jìn)一步解決。

c) Bogus資料同化方法用于遠(yuǎn)洋上的臺(tái)風(fēng)時(shí)，需注意臺(tái)風(fēng)模型與副高的配置關(guān)系。

[1] 瞿安詳, 麻素紅. 非對(duì)稱臺(tái)風(fēng)Bogus方案設(shè)計(jì)和初步試驗(yàn) [J]. 應(yīng)用氣象學(xué)報(bào), 2007, 18(3)： 380-387.

[2] Iwasaki T, Nakano H, Sugi M. The performance of a typhoon track prediction model [J]. J Meteor Soc Japan, 1987, 65： 555-570.

[3] Kurihara, Y, Bender M A, Ross R J. An initialization scheme of hurricane models by vortex specification [J]. Mon Wea Rev, 1993,121： 2030-2045.

[4] 章家琳, 隋世峰. 臺(tái)風(fēng)波浪數(shù)值預(yù)報(bào)的CHGS法Ⅱ： 臺(tái)風(fēng)風(fēng)場(chǎng)中梯度風(fēng)摩擦修正系數(shù)和風(fēng)向內(nèi)偏角的計(jì)算 [J]. 熱帶海洋, 1989, 8(1)： 58-66.

[5] 胡邦輝, 譚言科, 張學(xué)敏. 海面熱帶氣旋域內(nèi)風(fēng)速分布 [J]. 大氣科學(xué), 1999, 23(3)： 316-322.

[6] 房文鸞, 徐啟明, 章家琳, 等. 東海臺(tái)風(fēng)域內(nèi)的風(fēng)速分布 [J]. 海洋預(yù)報(bào), 1987, 4(3)： 1-14.

[7] 袁金南, 劉春霞. 軸對(duì)稱臺(tái)風(fēng)模型的改進(jìn)方案及其對(duì)0425號(hào)臺(tái)風(fēng)數(shù)值模擬的效果 [J]. 熱帶氣象學(xué)報(bào), 2007, 23(3)： 236-245.

[8] Lajoie F, Kevin W. A technique to determine the radius of maximum wind of a tropical cyclone [J]. Weather and Forecasting, 2008, 23(5)： 1 007-1 015.

[9] 陳孔沫, 秦曾灝. 臺(tái)風(fēng)海面最大風(fēng)速的計(jì)算 [J]. 海洋學(xué)報(bào), 1989, 11(1)： 31-41.

[10] 王國(guó)民, 王詩(shī)文, 李建軍. 一個(gè)人造臺(tái)風(fēng)方案及其在移動(dòng)套網(wǎng)格模式中的應(yīng)用 [J]. 熱帶氣象學(xué)報(bào), 1996, 12(1)： 9-15.

[11] 楊支中, 沙文鈺, 朱首賢, 等. 一種新型的非對(duì)稱臺(tái)風(fēng)海面氣壓場(chǎng)和風(fēng)場(chǎng)模型 [J]. 海洋通報(bào), 2005, 24(1)： 63-68.

[12] Lewis J M, Derber J C. The use of adjoint equation to solve a variational adjustment problem with advective conetraints [J]. Tellus, 1985, 37A：309-322.

[13] LeDemit F, Talagrand O. Variational algorithms for analysis and assimilation of meteorolagical observations： theoretical aspects [J]. Tellus, 1986, 38A：97-110.

[14] Talagrand O, Curtier P. Variational assimilaion of meteorological observations with the adjoint vorticity equation. I： Theory [J]. Quart J Roy Meteor Soc,1987, 113： 1 311-1 328.

[15] Xiaolei Zou, Qingnong Xiao. Studies on the initialization and simulation of a mature hurricane using a variational bogus data assimilation scheme [J]. J Atmos Sci, 1999, 57(6)： 836-860.

[16] Qingnong Xiao, Xiaolei Zou, Bin Wang. Initialization and simulation of a landfalling hurricane using a varational bogus data assimulation scheme [J].Mon Wea Rev, 2000, 128(7)： 2 252-2 269.

[17] Zhang Xiaoyan, Wang Bin, Ji Zhongzhen, et al. Initializations and simulation of a Typhoon Using 4-Dimensional Variational Data Assimilation-Research on Typhoon Herb (1996) [J]. Advances in Atmospheric Sciences, 2003, 20(4)： 612-622.

[18] Wang D L, Liang X D, Zhao Ying, et al. A comparison of two tropical cyclone bogussing scheme. Weather and Forecasting [J]. 2008, 23(1)： 194-204.

[19] 邢富沖. 一元三次方程求解新探 [J]. 中央民族大學(xué)學(xué)報(bào)(自然科學(xué)版), 2003, 12(3)： 207-218.

[20] Fajita T. Pressure Distribution in Typhoon [J]. Geophys Mag, 1952, 23：437.

[21] Myers V A. Maximum Hurricane Winds [J]. Bul1 Alner Metero Sco, 1957, 38(4)： 227.

[22] Holland J. An analytic model of the wind and pressure profiles in hurricanes [J]. Mon Wea Rev, 1980, 108： 1 212-1 218.

[23] 鈕學(xué)新. 臺(tái)風(fēng)的內(nèi)力 [J]. 大氣科學(xué), 1983, 7(1)： 42-49.

[24] Low-Nam S, Davis C. Development of a Tropical Cyclone Bogussing Scheme for the MM5 System [C]. Preprints of the Eleventh PSU/NCAR Mesoscale Model Users Workshop, Boulder, Colorado, 2001, 130-134.

[25] 黃小剛, 費(fèi)建芳, 陸漢城. 消去分析臺(tái)風(fēng)方法的對(duì)比研究[J]. 應(yīng)用氣象學(xué)報(bào), 2006, 17(1)： 81-86.

[26] 任素玲, 劉屹岷, 吳國(guó)雄. 西太平洋副熱帶高壓和臺(tái)風(fēng)相互作用的數(shù)值試驗(yàn)研究 [J]. 氣象學(xué)報(bào), 2007, 65(3)： 329-339.

Numerical simulation application of an asymmetrical typhoon bogus scheme

YUAN Bing, FEI Jian-fang, WANG Yun-feng, HAN Yue-qi, LU Qiang

(Institute of Meteorology, PLA University of Science and Technology, Nanjing 211101, China)

An axisymmetric bogus vortex at sea level was used in the traditional bogus data assimilation (BDA)scheme. This vortex could not accurately describe the specific characteristics of typhoon. The reasonable elements of the background field were discarded also. It is necessary to switch to a more sophisticated asymmetric bogus typhoon to make improvements for BDA. Thus, an asymmetrical typhoon bogus method with blend information from the analysis and the observation is brought forward while the impact of the subtropical anticyclone is considered as well.With the fifth-generation Penn State/NCAR Mesoscale Model (MM5) and its adjoint model, the four dimensional variatiaonal data assimilation (4DVAR) technique is employed to build a BDA scheme to assimilate this asymmetrical bogus vortex. The track and intensity of NO. 0613 and NO. 0704 typhoon are simulated. The results show that, the improvement of track simulation using new bogus data is more significant than that using Symmetrical Bogus data.The intensity forecasts are greatly improved by using both two kinds of bogus data, and the improvements are in a nearly same degree. Therefore, the non-symmetrical structure of typhoon, the rational element of the weak typhoon in the background field, the subtropical anticyclone, and other external environment field, play a key role in accurate simulation of typhoons. At the same time, when BDA method is used for ocean-going typhoon, it needs to pay attention to the configuration of the relationship between the subtropical and bogus typhoon model.

asymmetrical typhoon; 4DVAR; BDA scheme; track and intensity forecast

P457.8

A

1001-6932(2010)02-0187-07

2009-04-21；

2009-09-28

國(guó)家自然科學(xué)基金項(xiàng)目10871099、40805046、40830958，公益性行業(yè)(氣象)專項(xiàng)課題(GYH(QX)2007-6-15)，國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目(973)課題 (2009CB421502)

袁炳 (1982－)，男，博士，主要從事海氣數(shù)值模擬及資料同化。電子郵箱：yuanbing.mail@163.com