scholarly journals Predictability of Dry Convective Boundary Layers: An LES Study

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.

2006 ◽  
Vol 45 (9) ◽  
pp. 1224-1243 ◽  
Author(s):  
David Pino ◽  
Jordi Vilà-Guerau de Arellano ◽  
Si-Wan Kim

Abstract Dry convective boundary layers characterized by a significant wind shear on the surface and at the inversion are studied by means of the mixed-layer theory. Two different representations of the entrainment zone, each of which has a different closure of the entrainment heat flux, are considered. The simpler of the two is based on a sharp discontinuity at the inversion (zeroth-order jump), whereas the second one prescribes a finite depth of the inversion zone (first-order jump). Large-eddy simulation data are used to provide the initial conditions for the mixed-layer models, and to verify their results. Two different atmospheric boundary layers with different stratification in the free atmosphere are analyzed. It is shown that, despite the simplicity of the zeroth-order-jump model, it provides similar results to the first-order-jump model and can reproduce the evolution of the mixed-layer variables obtained by the large-eddy simulations in sheared convective boundary layers. The mixed-layer model with both closures compares better with the large-eddy simulation results in the atmospheric boundary layer characterized by a moderate wind shear and a weak temperature inversion. These results can be used to represent the flux of momentum, heat, and other scalars at the entrainment zone in general circulation or chemistry transport models.


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.


2014 ◽  
Vol 14 (23) ◽  
pp. 32491-32533 ◽  
Author(s):  
C. Darbieu ◽  
F. Lohou ◽  
M. Lothon ◽  
J. Vilà-Guerau de Arellano ◽  
F. Couvreux ◽  
...  

Abstract. We investigate the decay of planetary boundary layer (PBL) turbulence in the afternoon, from the time the surface buoyancy flux starts to decrease until sunset. Dense observations of mean and turbulent parameters were acquired during the Boundary Layer Late Afternoon and Sunset Turbulence (BLLAST) field experiment by several meteorological surface stations, sounding balloons, radars, lidars, and two aircraft flying extensively during the afternoon transition. We analyzed a case study based on some of those observations and Large-Eddy Simulation (LES) data focusing on the turbulent vertical structure throughout the afternoon transition. The decay of turbulence is quantified through the temporal and vertical evolution of (1) the turbulence kinetic energy (TKE), (2) the characteristic length scales of turbulence, (3) the shape of the turbulence spectra. A spectral analysis of LES data, airborne and surface measurements is performed in order to (1) characterize the variation of the turbulent decay with height and (2) study the distribution of turbulence over eddy size. This study points out the LES ability to reproduce the turbulence evolution throughout the afternoon. LES and observations agree that the afternoon transition can be divided in two phases: (1) a first phase during which the TKE decays with a low rate, with no significant change in turbulence characteristics, (2) a second phase characterized by a larger TKE decay rate and a change spectral shape, implying an evolution of eddy size distribution and energy cascade from low to high wavenumber. The changes observed either on TKE decay (during the first phase) or on the vertical wind spectra shape (during the second phase of the afternoon transition) occur first in the upper region of the PBL. The higher within the PBL, the stronger the spectra shape changes.


2016 ◽  
Vol 13 ◽  
pp. 63-67 ◽  
Author(s):  
Rachel Honnert

Abstract. Numerical weather prediction model forecasts at horizontal grid lengths in the range of 100 to 1 km are now possible. This range of scales is the "grey zone of turbulence". Previous studies, based on large-eddy simulation (LES) analysis from the MésoNH model, showed that some assumptions of some turbulence schemes on boundary-layer structures are not valid. Indeed, boundary-layer thermals are now partly resolved, and the subgrid remaining part of the thermals is possibly largely or completely absent from the model columns. First, some modifications of the equations of the shallow convection scheme have been tested in the MésoNH model and in an idealized version of the operational AROME model at resolutions coarser than 500 m. Secondly, although the turbulence is mainly vertical at mesoscale (>  2 km resolution), it is isotropic in LES (<  100 m resolution). It has been proved by LES analysis that, in convective boundary layers, the horizontal production of turbulence cannot be neglected at resolutions finer than half of the boundary-layer height. Thus, in the grey zone, fully unidirectional turbulence scheme should become tridirectional around 500 m resolution. At Météo-France, the dynamical turbulence is modelled by a K-gradient in LES as well as at mesoscale in both MésoNH and AROME, which needs mixing lengths in the formulation. Vertical and horizontal mixing lengths have been calculated from LES of neutral and convective cases at resolutions in the grey zone.


2011 ◽  
Vol 68 (1) ◽  
pp. 155-174 ◽  
Author(s):  
Alastair G. Williams ◽  
Wlodek Zahorowski ◽  
Scott Chambers ◽  
Alan Griffiths ◽  
Jörg M. Hacker ◽  
...  

Abstract Radon (222Rn) is a powerful natural tracer of mixing and exchange processes in the atmospheric boundary layer. The authors present and discuss the main features of a unique dataset of 50 high-resolution vertical radon profiles up to 3500 m above ground level, obtained in clear and cloudy daytime terrestrial boundary layers over an inland rural site in Australia using an instrumented motorized research glider. It is demonstrated that boundary layer radon profiles frequently exhibit a complex layered structure as a result of mixing and exchange processes of varying strengths and extents working in clear and cloudy conditions within the context of the diurnal cycle and the synoptic meteorology. Normalized aircraft radon measurements are presented, revealing the characteristic structure and variability of three major classes of daytime boundary layer: 1) dry convective boundary layers, 2) mixed layers topped with residual layers, and 3) convective boundary layers topped with coupled nonprecipitating clouds. Robust and unambiguous signatures of important atmospheric processes in the boundary layer are identifiable in the radon profiles, including “top-down” mixing associated with entrainment in clear-sky cases and strongly enhanced venting and subcloud-layer mixing when substantial active cumulus are present. In poorly mixed conditions, radon gradients in the daytime atmospheric surface layer significantly exceed those predicted by Monin–Obukhov similarity theory. In two case studies, it is demonstrated for the first time that a sequence of vertical radon profiles measured over the course of a single day can consistently reproduce major structural features of the evolving boundary layer.


2020 ◽  
Vol 77 (10) ◽  
pp. 3619-3630
Author(s):  
Jeremy A. Gibbs ◽  
Evgeni Fedorovich

AbstractWe extend our previous study, which dealt with structure functions of potential temperature fluctuations, and focus on the characteristics of second-order velocity structure functions and corresponding structure parameters in the atmospheric convective boundary layer. We consider the three previously reported methods to compute the structure parameters of turbulent velocity fields: the direct method, the true spectral method, and the approximate spectral method. The methods are evaluated using high-resolution gridded numerical data from large-eddy simulations of shear-free and shear-driven convective boundary layers. Results indicate that the direct and true spectral methods are more suitable than the approximate spectral method, which overestimates the structure parameters of velocity as a result of assuming the inertial-subrange shape of the velocity spectrum for all turbulence scales. Results also suggest that structure parameters of vertical velocity fluctuations are of limited utility because of violations of local isotropy, especially in shear-free convective boundary layers.


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.


2015 ◽  
Vol 15 (17) ◽  
pp. 10071-10086 ◽  
Author(s):  
C. Darbieu ◽  
F. Lohou ◽  
M. Lothon ◽  
J. Vilà-Guerau de Arellano ◽  
F. Couvreux ◽  
...  

Abstract. We investigate the decay of planetary boundary layer (PBL) turbulence in the afternoon, from the time the surface buoyancy flux starts to decrease until sunset. Dense observations of mean and turbulent parameters were acquired during the Boundary Layer Late Afternoon and Sunset Turbulence (BLLAST) field experiment by several meteorological surface stations, sounding balloons, radars, lidars and two aircraft during the afternoon transition. We analysed a case study based on some of these observations and large-eddy simulation (LES) data focusing on the turbulent vertical structure throughout the afternoon transition. The decay of turbulence is quantified through the temporal and vertical evolution of (1) the turbulence kinetic energy (TKE), (2) the characteristic length scales of turbulence and (3) the shape of the turbulence spectra. A spectral analysis of LES data, airborne and surface measurements is performed in order to characterize the variation in the turbulent decay with height and study the distribution of turbulence over eddy size. This study highlights the LES ability to reproduce the turbulence evolution throughout the afternoon. LESs and observations agree that the afternoon transition can be divided in two phases: (1) a first phase during which the TKE decays at a low rate, with no significant change in turbulence characteristics, and (2) a second phase characterized by a larger TKE decay rate and a change in spectral shape, implying an evolution of eddy size distribution and energy cascade from low to high wave number. The changes observed either in TKE decay (during the first phase) or in the vertical wind spectra shape (during the second phase of the afternoon transition) occur first in the upper region of the PBL. The higher within the PBL, the stronger the spectra shape changes.


Sign in / Sign up

Export Citation Format

Share Document