Rapid intensification and the bimodal distribution of tropical cyclone intensity

Abstract

The severity of a tropical cyclone (TC) is often summarized by its lifetime maximum intensity (LMI), and the climatological LMI distribution is a fundamental feature of the climate system. The distinctive bimodality of the LMI distribution means that major storms (LMI >96 kt) are not very rare compared with less intense storms. Rapid intensification (RI) is the dramatic strengthening of a TC in a short time, and is notoriously difficult to forecast or simulate. Here we show that the bimodality of the LMI distribution reflects two types of storms: those that undergo RI during their lifetime (RI storms) and those that do not (non-RI storms). The vast majority (79%) of major storms are RI storms. Few non-RI storms (6%) become major storms. While the importance of RI has been recognized in weather forecasting, our results demonstrate that RI also plays a crucial role in the TC climatology.

ARTICLE

Received 18 May 2015 | Accepted 5 Jan 2016 | Published 3 Feb 2016

The severity of a tropical cyclone (TC) is often summarized by its lifetime maximum intensity (LMI), and the climatological LMI distribution is a fundamental feature of the climate system. The distinctive bimodality of the LMI distribution means that major storms (LMI 496 kt) are not very rare compared with less intense storms. Rapid intensication (RI) is the dramatic strengthening of a TC in a short time, and is notoriously difcult to forecast or simulate. Here we show that the bimodality of the LMI distribution reects two types of storms: those that undergo RI during their lifetime (RI storms) and those that do not (non-RI storms). The vast majority (79%) of major storms are RI storms. Few non-RI storms (6%) become major storms. While the importance of RI has been recognized in weather forecasting, our results demonstrate that RI also plays a crucial role in the TC climatology.

DOI: 10.1038/ncomms10625 OPEN

Rapid intensication and the bimodal distribution of tropical cyclone intensity

Chia-Ying Lee1, Michael K. Tippett2,3, Adam H. Sobel2,4 & Suzana J. Camargo4

1 International Research Institute of Climate and Society, Columbia University, Palisades, New York 10964, USA. 2 Department of Applied Physics and Applied Mathematics, Columbia University, New York 10027, USA. 3 Center of Excellence for Climate Change Research, Department of Meteorology, King Abdulaziz University, Jeddah 21589, Saudi Arabia. 4 Division of Ocean and Climate Physics, Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York 10964, USA. Correspondence and requests for materials should be addressed to C.-Y.L. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:10625 | DOI: 10.1038/ncomms10625 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10625

The question of how climate change will affect tropical cyclone (TC) activity has drawn considerable attention in the past two decades14. The current expectation is that we

can expect a small increase in the global frequency of intense storms along with a small reduction in the total number of storms4. However, details of how the TC intensity distribution may change remain uncertain. This uncertainty reects both the difculty of simulating the most intense storms in climate change projections4, as well as an incomplete understanding of what determines the climatological intensity distribution in the current climate5. Improved understanding of the TC intensity distribution in the current climate seems necessary if we are to understand TC intensity changes in a warming climate. Among various metrics of TC intensity, lifetime maximum intensity (LMI) is an integrated statistic of TC intensication, and its distribution represents a fundamental property of the TC climatology68.

Several authors have noted that the probability density function (PDF) of global LMI is bimodal. The LMI PDF from International Best-Track Archive for Climate Stewardship for the period 19752007 has two local maxima around 40 and 100 kt, and a local minimum at 65 kt (ref. 9). The LMI PDF from a more temporally consistent global data for the period 19822009 shows the rst maximum at 50 kt and the second one around 120 kt (ref. 10). Considering individual basins, the LMI PDF has its rst peak at 50 kt for the North Atlantic, eastern and western North Pacic storms and its second maximum at 110 and 90 kt for eastern and western North Pacic storms11. The LMI distribution of Atlantic storms shows no clear secondary maximum, but substantial right skewness11. While the precise locations of the LMI PDF peaks vary with the data set and basin, the bimodal feature is quite robust and impossible to overlook. The existence of a second mode means that major storms (sustained winds 496 kt, that is, categories 35 TCs in Safr-Simpson Hurricane

Wind Scale) are much less rare than would be expected from the behaviour of the LMI distribution at lower values. The nonmonotonic behaviour of the LMI distribution is in contrast to that of the intensication rate distribution, which is exponential12.

Although no complete explanation has been offered in the literature for the bimodality of the LMI distribution, there have been suggestions and indications of what processes might be involved. Uncertainties in the best-track winds is one of them. For example, the widely used Dvorak technique for estimating TC intensity from satellite imagery has low resolution at higher intensities, with only one bin in the category 3 range of the Safr-Simpson scale. It has been argued that this may result in an articially low number of category 3 hurricanes (LMI between 96112 kt) in the Atlantic13. A recent study, on the other hand, proposed a parameterization of the ratio of surface exchange coefcients Ck/Cd, which appears in potential intensity theory14,15, as a function of wind speed with a local maximum around 115 kt, and speculated that such a maximum would be favourable for rapid intensication (RI) and might explain the bimodal distribution of LMI5. In addition to the slope change in Ck/Cd, observations have suggested an association between RI and eye formation16. Organized convective heating results in increased vortex efciency1720, which can potentially increase RI probability and cause intense systems. Perhaps most directly relevant, storms reaching the highest intensities in the western North Pacic21 and North Atlantic22 typically do so after undergoing RI at least once.

Here we extend these results to all basins, and show explicitly that RI explains the bimodality of the LMI distribution. We separate storms into two groups: those that undergo RI during their lifetime (RI storms) and those that do not (non-RI storms).

We nd that the bimodality of the LMI distribution reects the mixture of these two unimodal distributions, with the higher intensity peak consisting mostly of storms which have undergone RI at some point. In other words, the LMI distribution is unimodal when RI storms are excludedRI storms are responsible for the bimodality of the LMI distribution. Various thresholds have been used to dene RI23,24. Here we dene RI as an increase of at least 35 kt in the maximum sustained surface wind over a period of 24 h or less. We nd this denition of RI to be the most effective in separating the LMI distribution into two unimodal distributions. Sensitivity of the results to other RI thresholds, such as 30 kt, is discussed later.

ResultsThe bimodal LMI distribution. The global distribution of LMI for the period 19812012 has a peak at 45 kt and an indication of a secondary maximum around 120 kt (Fig. 1). Basin distributions of LMI are bimodal with local maxima around 45 and 120135 kt except for the Atlantic where there is only a hint of a second maximum (Fig. 2), similar to what is shown in the literature11. The observed LMI distribution shows that category 3 and 4 storms are more common than category 1 and 2 storms. Most types of natural hazards become more rare as they become more extreme, for example, earthquakes25 and tornadoes26. TCs, as measured by LMI, are unusual in having a range over which frequency increases with intensity. The second peak in the western North Pacic occurs at a higher intensity than in the other basins, consistent with the observation that the stronger storms globally occur more often in that basin27.

Relation to RI. The LMI distribution of the 2,303 non-RI storms in the global record is unimodal and forms the rst peak of the complete distribution (blue curve in Fig. 1). About 6% (141) of the non-RI storms are major storms (LMI 496 kt, categories 35), and the largest LMI value in this group is 145 kt (category 5). The second peak in the LMI distribution is formed by the 766 RI storms (red curve in Fig. 1), and 79% (603) of them are major storms. The same separation by RI in individual basins yields

12

10

8

6

4

2

0 0 50 100 150 200

LMI (kt)

PDF (%)

Global tropical cyclones (3,069)

RI 766 Non-RI 2,303

Figure 1 | Distributions of global tropical cyclone LMI. PDFs are calculated using 19812012 global tropical cyclone LMI. The grey bars show the raw data binned in 5 kt bins. The black, red and blue lines show the smoothed PDF for all storms, storms those undergo rapid intensication during their lifetime (RI storms), and those do not (non-RI storms), respectively. Smoothing is by moving average with window width of 15 kt. Total number of storms is listed in the title, while the numbers of RI and non-RI storms are given in the legend.

2 NATURE COMMUNICATIONS | 7:10625 | DOI: 10.1038/ncomms10625 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10625 ARTICLE

a

North Atlantic (487) Western North Pacific (956)

Eastern North Pacific (603) North Indian Ocean (163)

Southern Hemisphere basins (860)

224

636

b

12 10

86420 0 50 100 150 200

12 10

86420 0 50 100 150 200

12 10

864200 50 100 150 200

96

391

150

453

277

679

PDF (%)

PDF (%) PDF (%)

c d

12 10

8 6 4 2 0

0 50 100 150 200

19

144

PDF (%)

e

12 10

864200 50 100 150 200

PDF (%)

LMI (kt)

Figure 2 | Distributions of regional tropical cyclone LMI. PDFs are calculated using 19812012 tropical cyclone for individual basins: (a) North Atlantic, (b) Western North Pacic, (c) Eastern North Pacic, (d) North Indian Ocean, and (e) Southern Hemisphere basins. The black, red and blue lines show the smoothed PDFs for all, the subset of storms those undergo rapid intensication during their lifetime (RI storms), and those do not (non-RI storms), respectively. The raw data is binned in 5 kt bins and smoothing is by moving average with window width of 15 kt. The number of storms in each basin is given in the title, and the numbers of RI and non-RI storms are given in the legend.

similar results (Fig. 2). The intensication rate of 35 kt in 24 h (representing the 97th percentile of the intensication rate over that duration) is the optimal RI threshold for explaining the bimodality of the LMI distribution. Using other thresholds, such as 25, 30 or even 40 kt, does not separate the two peaks as clearly, and results a secondary maximum or a hint of it in either the RI or the non-RI LMI distributions (Supplementary Fig. 1). With the 35 kt denition, 8085% of major storms in the eastern and western North Pacic, North Indian Ocean and southern Hemisphere basins are RI storms, but only 70% of major storms in the Atlantic are RI storms. In the Atlantic, the LMI distribution of RI storms is quite different from that in other basins, being much less peaked (red lines in Fig. 2). The probability that an RI storm in the Atlantic will become a minor hurricane (categories 12) is close to the probability that it will become a major hurricane. Still, RI storms are responsible for the rightmost shoulder of the LMI distribution in the Atlantic.

Relation to storm lifetime. Storm lifetime is another factor related to LMI. The distribution of lifetime itself is unimodal

(Supplementary Fig. 2). Overall, LMI is positively correlated with storm lifetime with longer lived storms tending to have higher LMI values. The correlation of storm lifetime with LMI is 0.66 globally. Here the idea that the bimodality in LMI can be explained by variations in lifetime is tested by separating storms into two groups, longer and shorter lived, and examining the LMI distributions for the two groups, analogously to what was done for RI and non-RI storms above. Thresholds of lifetime from 7 to 13 days are tested for this purpose (Supplementary Fig. 3). The classication produces less completely unimodal LMI PDFs than were achieved using RI as the criterion, especially for storms in the North Indian Ocean and southern Hemisphere basins. The best threshold of lifetime is also basin dependent, and noticeable departures from uni-modality are found in even when that best value is used. We, therefore, conclude that while lifetime has an impact on LMI, RI is a better criteria for explaining the LMI bimodality.

DiscussionThe observed relation of RI with the bimodal LMI distribution is consistent with numerical simulation studies in which higher

NATURE COMMUNICATIONS | 7:10625 | DOI: 10.1038/ncomms10625 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10625

horizontal resolution results in both higher intensication rates and the appearance of the second LMI peak9,28. Such numerical simulations provide evidence that the observed LMI bimodality and its relationship to RI is not an artifact of the uncertainty associated with intensity estimation13,29,30. Our ndings are also in agreement with the hypothesis of ref. 5 that RI is responsible for the bimodal distribution of LMI, although we cannot comment on the link to the ratio of exchange coefcients. On the other hand, a recent study showed that RI has no apparent signature in the overall intensication rate statistics12. The intensication rate distributions of RI and non-RI storms do differ, even when the RI events are excluded, but no distinguishing features result when all storms are considered (Supplementary Fig. 4). A key difference between the distribution of intensication rates and that of LMI is that the LMI distribution accounts for the intensication behaviour over the storm lifetime, and thus is a measure of the joint (across different times) intensication rate statistics.

Another measure of the TC intensity climatology is LMI normalized by local potential intensity. The distribution of normalized LMI for storms whose peak intensities are not limited by declining potential intensity has a linear cumulative distribution with two slopes6, corresponding approximately to tropical storms (434 kt) and hurricanes (464 kt). This separation of the normalized LMI distribution appears to be distinct from the one here based on RI. The subset of storms not limited by declining potential intensity is a relatively small fraction of all storms, and the normalized LMI distribution over all storms (which we consider here) is not uniform (not shown; for example, Fig. 13 in ref. 6).

The basic physical mechanisms of RI are still not completely understood, and whether they are distinct from those of lower intensication rates remain an open question. We do not attempt to answer this question. The message here is that RI is relevant not only to short-term weather forecasting, but also to the relationship between TCs and climate. The most intense storms are those that undergo RI, and the storms that undergo RI are responsible for the observed bimodality of the LMI distribution. This nding suggests that a complete understanding of the most intense storms in either the current climate or future (or past) climates may need to include some understanding of RI. Our results also suggest that numerical models that do not simulate RI are likely to be incomplete in their representation of the LMI distribution and in the frequency of major storms. Therefore, an important research question is to what extent simulations and projections of the frequency of major storm occurrence can be accurate without either resolving RI or accounting for its absence.

Methods

Data. Best-track data from the National Hurricane Center (NHC)29,31 and the Joint Typhoon Warming Centers (JTWC)32 from 1981 to 2012 are used in this study. Best-track data include 1-min maximum sustained wind, minimum sea level pressure and location every 6 h.

Calculations. The LMI here is dened using maximum wind speed. Maximum wind speed itself is not an observed quantity, but rather estimated from in situ observations, remotely sensed estimates of winds or via satellite-based techniques29. The distribution of LMI is calculated globally, as well as for individual basins, following the denitions of the NHC and JTWC, that is, North Atlantic and eastern North Pacic (from NHC), and western North Pacic, North Indian Ocean and southern Hemisphere basins (from JTWC). Correlations between global LMI and storm lifetime are calculated using Spearmans rank correlation coefcient.

Terminology. We refer to storms with LMI 496 kt (categories 35 in Safr-Simpson Hurricane Wind Scale) as major storms.

References

1. Bengtsson, L., Botzet, M. & Esch, M. Will greenhouse-induced warming over the next 50 years lead to higher frequency and greater intensity of hurricanes? Tellus 48A, 5773 (1996).

2. Emanuel, K. A. Increasing destructiveness of tropical cyclones over the past 30 years. Nature 436, 686688 (2005).

3. Webster, P. J., Holland, G., Curry, J. A. & Chang, H.-R. Changes in tropical cyclone number, duration, and intensity in a warming environment. Science 309, 18441846 (2005).

4. Knutson, T. R. et al. Tropical cyclones and climate change. Nat. Geosci. 3, 157163 (2010).

5. Soloviev, A. V., Lukas, R., Donelan, M. A., Haus, B. K. & Ginis, I. The air-sea interface and surface stress under tropical cyclones. Sci. Rep. 4, 5306 (2014).

6. Emanuel, K. A. A statistical analysis of tropical cyclone intensity. Mon. Weather Rev. 128, 11391152 (2000).

7. Kossin, J. P., Emanuel, K. A. & Vecchi, G. A. The poleward migration of the location of tropical cyclone maximum intensity. Nature 509, 349352 (2014).

8. Park, D.-S. R., Ho, C.-H. & Kim, J.-H. Growing threat of intense tropical cyclones to East Asia over the period 1977-2010. Environ. Res. Lett. 9, 014008 (2014).

9. Manganello, J. V. et al. Tropical cyclone climatology in a 10-km global atmospheric GCM: toward weather-resolving climate modeling. J. Climate 25, 38673893 (2012).

10. Kossin, J. P., Olander, T. L. & Knapp, K. R. Trend analysis with a new global record of tropical cyclone intensity. J. Climate 26, 99609976 (2013).

11. Zhao, M., Held, I. M., Lin, S.-J. & Vecchi, G. A. Simulations of global hurricane climatology, interannual variability, and response to global warming using a 50-km resolution GCM. J. Climate 22, 66536678 (2009).

12. Kowch, R. & Emanuel, K. Are special processes at work in the rapid intensication of tropical cyclones? Mon. Weather Rev. 143, 878882 (2014).

13. Torn, R. D. & Snyder, C. Uncertainty of tropical cyclone best-track information. Weather Forecast. 27, 715729 (2012).

14. Bister, M. & Emanuel, K. A. Low frequency variability of tropical cyclone potential intensity 1. Interannual to interdecadal variability. J. Geophys. Res. 107, 4801 (2002).

15. Emanuel, K. A. Sensitivity of tropical cyclones to surface exchange coefcients and a revised steady-state model incorporating eye dynamics. J. Atmos. Sci. 52, 39693976 (1995).

16. Vigh, J. L., Knaff, J. A. & Schubert, W. H. A climatology of hurricane eye formation. Mon. Weather Rev. 140, 14051426 (2012).

17. Shapiro, L. J. & Willoughby, H. E. The response of balanced hurricanes to local-sources of heat and momentum. J. Atmos. Sci. 39, 378394 (1982).18. Hack, J. J. & Schubert, W. H. Nonlinear response of atmospheric vortices to heating by organized cumulus convection. J. Atmos. Sci. 43, 15591573 (1986).

19. Vigh, J. L. & Schubert, W. H. Rapid development of the tropical cyclone warm core. J. Atmos. Sci. 66, 33353350 (2009).

20. Nolan, D. S., Moon, Y. & Stern, D. P. Tropical cyclone intensication from asymmetric convection: energetics and efciency. J. Atmos. Sci. 64, 33773405 (2007).

21. Holliday, C. R. & Thompson, A. H. Climatological characteristics of rapidly intensifying typhoons. Mon. Weather Rev. 107, 10221034 (1979).

22. Kaplan, J. & DeMaria, M. Large-scale characteristics of rapidly intensifying tropical cyclones in the North Atlantic basin. Weather Forecast. 18, 10931108 (2003).

23. Knaff, J. et al. in Sixth International Workshop on Tropical Cyclones vol. 1.5 (WMO, Jose, Costa Rica., 2006).

24. Kaplan, J., DeMaria, M. & Knaff, J. A. A revised tropical cyclone rapid intensication index for the Atlantic and Eastern North Pacic basins. Weather Forecast. 25, 220241 (2010).

25. Hristopulos, D. T. & Mouslopoulou, V. Strength statistics and the distribution of earthquake interevent times. Physica A 392, 485496 (2013).

26. Feuerstein, B., Dotzek, N. & Grieser, J. Assessing a tornado climatology from global tornado intensity distributions. J. Climate 18, 585596 (2005).

27. Knapp, K. R., Kruk, M. C., Levinson, D. H., Diamond, H. J. & Neumann, C. J. The international best track archive for climate stewardship (IBTrACS). Bull. Am. Meteor. Soc. 91, 363376 (2010).

28. Murakami, H. et al. Future changes in tropical cyclone activity projected by the new high-resolution MRI-AGCM. J. Climate 25, 32373260 (2012).

29. Landsea, C. W. & Franklin, J. L. Atlantic hurricane database uncertainty and presentation of a new database format. Mon. Weather Rev. 141, 35763592 (2013).

30. Knaff, J. A., Brown, D. P., Courtney, J., Gallina, G. M. & Beven, J. L. An evaluation of Dvorak technique-based tropical cyclone intensity estimates. Weather Forecast. 25, 13621379 (2010).

31. Jarvinen, B. R., Neumann, C. J. & Davis, M. A. S. A Tropical Cyclone Data Tape for the North Atlantic Basin, 1886-1983: Contents, Limitations, and Uses pp 21 (NOAA Technical Memorandum NWS NHC 22, 1984).

32. Chu, J.-H., C. R. Sampson, A. S. L. & Fukada, E. The Joint Typhoon Warning Center tropical cyclone best-tracks, 1945-2000. Report No. NRL/MR/7540-02-16, 22pp (Washington, D. C., USA, 2002).

4 NATURE COMMUNICATIONS | 7:10625 | DOI: 10.1038/ncomms10625 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms10625 ARTICLE

Acknowledgements

The research is supported by Ofce of Naval Research under the research grant of MURI

(N00014-12-1-0911). We thank Dr Shuyi S. Chen from University of Miami for her

inuential suggestion on the possible role of RI on the LMI distribution during 2014

AGU Fall Meeting. Comments and suggestions from Dr John Knaff and the other

anonymous reviewer are appreciated.

Author contributions

The study was led by C.-Y.L. and M.K.T., calculations were carried out and manuscript

was drafted by C.-Y.L. All authors were involved with designing the research, analysing

the results and revising and editing the manuscript. The data were prepared by S.J.C.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/

http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests. Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions

Web End =http://npg.nature.com/

http://npg.nature.com/reprintsandpermissions

Web End =reprintsandpermissions /

How to cite this article: Lee, C-Y. et al. Rapid intensication and the bimodal

distribution of tropical cyclone intensity. Nat. Commun. 7:10625 doi: 10.1038/

ncomms10625 (2016).

This work is licensed under a Creative Commons Attribution 4.0

International License. The images or other third party material in this

article are included in the articles Creative Commons license, unless indicated otherwise

in the credit line; if the material is not included under the Creative Commons license,

users will need to obtain permission from the license holder to reproduce the material.

To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

NATURE COMMUNICATIONS | 7:10625 | DOI: 10.1038/ncomms10625 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5