周 才 揚，殷 成 團(tuán)，章 衛(wèi) 勝，張 金 善

(1.長江水利委員會水文局 長江口水文水資源勘測局，上海 200136； 2.南京水利科學(xué)研究院，江蘇 南京 210029)

近年來，強臺風(fēng)和超強臺風(fēng)登陸我國的頻數(shù)明顯增多[1]，其誘發(fā)的風(fēng)暴潮在河口海岸地區(qū)造成了巨大的經(jīng)濟(jì)損失。2017年臺風(fēng)“天鴿”在珠海登陸，最大風(fēng)速達(dá)45 m/s，中心氣壓950 hPa，造成最大增水2.79 m，直接經(jīng)濟(jì)損失51.54億元[2]；2018年“山竹”在廣東省臺山市登陸，最大風(fēng)速達(dá)48 m/s，中心氣壓955 hPa，造成最大增水3.39 m，直接經(jīng)濟(jì)損失24.57億元[3]；2019年臺風(fēng)“利奇馬”橫掃我國南北沿海地區(qū)，帶來大范圍降雨，直接經(jīng)濟(jì)損失超過500億元。長江口作為我國第一大河口屢受臺風(fēng)暴潮襲擾，造成嚴(yán)重經(jīng)濟(jì)損失。近70 a共計160場臺風(fēng)影響長江口地區(qū)，年均2.3場，而2018年臺風(fēng)數(shù)多達(dá)6.0場。

以往對熱帶氣旋活動變化趨勢的研究多集中在臺風(fēng)頻數(shù)及登陸頻數(shù)[4-6]，臺風(fēng)對區(qū)域的影響程度與臺風(fēng)強度、距離、最大風(fēng)速半徑、持續(xù)時間及登陸位置等密切相關(guān)，頻數(shù)只能反映熱帶氣旋對近岸影響能力的一部分，綜合評價仍需其他指標(biāo)予以補充。Emanuel采用氣旋能量指數(shù)PDI(Power Dissipation Index)研究了西太平洋熱帶氣旋與海表溫度的相關(guān)關(guān)系[7]，Bell等采用熱帶氣旋累計能量ACE(Accumulated Cyclone Energy)統(tǒng)計分析了1950～1999年熱帶風(fēng)暴和颶風(fēng)逐季變化特征[8]，Yu等在2009年將熱帶氣旋累積能力進(jìn)行了改進(jìn)(the Revised Accumulated Cyclone Energy Index)[9]，Lin等采用熱帶氣旋潛在影響指數(shù)PI(the Potential Intensity Index)建立了熱帶氣旋和海表溫度、拖曳力等變量之間的聯(lián)系[10]，Haig等采用氣旋活動指數(shù)(Cyclone Activity Indices)對熱帶氣旋活動特征變化進(jìn)行了評估[11]。上述指標(biāo)從臺風(fēng)內(nèi)部主要構(gòu)造特征參數(shù)對臺風(fēng)影響能力進(jìn)行了刻畫，其描述對象為大范圍區(qū)域如大洋乃至全球。 而臺風(fēng)在小范圍地區(qū)影響能力評估指標(biāo)尚不多見。影響長江口的熱帶氣旋活動特征尚不明晰，本文旨在從風(fēng)暴潮防災(zāi)減災(zāi)的角度對影響該地區(qū)的臺風(fēng)活動特征進(jìn)行長周期統(tǒng)計研究。

1 長江口臺風(fēng)影響圓假定

目前，評價熱帶氣旋對小區(qū)域的影響能力尚無較好的指標(biāo)。本文根據(jù)以往對風(fēng)暴潮的模擬研究和7級風(fēng)圈半徑大小，提出了長江口臺風(fēng)影響圓假定。即假定以長江口概化點(122°E，31.5°N)為中心，8倍最大風(fēng)速半徑Rmax為半徑(約400 km)的圓形區(qū)域為長江口臺風(fēng)影響圓。當(dāng)臺風(fēng)路徑與長江口臺風(fēng)影響圓相交或相切時，則認(rèn)為該臺風(fēng)對長江口造成影響。影響長江口的臺風(fēng)路徑概況及長江口臺風(fēng)影響圓如圖1所示，圖中灰色線為未對長江口造成影響的臺風(fēng)路徑。

圖1 1949～2018年影響長江口的臺風(fēng)路徑Fig.1 Tracks of typhoons affecting the Yangtze River Estuary from 1949 to 2018

2 長江口臺風(fēng)影響能力分析

2.1 長江口臺風(fēng)影響頻數(shù)分析

基于最近70 a來(1949～2018年)的熱帶氣旋最佳路徑資料，本節(jié)從統(tǒng)計學(xué)角度對影響我國各類熱帶氣旋尤其是近10 a的變化規(guī)律進(jìn)行分析。1949～2018年影響長江口地區(qū)臺風(fēng)共計160.0場，平均每年登陸2.3場，占西太平洋海區(qū)發(fā)生臺風(fēng)數(shù)的8.92%。影響長江口臺風(fēng)按路徑的不同可分為兩大類，第一類為“登陸北上型”臺風(fēng)，共計92場，占總數(shù)的57.86%；第二類為“外?；顒颖鄙闲汀迸_風(fēng)，占總數(shù)的42.14%?！暗顷懕鄙闲汀迸_風(fēng)中登陸點離長江口較近的“長江口直接登陸型臺風(fēng)”危害最大，按其風(fēng)向變化過程又可分為長江口南側(cè)登陸(87%)和北側(cè)登陸(13%)兩類。還有少數(shù)奇怪路徑的臺風(fēng)，因數(shù)量較少本文不做分析。

圖2為1949～2018年長江口及全國臺風(fēng)生成時間分布圖。由圖2可以看出：不同于全國統(tǒng)計情況，影響長江口的臺風(fēng)只發(fā)生在汛期，即每年的5～11月份；其中7～8月份最多，占臺風(fēng)總數(shù)的64.8%，同時7～8月份也是我國中東部地區(qū)雨量比較集中的月份，暴雨和上游洪水加劇了長江口河口風(fēng)暴增水過程。從生成日方面看(餅圖)，臺風(fēng)產(chǎn)生的時間為中旬(39%)>下旬(32%)>上旬(28%)，月內(nèi)差別較小。

圖2 1949～2018年影響中國臺風(fēng)生成時間分布Fig.2 Formation time of typhoons affecting China from 1949 to 2018

2.2 長江口氣旋活動強度指數(shù)分析

為量化臺風(fēng)給長江口造成影響，采用熱帶氣旋活躍指數(shù)(Cyclone Activity Index)、區(qū)域累計氣旋能量(Regional Accumulated Cyclone Energy)、年影響臺風(fēng)最大風(fēng)速Vmax和年影響臺風(fēng)頻數(shù)C四個指標(biāo)對長江口的氣旋影響強度進(jìn)行綜合評估。

Cyclone Activity Index(CAI)是Haig[11]等基于NOAA和其他氣象專家對氣旋強度和尺度研究進(jìn)行改進(jìn)，提出的一個新的氣旋活動指數(shù)。CAI從風(fēng)暴數(shù)量、強度、大小和壽命等方面描述了年均臺風(fēng)破壞性大小。其公式為

(1)

當(dāng)臺風(fēng)距長江口尚遠(yuǎn)如仍在西太平洋面活動時，除從大洋傳到近岸的先兆波外，臺風(fēng)氣旋對近岸水體影響很小。為描述長江口地區(qū)臺風(fēng)破壞能力影響能力，本節(jié)基于長江口臺風(fēng)影響圓假定定義一個新的指標(biāo)——區(qū)域累積氣旋能量(Regional Accumulated Cyclone Energy，RACE)。該指標(biāo)認(rèn)為只有臺風(fēng)眼進(jìn)入長江口影響圓范圍內(nèi)的臺風(fēng)，才會在河口地區(qū)造成較為嚴(yán)重的災(zāi)害，此類臺風(fēng)路徑在影響圓內(nèi)的部分便被計入RACE指標(biāo)。同時，RACE在風(fēng)速權(quán)重方面的考慮取了平方，更貼合能量的概念，其公式為

(2)

式中：i為每年進(jìn)入長江口臺風(fēng)影響圓內(nèi)臺風(fēng)路徑點數(shù)；V為最大風(fēng)速半徑處風(fēng)速，m/s；di為第i個路徑點到長江口的距離，km。

臺風(fēng)是一類中尺度天氣事件，其范圍遠(yuǎn)遠(yuǎn)大于河口等研究區(qū)域。CAI和RACE都是評估小研究區(qū)域的指標(biāo)，但是其側(cè)重點不同。CAI空間計算范圍是西太平洋所有發(fā)生的熱帶氣旋，時間計算范圍是熱帶氣旋從生成到消亡的整個生命周期。其側(cè)重點在于氣旋對區(qū)域的影響能力。但有些臺風(fēng)距離目標(biāo)地區(qū)較遠(yuǎn)，風(fēng)速不大，卻依然被計如參數(shù)計算范圍。而指標(biāo)RACE基于“長江口臺風(fēng)影響圓”假設(shè)提出，更強調(diào)研究地區(qū)的臺風(fēng)能量及災(zāi)害。2個指標(biāo)各有側(cè)重，計算結(jié)果可作相互補充。

年影響臺風(fēng)最大風(fēng)速Vmax為長江口影響圓內(nèi)出現(xiàn)過的年海表最大風(fēng)速，單位m/s，這個指標(biāo)主要是對臺風(fēng)強度進(jìn)行刻畫，尤其是強臺風(fēng)和超強臺風(fēng)。年影響臺風(fēng)頻數(shù)C為長江口臺風(fēng)影響圓內(nèi)出現(xiàn)過的臺風(fēng)個數(shù)，主要從年臺風(fēng)頻數(shù)方面進(jìn)行描述。

從氣旋活動指數(shù)(CAI)看，近70 a熱帶氣旋按其活躍度可分為兩部分：第一部分1949～1985年，CAI指數(shù)在平均線上下波動劇烈，表明熱帶氣旋影響能力極不穩(wěn)定，處于活躍期；第二部分1986～2018年，除2015年外，其余CAI指數(shù)在平均線以下小幅波動，表明該階段內(nèi)熱帶氣旋對長江口地區(qū)的影響能力較為穩(wěn)定。

區(qū)域累積氣旋能量(RACE)結(jié)果表明，長江口的臺風(fēng)影響年際變化較大。1976～1991年的5 a移動平均值(虛線)均超過多年平均，表明這15 a中氣旋影響能力較強；1991～2018年氣旋影響能力較為穩(wěn)定，僅部分年份大于均值。

年影響臺風(fēng)最大風(fēng)速Vmax歷史最值達(dá)到84 m/s(5 907超強臺風(fēng)Sarah )。2010年之后最大風(fēng)速都一般在50 m/s左右，稍高于平均，但未出現(xiàn)異常值。從本文指標(biāo)來看，近年登陸我國的強臺風(fēng)和超強臺風(fēng)雖然有所增多[1]，但是因為距離長江口較遠(yuǎn)，影響有限。年影響臺風(fēng)頻數(shù)C沒有比較明顯的長期規(guī)律，從短期角度看，自2009年之后一直高于歷年平均值，且保持上升趨勢，并在2018年達(dá)到歷史極值6場。

CAI，RACE和Vmax三個指標(biāo)的5 a滑動平均線均表明：近30 a來臺風(fēng)對長江口的影響程度并不算高(略低于均值)，但自2009年之后長江口臺風(fēng)影響能力一直保持增長趨勢，且近年臺風(fēng)影響頻數(shù)較大(見圖3)。

圖3 1949～2018年長江口臺風(fēng)影響能力年際變化Fig.3 Interannual variation of typhoon impact capacity in the Yangtze Estuary from 1949 to 2018

3 結(jié) 論

利用1949～2018年最佳臺風(fēng)路徑資料，在長江口臺風(fēng)影響圓假定的前提下，提出了臺風(fēng)對小區(qū)域影響能力的評估指標(biāo)，并采用多個指標(biāo)對長江口臺風(fēng)影響能力進(jìn)行了長周期的統(tǒng)計和研究，結(jié)果如下：

(1) 本文提出的指標(biāo)區(qū)域累積氣旋能量(RACE)，可以較好地評估臺風(fēng)對小區(qū)域的影響能力。

(2) 1949～1985年熱帶氣旋影響能力極不穩(wěn)定，處于活躍期；1986～2018年表明熱帶氣旋影響能力較低且強度穩(wěn)定，影響長江口的強臺風(fēng)和超強臺風(fēng)有所減少。

(3) 1991～2018年的近30 a來臺風(fēng)對長江口的影響程度并不算高(略低于均值)，但自2009年之后臺風(fēng)影響指數(shù)一直保持增長趨勢，且近年臺風(fēng)影響頻數(shù)較大。