skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Sensitivity of large eddy simulations of tropical cyclone to sub-grid scale mixing parameterization.
The surface wind structure and vertical turbulent transport processes in the eyewall of hurricane Isabel (2003) are investigated using six large-eddy simulations (LESs) with different horizontal grid spacing and three-dimensional (3D) sub-grid scale (SGS) turbulent mixing models and a convection permitting simulation that uses a coarser grid spacing and one-dimensional vertical turbulent mixing scheme. The mean radius-height distribution of storm tangential wind and radial flow, vertical velocity structure, and turbulent kinetic energy and momentum fluxes in the boundary layer generated by LESs are consistent with those derived from historical dropsonde composites, Doppler radar, and aircraft measurements. Unlike the convection permitting simulation that produces storm wind fields lacking small-scale disturbances, all LESs are able to produce sub-kilometer and kilometer scale eddy circulations in the eyewall. The inter-LES differences generally reduce with the decrease of model grid spacing. At 100-m horizontal grid spacing, the vertical momentum fluxes induced by the model-resolved eddies and the associated eddy exchange coefficients in the eyewall simulated by the LESs with different 3D SGS mixing schemes are fairly consistent. Although with uncertainties, the decomposition in terms of eddy scales suggests that sub-kilometer eddies are mainly responsible for the vertical turbulent transport within the boundary layer (~1 km depth following the conventional definition) whereas eddies greater than 1 km become the dominant contributors to the vertical momentum transport above the boundary layer in the eyewall. The strong dependence of vertical turbulent transport on eddy scales suggests that the vertical turbulent mixing parameterization in mesoscale simulations of tropical cyclones is ultimately a scale-sensitive problem.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Atmospheric research
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In a tropical cyclone (TC), turbulence not only exists in the planetary boundary layer (PBL) but also can be generated above the PBL by the cloud processes in the eyewall and rainbands. It is found that the Hurricane Analysis and Forecast System (HAFS), a new multi-scale operational model for TC prediction, fails to capture the intense turbulent mixing in eyewall and rainband clouds due to a poor estimation of static stability in clouds. The problem is fixed by including the effects of multi-phase water in the stability calculation. Simulations of 21 TCs and tropical storms in the North Atlantic basin of 2016–2019 hurricane seasons totaling 118 forecast cycles show that the stability correction substantially improves HAFS's skill in predicting storm track and intensity. Analyses of HAFS's simulations of Hurricane Michael (2018) show that the positive tendency of vortex's tangential wind resulting from the radially inward transport of absolute vorticity dominates the eddy correlation tendencies induced by the model-resolved asymmetric eddies and serves as a main mechanism for the rapid intensification of Michael. The sub-grid scale (SGS) turbulent transport above the PBL in the eyewall plays a pivotal role in initiating a positive feedback among the eyewall convection, mean secondary overturning circulation, vortex acceleration via the inward transport of absolute vorticity, surface evaporation, and radial convergence of moisture in the PBL. Without the SGS transport above the PBL, the model-resolved vertical transport alone may not be sufficient in initiating the positive feedback underlying the rapid intensification of TCs. 
    more » « less
  2. Abstract. The tropical tropopause layer (TTL) is a sea of vertical motions. Convectively generated gravity waves create vertical winds on scales of a few to thousands of kilometers as they propagate in a stable atmosphere. Turbulence from gravity wave breaking, radiatively driven convection, and Kelvin–Helmholtz instabilities stirs up the TTL on the kilometer scale. TTL cirrus clouds, which moderate the water vapor concentration in the TTL and stratosphere, form in the cold phases of large-scale (> 100 km) wave activity. It has been proposed in several modeling studies that small-scale (< 100 km) vertical motions control the ice crystal number concentration and the dehydration efficiency of TTL cirrus clouds. Here, we present the first observational evidence for this. High-rate vertical winds measured by aircraft are a valuable and underutilized tool for constraining small-scale TTL vertical wind variability, examining its impacts on TTL cirrus clouds, and evaluating atmospheric models. We use 20 Hz data from five National Aeronautics and Space Administration (NASA) campaigns to quantify small-scale vertical wind variability in the TTL and to see how it varies with ice water content, distance from deep convective cores, and height in the TTL. We find that 1 Hz vertical winds are well represented by a normal distribution, with a standard deviation of 0.2–0.4 m s−1. Consistent with a previous observational study that analyzed two out of the five aircraft campaigns that we analyze here, we find that turbulence is enhanced over the tropical west Pacific and within 100 km of convection and is most common in the lower TTL (14–15.5 km), closer to deep convection, and in the upper TTL (15.5–17 km), further from deep convection. An algorithm to classify turbulence and long-wavelength (5 km < λ < 100 km) and short-wavelength (λ < 5 km) gravity wave activity during level flight legs is applied to data from the Airborne Tropical TRopopause EXperiment (ATTREX). The most commonly sampled conditions are (1) a quiescent atmosphere with negligible small-scale vertical wind variability, (2) long-wavelength gravity wave activity (LW GWA), and (3) LW GWA with turbulence. Turbulence rarely occurs in the absence of gravity wave activity. Cirrus clouds with ice crystal number concentrations exceeding 20 L−1 and ice water content exceeding 1 mg m−3 are rare in a quiescent atmosphere but about 20 times more likely when there is gravity wave activity and 50 times more likely when there is also turbulence, confirming the results of the aforementioned modeling studies. Our observational analysis shows that small-scale gravity waves strongly influence the ice crystal number concentration and ice water content within TTL cirrus clouds. Global storm-resolving models have recently been run with horizontal grid spacing between 1 and 10 km, which is sufficient to resolve some small-scale gravity wave activity. We evaluate simulated vertical wind spectra (10–100 km) from four global storm-resolving simulations that have horizontal grid spacing of 3–5 km with aircraft observations from ATTREX. We find that all four models have too little resolved vertical wind at horizontal wavelengths between 10 and 100 km and thus too little small-scale gravity wave activity, although the bias is much less pronounced in global SAM than in the other models. We expect that deficient small-scale gravity wave activity significantly limits the realism of simulated ice microphysics in these models and that improved representation requires moving to finer horizontal and vertical grid spacing. 
    more » « less
  3. This talk presents results from the authors’ recent work on evaluating the role of turbulence and boundary-layer parameterizations on hurricane intensification. We show that observation-based modification of these physical parameterizations significantly improved the HWRF intensity forecast. Turbulent mixing in both the vertical and horizontal directions are found to be crucial for hurricane spin-up dynamics in 3D numerical simulations and HWRF forecasts. Vertical turbulent mixing regulates the inflow strength and the location of boundary-layer convergence that in turns regulates the distribution of deep convection and the intensification of the whole hurricane vortex. Convergence of angular momentum in the boundary layer that is a key component of the hurricane spin-up theory is also found to be regulated by vertical turbulent mixing in connection to the boundary layer inflow. Horizontal turbulent mixing, on the other hand, mainly influences the eddy momentum flux inside the radius of the maximum wind speed in the angular momentum budget. The effect of horizontal turbulent mixing on the convergence of angular momentum is on smoothing the radial gradient of the angular momentum when the horizontal mixing length is large. In a sheared storm, both the vertical and horizontal turbulent mixing affect vortex and shear interaction in terms of the evolution of vortex tilt and boundary-layer recovery processes. 
    more » « less
  4. Abstract

    Kilometer‐scale grid spacing is increasingly being used in regional numerical weather prediction and climate simulation. This resolution range is in the terra incognita, where energetic eddies are partially resolved and turbulence parameterization is a challenge. The Smagorinsky and turbulence kinetic energy 1.5‐order models are commonly used at this resolution range, but, as traditional eddy‐diffusivity models, they can only represent forward‐scattering turbulence (downgradient fluxes), whereas the dynamic reconstruction model (DRM), based on explicit filtering, permits countergradient fluxes. Here we perform large‐eddy simulation of deep convection with 100‐m horizontal grid spacing and use these results to evaluate the performance of turbulence schemes at 1‐km horizontal resolution. The Smagorinsky and turbulence kinetic energy 1.5 schemes produce large‐amplitude errors at 1‐km resolution, due to excessively large eddy diffusivities attributable to the formulation of the squared moist Brunt‐Väisälä frequency (). With this formulation in cloudy regions, eddy diffusivity can be excessively increased in “unstable” regions, which produce downward (downgradient) heat flux in a conditionally unstable environment leading to destabilization and further amplification of eddy diffusivities. A more appropriate criterion based on saturation mixing ratio helps eliminate this problem. However, shallow clouds cannot be simulated well in any case at 1‐km resolution with the traditional models, whereas DRM allows for countergradient heat flux for both shallow and deep convection and predicts the distribution of clouds and fluxes satisfactorily. This is because DRM employs an eddy diffusivity model that is dynamically adjusted and a reconstruction approach that allows countergradient fluxes.

    more » « less
  5. Abstract. The fundamental mechanism underlying tropical cyclone (TC) intensification may be understood from the conservation of absolute angular momentum, where the primary circulation of a TC is driven by the torque acting on air parcels resulting from asymmetric eddy processes, including turbulence. While turbulence is commonly regarded as a flow feature pertaining to the planetary boundary layer (PBL), intense turbulent mixing generated by cloud processes also exists above the PBL in the eyewall and rainbands. Unlike the eddy forcing within the PBL that is negative definite, the sign of eyewall/rainband eddy forcing above the PBL is indefinite and thus provides a possible mechanism to spin up a TC vortex. In this study, we show that the Hurricane Weather Research & forecasting (HWRF) model, one of the operational models used for TC prediction, is unable to generate appropriate sub-grid-scale (SGS) eddy forcing above the PBL due to lack of consideration of intense turbulent mixing generated by the eyewall and rainband clouds. Incorporating an in-cloud turbulent mixing parameterization in the PBL scheme notably improves HWRF's skills on predicting rapid changes in intensity for several past major hurricanes. While the analyses show that the SGS eddy forcing above the PBL is only about one-fifth of the model-resolved eddy forcing, the simulated TC vortex inner-core structure and the associated model-resolved eddy forcing exhibit a substantial dependence on the parameterized SGS eddy processes. The results highlight the importance of eyewall/rainband SGS eddy forcing to numerical prediction of TC intensification, including rapid intensification at the current resolution of operational models.

    more » « less