scholarly journals An Integrated TKE-Based Eddy Diffusivity/Mass Flux Boundary Layer Closure for the Dry Convective Boundary Layer

2011 ◽  
Vol 68 (7) ◽  
pp. 1526-1540 ◽  
Author(s):  
Marcin L. Witek ◽  
Joao Teixeira ◽  
Georgios Matheou

Abstract This study presents a new approach to the eddy diffusivity/mass flux (EDMF) framework for the modeling of convective boundary layers. At the root of EDMF lies a decomposition of turbulent transport mechanisms into strong ascending updrafts and smaller-scale turbulent motions. The turbulent fluxes can be therefore described using two conventional approaches: mass flux (MF) for the organized thermals and eddy diffusivity (ED) for the remaining turbulent field. Since the intensities of both MF and ED transports depend on the kinetic energy of the turbulent motions, it seems reasonable to formulate an EDMF framework based on turbulent kinetic energy (TKE). Such an approach allows for more physical and less arbitrary formulations of parameters in the model. In this study the EDMF–TKE coupling is achieved through the use of (i) a new parameterization for the lateral entrainment coefficient ɛ and (ii) the MF contribution to the buoyancy source of TKE. Some other important features of the EDMF parameterization presented here include a revised mixing length formulation and Monin–Obukhov stability scaling for the surface layer. The scheme is implemented in a one-dimensional (1D) model. Several cases of dry convective boundary layers (CBL) with different surface sensible heat fluxes in the free-convection limit are investigated. Results are compared to large-eddy simulation (LES). Good agreement between LES and the 1D model is achieved with respect to mean profiles, boundary layer evolution, and updraft characteristics. Some disagreements between the models are found to most likely relate to deficiencies in the TKE simulation in the 1D model. Comparison with other previously established ɛ parameterizations shows that the new TKE-based formulation leads to equally accurate, and in many respects better, simulation of the CBL. The encouraging results obtained with the proposed EDMF framework indicate that full integration of EDMF with higher-order closures is possible and can further improve boundary layer simulations.

2011 ◽  
Vol 68 (10) ◽  
pp. 2385-2394 ◽  
Author(s):  
Marcin L. Witek ◽  
Joao Teixeira ◽  
Georgios Matheou

Abstract In this study a new approach to the vertical transport of the turbulent kinetic energy (TKE) is proposed. The principal idea behind the new parameterization is that organized updrafts or convective plumes play an important role in transferring TKE vertically within convectively driven boundary layers. The parameterization is derived by applying an updraft environment decomposition to the vertical velocity triple correlation term in the TKE prognostic equation. The additional mass flux (MF) term that results from this decomposition closely resembles the features of the TKE transport diagnosed from the large-eddy simulation (LES) and accounts for 97% of the LES-diagnosed transport when the updraft fraction is set to 0.13. Another advantage of the MF term is that it is a function of the updraft vertical velocity and can be readily calculated using already existing parameterization. The new MF approach, combined with several eddy diffusivity (ED) formulations, is implemented into a simplified 1D TKE prognostic model. The 1D model results, compared against LES simulations of dry convective boundary layers, show substantial improvement in representing the vertical structure of TKE. The new combined ED–MF parameterization, as well as the MF term alone, surpasses in accuracy the ED parameterizations. The proposed TKE transport parameterization shows large potential of improving TKE simulations in mesoscale and global circulation models.


2007 ◽  
Vol 64 (4) ◽  
pp. 1230-1248 ◽  
Author(s):  
A. Pier Siebesma ◽  
Pedro M. M. Soares ◽  
João Teixeira

Abstract A better conceptual understanding and more realistic parameterizations of convective boundary layers in climate and weather prediction models have been major challenges in meteorological research. In particular, parameterizations of the dry convective boundary layer, in spite of the absence of water phase-changes and its consequent simplicity as compared to moist convection, typically suffer from problems in attempting to represent realistically the boundary layer growth and what is often referred to as countergradient fluxes. The eddy-diffusivity (ED) approach has been relatively successful in representing some characteristics of neutral boundary layers and surface layers in general. The mass-flux (MF) approach, on the other hand, has been used for the parameterization of shallow and deep moist convection. In this paper, a new approach that relies on a combination of the ED and MF parameterizations (EDMF) is proposed for the dry convective boundary layer. It is shown that the EDMF approach follows naturally from a decomposition of the turbulent fluxes into 1) a part that includes strong organized updrafts, and 2) a remaining turbulent field. At the basis of the EDMF approach is the concept that nonlocal subgrid transport due to the strong updrafts is taken into account by the MF approach, while the remaining transport is taken into account by an ED closure. Large-eddy simulation (LES) results of the dry convective boundary layer are used to support the theoretical framework of this new approach and to determine the parameters of the EDMF model. The performance of the new formulation is evaluated against LES results, and it is shown that the EDMF closure is able to reproduce the main properties of dry convective boundary layers in a realistic manner. Furthermore, it will be shown that this approach has strong advantages over the more traditional countergradient approach, especially in the entrainment layer. As a result, this EDMF approach opens the way to parameterize the clear and cumulus-topped boundary layer in a simple and unified way.


2013 ◽  
Vol 70 (7) ◽  
pp. 1929-1953 ◽  
Author(s):  
Kay Sušelj ◽  
João Teixeira ◽  
Daniel Chung

Abstract A single-column model (SCM) is developed for representing moist convective boundary layers. The key component of the SCM is the parameterization of subgrid-scale vertical mixing, which is based on a stochastic eddy-diffusivity/mass-flux (EDMF) approach. In the EDMF framework, turbulent fluxes are calculated as a sum of the turbulent kinetic energy–based eddy-diffusivity component and a mass-flux component. The mass flux is modeled as a fixed number of steady-state plumes. The main challenge of the mass-flux model is to properly represent cumulus clouds, which are modeled as moist plumes. The solutions have to account for a realistic representation of condensation within the plumes and of lateral entrainment into the plumes. At the level of mean condensation within the updraft, the joint pdf of moist conserved variables and vertical velocity is used to estimate the proportion of dry and moist plumes and is sampled in a Monte Carlo way creating a predefined number of plumes. The lateral entrainment rate is modeled as a stochastic process resulting in a realistic decrease of the convective cloudiness with height above cloud base. In addition to the EDMF scheme, the following processes are included in the SCM: a pdf-based parameterization of subgrid-scale condensation, a simple longwave radiation, and one-dimensional dynamics. Note that in this approach there are two distinct pdfs, one representing the variability of updraft properties and the other one the variability of thermodynamic properties of the surrounding environment. The authors show that the model is able to capture the essential features of moist boundary layers, ranging from stratocumulus to shallow-cumulus regimes. Detailed comparisons, which include pdfs, profiles, and integrated budgets with the Barbados Oceanographic and Meteorological Experiment (BOMEX), Dynamics and Chemistry of Marine Stratocumulus (DYCOMS), and steady-state large-eddy simulation (LES) cases, are discussed to confirm the quality of the present approach.


2005 ◽  
Vol 44 (9) ◽  
pp. 1436-1452 ◽  
Author(s):  
Wayne M. Angevine

Abstract A scheme is described that provides an integrated description of turbulent transport in free convective boundary layers with shallow cumulus. The scheme uses a mass-flux formulation, as is commonly found in cumulus schemes, and a 1.5-order closure, involving turbulent kinetic energy and eddy diffusivity. Both components are active in both the subcloud and cloud layers. The scheme is called “mass flux–diffusion.” In the subcloud layer, the mass-flux component provides nonlocal transport. The scheme combines elements from schemes that are conceptually similar but differ in detail. An entraining plume model is used to find the mass flux. The mass flux is continuous through the cloud base. The lateral fractional entrainment rate is constant with height, while the detrainment-rate profile reduces the mass flux smoothly to zero at the cloud top. The eddy diffusivity comes from a turbulent kinetic energy–length scale formulation. The scheme has been implemented in a simple one-dimensional (single column) model. Results of simulations of a standard case that has been used for other model intercomparisons [Atmospheric Radiation Measurement (ARM), 21 June 1997] are shown and indicate that the scheme provides good results. The model also simulates the profile of a conserved scalar; this capability is applied to a case from the 1999 Southern Oxidants Study Nashville (Tennessee) experiment, where it produces good simulations of vertical profiles of carbon monoxide in a cloud-topped boundary layer.


2016 ◽  
Vol 144 (3) ◽  
pp. 1161-1177 ◽  
Author(s):  
Hyeyum Hailey Shin ◽  
Jimy Dudhia

Abstract Planetary boundary layer (PBL) parameterizations in mesoscale models have been developed for horizontal resolutions that cannot resolve any turbulence in the PBL, and evaluation of these parameterizations has been focused on profiles of mean and parameterized flux. Meanwhile, the recent increase in computing power has been allowing numerical weather prediction (NWP) at horizontal grid spacings finer than 1 km, at which kilometer-scale large eddies in the convective PBL are partly resolvable. This study evaluates the performance of convective PBL parameterizations in the Weather Research and Forecasting (WRF) Model at subkilometer grid spacings. The evaluation focuses on resolved turbulence statistics, considering expectations for improvement in the resolved fields by using the fine meshes. The parameterizations include four nonlocal schemes—Yonsei University (YSU), asymmetric convective model 2 (ACM2), eddy diffusivity mass flux (EDMF), and total energy mass flux (TEMF)—and one local scheme, the Mellor–Yamada–Nakanishi–Niino (MYNN) level-2.5 model. Key findings are as follows: 1) None of the PBL schemes is scale-aware. Instead, each has its own best performing resolution in parameterizing subgrid-scale (SGS) vertical transport and resolving eddies, and the resolution appears to be different between heat and momentum. 2) All the selected schemes reproduce total vertical heat transport well, as resolved transport compensates differences of the parameterized SGS transport from the reference SGS transport. This interaction between the resolved and SGS parts is not found in momentum. 3) Those schemes that more accurately reproduce one feature (e.g., thermodynamic transport, momentum transport, energy spectrum, or probability density function of resolved vertical velocity) do not necessarily perform well for other aspects.


2011 ◽  
Vol 68 (12) ◽  
pp. 3112-3131 ◽  
Author(s):  
Rachel Honnert ◽  
Valéry Masson ◽  
Fleur Couvreux

Abstract Turbulence is well represented by atmospheric models at very fine grid sizes, from 10 to 100 m, for which turbulent movements are mainly resolved, and by atmospheric models with grid sizes greater than 2 km, for which those movements are entirely parameterized. But what happens at intermediate scales, Wyngaard’s so-called terra incognita? Here an original method is presented that provides a new diagnostic by calculating the subgrid and resolved parts of five variables at different scales: turbulent kinetic energy (TKE), heat and moisture fluxes, and potential temperature and mixing ratio variances. They are established at intermediate scales for dry and cumulus-topped convective boundary layers. The similarity theorem allows the determination of the dimensionless variables of the problem. When the subgrid and resolved parts are studied, a new dimensionless variable, the dimensionless mesh size , needs to be added to the Deardorff free convective scaling variables, where h is the boundary layer height and hc is the height of the cloud layer. Similarity functions for the subgrid and resolved parts are assumed to be the product of the similarity function of the total (subgrid plus resolved) variables and a “partial” similarity function that depends only on . In order to determine the partial similarity function form, large-eddy simulations (LES) of five dry and cloudy convective boundary layers are used. The resolved and subgrid parts of the variables at coarser grid sizes are then deduced from the LES fields. The evolution of the subgrid and resolved parts in the boundary layer with is as follows: fine grids mainly resolve variables. As the mesh becomes coarser, more eddies are subgrid. Finally, for very large meshes, turbulence is entirely subgrid. A scale therefore exists for which the subgrid and resolved parts are equal. This is obtained for in the case of TKE, 0.4 for the potential temperature variance, and 0.8 for the mixing ratio variance, indicating that the velocity structures are smaller than those for the potential temperature, which are smaller than those for the mixing ratio. Furthermore, boundary layers capped by convective clouds have structures larger than dry boundary layer ones as displayed by the scaling in the partial similarity functions. This new diagnostic gives a reference for evaluating current and future parameterizations at kilometric scales. As an illustration, the parameterizations of a mesoscale model are eventually evaluated at intermediate scales. In its standard version, the model produces too many resolved movements, as the turbulence scheme does not sufficiently represent the impact of the subgrid thermal. This is not true when a mass-flux scheme is introduced. However in this case, a completely subgrid thermal is modeled leading to an overestimation of the subgrid part.


2016 ◽  
Vol 73 (7) ◽  
pp. 2715-2727 ◽  
Author(s):  
Siddhartha Mukherjee ◽  
Jerôme Schalkwijk ◽  
Harmen J. J. Jonker

Abstract The predictability horizon of convective boundary layers is investigated in this study. Large-eddy simulation (LES) and direct numerical simulation (DNS) techniques are employed to probe the evolution of perturbations in identical twin simulations of a growing dry convective boundary layer. Error growth typical of chaotic systems is observed, marked by two phases. The first comprises an exponential error growth as , with δ0 as the initial error, δ(t) as the error at time t, and Λ as the Lyapunov exponent. This phase is independent of the perturbation wavenumber, and the perturbation energy grows following a self-similar spectral shape dominated by higher wavenumbers. The nondimensional error growth rate in this phase shows a strong dependence on the Reynolds number (Re). The second phase involves saturation of the error. Here, the error growth follows Lorenz dynamics with a slower saturation of successively larger scales. An analysis of the spectral decorrelation times reveals two regimes: an Re-independent regime for scales larger than the boundary layer height and an Re-dependent regime for scales smaller than , which are found to decorrelate substantially faster for increasing Reynolds numbers.


2010 ◽  
Vol 138 (7) ◽  
pp. 2895-2912 ◽  
Author(s):  
Wayne M. Angevine ◽  
Hongli Jiang ◽  
Thorsten Mauritsen

Abstract Comparisons between single-column (SCM) simulations with the total energy–mass flux boundary layer scheme (TEMF) and large-eddy simulations (LES) are shown for four cases from the Gulf of Mexico Atmospheric Composition and Climate Study (GoMACCS) 2006 field experiment in the vicinity of Houston, Texas. The SCM simulations were run with initial soundings and surface forcing identical to those in the LES, providing a clean comparison with the boundary layer scheme isolated from any other influences. Good agreement is found in the simulated vertical transport and resulting moisture profiles. Notable differences are seen in the cloud base and in the distribution of moisture between the lower and upper cloud layer. By the end of the simulations, TEMF has dried the subcloud layer and moistened the lower cloud layer more than LES. TEMF gives more realistic profiles for shallow cumulus conditions than traditional boundary layer schemes, which have no transport above the dry convective boundary layer. Changes to the formulation and its parameters from previous publications are discussed.


2011 ◽  
Vol 68 (4) ◽  
pp. 727-738 ◽  
Author(s):  
Robert van Driel ◽  
Harm J. J. Jonker

In this study the response of dry convective boundary layers to nonstationary surface heat fluxes is systematically investigated. This is relevant not only during sunset and sunrise but also, for example, when clouds modulate incoming solar radiation. Because the time scale of the associated change in surface heat fluxes may differ from case to case, the authors consider the generic situation of oscillatory surface heat fluxes with different frequencies and amplitudes and study the response of the boundary layer in terms of transfer functions. To this end both a mixed layer model (MLM) and a large-eddy simulation (LES) model are used; the latter is used to evaluate the predictive quality of the mixed layer model. The mixed layer model performs generally quite well for slow changes in the surface heat flux and provides analytical understanding of the transfer characteristics of the boundary layer such as amplitude and phase lag. For rapidly changing surface fluxes (i.e., changes within a time frame comparable to the large eddy turnover time), it proves important to account for the time it takes for the information to travel from the surface to higher levels of the boundary layer such as the inversion zone. As a follow-up to a 1997 study by Sorbjan, who showed that the conventional convective velocity scale is inadequate as a scaling quantity during the decay phase, this paper addresses the issue of defining, in (generic) transitional situations, a velocity scale that is solely based on the surface heat flux and its history.


Sign in / Sign up

Export Citation Format

Share Document