Sensitivity of Tropical Cyclone Intensity Variability to Different Stochastic Parameterization Methods

This study examines the dependence of tropical cyclone (TC) intensity forecast errors on track forecast errors in the Coupled Ocean–Atmosphere Mesoscale Prediction System for Tropical Cyclones (COAMPS-TC) model. Using real-time forecasts and retrospective experiments during 2015–18, verification of TC intensity errors conditioned on different 5-day track error thresholds shows that reducing the 5-day track errors by 50%–70% can help reduce the absolute intensity errors by 18%–20% in the 2018 version of the COAMPS-TC model. Such impacts of track errors on the TC intensity errors are most persistent at 4–5-day lead times in all three major ocean basins, indicating a significant control of global models on the forecast skill of the COAMPS-TC model. It is of interest to find, however, that lowering the 5-day track errors below 80 n mi (1 n mi = 1.852 km) does not reduce TC absolute intensity errors further. Instead, the 4–5-day intensity errors appear to be saturated at around 10–12 kt (1 kt ≈ 0.51 m s−1) for cases with small track errors, thus suggesting the existence of some inherent intensity errors in regional models. Additional idealized simulations under a perfect model scenario reveal that the COAMPS-TC model possesses an intrinsic intensity variation at the TC mature stage in the range of 4–5 kt, regardless of the large-scale environment. Such intrinsic intensity variability in the COAMPS-TC model highlights the importance of potential chaotic TC dynamics, rather than model deficiencies, in determining TC intensity errors at 4–5-day lead times. These results suggest a fundamental limit in the improvement of TC intensity forecasts by numerical models that one should consider in future model development and evaluation.

This study examines the potential limit in predicting tropical cyclone (TC) intensity under idealized conditions. Using the phase-space reconstruction method for TC intensity time series obtained from the CM1 idealized simulations, it is found that CM1 axisymmetric dynamics contain low-dimensional chaos at the maximum intensity equilibrium. Examination of several attractor invariants including the largest Lyapunov exponent, the Sugihara–May correlation, and the correlation dimension captures a consistent range of the chaotic attractor dimension between 4 and 5 for TC intensity at the maximum intensity equilibrium. In addition, the intensity error doubling time estimated from the largest Lyapunov exponent is roughly 1–3 h, which accords with the decay time obtained from the Sugihara–May correlation. Furthermore, the findings in this study reveal a relatively short TC intensity predictability limit for CM1, which is ∼3–9 h based on the maximum tangential wind but noticeably longer for the minimum central pressure (∼12–18 h) after reaching the mature stage. So long as the traditional metrics for TC intensity such as the maximum surface wind or the minimum central pressure is used for intensity forecast, our results support that TC intensity forecast errors will not be reduced indefinitely in any model, even in the absence of all model and observational errors. As such, the future improvement of TC intensity forecast should be based on different metrics beyond the absolute intensity errors that are currently used in real-time intensity verification. Significance Statement Using the phase-space reconstruction method for tropical cyclone (TC) intensity time series obtained from idealized axisymmetric simulations, we show that TC axisymmetric dynamics in CM1 possesses low-dimensional chaos at the maximum intensity equilibrium. This low-dimensional dynamics explains the long tradition of representing TC intensity by a few measures as in the current practice. The chaotic property of CM1 axisymmetric dynamics also suggests a relatively short predictability range for TC intensity at the maximum intensity equilibrium. The potential existence of low-dimensional chaos for TC intensity in CM1 idealized simulations as found in this study supports the use of different intensity verification metrics beyond the traditional absolute intensity errors currently used in operational model evaluation.

Climatological characteristics of simulated intense tropical cyclones (TCs) in the western North Pacific were explored with a 20-km-mesh atmospheric general circulation model (AGCM20) and a 5-km-mesh regional atmospheric nonhydrostatic model (ANHM5). From the AGCM20 climate runs, 34 intense TCs with a minimum central pressure (MCP) less than or equal to 900 hPa were sampled. Downscaling experiments were conducted with the ANHM5 for each intense TC simulated by the AGCM20. Only 23 developed into TCs with MCP ≤ 900 hPa. Most of the best-track TCs with an MCP ≤ 900 hPa underwent rapid intensification (RI) and attained maximum intensities south of 25°N. The AGCM20 simulated a similar number of intense TCs as the best-track datasets. However, the intense AGCM20 TCs tended to intensify longer and more gradually; only half of them underwent RI. The prolonged gradual intensification resulted in significant northward shifts of the location of maximum intensity compared with the location derived from two best-track datasets. The inner-core structure of AGCM20 TCs exhibited weak and shallow eyewall updrafts with maxima below an altitude of 6 km, while downscaling experiments revealed that most of the intense ANHM5 TCs underwent RI with deep and intense eyewall updrafts and attained their maximum intensity at lower latitudes. The altitudes of updraft maxima simulated by the AGCM20 descended rapidly during the phase of greatest intensification as midlevel warming markedly developed. The change in major processes responsible for precipitation in AGCM20 TCs before and after maximum intensification suggests close relationships between the large-scale cloud scheme and midlevel warming and prolonged gradual intensification.

In this study, the relationship between tropical cyclone (TC) intensity and rain rate over the ocean is investigated using a full-physics numerical model (WRF) and a physics-based TC rainfall model (TCR). TC intensity is found to be nearly linearly correlated with the average rain rate in the inner core [∼0.97 (mm h−1 m−2)/(m s−1)], while the correlation is weak at outer radii. This difference is induced because TC intensity is significantly correlated with both the vertical velocity and specific humidity in the inner core but is not significantly correlated with the vertical velocity in the outer radii. Further investigation shows that the intensity–rain-rate relationship at the outer radii is influenced by the TC evolution stage. The rain rate for the outer radii is positively correlated with TC intensity for nondecaying TCs, while this correlation is reduced for decaying TCs due to systematic downdrafts in the outer radii. In the context of climate change, the sensitivity of the TC rain rate to sea surface temperature (SST) is found to be +9% per 1 K increase of SST, roughly the product of the sensitivity of TC intensity to SST (+3%) and the Clausius-Clapeyron scaling (+7%). Coupled with synthetic storms, evolution of the TC rain rate over the twenty-first century under the SSP5-8.5 scenario is projected by the TCR (calibrated with the WRF simulations). The annual increase rates of averaged TC rain rate are 0.17% and 0.20% for the inner core and outer radii, respectively, larger than the annual increase rate of TC intensity (0.046%) but comparable to that of cube of intensity (0.18%).

Tropical cyclone (TC) intensity is controlled by both environmental factors and its internal structure. This paper investigates the influence of inner core features derived from satellite data on TC intensity changes and explores effective methods to improve TC intensity forecast accuracy. We construct models for forecasting TC intensity over the next 5 days using the light gradient boosting machine (LGBM) algorithm in the northwestern Pacific. The models integrate climatological and persistence predictors, large-scale environmental variables at the initial time, inner core features extracted from satellite data, and variables related to the forecasted TC location. The models are trained on TC data from 2005 to 2019, tested using best track data from 2020 to 2022 and real-time data in 2023. Results show that the inner core features derived from satellite data are helpful for intensity forecasts within 36 h, and variables related to the forecasted location are helpful for improving TC intensity forecast accuracy in most cases, especially for the forecasts beyond 48 h. Compared with the National Centers for Environmental Prediction (NCEP) and European Centre for Medium-Range Weather Forecasts (ECMWF) models, our models demonstrate 4%–48% and 19%–70% decreases in mean absolute errors (MAEs) in real-time intensity forecasts. The most important predictors identified include potential future intensity change, intensity changes during the previous 12 h, vertical wind shear, and sea–land ratio. The inner core feature extracted from satellite data is also among the top 10 important variables for short-time intensity forecasts. Significance Statement This study advances tropical cyclone (TC) intensity forecasting by integrating climatological and persistence predictors, large-scale environmental variables at the initial time, satellite-derived inner core features, and variables related to the forecasted TC location into the light gradient boosting machine (LGBM) models. Our models significantly improve forecast accuracy, reducing mean absolute errors by 4%–48% compared with the National Centers for Environmental Prediction (NCEP) and 19%–70% compared with the European Centre for Medium-Range Weather Forecasts (ECMWF). Results denote that the inner core features improve the intensity prediction in the short time, and variables related to the forecasted TC location improve the intensity forecasts in most cases, especially for the forecasts beyond 48 h. This study provides valuable insights for enhancing operational TC intensity forecasting.

The paper uses observational data from 1950 to 2014 to investigate rapid intensification (RI) variability of tropical cyclones (TCs) in the North Atlantic and its relationships with large-scale climate variations. RI is defined as a TC intensity increase of at least 15.4 m/s (30 knots) in 24 h. The seasonal RI distribution follows the seasonal TC distribution, with the highest number in September. Although an RI event can occur anywhere over the tropical North Atlantic (TNA), there are three regions of maximum RI occurrence: (1) the western TNA of 12°N–18°N and 60°W–45°W, (2) the Gulf of Mexico and the western Caribbean Sea, and (3) the open ocean southeast and east of Florida. RI events also show a minimum value in the eastern Caribbean Sea north of South America—a place called a hurricane graveyard due to atmospheric divergence and subsidence. On longer time scales, RI displays both interannual and multidecadal variability, but RI does not show a long-term trend due to global warming. The top three climate indices showing high correlations with RI are the June-November ENSO and Atlantic warm pool indices, and the January-March North Atlantic oscillation index. It is found that variabilities of vertical wind shear and TC heat potential are important for TC RI in the hurricane main development region, whereas relative humidity at 500 hPa is the main factor responsible for TC RI in the eastern TNA. However, the large-scale oceanic and atmospheric variables analyzed in this study do not show an important role in TC RI in the Gulf of Mexico and the open ocean southeast and east of Florida. This suggests that other factors such as small-scale changes of oceanic and atmospheric variables or TC internal processes may be responsible for TC RI in these two regions. Additionally, the analyses indicate that large-scale atmospheric and oceanic variables are not critical to TC genesis and formation; however, once a tropical depression forms, large-scale climate variations play a role in TC intensification.

This study examines the variation in tropical cyclone (TC) intensity for different phases of the El Niño–Southern Oscillation (ENSO) phenomenon in the Fiji, Samoa, and Tonga (FST) region. The variation in TC intensity is inferred from the accumulated cyclone energy (ACE), which is constructed from the 6-hourly Joint Typhoon Warning Center best-track data for the period 1985–2006. Overall, results suggest that ACE in the FST region is considerably influenced by the ENSO signal. A substantial contribution to this ENSO signal in ACE comes from the region equatorward of 15°S where TC numbers, lifetime, and intensity all play a significant role. However, the ACE–ENSO relationship weakens substantially poleward of 15°S where large-scale environmental variables affecting TC intensity are found to be less favorable during El Niño years than during La Niña years; in the region equatorward of 15°S, the reverse is true. Therefore, TCs entering this region poleward of 15°S are able to sustain their intensity for a longer period of time during La Niña years as opposed to TCs entering the region during El Niño years, when they decay more rapidly.

Using primitive equations with adiabatic, non‐friction non‐environmental steering, the dynamical equilibrium features of tropical cyclones (TC) are investigated. The parameter, dI/dt, was obtained and its dynamical balance formulation describing TC intensity change. The results of the we discussed the characteristics of TC motion and intensity change in barotropic and baroclinic model showed that the stationary motion of TC possess “equivalent‐geostrophic equilibrium” and isolines of potential vorticity and E(E = φ + (u2 + v2)/2) coincide with streamlines. The TC intensity will be weakened by “ventilation flow” in “beta (β) gyres” and “shear gyres” from the positive (negative) wind vertical shear at high‐ (low‐) level. The positive (negative) vorticity vertical convection of high‐ (low‐) level strengthens TC intensity. Generally, the TC center has the tendency to drift poleward and there is a gradient (ascendent) direction of vorticity vertical convection at high‐ (low‐) level. The TC motion is connected with the asymmetry structure of vertical shear and vertical motion, the “γ effect” (γ = ∂β/∂y) that is one of the dynamical mechanisms of “ventilation flow” steering TC center motion.

A tropical cyclone (TC) vortex is an immense, coherent, organized-convective system. Beneath this large-scale organization, is a litany of azimuthally asymmetric convective motions that exist on a spectrum of scales. These asymmetries are especially dominant during periods when the vortex undergoes critical transitions in its intensity and structure. However, the precise nature of influence of the organization of asymmetries on TC intensity change remains an enigma. The inherent difficulty in predicting their behavior is because asymmetries may arise due to different external or intrinsic sources and occur at different spatial and temporal scales while several complex mechanisms act near-simultaneously to dictate their evolution in time. As a result, multiple pathways are possible for a TC vortex that is influenced by these asymmetries. Our preliminary investigations using numerical models made it apparent that there wasn't a single, unifying way to address this problem. In this thesis, I outline multiple novel techniques of diagnosing and predicting which of the many pathways are likely for a TC vortex that is influenced by azimuthal asymmetries. First, using three-dimensional numerical simulations of a pair of sheared and non-sheared vortices, I demonstrate the diagnostic potential of the juxtaposition in the azimuthal phasing of: (i) the asymmetrically distributed vertical eddy flux of moist-entropy across the top of the boundary layer, and the radial eddy flux of moist-entropy within the boundary layer; and (ii) eddy relative vorticity, eddy moist-entropy, and vertical velocity throughout the depth of the vortex. Second, I introduce an energetics-based diagnostic framework that computes the energy transactions occurring at asymmetries across various length-scales in the wavenumber domain. By applying it to select cases, this thesis uncovers the relative importance of all the energy pathways that support or disrupt the growth of asymmetries within the vortex. Contrary to the traditional explanations of convective aggregation/disaggregation and axi/asymmetrization through barotropic mean-eddy transactions, my thesis reveals that the growth or disruption of asymmetries are predominantly due to (i) the baroclinic conversion from available potential to kinetic energy at individual scales of asymmetries and (ii) the transactions of kinetic energy across asymmetries of different length scales. Finally, this thesis introduces two further diagnostic frameworks targeted at tackling the problem of real-time forecasting of TC rapid intensity changes. The first is an empirical framework which examines symmetric and asymmetric convection and other state variables within the vortex, and in the environment across a suite of TCs and identifies a set of `important' variables that are significantly different during time periods that precede a rapid intensification as opposed to a rapid weakening. My framework then ranks the variables identified based on how significantly they influence a rapid intensity change in a TC and the amplification factor of any associated variability. We recommend that future observational, and consequent TC modeling and data assimilation efforts prioritize the highest ranked variables identified here. The second is a stochastic model wherein a scale-specific stochastic term is added to the equations describing the energy transactions within the TC vortex. By simulating a stochastic forcing that may arise from any scale, I compute the probability of the vortex transitioning into a rapidly intensifying or a rapidly weakening configuration across an ensemble of scenarios. In summary, this thesis introduces and applies a variety of diagnostic techniques that help determine the impact of azimuthal asymmetries on TC intensity evolution.

Several key issues in the simple time-dependent theories of tropical cyclone (TC) intensification developed in recent years remain, including the lack of a closure for the pressure dependence of saturation enthalpy at sea surface temperature (SST) under the eyewall and the definition of environmental conditions, such as the boundary layer enthalpy in TC environment and the TC outflow-layer temperature. In this study, some refinements to the most recent time-dependent theory of TC intensification have been accomplished to resolve those issues. The first is the construction of a functional relationship between the surface pressure under the eyewall and the TC intensity, which is derived using the cyclostrophic wind balance and calibrated using full-physics axisymmetric model simulations. The second is the definition of TC environment that explicitly includes the air–sea temperature difference. The third is the TC outflow-layer temperature parameterized as a linear function of SST based on global reanalysis data. With these refinements, the updated time-dependent theory becomes self-contained and can give both the intensity-dependent TC intensification rate (IR) and the maximum potential intensity (MPI) under given environmental thermodynamic conditions. It is shown that the pressure dependence of saturation enthalpy at SST can lead to an increase in the TC MPI and IR by about half of that induced by dissipative heating due to surface friction. Results also show that both MPI and IR increase with increasing SST, surface enthalpy exchange coefficient, environmental air–sea temperature difference, and decreasing environmental boundary layer relative humidity, but the maximum IR is insensitive to surface drag coefficient. Significance Statement A new advancement in the recent decade is the development of simple time-dependent theories of tropical cyclone (TC) intensification, which can provide quantitative understanding of TC intensity change. However, several key issues in these simple time-dependent theories remain, including the lack of a closure for the pressure dependence of saturation enthalpy at sea surface temperature under the eyewall and the definition of environmental conditions. These are resolved in this study with several refinements, which make the most recent time-dependent theory of TC intensification self-contained and practical.

The practical predictability of tropical cyclone (TC) intensity in terms of mean absolute forecast error with respect to different conditions at forecast initialization was explored through convection-permitting hindcasts of all Atlantic storms during the 2008–12 hurricane seasons using the Weather Research and Forecasting (WRF) Model. Averaged over a total of 2190 simulations, the day 1–5 performance of these WRF hindcasts was comparable to two operational regional-scale hurricane prediction models used by the National Hurricane Center (NHC) but was slightly inferior to the NHC official forecasts. It was found that the prediction accuracy of TC intensity, both at the initialization time and the targeted forecast hours, was strongly correlated with the TC intensity. On average, for both the WRF hindcasts and the NHC official forecasts, stronger intensities and larger intensity variations led to larger forecast errors. A number of synoptic-scale environmental parameters, such as vertical wind shear, sea surface temperature (SST), and the underlying surface condition (land vs sea), affected the intensity forecast errors of TCs, in part due to their influence on intensity changes, while other thermodynamic environmental parameters, such as moisture and instability, had relatively minor effects. The accuracy of the intensity prediction was also found to be sensitive to the translation speed of the TCs. A moderate TC translation speed of 11–15 knots (kt; 1 kt = 0.51 m s−1) corresponded to the largest intensity errors during forecast lead times less than 60 h, while the slowest translation speed (<7 kt) was associated with the largest errors after the 60-h forecast lead time.

A tropical cyclone intensity prediction model using conditional generative adversarial network

The relationship between sea surface temperature (SST) and tropical cyclone (TC) intensity change exhibits a strong dependence on the current TC intensity. Using western North Pacific TC observations from 1982 to 2018, a threshold SST (TSST) is identified as the SST required to maintain TC intensity. TSST increases with TC intensity, with TCs intensifying and weakening when SST is higher and lower than TSST, respectively. Across the dataset, mean TC intensity change is proportional to the difference between SST and TSST. This study also formulates an equation to quantify TC intensity change using SST and current TC intensity, which replicates 99.46% of the mean observed TC intensity changes. This equation could serve as an alternative to the linear regression‐based relationship between SST and TC intensity change that is widely used in statistical‐dynamical intensity models, thereby improving intensity forecasts.

Recent studies have shown that assimilating enhanced satellite-derived atmospheric motion vectors (AMVs) has improved mesoscale forecast of tropical cyclones (TC) track and intensity. The authors conduct data-denial experiments to understand where the TC analyses and forecasts benefit the most from the enhanced AMV information using an ensemble Kalman filter and the Weather Research and Forecasting Model. The Cooperative Institute for Meteorological Satellite Studies at the University of Wisconsin provides enhanced AMV datasets with higher density and temporal resolution using shorter-interval image triplets for the duration of Typhoon Sinlaku and Hurricane Ike (both 2008). These AMV datasets are then spatially and vertically subsetted to create six parallel cycled assimilation-forecast experiments for each TC: all AMVs; AMVs withheld between 100 and 350 hPa (upper layer), between 350 and 700 hPa (middle layer), and between 700 and 999 hPa (lower layer); and only AMVs within (interior) and outside (exterior) 1000-km radius of the TC center. All AMV subsets are found to be useful in some capacity. The interior and upper-layer AMVs are particularly crucial for improving initial TC position, intensity, and the three-dimensional wind structure along with their forecasts. Compared with denying interior or exterior AMVs, withholding AMVs in different tropospheric layers had less impact on TC intensity and size forecasts. The ensemble forecast is less certain (larger spread) in providing accurate TC track, intensity, and size when upper-layer AMVs or interior AMVs are withheld. This information could be useful to potential targeting scenarios, such as activating and focusing satellite rapid-scan operations, and decisions regarding observing system assessments and deployments.

Impacts of Climate Change on Tropical Cyclones in the Western North Pacific Basin. Part II: Late Twenty-First Century Projections