王平，陳葆德，曾智華

（1.中國氣象局上海臺風(fēng)研究所，上海 200030；2.中國氣象局臺風(fēng)數(shù)值預(yù)報重點(diǎn)實(shí)驗(yàn)室，上海 200030）

海洋飛沫對熱帶氣旋邊界層結(jié)構(gòu)的影響

王平1，2，陳葆德1，2，曾智華1

（1.中國氣象局上海臺風(fēng)研究所，上海 200030；2.中國氣象局臺風(fēng)數(shù)值預(yù)報重點(diǎn)實(shí)驗(yàn)室，上海 200030）

將海洋飛沫參數(shù)化引入到一個高分辨率、非靜力中尺度WRF模式中，對0908號熱帶氣旋Morakot進(jìn)行數(shù)值模擬，探討了海洋飛沫對熱帶氣旋Morakot邊界層結(jié)構(gòu)和強(qiáng)度的影響。模擬結(jié)果表明：采用新參數(shù)化后，對熱帶氣旋Morakot的強(qiáng)度預(yù)報有改進(jìn)，但對熱帶氣旋移動路徑改進(jìn)不大；其次，通過對邊界層過程的改進(jìn)，使得眼墻區(qū)域的平均徑向風(fēng)速、切向風(fēng)速、溫度、相對濕度、垂直風(fēng)速、熱通量，降水等物理量均有增強(qiáng)，各物理量的貢獻(xiàn)對熱帶氣旋Morakot強(qiáng)度和結(jié)構(gòu)變化的影響十分重要。

海洋飛沫；熱帶氣旋；Morakot；邊界層；參數(shù)化方案

1 引言

熱帶氣旋（簡稱TC）是一種強(qiáng)烈的高影響天氣系統(tǒng)。開闊洋面上，熱帶氣旋期間的強(qiáng)風(fēng)使波浪和氣泡破碎，在大氣邊界層低層產(chǎn)生大量的海洋飛沫。由于在強(qiáng)風(fēng)條件下的海上觀測存在很大的困難，因此在模式中采用參數(shù)化的方法來計算海洋飛沫對海氣間通量交換的貢獻(xiàn)。目前的數(shù)值預(yù)報模式中，海氣或陸氣間的通量的計算都是體現(xiàn)在行星邊界層參數(shù)化中，都沒考慮海洋飛沫的作用。

1982年，Anthes［1］提出海洋飛沫的蒸發(fā)會導(dǎo)致大氣下界面的冷卻，進(jìn)而提高海氣界面的感熱交換，使風(fēng)暴強(qiáng)度增強(qiáng)。1994年，F(xiàn)airall等［2］最先提出了海洋飛沫通量的總體參數(shù)化方案（FKH方案），他們把此參數(shù)化方案加入到模式中模擬熱帶氣旋，雖然沒有明確指出飛沫是否會對熱帶氣旋的強(qiáng)度產(chǎn)生影響，但模擬結(jié)果表明，加入飛沫的影響，模式模擬出的熱帶氣旋發(fā)展演變過程更接近觀測，海氣溫差與觀測結(jié)果也更加一致。通過擬合分析HEXOS（humidity exchange over the sea）試驗(yàn)海氣界面資料，Andreas等［3—11］給出了一種適合于較強(qiáng)風(fēng)速情況的海洋飛沫通量參數(shù)化方法，該參數(shù)化方案為研究海洋飛沫在海氣間熱量和動量傳輸中的作用及其對海洋上風(fēng)暴的影響提供了新途徑。Andreas和Emanuel［6］重新考慮了飛沫液滴（spray droplet）的微物理過程，利用Andreas和DeCosmo［8］的參數(shù)化方法進(jìn)行研究指出海洋飛沫的蒸發(fā)使熱帶氣旋的最大潛在強(qiáng)度增強(qiáng)。Bao等［12］將FKH參數(shù)化方案應(yīng)用到一個海氣耦合模式，對颶風(fēng)個例Opal進(jìn)行模擬，發(fā)現(xiàn)如果只有一部分飛沫蒸發(fā)，蒸發(fā)所需的熱量來自海洋，這部分蒸發(fā)會使颶風(fēng)強(qiáng)度明顯增強(qiáng)，而如果飛沫全部蒸發(fā)，則蒸發(fā)會從大氣吸取熱量，從而使大氣邊界層低層冷卻，對颶風(fēng)增強(qiáng)產(chǎn)生負(fù)面影響。Wang等［13］使用其高分辨率的熱帶氣旋模式（TCM3）進(jìn)行理想實(shí)驗(yàn)，他們的結(jié)果表明飛沫蒸發(fā)使模擬熱帶氣旋增強(qiáng)率減小而且對其最終強(qiáng)度影響也很小。黎偉標(biāo)等［14］在中尺度模式中引入海洋飛沫參數(shù)化方案，可使臺風(fēng)期間海氣界面的潛熱通量增加50%，10 m層風(fēng)速最大值增加30%，從而使模擬臺風(fēng)的強(qiáng)度明顯增加，使模擬結(jié)果更趨于合理。

綜上所述，以上海洋飛沫對熱帶氣旋影響的研究大都針對大西洋，關(guān)于西太平洋和南海熱帶氣旋期間海洋飛沫對其演變的影響規(guī)律和機(jī)理，目前尚缺乏系統(tǒng)深入的研究，因此，在對臺風(fēng)的數(shù)值模擬中，考慮大氣和海洋飛沫的動力、熱力相互作用過程是必要的。本文使用的NOAA／ESRL海洋飛沫參數(shù)化方案版本是基于Fairall等［15］研究設(shè)計的，該參數(shù)化方案基本恰當(dāng)?shù)胤磻?yīng)了海洋飛沫對海氣交換過程的影響，把海洋飛沫參數(shù)化方案引入到中尺度大氣模式WRF中，選取2009年8月的熱帶氣旋“莫拉克”（Morakot）作為研究個例，分析海洋飛沫對熱帶氣旋邊界層結(jié)構(gòu)和強(qiáng)度的影響。

2 個例選取和數(shù)值試驗(yàn)

2.1 莫拉克臺風(fēng)過程簡介

2009年第8號熱帶氣旋Morakot于8月4日02時在菲律賓東北方海面生成，生命期為7 d。莫拉克臺風(fēng)形成不久后，向西北偏西方向移動，并在8月5日5時起轉(zhuǎn)為偏西行進(jìn)，朝臺灣東部沿?？拷?日20時增強(qiáng)為中度臺風(fēng)，7日23時50分左右在花蓮市附近登陸，于8日14時左右在桃園附近出海，進(jìn)入臺灣海峽北部，并持續(xù)向西北緩慢移動，9日18時30分左右在馬祖北方進(jìn)入福建，強(qiáng)度減弱，10日早晨進(jìn)入浙閩交界，然后穿過浙江，于8月11日減弱為熱帶性低氣壓，經(jīng)過太湖流域，11日下午在江蘇鹽城東臺入海。

2.2 模式方案

采用中尺度數(shù)值模式WRF（weather research and forecast），模擬試驗(yàn)選取0.5°分辨率GFS的AVN資料作為初始場，模擬采用雙重嵌套網(wǎng)格，區(qū)域中心點(diǎn)位于（27°N，132°E），水平分辨率由粗到細(xì)為9 km、3 km，垂直層數(shù)為35層，模式層頂50 hPa。粗網(wǎng)格的起始時間為2009年8月6日00時（世界時，下同），積分96 h，細(xì)網(wǎng)格啟動時間與粗網(wǎng)格一致，采用渦旋追隨嵌套網(wǎng)格（vortex following nested grid）技術(shù)，該技術(shù)最大優(yōu)點(diǎn)是熱帶氣旋（或臺風(fēng)）始終處在移動的最高分辨率網(wǎng)格區(qū)的中心，因此它的模擬區(qū)域是隨著模式運(yùn)行自適應(yīng)調(diào)整的（詳細(xì)模式配置見表1）。

表1 模式設(shè)計Tab.1 The model configuration

2.3 海洋飛沫參數(shù)化方案

在WRF模式中，洋面上的湍流交換過程是用Monin-Obukhov相似理論描述的。計算如下：

式中，κ是von Karman常數(shù)，取為0.4；z0m、z0t和z0q分別是空氣動力、溫度和濕度的粗糙度長度。

本文采用Fairall等［15］和Bao等［16］改進(jìn)的海洋飛沫參數(shù)化方案。該方案的特點(diǎn)是考慮了海洋飛沫的反饋影響，他們認(rèn)為，海氣之間總的動量、感熱和潛熱傳輸包括了湍流傳輸和海洋飛沫的傳輸，表達(dá)式為：

式中，τtot、Hs，tot、Hl，tot分別為動量、感熱和潛熱的總通量，α為次網(wǎng)格反饋常數(shù)，τsp、Qs、Ql分別為海洋飛沫引起的動量、感熱和潛熱通量，它們分別由下面的公式計算：

式中，F(xiàn)a、Fv是飛沫質(zhì)量流量，m是飛沫質(zhì)量，Vf是液滴的平均降落速度，S為海洋飛沫的濃度分布，Ts是海洋溫度，Ta是空氣溫度，Tw是飛沫液滴的濕球溫度，Wb＝3.8×10－6U3.4z是白浪的面積，qs（Ta）是海表上的大氣飽和混合比，qs（Tw）是在濕球溫度Tw時的飽和混合比。詳細(xì)計算見Bao等［16］。

2.4 試驗(yàn)設(shè)計

為了評估海洋飛沫對熱帶氣旋“莫拉克”的邊界層結(jié)構(gòu)和強(qiáng)度的影響，設(shè)計了比對試驗(yàn)（表2），即在邊界層過程中不考慮飛沫作用和考慮飛沫作用，分別稱之為控制試驗(yàn)（Control）和飛沫試驗(yàn)（Spray），飛沫試驗(yàn)中邊界層過程采用飛沫參數(shù)化方案。

表2 試驗(yàn)設(shè)計Tab.2 The experiment configure

3 數(shù)值模擬結(jié)果

3.1 對路徑和強(qiáng)度的影響

本文采用的是中國氣象局上海臺風(fēng)研究所最佳路徑和強(qiáng)度（引用CMA－STI的2009年資料），這套資料是中國氣象局上海臺風(fēng)研究所按照《臺風(fēng)業(yè)務(wù)和服務(wù)規(guī)定》，根據(jù)收集到的常規(guī)和非常規(guī)氣象觀測資料，對熱帶氣旋的路徑和強(qiáng)度資料進(jìn)行整編，是一套后分析產(chǎn)品，包括了業(yè)務(wù)當(dāng)時所獲取的附加信息，可認(rèn)為產(chǎn)品中熱帶氣旋位置和強(qiáng)度是相對最好的，即“最佳路徑和強(qiáng)度”。

圖1表示臺風(fēng)“莫拉克”從8月6日00時到10日00時的最佳路徑和強(qiáng)度以及Control試驗(yàn)和Spray試驗(yàn)的模擬路徑和強(qiáng)度。從圖1a可以看到模擬開始后，“莫拉克”主要向西北方向移動，兩試驗(yàn)較好地模擬了“莫拉克”的路徑，對比兩模擬試驗(yàn)和觀測發(fā)現(xiàn)，兩試驗(yàn)?zāi)M的路徑移動都稍微快于最佳路徑，這主要是由于模式嵌套網(wǎng)格在初始時刻（8月6日00時）兩個試驗(yàn)的臺風(fēng)中心均位于實(shí)況的西北面。在6日00時到7日00時之前的路徑模擬中，Control試驗(yàn)和Spray試驗(yàn)的路徑相似，和臺風(fēng)最佳移動路徑的誤差不是很大。在8日7日12時至8日00時，兩個試驗(yàn)?zāi)M的路徑都在最佳路徑的南側(cè)，這主要是由于在8月7日12時登陸臺灣，受陸地摩擦耗散的影響，其強(qiáng)度逐漸減弱，以較慢的速度穿過臺灣島，隨后沿西北向移動。在8日18時以后，Control試驗(yàn)和Spray試驗(yàn)?zāi)M的路徑相對最佳路徑移動都偏快。熱帶氣旋的移動路徑主要也受引導(dǎo)氣流、地形和不同尺度系統(tǒng)相互作用等因子的影響，而海洋飛沫主要是通過改變邊界層內(nèi)的熱量交換對熱帶氣旋產(chǎn)生影響，這可能是兩個試驗(yàn)?zāi)M的移動路徑較為相似的原因。

在對熱帶氣旋強(qiáng)度的模擬上，Control試驗(yàn)和Spray試驗(yàn)?zāi)M的臺風(fēng)系統(tǒng)均偏弱。圖1b給出了考慮和不考慮海洋飛沫作用時“莫拉克”最低海平面氣壓模擬結(jié)果的對比，由圖可見，在48 h之前，兩個試驗(yàn)?zāi)M的中心海平面氣壓都大于實(shí)況，這是由初始場中初始臺風(fēng)的強(qiáng)度比實(shí)際偏弱造成的。由于臺風(fēng)發(fā)生在洋面上，觀測資料的稀缺，往往造成臺風(fēng)的初始化存在很大的誤差。此后臺風(fēng)繼續(xù)發(fā)展，積分12 h后，加入飛沫參數(shù)化模擬結(jié)果強(qiáng)于不考慮飛沫作用模式結(jié)果，中心氣壓與實(shí)況值更為接近（8月8日00—18時）。圖1c給出了最大10 m風(fēng)速的模擬結(jié)果，在模擬前期，Control試驗(yàn)?zāi)M風(fēng)場較強(qiáng)，Spray試驗(yàn)?zāi)M風(fēng)場相對偏弱，是因?yàn)镾pray試驗(yàn)中，近表層很不穩(wěn)定，近表層的風(fēng)很強(qiáng)，高層的溫度梯度分布是穩(wěn)定的，這意味著低層與高層很少有動量交換，導(dǎo)致模式最低層風(fēng)速降低，同樣也影響10 m層風(fēng)速。

3.2 對邊界層結(jié)構(gòu)和強(qiáng)度的影響

本文分析主要選取6日18時至7日06時內(nèi)12 h平均，即熱帶氣旋加強(qiáng)階段。

圖2是兩個試驗(yàn)6日18時至7日06時內(nèi)12 h平均的邊界層徑向風(fēng)、切向風(fēng)軸對稱結(jié)構(gòu)。從圖中可以發(fā)現(xiàn)，模式模擬的熱帶氣旋的0.6 km邊界層高度里，邊界層內(nèi)有強(qiáng)的徑向風(fēng)，說明此處入流發(fā)展旺盛，有強(qiáng)烈的徑向入流輻合，徑向風(fēng)的最大水平梯度區(qū)位于“莫拉克”的墻眼區(qū)（半徑135 km處），值得注意的是，近表層的流速向內(nèi)逐漸增大，然后再逐漸減小，然后在風(fēng)眼中心快速的減為0，顯然地，導(dǎo)致徑向風(fēng)快速減速的原因是超梯度風(fēng)的存在。Control試驗(yàn)（見圖2a）中最大徑向入流風(fēng)速為8 m／s，Spray試驗(yàn)（見圖2b）中最大徑向入流風(fēng)速為9 m／s，圖2b表明當(dāng)加入海洋飛沫時，邊界層結(jié)構(gòu)發(fā)生了顯著的變化，徑向風(fēng)速顯著增加，徑向風(fēng)最大水平梯度也明顯增強(qiáng)，說明邊界層中的輻合作用有所加強(qiáng)；切向風(fēng)大值區(qū)略有減小，表明引入海洋飛沫時會增加表面層的熱量，使垂直混合增加，湍流混合增加，邊界層摩擦增加，從而導(dǎo)致切向風(fēng)的大值區(qū)略有較小。值得注意的是，邊界層結(jié)構(gòu)的變化會造成熱帶氣旋結(jié)構(gòu)與強(qiáng)度的變化。

圖1 2009年8月6日00：00—10日00：00臺風(fēng)“Morakot”最佳移動路徑、Control試驗(yàn)和Spray試驗(yàn)的路徑（a），中心最低氣壓（b）和近中心最大風(fēng)速（c）（臺風(fēng)位置每隔6h標(biāo)識）Fig.1 The best moving track，tracks of Control and Spray experiments（a），minimum sea level pressure（b）and maximum surface wind（c）for typhoon“Morako”from 6 to 10 August 2009

圖2 2009年8月6日18時至7日06時12 h平均的邊界層切向風(fēng)速（陰影，單位：m／s）、和徑向風(fēng)速（等值線，單位：m／s）軸對稱垂直結(jié)構(gòu)Fig.2 12 h averaged axisymmetric structure of tangential wind（shading，unit：m／s）and radial wind（solid line，unit：m／s）within boundary layer

感熱通量和潛熱通量的交換對臺風(fēng)的發(fā)生發(fā)展具有十分重要的作用。圖3是Control試驗(yàn)和Spray試驗(yàn)積分18 h后，潛熱和感熱通量的空間分布。從圖3可以看到，潛熱通量和感熱通量的主要分布形式大致相似，呈現(xiàn)高度的非對稱性，從外向中心逐漸減小，達(dá)到最大值后，向中心又逐漸減小，在熱帶氣旋中心熱通量數(shù)值較低，潛熱通量的數(shù)值均為正值，且量值比感熱通量大得多，說明海洋飛沫的蒸發(fā)為臺風(fēng)發(fā)展提供了大量的水汽和熱量，潛熱通量在臺風(fēng)的發(fā)展和發(fā)展中作用更大。這與Jan等［17］研究觀點(diǎn)一致：隨著風(fēng)速的增加，潛熱通量的增加比感熱通量明顯，潛熱通量對熱帶氣旋的發(fā)生和維持作用尤為重要。從圖3a和圖3b看到，Control試驗(yàn)的潛熱通量的最大值為600 W／m2，Spray試驗(yàn)最大值比Control試驗(yàn)大100 W／m2，且臺風(fēng)中心附近通量大于200 W／m2的區(qū)域Spray試驗(yàn)比Control試驗(yàn)大，兩試驗(yàn)的感熱通量的最大值都為80 W／m2（圖3c、d）。熱通量的低值區(qū)域在臺風(fēng)中心附近，在臺風(fēng)“莫拉克”中心的南側(cè)和西南側(cè)的熱通量都較大。在整個Spray試驗(yàn)?zāi)M時段內(nèi)，潛熱通量最大值為1 200 W／m2。從圖3b可以看到，感熱通量的大值區(qū)分布在臺風(fēng)中心西南側(cè)，而臺風(fēng)中心的北側(cè)為低值區(qū)，這是因?yàn)楦袩嵬康拇笮∪Q于海氣溫差和風(fēng)速大小。

飛沫的蒸發(fā)使邊界層的結(jié)構(gòu)發(fā)生了改變，更重要的是直接影響了感熱通量和潛熱通量的重新分配。圖4是Control試驗(yàn)和Spray試驗(yàn)從6日18時到7日06時12小時平均的軸對稱平均感熱通量、潛熱通量、總熱量隨半徑的分布。由圖4a感熱通量的分布可以看出，通量的大值區(qū)位于距熱帶氣旋中心130～220 km處，從熱帶氣旋中心到距中心約135 km內(nèi)，Spray試驗(yàn)和Control試驗(yàn)的感熱通量都隨著半徑逐漸增大，加入飛沫后，感熱通量增加了22%，距熱帶氣旋中心135 km以外，Spray模擬的感熱通量低于Control試驗(yàn)?zāi)M的感熱通量，這主要是由于飛沫的蒸發(fā)使得邊界層冷卻；Control試驗(yàn)中，感熱通量的量值約是潛熱通量的5%（圖4a、4b），只占總熱通量的6%（圖4a、4c）；而在130 km外，加入飛沫后，感熱通量減小。從圖4b可以看到，加入飛沫后，從熱帶氣旋中心到450 km處，潛熱通量明顯增加。由于飛沫的蒸發(fā)引起的潛熱通量的明顯增加，因此總的熱通量有所增加（圖4c），增加的大值區(qū)位于眼墻附近，這也很好的解釋了Spray試驗(yàn)?zāi)M的最大風(fēng)速和海平面氣壓強(qiáng)于Control試驗(yàn)?zāi)M的最大風(fēng)速和海平面氣壓，同樣，從圖4c中也可以看到，熱通量存在兩次大值，分別位于180 km和405 km的眼墻區(qū)域內(nèi)，這與圖2有很好的對應(yīng)關(guān)系。

圖3 2009年8月6日18時熱通量、海平面氣壓和10 m風(fēng)速的分布Fig.3 Distribution of heat flux，sea level pressure and 10 m wind speed at 18∶00 on 6 August 2009

圖5反映的是兩個試驗(yàn)軸對稱平均總熱通量（enthalpy flux）的徑向半徑與時間關(guān)系。從圖5a與圖5b可以看出，兩個試驗(yàn)的平均熱通量最大值出現(xiàn)在熱帶氣旋的眼墻區(qū)域，在3～9 h時平均熱通量出現(xiàn)極大值，Control試驗(yàn)最大超過440 W／m2，加入海洋飛沫后，熱通量最大超過480 W／m2，此時也是熱帶氣旋風(fēng)速相對較大的時段（見圖1c）。從圖5b可以看出，在整個模擬期間的大多數(shù)時間段內(nèi)，海洋飛沫對熱通量的影響可以達(dá)到或超過320 W／m2，分別在3～9 h和21～27 h這兩個時間段內(nèi)達(dá)到最大，因此，我們得到結(jié)論，加入海洋飛沫后，熱通量的增加十分明顯。

上述的研究已經(jīng)說明了海洋飛沫改變了熱帶氣旋的邊界層結(jié)構(gòu)，從而使熱帶氣旋的強(qiáng)度增強(qiáng)。圖6是Spray試驗(yàn)和Control試驗(yàn)從6日18時至7日06時12 h平均溫度、相對濕度、垂直風(fēng)速之差的軸對稱垂直結(jié)構(gòu)。從圖6a中可以看到，海洋飛沫的蒸發(fā)，使眼墻附近的溫度降低，最大降低0.3℃，臺風(fēng)中心到眼墻附近的濕度加大，從邊界層低層一直可延伸到1.5 km以上，很明顯，由于輻合作用，濕度增加的最大值集中在熱帶氣旋的中心0.6 km高度以上，因此在此高度以上利于云的形成。邊界層的冷卻有兩種潛在的影響，首先，它能增加對流層低層的靜態(tài)穩(wěn)定度，那樣會抑制眼墻的對流發(fā)展，很可能會使熱帶氣旋的強(qiáng)度減弱；其次，海氣之間溫度差異的增加會使地表層不穩(wěn)定，使近表層湍流加強(qiáng)，從而加強(qiáng)了眼墻的對流使熱帶氣旋的強(qiáng)度增強(qiáng)［18］。在這種環(huán)境下，來自海表和海洋飛沫的水汽通量被眼墻附近的上升氣流輸送到邊界層高處（見圖6b），從圖6b我們可以看到，加入飛沫后，眼墻的上升氣流增強(qiáng)。

圖4 2009年8月6日18時至7日06時12 h平均的軸對稱平均感熱通量（a）、潛熱通量、（b）和熱通量（c）隨半徑的變化Fig.4 Radial distribution of 12 h averaged axisymmetric sensible heat flux（a），latent heat flux（b）and the enthalpy flux（c）

圖5 軸對稱平均熱通量的時間-半徑演變圖（單位：W／m2）Fig.5 Radial-time Hovmoller diagram of axisymmetric structure of the averaged surface enthalpy flux（W／m2）

圖6 試驗(yàn)Control與Spray試驗(yàn)6日18時至7日06時12 h平均的（a）溫度（等值線，單位：℃）、相對濕度（陰影，單位：%）、（b）垂直風(fēng)速（單位：m／s）軸對稱垂直結(jié)構(gòu)差異（Spray－Control）Fig.6 The difference 12 h averaged axisymmetric structure of（a）temperature（contour line，unit：℃）and relative humility（shaded，unit：%），and（b）vertical wind speed（unit：m／s）.（Spray minus Control）

海洋飛沫的蒸發(fā)導(dǎo)致潛熱通量的增加，從而導(dǎo)致試驗(yàn)Spray比試驗(yàn)Control在熱帶氣旋眼墻附近更多的降水（見圖7）。從圖中可以看出，降水隨半徑的變化與熱通量的分布相似，從氣旋中心到225 km，隨著半徑的增大，降水量逐漸增大，在225 km處，降水量達(dá)到最大，隨后隨著半徑的增加，降水量減少。很明顯，考慮海洋飛沫作用后，降水明顯增強(qiáng)，尤其在眼壁附近。

4 結(jié)論

將海洋飛沫參數(shù)化引入到一個高分辨率、非靜力中尺度WRF模式中，我們的數(shù)值模擬結(jié)果清楚的表明：

（1）采用海洋飛沫參數(shù)化后，模式對熱帶氣旋的強(qiáng)度預(yù)報有改進(jìn)，對熱帶氣旋移動路徑可能影響不大。

（2）引入海洋飛沫參數(shù)化后，熱帶氣旋邊界層內(nèi)眼墻區(qū)域的平均徑向風(fēng)速、切向風(fēng)速、均增強(qiáng)，表明飛沫對熱帶氣旋邊界層結(jié)構(gòu)變化的影響顯著。

（3）由于飛沫的蒸發(fā)，增加了表面熱通量，特別是潛熱通量，從而氣旋邊界層內(nèi)的垂直混合增大，使得熱帶氣旋的直接湍流增強(qiáng)，使得邊界層內(nèi)眼墻區(qū)域的垂直風(fēng)速、溫度、相對濕度、熱通量，降水等物理量均有增強(qiáng)，表明飛沫對熱帶氣旋強(qiáng)度變化的影響顯著。

圖7 2009年8月6日18時至7日06時12 h平均的軸對稱平均累積降水（單位：mm）隨半徑的變化Fig.7 Radial distribution of 12 h averaged axisymmetric of accumulated precipitation（unit：mm）

最后，需要注意的是，低風(fēng)速情況下的飛沫參數(shù)化的不確定性，還需要更多的研究和觀測事實(shí)來進(jìn)一步探討；其次，海洋飛沫對熱帶氣旋發(fā)展和演變的影響機(jī)制，以及海洋飛沫對臺風(fēng)影響的定量描述都將作為今后更進(jìn)一步的研究。

［1］Anthes R A.Tropical cyclones：Their evolution，structure and effects［J］.Meteor Monogr，1982，41：208.

［2］Fairall C W，Kepert J D，Holland G J.The effect of sea spray on the surface energy transports over the ocean［J］.The Global Atmos and Ocean Sys，1994，2：121—142.

［3］Andreas E L.Sea spray and the turbulent air-sea heat fluxes［J］.J Geophys Res，1992，97（C7）：11429—11441.

［4］Andreas E L.The temperature of evaporating sea spray droplets［J］.J Atmos Sci，1995，52：852—862.

［5］Andreas E L.A new sea spray generation for wind speed up to 32m／s［J］.J Phys Oceanogr，1998，28：2175—2184.

［6］Andreas E L，Emanuel K A.Effects of sea spray on tropical cyclone intensity［C］／／Preprints，23rdConf on Hurricanes and Tropical Meteorology.Dallas，TX，Amer Meteor Soc，1999：22—25.

［7］Andreas E L，Emanuel K A.Effects of sea spray on tropical cyclone intensity［J］.J Atmos Sci，2001，58：3741—3751.

［8］Andreas E L，DeCosmo J.The signature of sea spray in the HEXOS turbulent heat flux data［J］.Bound-Layer Meteor，2002，103：303—333.

［9］Andreas E L.An algorithm to predict the turbulent air-sea fluxes in high-wind，spray conditions［C］／／Preprints，12thConf On Interaction of the Sea and Atmosphere.Long Beach，CA，9－13 February，Amer Meteor Soc，2003，CD－ROM3.4，7.

［10］Andreas E L.Spray stress revisited［J］.J Phys Oceanogr，2004，34：1429—1440.

［11］Andreas E L.A bulk air-sea flux algorithm for high-wind，spray conditions，version 2.0［C］／／Preprints，13thConf On Interaction of the Sea and Atmosphere.Long Beach，CA，Amer Meteor Soc，2004，1.5.

［12］Bao J W，Wilczak J M，Choi J K，et al.Numerical simulation of air-sea interaction under high wind conditions using a coupled model：a study of hurricane development［J］.Mon Wea Rev，2000，128：2190—2209.

［13］Wang Y，Kepert J D，Holland G J.The effect of sea spray evaporation on tropical cyclone boundary layer structure and intensity［J］.Mon Wea Rev，2001，129：2481—2500.

［14］黎偉標(biāo)，何溪澄，唐潔.臺風(fēng)“森拉克”的數(shù)值模擬研究：海洋飛沫的作用［J］.熱帶海洋學(xué)報，2004，23（3）：58—65.

［15］Fairall C W，Bradley E E，Hare J E，et al.Bulk parameterization of sir-sea fluxes：Updates and verification for the COARE algorithm［J］.J Climate，2003，16：571—591.

［16］Bao J W，F(xiàn)airall C W，Michelson S A.Evaluation and improvement of spray-modified air-sea enthalpy and momentum flux parameterizations for operational hurricane prediction［R］.The First Semiannual Report of Project，2008.

［17］Jan E Meirink，Vladimir K Makin.The impact of sea spray evaporation in a numerical weather prediction model［J］.J Atmos Sci，2001，58（23）：3626—3638.

［18］Liu T，Katsaros K B，Businger J A.Bulk parameterization of sir-sea exchange of heat and water vapor including the molecular constants at the interface［J］.J Atmos Sci，1979，101：20771—20797.

［19］Guan Hao，Zhou Lin，Wang Hanjie，et al.A mesoscale atmosphere-ocean wave coupling model and numerical simulations on a strong typhoon process in South China Sea［J］.Acta Meteorologica Sinica，2008，66（3）：342—350.

［20］Liu Y，Zhang D L，Yau M K.A multiscale numerical study of hurricane Andrew（1992）PartⅠ.Explicit simulation and verification［J］.Mon Wea Rev，1997，125：169—176.

［21］Jeffrey S Gall，William M Frank.Effects of sea spray on tropical cyclones simulated under idealized conditions［J］.Mon Wea Rev，2008，136：1686—1705.

［22］Perrie W，Andreas E L，Zhang W，et al.Sea spray impacts on intensifying midlatitude cyclones［J］.J Atmos Sci，2005，62（6）：1867—1883.

［23］Uang C L.Impact of sea spray and oceanic response on the development of tropical cyclones［C］／／Preprints，23rdConf on Hurricanes and Tropical Meteorology.Dallas，TX，Amer Meteor Soc，1999：30—31.

［24］Zeng Z H，Chen L S，Wang Y，et al.A numerical simulation study of super Typhoon Saomei（2006）intensity and structure changes［J］.Acta Meteorologica Sinica，2009，67（5）：750—763.

［25］Zhang W Q，Perrie W.Impacts of waves and sea spray on midlatitude storm structure and intensity［J］.Mon Wea Rev，2006，134：2418—2442.

［26］Zhang W Q，Perrie W.The influence of air-sea roughness，sea spray，and storm translation speed on waves in North Atlantic storms［J］.J Phys Oceanogr，2008，38（4）：817—839.

［27］鄭靜，費(fèi)建方，王元，等.海洋飛沫對熱帶氣旋影響的數(shù)值試驗(yàn)［J］.熱帶氣象學(xué)報，2008，24（5）：467—474.

Effect of sea spray on tropical cyclone boundary layer structure

Wang Ping1，2，Chen Baode1，2，Zeng Zhihua1

（1.Shanghai Typhoon Institute of China Meteorological Administration，Shanghai 200030，China；2.Key Laboratory of Numerical Modeling for Tropical Cyclone，China Meteorological Administration，Shanghai 200030，China）

The sea spray flux parameterization scheme is incorporated into a high resolution and non-hydrostatic weather research and forecasting（WRF）model to simulate tropical cyclone（TC）Morakot（0908），and to study the effect of sea spray on boundary layer structure and intensity of Morakot.The results show that the intensity prediction of Morakot has been improved by using the new scheme，while there is no any impact on its track forecasting.Moreover，the averaged elements，such as radial wind，tangential wind，temperature，relative humidity，vertical wind，heat flux and precipitation in TCinner core region are increase by improving TC boundary layer processes.That means the contributions of all the elements have play important roles on Morakot structure and intensity change.

sea spray；tropical cyclones；morakot；boundary layer；parameterization scheme

P732.3

A

0253-4193（2014）09-0084-10

王平，陳葆德，曾智華.海洋飛沫對熱帶氣旋邊界層結(jié)構(gòu)的影響［J］.海洋學(xué)報，2014，36（9）：84—93，

10.3969／j.issn.0253-4193.2014.09.010

Wang Ping，Chen Baode，Zeng Zhihua.Effect of sea spray on tropical cyclone boundary layer structure［J］.Acta Oceanologica Sinica（in Chinese），2014，36（9）：84—93，doi：10.3969／j.issn.0253-4193.2014.09.010

2014-02-12；

2014-07-08。

國家自然科學(xué)基金項(xiàng)目——?dú)馊苣z-云微物理-降水相互作用對登陸我國臺風(fēng)的影響（41175094）；熱帶氣旋強(qiáng)度變化的環(huán)境因子影響研究（41275067）；近海熱帶氣旋邊界層過程對其結(jié)構(gòu)和強(qiáng)度變化的影響（40875039）；基于ADAS－WRF的RUC系統(tǒng)引進(jìn)及本地化改進(jìn)技術(shù)研究（GYHY201306014－2）。

王平（1986—），女，陜西省綏德縣人，主要從事數(shù)值模擬研究。E-mail：wangp＠m(xù)ail.typhoon.gov.cn