張曉慧 朱名名 鄭瑩

摘??要：本文采用中尺度WRF模式對(duì)2013年雙臺(tái)風(fēng)“Fitow”和“Danas”進(jìn)行了不同參數(shù)化方案敏感性試驗(yàn)，研究結(jié)果表明Kessler微物理方案和Kain-Fritsch積云參數(shù)化方案的組合對(duì)臺(tái)風(fēng)路徑模擬效果最好;YSU邊界層方案與Kain-Fritsch積云參數(shù)化方案的組合對(duì)臺(tái)風(fēng)強(qiáng)度模擬效果最好;綜合分析物理過程中的Kessler方案、Kain-Fritsch方案和YSU方案是最佳的組合，主要是對(duì)直接影響雙臺(tái)風(fēng)路徑、強(qiáng)度的西太平洋副熱帶高壓和西風(fēng)帶槽脊的位置和強(qiáng)度模擬效果好。

關(guān)鍵詞：數(shù)值模擬;WRF模式;雙臺(tái)風(fēng);物理方案

1?概述

兩個(gè)臺(tái)風(fēng)靠近時(shí)它們的路徑、強(qiáng)度等都會(huì)收到彼此的影響，早在1923年Fujiharal[1]就通過一系列的渦旋實(shí)驗(yàn)及觀測(cè)發(fā)現(xiàn)兩個(gè)接近的水旋渦繞著兩者連線的中心互相旋轉(zhuǎn)。1979年曾慶存[2]則從大氣動(dòng)力學(xué)理論上研究了兩個(gè)大氣渦旋的相互作用以及其移動(dòng)方式。王玉清等[3]采用正壓無輻散模式模擬發(fā)現(xiàn)兩個(gè)同等強(qiáng)度的渦旋“合并”后體現(xiàn)了西北向的平移，而不同強(qiáng)度的兩個(gè)渦旋“合并”后體現(xiàn)了陀螺運(yùn)動(dòng)和平移相疊加的形態(tài)，這與實(shí)際中雙臺(tái)風(fēng)的相互作用表現(xiàn)類似。張曉慧等[4-5]采用集合變分混合同化的方法研究了臺(tái)風(fēng)Fitow（1324）登陸前、后的強(qiáng)度變化與臺(tái)風(fēng)Danas（1323）相連的水汽輸送通道的變化、相互作用等直接相關(guān)。陳德輝等[6]研究發(fā)現(xiàn)對(duì)模式物理過程參數(shù)化方案改進(jìn)后能明顯提高臺(tái)風(fēng)路徑和強(qiáng)度的預(yù)報(bào)效果。Li?等[7]發(fā)現(xiàn)邊界層參數(shù)化方案對(duì)颶風(fēng)的早期發(fā)展非常敏感。所以有必要進(jìn)行單獨(dú)個(gè)例分析和研究，這樣有助于改進(jìn)雙臺(tái)風(fēng)的天氣預(yù)報(bào)水平，。，對(duì)具有自身明顯特點(diǎn)和造成較大災(zāi)害的，更是如此。

本文采用中尺度數(shù)值模式WRFV3.3，物理過程則使用三種微物理方案，三種邊界層方案和三種積云對(duì)流參數(shù)化方案的不同組合配置，來討論對(duì)雙臺(tái)風(fēng)“菲特”（Fitow）和“丹娜絲”（Danas）的路徑、強(qiáng)度和降水模擬的影響，以揭示?WRF?模式中物理過程的組合配置方案在該雙臺(tái)風(fēng)模擬中的作用，也為集合預(yù)報(bào)中各物理集合成員的選取打下基礎(chǔ)。

2?本文雙臺(tái)風(fēng)簡介

雙臺(tái)風(fēng)中的1323號(hào)臺(tái)風(fēng)“Fitow”于2013年9月30日12時(shí)開始編號(hào)（熱帶風(fēng)暴級(jí)）后，路徑沿西北偏北方向移動(dòng)，10月04日08時(shí)“Fitow”（強(qiáng)臺(tái)風(fēng)級(jí)）轉(zhuǎn)向西北偏西方向移動(dòng)（此時(shí)“Danas”編號(hào)后2h），10月6日15時(shí)強(qiáng)臺(tái)風(fēng)“Fitow”登陸福建福鼎市沙埕鎮(zhèn)，登陸后路徑稍向南折?！癉anas”生成編號(hào)到10月6日15時(shí)“Fitow”登陸時(shí)，“Danas”一直沿西北方向移動(dòng)，10月7日18時(shí)轉(zhuǎn)折北上，10月8日00時(shí)由向北運(yùn)動(dòng)轉(zhuǎn)為向東北方向移動(dòng)（圖1）。

3?所用資料和數(shù)值模擬方案設(shè)計(jì)

本文采用的數(shù)值模擬模式為非靜力的中尺度數(shù)值模式WRF。為揭示不同參數(shù)化方案對(duì)雙臺(tái)風(fēng)模擬的影響，利用WRF模式中自帶的三種微物理參數(shù)化方案、三種邊界層參數(shù)化方案以及三種積云對(duì)流參數(shù)化方案（在表1和下文中分別簡稱為微物理（方案）、邊界層（方案）、積云（方案））的不同的組合，設(shè)計(jì)了27組試驗(yàn)，具體見表1。模擬初始場(chǎng)和側(cè)邊界條件采用1°×?1°的FNL?（?Final?Global?Data?Assimilation?System）?再分析資料，時(shí)間間隔為6h，垂直方向26層，模式水平格距為30km，格點(diǎn)數(shù)為200×120，垂直方向分為28層，模式層頂為50hPa，時(shí)間步長為120s，模擬時(shí)段為：2013年10月04日00時(shí)（世界時(shí)，以下同）至2013年10月07日00時(shí)。

4?數(shù)值模擬結(jié)果分析

4.1對(duì)臺(tái)風(fēng)路徑的影響

為揭示微物理過程對(duì)雙臺(tái)風(fēng)路徑模擬的影響，在固定邊界層和積云方案后，對(duì)微物理過程選用Kessler方案較Lin方案和Ferrier方案模擬的誤差小。Kessler微物理方案是暖云方案，主要考慮雨水的產(chǎn)生、降落以及蒸發(fā)，云水的凝結(jié)等過程，而Lin/Ferrier方案部分的考慮了云冰等冷云過程。此次數(shù)值模擬過程中Kessler優(yōu)于Lin/Ferrier方案，這說明對(duì)低緯臺(tái)風(fēng)30Km分辨率的模擬，考慮冷云過程并不一定能提高預(yù)報(bào)效果。

固定微物理方案和邊界層方案，積云方案選用Kain-Fritsch則要較Betts-Miller-Janjic?以及Grell-Devenyi方案模擬的臺(tái)風(fēng)路徑誤差小，對(duì)于臺(tái)風(fēng)“Danas”路徑的模擬，幾種方案差異不大，但總體來講選用有利于熱對(duì)流混合層高度升高，使邊界層結(jié)構(gòu)更接近中性的Kain-Fritsch方案效果最好。

為揭示邊界層方案對(duì)臺(tái)風(fēng)路徑模擬效果的影響，采用相同的微物理方案和積云物理方案，不同的邊界層方案來模擬雙臺(tái)風(fēng)路徑并與實(shí)況資料作對(duì)比分析發(fā)現(xiàn)在微物理過程和積云方案固定的情況下，邊界層方案無論采用何種方案，對(duì)臺(tái)風(fēng)移動(dòng)路徑影響都不大;但對(duì)臺(tái)風(fēng)Danas的模擬，MRF方案模擬結(jié)果與實(shí)況更要接近些，總體來看MRF方案模擬效果更好些。

對(duì)比試驗(yàn)01（Kain-Fritsch）和試驗(yàn)13（Betts-Miller-Janji）發(fā)現(xiàn)，在積分72?h時(shí)（07日00時(shí)）模擬的臺(tái)風(fēng)路徑差異最顯著，此時(shí)試驗(yàn)01（Kain-Fritsch）方案模擬的副高中心位置、強(qiáng)度與實(shí)況較接近，受副高外圍引導(dǎo)氣流影響，模擬的“Fitow”移動(dòng)路徑較實(shí)況偏北。而臺(tái)風(fēng)“Danas”在04日-07日期間基本位于副高南側(cè)，故副高位置和強(qiáng)度對(duì)其影響不大。試驗(yàn)13方案模擬的副高脊線呈西北偏西向，故“fitow”和“Danas”受副高南側(cè)引導(dǎo)氣流的影響，移動(dòng)路徑較實(shí)況偏南。由于雙臺(tái)風(fēng)的副高引導(dǎo)氣流都較實(shí)況偏弱，所以模擬的雙臺(tái)風(fēng)移速都較實(shí)況偏慢?？傮w對(duì)比分析可見，副高中心位置，副高西伸脊點(diǎn)位置和副高強(qiáng)度都影響到雙臺(tái)風(fēng)的移動(dòng)方向和移動(dòng)速度。

以臺(tái)風(fēng)中心海平面氣壓作為衡量臺(tái)風(fēng)強(qiáng)度的標(biāo)準(zhǔn)，選用相同的邊界層方案和積云方案，不同的微物理方案組合對(duì)雙臺(tái)風(fēng)強(qiáng)度模擬研究發(fā)現(xiàn)不同的微物理方案模擬臺(tái)風(fēng)“Fitow”的強(qiáng)度差別不大，且都較實(shí)況強(qiáng)度偏弱。但對(duì)“Danas”強(qiáng)度的模擬差別較大?？傮w來看48h內(nèi)微物理過程選用Kessler方案較Lin方案和Ferrier方案模擬的臺(tái)風(fēng)強(qiáng)度誤差要小。

選用相同的邊界層方案和微物理方案，不同的積云方案對(duì)雙臺(tái)風(fēng)強(qiáng)度的模擬研究發(fā)現(xiàn)不同的微物理方案模擬“Fitow”?的強(qiáng)度，積云方案采用Kain-Fritsch方案要比Betts-Miller-Janjic、Grell-Devenyi方案模擬的臺(tái)風(fēng)強(qiáng)度誤差小。試驗(yàn)01方案模擬的副高和中緯度西風(fēng)帶槽脊的強(qiáng)度較試驗(yàn)13強(qiáng)，而模擬的雙臺(tái)風(fēng)強(qiáng)度也較試驗(yàn)13強(qiáng)（圖2）?？傮w對(duì)比分析可見，副高、西風(fēng)帶槽脊位置和強(qiáng)度變化對(duì)雙臺(tái)風(fēng)的強(qiáng)度有影響。

敏感性試驗(yàn)選用相同的微物理方案和積云方案，不同的邊界層方案組合對(duì)雙臺(tái)風(fēng)強(qiáng)度的模擬結(jié)果顯示所有方案模擬的“Fitow”臺(tái)風(fēng)強(qiáng)度都較實(shí)況偏弱，YSU方案較MRF方案、MYNN2.5方案方案模擬的效果要好。對(duì)“Danas”強(qiáng)度的模擬，MRF方案模擬的臺(tái)風(fēng)強(qiáng)度無論是大小還是臺(tái)風(fēng)減弱、增強(qiáng)的趨勢(shì)，都比YSU、MYNN2.5方案的效果要好。

結(jié)語

通過27種不同參數(shù)化方案組合對(duì)2013年雙臺(tái)風(fēng)“Fitow”和“Danas”進(jìn)行了數(shù)值模擬，得到以下主要結(jié)論：

（1）微物理方案的不同和積云方案的不同對(duì)臺(tái)風(fēng)路徑的影響明顯，其中Kessler微物理方案和Kain-Fritsch積云方案的組合，模擬的臺(tái)風(fēng)路徑與實(shí)況的平均偏差最小。

（2）邊界層方案與積云方案都對(duì)臺(tái)風(fēng)強(qiáng)度模擬效果影響明顯，YSU邊界層方案與Kain-Fritsch積云方案的組合模擬的臺(tái)風(fēng)強(qiáng)度變化與實(shí)況更接近。

（3）模式中不同物理過程的組合對(duì)模擬臺(tái)風(fēng)的路徑、強(qiáng)度等影響較大，物理過程中的Kessler方案、Kain-Fritsch方案和YSU方案是對(duì)雙臺(tái)風(fēng)“Fitow”和“Danas”模擬的最佳的組合方案。

（4）西太副高和西風(fēng)帶槽脊的位置和強(qiáng)度變化直接影響了雙臺(tái)風(fēng)“Fitow”和“Danas”的移動(dòng)路徑及強(qiáng)度。

參考文獻(xiàn)：

[1]?FujiwharalS，1923.?The?mutual?tendency?towards?symmetry?of?motion?and?its?application?as?a?principle?in?meteorology[J].Quart?J?Roy?Meteor?Soc，?1923，49（12）：287-93.

[2]曾慶存.數(shù)值天氣預(yù)報(bào)的數(shù)學(xué)物理基礎(chǔ)（第一卷）[M].北京：科學(xué)出版社，1979：421-440.

[3]王玉清，?朱永禔.?正壓無輻散模式中雙渦的相互作用[J].熱帶氣象，1989，5（2）：?105-115.

[4]王玉清，?朱永禔.?雙渦藤原效應(yīng)的數(shù)值模擬研究[J].應(yīng)用氣象學(xué)報(bào)，1989（1）：13-20.

[5]張曉慧，張立鳳，周海申，魏通峰.?雙臺(tái)風(fēng)相互作用及其影響[J].?應(yīng)用氣象學(xué)報(bào)，2019，30（04）：456-466.

[6]張曉慧，張立鳳，熊春暉，彭軍.?基于集合變分混合同化方法的雙臺(tái)風(fēng)數(shù)值模擬[J].?熱帶氣象學(xué)報(bào)，2015，31（04）：505-516.

[7]?Li?Y，?Ellingwood?B?R.?Hurricane?damage?to?residential?construction?in?the?US：?Importance?of?uncertainty?modeling?in?risk?assessment[J].?Engineering?Structures，?2006，28（7）：1009-1018.

基金項(xiàng)目：“雙重”兵種作戰(zhàn)《基于大數(shù)據(jù)的航空兵作戰(zhàn)指揮氣象輔助決策系統(tǒng)》，課題號(hào)：SZ063-BZ15-JY03

作者簡介：張曉慧（1989—??），女，甘肅天水人，大氣科學(xué)碩士，講師

研究方向：數(shù)值天氣預(yù)報(bào)與航空危險(xiǎn)天氣預(yù)報(bào)。