Abstract
It remains unclear how tropical cyclones (TCs) decay from their ocean lifetime maximum intensity (LMI) to landfall intensity (LI), yet this stage is of fundamental importance governing the socio-economic impact of TCs. Here we show that TCs decay on average by 25% from LMI to LI. A logistic decay model of energy production by ocean enthalpy input and surface dissipation by frictional drag, can physically connect the LMI to LI. The logistic model fits the observed intensity decay as well as an empirically exponential decay does, but with a clear physical foundation. The distance between locations of LMI and TC landfall is found to dominate the variability of the decay from the LMI to LI, whereas environmental conditions are generally less important. A major TC at landfall typically has a very large LMI close to land. The LMI depends on the heating by ocean warming, but the LMI location is also important to future landfall TC intensity changes which are of socio-economic importance.
Similar content being viewed by others
Tropical cyclones (TCs) are one of the major global natural hazards and of great concern to coastal regions. The TC intensity is traditionally recorded as the surface maximum sustained wind. The lifetime maximum intensity (LMI) is the intensity record that has been extensively studied. The LMI is the intensity that is closest to a theoretical upper limit, i.e., the maximum potential intensity, which can be predicted based on environmental thermodynamic conditions1,2. Previous studies3,4,5 showed that the LMI has been increasing in the past decades. The location of LMI is also migrating toward the coasts6 and the poles7, potentially due to the expansion of the tropics8. Regionally, this coastal poleward shift of LMI may change the TC threat to the coasts in the western North Pacific9,10,11 and elsewhere. However, the location of LMI is on average more than 700 km to the coastline6 and this distance is much larger than the typical TC wind radius of gale-force wind of about 200 km12.
It is the landfall intensity (LI) that dominates the destructive potential in coastal regions13,14. More than 80% of normalized TC-related damage in the US is caused by the landfall TCs with major hurricane intensities, that is, LI ≥ 50 m s−1. The intensity from the LMI to landfall have a wide range of behaviours. For example, Hurricane Dorian (2019) made landfall in the Bahamas just after a raid intensification15. However, Typhoon Kong-rey (2018) decayed from a category-5 intensity by about 50% to landfall in South Korea16. To date, there has been no systematic analysis of the intensity decay from the LMI to LI, which we provide here.
Internal and environmental factors have been proposed as the cause of the intensity decay. For example, simulations17 showed a progressive self-weaking of TC intensity after reaching the LMI, which is latter confirmed by observations18. The changes of coastal sea surface temperature19, vertical wind shear20 and entrainment of dry air21,22,23 may also modulate the intensity decay from LMI to landfall. From an energetic perspective a TC system is close to the balance of power generation by heat fluxes and surface frictional dissipation when its intensity reaches the LMI24. Any subsequent unfavourable environmental perturbation but before landfall may break the balance, reduce the power generation, and therefore lead to an intensity decay dominated by surface friction over oceans.
For the intensity reduction after landfall an empirically exponential decay model has been proposed25 and widely used26,27,28. Recently, a physically based algebraic model was also proposed for the decay after landfall29. However, there has been a lack of theoretical model for the decay from LMI to LI prior to landfall. In this study we will propose a simple physical model that connects LMI to LI over oceans. We will demonstrate that the intensity decay from LMI to LI can be understood with a physical logistic model.
Results
The intensity change can be considered as the residual between energy production by ocean enthalpy input and surface dissipation by frictional drag30. Our decay model is an extension of the approach proposed by ref31. We start with their Eq. (12):
where Vm is the maximum wind speed near the surface, Vmpi is the maximum potential intensity, CD is the drag coefficient, E is an efficiency taken to be a normalised inertial stability frequency31, and H is an “effective” depth of vortex. The height parameter H can be understood as the depth over which friction acts to spin down the cyclone, which has been shown to be roughly twice the depth of the boundary layer32. The first term on the r.h.s of Eq. (1) represents surface dissipation due to friction, and the second term corresponds to energy production.
Since we will apply the model to the intensity decay after LMI that is at least 33 m s−1 (i.e., category-1 TCs, see Methods) over a relatively short period of time before landfall (about 1.5 days on average), Eq. (1) can be further simplified to zero order as follows. First, Vmpi during this period may be assumed as stationary since it changes by less than about 5%33. Second, the radius of maximum wind (rm) during this short period may be also assumed as stationary for typical change of less than 15%18. Third, the Coriolis parameter (f) is a small term compared to Vm/rm. With these assumptions, the efficiency, E, is then proportional to Vm [Eq. (9) in Ref31, i.e., \(E={\left[(f+2{V}_{m}/{r}_{m})/(f+2{V}_{mpi}/{r}_{mpi})\right]}^{n}\), where rmpi is the radius of Vmpi and n = 1 as recommended31], and the second term on the r.h.s. in Eq. (1) can be simplified as αVm, where α is a constant for each storm decay. Thus, the simplified version of Eq. (1) can be written in the form of a logistic equation:
where \(\kappa =\frac{{C}_{D}}{H}\) and \(\alpha =\kappa \frac{{r}_{mpi}{V}_{mpi}}{{r}_{m}}\). The parameter κ is defined here as a decay parameter. Since the two terms on the r.h.s. represent frictional dissipation and energy production, we have \(\kappa >0\) and \(\alpha >0\).
Integration from the LMI (Vo) at time 0 to an intensity Vm after time t during decay yields
Equation (3) is the logistic decay model that will be tested against observations.
When \(\alpha t\) is small \({e}^{-\alpha t}\) can be approximated to \(1-\alpha t\), and therefore Eq. (3) can be simplified to \(1/{V}_{m}=1/{V}_{o}+(\kappa -\alpha /{V}_{o})t\). When \(\alpha /{V}_{o}\ll \kappa\) as will be shown in Fig. 2, Eq. (3) can be further simplified as:
Equation (4) represents an algebraic decay consistent with a previous vortex spin down model34,35 when frictional loss dominates. This algebraic decay has been recently validated over land29.
For comparison with the logistic decay [Eq. (3)], we also use a simple exponential decay approximation for the intensity decay from LMI to LI over oceans, which can be written as
where τ can be defined as a decay timescale with a unit of hr36.
Figure 1 shows that the observed intensity decays to the LI on average by about 25% of LMI. The decay time is defined from the last LMI to landfall. Due to the temporally discrete best-track records, the last TC centre record is over land. The observed intensity in Fig. 1 therefore decays more rapidly in the last 10% of decay duration, which reflects an abrupt enhancement of surface friction after landfall. The fit of the logistic decay model, i.e., Eq. (3), shows an excellent ability to quantitatively characterise the observed intensity decay from LMI to landfall with a mean coefficient of determination (r2) of 0.88 (Fig. 1). The same mean r2 is obtained if the exponential decay, i.e., Eq. (5), is used for the fit. The observed mean LI, LMI and decay duration of global TCs are 38 m s−1, 54 m s−1 and 35 h, respectively (Fig. S1). The fitted mean κ, α and τ are 2.2 × 10–7 m−1, 3.1 × 10–6 s−1 and 92 h, respectively (Fig. S1).
With the fitted κ and α, the relative importance of frictional dissipation (\(\kappa {{V}_{m}}^{2}\)) and energy production (\(\alpha {V}_{m}\)) in Eq. (2) is compared in terms of their ratio at LMI [\(\alpha /(\kappa {V}_{o})\)]. Since here we focus on the decay of intensity, by definition, dissipation is larger than production, i.e., all \(\alpha /(\kappa {V}_{o})<1\). For about 60% of decaying TCs, the energy production term is much less than 5% of the frictional dissipation term (the blue bar in Fig. 2a) and thus is negligible. Figure 2b shows that when the production term is important (although still less than the dissipation term), the logistical model performs very well (red lines), which gives confidence in the physical assumptions of the model.
The landfall intensity, LI, can be understood with the logistic decay model by rewriting Eq. (3) as
where T is the duration of the total decay from an initial Vo of LMI to LI, which can be further written as d/c, where d is the distance travelled from LMI to landfall and c is the mean translation speed of the TC during decay. Figure 3 shows that the relative intensity reduction, defined as the ratio of LI to LMI, is strongly and significantly correlated with the distance, d, travelled during the decay (r2 = 0.47, p < 0.05, Fig. 3a), and the decay duration, T (r2 = 0.46, p < 0.05, Fig. 3b). There is also a significant relationship between d and T themselves (r2 = 0.57, p < 0.05). The translation speed, c (i.e., d/T for each TC), however, appears to have no significant impact on the decay (p = 0.46, Fig. 3c). It is important to note that the strong correlation between intensity decay and decay distance d (and therefore decay duration T) is independent of the choice of decay models. By contrast, the fitted parameters (see Methods), i.e., κ and α in the logistic model and τ in the exponential decay, show surprisingly weak correlation with the relative intensity decay from LMI to LI (r2 < 0.05, Fig. 3d–f). For the TCs with major LIs (Fig. S2), the distance d and duration T also show stronger correlations with the decay than the other parameters examined.
We next examine the potential environmental control of the intensity decay parameters. Here three mean environmental conditions (see Methods) during decay are analysed: the potential intensity, vertical wind shear and steering flow approximated by TC translation speed. Figure 4 shows that none of the decay parameters is significantly associated with any of the environmental conditions. Similar results are found (not shown) for the TCs with the distance from LMI to LI above or below the median distance (500 km). The location of TC activity in the western North Pacific can be modulated by the El Niño-Southern Oscillation47. However, the detrended annual mean distance, d, in the northern West Pacific and the Niño3.4 index (August-October) are not significantly correlated (r2 = 0.1, p value = 0.06). The analysis in Figs. 3 and 4 points to the dominant role of the variability of d, rather than the fitted κ, α and τ, or other environmental conditions, as the primary cause of the LI variability.
The potential damage caused by TC landfall increases substantially for major TCs37. We next focus on the properties of these storms with major LIs (≥ 50 m s−1) in terms of the variables in the decay models. As expected, very high LMIs are found in these cases. The global median LMI of these TCs is ranked at about the 80th percentile of the LMIs of all the landfalling TCs (Table 1), which shows the connection between extreme LMIs and major LIs. Globally, relatively small distance (short duration) is important to achieving a weak decay from LMI to LI and hence high LI. The LMI of the TCs with major LIs occurs very close to land (122 km), less than half of the median distance to land for all global landfalling cases (282 km). However, the translation speed, c, is not unusual. All the variables in Table 1 as mentioned above are independent of decay approximations (logistic or exponential). Regarding the fitted parameters in two decay models, τ and α of the TCs with major LIs are close to the medians of all landfall TCs. The energy production term is therefore not unusual for these cases. Interestingly, a small κ appears to contribute to a high LI, which will be discussed in the next section.
Discussion and conclusions
The logistic intensity decay, i.e., Eq. (2), is a modification of an intensification theory31. Two assumptions based on observations18,33 are applied to derive the decay model, that is, the Vmpi and rm are stationary to first approximation from the LMI to LI. Since the derived model captures the intensity decay over oceans as shown in Fig. 1, these two zero-order assumptions give plausible results. A similar logistic intensification model has been used in operational forecast27, but developed from a statistical point of view. This study shows the validity of the logistical model for decay, but from a theoretical perspective.
Previous study38 points to the importance of ventilation of the boundary layer mass by deep convection. Once this ventilation ceases to be sufficient then the spin down mechanism will dominate. From an energetic perspective, TC intensification and decay can be understood as a competition between power production from ocean enthalpy input and surface frictional dissipation31. These two factors reach a balance at the time of LMI24. However, it has been found that even at the time of LMI the surface frictional dissipation rate at the radius of maximum wind is about 25% larger than the local energy production rate39. Thus, an intensity decay can be easily triggered post LMI with minor environmental influence if the excess frictional dissipation under the eyewall cannot be balanced by the energy production outside the eyewall in time, or the boundary layer ventilation by convection becomes ineffective. Our analysis shows that the decay from LMI to landfall can be mostly captured by frictional dissipation in about 60% of the TCs without the need to invoke a strong role for energy production by enthalpy fluxes.
The simplified algebraic decay [Eq. (4)] also predicts that the LI can be understood with two terms: an upper limit, 1/LMI, and the decay, κT. In line with this prediction, Table 1 suggests that TCs with major LI usually have relatively small κ. The median κT of the TCs with major LI is ranked at the 18th percentile of that for all landfall TCs, which is a more extreme percentile compared to the rankings in Table 1. It is nevertheless surprising that in general, the decay parameter, κ, does not have a strong predictive relationship with the intensity reduction when all landfall TCs are included in the statistics. We note that there may be some cancellation effects as κ depends on both the effective vortex depth H and the drag coefficient CD which may vary in the same direction when TC is relatively weak. However, there are no direct observations of either term. There are also no systematically significant relationships between κ and either environmental thermodynamic condition or vertical wind shear. Although κ appears on average to play a minor role in the climatology of the decay, we still expect environmental conditions to be important for individual cases40 and some subgroups of TCs, for example, those experiencing a rapid weakening41. We show that a small κ is a necessary condition for major TCs at landfall.
The decay parameter, κ, can be quantitatively broken into physically relevant components: the ratio of surface drag coefficient to an effective depth of vortex (i.e., κ = CD/H as defined in section “Results”). The mean fitted κ of 2.2 × 10–7 m−1 is also consistent with typical values of CD = 1.0 × 10–3 over oceans42 and an H = 4.5 km. Recall that H is defined as the “effective” vortex depth and is thus also the depth over which friction acts to spin down the cyclone. This estimation of 4.5 km is close to a boundary layer depth scale of 5 km in the eyewall as previously hypothesised32. We also note that it is approximately the height below which the buoyancy and Richardson number decrease rapidly in model simulations43.
We may also define a half-life of intensity decay with the simplified algebraic model as 1/(κLMI). This definition of half-life shows a self-regulating nature of intensity decay, that is, a more intense TC decays faster, whereas an exponential decay depends solely on the decay timescale τ36. Previous study28 found that the exponential decay approximation does not perform well when the initial decay intensity is high, which supports our finding here that the intensity decay should depend not only on the decay parameter (i.e., κ or τ), but also the initial decay intensity. This dependence on the initial decay intensity can lead to dampening of any positive LMI trend at landfall through an increasingly more rapid decay. For example, according to Eq. (4) the intensity difference of two TCs with LMIs of 70 and 50 m s−1 converges at landfall to a difference of only 5 m s−1 for T = 48 h and κ = 1 × 10–7 m−1. This may provide an explanation for the less significant LI trend44 compared to the LMI trend over oceans4,45,46, and fundamentally distinguishes the proposed physically based logistic model from exponential decay models25,26.
Major TCs at landfall are characterised by both high LMI and close proximity to land. The proposed decay model usefully points to the stringent and rare conditions required for a TC to make landfall as a major TC. Most climate models (e.g., Ref47) already struggle to generate the most intense storms in terms of LMI. This importance of distance creates another challenge to simulate plausible LIs given the model typical horizontal resolution of more than 100 km. Our findings may also have implications for weather forecasting and stochastic catastrophe models used in the insurance sector. Forecasting the most damaging storms with extreme intensities at landfall requires forecasts of both the LMI and its location close to land. For the stochastic models the TC evolution from genesis to landfall may be much less important than the ability to modelling the LMI, LMI location and subsequent decay. There has been much discussion in the literature and in the public domain on the role of global warming in changing TC intensity. Thermodynamic arguments point to the role of a warming ocean as the source of heat to drive intensity changes. Here we show that to understand the changes at landfall this is only part of the story. The changes in distance from LMI to landfall are also of critical importance to understand the past variability and predict future changes in tropical cyclone intensities at landfall.
Methods
Data
We consider landfalling TCs with an LMI of at least hurricane-force wind (LMI ≥ 33 m s−1) for the period 1982–2019. In this period we have the highest confidence in data quality and the completeness of global TC observations48. We take the TC best-track data from the International Best Track Archive for Climate Stewardship (IBTrACS) v04r0049, with the original data sources from the National Hurricane Center and the Joint Typhoon Warming Center. The TC positions and the other measures in the IBTrACS are interpolated from 6-h to 3-h intervals with splines and linear interpolations, respectively. In this study we use the IBTrACS best-track data at 3-h intervals, i.e., 00, 03, 06, 09, 12, 15, 18 and 21 Universal Time Coordinates. Only TCs within the 40oN/S latitude band are considered to reduce extra-tropical impacts.
Landfall criterion and decay period
The distance to the nearest land is available in the IBTrACS for each TC location. The smallest landmass used for the distance calculation is 1400 km2, equivalent to the area of Kauai, Hawaii. This distance is positive when the TC centre is over the ocean and drops to zero after landfall. We define a landfall at the time when the distance to land is reduced to zero from a positive distance. Some TCs made landfall more than once. In these cases, only the landfall with the highest landfall intensity is considered. A TC may also reach the same LMI more than once during its lifecycle. We define the decay period from the time of the last LMI to landfall. Seven percent of global landfalls happen at LMI, in which case the LMI and LI are equal and there is no intensity decay.
Parameter estimation
The parameters in the decay models [i.e., κ and α in Eq. (3) and τ in Eq. (5)] are estimated with the least-squares minimization fit to the entire intensity evolution from LMI to LI for each TC. There are on average 13 data points including LI and LMI for all the landfalling TCs with LMI ≥ 33 m s−1, and 6 data points for the TCs with major LI (LI ≥ 50 m s−1). For each fit we use the observed LMI for Vo and decay duration for t.
Environmental conditions
The hourly atmospheric reanalysis product, ERA550, is used to calculate the actual TC-related vertical wind shear and maximum potential intensity (MPI) for each TC during the period 1982–2019 rather than monthly means. For each TC the variables are calculated three-hourly along its track during decay, and we then take the mean to represent this environmental condition for the TC. Vertical wind shear is calculated as the magnitude of wind vector difference between 200- and 850-hPa pressure levels in the 2-to-8-degree-latitude annulus around TC centres. The MPI1 is calculated along the track three days before a TC arrives.
Data availability
Tropical cyclone best track data can be downloaded from the National Centers for Environmental Information website (https://www.ncei.noaa.gov/data/international-best-track-archive-for-climate-stewardship-ibtracs/v04r00/access/csv/ibtracs.ALL.list.v04r00.csv). The ERA5 reanalysis data is available at the European Centre for Medium-Range Weather Forecasts (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=form).
Code availability
All codes used to analyse and plot the data are available at https://data.hpc.imperial.ac.uk/resolve/?doi=8721&access= .
Change history
04 August 2022
A Correction to this paper has been published: https://doi.org/10.1038/s41598-022-17834-4
References
Emanuel, K. Sensitivity of tropical cyclones to surface exchange coefficients and a revised steady-state model incorporating eye dynamics. J. Atmos. Sci. 52, 3969–3976 (1995).
Holland, G. J. The maximum potential intensity of tropical cyclones. J. Atmos. Sci. 54, 2519–2541 (1997).
Emanuel, K. Increasing destructiveness of tropical cyclones over the past 30 years. Nature 436, 686–688 (2005).
Wang, S., Rashid, T., Throp, H. & Toumi, R. A shortening of the life cycle of major tropical cyclones. Geophys. Res. Lett. 47, 1–8 (2020).
Wang, S. & Toumi, R. More tropical cyclones are striking coasts with major intensities at landfall. Sci. Rep. Under review. (2022).
Wang, S. & Toumi, R. Recent migration of tropical cyclones toward coasts. Science (80-) 371, 514–517 (2021).
Kossin, J. P., Emanuel, K. A. & Vecchi, G. A. The poleward migration of the location of tropical cyclone maximum intensity. Nature 509, 349–352 (2014).
Sharmila, S. & Walsh, K. J. E. Recent poleward shift of tropical cyclone formation linked to Hadley cell expansion. Nat. Clim. Chang. 8, 730–736 (2018).
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).
Kossin, J. P., Emanuel, K. A. & Camargo, S. J. Past and projected changes in western north pacific tropical cyclone exposure. J. Clim. 29, 5725–5739 (2016).
Wang, S. et al. Is the tropical cyclone surge in Shanghai more sensitive to landfall location or intensity change?. Atmos. Sci. Lett. 22, 1–8 (2021).
Chan, K. T. F. & Chan, J. C. L. Global climatology of tropical cyclone size as inferred from QuikSCAT data. Int. J. Climatol. 35, 4843–4848 (2015).
Wang, S. & Toumi, R. On the relationship between hurricane cost and the integrated wind profile. Environ. Res. Lett. 11, 114005 (2016).
Klotzbach, P. P. J., Bowen, S. G., PielKe, R. G. R. & Bell, M. Continental U.S. hurricane landfall frequency and associated damage: Observations and future risks. Bull. Am. Meteorol. Soc. 99, 1359–1376 (2018).
Ryglicki, D. R., Velden, C. S., Reasor, P. D., Hodyss, D. & Doyle, J. D. Observations of atypical rapid intensification characteristics in hurricane dorian (2019). Mon. Weather Rev. 149, 2131–2150 (2021).
Wada, A. The impacts of a cold eddy induced by Typhoon Trami (2018) on the intensity forecast of Typhoon Kong-Rey (2018). CAS/JSC WGNE Res Act Atm Ocean Model 49, 7–9 (2019).
Kilroy, G., Smith, R. K. & Montgomery, M. T. Why do model tropical cyclones grow progressively in size and decay in intensity after reaching maturity?. J. Atmos. Sci. 73, 467–503 (2016).
Wang, S. & Toumi, R. A historical analysis of the mature stage of tropical cyclones. Int. J. Climatol. 38, 2490–2505 (2018).
Mei, W. & Xie, S. Intensification of landfalling typhoons over the northwest Pacific since the late 1970s. Nat. Geosci. 9, 753–757 (2016).
Kossin, J. P. Hurricane intensification along United States coast suppressed during active hurricane periods. Nature 541, 390–393 (2017).
Tang, B. & Emanuel, K. Midlevel ventilation’s constraint on tropical cyclone intensity. J. Atmos. Sci. 67, 1817–1830 (2010).
Corsaro, C. M. & Toumi, R. A self-weakening mechanism for tropical cyclones. Q. J. R. Meteorol. Soc. 143, 2585–2599 (2017).
Wang, S. & Toumi, R. Impact of dry midlevel air on the tropical cyclone outer circulation. J. Atmos. Sci. 76, 1809–1826 (2019).
Emanuel, K. Some aspects of hurricane inner-core dynamics and energetics. J. Atmos. Sci. 54, 1014–1026 (1997).
Kaplan, J. & DeMaria, M. A simple empirical model for predicting the decay of tropical cyclone winds after landfall. J. Appl. Meteorol. 34, 2499–2512 (1995).
Kaplan, J. & DeMaria, M. On the decay of tropical cyclone winds after landfall in the New England Area. J. Appl. Meteorol. 40, 280–286 (2001).
DeMaria, M. A simplified dynamical system for tropical cyclone intensity prediction. Mon. Weather Rev. 137, 68–82 (2009).
Jing, R. & Lin, N. Tropical cyclone intensity evolution modeled as a dependent hidden Markov process. J. Clim. 32, 7837–7855 (2019).
Phillipson, L. M. & Toumi, R. A physical interpretation of recent tropical cyclone post-landfall decay. Geophys. Res. Lett. 48, 1–9 (2021).
Emanuel, K. The theory of hurricanes. Ann Rev Fluid Mech 23, 179–196 (1991).
Wang, Y., Li, Y., Xu, J., Tan, Z.-M. & Lin, Y. The intensity-dependence of tropical cyclone intensification rate in a simplified energetically based dynamical system model. J. Atmos. Sci. https://doi.org/10.1175/jas-d-20-0393.1 (2021).
Emanuel, K. Self-stratification of tropical cyclone outflow. Part II: implications for storm intensification. J. Atmos. Sci. 69, 988–996 (2012).
Emanuel, K. A statistical analysis of tropical cyclone intensity. Mon. Weather Rev. 128, 1139–1152 (2000).
Eliassen, A. & Lystad, M. The Ekman layer of a circular vortex. A numerical and theoretical study. Geophys. Norvegica 31, 1–16 (1977).
Montgomery, M. T., Snell, H. D. & Yang, Z. Axisymmetric spindown dynamics of hurricane-like vortices. J. Atmos. Sci. 58, 421–435 (2001).
Li, L. & Chakraborty, P. Slower decay of landfalling hurricanes in a warming world. Nature 587, 230–234 (2020).
Pielke, R. A. & Landsea, C. W. Normalized hurricane damages in the United States: 1925–95. Weather Forecast. 13, 621–631 (1998).
Smith, R. K., Kilroy, G. & Montgomery, M. T. Tropical cyclone life cycle in a three-dimensional numerical simulation. Q. J. R. Meteorol. Soc. 147, 3373–3393 (2021).
Wang, Y. & Xu, J. Energy production, frictional dissipation, and maximum intensity of a numerically simulated tropical cyclone. J. Atmos. Sci. 67, 97–116 (2010).
Colomb, A., Kriat, T. & Leroux, M. D. On the rapid weakening of very intense tropical Cyclone Hellen (2014). Mon. Weather Rev. 147, 2717–2737 (2019).
Song, J., Klotzbach, P. J., Duan, Y. & Guo, H. Recent increase in tropical cyclone weakening rates over the western north pacific. Geophys. Res. Lett. 47, 1–11 (2020).
Holton, J. R. & Hakim, G. J. An introduction to dynamic meteorology: Fifth edition. An Introd. to Dyn. Meteorol. Fifth Ed. 9780123848, 1–532 (2012).
Zhu, P., Hazelton, A., Zhang, Z., Marks, F. D. & Tallapragada, V. The Role of Eyewall Turbulent Transport in the Pathway to Intensification of Tropical Cyclones. J. Geophys. Res. Atmos. 126, (2021).
Weinkle, J., Maue, R. & Pielke, R. Historical global tropical cyclone landfalls. J. Clim. 25, 4729–4735 (2012).
Kossin, J. P., Olander, T. L. & Knapp, K. R. Trend analysis with a new global record of tropical cyclone intensity. J. Clim. 26, 9960–9976 (2013).
Kossin, J. P., Knapp, K. R., Olander, T. L. & Velden, C. S. Global increase in major tropical cyclone exceedance probability over the past four decades. Proc. Natl. Acad. Sci. U. S. A. 117, 11975–11980 (2020).
Roberts, M. J. et al. Impact of model resolution on tropical cyclone simulation using the HighResMIP-PRIMAVERA multimodel ensemble. J. Clim. 33, 2557–2583 (2020).
Knapp, K. R. & Kruk, M. C. Quantifying interagency differences in tropical cyclone best-track wind speed estimates. Mon. Weather Rev. 138, 1459–1473 (2010).
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. Meteorol. Soc. 91, 363–376 (2010).
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146, 1999–2049 (2020).
Acknowledgements
We acknowledge that this research is jointly supported by the Natural Environment Research Council/UKRI (NE/V017756/1), and the UK-China Research and Innovation Partnership Fund through the Met Office Climate Science for Service Partnership (CSSP) China as part of the Newton Fund.
Author information
Authors and Affiliations
Contributions
S.W. and R.T. conceived the study. S.W. performed the analysis. Both the authors discussed the results and jointly contributed to writing the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this Article was revised: The original version of this Article contained an error in the Results section below Equation 2, where the Greek letter κ was omitted in the explanation of variable α. Full information regarding the corrections made can be found in the correction for this Article.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, S., Toumi, R. On the intensity decay of tropical cyclones before landfall. Sci Rep 12, 3288 (2022). https://doi.org/10.1038/s41598-022-07310-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-07310-4
This article is cited by
-
Sea surface temperature warming to inhibit mitigation of tropical cyclone destructiveness over East Asia in El Niño
npj Climate and Atmospheric Science (2024)
-
The Imperial College Storm Model (IRIS) Dataset
Scientific Data (2024)
-
Temporal and spatial distribution of extreme rainfall from tropical storms in the Gulf of Mexico from 1979 to 2021
Stochastic Environmental Research and Risk Assessment (2024)
-
More tropical cyclones are striking coasts with major intensities at landfall
Scientific Reports (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.