An official website of the United States government
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
A Generalized Mixing Length Closure for Eddy‐Diffusivity Mass‐Flux Schemes of Turbulence and Convection
Abstract Because of their limited spatial resolution, numerical weather prediction and climate models have to rely on parameterizations to represent atmospheric turbulence and convection. Historically, largely independent approaches have been used to represent boundary layer turbulence and convection, neglecting important interactions at the subgrid scale. Here we build on an eddy‐diffusivity mass‐flux (EDMF) scheme that represents all subgrid‐scale mixing in a unified manner, partitioning subgrid‐scale fluctuations into contributions from local diffusive mixing and coherent advective structures and allowing them to interact within a single framework. The EDMF scheme requires closures for the interaction between the turbulent environment and the plumes and for local mixing. A second‐order equation for turbulence kinetic energy (TKE) provides one ingredient for the diffusive local mixing closure, leaving a mixing length to be parameterized. Here, we propose a new mixing length formulation, based on constraints derived from the TKE balance. It expresses local mixing in terms of the same physical processes in all regimes of boundary layer flow. The formulation is tested at a range of resolutions and across a wide range of boundary layer regimes, including a stably stratified boundary layer, a stratocumulus‐topped marine boundary layer, and dry convection. Comparison with large eddy simulations (LES) shows that the EDMF scheme with this diffusive mixing parameterization accurately captures the structure of the boundary layer and clouds in all cases considered.
more »
« less
- Award ID(s):
- 1835860
- PAR ID:
- 10452250
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 12
- Issue:
- 11
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract In the inner core of a tropical cyclone, turbulence not only exists in the boundary layer (BL) but can also be generated above the BL by eyewall and rainband clouds. Thus, the treatment of vertical turbulent mixing must go beyond the conventional scope of the BL. The turbulence schemes formulated based on the turbulent kinetic energy (TKE) are attractive as they are applicable to both deep and shallow convection regimes in the tropical cyclone (TC) inner core provided that the TKE production and dissipation can be appropriately determined. However, TKE schemes are not self-closed. They must be closed by an empirically prescribed vertical profile of mixing length. This motivates this study to investigate the sensitivity of the simulated TC intensification to the sloping curvature and asymptotic length scale of mixing length, the two parameters that determine the vertical distribution of a prescribed mixing length. To tackle the problem, both idealized and real-case TC simulations are performed. The results show that the simulated TC intensification is sensitive to the sloping curvature of mixing length but only exhibits marginal sensitivity to the asymptotic length scale. The underlying reasons for such sensitivities are explored analytically based on the Mellor and Yamada level-2 turbulence model and the analyses of azimuthal-mean tangential wind budget. The results highlight the uncertainty and importance of mixing length in the numerical prediction of TCs and suggest that future research should focus on searching for physical constraints on mixing length, particularly in the low- to midtroposphere, using observations and large-eddy simulations. Significance StatementThe parametric representation of subgrid-scale turbulent mixing is one of the major sources of uncertainty in numerical predictions of tropical cyclones (TCs). This study investigates how the numerical prediction of TC intensification is affected by the turbulent mixing length, a length scale that is required to close a turbulence scheme formulated based on the turbulent kinetic energy (TKE). The research highlights the uncertainty and importance of mixing length in numerical prediction of TCs and suggests that future research should focus on searching for physical constraints on the mixing length, particularly in the low- to midtroposphere, using observations and large-eddy simulations.more » « less
-
Abstract Turbulence parameterizations for convective boundary layer in coarse‐scale atmospheric models usually consider a combination of the eddy‐diffusive transport and a non‐local transport, typically in the form of a mass flux term, such as the widely adopted eddy‐diffusivity mass‐flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large‐eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass‐flux transport can be used to approximate the eddy‐diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy‐diffusivity and mass‐flux transport theories in convective conditions, which also quantifies departure from the Monin‐Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS.more » « less
-
Abstract Subgrid-scale turbulence in numerical weather prediction models is typically handled by a PBL parameterization. These schemes attempt to represent turbulent mixing processes occurring below the resolvable scale of the model grid in the vertical direction, and they act upon temperature, moisture, and momentum within the boundary layer. This study varies the PBL mixing strength within 4-km WRF simulations of a 26–29 January 2015 snowstorm to assess the sensitivity of baroclinic cyclones to eddy diffusivity intensity. The bulk critical Richardson number for unstable regimes is varied between 0.0 and 0.25 within the YSU PBL scheme as a way of directly altering the depth and magnitude of subgrid-scale turbulent mixing. Results suggest that varying the bulk critical Richardson number is similar to selecting a different PBL parameterization. Differences in boundary layer moisture availability, arising from reduced entrainment of dry, free tropospheric air, lead to variations in the magnitude of latent heat release above the warm frontal region, producing stronger upper-tropospheric downstream ridging in simulations with less PBL mixing. The more amplified flow pattern impedes the northeastward propagation of the surface cyclone and results in a westward shift of precipitation. In addition, trajectory analysis indicates that ascending parcels in the less-mixing simulations condense more water vapor and terminate at a higher potential temperature level than do ascending parcels in the more-mixing simulations, suggesting stronger latent heat release when PBL mixing is reduced. These results suggest that spread within ensemble forecast systems may be improved by perturbing PBL mixing parameters that are not well constrained.more » « less
-
Abstract Modeling experiments and field campaigns have evaluated shallow convective mixing as a potential constraint on the low‐cloud climate feedback, which is critical for establishing climate sensitivity. Yet the apparent relationship between low‐cloud fraction and shallow convective mixing differs substantially among general circulation models (GCMs), large eddy simulations, and both remote sensing and in situ observations. Here, we consider how changes in GCMs' representations of subgrid‐scale vertical moist fluxes can alter the cloud‐mixing relationship. Using vertical profiles of water vapor isotope ratios (δD) to characterize the strength of shallow convective mixing in a manner that can be compared directly to satellite observations, we evaluate the cloud‐mixing relationship produced in tiered experiments with the Community Atmosphere Model (CAM). From versions 5 to 6 of CAM, the most notable physics change is CLUBB, a scheme that unifies the representation of shallow convection and boundary layer turbulence through a joint probability density function (PDF) for subgrid velocity and moisture. CLUBB reduces the covariance between low‐cloud fraction and shallow convective mixing, producing a bivariate distribution that is more similar in character to monthly averaged satellite observations. Using parameter sensitivity experiments, we argue that CLUBB's ability to simulate skewness in the distribution of vertical velocity produces more isolated but stronger moist updrafts, which reduce the grid‐mean low‐cloud fraction while maintaining efficient hydrological connectivity between the boundary layer and the free troposphere. These results suggest that mixing is not an effective predictor of low‐cloud feedback in GCMs with PDF closure schemes.more » « less
