Marine Meteorology

Climatic characteristics and interannual variability of tropical cyclone rapid intensification in the South Indian Ocean

  • WANG Hao , 1, 2 ,
  • WANG Jing , 1 ,
  • ZHENG Jiayu 2
Expand
  • 1. School of Geography and Planning, Sun Yat-sen University, Guangzhou 510275, China
  • 2. State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China
WANG Jing. email:

Copy editor: YAO Yantao

Received date: 2021-02-20

  Revised date: 2021-05-04

  Online published: 2021-05-10

Supported by

National Key Research and Development Program of China(2019YFA0606701)

National Natural Science Foundation of China(41731173)

Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou)(GML2019ZD0306)

National Program on Global Change and Air-Sea Interaction under grant(GASI-IPOVAI-03)

Innovation Academy of South China Sea Ecology and Environmental Engineering, the Chinese Academy of Sciences(ISEE2021ZD01)

Abstract

Using tropical cyclone (TC) data from the Joint Typhoon Warning Center, the Extended Reconstructed Sea Surface Temperature version 5 from the National Oceanic and Atmospheric Administration, the temperature and salinity data from the Global Ocean Data Assimilation System, and the reanalysis data from the National Centers for Environmental Prediction/National Center for Atmospheric Research, we investigated climatic characteristics and interannual variability of TC rapid intensification (RI) in the South Indian Ocean (SIO) during the period of 1981~2019. The occurrence frequency of RI events in the SIO presents a single peak distribution, which ranges from October to June of the next year. The occurrence of SIO RI events is distributed in the subtropical area in a belt shape, with three high-value centers of RI events located in the northeast of Madagascar, the middle of the SIO and the northwest of Australia. The spatial distribution of SIO RI events is mainly modulated by two large-scale environmental variables: tropical cyclone heat potential (TCHP) and vertical wind shear (VWS). We also found that the influence of El Niño-Southern Oscillation on the occurrence frequency of SIO RI events is asymmetric: the occurrence frequencies of RI events during El Niño and La Niña years are both reduced. During El Niño years, the decrease of RI events is mainly due to higher VWS. During La Niña years, the decrease of RI events is mainly due to lower TCHP, and SST; VWS and RHUM also contribute to the decrease of RI events.

Cite this article

WANG Hao , WANG Jing , ZHENG Jiayu . Climatic characteristics and interannual variability of tropical cyclone rapid intensification in the South Indian Ocean[J]. Journal of Tropical Oceanography, 2022 , 41(1) : 94 -105 . DOI: 10.11978/2021021

热带气旋(Tropical Cyclone, TC)是生成在热带暖洋面上具有暖心结构的强烈气旋性环流, 其带来的强风和强降水对所经过地区的人类生命和财产安全造成了巨大的威胁(Zhang et al, 2009; 张娇艳 等, 2011)。因此, 准确地预测TC的生成及发展至关重要。目前, 在TC路径预报方面, 前人已经做过大量研究, 取得了显著的进展。Chan(2005)指出对于西北太平洋TC的生成高峰季(7—9月), 西北太平洋TC在El Niño年倾向于弯曲的移动路径, 在La Niña年倾向于在偏北位置形成, 并向北移动; 对于10—12月来说, 西北太平洋TC在El Niño年仍然倾向于弯曲的移动路径, 但在La Niña年以西行路径为主。McAdie等(2000)利用持久性模型并结合大西洋TC的气候特征, 提出了一种新的TC路径预测方法, 减少了之前路径预测存在的诸多误差。然而, TC强度预测的准确性要远低于路径预测, 对TC强度特别是快速增强(rapid intensification, RI)事件的准确预测仍旧是亟待解决的难题之一(Elsberry et al, 2007; Rappaport et al, 2009; DeMaria et al, 2014)。无论是在北大西洋或是西北太平洋, 强TC在其完整生命周期中都经历了至少一次的RI事件(Kaplan et al, 2003; Wang et al, 2008)。因此, 增进对RI事件变异和物理机制的认识对于理解TC强度变异和提高预测准确度均具有重要的意义。
关于西北太平洋、北大西洋、北印度洋RI事件变率的研究已有很大的进展(Wang et al, 2008, 2015, 2017; Zhao et al, 2018; Gao et al, 2020)。在西北太平洋TC生成高峰季(7—9月)期间, 厄尔尼诺年的RI事件比例(53%)远高于拉尼娜年同期(37%)(Wang et al, 2008)。Wang等(2015)的研究则表明, 与TC密切相关的大尺度环境因子, 如海表面温度(sea surface temperature, SST)、TC热潜(tropical cyclone heat potential, TCHP)、垂直风切变(vertical wind shear, VWS)、相对湿度(relative humidity, RHUM)等, 对RI事件的发生有着显著的影响。在北大西洋, Kaplan(2003)等指出产生RI事件的TC相较于没有产生RI事件的, 大多形成于偏南及偏西的区域, 且有西移的趋势, 并提出了一种估算RI事件概率的简单方法。Wang等(2017)的研究发现, 较低的VWS和较高的TCHP是北大西洋TC主要发展区RI事件多发的首要原因, 而北大西洋东部RI事件主要受RHUM所调控。此外, 该研究强调大尺度的大气和海洋变率对北大西洋TC的生成并不重要, 但在热带低压形成之后, 对TC强度的影响则十分明显。Ng等(2020)的研究进一步指出, 在热带北大西洋中心区域RI事件的年际变率受季节性大尺度环境变量影响, 而墨西哥湾和加勒比海则不存在这种关系, 主要由于前一区域VWS的季节异常与潜在强度之间表现出强烈的负相关。在北印度洋, 孟加拉湾的RI事件在拉尼娜年的产生频率远高于厄尔尼诺年的产生频率, 且RHUM和VWS是该海盆调制RI事件变化的主要环境变量(Girishkuma et al, 2015)。
以往的研究大多集中于西北太平洋、北大西洋等海域, 对于印度洋特别是南印度洋TC强度变化的研究是相对缺乏的。已有的关于印度洋TC的研究也主要是侧重于对TC的生成频率和登陆频率, 以及受大气波动影响等方面的分析(Bessafi et al, 2006; Mavume et al, 2009; Kuleshov et al, 2010), 关于RI事件的研究则更加匮乏。本文拟对南印度洋RI事件的气候特征和年际变率进行研究, 分析大尺度环境变量(SST、TCHP、VWS和RHUM)对南印度洋RI事件气候特征和年际变率的影响。

1 数据和方法

1.1 数据

本文的TC数据采用1981—2019年联合台风预警中心(Joint Typhoon Warning Center, JTWC)的TC最佳路径数据集(https://www.usno.navy.mil/JTWC), 该数据集包含了历次TC每隔6h的中心位置、中心最大持续风速、中心最低气压等信息。由于1981年以前缺少卫星数据的修正, JTWC数据集中南印度洋部分的资料存在较大缺失, 故本文研究的时间从1981年开始。1981年至今, 数据的各项参数均记录正常, 但也存在部分TC生成点位数据的重复记录, 在本文的研究过程中对这些重复记录数据已进行剔除。根据萨菲尔-辛普森风力等级量表, 最大中心持续风速为10~17m•s-1(20~33kn)的TC被归类为热带低压; 若热带低压达到18m•s-1(34kn)的风速, 则被列为热带风暴; 风速继续增强直至33m•s-1(64kn)则被归类为飓风。飓风又被划分为1~5级, 对应的最大中心持续风速分别为33~42m•s-1(64~82kn)、43~49m•s-1(83~95kn)、50~58m•s-1(96~112kn)、59~69m•s-1(113~136kn)和70m•s-1(137kn)。
在前人的研究中, 通常将TC间隔24h最大中心持续风速的变化(ΔV24=V24h-V0h, 其中V0h代表TC事件中某个记录时刻的风速值, V24h代表该TC经过24h后新记录时刻的风速值)的所有样本的第95个百分位数作为定义RI事件的阈值(Kaplan et al, 2003)。在本研究中, 南印度洋所有样本ΔV24的第95个百分位数为15.4m•s-1(30kn)。因此, 在一次TC事件中, 若某个数据点相对于24h前的最大中心持续风速增量ΔV24大于15.4m•s-1, 便将该数据点定义为一个RI事件。
参照前人的研究(Wang et al, 2015, 2017), 对影响RI事件产生的大环境变量主要选取了4个, 分别为SST、TCHP、VWS和RHUM。其中, 2个海洋变量SST和TCHP使用的数据如下: SST数据采用美国国家海洋和大气管理局(National Oceanic and Atmospheric Administration, NOAA)的扩展重构月平均海表温度数据第5版本(the Extended Reconstructed Sea Surface Temperature version 5) (https://psl.noaa.gov/data/gridded/data.noaa.ersst.v5.html), 水平分辨率为2°×2°; 计算TCHP变量使用全球海洋数据同化系统(Global Ocean Data Assimilation System)的月平均温度和盐度数据(https://psl.noaa.gov/data/gridded/data.go-das.html), 水平分辨率为1°×0.33°, 所有变量数据层数均为40层。VWS和RHUM这2个大气变量使用美国国家环境预报中心/国家大气研究中心(the National Centers for Environmental Prediction/National Center for Atmospheric Research)的再分析数据(https://psl.noaa.gov/data/gridded/data.ncep.reanalysis.pressure.html)得到, 分辨率为2.5°×2.5°。上述所有数据的时间跨度均为1981—2019年。
海洋尼诺指数(Oceanic Niño Index, ONI)是NOAA定义的热带太平洋厄尔尼诺(暖)和拉尼娜(冷)事件的标准(https://ggweather.com/enso/oni.htm), 它表示Niño3.4区域(即5°N—5°S, 120°—170°W)3个月的滑动平均SST异常。厄尔尼诺事件定义为连续5个月ONI大于或等于0.5℃, 拉尼娜事件定义为连续5个月ONI小于或等于-0.5℃。SST异常阈值又进一步细分为弱(0.5~0.9℃)、中等(1.0~1.4℃)、强(1.5~1.9℃)和非常强(≥2.0℃)事件。NOAA的报告中指出, 对于被归类为弱、中等、强或非常强的事件, 必须至少连续3个月达到或超过阈值。

1.2 相关计算方法

1.2.1 TCHP计算

TCHP的概念最早来源于飓风热潜(Hurricane Heat Potential), 后来的研究让其概念得到延伸(Leipper et al, 1972; Goni et al, 2007), TCHP代表的是26℃等温层面至海表面之间暖水所包含的海洋热含量(Leipper, 1967; Leipper et al, 1972), 其计算公式如下:
$\text{TCHP=}{{c}_{p}}\int_{{{D}_{26}}}^{0}{\rho (z)[T(z)-26]dz}$
式中, cp表示恒定压力下的比热, 固定值为3.9kJ•kg-1•K-1, D26表示26℃等温线的深度, ρ(z)代表实测密度, T(z)表示实测温度。

1.2.2 相关分析及其检验

线性相关分析是用来衡量两个变量之间线性相关程度的量, 一般用相关系数r表示。r的取值范围为[-1, 1], r的绝对值越大, 说明两个变量之间的线性相关程度越大。当r > 0时, 表明两个变量呈正相关; 当r < 0时, 表明两个变量呈负相关。相关系数r的计算公式为:
$r\text{(}X\text{, }Y\text{) = Cov(}X\text{, }Y\text{)}\cdot {{\left( \sqrt{Var(X)\cdot Var(Y)} \right)}^{-1}}$
式中, Cov(X, Y)为XY的协方差, Var(X)为X的方差, Var(Y)为Y的方差。
相关性检验是对变量之间是否相关以及相关的程度所进行的统计检验, 一般基于t分布进行检验。在假设总体相关系数ρ=0成立的条件下, 相关系数r的概率密度函数正好是t分布的密度函数, 可以用t检验对r进行显著性检验。检验统计量为:
$t \text{= }r\cdot \left( \sqrt{n-2} \right)\cdot {{\left( \sqrt{1-{{r}^{2}}} \right)}^{-1}}$
式中, t为检验统计量, n为样本量, 统计量遵从自由度(n-2)的t分布。再根据给定显著性水平的αt分布表, 若t > tα(tα代表显著性水平为α的统计量的值), 则拒绝原假设, 认为相关系数是显著的。

1.2.3 合成分析

合成分析是气象学分析与预报中经常使用的方法, 是将两种不同特征或者状态的气象变量进行合成, 即求不同状态下的某气象要素的平均值, 以揭示各状态之间是否存在明显差异(施能, 2009)。
本文在合成分析中所使用的显著性检验方法为t检验(t test),在总体方差σ2未知的情况下, 利用样本方差s2构造检验总体均值的t统计量:
$t \text{= (}\bar{x}-{{\mu }_{0}}\text{)}\cdot {{s}^{-1}}\cdot \sqrt{n}$
式中, $\bar{x}$s分别代表样本均值和标准差, μ0为总体均值, n为样本量。在确定显著性水平α之后, 根据自由度v=n-1查t分布表, 若 $\left| \text{ }t\text{ } \right|$tα, 则认为通过了显著性水平为α的检验(魏凤英, 1999)。

2 结果分析

2.1 南印度洋RI事件的气候特征

2.1.1 时空分布特征

在1981—2019年间, 南印度洋共生成了602个TC, 其中热带低压17个。根据萨菲尔-辛普森风力等级量表, 热带低压的最大中心持续风速小于33kn(17m•s-1), 且在最佳路径资料中, 一个TC的起始记录风速最小不低于15kn(约7.1m•s-1)。因此, 无论在实际意义上还是RI的定义上, 热带低压都没有达到研究的标准, 故在本文研究中将热带低压剔除, 只选取达到热带风暴级别以上的TC进行研究。在剩余的585个TC中, 共有215个TC产生了RI事件, 在这些TC的完整周期中共产生了870个RI事件, 产生RI事件的TC占南印度洋TC生成总数的37%。在这215个TC中, 平均每个TC产生了4个RI事件, 其中1989年编号为28号的TC产生了最多的RI事件(10个)。
图1显示了南印度洋1981—2019年TC生成数量和RI事件产生数量的逐月分布特征。从图1a中看可出南印度洋TC生成的高峰季为11月至次年4月, 此时间范围内生成的TC数量占南印度洋全年TC生成总数量的81%; 而RI事件的主要产生季节是12月至次年4月(图1b), 占南印度洋全年RI事件总产生数量的92%。从图中还可以看出, 从10月(南半球夏末)开始, RI事件的数量逐渐增加, 在次年3月(南半球秋季)达到全年的最大值, 数量为203个, 4—6月开始逐步降低, 每年的7—9月则没有RI事件产生。图1c展示了逐月平均的RI事件产生比率(RI事件产生比率=RI事件产生数量/TC生成数量), 对比TC数量和RI事件产生比率的全年分布, 可以发现两者的峰值月份没有同步, 这一特征不同于其他海盆如北大西洋和西北太平洋, 这两个海盆的TC数量和RI事件产生数量在逐月变化上是高度同步的。控制TC生成与RI事件的环境变量类似, 但也有着不同之处, TC生成与RI事件产生的峰值月份不同, 可能是不同海盆主导的环境变量不同所致。南印度洋RI事件主要产生在12月至次年4月, 本研究将12月至次年4月这5个月定义为南印度洋RI事件的活跃季节, 下文将主要研究该时间段RI事件的分布特征及其年际变率。
图1 1981—2019年南印度洋TC生成数量(a)、RI事件产生数量(b)和RI事件产生比率的逐月分布特征(c)

Fig. 1 Monthly distribution of (a) TC genesis number, (b) occurrence number of RI events and (c) occurrence ratio of RI events in the SIO averaged from 1981 to 2019

为了探究南印度洋RI事件产生的空间分布特征, 本文将南印度洋(20°—130°E, 0°—40°S)划分为8×22个空间分辨率为5°×5°的网格, 再对每个网格中的RI事件进行统计。图2显示了南印度洋1981—2019年RI事件的空间分布特征, 可以看出南印度洋RI事件主要产生在10°—20°S的狭长带状区域, 并且存在多个中心, 较为显著的中心有3个, 分别是马达加斯加岛东北海域(66°E, 11°S)、南印度洋中部海域(85°E, 11°S), 以及澳大利亚西北海域(120°E, 18°S)。这3个中心所产生的RI事件数量约占南印度洋RI事件产生总数的67%。其中, RI事件最大值中心区为马达加斯加岛东北海域, 在以此为中心的高值区共产生了205个RI事件, 平均每年5.39个。
图2 1981—2019年南印度洋RI事件产生数量的平均空间分布

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作

Fig. 2 Spatial distribution of RI events in the SIO averaged from 1981 to 2019

此外, 本文分析了南印度洋TCRI事件年际变率的方差分布, 相应的分布图见图3。由图可知, 南印度洋TCRI事件数量的方差分布也主要呈带状分布在南印度洋副热带地区, 方差分布图中的3个高值中心与图2所示的TCRI事件的3个高值中心相互吻合。这表明TCRI事件产生数量较高的区域, 其TCRI事件变化程度也较大。
图3 1981—2019年南印度洋RI事件产生数量的方差分布

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作

Fig. 3 Variance distribution of TCRI events in the SIO from 1981 to 2019

2.1.2 大尺度环境变量对南印度洋RI事件空间分布的影响

海洋到大气的热通量是TC的主要能量来源。强劲的表面风和TC中持续下降的表面压力推动了热通量的传播(Schade, 2000)。SST和海洋热含量的相对细微变化也可导致TC的强度短时间内在萨菲尔-辛普森等级量表上变化数个类别(Sun et al, 2007), TCHP对TC强度的变化同样也有着不可忽视的影响, 在全球范围内TC生成数最多的西北太平洋, RI事件大多产生在TCHP相对较高的区域(Wada et al, 2007)。VWS也是影响TC强度的一个重要因子, VWS可以改变TC的内部位涡及相应的热力环境来影响其强度变化, 对流层VWS越小, 越有利于热力气柱的形成和CISK的运转, 从而引起RI事件的产生(DeMaria, 1996; Walsh et al, 2000; Emanuel et al, 2004)。此外, 充足的水汽可以为TC的发展提供潜热(Wu et al, 2012; 韩翔 等, 2018), 进而TC强度增加而形成RI事件, 所以RHUM也可能影响RI事件的产生。
图4显示了在南印度洋RI事件活跃季节SST、TCHP、VWS和500-hPa RHUM的气候态分布。从图4a可以看出南印度洋RI事件产生区域, 特别是马达加斯加岛东北海域、南印度洋中部海域和澳大利亚西北海域这3个RI事件高值区域, 大部分位于SST接近30℃的区域, 较高的SST可以促进TC强度在短时间内提升, 进而引起RI事件的产生。图4b显示了南印度洋TCHP的空间分布, TCHP在西南印度洋较低, 这与前人指出的西南印度洋温跃层较浅相一致(Du et al, 2009)。在RI事件空间分布中, 除马达加斯加岛东北海域这一RI事件高值中心的TCHP呈现出较低值外, 其余两个高值中心都位于TCHP值较高的海域, 这说明较高的TCHP可能也是这两个区域RI事件产生较多的环境因素。图4c显示了南印度洋VWS的空间分布, 可以看出马达加斯加岛东北海域和南印度洋中部海域这两个高值中心均位于VWS最低值的中心区域, 而澳大利亚西北海域这一RI事件高值中心的VWS则高于前两个中心。在图4d显示的500-hPa RHUM的空间分布中, 3个RI事件高值中心的RHUM均较高, 但与RI事件分布的联系相比于其他因子均较弱。
图4 1981—2019年南印度洋SST(a)、TCHP(b)、VWS(c)及500-hPa RHUM(d)的气候态

Fig. 4 Climatology of (a) SST, (b) TCHP, (c) VWS, and (d) 500-hPa RHUM in the SIO from 1981 to 2019

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作。图中黑色等值线为1981—2019年南印度洋RI事件数量的平均空间分布综上所述, 南印度洋RI事件产生区域大致位于SST接近30℃的区域, 而控制RI事件高值中心产生的大尺度环境因子主要是TCHP和VWS, 但3个高值中心的主要影响因子略有不同。在马达加斯加岛东北海域RI事件高值中心, 主要是气候态上较弱的VWS导致了RI事件频繁产生; 在南印度洋中部海域这一中心区域附近, 较高的TCHP以及较弱的VWS对RI事件产生均有贡献; 而在澳大利亚西北海域, RI事件的频繁产生主要是较高的TCHP所致。

2.2 RI事件年际变化特征

2.2.1 RI事件的年际变率

图5显示了南印度洋TC数量及RI事件产生数量的时间序列。在1981—2019年间, 南印度洋TC生成最少的年份是2011年, 数量为8个, TC生成最多的年份是1985年, 数量为20个, 其主要呈现年际变化, 未出现较为明显的下降或上升趋势。而在南印度洋RI活跃季节内RI事件的变化趋势则与TC的变化趋势有着显著的差异(图5), 可以看出南印度洋RI事件产生数量呈现年际和年代际的变化特征, 并且总体上呈现上升的趋势。在此期间, 产生RI事件最少的年份是1982年, 数量为0个, 产生RI事件最多的年份是2014年, 数量为44个。
图5 1981—2019年南印度洋快速增强(RI)事件数量(黑线)和热带气旋(TC)数量(红线)的时间序列

Fig. 5 Time series of RI events number (black) and TC number (red) in the SIO from 1981 to 2019

在年际时间尺度上, ENSO是现有的观测记录中热带气候主导模态(Xie et al, 2009)。前人的研究发现, 印度洋气候(包含海表温度)主要受ENSO的强迫, 其自身的自然变率较弱(Xie et al, 2002, 2009; Zhang et al, 2015)。图6显示了去趋势后RI事件数量异常和12月至次年4月的Niño3.4 (5°S—5°N, 120°—170°W)指数的时间序列, 对二者进行相关性分析, 相关系数为0.2, 未通过相关性显著检验。但是从图6可以看出, 当Niño3.4指数较高, 即厄尔尼诺事件发生时, 南印度洋RI事件产生数量较低; 而当Niño3.4较低, 即拉尼娜事件发生时, 绝大多数年份中南印度洋RI事件产生数量同样较少, 但也存在个别年份RI事件产生数量呈正异常。虽然从整体来看ENSO事件对南印度洋RI事件产生数量的影响不显著, 但这可能是由于影响的不对称造成的。为了更具体地考察ENSO事件对南印度洋热带气旋RI事件的影响, 下文将通过合成分析加以探讨。
图6 1981—2019年去趋势后的南印度洋快速增强(RI)事件异常值和Niño3.4指数的时间序列

Fig. 6 Time series of detrended RI events number anomaly (black) and Niño3.4 index (red) from 1981 to 2019

前人的研究指出, ENSO虽然是热带主导模态, 但是需要一定强度的ENSO异常事件才能强迫出印度洋的响应(Yu et al, 1999)。所以在本研究中, 为了突出ENSO信号, 选取ENSO事件强度等级为强以上的年份进行合成。根据ONI指数, 将厄尔尼诺事件强度等级为强和超强的年份及拉尼娜事件强度等级为强的年份挑选出来(见表1), 并分别与RI事件产生数量的异常进行合成(图7)。可以发现, 与图6中的结论类似, ENSO对南印度洋RI事件产生数量的调制是非对称的。在强和超强厄尔尼诺事件发生的年份(图7a), 南印度洋RI活跃季节的RI事件产生数量在整个海盆上呈现显著的负异常。这些年份出现的负异常值较为均匀地分布在南印度洋海盆中部10°S附近, 在60°E、75°E和85°E这3个区域出现较高的负异常值, 而在这条负异常线的两侧分别出现了两个小面积的正异常区域(95°E、5°S和85°E、15°S)。在强拉尼娜事件发生的年份(见图7b), 南印度洋RI活跃季节的RI事件产生数量同样呈现出大面积的负异常区域, 负异常值的3个中心点与厄尔尼诺事件发生年份的负异常中心位置大致相同, 不同的是在强拉尼娜事件发生年份澳大利亚西北海域存在一个数值较高的正异常中心。
表1 1981—2019年厄尔尼诺和拉尼娜事件发生的年份及对应强度

Tab. 1 El Niño and La Niña years and intensities from 1981 to 2019

厄尔尼诺事件 拉尼娜事件
中等 超强 中等
2004—2005 1986—1987 1987—1988 1982—1983 1983—1984 1995—1996 1988—1989
2006—2007 1994—1995 1991—1992 1997—1998 1984—1985 2011—2012 1998—1999
2014—2015 2002—2003 2015—2016 2000—2001 1999—2000
2009—2010 2005—2006 2007—2008
2008—2009 2010—2011
图7 南印度洋RI事件产生数量异常在强和超强厄尔尼诺年(a)以及强拉尼娜年(b)的合成

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作。图中白色“+”表示异常值通过90%置信度的t检验

Fig. 7 Composite anomalies of RI events number in the SIO in the years of strong and super strong El Nino events (a) and strong La Nina events (b)

2.2.2 大尺度环境因子对南印度洋RI事件年际变化的影响

为了考察ENSO对南印度洋RI事件的调制, 接下来将与RI事件联系密切的4个大尺度的环境因子(SST、TCHP、VWS、RHUM)在强和超强厄尔尼诺事件发生的年份和强拉尼娜事件发生的年份进行合成。如图8a所示, 在厄尔尼诺事件发生的年份, 印度洋的SST呈现显著的正异常, 这与前人的结论相一致(Xie et al, 2002, 2009; Du et al, 2009)。在厄尔尼诺发展时期, 海洋性大陆(Maritime Continent, MC)会产生沃克环流下沉异常, 下沉气流会在印度洋强迫暖的海洋罗斯贝波向西传播, 进而导致西南印度洋增暖, 加深西印度洋的温跃层(Xie et al, 2002)。相应地, 在厄尔尼诺年西南印度洋TCHP显著增加(图8b), 南印度洋VWS(图8c)和RHUM(图8d)均为显著正异常。虽然SST、TCHP、RHUM等因子有利于RI事件的产生, 但是RI事件产生数量在整个南印度洋洋盆中心区域的大部分海域仍表现为连续的纬度较为集中的负异常。而在对应区域仅VWS存在较高的正异常, 不利于RI事件的产生。因此, 在强和超强的厄尔尼诺事件发生年份, VWS是调制南印度洋RI事件异常最主要的环境因子, 因为它削弱了潜热释放对热带气旋风速的增强作用(潘静 等, 2010)。
图8 南印度洋SST(a)、TCHP(b)、VWS(c)和500-hPa RHUM(d)异常值在强及超强厄尔尼诺年的合成

Fig. 8 Composite anomalies of (a) SST, (b) TCHP, (c) VWS, and (d) RHUM in the SIO in the years of strong and super strong El Niño events

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作。图中黑色等值线表示南印度洋RI事件产生数量异常值在强和超强厄尔尼诺年的合成, 白色“+”表示SST、TCHP、VWS和RHUM异常值通过90%置信度的t检验与厄尔尼诺年相反, 拉尼娜年南印度洋SST(图9a)和TCHP(图9b)呈现出负异常, 并且以RI事件产生数量负异常中心处的TCHP负异常值最强。另外可以看到, 在部分区域VWS正异常(图9c)与RHUM(图9d)负异常对拉尼娜年南印度洋RI事件产生较少也起到一定的贡献作用。在澳大利亚的西北海域存在RI事件产生数量的正异常, 对应的3个热力学因素(SST、TCHP和RHUM)也都分别呈现出一定的正异常。因此, 不同于厄尔尼诺年, 拉尼娜年南印度洋RI事件产生数量减少主要是因为TCHP变弱, 此外SST、VWS、RHUM也有一定的贡献。
图9 南印度洋SST(a)、TCHP(b)、VWS(c)和500-hPa RHUM(d)异常值在强拉尼娜年的合成

Fig. 9 Composite anomalies of (a) SST, (b) TCHP, (c) VWS, and (d) RHUM in the SIO in the years of strong La Niña events

该图基于国家测绘地理信息局标准地图服务网站下载的审图号为GS(2016)1665号的标准地图制作。图中黑色等值线表示南印度洋RI事件产生数量异常值在强拉尼娜年的合成, 白色“+”表示SST、TCHP、VWS和RHUM异常值通过90%置信度的t检验以上的现象表明了在年际尺度上, ENSO对南印度洋RI事件的影响是不对称的, 即厄尔尼诺年与拉尼娜年南印度洋RI事件产生数量均减少。但是使其减少的主要因素在厄尔尼诺年与拉尼娜年却是不同的, 厄尔尼诺年南印度洋较少RI事件产生主要是由较高的VWS导致, 而拉尼娜年南印度洋较少RI事件产生主要是由于TCHP的降低。

3 结论

本文系统地研究了1981—2019年南印度洋RI事件的气候特征和年际变率, 通过对南印度洋的RI事件产生数量和大尺度环境变量进行分析, 主要得到以下主要结论:
1) 南印度洋RI事件产生数量的逐月变化呈现出单峰性, 主要产生在10月至次年6月这9个月(南半球夏季和秋季), 而7—9月这3个月(南半球冬季和春季)则无RI事件产生。
2) 南印度洋RI事件的空间分布主要集中于5°—15°S这一纬度带上, 且出现了3个高值中心, 分别是马达加斯加岛东北海域(66°E, 11°S)、南印度洋中部海域(85°E, 11°S)以及澳大利亚西北海域(120°E, 18°S)。其中, 马达加斯加岛东北海域这一高值中心所产生的RI事件数量最多, 换而言之, 是整个南印度洋最容易产生RI事件的区域, 且此中心靠近马达加斯加岛。马达加斯加岛(尤其是其东北部)常年遭受TC的侵扰, 有不少TC途经此中心区域而产生强度增强, 进而对该地区人民的生命和财产安全造成威胁。
3) 气候态上, 南印度洋RI事件产生数量的空间分布主要受到TCHP和VWS这两个大尺度环境变量调制, 且不同的RI高值中心区域所主导调制的因子不同。在马达加斯加岛东北海域RI事件高值中心, 主要是气候态上较弱的VWS导致了RI事件频繁产生; 而在南印度洋中部海域这一中心区域附近, 较高的TCHP以及较弱的VWS对RI事件的产生均有贡献; 在澳大利亚西北海域, RI事件的频繁产生主要是较高的TCHP所致。
4) 在年际变率上, ENSO对南印度洋RI事件的影响是不对称的, 即在厄尔尼诺年与拉尼娜年RI事件的发生数量均呈现减少的状况, 但使其减少的主要因子却截然不同。在厄尔尼诺年, 过高的VWS导致了RI事件的减少; 而在拉尼娜年, RI事件的减少主要是由于TCHP降低所致, 其次SST、TCHP、RHUM对RI事件的减少也有贡献。
本文研究了ENSO事件对南印度洋RI事件的影响及机制, 但不排除其他海温模态对南印度洋RI事件的可能影响。比如, 前人的研究指出ENSO可以通过影响沃克环流来改变海洋性大陆的对流运动(Leung et al, 2016, 2017, 2019; Zhang et al, 2019), 而ENSO对海洋性大陆对流活动的影响受到南印度洋偶极子(Southern Indian Ocean Dipole, SIOD)的调控(Gong et al, 2019; Xiao et al, 2020)。这些都表明南印度洋RI事件的产生可能还受到其他海温模态的调制作用。目前本文主要从环境因子的角度出发, 分析了ENSO对南印度洋RI事件的影响, 而其他海温模态对南印度洋RI事件的影响将在下一步工作中开展深入的研究与探讨。

*感谢3位审稿人认真的审阅和建设性的意见, 这对提高论文水平很有帮助。

[1]
韩翔, 赵海坤, 孙齐, 2018. 夏季热带大气准双周振荡对西北太平洋台风生成的影响[J]. 热带气象学报, 34(4): 524-534.

HAN XIANG, ZHAO HAIKUN, SUN QI, 2018. Effects of Quasi-Biweekly Oscillation of tropical atmosphere on typhoon formation in Northwest Pacific[J]. Journal of Tropical Meteorology, 34(4): 524-534. (in Chinese with English abstract)

[2]
潘静, 李崇银, 宋洁, 2010. 热带大气季节内振荡对西北太平洋台风的调制作用[J]. 大气科学, 34(6): 1059-1070.

PANG JING, LI CHONGYIN, SONG JIE, 2010. The modulation of Madden-Julian oscillation on typhoons in the northwestern Pacific Ocean[J]. Chinese Journal of Atmospheric Sciences, 34(6): 1059-1070. (in Chinese with English abstract)

[3]
施能, 2009. 气象统计预报[M]. 北京: 气象出版社: 235-236. (in Chinese)

[4]
魏凤英, 1999. 现代气候统计诊断与预测技术[M]. 北京: 气象出版社: 27-28. (in Chinese)

[5]
张娇艳, 吴立广, 张强, 2011. 全球变暖背景下我国热带气旋灾害趋势分析[J]. 热带气象学报, 27(4): 442-454.

ZHANG JIAOYAN, WU LIGUANG, ZHANG QIANG, 2011. Tropical cyclone damages in China under the background of global warming[J]. Journal of Tropical Meteorology, 27(4): 442-454. (in Chinese with English abstract)

[6]
BESSAFI M, WHEELER M C, 2006. Modulation of South Indian Ocean tropical cyclones by the Madden-Julian oscillation and convectively coupled equatorial waves[J]. Monthly Weather Review, 134(2): 638-656.

DOI

[7]
CHAN J C L, 2005. Interannual and interdecadal variations of tropical cyclone activity over the western North Pacific[J]. Meteorology and Atmospheric Physics, 89(1-4): 143-152.

DOI

[8]
DEMARIA M, 1996. The effect of vertical shear on tropical cyclone intensity change[J]. Journal of the Atmospheric Sciences, 53(14): 2076-2088.

DOI

[9]
DEMARIA M., SAMPSON C R, KNAFF J A, et al, 2014. Is tropical cyclone intensity guidance improving[J]. Bulletin of the American Meteorological Society, 95(3): 387-398.

DOI

[10]
DU YAN, XIE SHANGPING, HUANG GANG, et al, 2009. Role of air-sea interaction in the long persistence of el Niño-induced North Indian Ocean warming[J]. Journal of Climate, 22(8): 2023-2038.

DOI

[11]
ELSBERRY R L, LAMBERT T D B, BOOTHE M A, 2007. Accuracy of Atlantic and eastern North Pacific tropical cyclone intensity forecast guidance[J]. Weather and Forecasting, 22(4): 747-762.

DOI

[12]
EMANUEL K, DESAUTELS C, HOLLOWAY C, et al, 2004. Environmental control of tropical cyclone intensity[J]. Journal of the Atmospheric Sciences, 61(7): 843-858.

DOI

[13]
GAO JUN, ZHAO HAIKUN, KLOTZBACH P J, et al, 2020. Possible influence of tropical Indian Ocean sea surface temperature on the proportion of rapidly intensifying western North Pacific tropical cyclones during the extended boreal summer[J]. Journal of Climate, 33(21): 9129-9143.

DOI

[14]
GIRISHKUMAR M S, SUPRIT K, VISHNU S, et al, 2015. The role of ENSO and MJO on rapid intensification of tropical cyclones in the Bay of Bengal during October- December[J]. Theoretical and Applied Climatology, 120(3): 797-810.

DOI

[15]
GONG HAINAN, ZHOU WEN, CHEN WEN, et al, 2019. Modulation of the southern Indian Ocean dipole on the impact of El Niño-Southern Oscillation on Australian summer rainfall[J]. International Journal of Climatology, 39(4): 2484-2490.

DOI

[16]
GONI G J, KNAFF J, LIN I I, 2007. Tropical cyclone heat potential[J]. State of the Climate in, 43-45.

[17]
GRAY W M, 1979. Hurricanes: their formation, structure and likely role in the tropical circulation[M]//SHAW D B. Supplement to meteorology over the tropical oceans. Bracknell: James Glaisher House: 155-218.

[18]
KAPLAN J, DEMARIA M, 2003. Large-scale characteristics of rapidly intensifying tropical cyclones in the North Atlantic basin[J]. Weather and Forecasting, 18(6): 1093-1108.

DOI

[19]
KULESHOV Y, FAWCETT R, QI L, et al, 2010. Trends in tropical cyclones in the South Indian Ocean and the South Pacific Ocean[J]. Journal of Geophysical Research: Atmospheres, 115(D1): D01101.

[20]
LEIPPER D F, 1967. Observed ocean conditions and Hurricane Hilda, 1964[J]. Journal of the Atmospheric Sciences, 24(2): 182-186.

DOI

[21]
LEIPPER D F, VOLGENAU D, 1972. Hurricane heat potential of the Gulf of Mexico[J]. Journal of Physical Oceanography, 2(3): 218-224.

DOI

[22]
LEUNG M Y T, ZHOU WEN, 2016. Direct and indirect ENSO modulation of winter temperature over the Asian-Pacific-American region[J]. Scientific Reports, 6: 36356.

DOI

[23]
LEUNG M Y T, CHEUNG H H N, ZHOU WEN, 2017. Meridional displacement of the East Asian trough and its response to the ENSO forcing[J]. Climate Dynamics, 48(1-2): 335-352.

DOI

[24]
LEUNG M Y T, ZHOU WEN, CHEUNG K Y, et al, 2019. Enhancement of lower tropospheric winter synoptic temperature variations in Southwest China and the northern Indochina Peninsula after 2010[J]. Climate Dynamics, 53(3-4): 2281-2294.

DOI

[25]
MAVUME A F, RYDBERG L, ROUAULT M, et al, 2009. Climatology and landfall of tropical cyclones in the south-west Indian Ocean[J]. Western Indian Ocean Journal of Marine Science, 8(1): 15-36.

[26]
MCADIE C J, LAWRENCE M B, 2000. Improvements in tropical cyclone track forecasting in the Atlantic Basin, 1970-98[J]. Bulletin of the American Meteorological Society, 81(5): 989-998.

DOI

[27]
NG C H J, VECCHI G A, 2020. Large-scale environmental controls on the seasonal statistics of rapidly intensifying North Atlantic tropical cyclones[J]. Climate Dynamics, 54(9-10): 3907-3925.

DOI

[28]
RAPPAPORT E N, FRANKLIN J L, AVILA L A, et al, 2009. Advances and challenges at the National Hurricane Center[J]. Weather and Forecasting, 24(2): 395-419.

DOI

[29]
SCHADE L R, 2000. Tropical cyclone intensity and sea surface temperature[J]. Journal of the Atmospheric Sciences, 57(18): 3122-3130.

DOI

[30]
SUN DONGLIAN, KAFATOS M, CERVONE G, et al, 2007. Satellite microwave detected SST anomalies and hurricane intensification[J]. Natural Hazards, 43(2): 273-284.

DOI

[31]
WADA A, USUI N, 2007. Importance of tropical cyclone heat potential for tropical cyclone intensity and intensification in the Western North Pacific[J]. Journal of Oceanography, 63(3): 427-447.

DOI

[32]
WALSH K J E, RYAN B F, 2000. Tropical cyclone intensity increase near Australia as a result of climate change[J]. Journal of Climate, 13(16): 3029-3036.

DOI

[33]
WANG BIN, ZHOU XIAQIONG, 2008. Climate variation and prediction of rapid intensification in tropical cyclones in the western North Pacific[J]. Meteorology and Atmospheric Physics, 99(1-2): 1-16.

DOI

[34]
WANG CHUNZAI, WANG XIDONG, WEISBERG R H, et al, 2017. Variability of tropical cyclone rapid intensification in the North Atlantic and its relationship with climate variations[J]. Climate Dynamics, 49(11-12): 3627-3645.

DOI

[35]
WANG XIDONG, WANG CHUNZAI, ZHANG LIPING, et al, 2015. Multidecadal variability of tropical cyclone rapid intensification in the western North Pacific[J]. Journal of Climate, 28(9): 3806-3820.

DOI

[36]
WU LONGTAO, SU HUI, FOVELL R G, et al, 2012. Relationship of environmental relative humidity with North Atlantic tropical cyclone intensity and intensification rate[J]. Geophysical Research Letters, 39(20): L20809.

[37]
XIAO FUAN, WANG DONGXIAO, LEUNG M Y T, 2020. Early and extreme warming in the South China Sea during 2015/2016: role of an unusual Indian Ocean dipole event[J]. Geophysical Research Letters, 47(17): e2020GL089936.

[38]
XIE SHANGPING, ANNAMALAI H, SCHOTT F A, et al, 2002. Structure and mechanisms of South Indian Ocean climate variability[J]. Journal of Climate, 15(8): 864-878.

DOI

[39]
XIE SHANGPING, HU KAIMING, HAFNER J, et al, 2009. Indian ocean capacitor effect on indo-western pacific climate during the summer following El Niño[J]. Journal of Climate, 22(3): 730-747.

DOI

[40]
YU LISAN, RIENECKER M M, 1999. Mechanisms for the Indian Ocean warming during the 1997-98 El Niño[J]. Geophysical Research Letters, 26(6): 735-738.

DOI

[41]
ZHANG QIANG, WU LIGUANG, LIU QIUFENG, 2009. Tropical cyclone damages in China 1983-2006[J]. Bulletin of the American Meteorological Society, 90(4): 489-496.

DOI

[42]
ZHANG WENJUN, WANG YALAN, JIN FEIFEI, et al, 2015. Impact of different El Niño types on the El Niño/IOD relationship[J]. Geophysical Research Letters, 42(20): 8570-8576.

DOI

[43]
ZHANG YUE, ZHOU WEN, LEUNG M Y T, 2019. Phase relationship between summer and winter monsoons over the South China Sea: Indian Ocean and ENSO forcing[J]. Climate Dynamics, 52(9-10): 5229-5248.

DOI

[44]
ZHAO HAIKUN, DUAN XINGYI, RAGA G B, et al, 2018. Changes in characteristics of rapidly intensifying western North Pacific tropical cyclones related to climate regime shifts[J]. Journal of Climate, 31(19): 8163-8179.

DOI

Outlines

/