周冠博，林青，高守亭，聶高臻

(1. 國家氣象中心，北京100081；2. 福建省氣象臺(tái)，福建 福州350001 3. 中國科學(xué)院大氣物理研究所云降水物理與強(qiáng)風(fēng)暴實(shí)驗(yàn)室，北京100029)

1 引 言

臺(tái)風(fēng)暴雨的成因非常復(fù)雜，影響因素眾多。中央氣象臺(tái)已有的格點(diǎn)化降水預(yù)報(bào)產(chǎn)品雖然能夠提供較高時(shí)空分辨率的降水預(yù)報(bào)產(chǎn)品，但其并沒有考慮到臺(tái)風(fēng)特有的非對(duì)稱結(jié)構(gòu)演變、中尺度螺旋雨帶發(fā)展、地形作用、客觀降水預(yù)報(bào)訂正及臺(tái)風(fēng)路徑預(yù)報(bào)誤差等諸多影響因素，因此臺(tái)風(fēng)暴雨的預(yù)報(bào)準(zhǔn)確率還有待進(jìn)一步提高，預(yù)報(bào)產(chǎn)品也不足以滿足實(shí)際的預(yù)報(bào)服務(wù)的需求。對(duì)以往的若干臺(tái)風(fēng)個(gè)例進(jìn)行分析發(fā)現(xiàn)，臺(tái)風(fēng)中的斜壓渦度變化與臺(tái)風(fēng)暴雨的落區(qū)密切相關(guān)，臺(tái)風(fēng)暴雨雨帶的移動(dòng)往往落后于斜壓渦度的變化，并沿著斜壓渦度擾動(dòng)的變化趨勢(shì)方向移動(dòng)，因此提出一個(gè)問題：斜壓渦度變化的異常值是否會(huì)對(duì)實(shí)際臺(tái)風(fēng)暴雨的預(yù)報(bào)具有一定的參考價(jià)值？而以往在這方面的探討與研究還很少，因此本文的研究內(nèi)容是臺(tái)風(fēng)中的斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究。

簡要回顧可知，很多的氣象學(xué)家雖然已經(jīng)對(duì)登陸臺(tái)風(fēng)做了大量的研究工作，但是對(duì)臺(tái)風(fēng)登陸時(shí)很多過程認(rèn)識(shí)不足，這就導(dǎo)致我們對(duì)登陸臺(tái)風(fēng)強(qiáng)降水的預(yù)報(bào)水平不高。例如，許多專家學(xué)者引入了各種形式的Ertel 位渦[1-4]。針對(duì)位渦在赤道地區(qū)應(yīng)用時(shí)存在一定的局限性，Gao 等[5-6]又提出了兩維對(duì)流渦度矢量的概念，可以有效地表征赤道地區(qū)深對(duì)流系統(tǒng)的發(fā)展演變。

在流體力學(xué)中渦度表示流體質(zhì)塊的旋轉(zhuǎn)程度和旋轉(zhuǎn)方向，進(jìn)入20 世紀(jì)人們認(rèn)識(shí)到渦度場(chǎng)對(duì)臺(tái)風(fēng)的影響很大。Ertel[1]提出斜壓大氣中的位渦，又稱Ertel 位渦，它是絕對(duì)渦度矢量與位溫梯度矢量的點(diǎn)積，在絕熱無摩擦大氣中位渦具有守恒性，是一個(gè)既包含熱力因子又包含動(dòng)力因子的物理量[7]。Hoskins 等[8]系統(tǒng)地分析了位渦的性質(zhì)及其可反演性，并提出了等熵位渦的概念，使得位渦理論得到進(jìn)一步發(fā)展。然而，在伴有潛熱釋放等非絕熱過程的大氣中，位渦失去了保守性。為了使其在濕絕熱過程中仍具有保守性，Bennetts 等[9]利用相當(dāng)位溫替換位溫，定義了相當(dāng)位渦（下稱濕位渦），相當(dāng)位溫的引入，拓展了位渦作為診斷量的應(yīng)用范圍[10-11]，使其被廣泛應(yīng)用到熱帶氣旋、溫帶氣旋、颮線、鋒面雨帶等系統(tǒng)的診斷分析中[12-16]。經(jīng)典的濕位渦主要應(yīng)用于臺(tái)風(fēng)暴雨的診斷分析當(dāng)中[17-18]，并揭示了對(duì)暴雨落區(qū)的指示意義[19]，但是實(shí)際大氣并不滿足處處飽和，為了體現(xiàn)實(shí)際大氣非均勻飽和的狀態(tài)，Gao 等[20]定義了廣義的濕位渦，并將其應(yīng)用于暴雨系統(tǒng)的診斷分析當(dāng)中，取得了一些有意義的結(jié)果[21-25]。

Dutton[26]利用拉格朗日坐標(biāo)系將單個(gè)氣團(tuán)的渦度分解為正壓渦度和斜壓渦度。Zhang[27]在Dutton 的基礎(chǔ)上進(jìn)一步推導(dǎo)出拉格朗日中的次級(jí)渦度概念。斜壓渦度在臺(tái)風(fēng)研究中很少運(yùn)用到。

風(fēng)場(chǎng)可以分解為輻合和輻散兩個(gè)物理量表示有旋部分和無旋部分，渦度從等熵面與垂直速度的夾角不同也可以分解為PV（位渦）以及CVV（對(duì)流渦度矢量）。渦度從大氣的正壓性以及斜壓性也能分解出正壓渦度和斜壓渦度兩個(gè)物理量。

雖然臺(tái)風(fēng)常被用對(duì)稱結(jié)構(gòu)研究，但從實(shí)際大氣狀況來看，其非對(duì)稱更明顯，所以有必要研究斜壓性渦度在臺(tái)風(fēng)中的表現(xiàn)。

2 斜壓渦度的公式推導(dǎo)

在笛卡爾坐標(biāo)下的渦度方程中，其方程右側(cè)會(huì)出現(xiàn)力管項(xiàng)，這個(gè)力管項(xiàng)通常被稱為是大氣斜壓項(xiàng)，也是產(chǎn)生渦度的源項(xiàng)?？梢姕u度本身應(yīng)含有兩部分，即正壓部分和斜壓部分。那么這兩部分是否可以分解出來呢？在過去很少人關(guān)注這一問題，也沒有多少學(xué)者去討論斜壓渦度在哪些方面有用處，為了探討斜壓渦度在臺(tái)風(fēng)及其暴雨中的表現(xiàn)及應(yīng)用，本節(jié)給出斜壓渦度的數(shù)學(xué)表達(dá)式，以便定量計(jì)算。

在歐拉坐標(biāo)系中，通常梯度的數(shù)學(xué)表達(dá)式為：

而在拉格朗日坐標(biāo)系中，其梯度表達(dá)式又可寫為：

這里I=1，2，3，i=1，2，3，分別表示拉格朗日坐標(biāo)系下的分量及歐拉坐標(biāo)系下的分量符號(hào)。

在拉格朗日坐標(biāo)系下有：

在濕大氣中，熵的表達(dá)式可寫為：

其中C 是常數(shù)。

于是有：

由式（5）可以看出氣壓梯度可表示為：

其中H=CpTe是焓。

于是運(yùn)動(dòng)方程可以寫為：

其中? 為重力位勢(shì)，忽略了摩擦力。

由式（3）可知：

則對(duì)式（8）進(jìn)行拉格朗日坐標(biāo)系下的時(shí)間積分后有：

這里利用等熵過程，即在短時(shí)間內(nèi)熵是不變的。

對(duì)方程（10）兩邊求旋度，則有：

從式（11）中可以看出，科氏力產(chǎn)生的渦度與大氣中的斜壓性無關(guān)是正壓部分，而其斜壓渦度為：

下面我們便可以利用觀測(cè)資料或數(shù)值模式模擬的資料來具體計(jì)算斜壓渦度了。

3 高分辨率數(shù)值模擬

本文選取2009年第8 號(hào)臺(tái)風(fēng)“莫拉克”(Morakot) 作為個(gè)例進(jìn)行數(shù)值模擬，“莫拉克”于2009年8月4日凌晨（北京時(shí)，下同）生成，5日加強(qiáng)為臺(tái)風(fēng)，7日23:45 在臺(tái)灣省花蓮登陸，9日16:20 在福建省霞浦再次登陸，9日晚上在福建省境內(nèi)減弱為強(qiáng)熱帶風(fēng)暴，10日凌晨減弱為熱帶風(fēng)暴，12日02：00 停止編號(hào)。從臺(tái)風(fēng)“莫拉克”生成到結(jié)束9 天時(shí)間里給我國多省(市)帶來嚴(yán)重創(chuàng)傷, 其中臺(tái)灣受創(chuàng)最為嚴(yán)重。近年來，許多專家和學(xué)者都對(duì)“莫拉克”產(chǎn)生的降水過程做了詳細(xì)的模擬和分析[28-32]。

林青等[31]利用WRF v3.2.1 模式對(duì)臺(tái)風(fēng)“莫拉克”的精細(xì)結(jié)構(gòu)開展了高分辨率的數(shù)值模擬，模擬時(shí)間范圍為2009年8月2日12:00—10日12:00，物理過程采用Goddard 方案和Kain-Fritsch 積云對(duì)流參數(shù)化方案，同時(shí)采用YSU 邊界層方案、Dudhia 短波輻射和RRTM 長波輻射方案。3 個(gè)模擬區(qū)域水平分辨率分別為18 km、6 km、2 km，模式積分區(qū)域中心為120.5 °E，26.0 °N，格點(diǎn)數(shù)分別為520×520、378×378、520×492。內(nèi)二層使用移動(dòng)嵌套。模式垂直方向取41 層。模擬結(jié)果對(duì)臺(tái)風(fēng)的路徑以及強(qiáng)度的模擬與觀測(cè)總體接近，同時(shí)能夠模擬出臺(tái)風(fēng)“莫拉克”在登陸內(nèi)陸期間的精細(xì)化風(fēng)雨結(jié)構(gòu)，本文應(yīng)用該模擬資料進(jìn)行斜壓渦度的計(jì)算與分析。

圖1 是2009年8月9日08：00 和14：00 的逐6 小時(shí)降水量的實(shí)況和模擬的對(duì)比，基本上模擬的雨帶分布和實(shí)況一致，沿閩浙沿海呈現(xiàn)東北-西南走向，而強(qiáng)降水區(qū)域主要在福鼎一帶，整體模擬雨量高于實(shí)況雨量，且位于福建中部沿海即莆田、泉州一帶的降水量級(jí)略偏小，但整體看來，該模擬結(jié)果能夠再現(xiàn)出臺(tái)風(fēng)登陸閩浙過程的降水變化精細(xì)過程，其結(jié)果比較可靠，可以在此基礎(chǔ)上開展進(jìn)一步的診斷分析。

圖1 2009年8月9日08：00（a、b）和14：00（c、d）6 小時(shí)觀測(cè)累積降水（a、c）和模擬累積降水（b、d）

4 計(jì)算結(jié)果與分析

下面采用以上高分辨率的數(shù)值模擬結(jié)果剖析斜壓渦度的時(shí)空分布特點(diǎn)。圖2 和圖3 是斜壓渦度和MPV（濕位渦，）對(duì)比，從圖中可以看出在臺(tái)風(fēng)“莫拉克”登陸前，很明顯低層的斜壓渦度有更強(qiáng)的異常信號(hào)，圍繞臺(tái)風(fēng)內(nèi)核呈現(xiàn)出正負(fù)正的位相特點(diǎn)，登陸后也是斜壓渦度的異常值信號(hào)更強(qiáng)于MPV。

由圖2 可以看到，在臺(tái)風(fēng)“莫拉克”登陸臺(tái)灣島以前，斜壓渦度異常值區(qū)主要位于臺(tái)風(fēng)眼區(qū)和外圍西北部象限的大部分地區(qū)，其分布形勢(shì)與對(duì)流的分布狀況相近。而濕位渦的異常值區(qū)主要位于臺(tái)風(fēng)中心，其濕位渦的異常值區(qū)分布比較分散，與臺(tái)風(fēng)結(jié)構(gòu)的分布情況不能很好地對(duì)應(yīng)。

（二）培養(yǎng)學(xué)生數(shù)學(xué)學(xué)習(xí)中的創(chuàng)新思維。創(chuàng)新性思維是發(fā)散性思維的核心。發(fā)散性思維是指擴(kuò)散思維或求異思維，是指大腦在思維時(shí)呈現(xiàn)出的一種擴(kuò)散模式，具體表現(xiàn)為思維視野的廣闊，在數(shù)學(xué)教學(xué)中的體現(xiàn)就是“一題多解”。因此，培養(yǎng)學(xué)生的發(fā)散性思維不僅有利于數(shù)學(xué)知識(shí)的學(xué)習(xí)，還有利于其他方面的創(chuàng)造力。在數(shù)學(xué)教學(xué)中教師要巧妙地利用“問題”“障礙”“趣味”等方式引入以激發(fā)起學(xué)生潛在的強(qiáng)烈求知欲。另外，因?yàn)樾W(xué)生的身心特點(diǎn)，他們的思維往往是以形象思維為主的，故而，教師應(yīng)多方位、多角度地培養(yǎng)學(xué)生的抽象思維。

由圖3 可以看出，在臺(tái)風(fēng)“莫拉克”登陸臺(tái)灣島以后，斜壓渦度異常值區(qū)主要位于臺(tái)風(fēng)眼區(qū)以及其中心的南側(cè)，眼區(qū)內(nèi)呈現(xiàn)正負(fù)渦度替換的趨勢(shì)，其分布形勢(shì)與對(duì)流的分布狀況相近。

圖3 2009年8月8日00：00700 hPa (a、b)、850 hPa (c、d)濕位渦(a、c)和斜壓渦度(b、d)的平面分布

濕位渦在臺(tái)風(fēng)中心的信號(hào)比較弱，而斜壓渦度較強(qiáng)，這是由于斜壓渦度為，而濕位渦中的斜壓部分為，即因?yàn)棣哩尅摹哩孲·▽?duì)萫=0，所以濕位渦中的斜壓部分為0，沒有這部分信息。

從沿著臺(tái)風(fēng)中心時(shí)間剖面圖4a～4c 可以看出，登陸前斜壓渦度低層多為負(fù)正負(fù)的位相，而且隨著時(shí)間的推移，斜壓渦度有從大氣的高層向臺(tái)風(fēng)的移動(dòng)中心傳遞的趨勢(shì)，即在臺(tái)風(fēng)即將到達(dá)時(shí)原先的正渦度被替換為負(fù)的渦度。所以斜壓渦度對(duì)其移動(dòng)有一定指示意義。

由圖5 可以看出，在臺(tái)風(fēng)“莫拉克”過臺(tái)灣島時(shí)，會(huì)以次中心誘生過島，斜壓渦度體現(xiàn)為負(fù)渦度消失，在山地附近有正渦度生成，完成過島，臺(tái)風(fēng)中心被替換?？梢钥闯鲂眽簻u度對(duì)于臺(tái)風(fēng)以次生中心過島的移動(dòng)同樣具有一定的指示意義。

圖4 2009年8月斜壓渦度沿著臺(tái)風(fēng)中心的時(shí)間剖面分布a. 7日12：00 的臺(tái)風(fēng)中心位于122.2 °E，24.0 °N；b. 7日21：00 的臺(tái)風(fēng)中心位于121.4 °E，24.0 °N；c. 8日02：00 的臺(tái)風(fēng)中心位于120.8 °E，24.0 °N。

圖5 2009年8月7日23：00 700 hPa 斜壓渦度的平面分布（a）、斜壓渦度沿著臺(tái)風(fēng)中心121.9 °E，23.8 °N（b）和斜壓渦度沿著臺(tái)風(fēng)中心121.4 °E，23.7 °N（c）的時(shí)間剖面分布

大氣運(yùn)動(dòng)是三維的，對(duì)于導(dǎo)致臺(tái)風(fēng)暴雨的中尺度系統(tǒng)除了平行于位溫梯度方向的渦度分量，還要考慮垂直于位溫梯度方向的渦度分量。而且垂直于位溫梯度方向的渦度分量對(duì)于中尺度系統(tǒng)的發(fā)生發(fā)展有十分重要的影響。為了不遺漏大氣運(yùn)動(dòng)的動(dòng)力信息，考慮了渦度的三維分量，在位渦的基礎(chǔ)上又引入一個(gè)新的物理量——斜壓渦度。如圖6a、6b 所示，斜壓渦度的異常值主要位于中低層大氣時(shí)，一般會(huì)產(chǎn)生較強(qiáng)的臺(tái)風(fēng)降水。

圖6 2009年8月9日20：00（a）和2009年8月10日02：00（b）斜壓渦度沿著119 °E 的剖面 其中直方圖代表降水。

5 臺(tái)風(fēng)“海棠”與“鳳凰”的個(gè)例分析

下面采用1 °×1 °的NCEP/FNL 粗網(wǎng)格再分析資料初步分析斜壓渦度在2005年臺(tái)風(fēng)“海棠”和2008年臺(tái)風(fēng)“鳳凰”中的表現(xiàn)。

從圖7a 可以看到，在臺(tái)風(fēng)“海棠”過臺(tái)灣島前，斜壓渦度的異常大值區(qū)主要位于臺(tái)風(fēng)眼區(qū)和外圍眼墻的大部分地區(qū)，其分布形勢(shì)與對(duì)流的分布狀況非常相近。從圖7b 可以看出，登陸后斜壓渦度在眼區(qū)內(nèi)呈現(xiàn)出正負(fù)渦度替換的趨勢(shì)，其分布形勢(shì)與對(duì)流的分布狀況也非常相近。從剖面圖7c 和7d 上看，其中心自高到低都是斜壓渦度異常區(qū)，正負(fù)渦度相間，登陸前后的位相明顯不同。

從圖8a～8d 可以看出，斜壓渦度在臺(tái)風(fēng)“鳳凰”個(gè)例中也有類似的表現(xiàn)，由此可以看出斜壓渦度可以作為臺(tái)風(fēng)追蹤的一個(gè)示蹤量，并且斜壓渦度與臺(tái)風(fēng)內(nèi)部的對(duì)流分布也有一定關(guān)系，其正負(fù)位相變化與臺(tái)風(fēng)的結(jié)構(gòu)變化及登陸與否都存在一定的聯(lián)系。

圖7 斜壓渦度（850 hPa）2005年8月17日14：00(a、c)及18日20：00(b、d)的平面分布(a、b)和過臺(tái)風(fēng)中心的剖面分布(c、d)

圖8 斜壓渦度（850 hPa）2008年8月28日14：00(a、c)及29日02：00(b、d)的平面分布(a、b)和過臺(tái)風(fēng)中心的剖面分布(c、d)

6 結(jié) 論

（1）斜壓渦度和MPV（濕位渦）對(duì)比，可以看出在登陸前和登陸后，明顯低層斜壓渦度有更強(qiáng)的異常信號(hào)，圍繞臺(tái)風(fēng)內(nèi)核呈現(xiàn)正負(fù)正的位相特點(diǎn)。

（2）從沿著臺(tái)風(fēng)中心時(shí)間剖面可以看出，登陸前斜壓渦度低層多為負(fù)正負(fù)的位相，而且隨時(shí)間有從上層向臺(tái)風(fēng)的移動(dòng)中心傳遞的趨勢(shì)，即在臺(tái)風(fēng)即將到達(dá)時(shí)原先的正渦度被替換為負(fù)的渦度。所以對(duì)其移動(dòng)有一定指示意義。

（3）登陸后，在臺(tái)風(fēng)“莫拉克”過臺(tái)灣島時(shí)，以次中心誘生過島，其斜壓渦度體現(xiàn)為負(fù)渦度消失，在山地附近正渦度生成，完成過島，臺(tái)風(fēng)中心被替換。

（4）為了不遺漏大氣運(yùn)動(dòng)的動(dòng)力信息，考慮了渦度的三維分量，在位渦的基礎(chǔ)上又引入一個(gè)新的物理量——斜壓渦度。當(dāng)斜壓渦度的異常值主要位于中低層大氣時(shí)，一般會(huì)產(chǎn)生較強(qiáng)的臺(tái)風(fēng)降水。

目前在該領(lǐng)域的研究工作值得繼續(xù)深入進(jìn)行，今后將在本文已取得的研究成果基礎(chǔ)上，繼續(xù)開展“斜壓渦度”在不同臺(tái)風(fēng)個(gè)例中的應(yīng)用研究工作，深入研究臺(tái)風(fēng)中的斜壓渦度的變化與臺(tái)風(fēng)暴雨關(guān)系的普適性，并希望在臺(tái)風(fēng)暴雨的實(shí)際預(yù)報(bào)工作中取得更好的預(yù)報(bào)效果，同時(shí)提供更科學(xué)的理論依據(jù)。