華東沿海地形對登陸熱帶氣旋運動影響的理想數值研究

陳小宇 吳立廣 吳玉琴

摘要利用中尺度數值模式設計一組高分辨率理想試驗，采用位渦趨勢方法定量診斷分析熱帶氣旋在登陸我國華東沿海地形時，其運動發(fā)生的精細化變化以及不同因子的貢獻。結果表明，平地的存在使得登陸熱帶氣旋移速相對更快，當華東沿海地形存在時，熱帶氣旋移速顯著增大，這種增速現象主要是由于平地和地形所引起的非對稱氣流以及相應的引導氣流變化所致，這很可能是導致預報路徑誤差的一個重要原因。平地試驗中，陸地在熱帶氣旋低層激發(fā)出中小尺度的非對稱氣流，與之不同的是，實際地形的加入激發(fā)出更大尺度并且更強的非對稱偏南氣流。位渦趨勢方法的診斷結果表明，非引導效應總體而言對熱帶氣旋運動貢獻較小，這是因為這些因子相互抵消，但在不同的垂直層次上，不同的非引導因子貢獻存在明顯的差異。

關鍵詞登陸熱帶氣旋運動；地形；非對稱氣流；位渦趨勢方法

熱帶氣旋（TC）的路徑和運動一直以來都是熱帶氣旋預報和研究中最為重要的方面之一，路徑預報的準確度直接關系著強度、降水以及災害等方面的預報結果。早期的觀點大多認為熱帶氣旋的運動主要受大尺度環(huán)境氣流的引導效應所影響（Adem and Lezama，1960；George and Gray，1976；Neumanm，1979；Chan and Gray，1982）。但許多學者在研究中發(fā)現，觀測中熱帶氣旋的運動方向和環(huán)境氣流的方向之間始終存在一定的偏差，并且這種偏差并不是由于引導氣流計算方式和范圍的差別所引起的（George and Gray，1976；DeMaria，1985；Chan and Williams，1987；Fiorino and Elsberry，1989）。早在20世紀80年代初，Holland（1983）就提出了β漂移（βdrift）理論，對這一現象的機制進行了解釋，發(fā)現行星渦度的經向梯度在熱帶氣旋內部激發(fā)出一對渦旋對，渦旋對之間的向極向西的通風流是導致熱帶氣旋運動與環(huán)境引導氣流存在偏差的原因。這一理論得到后續(xù)大量的觀測和數值研究的驗證（Holland，1983；Chan and Williams，1987；Fiorino and Elsberry，1989；Willoughby，1990，1992；Peng and Williams，1990；Li and Wang，1994；Wang et al.，1997）。

除了環(huán)境導引和β漂移之外，Wu and Wang（2000）發(fā)現熱帶氣旋運動還受到非引導氣流效應的影響。他們提出了一種新的診斷方法，把熱帶氣旋看作具有較強正位勢渦度（PV）異常的系統(tǒng)，因此將熱帶氣旋運動與位勢渦度的趨勢變化聯(lián)系起來，這種方法被稱為位渦趨勢（PVT）方法。該方法不僅考慮傳統(tǒng)的引導氣流和β漂移的作用，還包含了垂直運動、非絕熱加熱以及摩擦等非引導效應對于熱帶氣旋運動的貢獻。Wu and Wang（2001a，2001b）利用這一方法進一步分析了絕熱加熱和非絕熱加熱情形下，引導和非引導效應對于熱帶氣旋運動的貢獻。Chan et al.（2002）采用該方法診斷觀測中熱帶氣旋路徑變化時，發(fā)現不同的因子對于各種類型的熱帶氣旋路徑貢獻有所差別。此外，一些研究在分析熱帶氣旋登陸的理想試驗和個例模擬中利用該診斷方法區(qū)別不同物理過程的作用（Wong and Chan，2006；吳玉琴等，2015）。這些結果表明，位渦趨勢方法相比于傳統(tǒng)的引導氣流方法，可以更好地應用于熱帶氣旋運動的診斷分析中。

登陸熱帶氣旋給我國沿海地區(qū)的人民生命財產安全造成了巨大的損失（Zhang et al.，2009；張嬌艷等，2011；曹劍等，2012；李肖雅等，2014），因此對于登陸熱帶氣旋路徑的研究直接關系著災害損失和防臺減災工作的展開。Tuleya and Kurihara（1978）最早設計了理想試驗以研究登陸熱帶氣旋的特征變化。在這之后，許多研究針對島嶼和山脈地形（例如臺灣島）進行了大量的觀測和數值研究（Chang，1982；Bender et al.，1987；Yeh and Elsberry，1993a，1993b；Lin et al.，2005；Jian and Wu，2008；Yang et al.，2008；Huang et al.，2011；Wang et al.，2012；Xie and Zhang，2012）。與之相對應的是，較少有研究涉及大陸地形對熱帶氣旋路徑的影響。針對這一問題，Wong and Chan（2006）設計了幾組理想試驗，研究在沒有背景氣流和地形的影響下，熱帶氣旋登陸f平面的平坦陸地過程中運動的變化。他們發(fā)現陸地的存在激發(fā)熱帶氣旋環(huán)流產生不對稱的結構，這種不對稱性逐漸發(fā)展，最終使得熱帶氣旋產生向陸地的運動趨勢，運動速度約為1 m/s。這些結果表明，陸地使得熱帶氣旋的不對稱結構發(fā)展從而可能導致登陸過程中移速增加。后續(xù)的數值研究進一步驗證，在β平面和三角洲地形影響下這種運動加速現象的存在（AuYeung and Chan，2010；Szeto and Chan，2010），但是這種研究中同樣沒有考慮環(huán)境氣流的影響。針對這一問題，Li et al.（2013）在試驗中加入了均勻的環(huán)境氣流，發(fā)現了類似的現象，但是他們在研究中著重關注熱帶氣旋登陸過程中降水的變化，缺乏對于運動變化的定量性分析。在熱帶氣旋登陸大陸地形的過程中，較少有研究考慮環(huán)境氣流的影響，尤其當環(huán)境氣流和地形高度同時存在時，會對登陸熱帶氣旋運動產生何種影響，仍然缺乏更深入的認識。

近三十年間，登陸我國的熱帶氣旋陸上持續(xù)時間和移動速度都出現了明顯的變化（Chen et al.，2011），一些研究結果也表明，登陸我國大陸與海南島的熱帶氣旋造成了重大的直接經濟損失和人員傷亡（Zhang et al.，2009；張嬌艷等，2011），但是對于登陸我國大陸沿海地區(qū)的熱帶氣旋的運動研究仍然偏少。因此，本文設計高分辨率的理想試驗模擬熱帶氣旋登陸我國華東沿海地區(qū)，通過對比試驗結果定量地診斷背景氣流和華東地區(qū)實際地形對于登陸熱帶氣旋運動的影響，以及不同的物理過程對于運動變化的貢獻。

1資料與試驗設計

使用WRFARW V 221模式，試驗設計采用了五層網格雙向嵌套設置，最外層網格為水平分辨率27 km，格點數為230×210，網格中心點為1325°E、300°N。四層內網格的水平分辨率分別為9、3、1、1/3 km，格點數分別為432×399、333×333、501×501以及720×720。模式的垂直層次為40層，模式頂氣壓為50 hPa。最外兩層網格為固定網格，其范圍設置足夠包含了西北太平洋海域的主要系統(tǒng)的發(fā)展和演變；最內三層網格為移動網格，伴隨著模擬的熱帶氣旋運動而移動，主要覆蓋了逐個時次上的熱帶氣旋影響范圍。

選取2005年典型熱帶氣旋“麥莎”活動并登陸的2005年8月5—9日這一時段，利用Lanczos時間濾波方法（Duchon，1979）得出該時段的20 d以上低頻背景場作為試驗的背景場（圖1）。8月5號00時（世界時，下同），“麥莎”位于臺灣東北側洋面之上，此時該熱帶氣旋的近地面風速約為45 m·s-1，在后續(xù)的3 d中，“麥莎”沿著環(huán)境場氣流向西北方向移動并于8月5日19時40分在我國東南沿海地區(qū)登陸，登陸之后強度開始下降并逐漸消散。初始時刻，該低頻背景場呈現出典型的西北太平洋夏季的環(huán)流形式，主要表現為較為明顯的季風槽和副高的分布特征。需要指出的是由于環(huán)境場已經過濾波處理（已去除20 d以下頻段的系統(tǒng)），因此后續(xù)登陸試驗中的TC路徑與強度均與臺風“麥莎”存在明顯區(qū)別。

最外層網格的微物理參數化方案和積云對流方案分別采用WSM3簡單冰方案和KainFritropical cycloneh方案（Dudhia，1989；Kain and Fritch，1993），第二至第五層網格采用WSM6微物理方案并關閉積云對流方案（Hong and Lim，2006）。邊界層過程選用的是Yonsei University 邊界層方案（Noh et al.，2003），輻射過程選擇了Dudhia 短波方案（Dudhia，1989）、Rapid Radiative Transfer Model（RRTM）長波方案（Mlawer et al.，1997），陸面過程采用的是Noah陸面過程模式。模式初始化選用美國國家環(huán)境預報中心（National Centers for Environmental Prediction，NCEP）發(fā)布的分辨率為1°×1°的FNL全球分析資料（Final Operational Global Analysis），該資料時間間隔為6 h。試驗中對最外層網格進行張弛逼近（Nudging），以便積分過程的背景場更為穩(wěn)定。最外層網格的積分時間步長為90 s，總積分時間為96 h。

在該低頻背景場中，放入人工構造的理想軸對稱渦旋。為了進一步改進人造渦旋的結構和風壓關系，本研究中，預先將構造的理想軸對稱渦旋放入一個f平面無環(huán)境場的全海洋試驗中，讓渦旋在該環(huán)境中進行初旋（Spinup），使其結構和風壓關系得到模式動力與熱力過程的調整和改善。進行初旋之前，構造的理想渦旋的海平面最低氣壓為994 hPa，近地面最大風速約30 m·s-1。經過12 h的初旋后，該渦旋結構依然保持對稱，海平面最低氣壓下降至965 hPa，近地面最大風速達到40 m·s-1，已達到了臺風強度（圖2）。之后從初旋試驗結果中計算得出該渦旋的風場、溫度場、濕度場以及氣壓場等物理變量，加入圖1所示低頻背景場之中與熱帶氣旋“麥莎”（8月5日00時）相近的位置（1226°E，251°N），使其沿著季風槽與副高之間的東南氣流向西北方向運動，之后登陸我國華東沿海地區(qū)。

在低頻背景場和理想渦旋的基礎上，設計一組理想試驗，以此對比分析在高分辨率的情況下，熱帶氣旋在登陸過程中運動的精細化變化。三個試驗分別為全海洋試驗（OCN），平地試驗（FLT）和地形試驗（TOP）。全海洋試驗中，將整個區(qū)域的下墊面和表層屬性全部設置成海洋屬性；平地試驗中，去除所有的地形高度，將其全部設置為0；而在地形試驗中，包括真實的地形高度和分布。同時，在所有試驗中，海平面溫度統(tǒng)一設置為29 ℃，在平地和地形試驗中，地表溫度也同樣設置為29 ℃，這么設置是為了排除海平面和地表溫度的差異所可能造成的熱帶氣旋登陸過程中的差異。同時，在平地和地形試驗中，陸表屬性也設置為均一的類型，本研究中選擇的陸表類型是我國華中華北地區(qū)常見的灌溉型農田與牧場（irrigated cropland and pasture），該類型的地表摩擦層厚度為007 m。

2模擬結果

圖3給出了模擬的海洋試驗、平地試驗以及地形試驗中熱帶氣旋的積分前48 h移動路徑，其中熱帶氣旋中心定義為最大方位角平均的切向風所對應的中心（Wu et al.，2006）。為了剔除地形對低層氣旋環(huán)流的影響，三組實驗中將700 hPa氣旋中心的移動軌跡作為氣旋整體的移動路徑。從圖3可以看出，三組試驗中熱帶氣旋的運動方向一致但是移速存在區(qū)別，其中海洋試驗、平地實驗、地形試驗中熱帶氣旋分別在積分42 h、39 h、24 h到達海岸線，這表明在平地試驗、地形試驗中熱帶氣旋的移動速度較海洋試驗中更快。需要指出的是，在海洋試驗中去除了陸地和地形，因此在該試驗的分析中將熱帶氣旋移動至原陸地位置的時刻定義為該試驗的登陸時刻。

從強度模擬結果（圖4）看，平地實驗、地形試驗中熱帶氣旋登陸后強度迅速減弱，48 h后，兩個試驗中熱帶氣旋近地面風速降為熱帶風暴級別。在海洋試驗中，渦旋在積分前12 h強度增加，并維持在40 m/s。與海洋試驗類似，平地試驗和地形試驗中氣旋在海洋中移動強度增加，在登陸前3 h氣旋強度開始減弱。熱帶氣旋登陸時，平地試驗、地形試驗中熱帶氣旋強度為310、313 m/s，登陸9 h后兩個試驗中強度減弱為194、176 m/s。積分48 h后，三個試驗中熱帶氣旋強度分別為：383、191、148 m/s。海平面氣壓場的變化表明，在登陸前24 h，三個試驗中熱帶氣旋的海平面最低氣壓接近，但這之后，海洋試驗熱帶氣旋繼續(xù)增強，氣壓持續(xù)下降。與之相反的是，平地和地形試驗中的熱帶氣旋在登陸前氣壓即開始上升，登陸后12 h，海洋試驗、平地試驗和地形試驗中的熱帶氣旋海平面最低氣壓分別為950、976、989 hPa，同樣表明了地形試驗中熱帶氣旋強度下降的速率更大。

吳玉琴等（2015）利用WRF模式模擬熱帶氣旋“莫拉克”（2008），通過考慮熱帶氣旋在移動過程中自身的PV一波發(fā)展改進了原有的PVT診斷方法，并通過對比PVT估算速度、大尺度環(huán)境引導氣流與模擬的熱帶氣旋移速，發(fā)現改進后的PVT估算速度較引導氣流能更好地表征熱帶氣旋的運動，特別是在熱帶氣旋運動出現振蕩時刻。因此，本文利用改進后的PVT方法診斷理想試驗結果，重點分析熱帶氣旋登陸過程中引導效應與非對稱結構變化之間的關系。

為了評估大尺度環(huán)境引導氣流與PVT估算速度在表征熱帶氣旋移動上的差異，在此分別計算了海洋、平地、地形數值試驗中熱帶氣旋移速、PVT估算速度、引導氣流的對比（圖5）。參考前人計算方法（Wu and Wang，2000；吳玉琴等，2015），計算PVT各項估算速度和引導氣流時采用第二層網格（D02）的固定坐標系數據，熱帶氣旋移速采用第五層高分辨率網格中的氣旋中心計算得出。此外，三個數值試驗中各垂直層上PVT估算速度均為以熱帶氣旋中心為中心，直徑為540 km的圓形區(qū)域計算，垂直方向上采用850～300 hPa多層平均，各層引導氣流采用相同的計算范圍和垂直層次所得，熱帶氣旋移速為前后兩小時位置差計算所得。

32熱帶氣旋移速的差異

三個數值試驗中，熱帶氣旋的移速出現明顯的差異。海洋、平地和地形試驗中氣旋在登陸前后的37 h平均移速分別為267、313、462 m/s，同時段平均引導氣流大小分別為251、271、423 m/s。這表明在平地試驗中即出現了熱帶氣旋登陸增速的現象，與Wong and Chan （2006）結果類似，但僅約05 m/s的增速幅度明顯小于他們所發(fā)現的1 m/s的增幅，這可能是由于本研究中所采用的地表摩擦層厚度僅為007 m，小于Wong and Chan （2006）的試驗設置。更值得注意的是，在地形試驗中，熱帶氣旋移速相對于海洋試驗中的結果增加了約2 m/s，表明實際地形的存在使得熱帶氣旋在登陸過程中移速顯著增大，并且這種現象很可能是由于引導氣流的增大所導致。楊瓊瓊和吳立廣（2015）在分析臺風路徑預報誤差時發(fā)現大多數的預報誤差來源于臺風移速的預報偏差而非移動方向的偏差，因此在環(huán)境氣流和地形共同影響下，熱帶氣旋登陸時運動所出現的變化很可能是導致路徑預報誤差的一個重要原因。

分析試驗中熱帶氣旋移速的經向和緯向分量后發(fā)現，PVT估算速度較引導氣流更接近熱帶氣旋的移速（圖5），特別是在熱帶氣旋移動過程中出現速度突變情況下，例如海洋試驗和平地試驗中熱帶氣旋登陸前后6 h經向和緯向速度均出現短時振蕩特征（圖3a—d），PVT估算速度很好地表征了熱帶氣旋移動中的速度變化特征，而傳統(tǒng)的引導氣流計算難以表征這一現象。熱帶氣旋登陸前24 h到登陸后12 h內，平地試驗中平均經向、緯向速度分別為-123、286 m/s，平均PVT估算速度為-116、262 m/s，相應的平均引導氣流為-057、263 m/s。PVT估算速度與熱帶氣旋移速之間均方根誤差分別為013、029 m/s，引導氣流與熱帶氣旋移速之間均方根誤差為066、04 m/s，從計算的平均移速、均方根誤差數值上都可以看出PVT估算速度相較于引導氣流能更好地表征熱帶氣旋的移動，并且在海洋試驗和地形試驗中存在類似的結果。這是由于PVT診斷方法不僅包含了大尺度環(huán)境引導氣流對熱帶氣旋運動的作用，還包含了熱帶氣旋結構變化、非絕熱加熱、PV垂直平流等對熱帶氣旋運動的作用。需特別指出的是，在地形試驗中的熱帶氣旋經向速度分析中，PVT速度與熱帶氣旋移速的均方根誤差略大于引導氣流的偏差，這主要是由于該試驗中熱帶氣旋強度的迅速下降所導致的，但PVT方法仍然可以更好地再現熱帶氣旋在登陸前后的緯向速度以及經向速度的振蕩。結合上述結果，可利用PVT診斷方法診斷三個數值試驗中熱帶氣旋運動的差異。

Wu and Wang（2000）指出局地PV一波變化由PV的水平平流、PV垂直平流、非絕熱加熱、摩擦四項組成，為了評估各物理過程對熱帶氣旋運動的作用，在此分別計算三組數值試驗中PV的水平平流（簡稱HA）、PV垂直平流（簡稱VA）、非絕熱加熱（簡稱DH）、摩擦（簡稱FR）對應的估算速度與PVT估算速度的對比（圖6）。

從圖6可以看出，三組數值試驗中PV水平平流對應的估算速度最接近PVT估算速度，在海洋數值試驗中，PVT、HA、VA、DH、FR各項所對應的緯向估算速度為-087、-021、-07、005、-001 m/s，經向速度為242、261、038、-043、-014 m/s，相應的引導氣流為-052、243 m/s。從數值上可以看出，HA對應的估算速度接近引導氣流的大小，從長時間平均來看，HA對熱帶氣旋運動的作用要遠大于其他三項的作用總和，此外，VA與DH對熱帶氣旋運動的作用是相互抵消的，這也導致了熱帶氣旋的移動基本由HA決定。與海洋試驗相同，平地試驗與地形試驗中HA對PVT的貢獻最大，但HA與引導氣流之間的差異相比海洋試驗增加，這可能是兩個試驗登陸過程中由熱帶氣旋非對稱結構的發(fā)展引起的非引導效應所導致的。總體而言，PVT方法診斷的結果表明，HA項對于三個試驗熱帶氣旋的運動起著最重要的貢獻，同時可以表征三個熱帶氣旋移速的明顯差異（圖6）。Wu and Wang （2001a）指出，HA項的貢獻主要表示了熱帶氣旋中的引導效應，但與傳統(tǒng)意義的引導氣流有所區(qū)別，不僅包括了大尺度環(huán)境引導氣流和β漂移的作用，還包含了熱帶氣旋內垂直切變、垂直運動、絕熱和非絕熱加熱等過程所造成的非對稱氣流對熱帶氣旋運動的直接貢獻。

33登陸熱帶氣旋的非對稱結構

為了進一步分析產生登陸熱帶氣旋運動變化的原因，分別計算三組試驗中熱帶氣旋登陸時刻降水率與700 hPa的風場進行對比（圖7）。從圖7可以看出，海洋試驗中熱帶氣旋登陸時刻降水率與700 hPa風場基本呈對稱結構，降水主要集中在距離熱帶氣旋中心120 km區(qū)域內，且熱帶氣旋眼墻北側區(qū)域內降水最大（圖7a，b）。平地實驗中，登陸時刻熱帶氣旋的風場、降水速率結構上非對稱性增強，大值區(qū)均在氣旋東北側（圖7c，d），地形試驗中風場、降水分布與平地試驗存在相似結構，但強度上地形試驗的非對稱性要略弱于平地試驗。三個試驗中的降水大值均分布在順切變左側區(qū)域，與之前的研究結果一致，表明強對流更傾向于發(fā)生在順切變左側區(qū)域（Frank and Ritchie，2001；Chen and Yau，2003）。

為了進一步了解地形對熱帶氣旋非對稱結構的影響，分別計算平地試驗和海洋試驗中不同時間段內中低層熱帶氣旋的一波非對稱風場差值以及降水率的差值（圖8a—d）。為了更深入地研究登陸前熱帶氣旋結構的變化，在此將熱帶氣旋登陸前的24 h分為兩個時間段：第一個為登陸前24～13 h（b1時間段），第二個為登陸前12 h（b2時間段），下同。從圖中可以看出，平地試驗中b1時間段內，700 hPa一波非對稱風場熱帶氣旋中心處有弱西北氣流，同時在熱帶氣旋西北側（陸地的區(qū)域）出現了氣旋式的風場異常，平地試驗和海洋試驗中的降水速率分布差異是這種低層地形導致輻合氣流的結果。熱帶氣旋內核區(qū)域的西南向的通風流與兩個試驗中水平速度的差異相當。b2時間段內，低層非對稱氣流的增強與地形試驗中熱帶氣旋在登陸之前的加速度是一致的，700 hPa上出現了較為明顯的南風氣流，氣流左側的氣旋性環(huán)流異常與陸地的分布相對應（圖6d）。但是這種非對稱結構在中層難以被識別出，500 hPa上的南風非對稱氣流特征明顯偏弱，且在熱帶氣旋內核區(qū)域出現了北風異常，與熱帶氣旋的移速增速趨勢相反。從中低層對比可以看出，熱帶氣旋在平地試驗的登陸過程中的增速現象主要來源于700 hPa上非對稱南風氣流的貢獻，說明平地試驗中地形所激發(fā)的非對稱結構主要表現在低層。

在地形試驗中，b1時間段內，相較于海洋試驗和平地試驗，地形導致的非對稱氣流強度更強；非對稱的南風氣流在低層到高層上都比平地試驗中的結果更顯著（圖9a—c）。另一個顯著特征是，非對稱氣流的半徑超過熱帶氣旋的半徑，達到約1 000 km半徑的范圍，雨帶的分布與平地試驗類似。在b2時間段內，700 hPa上非對稱氣流最顯著，這可能是因為熱帶氣旋低層環(huán)流受地形的影響最為顯著。結合上述分析：在地形試驗中，非對稱南風氣流與平地試驗中的結果相比明顯增強，非對稱南風氣流不僅出現在中低層，在高層也有明顯的分布特征，同時氣流強度和覆蓋范圍也明顯增大，因此導致地形試驗中熱帶氣旋的移速最快。此外值得一提的是，非對稱氣流的半徑明顯大于熱帶氣旋的尺度，這一結果表明，地形不僅導致了熱帶氣旋出現了非對稱結構，還顯著地影響了環(huán)境風場。這一現象體現了地形與平地對于登陸熱帶氣旋運動的不同影響程度，平地只能對中低層的熱帶氣旋環(huán)流產生影響，而地形可以對熱帶氣旋和環(huán)境氣流同時造成非常重要的影響，使得氣旋在登陸過程中產生明顯的增速現象。

34非引導效應的作用

吳玉琴等（2015）曾指出非引導效應對于臺風“莫拉克”的移速振蕩有著非常重要的作用。那么在登陸過程中，其余非引導效應起著何種作用？為什么熱帶氣旋運動主要受到引導氣流的影響？針對這些問題，為了進一步定量地分析非引導效應的作用，分別計算三個試驗中熱帶氣旋在b1時間段內、b2時間段內各非引導效應的平均貢獻，結果如圖10所示。

Wu and Wang（2001a）進一步將PV水平平流項分為兩部分：

∧1（-V·P）≈-V1·Ps-Vs·P1。

其中：腳標s和1表示對稱分量和一波分量；等式右端第一項為非對稱氣流引起的對稱PV分量的輸送，簡稱ASAF，包含大尺度引導氣流和通風流對熱帶氣旋運動的影響，與一波風場的區(qū)域平均引導氣流對應較好；第二項為對稱氣流引起的非對稱PV分量的平流，簡稱AASF，包含了β渦對、垂直切變及非對稱加熱產生的非對稱氣流等物理過程對熱帶氣旋運動的作用。本研究中將引導氣流以外的各物理過程統(tǒng)稱為非引導效應，包括了AASF項，PV垂直輸送項（VA），非絕熱加熱項（DH）以及摩擦項（FR）。

圖10海洋試驗（a，d）、平地試驗（b，e）、地形試驗（c，f）中登陸前24～13 h（a，b，c）和12～1 h（d，e，f）平均的多層平均位渦趨勢各非引導項估算移速（VA項估算速度，藍色；DH項估算速度，紅色；FR項目估算速度，灰色；AASF項估算速度，紫色；上述各項的矢量和，棕色）

Fig.10Averaged nonsteering terms including vertical advection（VA，blue），diabatic heating（DH，red），friction（FR，gray），advection of the asymmetric PV structure by the symmetric flow（AASF，purple） and vector sum of nonsteering terms（brown） during （a，b，c）24—13 h and （d，e，f）12—1 h before landfall in （a，d）OCN，（b，e）FLT and （c，f）TOP experiments

由圖10可以看出，三組試驗中，b1、b2時間段內，去除引導氣流后PVT各分項對熱帶氣旋運動的作用不一致，PVT各分量的矢量和小于引導氣流作用。海洋試驗中b1、b2時間段內非引導效應總貢獻的矢量和分別為052、022 m/s，平地試驗中兩個時間段矢量和分別為069、022 m/s（圖10中棕色矢量），而在兩組試驗中熱帶氣旋移動的平均速度大約3 m/s，因此海洋試驗、平地試驗中登陸時段內非引導效應對熱帶氣旋運動的貢獻遠小于引導氣流的貢獻。在地形試驗中，b1、b2時間段內非引導效應的總貢獻分別為116、138 m/s，明顯大于海洋和平地試驗的結果，但在該試驗中熱帶氣旋移速明顯增大，因此非引導效應對于移速的整體貢獻依然偏小。值得注意的是，由圖10可以看出：VA、DH項各自的估算速度數值上均大于或者接近非引導項的總貢獻，但矢量方向上存在一定程度的反位相分布特征，這說明在登陸過程中，非引導效應各項之間存在一定程度相互抵消，因此對熱帶氣旋運動的作用整體貢獻偏小，熱帶氣旋在登陸過程中總體仍受到引導效應的作用。此外，非引導效應導致熱帶氣旋有向西運動的趨勢，同時在地形試驗中，由于地形的存在使得熱帶氣旋登陸過程中FR項加大。

圖10的結果表明垂直方向上多層平均的非引導效應對熱帶氣旋運動的總體作用較小，但在各垂直層上是否有同樣的作用呢？Wu and Wang（2001b）在數值研究中發(fā)現，加熱場在熱帶氣旋高低層呈現相反的作用。為了進一步分析各物理因子在不同垂直層次上的不同貢獻，分別診斷了三個試驗中，熱帶氣旋登陸前12 h內高中低層PVT各分項對熱帶氣旋運動的貢獻（圖11）。從圖11可以看出：熱帶氣旋登陸前12 h內，三個試驗中PVT方程各項的高低層貢獻存在較大差異，例如HA項和ASAF項高低層（700 hPa與300 hPa）貢獻相反，DH項與ASAF項在同一層次貢獻相反，這些發(fā)現與Wu and Wang（2001b）的結論一致。需要指出的是，海洋和平地試驗中ASAF項在中層貢獻較小，而在高層對熱帶氣旋運動的作用均為負貢獻，但在地形試驗的ASAF項中高層均為較為顯著的正貢獻，這一結果與該試驗中引導氣流的顯著增大對應。此外，兩個登陸試驗的低層DH項與海洋試驗結果存在明顯差異，這說明熱帶氣旋在登陸過程中受陸地和地形的影響，低層非絕熱加熱以及相應的非對稱對流、降水等結構發(fā)生了相應的變化，例如圖7中的降水率差別所示。以上結果表明，雖然非引導效應對于熱帶氣旋運動的整體貢獻相對于引導效應偏小，但在不同的垂直層次上各項的作用存在明顯差異，尤其在登陸試驗中，DH項發(fā)生了明顯的變化。

圖11海洋試驗（a，b，c）、平地試驗（d，e，f）、地形試驗（h，i，j）中登陸前12～1 h平均的300 hPa（a，d，h）、500 hPa（b，e，i）和700 hPa（c，f，j）位渦趨勢各項估算移速（VA項估算速度，藍色；DH項估算速度，紅色；FR項目估算速度，灰色；AASF項估算速度，紫色；ASAF估算速度，黑色）

Fig.11Averaged PVT terms including vertical advection（VA，blue），diabatic heating（DH，red），friction（FR，gray），advection of the asymmetric PV by the symmetric flow（AASF，purple） and advection of the symmetric PV by the asymmetric flow（ASAF，black） at （a，d，h）300 hPa，（b，e，i）500 hPa and （c，f，j）700 hPa during 12—1h before landfall in （a，b，c）OCN，（d，e，f）FLT and （h，i，j）TOP experiments

4結論與討論

對于熱帶氣旋登陸過程中的運動變化一直以來都是備受關注的一個問題，但是長期以來，更多的研究關注于特殊的島嶼地形（例如臺灣島）對于路徑的影響，僅有少量研究涉及大陸地形的作用（Wong and Chan，2006；AuYeung and Chan，2010；Szeto and Chan，2010；Li et al.，2013）。然而這些研究中，通常又沒有考慮地形高度和環(huán)境氣流的影響。因此，本研究針對這一問題，設計了一組高分辨率的數值理想試驗，通過試驗結果對比以及位渦趨勢方法的診斷結果，分析研究登陸我國華東沿海地區(qū)的熱帶氣旋在登陸過程中的運動變化，以及不同的物理因子對于這種變化的不同貢獻。

分析結果表明：當典型的低頻背景場存在時，在熱帶氣旋登陸前24 h至登陸后12 h這一段時間內，熱帶氣旋在平地試驗中的移速相比在海洋試驗中的移速增加了約05 m/s，在具有實際地形高度的試驗中，熱帶氣旋移速存在約2 m/s的增幅。進一步的分析發(fā)現，這些移速增加的現象主要是由地形所引起的非對稱氣流以及相應的引導效應變化所引起的。在熱帶氣旋登陸前，平地下墊面可以激發(fā)出中低層的非對稱氣流，這種熱帶氣旋尺度的非對稱氣流導致了引導氣流的增大，從而產生移速的加快。當地形高度存在時，地形可以激發(fā)出更大尺度的非對稱南風氣流，同時這些非對稱氣流更強，從熱帶氣旋的低層發(fā)展到高層，因此造成了地形試驗中顯著增大的熱帶氣旋移速。需要指出的是，這種地形所造成的非對稱氣流同時包括了地形對于熱帶氣旋環(huán)流和對于環(huán)境氣流的分別貢獻，如何定量化地區(qū)別兩者的不同貢獻，有待今后的研究中進一步的分析。

除了引導效應之外，非引導效應對于登陸過程中的熱帶氣旋運動的整體貢獻相對較小，但是其作用不可忽視。位渦趨勢方法的診斷結果表明，這些非引導因子整體存在相互抵消的現象，因此對熱帶氣旋運動總貢獻偏小，這與前人研究結果一致。此外，在不同的垂直層次上，不同的非引導效應的貢獻存在著明顯的差異，這些差異與熱帶氣旋登陸過程中強度以及結構等方面的變化有著密切聯(lián)系，表明熱帶氣旋強度和結構的變化同樣可能在一定程度會對其運動產生影響。

綜合以上結論，可以看出我國大陸華東沿海地區(qū)的地形可以使得熱帶氣旋在登陸該區(qū)域時移速發(fā)生顯著的增大，增幅明顯大于前人在平地試驗中發(fā)現的結果，這種現象在前人研究中未有涉及，可能是導致臺風路徑預報誤差的一個原因，值得更深入的研究和探討。

參考文獻（References）

Adem J，Lezama P，1960.On the motion of cyclone embedded in a uniform flow[J].Tellus，12：255258.

AuYeung A Y M，Chan J C L，2010.The effect of a river delta and coastal roughness variation on a landfalling tropical cyclone[J].J Geophys Res，115（D19）.

Bender M A，Tuleya R E，Kurihara Y，1987.A numerical study of the island terrain on tropical cyclones[J].Mon Wea Rev，115：130155.

曹劍，吳立廣，潘維玉，20122006年7—9月西北太平洋熱帶氣旋季節(jié)活動的數值模擬[J].大氣科學學報，35（2）：148162.Cao J，Wu L G，Pan W Y，2012.Simulated seasonal activity of tropical cyclone over the western North Pacific during JulySeptember 2006[J].Trans Atmos Sci，35（2）：148162.（in Chinese）.

Chan J C L，Gray W M，1982.Tropical cyclone motion and surrounding flow relationship[J].Mon Wea Rev，110：13541374.

Chan J C L，Williams R T，1987.Analytical and numerical studies of betaeffect in tropical cyclone motion.Part Ⅰ：zero mean flow[J].J Atmos Sci，44：12571265.

Chan J C L，Ko F M F，Lei Y M，2002.Relationship between potential vorticty tendency and tropical cyclone motion[J].J Atmos Sci，59：13171336.

Chang S W J，1982.The orographic effects induced by an island mountain range on propagating tropical cyclones[J].Mon Wea Rev，110：12551268.

Chen X，Wu L，Zhang J，2011.Increasing duration of tropical cyclones over China[J].Geophys Res Lett，38：L02708.doi：101029/2010GL046137.

Chen Y，Yau M K，2003.Asymmetric structures in a simulated landfalling hurricane[J].J Atmos Sci，60（18）：22942312.

DeMaria M，1985.Tropical cyclone motion in a nondivergent barotropic model[J].Mon Wea Rev，113：11991210.

Duchon C E，1979.Lanczos filtering in one and two dimensions[J].J Appl Meteor，18：10161022.

Dudhia J，1989.Numerical study of convection observed during the winter monsoon experiment using a mesoscale twodimensional model[J].J Atmos Sci，46：30773107.

Fiorino M J，Elsberry R L，1989.Some aspects of vortex structure related to tropical cyclone motion[J].J Atmos Sci，46：975990.

Frank W M，Ritchie E A，2001.Effects of vertical wind shear on the intensity and structure of numerically simulated hurricanes[J].Mon Wea Rev，129：22492269.

George J E，Gray W M，1976.Tropical cyclone motion and surrounding parameter relationships[J].J Appl Meteor，15：12521264.

Holland G J，1983.Tropical cyclone motion：environmental interaction plus a beta effect[J].J Atmos Sci，40：328342.

Hong S Y，Lim J O J，2006.The WRF singlemoment 6class microphysics scheme （WSM6）[J].J Korean Meteor Soc，42：129151.

Huang Y H，Wu C C，Wang Y，2011.The influence of island topography on typhoon track deflection[J].Mon Wea Rev，139：17081727.

Jian G J，Wu C C，2008.A numerical study of the track deflection of super typhoon Haitang （2005） prior to landfall in Taiwan[J].Mon Wea Rev，136：598615.

Kain J S，Fritch J M，1993.Convective parameterization for mesoscale models：the KainFritch scheme.The representation of cumulus convection in numerical models Meteor Monogr[J].Amer Meteor Soc，46：165170.

Li X，Wang B，1994.Acceleration of the hurricane beta drift by shear strain rate of an environmental flow[J].J Atmos Sci，53：327334.

李肖雅，吳立廣，宗慧君，2014.季風渦旋影響西北太平洋臺風生成初步分析[J].大氣科學學報，37（5）：653664.Li X Y，Wu L G，Zong H J，2014.Analysis of influence of monsoon gyres on tropical cyclogenesis over the western North Pacific[J].Trans Atmos Sci，37（5）：653664.（in Chinese）.

Li Y，Cheung K K W，Chan J C L，2013.Numerical study on the development of asymmetric convection and vertical wind shear during tropical cyclone landfall[J].Quart J Roy Meteor Soc，140：18661877.

Lin Y L，Chen S Y，Hill C M，et al.，2005.Control parameters for the influence of a mesoscale mountain range on cyclone track continuity and deflection[J].J Atmos Sci，62：18491866.

Mlawer E J，Taubman S J，Brown P D，1997.Radiative transfer for inhomogeneous atmosphere：RRTM，a validated correlatedk model for the longwave[J].J Geophys Res，102（D14）：1666316682.

Neumann C J，1979.A guide to Atlantic and Eastern Pacific models for the prediction of tropical cyclone motion[C]//NOAA Tech Memo NWS NHC11.

Noh Y，Cheon W G，Hong S Y，2003.Improvement of the Kprofile model for the planetary boundary layer based on largeeddy simulation data[J].Bound Layer Meteor，107：401427.

Peng M S，William R T，1990.Dynamics of vortex asymmetries and their influence on vortex motion on a betaplane[J].J Atmos Sci，47：19872003.

Szeto K C，Chan J C L，2010.Structural changes of a tropical cyclone during landfall：βplane simulations[J].Adv Atmos Sci，27：11431150.

Tuleya R E，Kurihara Y，1978.A numerical simulation of the landfall of tropical cyclones[J].J Atmos Sci，35：242257.

Wang B，Li X，Wu L，1997.Hurricane beta drift direction in horizontally sheared flows[J].J Atmos Sci，54：14621471.

Wang C C，Kuo H C，Chen Y H，et al.，2012.Effects of asymmetric latent heating on typhoon movement crossing Taiwan：the case of Morakot （2009） with extreme rainfall[J].J Atmos Sci，69：31723196.

Willoughby H E，1990.Temporal changes of the primary circulation in tropical cyclones[J].J Atmos Sci，47：242264.

Willoughby H E，1992.Linear motion of a shallow water barotropic vortex as an initialvalue problem[J].J Atmos Sci，47：242264.

Wong M L M，Chan J C L，2006.Modeling the effects of landsea roughness contrast on tropical cyclone winds[J].J Atmos Sci，64：32493264.

Wu L，Wang B，2000.A potential vorticity tendency diagnostic approach for tropical cyclone motion[J].Mon Wea Rev，128：18991911.

Wu L，Wang B，2001a.Movement and vertical coupling of adiabatic baroclinic tropical cyclones[J].J Atmos Sci，58：18011814.

Wu L，Wang B，2001b.Effects of convective heating on movement and vertical coupling of tropical cyclones：a numerical study[J].J Atmos Sci，58：36393649.

Wu L，Braun S A，Halverson J，et al.，2006.A numerical study of Hurricane Erin （2001）.Part I：model verification and storm evolution[J].J Atmos Sci，63：6586.

吳玉琴，吳立廣，梁佳，2015.利用位渦趨勢法診斷臺風“莫拉克”（2009）的運動[J].氣象學報，73（2）：236246.Wu Y Q，Wu L G，Liang J，2015.Diagnosi of the movement of Typhoon Morakot （2009） using the potential vorticity tendency approach[J].Acta Meteor Sinica，73（2）：236246.（in Chinese）.

Xie B，Zhang F，2012.Impacts of typhoon track and island topography on the heavy rainfalls in Taiwan associated with Morakot （2009）[J].Mon Wea Rev，140：33793394.

Yang M J，Zhang D L，Huang H L，2008.A modeling study of typhoon Nari （2001） at landfall.Part I：topography effects[J].J Atmos Sci，65：30953115.

楊瓊瓊，吳立廣，2015.中國24 h臺風路徑預報難點及其大尺度環(huán)流分析[J].氣象科學，35（4）：389397.Yang Q Q，Wu L G，2015.Analysis on difficulty of 24 h typhoon track forecast in China and associated largescale circulation[J].J Meteor Sci，35（4）：389397.（in Chinese）.

Yeh T C，Elsberry R L，1993a.Interaction of typhoons with the orography.Part I：upstream track deflections[J].Mon Wea Rev，121：31933212.

Yeh T C，Elsberry R L，1993b.Interaction of typhoons with the orography.Part II：continuous and discontinuous tracks across the island[J].Mon Wea Rev，121：32133233.

張嬌艷，吳立廣，張強，2011.全球變暖背景下我國熱帶氣旋災害趨勢分析[J].熱帶氣象學報，27（4）：442454.Zhang J Y，Wu L G，Zhang Q，2011.Tropical cyclone damages in China under the background of global warming[J].J Tropical Meteor，27（4）：442454.（in Chinese）.

Zhang Q，Liu Q，Wu L，2009.Tropical cyclone damages in China 1983—2006[J].Bull Amer Meteor Soc，90：489495.

Idealized numerical study on influence of topography along East China coast on landfall tropical cyclone motion