2
years 1979 to present, the opportunity to extend the F16
analyses to additional years became possible. The goal of
this study is to determine the representativeness of F16’s
findings compared to a longer period (2008-2017) and
a different, and presumably improved, reanalysis dataset.
Utilizing this longer analysis period, we investigate vari-
ability of vortex formation and its relationship to various
large-scale signals such as the MJO, IOD, and ENSO sig-
nals. Finally, we consider some long-term IO TC statistics
to examine the relationship of TCs to wind regimes and
large-scale indices.
2. Data and Methodology
a. Datasets
ECMWF Reanalysis 5th Generation (or ERA5), which
replaces the ERA-Interim reanalysis, is based on 4D-Var
data assimilation using Cycle 41r2 of the Integrated Fore-
casting System (IFS), which was operational at ECMWF in
2016. A detailed description of the ERA5 configuration,
which benefits from a decade of developments in model
physics, core dynamics, and data assimilation relative to
ERA-Interim, can be found in Hersbach et al. (2020). For
this study we used both vorticity and zonal winds from
ERA5.
The Madden-Julian Oscillation (MJO) is the major fluc-
tuation in tropical weather on weekly to monthly timescales
(Madden and Julian 1971). The MJO can be character-
ized as an eastward moving ’pulse’ of cloud and rainfall
near the equator that typically recurs every 30 to 60 days.
The location and strength of the MJO are given by the
Real-time Multivariate MJO (RMM) index developed by
Wheeler and Hendon (2004). This index is based on a pair
of empirical orthogonal functions (EOFs) of the combined
fields of near-equatorially averaged 850-hPa zonal wind,
200-hPa zonal wind, and satellite-observed outgoing long-
wave radiation (OLR) data. Daily values of the RMM
index provide the amplitude and phase of the MJO. For
this study, classification of the vortices by MJO phase was
only done when the MJO signal had a significantly robust
amplitude defined here as when the RMM amplitude was
> 1.
To represent the variability of the El Niño/Southern Os-
cillation (ENSO) we use the Multivariate ENSO Index
Version 2 (MEI.v2) which combines both oceanic and at-
mospheric variables into a single assessment of the state of
ENSO (Wolter and Timlin 2011; Zhang et al. 2019). Pos-
itive (negative) values of this index imply warm, El Niño
(cool, La Niña) conditions across the east-central equato-
rial Pacific. The MEI.v2 is available on a monthly basis.
Variability in the IO region is often characterized by
changes in the east-west sea surface gradient across the IO
basin referred as the Indian Ocean Dipole (IOD) mode (Saji
et al. 1999). The IOD is commonly measured by the Dipole
Mode Index (DMI), that is, the difference between sea
surface temperature (SST) anomalies between the western
(50◦-70◦E) and eastern (90◦-110◦E) tropical (10◦S-10◦N)
IO. A positive (negative) IOD period is characterized by
cooler (warmer) than average water in the tropical eastern
Indian Ocean and warmer (cooler) than average water in
the tropical western Indian Ocean.
Information on IO tropical cyclone tracks and inten-
sity was obtained from the Joint Typhoon Warning Center
(JTWC) best-track archive which goes back to 1945 (Chu
et al. 2002). For this study TC data were used from 1980
to 2017, although Chu et al. (2002) state that the years
after 1984 have the best data quality. The best-track data
contains the storm center locations and intensities (i.e., the
maximum 1-minute mean sustained 10-meter wind speed)
at six-hour intervals. Storms are only considered here that
(1) reach TC intensity (i.e., maximum sustained winds of
35 knots or 39 mph) and (2) formed over the IO (i.e., ex-
tending from Africa to 105◦E). Storms with missing wind
speeds, which are most prevalent in the earlier years of the
archive, were not considered.
b. Vortex tracking
As in F16, identification and tracking of low-level lee
vortices was carried out based on the relative vorticity field
using the objective feature tracking code of Hodges (1995,
1999). To facilitate the use of this software, relative vor-
ticity at 6 h and 0.25◦horizontal resolution from ERA5
analyses were vertically averaged using data at 850, 875,
900, 925 hPa from 50◦-110◦E, 20◦N-20◦S. The F16 anal-
yses had only three vertical levels in the 850 to 925 hPa
layer. After smoothing the vorticity field to retain spatial
scales greater than 450 km, cyclonic vortex features were
tracked if they maintained an amplitude greater than 1.0 ×
10−5s−1for longer than 2 days.
To focus on wake vortices with the potential to develop
into TCs, we restrict our analyses to cases where the vortex
moved westward over the Indian Ocean. As in F16 these
shed vortices are defined as those with 1) a final location
over the Indian Ocean that is > 500 km from Sumatra, and
either 2) a final minus initial displacement from Sumatra
> 250 km, or 3) their average speed away from Sumatra
is > 0.5 ms−1. The first condition ensures that the shed
vortex at the end of its track is some critical distance from
Sumatra while conditions 2 or 3 guarantee that the vortex
is moving away from its generating landmass.
Application of the tracking code to ten years of ERA5
data yields information on both the tracks and locations
of cyclonic lee vortices. Because of seasonal changes in
the flow regime, genesis locations are shown in Fig. 2 for
the boreal cold season (November to April) and the boreal
warm season (May to October). Overall, these analyses
based on 10 years of statistics are remarkably similar to
those shown in F16. For example, during the boreal win-
ter season (Fig. 2a), northeasterly flow across the Malay