skip to main content

Title: Key Elements of Turbulence Closures for Simulating Deep Convection at Kilometer‐Scale Resolution

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
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Page Range / eLocation ID:
p. 818-838
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Turbulence parameterization plays a critical role in the simulation of many weather regimes. For challenging cases such as the stratocumulus-capped boundary layer (SCBL), traditional schemes can produce unrealistic results even when a fine large-eddy-simulation (LES) resolution is used. Here we present an implicit generalized linear algebraic subfilter-scale model (iGLASS) to better represent unresolved turbulence in the simulation of the atmospheric boundary layer, at both standard LES and so-called terra incognita (TI) resolutions. The latter refers to a range of model resolutions where turbulent eddies are only partially resolved, and therefore the simulated processes are sensitive to the representation of unresolved turbulence. iGLASS is based on the truncated conservation equations of subfilter-scale (SFS) fluxes, but it integrates the full equations of the SFS turbulence kinetic energy and potential energy to retain “memory” of the SFS turbulence. Our evaluations suggest iGLASS can perform significantly better than traditional eddy-diffusivity models and exhibit skills comparable to the dynamic reconstruction model (DRM). For a neutral boundary layer case run at LES resolution, the simulation using iGLASS exhibits a wind profile that reasonably matches the similarity-theory solution. For an SCBL case with 5-m vertical resolution, iGLASS maintains more realistic cloud water profiles and boundary layer structure than traditional schemes. The SCBL case is also tested at TI resolution, and iGLASS also exhibits superior performance. iGLASS permits significant backscatter, whereas traditional models allow forward scatter (diffusion) only. As a physics-based approach, iGLASS appears to be a viable alternative for turbulence parameterization.

    more » « less
  2. 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
  3. Abstract

    While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement.

    Significance Statement

    In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

    more » « less
  4. Abstract

    Depth‐averaged eddy buoyancy diffusivities across continental shelves and slopes are investigated using a suite of eddy‐resolving, process‐oriented simulations of prograde frontal currents characterized by isopycnals tilted in the opposite direction to the seafloor, a flow regime commonly found along continental margins under downwelling‐favorable winds or occupied by buoyant boundary currents. The diagnosed cross‐slope eddy diffusivity varies by up to three orders of magnitude, decaying fromin the relatively flat‐bottomed region toover the steep continental slope, consistent with previously reported suppression effects of steep topography on baroclinic eddy fluxes. To theoretically constrain the simulated cross‐slope eddy fluxes, we examine extant scalings for eddy buoyancy diffusivities across prograde shelf/slope fronts and in flat‐bottomed oceans. Among all tested scalings, the GEOMETRIC framework developed by D. P. Marshall et al. (2012, and a parametrically similar Eady scale‐based scaling proposed by Jansen et al. (2015, most accurately reproduce the diagnosed eddy diffusivities across the entire shelf‐to‐open‐ocean analysis regions in our simulations. This result relies upon the incorporation of the topographic suppression effects on eddy fluxes, quantified via analytical functions of the slope Burger number, into the scaling prefactor coefficients. The predictive skills of the GEOMETRIC and Eady scale‐based scalings are shown to be insensitive to the presence of along‐slope topographic corrugations. This work lays a foundation for parameterizing eddy buoyancy fluxes across large‐scale prograde shelf/slope fronts in coarse‐resolution ocean models.

    more » « less
  5. Abstract

    We demonstrate that an extended eddy‐diffusivity mass‐flux (EDMF) scheme can be used as a unified parameterization of subgrid‐scale turbulence and convection across a range of dynamical regimes, from dry convective boundary layers, through shallow convection, to deep convection. Central to achieving this unified representation of subgrid‐scale motions are entrainment and detrainment closures. We model entrainment and detrainment rates as a combination of turbulent and dynamical processes. Turbulent entrainment/detrainment is represented as downgradient diffusion between plumes and their environment. Dynamical entrainment/detrainment is proportional to a ratio of a relative buoyancy of a plume and a vertical velocity scale, that is modulated by heuristic nondimensional functions which represent their relative magnitudes and the enhanced detrainment due to evaporation from clouds in drier environment. We first evaluate the closures offline against entrainment and detrainment rates diagnosed from large‐eddy simulations (LES) in which tracers are used to identify plumes, their turbulent environment, and mass and tracer exchanges between them. The LES are of canonical test cases of a dry convective boundary layer, shallow convection, and deep convection, thus spanning a broad range of regimes. We then compare the LES with the full EDMF scheme, including the new closures, in a single column model (SCM). The results show good agreement between the SCM and LES in quantities that are key for climate models, including thermodynamic profiles, cloud liquid water profiles, and profiles of higher moments of turbulent statistics. The SCM also captures well the diurnal cycle of convection and the onset of precipitation.

    more » « less