A Modelling Study of Flux Imbalance and the Influence of Entrainment in the Convective Boundary Layer

It is frequently observed in field experiments that the eddy covariance heat fluxes are systematically underestimated as compared to the available energy. The flux imbalance problem is investigated using the NCAR’s large-eddy simulation (LES) model imbedded with an online scheme to calculate Reynolds-averaged fluxes. A top–down and a bottom–up tracer are implemented into the LES model to quantify the influence of entrainment and bottom–up diffusion processes on flux imbalance. The results show that the flux imbalance follows a set of universal functions that capture the exponential decreasing dependence on u _*/w _*, where u _* and w _* are friction velocity and the convective velocity scale, respectively, and an elliptic relationship to z/z _i , where z _i is the mixing-layer height. The source location in the boundary layer is an important factor controlling the imbalance magnitude and its horizontal and vertical distributions. The flux imbalance of heat and the bottom–up tracer is tightly related to turbulent coherent structures, whereas for the top–down diffusion, such relations are weak to nonexistent. Our results are broadly consistent with previous studies on the flux imbalance problem, suggesting that the published results are robust and are not artefacts of numerical schemes.     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

Statistics of Scalar Fields in the Atmospheric Boundary Layer Based on Large-Eddy Simulations. Part II: Forced Convection

Forced convection in a quasi-steady atmospheric boundary layer is investigated based on a large-eddy simulation (LES) model. The performed simulations show that in the upper portion of the mixed layer the dimensionless (in terms of mixed layer scales) vertical gradients of temperature, humidity, and wind velocity depend on the dimensionless height z/z _i and the Reech number Rn. The peak values of variances and covariances at the top of the mixed layer, scaled in terms of the interfacial scales, are functions of the interfacial Richardson number Ri. As a result expressions for the entrainment rates, in the case when the interfacial layer has a finite depth, and a condition for the presence of moistening or drying regimes in the mixed layer, are derived. Profiles of dimensionless scalar moments in the mixed layer are proposed to be expressed in terms of two empirical similarity functions F _m and F _i , dependent on dimensionless height z/z _i , and the interfacial Richardson number Ri. The obtained similarity expressions adequately approximate the LES profiles of scalar statistics, and properly represent the impact of stability, shear, and entrainment. They are also consistent with the parameterization proposed for free convection in the first part of this paper.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

A Study in Search of Interconnection between Surface Parameters and Surrounding Synoptic And Subsynoptic Features

The paper reveals that the variations in parameters like u*, the scaling velocity and θ*. The scaling tempera-ture during the various phases of monsoon might be linked with subsynoptic features. The rise in u* is mainly connected with the presence of lower tropospheric cyclonic vorticity over a subsynoptic scale of the site. However the variations in θ* is mainly linked with the various phases of monsoon and θ* shows a sharp rise in presence of low level convective cloud.Besides the correlation studies of u and u*, θv and θv* , θv-θv0 and θv* are undertaken. The correlation be?tween θv and θv* is poor. In other two cases correlations are good. Besides u/u* , has shown good coefficient of variation values within the ζ range.     

Scaling Variances of Scalars in a Convective Boundary Layer Under Different Entrainment Regimes

For the presentation and analysis of atmospheric boundary-layer (ABL) data, scales are used to non-dimensionalise the observed quantities and independent variables. Usually, the ABL height, surface sensible heat flux and surface scalar flux are used. This works well, so long as the absolute values of the entrainment ratio for both the scalar and temperature are similar. The entrainment ratio for temperature naturally ranges from −0.4 to −0.1. However, the entrainment ratio for passive scalars can vary widely in magnitude and sign. Then the entrainment flux becomes relevant as well. The only customary scalar scale that takes into account both the surface flux and the entrainment flux is the bulk scalar scale, but this scale is not well-behaved for large negative entrainment ratios and for an entrainment ratio equal to −1. We derive a new scalar scale, using previously published large-eddy simulation results for the convective ABL. The scale is derived under the constraint that scaled scalar variance profiles are similar at those heights where the variance producing mechanisms are identical (i.e., either near the entrainment layer or near the surface). The new scale takes into account that scalar variance in the ABL is not only related to the surface flux of that scalar, but to the scalar entrainment flux as well. Furthermore, it takes into account that the production of variance by the entrainment flux is an order of magnitude larger than the production of variance by the surface flux (per unit flux). Other desirable features of the new scale are that it is always positive (which is relevant when scaling standard deviations) and that the scaled variances are always of order 1–10.     

Parameterization of Sheared Convective Entrainment in the First-Order Jump Model: Evaluation Through Large-Eddy Simulation

In this note, two different approaches are used to estimate the entrainment-flux to surface-flux ratio for a sheared convective boundary layer (CBL); both are derived under the framework of the first-order jump model (FOM). That suggested by Sun and Wang (SW approach) has the advantage that there is no empirical constant included, though the dynamics are described in an implicit manner. The second, which was proposed by Kim et al. and Pino et al. (KP approach), explicitly characterizes the dynamics of the sheared entrainment, but uncertainties are induced through the empirical constants. Their performances in parameterizing the CBL growth rate are compared and discussed, and a new value of the parameter A ₃ in the KP approach is suggested. Large-eddy simulation (LES) data are employed to test both approaches: simulations are conducted for the CBL growing under varying conditions of surface roughness, free-atmospheric stratification, and wind shear, and data used when the turbulence is in steady state. The predicted entrainment rates in each case are tested against the LES data. The results show that the SW approach describes the evolution of the sheared CBL quite well, and the KP approach also reproduces the growth of the CBL reasonably, so long as the value of A ₃ is modified to 0.6.     

对流边界层大涡模式的改进及对夹卷速度的研究

对已建的对流边界层(CBL)大涡模式进行了改进，将次网格闭合方案改为次网格能量闭合，并考虑了水汽的源汇项和水汽相变潜热的作用。通过对均匀下垫面上由热扰动发展的对流边界层的模拟及与实验结果的比较表明，模式较好地模拟了对流边界层的主要物理结构，较好地反映了各物理量之间的对应关系。本文在一定的对流理查森数(Ri*)范围内给出了一些算例，对无量纲夹卷速度(We/W*)进行了研究。结果表明，无量纲夹卷速度随地表热通量(Qs)的增大而增大，随对流边界层上部温度递减率（γ)的增大而减小。当9.06≤Ri*≤45.29时，无量纲夹卷速度We/W*可以拟合成A(Ri*)-1的形式，其中A=0.226。并且与我们的对流槽实验结果，Sullivan等人的大涡模拟结果以及Deardorff等人的对流槽实验结果作了比较，四者吻合较好。     

A Simple Model for the Vertical Transport of Reactive Species in the Convective Atmospheric Boundary Layer

We have developed a simple, steady-state, one-dimensional second-order closure model to obtain continuous profiles of turbulent fluxes and mean concentrations of non-conserved scalars in a convective boundary layer without shear. As a basic tool we first set up a model for conserved species with standard parameterizations. This leads to formulations for profiles of the turbulent diffusivity and the ratio of temperature-scalar covariance to the flux of the passive scalar. The model is then extended to solving, in terms of profiles of mean concentrations and fluxes, the NO _x –O₃ triad problem. The chemical reactions involve one first-order reaction, the destruction of NO₂ with decay time τ, and one second-order reaction, the destruction of NO and O₃ with the reaction constant k. Since the fluxes of the sum concentrations of NO _x = NO + NO₂ and O₃ + NO₂ turn out to be constant throughout the boundary layer, the problem reduces to solving two differential equations for the concentration and the flux of NO₂. The boundary conditions are the three surface fluxes and the fluxes at the top of the boundary layer, the last obtained from the entrainment velocity, and the concentration differences between the free troposphere and the top of the boundary layer. The equations are solved in a dimensionless form by using 1/(kτ) as the concentration unit, the depth h of the boundary layer as the length unit, the convective velocity scale w _* as the velocity unit, and the surface temperature flux divided by w _* as the temperature unit. Special care has been devoted to the inclusion of the scalar–scalar covariance between the concentrations of O₃ and NO. Sample calculations show that the fluxes of the reactive species deviate significantly from those of non-reactive species. Further, the diffusivities, defined by minus the flux divided by the concentration gradient may become negative for reactive species in contrast to those of non-reactive species, which in the present model are never negative.     

Local Scales of Turbulence in the Stable Boundary Layer

Local, gradient-based scales, which contain the vertical velocity and temperature variances, as well as the potential temperature gradient, but do not include fluxes, are tested using data collected during the CASES-99 experiment. The observations show that the scaling based on the temperature variance produces relatively smaller scatter of empirical points. The resulting dimensionless statistical moments approach constant values for sufficiently large values of the Richardson number Ri. This allows one to derive predictions for the Monin–Obukhov similarity functions φ _m and φ _h , the Prandtl number Pr and the flux Richardson number Rf in weak turbulence regime.     

The role of molecular diffusion in the deepening of the mixed layer

Turbulent mixing across heat-stratified density interfaces was studied in the laboratory using oscillating-grid generated turbulence. The aim was to study the transition between the entrainment regimes dominated by interfacial wave-breaking and molecular diffusion, and to study the characteristics of the latter. It was observed that, above a critical Richardson number Ri_c, which depends on the Peclet number Pe, the mixing due to wave breaking disappears and that Ri_c Pe⁻ⁿ, where the mean value of the exponent n is approximately . Above Ri_c, the entrainment is molecular-diffusion dominated and takes place through a sequence of events: the buoyancy gradient of the initially sharp density interface is weakened by molecular diffusion until the mixed-layer eddies can engulf a portion of the interfacial layer wherefore the interface sharpens again. Thus, the entrainment events are recurrent with a rate-controlling diffusion stage between them. An entrainment law of the form E Ri⁻²Pe⁻², where E is the entrainment coefficient and Ri is the Richardson number, is suggested for the diffusion-dominated entrainment regime.     

Laboratory Experiments on Convective Entrainment Using a Saline Water Tank

Entrainment fluxes in a shear-free convective boundary layer have been measured with a saline water tank set-up. The experiments were targeted towards measuring the entrainment behaviour for medium to high Richardson numbers and use a two-layer design, i.e. two stacked non-stratified (neutral) layers with different densities. With laser induced fluorescence (LIF), the entrainment flux of a fluorescent dye is measured for bulk Richardson numbers in the range 30–260. It is proposed that a carefully chosen combination of top-down and bottom-up processes improves the accuracy of LIF-based entrainment observations. The observed entrainment fluxes are about an order of magnitude lower than reported for thermal water tanks: the derived buoyancy entrainment ratio, $A$ , is found to be $A \approx 0.02$ , which is to be compared with $A\approx 0.25$ for a thermal convection tank (Deardorff et al., J Fluid Mech 100:41–64, 1980). An extensive discussion is devoted to the influence of the Reynolds and Prandtl numbers in laboratory experiments on entrainment.     

An Evaluation of the Flux–Gradient Relationship in the Stable Boundary Layer

Data collected during the SHEBA and CASES-99 field programs are employed to examine the flux–gradient relationship for wind speed and temperature in the stably stratified boundary layer. The gradient-based and flux-based similarity functions are assessed in terms of the Richardson number Ri and the stability parameter z/Λ_*, z being height and Λ_* the local Obukhov length. The resulting functions are expressed in an analytical form, which is essentially unaffected by self-correlation, when thermal stratification is strong. Turbulence within the stably stratified boundary layer is classified into four regimes: “nearly-neutral” (0 < z/Λ_* < 0.02), “weakly-stable” (0.02 < z/Λ_* < 0.6), “very-stable” (0.6 < z/Λ_* < 50), and “extremely-stable” (z/Λ_* > 50). The flux-based similarity functions for gradients are constant in “nearly-neutral” conditions. In the “very-stable” regime, the dimensionless gradients are exponential, and proportional to (z/Λ_*)^3/5. The existence of scaling laws in “extremely-stable” conditions is doubtful. The Prandtl number Pr decreases from 0.9 in nearly-neutral conditions and to about 0.7 in the very-stable regime. The necessary condition for the presence of steady-state turbulence is Ri < 0.7.     

Parameterization of sheared entrainment in a well-developed CBL. Part I: Evaluation of the scheme through large-eddy simulations

The entrainment flux ratio A _e and the inversion layer (IL) thickness are two key parameters in a mixed layer model. A _e is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of A _e is derived from the TKE budget in the firstorder model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness. The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized A _e and IL thickness agree well with the LES outputs.     

Nearly steady convection and the boundary-layer budgets of water vapor and sensible heat

The budgets of water vapor and sensible heat in the convective atmospheric boundary (mixed) layer are analyzed by means of a simple slab approach adapted to steady large-scale advective conditions with radiation and cloud activity. The entrainment flux for sensible heat is assumed to be a linear function of the surface flux. The flux of water vapor at the top of the mixed layer is parameterized by extending the first-order Betts-Deardorff approach, i.e., by adopting linear changes for both the specific humidity and the flux across the mixed layer and across the inversion layer of finite thickness. In this way the dissimilarity of sensible heat and water vapor transport in the mixed layer can be taken into account. The experimental data were obtained from the Air Mass Transformation Experiment (AMTEX). The entrainment constant for sensible heat at the top of the mixed layer was found to have values similar to those observed in other weakly convective situations, i.e., around 0.4 to 0.6. This appears to indicate that the effect of mechanical turbulence was not negligible; however, the inclusion of this effect in the formulation did not improve the correlation. In contrast to the first-order approach, the zero-order approach, i. e., the jump equation commonly used for the flux of a scalar at the inversion, (ovwc ) _h = w_e c (where w _e is the entrainment velocity and c the concentration jump across the inversion), was found to be invalid and incapable of describing the data.     

An Alternative Eddy-Viscosity Model for the Horizontally Uniform Atmospheric Boundary Layer

This paper explores the utility of specifying the eddy viscosity for the horizontally uniform boundary layer as the product of the variance of vertical velocity and an empirical time scale τ _w , as opposed to the more usual formulation where k is the turbulent kinetic energy (TKE), λ _k is a length scale and α is a dimensionless coefficient. Simulations were compared with the observations on Day 33 of the Wangara experiment, and with a plausible specification of τ _w (or λ _k ) each model simulated convective boundary-layer development reasonably well, although the closure produced a more realistic width for the entrainment layer. Under the light winds of Day 33, and with the onset of evening cooling, an excessively shallow and strongly-stratified nocturnal inversion developed, and limited its own further deepening. Boundary-layer models that neglect radiative heat transport and parametrize convective transport by eddy viscosity closure are prone to this runaway (unstable) feedback when forced by a negative (i.e. downward) surface flux of sensible heat.     

Parameterization of Sheared Entrainment in a Well-developed CBL.Part Ⅱ:A Simple Model for Predicting the Growth Rate of the CBL

Following the parameterization of sheared entrainment obtained in the companion paper, Liu et al. (2016), the present study aims to further investigate the characteristics of entrainment, and develop a simple model for predicting the growth rate of a well-developed and sheared CBL. The relative stratification, defined as the ratio of the stratification in the free atmosphere to that in the entrainment zone, is found to be a function of entrainment flux ratio (A _e). This leads to a simple expression of the entrainment rate, in which A _e needs to be parameterized. According to the results in Liu et al. (2016), A _e can be simply expressed as the ratio of the convective velocity scale in the sheared CBL to that in the shear-free CBL. The parameterization of the convective velocity scale in the sheared CBL is obtained by analytically solving the bulk model with several assumptions and approximations. Results indicate that the entrainment process is influenced by the dynamic effect, the interaction between mean shear and environmental stratification, and one other term that includes the Coriolis effect. These three parameterizations constitute a simple model for predicting the growth rate of a well-developed and sheared CBL. This model is validated by outputs of LESs, and the results show that it performs satisfactorily. Compared with bulk models, this model does not need to solve a set of equations for the CBL. It is more convenient to apply in numerical models.     

Vertical Structure of Selected Turbulence Characteristics above an Urban Canopy

Summary In this paper the results of an urban measurement campaign are presented. The experiment took place from July 1995 to February 1996 in Basel, Switzerland. A total of more than 2000 undisturbed 30-minute runs of simultaneous measurements of the fluctuations of the wind vector u′, v′, w′ and the sonic temperature θ _s ′ at three different heights (z=36, 50 and 76 m a.g.l.) are analysed with respect to the integral statistics and their spectral behaviour. Estimates of the zero plane displacement height d calculated by the temperature variance method yield a value of 22 m for the two lower levels, which corresponds to 0.92 h (the mean height of the roughness elements). At all three measurement heights the dimensionless standard deviation σ _w /u _* is systematically smaller than the Monin-Obukhov similarity function for the inertial sublayer, however, deviations are smaller compared to other urban turbulence studies. The σ_θ/θ_* values follow the inertial sublayer prediction very close for the two lowest levels, while at the uppermost level significant deviations are observed. Profiles of normalized velocity and temperature variances show a clear dependence on stability. The profile of friction velocity u _* is similar to the profiles reported in other urban studies with a maximum around z/h=2.1. Spectral characteristics of the wind components in general show a clear dependence on stability and dimensionless measurement height z/h with a shift of the spectral peak to lower frequencies as thermal stability changes from stable to unstable conditions and as z/h decreases. Velocity spectra follow the −2/3 slope in the inertial subrange region and the ratios of spectral energy densities S _w (f)/S _u (f) approach the value of 4/3 required for local isotropy in the inertial subrange. Velocity spectra and spectral peaks fit best to the well established surface layer spectra from Kaimal et al. (1972) at the uppermost level at z/h=3.2. Received September 26, 1997 Revised February 15, 1998     

Entrainment Processes in the Convective Boundary Layer with Varying Wind Shear

Large-eddy simulations (LES) are performed to investigate the entrainment andthe structure of the inversion layer of the convective boundary layer (CBL) withvarying wind shears. Three CBLs are generated with the constant surface kinematicheat flux of 0.05 K m s^-1 and varying geostrophic wind speeds from 5 to 15m s^-1. Heat flux profiles show that the maximum entrainment heat flux as afraction of the surface heat flux increases from 0.13 to 0.30 in magnitude withincreasing wind shear. The thickness of the entrainment layer, relative to the depthof the well-mixed layer, increases substantially from 0.36 to 0.73 with increasingwind shear. The identification of vortices and extensive flow visualizations nearthe entrainment layer show that concentrated vortices perpendicular to the meanboundary-layer wind direction are identified in the capping inversion layer for thecase of strong wind shear. These vortices are found to develop along the mean winddirections over strong updrafts, which are generated by convective rolls and to appearas large-scale wavy motions similar to billows generated by the Kelvin–Helmholtzinstability. Quadrant analysis of the heat flux shows that in the case of strong windshear, large fluctuations of temperature and vertical velocity generated by largeamplitude wavy motions result in greater heat flux at each quadrant than that inthe weak wind shear case.