THERMODYNAMIC PROPERTIES AND ROLE OF THERMALS IN THE STRATUS-TOPPED BOUNDARY LAYER*

From results of two large-eddy simulation of stratus-topped boundary layer,the structure,thermodynamic properties and role of thermals are investigated by using conditional sampling methods,which divided the thermals into the warm/moist,cool/dry,warm/dry and cool/moist events.The results show that the main turbulent circulation in the stratus-topped boundary layer is composed of the warm/moist updraft and cool/dry downdraft.Below entrainment region,the warm/moist updrafts and cool/dry downdrafts are,respectively,positively and negatively buoyant and contribute most to total fluxes and variances.Evaporative cooling has important effect on the structure and thermodynamic properties of thermals in stratus-topped boundary layer.     

THERMODYNAMIC PROPERTIES AND ROLE OF THERMALS IN THE STRATUS-TOPPED BOUNDARY LAYER

From results of two large-eddy simulation of stratus-topped boundary layer,the structure,thermodynamic prop-erties and role of thermals are investigated by using conditional sampling methods,which divided the thermals into thewarm/moist,cool/dry,warm/dry and cool/moist events.The results show that the main turbulent circulation in thestratus-topped boundary layer is composed of the warm/moist updraft and cool/dry downdraft.Below entrainmentregion,the warm/moist updrafts and cool/dry downdrafts are,respectively,positively and negatively buoyant andcontribute most to total fluxes and variances.Evaporative cooling has important effect on the structure and thermodynamic properties of thermals instratus-topped boundary layer.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part II: Continuous Formulations of Mixing Rates for Mass-Flux Schemes

The conditional sampling of coherent structures in large-eddy simulations of the convective boundary layer (Couvreux et al. Boundary-layer Meteorol 134:441–458, 2010) is used to propose and evaluate formulations of fractional entrainment and detrainment rates for mass-flux schemes. The proposed formulations are physically-based and continuous from the surface to the top of clouds. Entrainment is related to the updraft vertical velocity divergence, while detrainment depends on the thermal vertical velocity, on buoyancy and on the moisture contrast between the mean plume and its environment. The proposed formulations are first directly evaluated in simulations of shallow clouds. They are then tested in single-column simulations with the thermal plume model, a mass-flux representation of boundary-layer thermals.     

Bivariate conditional sampling of buoyancy flux during an intense cold-air outbreak

Buoyancy fluxes in the marine atmospheric boundary layer (MABL) for the cloud street regime, observed during the Genesis of Atlantic Lows Experiment (GALE), have been analyzed using the technique of joint frequency distribution. For the lower half of the MABL, the results suggest that the buoyancy flux is mainly generated by the rising thermals and the sinking compensating ambient air, and is mainly consumed by the entrainment and detrainment of thermals, penetrative convection, and the entrainment from the MABL top.The results are compared to those from previous studies of mesoscale cellular convection (Air-Mass Transformation Experiment, AMTEX), the dry convective boundary layer, and the trade-wind MABL. For the lower MABL, the quadrant buoyancy fluxes, fractional coverages, and flux intensities are in good agreement with those of mesoscale cellular convection (AMTEX) and the dry convective boundary layer. The results suggest that, if the buoyancy flux is primarily driven by the temperature flux, the physical processes for generating buoyancy flux mentioned above are about the same for the lower boundary layers over land and ocean, even with different convective regimes. For the trade-wind MABL, the buoyancy flux is mainly driven by the moisture flux; the quadrant flux intensities are stronger than those of the other three studies except for the buoyant updrafts (thermals). These results suggest that the entrainment and detrainment of thermals are more effective in the trade-wind MABL than in the boundary layers driven by the temperature flux.Scale analysis of the buoyancy flux is in good agreement with that of AMTEX. For the lower half of the MABL, the buoyancy flux is mainly generated by the intermediate scale (200 m to 2 km), which includes the dominant convective thermals in the surface layer and the mixed layer. The scale smaller than 200 m is important only in the surface layer. The scale larger than 2 km, which includes the roll vortices, increases its significance upward. While most of the positive and negative fluxes are associated with the updrafts for the intermediate scale, the downdrafts are as important as updrafts for the larger scale.ST Systems Corporation, Lanham, MD, 20706, U.S.A.     

Airborne measurements of turbulent trace gas fluxes and analysis of eddy structure in the convective boundary layer over complex terrain

Using the new high-frequency measurement equipment of the research aircraft DO 128, which is described in detail, turbulent vertical fluxes of ozone and nitric oxide have been calculated from data sampled during the ESCOMPTE program in the south of France. Based on airborne turbulence measurements, radiosonde data and surface energy balance measurements, the convective boundary layer (CBL) is examined under two different aspects. The analysis covers boundary-layer convection with respect to (i) the control of CBL depth by surface heating and synoptic scale influences, and (ii) the structure of convective plumes and their vertical transport of ozone and nitric oxides. The orographic structure of the terrain causes significant differences between planetary boundary layer (PBL) heights, which are found to exceed those of terrain height variations on average. A comparison of boundary-layer flux profiles as well as mean quantities over flat and complex terrain and also under different pollution situations and weather conditions shows relationships between vertical gradients and corresponding turbulent fluxes. Generally, NO_x transports are directed upward independent of the terrain, since primary emission sources are located near the ground. For ozone, negative fluxes are common in the lower CBL in accordance with the deposition of O₃ at the surface.The detailed structure of thermals, which largely carry out vertical transports in the boundary layer, are examined with a conditional sampling technique. Updrafts mostly contain warm, moist and NO_x loaded air, while the ozone transport by thermals alternates with the background ozone gradient. Evidence for handover processes of trace gases to the free atmosphere can be found in the case of existing gradients across the boundary-layer top. An analysis of the size of eddies suggests the possibility of some influence of the heterogeneous terrain in mountainous area on the length scales of eddies.     

Characteristics of secondary circulations over an inhomogeneous surface simulated with large-eddy simulation

Large-eddy simulation is used to study secondary circulations in the convective boundary layer modulated as a result of horizontally varying surface properties and surface heat fluxes over flat terrain. The presence of heat flux heterogeneity and its alignment with respect to geostrophic wind influences the formation, strength and orientation of organized thermals. Results show boundary-attached roll formation along heat flux maxima in the streamwise direction. The streamwise organization of the updrafts and downdrafts formed downwind of heterogeneities leads to counter-rotating secondary circulations in the crosswind plane. The distribution of resolved-scale pressure deviations shows large pressure gradients in the crosswind plane. Spanwise and vertical velocity variances and heat flux profiles depict considerable spatial variability compared to a homogeneous forest simulation. Secondary circulations are observed for various ambient wind scenarios parallel and perpendicular to heterogeneities. In the presence of increased wind speed, thermals emerging from the heat flux heterogeneity are elongated, and organize along and downwind of large-scale heterogeneity in the streamwise direction. Simulation with a reduced heat flux shows a shallower circulation with a lower aspect ratio. Point measurements of heat flux inside the roll circulation could be overestimated by up to 15–25% compared to a homogeneous case.     

涡旋相关法测定湍流通量偏低的研究

针对野外实验所发现的不同观测法测定地表能通量不平衡问题,进行了均匀加热大气边界层的大涡模拟实验.用模拟的湍流风、温度和湿度涨落的时间序列证实,对流边界层低频涡普遍存在,并经常以一簇一簇热泡的形式出现.风速较小时,有限时长的取样不足以捕捉低频涡的贡献,可造成涡旋相关法测量的统计量异常偏低.仿照涡旋相关法的步骤进行数据处理发现,经去除平均或趋势计算的温度和湿度通量偏低程度在边界层下部随观测高度的增高而显著,其中尤以湿度通量为甚.其结果在一定程度上可以解释低风速条件下地表能通量测量的不闭合问题,但是尚不能完全解释诸如青藏高原实验出现的严重不闭合.文中对此作了探讨性的讨论.     

Mean Profiles of Moisture Fluxes in Snow-Filled Boundary Layers

Profiles of moisture fluxes have been examined for convective boundary layers containing clouds and snow, using data derived from aircraft measurements taken on four dates during the 1983/1984 University of Chicago lake-effect snow project. Flux profiles were derived from vertical stacks of aircraft cross-wind flight legs taken at various heights over Lake Michigan near the downwind shore. It was found that, if ice processes are taken into account, profiles of potential temperature and water content were very similar to those presented in past studies of convective boundary layers strongly heated from below. Profiles of total water content and equivalent potential temperature adjusted for ice were nearly invariant with height, except very near the top of the boundary layer, suggesting that internal boundary-layer mixing processes were rapid relative to the rates at which heat and vapour were transported into the boundary layer through entrainment and surface fluxes. Ice was found to play a significant, measurable role in boundary-layer moisture fluxes. It was estimated that 40 to 57% of the upward vapour flux was returned to the surface in the form of snow, converting about 45 to 64% of the surface latent heat flux into sensible heat in the snow-producing process. Assuming advective fluxes are relatively small (thought to be appropriate after the first few tens of km over the lake as suggested by past studies), the boundary layer was found to warm at a rate faster than could be explained by surface heat fluxes and latent heat releases alone, the remainder of the heating presumably coming from radiational processes and entrainment. Discussions of moisture phase change processes throughout the boundary layer and estimates of errors of these flux measurements are presented.     

Impact of coherent eddies on airborne measurements of vertical turbulent fluxes

During the Hydrological-Atmospheric Pilot Experiment (HAPEX)-Sahel, which took place in Niger in the transitional period between the wet and dry seasons, two French aircraft probed the Sahelian boundary layer to measure sensible and latent heat fluxes. The measurements over the Niamey area often revealed organised structures of a few km scale that were associated with both thermals and dry intrusions. We study the impact of these coherent structures using a single day’s aircraft-measured fluxes and a numerical simulation of that day with a mesoscale model. The numerical simulation at high horizontal resolution (250 m) contains structures that evolve from streaks in the early morning to cells by noon. This simulation shows distribution, variance and skewness similar to the observations. In particular, the numerical simulation shows dry intrusions that can penetrate deeply into the atmospheric boundary layer (ABL), and even reach the surface in some cases, which is in accordance with the observed highly negatively skewed water vapour fluctuations. Dry intrusions and thermals organised at a few km scale give skewed flux statistics and can introduce large errors in measured fluxes. We use the numerical simulation to: (i) evaluate the contribution of the organised structures to the total flux, and (ii) estimate the impact of the organised structures on the systematic and random errors resulting from the 1D sampling of the aircraft as opposed to the 2D numerical simulation estimate. We find a significant contribution by the organised structures to the total resolved fluxes. When rolls occur, and for a leg length of about 30 times the ABL depth, the 1D sampled flux is shown to be sometimes 20% lower than the corresponding 2D flux when the 1D sampling direction is the same as the main axis of the rolls, whereas the systematic error is much lower when the direction of the leg is transverse to the rolls. In the case of cells, an underestimate of around 10% can still be observed with the 1D approach independent of direction, due to poor sampling of the energy-containing scales.     

Analyses of variance and flux budgets in cumulus-topped boundary layers

A large-eddy simulation (LES) of a cumulus-topped boundary layer (CTBL) has been made, based on observations gathered on 10 September 1974 during the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE). On this day there were light winds, making buoyancy production the only important mechanism for generating turbulent kinetic energy. Cumulus clouds with a depth of a few hundred meters covered about 10% of the sky. A comparison between LES results and observed variances and fluxes will be made.In order to understand the dynamics of the cumulus, we will investigate the physical mechanisms by which variances and fluxes are produced in CTBL's. This will be done by analyzing the budget equations which describe the time evolution of the variances and fluxes. The production, transport and destruction terms, that appear in these equations, will be calculated explicitly with the LES model.Budgets are calculated from the GATE simulation and two slightly modified GATE simulations, in which only the initial water vapor specific humidity is increased. As a result in these two modified GATE simulations the vertical extent of the clouds and the cloud cover is larger. In this way, the influence of the cumulus clouds on the budgets can be investigated more thoroughly. Finally, the budgets are compared with the results obtained from an earlier simulation based on the Puerto Rico Field Experiment and other results found in the literature.     

Boundary Conditions for Scalar (Co)Variances over Heterogeneous Surfaces

The problem of boundary conditions for the variances and covariances of scalar quantities (e.g., temperature and humidity) at the underlying surface is considered. If the surface is treated as horizontally homogeneous, Monin–Obukhov similarity suggests the Neumann boundary conditions that set the surface fluxes of scalar variances and covariances to zero. Over heterogeneous surfaces, these boundary conditions are not a viable choice since the spatial variability of various surface and soil characteristics, such as the ground fluxes of heat and moisture and the surface radiation balance, is not accounted for. Boundary conditions are developed that are consistent with the tile approach used to compute scalar (and momentum) fluxes over heterogeneous surfaces. To this end, the third-order transport terms (fluxes of variances) are examined analytically using a triple decomposition of fluctuating velocity and scalars into the grid-box mean, the fluctuation of tile-mean quantity about the grid-box mean, and the sub-tile fluctuation. The effect of the proposed boundary conditions on mixing in an archetypical stably-stratified boundary layer is illustrated with a single-column numerical experiment. The proposed boundary conditions should be applied in atmospheric models that utilize turbulence parametrization schemes with transport equations for scalar variances and covariances including the third-order turbulent transport (diffusion) terms.     

A simple model of the atmospheric boundary layer; sensitivity to surface evaporation

A simple formulation of the boundary layer is developed for use in large-scale models and other situations where simplicity is required. The formulation is suited for use in models where some resolution is possible within the boundary layer, but where the resolution is insufficient for resolving the detailed boundary-layer structure and overlying capping inversion. Surface fluxes are represented in terms of similarity theory while turbulent diffusivities above the surface layer are formulated in terms of bulk similarity considerations and matching conditions at the top of the surface layer. The boundary-layer depth is expressed in terms of a bulk Richardson number which is modified to include the influence of thermals. Attention is devoted to the interrelationship between predicted boundary-layer growth, the turbulent diffusivity profile, countergradient heat flux and truncation errors.The model predicts growth of the convectively mixed layer reasonably well and is well-behaved in cases of weak surface heat flux and transitions between stable and unstable cases. The evolution of the modelled boundary layer is studied for different ratios of surface evaporation to potential evaporation. Typical variations of surface evaporation result in a much greater variation in boundary-layer depth than that caused by the choice of the boundary-layer depth formulation.     

The role of thermals in the convective boundary layer

Detailed measurements of the structure of thermals throughout the convective boundary layer were obtained from the NCAR Electra aircraft over the ocean during the Air Mass Transformation Experiment (AMTEX). Humidity was used as an indicator of thermals. The variables were first high-pass filtered with a 5 km cutoff digital filter to eliminate mesoscale variations. Segments of the 5 min (30 km length) horizontal flight legs with humidity greater than half the standard deviation of humidity fluctuations for that leg were defined as thermals. This was found to be a better indicator of thermals than temperature in the upper part of the boundary layer since the temperature in a thermal is cooler than its environment in the upper part of the boundary layer. Using mixed-layer scaling, the normalized length and number of thermals were found to scale with the 1/3 and -1/3 powers, respectively, of normalized height, while vertical velocity and temperature scaled according to similarity predictions in the free convection region of the surface layer. The observational results presented here extend throughout the entire mixed layer. Using these results in the equation for mean updraft velocity of a field of thermals, the sum of the vertical pressure gradient and edge-effect terms can be estimated. This residual term is found to be important throughout most of the boundary layer. The magnitude of the divergence of vertical velocity variance within a thermal is found to be larger than the magnitude of the mean updraft velocity term throughout most of the mixed layer.Part of this work was completed while visiting Risø National Laboratory, Denmark.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Three-dimensional numerical study of turbulence in an entraining mixed layer

The mean structure calculated by a three-dimensional numerical model of a heated planetary boundary layer, in simulation of DAY 33 of the Australian Wangara data, has been previously described. The present study supplements it by describing properties of the calculated turbulence.A major finding is the importance of entrainment upon turbulence statistics relating to specific humidity, relative to those for potential temperature. The variances, skewness and spectra of velocity, temperature and humidity are presented, as are budget equations for kinetic energy, temperature and humidity variances and heat/moisture fluxes. These are interpreted with regard to the relative importance of the surface flux vs the flux due to entrainment at the top of the mixed layer, and in regard to the structure which would occur if the entrainment were to vanish.The Rotte-type closure assumption is tested for the correlation between the pressure fluctuation and the vertical gradient of vertical velocity, potential temperature, or specific humidity, and found to be qualitatively correct except near the top of the mixed layer.NCAR is sponsored by the National Science Foundation (U.S.A.).     

Atmosphere-ocean interaction in mid-latitude storms

Summary Surface fluxes of heat, latent heat, and momentum, and entrainment fluxes and vertical motion at the top of the boundary layer have been calculated for limited regions of several mid-latitude ocean storms. Results have been combined to describe distributions of boundary layer processes which are characteristic of such storms. Surface heat fluxes have important effects in the region west of cold or occluded fronts and are relatively unimportant within a band of about 200 km width east of fronts. Entrainment in pre-frontal regions is driven largely by vertical shear at the top of the boundary layer, while in post-frontal regions it is driven largely by surface heat flux. Boundary layers are well defined in regions more than roughly 200 km east or west of fronts; but closer to fronts boundary layers are not well defined due to the combined effects of entrainment, condensation, and vertical motion associated with the distribution of surface stress.With 12 Figures     

Boundary-layer evolution and nocturnal inversion dispersal—Part II

Observations from an instrumented aircraft are used to study the small-scale structure of turbulence and convection in well-mixed boundary layers and the erosion processes in the nocturnally-formed inversions above them. The ways in which turbulence statistics for temperature, humidity and vertical velocity scale with height in the mixed layer are compared with the results of a three-dimensional model by Deardorff (1974a, b), and agreement is found in many aspects. Conditional sampling enables the statistics of thermals and their environment to be considered separately and, in particular, shows that the mode of the vertical velocity in thermals markedly decreases with height in the upper half of the mixed layer. Thermals may be recognized equally readily by either their excess of temperature or humidity. Transfers of heat and moisture through the nocturnal inversions influence the structure of the upper region of the mixed layer and there is strong evidence that these transfer processes are turbulent and not organized on scales similar to convective thermals.     

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).     

Air mass modification over Lake Michigan

The modification of a relatively cold air mass over the warm water of Lake Michigan is studied by using a two-dimensional nonlinear mesoscale model. Considerable amounts of heat and water vapor are supplied from the water surface to the lower atmosphere by turbulent eddies. A convective mixed layer develops and grows toward the downwind region with stratocumulus clouds over the lake.The model simulates the warming and moistening of the mixed layer, the development of a boundary layer, the divergence and convergence of wind near the coastlines, and the turbulent fluxes.The model warming of the mixed layer across the lake was about 2.2 °K and the moistening of the mixed layer was about 0.8 g kg^–1, which are comparable to 2.7 °K and 0.8 g kg^–1 observed by Lenschow (1973). The convective boundary layer, which includes the cloud layer, subcloud layer, and superadiabatic layer near the water surface, is well simulated. The tilt of the inversion which coincides with the cloud top is also well reproduced. When a prescribed cooling rate is applied at the cloud top, stronger turbulence and a deeper cloud layer are generated. Without the cooling, the cloud is shallow and the shape of the cloud base is determined by surface conditions. The rise of the inversion is due to upward vertical motion, and deepening of the convective layer in the downwind region.     

Observations Of The Turbulence Structure Within Two Stratocumulus-Topped, Marine Boundary Layers

Two situations observed during the second Aerosol Characterization Experiment (ACE-2) are analysed from aircraft measurements in the broken stratocumulus (Sc)-topped marine boundary layer. The first one (26 June 1997), characterized by a non-polluted, oceanic air mass, presents a decoupling between the Sc layer (1400–1520 m) and the turbulent mixed layer, this latter extending from the surface up to 580 m. In contrast, the second case (9 July 1997), during which continental air had been advected over the experimental area, presents a well-coupled layer extending from the surface up to the top of the Sc layer(910 m). This coupling, uncommon in this area in the middle of the day, isrelated to the relative shallowness of the boundary layer. For both situations,it is shown that the turbulent fluxes can be computed with reasonably goodaccuracy (better than 10 %), taking into account both the random and thesystematic errors involved in the eddy-correlation technique. Estimationof random error is based on the computation of the integral scale of thecovariance, and systematic error is estimated from the parameterizationof Mann and Lenschow. The fluxes show that the buoyancy, as a sourceof turbulence, is due to latent heat flux rather than sensible heat flux,with values comparable to previous experiments in the Azores-Canariesbasin. In addition, we propose a method to analyse, for coupled situations,the relationship between the fractional cloudiness and the organization ofthe turbulent field below the clouds. This method is based on a conditionalsampling technique. It is shown that this organization cannot be deducedfrom the analysis of the velocity signal, which is dominated by turbulence.However, when the signals are conditionally sampled according to thepresence or absence of clouds, a weak cloud-related organization can beshown, and the cloud-related transports quantified; the values found areof the order of 10 % of the total transfers, i.e. the same order of magnitude asthe errors on the total flux computation. The method developed is thereforepromising, provided that the uncertainties can be reduced by analyzing a highamount of data.     

On the structure of convection

Conditional sampling methods are used to investigate the structure of thermals in a dry atmospheric convection layer. Temperature and vertical velocity data obtained from an aircraft flying above the surface layer are analysed. The thermal structure is found to be consistent with free-convection similarity theory.