駱方露，朱敏，李江南

（1.中山大學(xué)大氣科學(xué)學(xué)院/廣東省氣候變化與自然災(zāi)害研究重點(diǎn)實(shí)驗(yàn)室，廣東珠海519082；2.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室（珠海），廣東珠海519082）

1 引言

隨著全球變暖形勢(shì)愈發(fā)嚴(yán)峻，全球極端天氣氣候事件及其造成的災(zāi)害頻率呈現(xiàn)出增加的趨勢(shì)。熱帶氣旋（Tropical Cyclone，TC）作為一種破壞性極大的全球性氣象災(zāi)害，其活動(dòng)與全球平均氣溫密切相關(guān)。多年來(lái)，許多國(guó)內(nèi)外研究學(xué)者致力于TC的研究，剖析其發(fā)生機(jī)理和時(shí)空變化規(guī)律，這對(duì)于提高TC預(yù)報(bào)水平以及減少TC的發(fā)生風(fēng)險(xiǎn)和災(zāi)害損失具有重要意義。

過(guò)去對(duì)TC活動(dòng)趨勢(shì)的研究主要集中在頻數(shù)和強(qiáng)度上。Chan等[1]發(fā)現(xiàn)20世紀(jì)70年代中期后，西北太平洋TC發(fā)生數(shù)量呈現(xiàn)增長(zhǎng)趨勢(shì)，而北大西洋的變化趨勢(shì)恰恰相反，這一結(jié)論與Landsea等[2]一致，上述研究表明全球總TC活動(dòng)可能保持相對(duì)穩(wěn)定，一個(gè)海區(qū)TC活動(dòng)的增加可能是另一個(gè)海區(qū)活動(dòng)減少所補(bǔ)償產(chǎn)生的。Camargo等[3]發(fā)現(xiàn)在厄爾尼諾年，TC強(qiáng)度比拉尼娜年更大，持續(xù)時(shí)間也更長(zhǎng)。楊亞新等[4]指出2006—2015年間TC發(fā)生頻數(shù)呈減小趨勢(shì)，7—10月是TC發(fā)生的活躍期。對(duì)于西北太平洋海域2018年8月出現(xiàn)的TC頻數(shù)異常，張潤(rùn)宇等[5]研究發(fā)現(xiàn)其與850 hPa相對(duì)渦度和600 hPa相對(duì)濕度之間存在良好的正相關(guān)關(guān)系，這兩個(gè)因子正異常的出現(xiàn)，為TC的生成提供了有利條件。較大的垂直切變會(huì)直接導(dǎo)致北大西洋熱帶風(fēng)暴和颶風(fēng)頻數(shù)的減少[6]。

然而，Landsea等[2]發(fā)現(xiàn)1944—1995年北大西洋出現(xiàn)了強(qiáng)颶風(fēng)頻數(shù)減少和TC強(qiáng)度減弱的趨勢(shì)，但在這期間TC造成的損失卻是巨大的，這說(shuō)明了TC頻數(shù)和強(qiáng)度若用于評(píng)估TC造成的損失，還存在很大的局限性。1969年，土木工程師Saffir提出薩菲爾-辛普森颶風(fēng)等級(jí)（Saffir-Simpson Hurricane Scale，SSHS）這個(gè)概念，它利用最大持續(xù)風(fēng)速評(píng)定颶風(fēng)強(qiáng)度，后來(lái)被廣泛應(yīng)用于颶風(fēng)預(yù)警。Kantha[7]提出了連續(xù)的颶風(fēng)強(qiáng)度指數(shù)（Hurricane Intensity Index，HII）和考慮颶風(fēng)尺度的颶風(fēng)危害指數(shù)（Hurricane Hazard Index，HHI），通過(guò)對(duì)比這兩個(gè)指數(shù)和SSHS評(píng)定的颶風(fēng)強(qiáng)度，發(fā)現(xiàn)SSHS的評(píng)定結(jié)果與颶風(fēng)實(shí)際造成的損失相差很大，因?yàn)樗鼪]有考慮颶風(fēng)尺度，忽視了強(qiáng)度弱但尺度大的TC可能比強(qiáng)度大但尺度小的TC造成更大破壞的情況。因此，為了更好地描述TC的破壞力，很多研究學(xué)者一直嘗試開發(fā)各種新的影響指數(shù)，比如Bell等[8]提出了累積氣旋能量指數(shù)（Accumulated Cyclone Energy，ACE），定義為整個(gè)TC壽命期間最大風(fēng)速的平方和；Emanuel[9]基于能量耗散概念提出了類似ACE的能量耗散指數(shù)（Power Dissipation Index，PDI），定義為TC壽命內(nèi)最大持續(xù)風(fēng)速立方的時(shí)間積分；Powell等[10]基于TC的綜合動(dòng)能（Integrated Kinetic Energy，IKE）提出了風(fēng)破壞力評(píng)級(jí)指數(shù)（Wind Destructive Potential，WDP）和風(fēng)暴潮與波浪破壞力評(píng)級(jí)指數(shù)（Storm surge and wave Destructive Potential，SDP），這兩個(gè)評(píng)級(jí)指標(biāo)中納入了TC的半徑信息，充分考慮了尺度在TC破壞力中的貢獻(xiàn)。上述指數(shù)的提出綜合考慮了TC的頻數(shù)、持續(xù)時(shí)間和強(qiáng)度，彌補(bǔ)了單獨(dú)研究頻數(shù)和強(qiáng)度的局限性[11]。王功錄[12]對(duì)1949—2011年登陸我國(guó)TC的ACE進(jìn)行分析，發(fā)現(xiàn)基本呈減弱趨勢(shì)。1982—2014年西北太平洋TC的ACE年際變化在整體上也呈現(xiàn)緩慢下降的趨勢(shì)[13]。屠佳雨等[14]研究發(fā)現(xiàn)，1949—2014年影響海南島的西北太平洋臺(tái)風(fēng)年際總PDI變化趨勢(shì)表現(xiàn)為緩慢減小。季倩倩[15]的研究結(jié)果表明1979—2018年西北太平洋登陸TC的PDI表現(xiàn)為顯著的增長(zhǎng)趨勢(shì)，且在研究時(shí)間范圍內(nèi)增大了33%。喬守文等[16]發(fā)現(xiàn)TC頻數(shù)和平均最大風(fēng)速的增加促進(jìn)了累積PDI的增長(zhǎng)。

海洋是TC的主要能量來(lái)源，更高的海表溫度（Sea Surface Temperature，SST）降低了大氣穩(wěn)定性，增加了渦的穿透深度，使得TC能夠抵抗垂直風(fēng)切變繼續(xù)發(fā)展[17]。一般來(lái)說(shuō)，SST達(dá)26.5℃以上是TC發(fā)展的必要條件。Emanuel[18]和Knutson等[19]預(yù)測(cè)TC強(qiáng)度將隨著全球平均溫度的升高而增加。赤道東太平洋SST的年際變化與西北太平洋中、西部臺(tái)風(fēng)頻數(shù)之間有很強(qiáng)的滯后相關(guān)關(guān)系[20]。Goldenberg等[17]認(rèn)為1995—2000年以來(lái)北大西洋颶風(fēng)的異?；钴S可能是北大西洋SST和垂直風(fēng)切變的多年代際變化，以及全球變暖導(dǎo)致SST進(jìn)一步增加的共同作用的結(jié)果。Emanuel[9]發(fā)現(xiàn)自20世紀(jì)70年代中期以來(lái)TC能量耗散明顯增加，這是由于TC壽命延長(zhǎng)和強(qiáng)度增大所致。TC能量耗散與SST變化之間存在很強(qiáng)的相關(guān)關(guān)系，但觀測(cè)到的能量耗散增加僅部分是由SST增加直接引起。

綜上所述，前人在TC氣候特征方面做了大量的研究工作，其中關(guān)于頻率和強(qiáng)度的統(tǒng)計(jì)研究很多，關(guān)于破壞力的相對(duì)較少，并且以往的研究多集中于西太平洋和北大西洋，很少涉及到全球范圍內(nèi)的其他海區(qū)，影響指數(shù)與SST之間的相關(guān)關(guān)系也較少被提及。因此，本研究旨在應(yīng)用PDI、WDP和SDP對(duì)全球范圍內(nèi)的TC進(jìn)行統(tǒng)計(jì)，分析全球TC破壞力的時(shí)間變化和空間分布特征，并探討各影響指數(shù)與SST變化之間的聯(lián)系。

2 數(shù)據(jù)和指標(biāo)的定義

2.1 數(shù)據(jù)和研究方法

本文使用的TC資料來(lái)自國(guó)際熱帶氣旋最佳路徑資料集（International Best Track Archive for Climate Stewardship，IBTrACS），它集合了來(lái)自全球多個(gè)區(qū)域?qū)I(yè)氣象中心和其他機(jī)構(gòu)的TC數(shù)據(jù)集信息，包含間隔3 h的TC位置、強(qiáng)度和其他信息。全球海區(qū)分為北大西洋（North Atlantic，NA）、東太平洋（East Pacific ocean，EP）、西太平洋（West Pacific ocean，WP）、北印度洋（North India，NI）和南半球海區(qū)（Southern Hemisphere，SH，包含南印度洋、南太平洋和南大西洋），各海區(qū)的SST平均范圍選自2000—2010年歷史TC（≥18 m/s）路徑覆蓋的范圍：北大西洋（0°～60°N，0°～100°W）、東太平洋（0°～35°N，100°～180°W）、西太平洋（0°～60°N，100°～180°E）、北印度洋（0°～25°N，50°～100°E）、南印度洋（0°～45°S，30°～90°E）、南太平洋（0°～60°S，160°E～120°W）和南大西洋（20°～40°S，30°～60°W）[21]。

SST資料來(lái)自Hadley中心海冰和海表溫度（Hadley center sea Ice and Sea Surface Temperature，HadISST）資料集[22]，該數(shù)據(jù)集由全球月平均海溫?cái)?shù)據(jù)和海冰數(shù)據(jù)組成，分辨率為1°×1°，覆蓋范圍為89.5°S～89.5°N，179.5°W～179.5°E，時(shí)間跨度為1870年至今。

為分析指數(shù)的年際變化序列，本文采用一元線性回歸的統(tǒng)計(jì)學(xué)方法，對(duì)回歸系數(shù)進(jìn)行t檢驗(yàn)，并計(jì)算皮爾遜相關(guān)系數(shù)。文中對(duì)部分時(shí)間序列作了功率譜分析。此外，在對(duì)去除了線性趨勢(shì)的全球海表溫度異常（Sea Surface Temperature Anomaly，SSTA）進(jìn)行統(tǒng)計(jì)分析時(shí)，利用了經(jīng)驗(yàn)正交函數(shù)（Empirical Orthogonal Functio，EOF）分解，并計(jì)算特征值誤差范圍對(duì)主分量進(jìn)行顯著性檢驗(yàn)[23]。

2.2 指標(biāo)定義

2.2.1 能量耗散指數(shù)

本文將應(yīng)用Emanuel[9]定義的TC能量總耗散指數(shù)PDI，計(jì)算公式如下：

式中，VMS是最大持續(xù)風(fēng)速（單位：m/s）；τ是TC壽命（單位：s）；PDI單位：m3/s2。

對(duì)于PDI各大海區(qū)的統(tǒng)計(jì)時(shí)間分別是：北大西洋、東太平洋和西太平洋為1949—2019年，北印度洋和南半球?yàn)?969—2019年（見表1）。

2.2.2 破壞力指數(shù)

Powell等[10]提出用IKE作為破壞力的指標(biāo)，證明了IKE與風(fēng)破壞力以及風(fēng)暴潮與波浪破壞力之間的相關(guān)性，并由此提出了評(píng)定兩種破壞力的連續(xù)評(píng)級(jí)方法——WDP和SDP，分別從風(fēng)荷載和海表應(yīng)力的物理基礎(chǔ)上傳遞破壞力。

本文將應(yīng)用WDP和SDP指數(shù)對(duì)全球TC進(jìn)行分級(jí)。方法是通過(guò)Powell等[10]提供的公式，使用IBTrACS中風(fēng)半徑和最大風(fēng)速數(shù)據(jù)進(jìn)行計(jì)算。

WDP的計(jì)算公式如下：

當(dāng)TC的最大風(fēng)速＜55 m/s時(shí)：

當(dāng)TC的最大風(fēng)速≥55 m/s時(shí)：

SDP的計(jì)算公式如下：

式 中，R18、R33和Rmax分別 是TC風(fēng)速 達(dá)到18 m/s、33 m/s和最大風(fēng)速時(shí)的風(fēng)半徑（單位：km）；VMS是最大持續(xù)風(fēng)速（單位：m/s）。

由于風(fēng)半徑數(shù)據(jù)自2004年后才有完整記錄，對(duì)各大海區(qū)的WDP和SDP的分級(jí)統(tǒng)計(jì)時(shí)間為2004—2019年（見表1）。

表1 全球各海區(qū)指數(shù)的統(tǒng)計(jì)時(shí)間

3 全球TC破壞力的時(shí)空變化特征

3.1 破壞力的月際變化特征

統(tǒng)計(jì)結(jié)果表明（見圖1），北大西洋、東太平洋和西太平洋的PDI均呈現(xiàn)顯著的單峰分布，主要集中在8—10月，峰值均出現(xiàn)在9月。北印度洋TC生成數(shù)量少，且大多較弱、壽命短，PDI數(shù)值相較其他海區(qū)要小得多，分布呈雙峰型，主要集中在5—6月和10—11月，峰值分別在5月和11月。北半球月均累積PDI最大的海區(qū)為西太平洋，其次為東太平洋和北大西洋，北印度洋最小。南半球海區(qū)的PDI主要集中在1—3月，峰值出現(xiàn)在3月。

圖1 全球海域累積PDI的月際分布

圖2 展示了全球海域WDP和SDP累積評(píng)級(jí)的月際分布。由圖可見，北大西洋、東太平洋、西太平洋和南半球的WDP和SDP均表現(xiàn)為顯著的單峰型結(jié)構(gòu)。北大西洋WDP和SDP的活躍期為8—10月，峰值均出現(xiàn)在9月。東太平洋WDP和SDP主要集中在7—10月，峰值出現(xiàn)在8月。西太平洋WDP和SDP活躍在7—10月，峰值在9月。南半球的WDP和SDP主要集中在夏秋季，即1—3月，峰值出現(xiàn)在3月。與其他海區(qū)不同，北印度洋WDP和SDP呈雙峰型分布，主要集中在4—5月和10—11月，峰值分別為5月和10月。

圖2 2004—2019年全球海域WDP和SDP累積評(píng)級(jí)的月際分布

上述結(jié)果表明，北大西洋、東太平洋、西太平洋和南半球的PDI、WDP和SDP均活躍在夏秋季節(jié)，北半球海區(qū)主要集中在7—10月，峰值出現(xiàn)在8—9月，南半球季節(jié)與北半球相反，破壞力主要集中在1—3月，峰值出現(xiàn)在3月。與其他海區(qū)不同，北印度洋表現(xiàn)出特殊的雙峰型結(jié)構(gòu)，主要集中在4—6月和10—11月，峰值分別出現(xiàn)在5月和10—11月。梁梅等[24]在對(duì)北印度洋TC的頻數(shù)特征進(jìn)行統(tǒng)計(jì)時(shí)也發(fā)現(xiàn)類似結(jié)果，即在該海區(qū)，春末夏初和秋季是最有利于TC發(fā)生發(fā)展的時(shí)期；此外，作者還發(fā)現(xiàn)雙峰型結(jié)構(gòu)的成因是兩個(gè)高峰期時(shí)，北印度洋SST較高，且正值冬、夏季風(fēng)轉(zhuǎn)換期間，因此為TC的生成和發(fā)展提供了有利條件。

由WDP和SDP平均評(píng)級(jí)的月際變化可知（見圖3），北大西洋的WDP和SDP表現(xiàn)為單峰型分布，主要集中在7—11月，峰值分別出現(xiàn)在9月和10月，SDP的全年最大值出現(xiàn)在1月。據(jù)美國(guó)國(guó)家颶風(fēng)中心（National Hurricane Center，NHC）整理的歷史TC信息可知，2016年大西洋颶風(fēng)季相當(dāng)活躍，TC的活動(dòng)強(qiáng)度高于1981—2010年的平均水平，1月出現(xiàn)的異常值正是由發(fā)生于該年的颶風(fēng)“Alex”造成的，它貢獻(xiàn)了全年最高的SDP值。此外，SDP在5月也存在一個(gè)異常高值。與11月相比，在研究時(shí)間范圍內(nèi)，5月產(chǎn)生了5個(gè)熱帶風(fēng)暴級(jí)的TC，11月產(chǎn)生了5個(gè)熱帶風(fēng)暴和4個(gè)颶風(fēng)，但5月和11月的平均SDP值相差不大，分別為2.85和2.73。這說(shuō)明在某些TC生成數(shù)量較少且強(qiáng)度較弱的月份，也可能具有較強(qiáng)的風(fēng)暴潮破壞力，若僅對(duì)頻數(shù)或強(qiáng)度進(jìn)行單一的統(tǒng)計(jì)研究，則難以準(zhǔn)確預(yù)估TC造成的影響。在東太平洋，WDP和SDP呈現(xiàn)雙峰型分布，前者峰值分別出現(xiàn)在6月和9月，后者分別在4月和9月。考慮SDP在4月出現(xiàn)的異常高值，它是由2004年的超強(qiáng)臺(tái)風(fēng)“Sudal”產(chǎn)生的，與颶風(fēng)“Alex”相似，臺(tái)風(fēng)“Sudal”作為研究時(shí)間范圍內(nèi)唯一在4月產(chǎn)生的TC，造成了SDP的最大值；不同的是，它的強(qiáng)度和破壞程度都比颶風(fēng)“Alex”大得多。臺(tái)風(fēng)“Sudal”在其生命史中共造成了約1 400萬(wàn)美元的損失，其中遭受損失最大的區(qū)域是密克羅尼西亞聯(lián)邦雅蒲州。臺(tái)風(fēng)“Sudal”經(jīng)過(guò)雅蒲州時(shí)，給島上帶來(lái)狂風(fēng)暴雨，摧毀大量建筑物，島上的供水、供電和通訊系統(tǒng)遭到重創(chuàng)。沿岸受巨浪和風(fēng)暴潮影響，船只沉沒，海堤被破壞，風(fēng)暴潮還造成海水入侵，地下水被污染，災(zāi)難性地摧毀島上的農(nóng)作物。對(duì)比颶風(fēng)“Alex”和臺(tái)風(fēng)“Sudal”的WDP和SDP，可以發(fā)現(xiàn)兩者的SDP相差不大，但臺(tái)風(fēng)“Sudal”的WDP比颶風(fēng)“Alex”高一倍多，這說(shuō)明強(qiáng)烈的風(fēng)破壞力是導(dǎo)致臺(tái)風(fēng)“Sudal”造成的影響大于颶風(fēng)“Alex”的關(guān)鍵原因。西太平洋的WDP沒有顯著峰值，除1月以外，其他月份的WDP值都比較高；SDP表現(xiàn)為單峰型分布，峰值出現(xiàn)在9月，且全年均呈現(xiàn)較活躍的狀態(tài)。北印度洋的WDP和SDP呈現(xiàn)雙峰型分布，前者峰值分別出現(xiàn)在6月和10月，后者在5月和9月。在南半球，WDP和SDP除6—9月外均處于較活躍的狀態(tài)，無(wú)顯著峰值。

圖3 2004—2019年全球海域WDP和SDP平均評(píng)級(jí)的月際分布

由上述分析可知，全球各海區(qū)WDP和SDP平均評(píng)級(jí)與累積評(píng)級(jí)的月際分布特征具有顯著差異。相較而言，平均WDP和SDP具有更長(zhǎng)的活躍期，這是受到了TC頻數(shù)月際變化的影響。平均WDP和SDP會(huì)凸顯出頻數(shù)小但破壞力大的月份，這解釋了一些累積評(píng)級(jí)較低的月份在平均評(píng)級(jí)中反而較高的現(xiàn)象，也強(qiáng)調(diào)了僅使用頻數(shù)來(lái)評(píng)估TC的局限性。此外，還有一些風(fēng)破壞力較小的月份，風(fēng)暴潮與波浪破壞力卻不低，如北大西洋的1月、5—6月和12月，西太平洋的1月，北印度洋的8—9月和12月，這些月份生成的TC在以風(fēng)力為標(biāo)準(zhǔn)的傳統(tǒng)分級(jí)方式中評(píng)級(jí)較低，但實(shí)際造成的破壞會(huì)高于其所獲得的評(píng)級(jí)，從而產(chǎn)生超出預(yù)期的損失，為TC預(yù)警和防災(zāi)政策的制定帶來(lái)相當(dāng)大的阻力。

3.2 破壞力的年際變化特征

由圖4可知，北大西洋PDI整體呈上升趨勢(shì)，上升速率為7.28×1012m3/（s2·a），在96%的置信水平上顯著；其年際分布表現(xiàn)為一個(gè)相對(duì)平穩(wěn)階段和一個(gè)起伏階段，分別是1949—1965年和1966—2019年，其中，1972—1994年為最不活躍的時(shí)期；對(duì)時(shí)間序列作功率譜分析（見圖5a），發(fā)現(xiàn)北大西洋的PDI具有顯著的3 a左右的周期。從長(zhǎng)期趨勢(shì)上看，東太平洋的PDI顯示出上升趨勢(shì)，上升速率為1.92×1010m3/（s2·a），在99%的置信水平上顯著；從年際分布上看，它具有明顯的波浪式起伏變化特征，這主要表現(xiàn)在兩個(gè)起伏階段，分別是1949—1977年和1978—2010年，此外，自2011年至今，它又表現(xiàn)出上升趨勢(shì)；在該海區(qū)，PDI具有顯著的2 a左右的周期（見圖5b）。西太平洋的PDI未顯示出顯著的線性變化趨勢(shì)，其年際分布大致呈現(xiàn)為兩個(gè)起伏階段，分別是1949—1977年和1978年至今，它具有顯著的3 a左右周期和35 a左右的年代際周期（圖略）。北印度洋PDI呈現(xiàn)上升趨勢(shì)，上升速率為2.97×109m3/（s2·a），在92%的置信水平上顯著；該海區(qū)PDI的年際變化趨勢(shì)大致表現(xiàn)為，1993—2009年為減弱期，1969—1992年和2009年至今為增長(zhǎng)期，由圖5d可知，該海區(qū)PDI具有顯著的7 a左右周期。在南半球海區(qū)，PDI整體呈現(xiàn)為上升趨勢(shì)，上升速率為4.63×1010m3/（s2·a），在99%的置信水平上顯著，其年際分布表現(xiàn)為單峰型結(jié)構(gòu)，1969—1996年為增長(zhǎng)期，1997年至今為減弱期，圖5e顯示該海區(qū)PDI具有顯著的3 a左右周期。

圖4 全球海域累積PDI的年際分布（虛線為線性趨勢(shì)擬合）

下面對(duì)平均WDP的年際變化進(jìn)行分析（見圖6）。北大西洋WDP呈現(xiàn)出不顯著的上升趨勢(shì)，上升速率約為0.024/a，其中，2006年、2012年和2013年是16 a間最不活躍的年份。東太平洋WDP未顯示出明顯的線性變化趨勢(shì)，但WDP值總體保持在一個(gè)較高水平。西太平洋的WDP顯示出緩慢的上升趨勢(shì)，上升速率為0.02/a，在96%的置信水平上顯著。北印度洋WDP整體上呈顯著下降趨勢(shì)，下降速率約為-0.12/a，在93%的置信水平上顯著。WDP最大值出現(xiàn)在2004年，由當(dāng)年第25號(hào)臺(tái)風(fēng)“梅花”產(chǎn)生。值得注意的是，2004年出現(xiàn)的異常高值和近4 a出現(xiàn)的異常低活躍現(xiàn)象是線性變化趨勢(shì)呈現(xiàn)為快速下降的主要原因，因此，若在更長(zhǎng)的時(shí)間尺度上進(jìn)行研究，可能會(huì)得到不同的結(jié)論。南半球海區(qū)WDP整體上呈現(xiàn)緩慢上升趨勢(shì)，上升速率約為0.019/a，在94%的置信水平上顯著。

圖6 2004—2019年全球海域WDP平均評(píng)級(jí)的年際分布（柱形數(shù)據(jù)為各海區(qū)年際平均WDP，虛線為線性趨勢(shì)擬合）

圖7 a顯示，北大西洋平均SDP的時(shí)間序列沒有明顯的線性變化趨勢(shì)，但呈現(xiàn)出明顯的年際變化，在2019年SDP達(dá)到最大值3.04。東太平洋SDP存在明顯的下降趨勢(shì)，下降速率約為-0.031/a，在93%的置信水平上顯著。近16 a來(lái)，西太平洋的SDP呈現(xiàn)顯著的上升趨勢(shì)，上升速率為0.029/a，在98%的置信水平上顯著。在2011年以前，SDP存在起伏變化，2005—2007年為增長(zhǎng)期，2008—2010年為減弱期；在此之后，SDP基本保持穩(wěn)定的增長(zhǎng)趨勢(shì)。從長(zhǎng)期趨勢(shì)上看，北印度洋SDP的年際變化表現(xiàn)為緩慢的下降趨勢(shì)，下降速率為-0.016/a，但未通過(guò)90%置信度檢驗(yàn)。北印度洋的年際分布大致呈現(xiàn)為先減少后增加的U字型變化，2005—2010年為減弱期，2013年至今為增長(zhǎng)期。南半球海區(qū)的SDP整體表現(xiàn)為上升趨勢(shì)，上升速率為0.039/a，在94%的置信水平上顯著。南半球的SDP在2011年以前大致經(jīng)歷了一段增長(zhǎng)期，此后維持減弱趨勢(shì)到2016年，直至近3 a又開始呈現(xiàn)增長(zhǎng)趨勢(shì)。

圖7 2004—2019年全球海域SDP平均評(píng)級(jí)的年際分布（柱形數(shù)據(jù)為各海區(qū)年際平均SDP，虛線為線性趨勢(shì)擬合）

綜上可知，全球各海區(qū)PDI年際變化顯示出明顯的波動(dòng)性變化特征，除西太平洋外均呈現(xiàn)顯著的線性增長(zhǎng)趨勢(shì)。此外，PDI年際分布還具有顯著的周期，北大西洋和南半球存在3 a左右的振蕩周期，東太平洋顯示出2 a左右的年際變化周期，北印度洋表現(xiàn)出7 a左右的變化周期，西太平洋則具有3 a左右的年際變化周期和35 a左右的年代際變化周期。平均WDP和SDP變化緩慢，基本圍繞2級(jí)浮動(dòng)，且全球海域變化趨勢(shì)不一致：WDP在西太平洋和南半球表現(xiàn)為上升趨勢(shì)，在北印度洋為下降趨勢(shì)，其余海區(qū)無(wú)明顯變化趨勢(shì)；SDP在西太平洋和南半球呈上升趨勢(shì)，在東太平洋為下降趨勢(shì)，其余海區(qū)變化趨勢(shì)不顯著。

4 SST變化對(duì)TC破壞力的影響

4.1 全球SST的時(shí)空分布及變化特征

圖8 對(duì)去趨勢(shì)的SSTA進(jìn)行EOF分解，從而分析近70 a來(lái)全球SST的時(shí)空變化特征。EOF分解結(jié)果顯示，第一模態(tài)的解釋方差為21.92%，且根據(jù)North等[23]的顯著性檢驗(yàn)方法，該模態(tài)是顯著的。EOF1的空間分布顯示（圖略），太平洋海區(qū)赤道中東太平洋表現(xiàn)為負(fù)異常，且存在全球海區(qū)負(fù)異常的高值中心。南北太平洋的中緯度海區(qū)呈現(xiàn)正異常，全球海區(qū)正異常的高值中心位于北太平洋。印度洋海區(qū)的異常情況與赤道中東太平洋的一致，呈現(xiàn)負(fù)異常。北大西洋整體上表現(xiàn)為正異常。

圖8 1949—2019年全球SSTA的EOF分解第一模態(tài)的時(shí)間系數(shù)

分析該模態(tài)的時(shí)間系數(shù)，發(fā)現(xiàn)其在1977—1998年間以負(fù)異常為主，1976年之前及1999年以后以正異常為主，即在前一時(shí)段，全球SSTA的分布與EOF1是反位相的，而在后兩個(gè)時(shí)段則表現(xiàn)為同位相。對(duì)時(shí)間系數(shù)進(jìn)行功率譜分析的結(jié)果顯示（見圖9），它具有顯著的3～5 a變化周期，中心周期為5 a。

圖9 PC1的功率譜分析（紅線表示95%的置信水平）

4.2 破壞力與SST變化的相關(guān)關(guān)系

由PDI與SST相關(guān)系數(shù)的年際變化可知（見圖10），在研究時(shí)間范圍內(nèi)，PDI與北大西洋、東太平洋、西太平洋和南半球海區(qū)SST均呈現(xiàn)顯著的正相關(guān)關(guān)系（通過(guò)95%置信度檢驗(yàn)），這說(shuō)明東太平洋的PDI在1977—1998年隨著SST的升高表現(xiàn)為增長(zhǎng)的趨勢(shì)，而在1976年以前及1999年以后隨SST的下降呈現(xiàn)減小的趨勢(shì)。其余3個(gè)海區(qū)的變化情況類似，1977—1998年P(guān)DI隨SST的降低而減小，其余時(shí)段則隨SST的升高而增大。北印度洋PDI與SST相關(guān)系數(shù)的年際變化大致可分為兩個(gè)階段，1969—1989年為負(fù)相關(guān)，1990年以后為正相關(guān)。這說(shuō)明北印度洋PDI在1976年以前隨著SST的下降呈現(xiàn)增大趨勢(shì)，1977—1989年間PDI隨著SST的升高呈減小趨勢(shì)，1990—1998年P(guān)DI隨SST升高而增大，1999年至今PDI隨SST下降而減小。綜合PDI與SST的顯著正相關(guān)性說(shuō)明，在全球變暖的背景下，TC能量耗散會(huì)隨著SST的升高而增大，從而產(chǎn)生更大的破壞力。這與Emanuel[9]的研究結(jié)果一致。

圖10 全球海域累積PDI與SST相關(guān)系數(shù)的年際分布（虛線為0.05置信水平）

由圖11可知，除北印度洋外，其余各海區(qū)WDP與SST整體上均呈顯著正相關(guān)。這說(shuō)明在東太平洋，WDP隨著SST的降低而減小，而在北大西洋、西太平洋和南半球，WDP隨SST的升高呈現(xiàn)增大趨勢(shì)。在北印度洋，由于每年生成TC的數(shù)量很少，導(dǎo)致圖中缺少大量相關(guān)系數(shù)的數(shù)據(jù)，因此不具備進(jìn)行長(zhǎng)期趨勢(shì)分析的條件。在該海區(qū)，最大相關(guān)系數(shù)出現(xiàn)在2017年，為-0.99，在95%的置信水平上顯著。

圖11 2004—2019年全球海域平均WDP與SST相關(guān)系數(shù)的年際分布（虛線為0.05置信水平）

圖12 展示了SDP與SST相關(guān)系數(shù)的年際分布。由圖可知，北大西洋的SDP與SST在2014年以前主要呈現(xiàn)為顯著正相關(guān)，自2014年起呈現(xiàn)顯著負(fù)相關(guān)。根據(jù)全球SST的時(shí)空分布可知，2014年前北大西洋大部分海區(qū)表現(xiàn)為增溫，SDP隨著SST的升高而增大，2014年后北大西洋主要表現(xiàn)為降溫，SDP隨著SST的降低而增大。在東太平洋、西太平洋和南半球海區(qū)，SDP整體上與SST均表現(xiàn)為顯著的正相關(guān)關(guān)系。在東太平洋，2004—2014年間SST呈現(xiàn)負(fù)異常，則SDP隨SST的降低而減小，而2014年至今主要表現(xiàn)為正異常，隨著SST的升高，SDP呈現(xiàn)增大趨勢(shì)。與東太平洋恰恰相反，西太平洋的SST在2004—2014年呈現(xiàn)正異常，2014年至今呈現(xiàn)負(fù)異常，因此，在前一時(shí)段，SDP隨SST的升高而增大，后一時(shí)段，SDP隨SST的下降而減小。南半球海區(qū)的變化情況與西太平洋大致相同。北印度洋的情況如前文所述，由于有效數(shù)據(jù)量少，不進(jìn)行長(zhǎng)期趨勢(shì)分析。該海區(qū)的最大相關(guān)系數(shù)出現(xiàn)在2017年，為-0.94，在95%的置信水平上顯著。

圖12 2004—2019年全球海域平均SDP與SST相關(guān)系數(shù)的年際分布（虛線為0.05置信水平）

通過(guò)上述對(duì)各指數(shù)與SST的相關(guān)分析可知，盡管在16 a間有少數(shù)年份顯示出顯著負(fù)相關(guān)，以及北大西洋SDP在2014年由正相關(guān)轉(zhuǎn)變?yōu)樨?fù)相關(guān)，但整體上指數(shù)與SST之間仍呈現(xiàn)為普遍的顯著正相關(guān)。這種相關(guān)性說(shuō)明，在全球變暖的背景下，TC的破壞力會(huì)隨著SST的升高而增大。此外，在全球海域中，北印度洋的情況較為特殊，該海區(qū)的PDI與SST在1969—1989年為負(fù)相關(guān)，1990年以后轉(zhuǎn)變?yōu)檎嚓P(guān)，而WDP和SDP由于相關(guān)系數(shù)的有效數(shù)據(jù)量少，無(wú)法進(jìn)行長(zhǎng)時(shí)間趨勢(shì)分析，兩個(gè)指數(shù)與SST的最大相關(guān)系數(shù)均為負(fù)，出現(xiàn)在2017年。

5 結(jié)論和討論

本文基于IBTrACS的歷史TC資料和HadISST數(shù)據(jù)集的SST數(shù)據(jù)，對(duì)1949—2019年全球海域TC的能量耗散指數(shù)PDI、風(fēng)破壞力評(píng)級(jí)指數(shù)WDP和風(fēng)暴潮與波浪破壞力評(píng)級(jí)指數(shù)SDP進(jìn)行統(tǒng)計(jì)，揭示了近70 a TC破壞力的時(shí)空變化規(guī)律，并進(jìn)一步研究了各指數(shù)與全球SST變化之間的相關(guān)關(guān)系，得到的主要結(jié)論如下：

（1）北大西洋、東太平洋和西太平洋的累積能量耗散、風(fēng)破壞力和風(fēng)暴潮與波浪破壞力主要集中在7—10月，峰值在8—9月，南半球海區(qū)集中在1—3月，峰值在3月。北印度洋的累積破壞力呈現(xiàn)雙峰型分布，活躍期為4—6月和10—11月，峰值分別在5月和10—11月。在北半球海域中，破壞力月際累積最大的是西太平洋，其次為東太平洋和北大西洋，北印度洋為全球海域最小。

（2）全球各海區(qū)WDP和SDP平均評(píng)級(jí)與累積評(píng)級(jí)的月際分布特征不同，平均評(píng)級(jí)具有更長(zhǎng)的活躍期，這是因?yàn)槠骄u(píng)級(jí)凸顯了TC生成數(shù)量少但破壞力大的月份。此外，還有一些月份風(fēng)破壞力小，但風(fēng)暴潮與波浪破壞力大，這說(shuō)明若僅用頻數(shù)或強(qiáng)度進(jìn)行評(píng)級(jí)，易造成預(yù)期之外的破壞。

（3）從年際變化的長(zhǎng)期趨勢(shì)上看，全球各海區(qū)累積PDI除西太平洋均表現(xiàn)為顯著的線性增長(zhǎng)趨勢(shì)。從年際分布上看，各海區(qū)PDI顯示出明顯的周期性變化，北大西洋和南半球存在3 a左右的變化周期，東太平洋為2 a，北印度洋為7 a，西太平洋具有3 a左右的年際變化周期和35 a左右的年代際變化周期。平均WDP和SDP年際變化緩慢，波動(dòng)不大，且未在全球各海區(qū)表現(xiàn)出較為一致的變化趨勢(shì)。

（4）TC破壞力變化與其生成海域的SST之間具有密切聯(lián)系。PDI與北大西洋、東太平洋、西太平洋和南半球SST表現(xiàn)為顯著的強(qiáng)相關(guān)，北印度洋在1990年由負(fù)相關(guān)轉(zhuǎn)變?yōu)榱苏嚓P(guān)。WDP和SDP在全球海域普遍表現(xiàn)為正相關(guān)。該結(jié)論表明，TC破壞力將在全球變暖的背景下顯示出顯著增長(zhǎng)的趨勢(shì)。

另外值得注意的是，在對(duì)指數(shù)進(jìn)行年際變化分析時(shí)，只有PDI在全球海域中表現(xiàn)出了較為一致的上升趨勢(shì)，而WDP和SDP在各海區(qū)的年際變化情況是相當(dāng)不一致的。以前的TC資料忽略了TC的非對(duì)稱性結(jié)構(gòu)，風(fēng)半徑值自2003年起才開始由對(duì)稱轉(zhuǎn)為按象限記錄，且過(guò)去用于計(jì)算WDP和SDP的R33觀測(cè)記錄也比較少[25]，故本文對(duì)WDP和SDP的研究時(shí)間范圍選擇了有完整記錄的2004—2019年，但時(shí)間范圍的縮短意味著統(tǒng)計(jì)結(jié)果在長(zhǎng)時(shí)間線性趨勢(shì)上的代表性較差，與實(shí)際情況之間可能存在較大的誤差。對(duì)TC影響指數(shù)的研究結(jié)果也顯示，國(guó)內(nèi)外對(duì)ACE開展的統(tǒng)計(jì)工作最多，關(guān)于PDI的研究在近幾年開始呈現(xiàn)發(fā)展態(tài)勢(shì)，而關(guān)于WDP和SDP的研究始終很少，風(fēng)半徑數(shù)據(jù)的缺失在一定程度上解釋了這一趨勢(shì)的產(chǎn)生原因。因此，未來(lái)關(guān)于WDP和SDP統(tǒng)計(jì)分析的研究發(fā)展，亟需多年歷史風(fēng)半徑數(shù)據(jù)的支持。

致謝：感謝江蘇省氣候變化協(xié)同創(chuàng)新中心對(duì)本研究的資助。