Acceleration of Tropical Cyclones As a Proxy For Extratropical Interactions: Synoptic-Scale Patterns and Long-Term Trends
AAcceleration of Tropical Cyclones As a Proxy For ExtratropicalInteractions: Synoptic-Scale Patterns and Long-Term Trends
Anantha Aiyyer and Terrell Wade Department of Marine, Earth and Atmospheric Sciences, North Carolina State University
Correspondence:
A. Aiyyer ([email protected])
Abstract.
It is well known that rapid changes in tropical cyclone motion occur during interaction with extratropical waves.While the translation speed has received much attention in the published literature, acceleration has not. Using a large datasample of Atlantic tropical cyclones, we formally examine the composite synoptic-scale patterns associated with tangential and curvature components of their acceleration. During periods of rapid tangential acceleration, the composite tropical cy-clone moves poleward between an upstream trough and downstream ridge of a developing extratropical wavepacket. The twosystems subsequently merge in a manner that is consistent with extratropical transition. During rapid curvature acceleration,a prominent downstream ridge promotes recurvature of the tropical cyclone. In contrast, during rapid tangential or curva-ture deceleration, a ridge is located directly poleward of the tropical cyclone. Locally, this arrangement takes the form of acyclone-anticyclone vortex pair somewhat akin to a dipole block. On average, the tangential acceleration peaks 18 hours priorto extratropical transition while the curvature acceleration peaks at recurvature. These findings confirm that rapid accelerationof tropical cyclones is mediated by interaction with extratropical baroclinic waves. Furthermore, The tails of the distribution ofacceleration and translation speed show a robust reduction over the past 5 decades. We speculate that these trends may reflectthe poleward shift and weakening of extratropical Rossby waves.
The track and movement of tropical cyclones are known to be governed by the background environment (e.g., Hodanishand Gray, 1993). It was recognized early on that the translation speed of a tropical cyclone can be approximated by thesurrounding wind field (Emanuel, 2018). A tropical cyclone, however, is not an isolated vortex that is passively carried by thecurrent. The background environment is comprised of synoptic and large-scale circulation features, with attendant gradientsof potential vorticity, moisture, and deformation. The tropical cyclone actively responds to these external stimuli. The largescale environment is also impacted by the interaction. For example, the generation of β - gyres – that influence tropical cyclonemotion – is a response to the background potential vorticity gradient (e.g., Chan and Williams, 1987; Wu and Emanuel, 1993).Similarly, vertical wind shear can limit the intensification of tropical cyclones by importing low entropy air (e.g., Tang andEmanuel, 2010). This, in turn, can impact subsequent tropical cyclone motion. On the other hand, a tropical cyclone caninfluence the large-scale flow by exciting waves in the extratropical stormtrack, leading to rapid downstream development and a r X i v : . [ phy s i c s . g e o - ph ] J a n earrangement of the flow (e.g., Jones et al., 2003). It can be contended that a tropical cyclone is always interacting with itsenvironment, and the interaction is partly manifested in its motion.Track and translation speed are two aspects of tropical cyclone motion that are particularly important for operational fore-casts. The track garners much attention for obvious reasons – it informs potential locations that will be affected by a storm.The translation speed impacts intensity change, storm surge and local precipitation amount. There is a large body of publishedliterature and review articles dealing with various research and operational problems related to tropical cyclone track and speed(e.g., Chan, 2005; Emanuel, 2018). When tropical cyclones move poleward, they often encroach upon the extratropical storm-track. This leads to their interaction with baroclinic eddies of the stormtrack. The outcome of the interaction is varied. Sometropical cyclones weaken and dissipate while others strengthen and retain their tropical nature for an additional period. Inthe North Atlantic, around 50% of tropical cyclones eventually experience extratropical transition - a complex process duringwhich the warm-core tropical cyclone evolves into a extratropical cyclone and becomes part of an extratropical stormtrack (e.g.,Hart and Evans, 2001). Several recent reviews have extensively documented the research history and dynamics of extratropicaltransition (e.g., Evans et al., 2017; Keller et al., 2019). Forecasters have long known that tropical cyclones accelerate forwardduring extratropical transition, but relatively less attention has been devoted to the details of this acceleration in the researchliterature.This paper has two main themes. The first concerns the synoptic-scale flow associated with rapid acceleration and decelera-tion of tropical cyclones. This is addressed via composites of global reanalysis fields and observed tropical cyclone tracks from1980–2016. We consider both tangential and curvature accelerations. To our knowledge, such formal delineation of extratrop-ical interaction based on categories of tropical cyclone acceleration has not been presented in the literature. Of note, however,is the recent work of Riboldi et al. (2019) that is relevant to this paper. That study examined the interaction of acceleratingand decelerating upper-level troughs and recurving western North Pacific typhoons. Their key findings are: (a) In the major-ity of cases, a recurving tropical cyclone is associated with a decelerating upper-level trough that remains upstream; (b) Theupper-level trough appears to phase lock with the tropical cyclone; and (c) Recurvatures featuring such trough deceleration arefrequently associated with downstream atmospheric blocking. As we shown in subsequent sections, many of these results canbe recovered independently when we approach the problem from the perspective of tropical cyclone motion. As such, part ofour work complements the findings of Riboldi et al. (2019).The second theme concerns the identification of long-term trends in tropical cyclone motion. While we focus on acceleration,we begin with translation speed to place our results within the context of related recent work. Kossin (2018) reported that thetranslation speed of tropical cyclones over the period 1949–2016 has reduced by about 10% over the globe, and about 6%over the North Atlantic. Without directly attributing it, Kossin (2018) noted that the trend is consistent with the observedslow-down of the atmospheric circulation due to anthropogenic climate change. In addition to the circulation changes, thegeneral warming of the atmosphere is associated with an increase in water vapor content per the Clausius-Clapeyron scaling(CCS). This translates to increase in precipitation rates that can locally exceed the CCS (e.g., Nie et al., 2018). Kossin (2018)made the point that the slowing of tropical cyclones may compound the problem of heavy precipitation in a warmer climate.However, Moon et al. (2019) and Lanzante (2019) argued that the historical record of tropical cyclone track data, particularly rior the advent of the satellite era around the mid-1960s, is likely incomplete. They also showed that annual-mean tropicalcyclone translation speed exhibits step-like changes and questioned the existence of a true monotonic trend. They attributedthese discrete changes to natural factors (e.g., regional climate variability) as well as omissions and errors due to lack of directobservation prior to the availability of extensive remote sensing tools.We revisit the issue of trends in tropical cyclone motion, but restrict our attention to the Atlantic basin and the years since1966. This year is considered to be the beginning of the satellite era , at least for the Atlantic (Landsea, 2007). In contrast,global coverage by geostationary satellites began much later, in 1981 (e.g., Schreck et al., 2014). Lanzante (2019) found achange point in 1965 for the timeseries of annual average tropical cyclone speed for the North Atlantic basin. Lanzante (2019)also reported that accounting for the change point dramatically reduces the trend from the value reported by Kossin (2018). Inlight of the issues raised about the reliability of tropical cyclone track data prior to the satellite era , we use 1966 as our startingpoint and seek to ascertain whether robust trends exist in the observed record of tropical cyclone acceleration. We use the IBTrACS v4 (Knapp et al., 2010) for tropical cyclone locations. The information in this database typically spansthe tropical depression stage to extratropical cyclone and/or dissipation. Kossin (2018), used all locations in the IBTrACs aslong as a storm lasted more than 3 days and did not apply a threshold for maximum sustained winds. We follow the samemethod with one major difference. We only consider those instances of the track information wherein a given storm was stillclassified as tropical. We recognize the fact that the nature of a storm is fundamentally altered after it loses its surface enthalpyflux-driven, warm-core structure. After extratropical transition, evolution is governed by stormtrack dynamics of barocliniceddies. We wish to avoid conflating the motions of two dynamically distinct weather phenomena. For the same reason, wealso omit subtropical storms which are catalogued in the track database. We argue that, by restricting our track data to onlyinstances that were deemed to be tropical in nature, we can paint a more appropriate picture of the composite environment andtrends associated with the rapid acceleration of tropical cyclones.To ascertain whether a tropical storm underwent extratropical transition, we use the nature designation in the IBTrACSdatabase that relies on the judgement of the forecasters from one or more agencies responsible for an ocean basin. In IBTrACS,the nature flag is set to "ET" after the transition is complete. Admittedly, there will be some subjectivity in this designation. Analternative would be to employ a metric such as the cyclone phase space (Hart, 2003), usually calculated using meteorologicalfields from global numerical weather prediction models (e.g., reanalysis). As we wish to remain independent of modeledproducts to characterize the storms, we rely on the forecaster designated storm nature. Furthermore, Bieli et al. (2019) foundthat the phase space calculation was sensitive to the choice of the reanalysis model, which would add another source ofuncertainty. Occasionally, the nature of the storm is not recorded (NR) or, if there is an inconsistency among the agencies, itis designated as mixed (MX). In the north Atlantic, only a small fraction of track data ( ≈ . ) in the IBTrACs is designatedas NR or MX. This fraction is higher in other basins (e.g., ≈ in the western North Pacific). This is another reason why e restrict our attention to the Atlantic basin. Henceforth, we will collectively refer to the designations NR, MX and ET as non-tropical , and unless explicitly stated, omit the associated track data in our calculations.For the trends and basic statistics, we focus on the years 1966–2019, within the ongoing satellite era for the North Atlanticbasin (e.g. Landsea, 2007; Vecchi and Knutson, 2008). Compared to years prior, tropical cyclone data is deemed more reli-able once satellite-based observations became available. For the composites, we use the European Center for Medium RangeWeather Forecasting (ECMWF) ERA-Interim (ERAi) reanalysis (Dee et al., 2011) for the period 1981–2016. In subsequentsections, we focus on the region 20–50 o N, wherein tropical cyclones are more likely to interact with extratropical baroclinicwaves.
The acceleration of a hypothetical fluid element moving with the center of a tropical cyclone can be written as: a = dVdt ˆs + V R ˆn (1)where V is the forward speed and R is the radius of curvature of the track at a given location. Here, ˆs and ˆn are orthogonalunit vectors in the so-called natural coordinate system. The former is directed along the tropical cyclone motion. The latter isdirected along the radius of curvature of the track. The first term in eq. 1 is the tangential acceleration and the second is thecurvature or normal acceleration. The speed V j at any track location (given by index j) is calculated as: V j = ( D j,j − + D j,j +1 ) δt (2)where D refers to the distance between the two consecutive points, indexed as shown above, along the track. Since we used3-hourly reports from the IBTrACS, δt = 6 h. The tangential acceleration was calculated from the speed using centered differ-ences.To calculate the curvature acceleration, it is necessary to first determine the radius of curvature, R. A standard approach tocalculating R, given a set of discrete points along a curve – in our case, a tropical cyclone track – is to fit a circle through threeconsecutive points. For a curved line on a sphere, it can be shown that: R = R e sin − (cid:32)(cid:115) d d d ( d + d + d ) − d + d + d ) (cid:33) (3.1)where R is the radius of curvature, R e is the radius of the Earth, and the d terms are expressed as follows: d = 1 − (cos T cos T cos ( N − N ) + sin T sin T ) (3.2) d = 1 − (cos T cos T cos ( N − N ) + sin T sin T ) (3.3) d = 1 − (cos T cos T cos ( N − N ) + sin T sin T ) (3.4) here, T , T , and T are the latitudes of the 3 points while N , N , and N are the longitudes. The center of the circle is givenby the coordinates: tan T = ± cos T cos T sin ( N − N ) + cos T cos T sin ( N − N ) + cos T cos T sin ( N − N ) (cid:112) α + β (4.1) tan N = − αβ (4.2)Where T and N are the latitude and longitude, respectively of the circle’s center, and α and β are obtained using: α = cos T (sin T − sin T ) cos N + cos T (sin T − sin T ) cos N + cos T (sin T − sin T ) cos N (4.3) β = cos T (sin T − sin T ) sin N + cos T (sin T − sin T ) sin N + cos T (sin T − sin T ) sin N (4.4) Figure 1.
Illustration of the circle-fit and radius of curvature calculations at five selected locations along the track of hurricanes Katrina(2005) and Irma (2017).
While the expressions above were derived independently for this work, we make no claim of originality since they are basedon elementary principles of geometry. Figure 1 shows some examples of the radius of curvature calculation for hurricanesKatrina (2005) and Irma (2017).
Over the years 1966–2019, 689 storms in the North Atlantic meet the 3-day threshold. Figure 2 shows some basic statisticsfor speed and accelerations as a function of latitude for Atlantic storms. Tropical cyclone locations are sorted in ◦ -widelatitude bins and all instances classified as non-tropical were excluded. Table 1 provides the same data. The average speed ofall North Atlantic tropical cyclones (including over-land) is about 21 km/hr. As expected, tropical cyclone speed clearly varies igure 2. Distribution of (a) Speed(km hr − ) ; (b) Tangential acceleration (km hr − day − ) and (c) Curvature acceleration (km hr − day − )of Atlantic TCs as a function of latitude. Storm instances classified as ET or NR were excluded. Data from 1966–2019 was binned within10 o -wide overlapping latitude bins. Statistics shown are: median (horizontal line within the box), mean (dot), and 10th, 25th, 75th and 90thpercentiles. with latitude. It is lower in the subtropics as compared to other latitudes. It increases sharply in the vicinity of o N. Thetangential acceleration, as defined in eq. 1, can be positive or negative. The mean and median tangential acceleration remainsnear-zero equatorward of ◦ N. This suggests that tropical cyclones in this region tend to translate steadily, or are equallylikely to accelerate and decelerate. The tangential acceleration is substantially positive poleward of ◦ N. Tropical cyclonesin these latitudes are subject to rapid acceleration. For example, the mean tangential acceleration in the ◦ − ◦ N latitudeband is about 5.0 km/hr day − . The curvature acceleration, by our definition, takes only positive values and steadily increaseswith latitude. The distributions of tropical cyclone speed and tangential acceleration are relatively symmetric about the medianvalue as compared to the curvature acceleration. We now examine the flow pattern associated with rapid tangential and curvature acceleration of tropical cyclones. For this,storm-relative ensemble average fields are constructed using ERA-interim reanalysis over the period 1980–2016. We use themethod outlined as follows. – All tropical cyclone track locations are binned into 10 ◦ wide latitude strips (e.g., 10-20 ◦ N, 20-30 ◦ N, 30-40 ◦ N).Instances where storms are classified as non-tropical (c.f. section 2) are excluded. For brevity, only results for the 30-40 ◦ N bin are discussed here. A total of 3515 track points are identified for this latitude bin. Note that a particular tropicalcyclone could appear more than once in a latitude bin at different times. – The tangential accelerations in each bin are separated into two categories: rapid acceleration and rapid deceleration. Therapid acceleration composite is calculated using all instances of acceleration exceeding a threshold: i.e., a ≥ a τ , where able 1. Speed (km hr − ), tangential and curvature acceleration (km hr − day − ) of Atlantic TCs in the IBTRaCs database as a function oflatitude. Storm instances classified as ET or NR were excluded. N refers to number of 3-hourly track positions in each latitude-bin over theperiod 1966–2019. Speed Tang. Accel Curv. AccelLatitude N Mean Median Std Dev. Mean Median Std Dev. Mean Median Std Dev.Full Basin 38822 20.85 18.87 12.2 2.27 0.68 19.5 16.42 10.88 19.05–15 5870 23.53 23.4 9.4 0.18 0.0 14.7 11.65 8.2 13.110–20 13087 21.10 20.6 9.3 -0.13 -0.0 15.3 12.29 8.5 13.315–25 13656 18.66 18.1 8.6 0.26 -0.0 16.0 13.90 9.4 15.620–30 13074 17.21 16.4 8.6 1.40 0.4 17.1 16.48 11.2 19.025–35 13905 17.74 16.2 10.2 2.74 1.3 19.7 18.04 12.2 20.630–40 10815 21.37 18.6 13.5 5.33 3.0 22.8 19.34 13.4 21.035–45 5272 29.41 26.7 18.0 8.27 4.8 27.4 22.41 15.8 22.940–50 1607 42.37 41.0 20.8 8.75 6.4 34.9 28.45 19.9 29.045–55 329 55.69 53.9 19.4 6.81 6.3 37.8 36.50 25.9 38.9 τ refers to a specified quantile of the acceleration distribution within the latitude bin. Similarly, the rapid tangentialdeceleration composite is based on all instances where a ≤ a τ . We tried a variety of thresholds for τ (e.g., . − . for rapid tangential acceleration, and . − . for rapid tangential deceleration). Our conclusions for the compositesare not sensitive to the exact choice of the threshold values as long as they are sufficiently far from the median. – We also use the same method to create categories of curvature accelerations. For this, since we only have positive values,we interpret these two quantile-based categories as rapid acceleration and near-zero acceleration respectively. – For each category, we compute an ensemble average composite field of the geopotential field at selected isobaric levels.The composites are calculated after shifting the grids such that the centers of all storms were coincident. The centroidposition of storms was used for the composite storm center, and the corresponding time was denoted as Day-0. Lag-composites were created by shifting the dates backward and forward relative to Day-0. – For anomaly fields, we subtract the long-term synoptic climatology from the total field. The climatology is calculatedfor each day of the year by averaging data for that day over the years 1980-2015. This is followed by a 7-day runningmean smoother. To account for diurnal fluctuations, the daily climatology is calculated for each available synoptic hourin the ERA-interim data (00, 06, 12, and 18 UTC). – For brevity, we only show the results for τ = 0 . for rapid acceleration and τ = 0 . for rapid deceleration. For reference,within 30–40 o N latitude range, τ = 0 . corresponds to a = 32 km/hr day − and the τ = 0 . corresponds a = − km/hrday − . For curvature acceleration, they are, respectively, 48 and 32 km/hr day − . The sample size of each composite was
52 ( ≈ of the total number of track points in this latitude range. These correspond to 196 and 168 unique stormsfor rapid acceleration and rapid deceleration respectively. – Statistical significance of anomaly fields shown in this section are evaluated by comparing them against 1000 compositescreated by randomly drawing 352 dates for each composite from the period July–October, 1980–2015. A two-tailedsignificance is evaluated at the 95% confidence level with the null hypothesis being that the anomalies could haveresulted from a random draw.
Figure 3.
Storm-relative composite average geopotential heights (thick orange lines) and anomalies (color shaded) for all TCs located in thelatitude bin 30-40 o N over the Atlantic. The composite fields are shown for three levels – 300 hPa, 500 hPa and 850 hPa. In each panel, thecomposite 1000 hPa anomalous geopotential is shown using thin black contours. All anomalies are defined relative to a long-term synopticclimatology. The contour intervals are: 12 dam, 6 dam and 3 dam for the three levels respectively. The shading interval in dam for the 300hPa anomaly fields is shown in the label-bar. It is half of the value for the other two levels. The left column is for rapid tangential accelerationand the right column is for rapid tangential deceleration. igure 3 shows storm-centered composite geopotential heights (thick orange lines) and their anomalies (color shaded) forall Atlantic tropical cyclones located within 30-40 ◦ N. Two categories of tangential acceleration are shown: rapid acceleration(left column) and rapid deceleration (right column). The fields are shown at three levels – 300 hPa, 500 hPa and 850 hPa. Ineach panel, the anomalous 1000 hPa geopotential is shown using thinner black contours. It highlights the composite tropicalcyclone and the surface development within the extratropical stormtrack.The ensemble average for rapid tangential acceleration (left column of Fig. 3) shows the composite tropical cyclone inter-acting with a well defined extratropical wavepacket. The tropical cyclone is straddled by an upstream trough and a downstreamridge. At 500 hPa, the geopotential anomalies of the tropical cyclone and the upstream trough are close and, consequently,appear to be connected. This yields a negative tilt in the horizontal and indicates the onset of cyclonic wrap-up of the trougharound the tropical cyclone. The 1000-hPa geopotential anomaly field is dominated by the composite tropical cyclone. It alsoshows the relatively weaker near-surface cyclones and anticyclones of the extratropical stormtrack. The entire wavepacketshows upshear tilt of geopotential anomalies with height, indicating baroclinic growth. This arrangement of the tropical cy-clone and the extratropical wavepacket is consistent with the synoptic-scale flow that is typically associated with extratropicaltransition (e.g., Bosart and Lackmann, 1995; Klein et al., 2000; McTaggart-Cowan et al., 2003; Riemer et al., 2008; Riemerand Jones, 2014; Keller et al., 2019). At this point, all storms in the ensemble were still classified as tropical . Thus, we interpretthis composite as pre-extratropical transition completion state.The ensemble average for rapid tangential deceleration cases (right column of Fig. 3) shows an entirely different synoptic-scale pattern. The extratropical wavepacket is substantially poleward, with a ridge immediately north of the composite tropicalcyclone. The geopotential anomalies of the extratropical wavepacket and the composite tropical cyclone appear to be distinctat all three levels, with no evidence of merger. The prominent synoptic structure is the cyclone-anticyclone dipole formed bythe tropical cyclone and the extratropical ridge.
To get a sense of the temporal evolution of the entire system, we show lag composites for Day-2 to Day+2 in Fig. 4. As inthe previous figure, the two categories of acceleration are arranged in the respective columns. The rows now show 500-hPageopotential height (thick contours) and anomalies (color shaded). In each panel, the corresponding 1000-hPa geopotentialheight anomalies are shown by thin black contours.The ensemble average for rapid tangential acceleration (left column of Fig. 4) shows a tropical cyclone moving rapidlytowards an extratropical wavepacket. At day-2, the tropical cyclone circulation is relatively symmetric as depicted by thecontours of 1000 hPa geopotential anomalies. The downstream extratropical ridge is prominent, but the upstream trough ismuch weaker at this time. On Day-1, the entire extratropical wavepacket has amplified and the 500 hPa geopotential anomaliesof the tropical cyclone and a developing upstream trough have merged. This process continues through Day 0. By Day+1,the composite storm has moved further poleward and eastward and is now located between the upper-level upstream troughand downstream ridge in a position that is optimal for further baroclinic development. The 1000 hPa geopotential field is nowasymmetric with a characteristic signal of a cold front. igure 4. Storm-relative average 500-hPa geopotential heights (thick orange lines) and anomalies (color shaded) for all TCs located in thelatitude bin 30-40 o N over the Atlantic. The fields are shown for lags Day-2 to Day+2. In each panel, the composite 1000 hPa anomalousgeopotential is shown using thin black contours. All anomalies are defined relative to a long-term synoptic climatology. The contour intervalis 6 dam and shading interval in dam is shown in the label-bar. The plus symbol shows the location of the composite TC at Day 0 and thehurricane symbol shows the approximate location at each lags. The left column is for rapid tangential acceleration and the right column isfor rapid tangential deceleration igure 5. Track of the composite tropical cyclone (blue) and the downstream 500-hPa extratropical ridge (orange) from Day-2 to Day +2 for(a) rapid tangential acceleration, and (b) rapid tangential deceleration. The composites are based on all TC tracks locations within 30–40 o N.Day 0 is the reference day for the composites in Fig. 4
The picture that is evident from these 500-hPa composite fields is that, over the course of the 4 days, downstream ridge-troughcouplet amplifies while simultaneously propagating eastward. The upstream trough cyclonically wraps around the tropicalcyclone and the two have merged by Day 1. The geopotential gradient poleward of the storm is also enhanced, indicating astrengthening jet streak. These features are consistent with the process of extratropical transition (e.g., Keller et al., 2019). Thepoleward moving tropical cyclone may either interact with an existing wavepacket or perturb the extratropical flow and excitea Rossby wavepacket that disperses energy downstream (e.g., Riemer and Jones, 2014). The outflow of the tropical cyclone isa source of low potential vorticity (PV) air that further reinforces the downstream ridge (e.g., Riemer et al., 2008).To further illustrate the interaction, the tracks of the tropical cyclone and the 500-hPa ridge of the extratropical wavepacketare presented in Fig. 5a. It can be clearly seen that the tropical cyclone merges with the extratropical stormtrack. Furthermore,the 500 hPa ridge has rapidly moved downstream during the 4 day period, indicating a very progressive pattern. The eastwardphase speed of the extratropical wavepacket, as inferred from the track of the ridge, is ≈ ms − . The tropical cyclone speed,averaged over the same 4-day period, is ≈ ms − . The close correspondence between the two and the merger of the tracksfurther supports the notion that the synoptic-scale evolution during rapid acceleration cases is consistent with the canonicalpattern associated with extratropical transition.On the other hand, during rapid deceleration (right column Fig. 4), the composite tropical cyclone remains equatorward ofthe extratropical wavepacket and maintains a nearly symmetric structure throughout the period. The arrangement of the tropicalcyclone and the extratropical ridge is akin to a vortex dipole. The extratropical wavepacket is not as progressive as in the rapidacceleration case. This is seen clearly from the tracks of the tropical cyclone and the ridge (Fig. 5b). The phase speed of theextratropical wavepacket is ≈ ms − while the tropical cyclone speed is ≈ . ms − . The phasing of the tropical cyclone andthe extratropical wavepacket has led to the formation of a cyclone-anticyclone vortex dipole. We return to this point and relateit to similar findings in Riboldi et al. (2019) in a later section. igure 6. As in Fig. 3, but for rapid and near-zero curvature acceleration
As in the previous section, Figure 6 shows storm-centered ensemble averages, but this time for the two categories of curvatureacceleration. The composite for rapid acceleration (left column) shows a tropical cyclone that is primarily interacting withan extratropical ridge that is poleward and downstream of it. The upstream trough in the extratropics is weaker and fartherwestward as compared to the rapid tangential acceleration composite (Fig. 6a). Furthermore, instead of the upstream troughwrapping cyclonically, in this case, we see the downstream ridge wrapping anticyclonically around the tropical cyclone. Thisis similar to the composite 500-hPa fields that were based on recurving tropical cyclones as shown in Fig. 5 of Aiyyer (2015).Thus, in an ensemble average sense, rapid curvature acceleration appears to mark the point of the recurvature of tropicalcyclones.The composite for near-zero curvature acceleration (right column) is quite similar to the composite for rapid tangentialdeceleration (Fig. 6d–f). The extratropical wavepacket is poleward and the tropical cyclone-ridge system appears as a vortexdipole. igure 7. As in Fig. 4, except for rapid (left column) and near-zero (right column) curvature acceleration.
The temporal evolution of the entire system for the two categories of curvature acceleration is shown in Fig. 7. For rapid cur-vature acceleration, we see a tropical cyclone that is moving poleward towards an extratropical ridge. During the subsequentdays, the ridge moves eastward initially and begins to wrap around the tropical cyclone. This arrangement promotes the re- urving of the tropical cyclone. By Day+2, the anticyclonic wrapping and thinning has resulted in a significantly weaker ridgeas compared to a few days prior. For the near-zero curvature cases ( Fig. 7f–g), the initial movement of the tropical cycloneis also directly poleward towards the extratropical ridge. However, in this case, the ridge remains poleward of the tropicalcyclone. There is also significantly less anticyclonic wrapping of ridge. The tropical cyclone-ridge system takes the form of acyclonic-anticyclonic vortex pair similar to the rapid tangential deceleration composite.The tracks in Fig. 8 clearly show how the tropical and extratropical systems propagate. For rapid curvature acceleration (Fig.8a), the tropical cyclone track shows a recurving tropical cyclone. The track of the ridge confirms the initial eastward motion,followed by a poleward shift after parts of it wrap around the tropical cyclone as noted from Fig. 7a-e. By Day+2, we do notobserve a merger of the tracks that happens in the case of rapid tangential acceleration (Fig. 5a).The tracks for near-zero curvature acceleration (Fig. 8a) are somewhat similar to the rapid tangential deceleration (Fig. 5b).The key point here is that, although the tropical cyclone moves poleward, tropical cyclone-ridge system acts like a vortex dipoleand is nearly stationary in the zonal direction. This arrangement of the tropical cyclone and the extratropical wavepacket issimilar to the composite fields in Fig. 10 of Riboldi et al. (2019), where they show upper-level potential vorticity (PV), 850-hPapotential temperature and sea-level pressure. The difference is that their composite was conditioned on the acceleration of theupstream trough for recurving western Pacific typhoons. It is, however, not surprising that we can recover a similar patternwhen we condition our composites on the basis of tropical cyclone acceleration. Since the the extratropical wavepacket and thetropical cyclone are actively interacting, they influence each other’s motion. Riboldi et al. (2019) referred to this as a phase-locking of the upstream trough and the tropical cyclone while we have viewed this as a phase-lock between the ridge and thetropical cyclone. The two are not mutually exclusive since the trough and the ridge are part of the same wavepacket. Figure 8.
Track of the composite tropical cyclone (blue) and the downstream 500-hPa extratropical ridge (orange) from Day-2 to Day +2for (a) rapid curvature acceleration, and (b) near-zero curvature acceleration. The composites are based on all TC tracks locations within30–40 o N. Day 0 is the reference day for the composites in Fig. 7
14 Extratropical Transition
In the previous section, we showed that the composite synoptic-scale flow associated with rapid tangential acceleration re-sembles a pattern that is favorable for extratropical transition. However, this does not imply that all storms in the compositeunderwent extratropical transition. Some tropical cyclones may begin the process of extratropical transition but dissipate beforeits completion (e.g., Kofron et al., 2010). We now consider tropical cyclone motion from a different perspective by consideringonly those storms that completed the transformation from being tropical to extratropical. Hart et al. (2006) found that the timetaken for extratropical transition completion can vary considerably. For the storms that they examined, this ranged from 12–168hours. To get a sense of the temporal evolution relative to extratropical transition completion, we examine composite tropicalcyclone speed and acceleration as a function of time. For this, we only considered those Atlantic storms during 1966–2019 thatwere classified as tropical at some time and subsequently underwent extratropical transition. Of the 689 candidate storms thatpassed the three-day threshold, 18 storms were never classified as tropical . Of the remaining 671 storms, 274 were eventuallyclassified as extratropical . This yields a climatological extratropical transition fraction of 41%. However, in the data record, afew instances exist where a storm was flagged as extratropical earlier than tropical . If we remove these instances, the extra-tropical transition fraction slightly reduces to ≈ . These estimates are lower than the fraction of 44% during 1979–2017in Bieli et al. (2019), and 46% during 1950-1993 in Hart and Evans (2001). The mean and median latitude of extratropicaltransition completion in our data set were, respectively, 40.5 ◦ N and 41.5 ◦ N. This is consistent with Hart and Evans (2001)who found that the highest frequency of extratropical transition in the Atlantic occurs between the latitudes of 35 o N–45 o N. Figure 9.
Composite speed and accelerations relative to time of (a) Extratropical transition; and (b) Recurvature. A single pass of 5-pointrunning average was applied to the speed and tangential acceleration curves. Two passes of the same filter were applied to the curvatureacceleration. ig. 9a shows the composite accelerations and speed relative to the time of extratropical transition. Hour 0 is defined as thefirst instance in the IBTrACS where the storm nature is designated as extratropical transition. We interpret this as the nearesttime after extratropical transition has been completed. In an ensemble-averaged sense, the forward speed of transitioningtropical storms is seen to reach its peak around the time of extratropical transition completion. The tangential acceleration peaksabout 18 hours prior to that. The curvature acceleration appears to steadily increase up to the time of extratropical transitionand stabilizes thereafter. The point here is that the peak tangential acceleration of tropical cyclones precedes extratropicaltransition completion. The rapid increase in the speed prior to extratropical transition completion time is a direct outcome ofthe interaction with the extratropical baroclinic wavepacket. In the previous section, we found that the composite synoptic-scale flow associated with rapid curvature acceleration closelymatches the pattern associated with recurving of tropical cyclones (Aiyyer, 2015). To further explore this connection, Fig.9b shows the acceleration and speed composite timeseries relative to recurvature. We follow the method described in Aiyyer(2015) to determine the location of recurvature. A total of 653 recurvature points were found for Atlantic tropical storms overthe period 1966–2019. Note that a given storm could have more than one instance of recurvature. Fig 9b confirms that, in anensemble average sense, the time of recurvature is associated with the highest curvature acceleration. Furthermore, it is alsoassociated with the lowest forward speed and a period of rapid increase in tangential acceleration.
We first examine the trends in the annual-mean translation speed to place our results within the context of recent studies oftropical cyclone motion. As noted in the introduction, Kossin (2018) found a decreasing trend in annual-mean tropical cyclonespeed during 1949–2018 over most of the globe. That study considered all storms in the IBTrACS dataset as long as theysurvived at least three days. We revisit this for the Atlantic and test the sensitivity of the trend when we exclude non-tropicalsystems. The rationale for this was discussed earlier in section 2. Figure 10 shows the annual-mean speed of tropical cyclonesfor two categories: All storms (grey) and storms excluding NR and extratropical transition designations (orange). Panel (a)shows this for the entire Atlantic and Panel (b) for the 20–40 o N band.Table 2 contains the trends calculated using linear regression and the Thiel-Sen estimate. We also include the trends for1949–2016 to compare our calculations with Kossin (2018). For 1949–2016, when we consider all tropical cyclones, the linearregression and Thiel-Sen estimates of the trends in annual-mean speed are − . and − . km hr − year − . These arepractically identical to the value of − . km hr − year − reported by Kossin (2018). However, the trend for the satelliteera over the entire basin switches to a positive value of ≈ . km hr − year − . The sensitivity to the choice of the yearsis consistent with Lanzante (2019) who showed that the negative trend of the annual-mean speed over the longer period1949–2016 was reduced in magnitude by accounting for the change points such as those associated with the advent of satellite- igure 10. Annual-mean speed and linear trend for (a) The entire Atlantic; and (b) 20–50 o N latitude band. The grey curve is for all stormsin the IBTRaCS dataset while the orange curve excludes instances when the storm was classified as ET or NR.
Table 2.
Trends in Speed (km hr − year − ). All storms refers to all instances of a system recorded in the IBTraCs. ET refers to storm naturedesignated as extratropical, while NR refers to instances when the storm nature was not recorded.1966–2019 1949–2019 1949–2016LR MK-TS LR MK-TS LR MK-TSTrend p-value Trend p-value Trend p-value Trend p-value Trend p-value Trend p-value
Atlantic (All storms)
Full basin 0.029 0.19 0.028 0.15 -0.016 0.28 -0.016 0.32 -0.019 0.24 -0.021 0.2520–50 0.041 0.15 0.035 0.11 -0.011 0.56 -0.019 0.29 -0.012 0.55 -0.023 0.26
Atlantic (Excluding ET,NR)
Full basin -0.007 0.70 -0.008 0.62 -0.004 0.77 -0.007 0.48 -0.002 0.90 -0.006 0.6320–50 o N 0.008 0.76 0.002 0.93 <0.001 1.00 -0.009 0.52 0.005 0.79 -0.007 0.72 based weather monitoring. When we remove non-tropical data, the trends for various periods and regions are generally lower.Furthermore, none of the trends shown in Table 2 can be deemed significant if we use a p-value of 0.05 as the cut-off. Wereturn to this point in the following section. igure 11. Cumulative distributions of tangential acceleration (km hr − day − ) for July-October, 1966–2019 for (a) The entire Atlantic; (b)0–20 o N; and (c) 20–50 o N Figure 12.
Cumulative distributions of curvature acceleration (km hr − day − ) for July-October, 1966–2019 for (a) The entire Atlantic; (b)0–20 o N; and (c) 20–50 o N Figure 11 compares the cumulative probability distribution (CDF) of tangential accelerations over three 10-year periods: 1966–1975, 1988-1997, and 2010–2019. The data covers the peak hurricane season: July–October (JASO). Three Atlantic regionsare shown - the entire basin, 0–20 o N, and 20–50 o N. In all three CDFs, the lower and upper-tail probability of the distributionappear to show a shift towards the median. The direction of the shift indicates a reduction in the frequency of both rapidtangential acceleration ( a ≥ km/hr day − ) and rapid tangential deceleration ( a ≤ − km/hr day − ) from the earlier torecent decades. This is most pronounced over the 20–50 o N latitude band. The CDFs for curvature acceleration (bottom row ofFig. 12) show a similar shift towards less frequent rapid acceleration. The CDF over the entire year shows similar shifts. Whenwe consider a smaller subset of months, we find that the shifts are more pronounced when we omit October and November(not shown). n the preceding sections, we showed that rapid acceleration or deceleration of tropical cyclones are typically associatedwith interactions with the extratropical baroclinic stormtrack. The attendant synoptic-scale pattern are distinct in the phasingof the tropical cyclone and the extratropical wavepacket. It is of interest to determine if the shift in CDFs of acceleration (Figs.11, 12) are related to long-term trends. The motivation being that it can inform us about potential changes in the nature oftropical cyclone-baroclinic stormtrack interaction. Given that we are interested in the long-term trends of rapid accelerationand deceleration – i.e., the tails of the probability distribution – we use quantile regressions (QR) as developed by Koenkerand Bassett (1978). QR is a useful tool to model the behavior of the entire probability distribution and has been used in diversefields and applications (e.g., Koenker and Hallock, 2001). In atmospheric sciences, QR has been applied to examine trendsin extreme precipitation and temperature (e.g., Koenker and Schorfheide, 1994; Barbosa et al., 2011; Gao and Franzke, 2017;Lausier and Jain, 2018; Passow and Donner, 2019).The standard form of the simple linear regression model for a response variable Y in terms of its predictor X is written as: µ { Y | X } = β o + β X (5)Where, µ { Y | X } is the conditional mean of Y given a variable X, and β o and β are, respectively, the intercept and the slope.This linear model fits the mean of a response variable under the assumption of constant variance over the entire range of thepredictor. However, when the data is heteroscedastic, and there is interest in characterizing the entire distribution - and not justthe mean - QR is more appropriate and insightful. The standard form of QR is written as follows (e.g., Lausier and Jain, 2018): Y ( τ | X ) = β ( τ ) o + β ( τ )1 X + (cid:15) ( τ ) (6)where Y ( τ | X ) denotes the conditional estimate of Y at the quantile τ for a given X. By definition, < τ < . In our case,Y is the timeseries vector of either acceleration or speed, and X is the vector comprised of the dates of the individual stormpositions. Here, β ( τ ) o and β ( τ )1 denote the intercept and slope, while the (cid:15) ( τ ) denotes the error. As noted in previous studiescited above, QR does not make any assumption about the distribution of parameters and is known to be relatively robust tooutliers. To determine the trends, we fit the quantile regression model for a range of quantiles between .05 and .95. Insteadof calculating annual averages to get one value of acceleration or speed per year, we retain all of the individual values for thetropical cyclones. The corresponding time for each data point is assigned as a fractional year.Figure 13 shows the results of QR for tangential acceleration. The panels on the left show the acceleration (light bluecircles) at individual track locations from 1966-2019 (ET and NR excluded). The dashed magenta lines are the linear fits forthe quantiles ranging from 0.5 to 0.95. The panels on the right show the slope (trend; km/hr day − year − ) of the linear line asa function of the quantile. These figures include the best fit using ordinary least squares (OLS; red line) that models the meanof the distribution. Also included are the associated 95% confidence bounds. The top row includes data from all months overthe entire Atlantic. The middle row is for 20–50 o N over the peak tropical cyclone months (July-October), and the bottom rowrestricts the data to August-September. The latter two illustrate some of the sensitivity to the choice of domain and monthsof analysis. They also focus our attention on the region where tropical cyclones are most likely to interact with extratropical igure 13. Quantile regressions of tangential acceleration for regions and months shown on the panels: The left columns show the accel-eration (light blue circles) for all TCs (NR and ET excluded) as a function of time. The dotted magenta lines show the linear fits quantilesranging from 0.05 to 0.95. The red line shows the ordinary least square fit for the mean. The right columns show the estimates of the slope(i.e., the trend in km/hr day − year − ) for each quantile, along with the 95% confidence band (dotted magenta). Also shown is the ordinaryleast square estimate of trend (red line) and its 95% confidence band for each quantile. systems. The corresponding numerical values are shown in Table 3. Recall from Table 1 that the median value ( τ = . ) of able 3. Quantile trends of tangential acceleration over 1966-2019 (km hr − day − year − ) for months and regions labeled below. Trendsof magnitude below . are not reported.Entire Atlantic 20–50 o N 20–50 o NAll Months Jul-Oct Aug-Sep τ Trend Change p 95% Conf Trend Change p 95% conf. Trend Change p 95% conf.OLS 0.01 – 0.01 0.00, 0.03 0.01 10 0.50 -0.01, 0.02 , – 4 0.83 -0.01, 0.020.05 0.12 23 <0.01 0.09, 0.16 0.07 13 <0.01 0.02, 0.11 0.14 25 <0.01 0.09, 0.200.10 0.06 18 <0.01 0.04, 0.08 0.03 10 0.03 0.00, 0.07 0.07 19 <0.01 0.03, 0.100.15 0.03 12 <0.01 0.01, 0.05 – 0 0.98 -0.02, 0.02 0.01 5 0.35 -0.02, 0.040.20 0.02 9 0.02 <0.01, 0.03 – -2 0.76 -0.02, 0.02 – 2 0.70 -0.02, 0.030.30 0.01 8 0.19 -0.00, 0.02 — 4 0.60 -0.01, 0.02 0.01 6 0.53 -0.01, 0.020.50 0.03 – <0.01 0.02, 0.03 0.03 – <0.01 0.02, 0.04 0.01 – 0.04 0.00, 0.030.70 0.02 18 <0.01 0.01, 0.03 0.03 19 <0.01 0.01, 0.05 – 2 0.75 -0.02, 0.020.80 0.02 9 0.01 0.00, 0.04 0.02 9 0.04 0.00, 0.05 -0.01 -5 0.40 -0.04, 0.020.85 – -2 0.73 -0.02, 0.02 0.01 1 0.64 -0.02, 0.04 -0.04 -11 0.05 -0.07, -0.000.90 -0.02 -5 0.13 -0.04, 0.01 -0.03 -7 0.09 -0.07, 0.00 -0.07 -15 <0.01 -0.12, -0.020.95 -0.11 -18 <0.01 -0.15, -0.07 -0.14 -18 <0.01 -0.20, -0.07 -0.17 -23 <0.01 -0.25, -0.10 tangential acceleration is a small positive number. As such, τ < . is indicative of deceleration, while τ ≥ . is indicative ofacceleration.From Fig. 13 and Table 3, we note that the OLS estimate of the trend is weakly positive when data from all months over theentire Atlantic are considered. However, the OLS estimate is not statistically significant for the 20–50 o N region. As expected,the regression quantiles show a fuller picture. The slopes of the individual quantiles provide an estimate of the trends of thespecific portions of the probability distribution. The key finding here is that the magnitudes of both rapid deceleration and rapidacceleration show a statistically-significant reducing trend. This is reflected in the positive slope for τ ≤ . and negative slopefor τ ≥ . (Table 3 and left columns of Fig. 13 ). It also appears that the trends for τ < . are generally positive, implyinga reduction in the magnitude of tangential deceleration at all quantile thresholds. On the other hand, the positive slopes seenfor . < τ < . suggest that there is an increasing trend in the values of tangential acceleration that are closer to the median.This shift towards less-extreme acceleration is noted in all three regional categories, albeit with varying degree of statisticalsignificance. The trends in these quantiles are, however, weaker compared with those of the tails. From the right column ofFig. 13, it is clear that the trends of the tails of the distribution are significant and fall outside the 95% confidence bounds ofthe OLS estimate of the trend.The QR results for curvature acceleration (Fig. 14 and Table 4) show statistically significant, weak positive trends for <τ < . . However, the trends switch to increasingly negative values above the median. As in the case of tangential acceleration, igure 14. As in Fig. 13, but for curvature acceleration. the decelerating trends in the upper quantiles of the distribution ( τ ≥ . are statistically significant and outside the 95%confidence bounds of the OLS estimate of the trend in the mean.For completeness, we also show the corresponding QR results for translation speed (Fig. 15 and Table 5). The OLS estimateof the trend is nearly the same value as it was for the annual-mean speeds. However, it is now statistically significant. Thechange in the p-value reflects the fact that the sample size is much higher since the data is not averaged annually. When we able 4. Quantile trends of curvature acceleration over 1966-2019 (km hr − day − year − ) for months and regions labeled below. Trendsof magnitude below . are not reported.Entire Atlantic 20–50 o N 20–50 o NAll Months Jul-Oct Aug-Sep τ Trend Change p 95% Conf Trend Change p 95% conf. Trend Change p 95% conf.OLS -0.02 -5 <0.01 [-0.03, -0.01] -0.04 -11 <0.01 [-0.06, -0.02] -0.06 -15 <0.01 [-0.08, -0.04]0.05 0.01 42 <0.01 0.01, 0.01 0.01 40 <0.01 0.00, 0.01 0.01 27 0.04 0.00, 0.010.10 0.01 32 <0.01 0.01, 0.02 0.01 25 <0.01 0.00, 0.02 0.01 16 0.05 0.00, 0.010.15 0.01 23 <0.01 0.01, 0.02 0.01 14 0.01 0.00, 0.01 <0.01 7 0.25 -0.00, 0.010.20 0.01 17 <0.01 0.01, 0.02 0.01 13 <0.01 0.00, 0.02 0.01 7 0.22 -0.00, 0.010.30 0.02 15 <0.01 0.01, 0.02 0.01 10 0.01 0.00, 0.02 0.01 7 0.11 -0.00, 0.020.50 0.02 8 <0.01 0.01, 0.03 <0.01 1 0.55 -0.01, 0.02 -0.01 -4 0.37 -0.02, 0.010.70 – 0 0.72 -0.01, 0.02 -0.03 -7 0.01 -0.05, -0.01 -0.04 -11 <0.01 -0.06, -0.020.80 -0.03 -7 <0.01 -0.05, -0.01 -0.06 -11 <0.01 -0.08, -0.03 -0.08 -15 <0.01 -0.11, -0.050.85 -0.06 -11 <0.01 -0.08, -0.03 -0.09 -14 <0.01 -0.13, -0.05 -0.12 -20 <0.01 -0.16, -0.080.90 -0.09 -13 <0.01 -0.12, -0.06 -0.12 -15 <0.01 -0.16, -0.07 -0.18 -22 <0.01 -0.23, -0.120.95 -0.15 -15 <0.01 -0.20, -0.09 -0.22 -20 <0.01 -0.31, -0.13 -0.33 -29 <0.01 -0.43, -0.22
Table 5.
Quantile trends of translation speed over 1966-2019 (km hr − year − ) for months and regions labeled below. Trends of magnitudebelow . are not reported.Entire Atlantic 20–50 o N 20–50 o NAll Months Jul-Oct Aug-Sep τ Trend Change p 95% Conf Trend Change p 95% conf. Trend Change p 95% conf.OLS -0.01 -2 0.05 -0.01, 0.0 -0.01 -3 0.04 -0.02, -0.0 -0.04 -9 <0.01 -0.05, -0.020.05 – -13 <0.01 -0.02, -0.01 -0.02 -17 <0.01 -0.03, -0.01 -0.02 -17 <0.01 -0.03, -0.010.10 – -6 0.02 -0.01, -0.00 -0.01 -8 0.02 -0.02, -0.00 -0.01 -9 0.04 -0.02, -0.000.15 – -2 0.46 -0.01, 0.00 – -1 0.98 -0.01, 0.01 – -1 0.93 -0.01, 0.010.20 – -2 0.52 -0.01, 0.01 – -2 0.66 -0.01, 0.01 – -2 0.73 -0.01, 0.010.30 – -1 0.61 -0.01, 0.01 -0.01 -3 0.25 -0.02, 0.00 – -7 0.01 -0.03, -0.000.50 – -1 0.99 -0.01, 0.01 -0.01 -5 0.01 -0.02, -0.00 -0.04 -12 <0.01 -0.05, -0.030.70 – -2 0.34 -0.01, 0.01 -0.01 -3 0.18 -0.02, 0.00 -0.04 -10 <0.01 -0.06, -0.030.80 -0.01 -2 0.32 -0.02, 0.01 -0.02 -5 0.02 -0.04, -0.00 -0.07 -13 <0.01 -0.09, -0.050.85 -0.03 -5 <0.01 -0.04, -0.01 -0.03 -5 0.01 -0.05, -0.01 -0.09 -15 <0.01 -0.12, -0.070.90 -0.03 -4 <0.01 -0.04, -0.01 – – 0.96 -0.03, 0.03 -0.09 -14 <0.01 -0.13, -0.060.95 0.02 2 0.20 -0.01, 0.05 -0.01 -2 0.70 -0.07, 0.04 -0.08 -10 <0.01 -0.14, -0.03 igure 15. As in Fig. 13, but for translation speed. consider the entire Atlantic, estimated trends from QR are nearly the same as the OLS trend with the exception of τ = 0 . .However, when we consider the subsets of the data for 20–50 o N , there are some notable differences. In particular, for August-September, the trends in the fastest translation speeds are even more negative.Tables 3, 4 and 5 also include the percent changes defined using the first and last value in the linear fit over the period1966–2019. If we subjectively assume that a statistically significant change of magnitude at least 10% over the past 54 years an be deemed robust , then the key outcome of the QR is the following: The trends in the tails of the distribution of theaccelerations and speeds are most robust for the August-September months. For the extratropical region (20–50 o N ), bothrapid tangential acceleration and deceleration show robust reductions. This indicates a general narrowing of the tangentialacceleration distribution over time. The curvature acceleration shows an increase for the lower quantiles and reduction forthe upper. This suggests a shift in the curvature acceleration towards smaller values, consistent with the results for tangentialacceleration. Forward speed shows mostly reducing trends for both upper and lower tails, indicating that extremes in speedsare reducing over time. Ensemble-average composites of atmospheric fields show distinct synoptic-scale patterns when they are categorized on thebasis of the acceleration of tropical cyclones. The composites for rapid tangential acceleration outside the deep tropics depicta synoptic-scale pattern that is consistent with the extratropical transition of tropical cyclones. This is unsurprising since it isgenerally known that tropical storms speed up during extratropical transition. The novel aspect here is that we have recoveredthe signal of extratropical transition from the perspective of acceleration. The composites show a poleward moving tropicalcyclone that is straddled by an upstream trough and a downstream ridge. Subsequently, the tropical cyclone merges with theextratropical wavepacket ahead of the trough in an arrangement that is conducive for further baroclinic development. Featurescommonly associated with extratropical transition such as the downstream ridge-building, amplification of the upper-level jetstreak and downstream development can be clearly seen in the composite maps (Fig. 3 and Fig. 4a-e). The composites forrapid curvature acceleration also show the impact of the phasing of the tropical cyclone and the extratropical wavepacket. Forthis category, we recover a synoptic-scale pattern that is similar to the one obtained in composites based on the recurvature oftropical cyclones (Aiyyer, 2015).In contrast, the composite fields for rapid tangential and curvature deceleration show a tropical cyclone that approaches anextratropical ridge. The upstream trough remains at a distance and the tropical cyclone does not merge with the extratropicalwavepacket, but instead remains equatorward of it over the following few days. The tropical cyclone and the extratropicalridge – at least locally – can be viewed as a vortex dipole. The combined system remains relatively stationary compared tothe progressive pattern for rapid acceleration. This arrangement qualitatively resembles a dipole block — an important modeof persistent anomalies in the atmosphere (e.g., McWilliams, 1980; Pelly and Hoskins, 2003). The canonical dipole block isdepicted as a vortex pair comprised of a warm anticyclone and a low-latitude cut-off cyclone (e.g., Haines and Marshall, 1987;McTaggart-Cowan et al., 2006). The dynamics of blocked flows are rich and the subject of a variety of theories that are farfrom settled (e.g. Woollings and Barriopedro, 2018). In the present case, the slowly propagating cyclone-anticyclone pair islikely an outcome of a fortuitous phasing of the tropical cyclone and the extratropical ridge.Our study takes a complementary view of the extratropical interaction described in Riboldi et al. (2019). In particular, wenote the close correspondence between their composite sea-level pressure and potential vorticity composites (their Fig. 10) fordecelerating troughs and our geopotential composites for rapidly decelerating tropical cyclones (right columns of Fig. 3 and ig. 6). The vortex dipole can be seen in all three figures. Riboldi et al. (2019) conditioned their composites on the basis ofthe acceleration of the upstream trough interacting with recurving typhoons in the western North Pacific. We recover the samepattern when we condition the composites based on rapidly decelerating tropical cyclones in the North Atlantic.The interaction between the tropical cyclone and the extratropical wavepacket that leads to the deceleration of the entiresystem can be viewed from a PV perspective. As noted by Riboldi et al. (2019), both adiabatic and diabatic pathways areactive in this interaction. In the former, the induced flow from the cyclonic vortex (tropical cyclone) will be westward withinthe poleward anticyclonic vortex (ridge). The induced flow from the ridge will also be westward at the location of the tropicalcyclone. The combined effect will be a mutual westward advection, and thus reduced eastward motion in the earth-relativeframe. The latter, diabatic pathway relies on the amplification of the ridge through the action of precipitating convection in thevicinity of the tropical cyclone. The negative vertical gradient of diabatic heating in upper levels of the tropical cyclone impliesthat its anticyclonic outflow is a source of low PV air (e.g., Wu and Emanuel, 1993). The advection of this low PV air by theirrotational component of the outflow and its role in ridge-building has been extensively documented (e.g., Atallah et al., 2007;Riemer et al., 2008; Grams et al., 2013; Archambault et al., 2015). Riboldi et al. (2019) also showed that the ridge is moreamplified for rapidly decelerating troughs as compared to accelerating troughs. They implicated stronger latent heating andirrotational outflow for this difference. Our composites also show a stronger ridge for rapid tropical cyclone deceleration ascompared to rapid acceleration. This can be noted by comparing the left and right columns of Figs. 4 and 7. The amplificationof the ridge also results in the slow-down of the upstream trough, and as shown by Riboldi et al. (2019) yields frequentdownstream atmospheric blocking events.The tracks of TC in the vicinity of extratropical wavetrains and subtropical ridges are sensitive to the existence of bifurcationpoints (e.g., Grams et al., 2013; Riemer and Jones, 2014), and small shifts in positions can yield different outcomes for motionand extratropical transition. Bifurcation points also exist in the case of tropical-cyclone cut-off low interactions (Pantillon et al.,2016). Our acceleration-based composites have further highlighted the impact of the phasing of the tropical cyclone and theextratropical wavepacket in mediating the interactions between them.While there is some storm-to-storm variability, in an average sense the tangential acceleration peaks 18 hours prior tocompletion of extratropical transition. Interestingly, the forward speed peaks around the time of completion of extratropicaltransition. Curvature acceleration increases rapidly prior to extratropical transition and remains nearly steady after this time.Composite time series also show that the curvature acceleration peaks at track recurvature while the forward speed is nearlyat its minimum. The tangential acceleration shows a sharp, steady increase around recurvature. This is consistent with theobservations that extratropical transition is typically completed within 2-4 days of recurvature. A related relevant questionis the following: Is rapid tangential acceleration a sign of imminent extratropical transition? We found that ≈ of thestorms that comprised the composites for rapid acceleration (Left column: Fig. 3) completed the transition within 3 days ofthe reference time (day 0). This is substantially higher than the climatological fraction of ≈ for extratropical transitionover the entire basin and storm lifetime. On the other hand, only ≈ of the storms that comprised the composites for rapiddeceleration (Right column: Fig. 3) completed the transition within a similar time range. This fraction is substantially lowerthan the climatological fraction. It is, however, consistent with the observation made earlier that the synoptic-scale pattern or rapid deceleration promotes recurvature rather than extratropical transition. Furthermore, not all recurving storms becomeextratropical (Evans et al., 2017).The composites show that rapid acceleration and deceleration can be viewed as a proxy for distinct types of tropical cyclone-extratropical interactions. It is of interest, therefore, to ascertain whether any trends in tropical cyclone accelerations can befound in the track data. For this, we use quantile regression to examine the linear trends over the entire distribution. Thetails of acceleration and translation speed show statistically significant trends. The trends for the extratropical regions of theAtlantic (20–50 o N) are most robust for August-September. Both rapid tangential acceleration and deceleration have reducedover the past five decades. The trends for curvature acceleration show increases for the lower quantiles and decreases forthe upper quantiles. The forward speed, particularly values above the median, also shows robust decreases for the August-September months. This supports the general conclusion of Kossin (2018) that tropical cyclones have slowed down in the pastfew decades.We have not explored the physical basis for the trends discussed above. We, however, speculate that they are indicativeof systematic changes in the interaction between tropical cyclones and extratropical waves. It is, however, unclear from ourpreliminary examination whether the trends reflect changes in the frequency or some measure of the strength of the interactions.We also recognize that the notion of strength of interaction needs a firm and objective definition. Nevertheless, there are somerecent findings related to the atmospheric general circulation that may be relevant to this point. First, there is growing evidencefor a poleward shift in the stormtrack of extratropical baroclinic eddies. As noted by Tamarin and Kaspi (2017) and referencestherein, this shift has been found in both reanalysis data and climate model simulations. Separately, Coumou et al. (2015) founda robust decrease of the zonal wind and the amplitude of synoptic-scale Rossby waves in the ERA-interim reanalysis over theNorthern hemisphere during the months June-August. We hypothesize that the poleward shift of the extratropical waves andtheir weakening could potentially account for the acceleration trends reported here. This, however, needs to be examined furtherif any robust conclusion regarding attribution to climate change is to be made.
10 Conclusions
When we separate tropical cyclones on the basis of their acceleration and consider ensemble average composites of atmosphericfields (e.g., geopotential), we get two broad sets of synoptic-scale patterns. The composite for rapid tangential accelerationshows a poleward moving tropical cyclone straddled by an upstream trough and a downstream ridge. The subsequent mergerof the tropical cyclone and the developing extratropical wavepacket is consistent with the process of extratropical transition.The composite for rapid curvature acceleration shows a prominent downstream ridge that promotes recurvature. On the otherhand, the synoptic-scale pattern for rapid tangential as well as curvature deceleration takes the form of a cyclone-anticyclonedipole with a ridge directly poleward of the tropical cyclone. We note the qualitative resemblance of this arrangement withthe canonical dipole block. Some of our findings from the perspective of tropical cyclone acceleration closely match those ofRiboldi et al. (2019), who conditioned their diagnostics on the basis of the acceleration of the upstream trough and recurvingwestern north Pacific typhoons. ccelerations and speed show robust trends in the tails of their distribution. For the extratropical region of the Atlantic(20–50 o N), and particularly for the months August-September, peak acceleration/deceleration, as well as speeds of tropicalcyclones, have reduced over the past 5 decades. The reduction in the tails of the speed distribution provide complementaryevidence for a general slowing trend of tropical cyclones reported by Kossin (2018). We also suggest that the robust reductionin the tails of the acceleration distribution is indicative of a systematic change in the interaction of tropical cyclones withextratropical baroclinic waves. We have not, however, examined the underlying processes. We speculate that poleward shiftand decreasing amplitude of extratropical Rossby waves found in other studies may account for the acceleration trends. Detailedmodeling and observational studies are needed to better understand the source of these trends.
11 Data and code availability
12 Author contributions
AA wrote the computer code for all analysis and visualization, and wrote the text of the paper. TW derived the expression forthe radius of curvature, and assisted with editing the text and interpretation of the results.
13 Competing interests
The authors declare that they have no conflict of interest.
14 Financial Support
This work was supported by NSF through award
Aiyyer, A.: Recurving western North Pacific tropical cyclones and midlatitude predictability, Geophysical Research Letters, 42, 7799–7807,https://doi.org/10.1002/2015GL065082, 2015.Archambault, H. M., Keyser, D., Bosart, L. F., Davis, C. A., and Cordeira, J. M.: A Composite Perspective of the Extratropical Flow Responseto Recurving Western North Pacific Tropical Cyclones, Monthly Weather Review, 143, 1122–1141, https://doi.org/10.1175/MWR-D-14-00270.1, 2015.Atallah, E., Bosart, L. F., and Aiyyer, A. R.: Precipitation Distribution Associated with Landfalling Tropical Cyclones over the EasternUnited States, Mon. Wea. Rev., 135, 2185–+, https://doi.org/10.1175/MWR3382.1, 2007.Barbosa, S. M., Scotto, M. G., and Alonso, A. M.: Summarising changes in air temperature over central europe by quantile regression andclustering, Nat. Hazards Earth Syst. Sci., 11, 3227–3233, https://doi.org/10.5194/nhess-11-3227-2011, 2011.Bieli, M., Camargo, S. J., Sobel, A. H., Evans, J. L., and Hall, T.: A Global Climatology of Extratropical Transition. Part I: Characteristicsacross Basins, Journal of Climate, 32, 3557–3582, https://doi.org/10.1175/JCLI-D-17-0518.1, 2019.Bosart, L. F. and Lackmann, G. M.: Postlandfall Tropical Cyclone Reintensification in a Weakly Baroclinic Environment:A Case Study of Hurricane David (September 1979), Monthly Weather Review, 123, 3268, https://doi.org/10.1175/1520-0493(1995)123<3268:PTCRIA>2.0.CO;2, 1995.Chan, J. C.: THE PHYSICS OF TROPICAL CYCLONE MOTION, Ann. Rev. of Fluid Mech., 37, 99–128,https://doi.org/10.1146/annurev.fluid.37.061903.175702, https://doi.org/10.1146/annurev.fluid.37.061903.175702, 2005.Chan, J. C. L. and Williams, R. T.: Analytical and Numerical Studies of the Beta-Effect in Tropical Cyclone Motion. Part I: Zero Mean Flow.,Journal of Atmospheric Sciences, 44, 1257–1265, https://doi.org/10.1175/1520-0469(1987)044<1257:AANSOT>2.0.CO;2, 1987.Coumou, D., Lehmann, J., and Beckmann, J.: The weakening summer circulation in the Northern Hemisphere mid-latitudes, Science, 348,324–327, https://doi.org/10.1126/science.1261768, 2015.Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P.,Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L.,Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Mor-crette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configurationand performance of the data assimilation system, Quart. J. Roy. Meteor. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011.Emanuel, K.: 100 Years of Progress in Tropical Cyclone Research, vol. 59, pp. 15.1–15.68, https://doi.org/10.1175/AMSMONOGRAPHS-D-18-0016.1, 2018.Evans, C., Wood, K. M., Aberson, S. D., Archambault, H. M., Milrad, S. M., Bosart, L. F., Corbosiero, K. L., Davis, C. A., Dias Pinto,J. R., Doyle, J., Fogarty, C., Galarneau, Thomas J., J., Grams, C. M., Griffin, K. S., Gyakum, J., Hart, R. E., Kitabatake, N., Lentink,H. S., McTaggart-Cowan, R., Perrie, W., Quinting, J. F. D., Reynolds, C. A., Riemer, M., Ritchie, E. A., Sun, Y., and Zhang, F.: TheExtratropical Transition of Tropical Cyclones. Part I: Cyclone Evolution and Direct Impacts, Monthly Weather Review, 145, 4317–4344,https://doi.org/10.1175/MWR-D-17-0027.1, 2017.Gao, M. and Franzke, C. L. E.: Quantile Regression-Based Spatiotemporal Analysis of Extreme Temperature Change in China, J. Climate,30, 9897–9914, https://doi.org/10.1175/JCLI-D-17-0356.1, 2017.Grams, C. M., Jones, S. C., and Davis, C. A.: The impact of Typhoon Jangmi (2008) on the midlatitude flow. Part II: Downstream evolution,Quarterly Journal of the Royal Meteorological Society, 139, 2165–2180, https://doi.org/10.1002/qj.2119, 2013. aines, K. and Marshall, J.: Eddy-forced coherent structures as a prototype of atmospheric blocking, Quarterly Journal of the Royal Meteo-rological Society, 113, 681–704, https://doi.org/10.1002/qj.49711347613, 1987.Hart, R. E.: A Cyclone Phase Space Derived from Thermal Wind and Thermal Asymmetry, Monthly Weather Review, 131, 585,https://doi.org/10.1175/1520-0493(2003)131<0585:ACPSDF>2.0.CO;2, 2003.Hart, R. E. and Evans, J. L.: A Climatology of the Extratropical Transition of Atlantic Tropical Cyclones., Journal of Climate, 14, 546–564,https://doi.org/10.1175/1520-0442(2001)014<0546:ACOTET>2.0.CO;2, 2001.Hart, R. E., Evans, J. L., and Evans, C.: Synoptic Composites of the Extratropical Transition Life Cycle of North Atlantic Tropical Cyclones:Factors Determining Posttransition Evolution, Monthly Weather Review, 134, 553, https://doi.org/10.1175/MWR3082.1, 2006.Hodanish, S. and Gray, W. M.: An Observational Analysis of Tropical Cyclone Recurvature, Monthly Weather Review, 121, 2665,https://doi.org/10.1175/1520-0493(1993)121<2665:AOAOTC>2.0.CO;2, 1993.Jones, S. C., Harr, P. A., Abraham, J., Bosart, L. F., Bowyer, P. J., Evans, J. L., Hanley, D. E., Hanstrum, B. N., Hart, R. E.,Lalaurette, F., Sinclair, M. R., Smith, R. K., and Thorncroft, C.: The Extratropical Transition of Tropical Cyclones: ForecastChallenges, Current Understanding, and Future Directions, Weather and Forecasting, 18, 1052–1092, https://doi.org/10.1175/1520-0434(2003)018<1052:TETOTC>2.0.CO;2, 2003.Keller, J. H., Grams, C. M., Riemer, M., Archambault, H. M., Bosart, L., Doyle, J. D., Evans, J. L., Galarneau, Thomas J., J., Griffin, K.,Harr, P. A., Kitabatake, N., McTaggart-Cowan, R., Pantillon, F., Quinting, J. F., Reynolds, C. A., Ritchie, E. A., Torn, R. D., and Zhang, F.:The Extratropical Transition of Tropical Cyclones. Part II: Interaction with the Midlatitude Flow, Downstream Impacts, and Implicationsfor Predictability, Monthly Weather Review, 147, 1077–1106, https://doi.org/10.1175/MWR-D-17-0329.1, 2019.Klein, P. M., Harr, P. A., and Elsberry, R. L.: Extratropical Transition of Western North Pacific Tropical Cyclones: AnOverview and Conceptual Model of the Transformation Stage, Weather and Forecasting, 15, 373–396, https://doi.org/10.1175/1520-0434(2000)015<0373:ETOWNP>2.0.CO;2, 2000.Knapp, K. R., Kruk, M. C., Levinson, D. H., Diamond, H. J., and Neumann, C. J.: The International Best Track Archive for ClimateStewardship (IBTrACS), Bulletin of the American Meteorological Society, 91, 363–376, https://doi.org/10.1175/2009BAMS2755.1, 2010.Koenker, R. and Bassett, G.: Regression quantiles, Econometrica, 46, 33–50, https://doi.org/doi: 10.2307/1913643, 1978.Koenker, R. and Hallock, K.: Quantile Regression, Journal of Economic Perspectives—Volume 15, Number 4—Fall 2001—Pages 143–156,15, 143–156, 2001.Koenker, R. and Schorfheide, F.: Quantile spline models for global temperature change, Climatic Change, 28, 395–404,https://doi.org/10.1007/BF01104081, 1994.Kofron, D. E., Ritchie, E. A., and Tyo, J. S.: Determination of a Consistent Time for the Extratropical Transition of Tropical Cyclones. PartII: Potential Vorticity Metrics, Monthly Weather Review, 138, 4344–4361, https://doi.org/10.1175/2010MWR3181.1, 2010.Kossin, J. P.: A global slowdown of tropical-cyclone translation speed, nature, 558, 104–107, https://doi.org/10.1038/s41586-018-0158-3,2018.Landsea, C. W.: Counting Atlantic Tropical Cyclones Back to 1900, EOS Transactions, 88, 197–202, https://doi.org/10.1029/2007EO180001,2007.Lanzante, J. R.: Uncertainties in tropical-cyclone translation speed, nature, 570, E6–E15, https://doi.org/10.1038/s41586-019-1223-2, 2019.Lausier, A. M. and Jain, S.: Overlooked Trends in Observed Global Annual Precipitation Reveal Underestimated Risks, Scientific Reports,8, 16746, https://doi.org/10.1038/s41598-018-34993-5, 2018.