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.     

Resolved Versus Parametrized Boundary-Layer Plumes. Part I: A Parametrization-Oriented Conditional Sampling in Large-Eddy Simulations

A conditional sampling based on the combination of a passive tracer emitted at the surface and thermodynamic variables is proposed to characterise organized structures in large-eddy simulations of cloud-free and cloudy boundary layers. The sampling is evaluated against more traditional sampling of dry thermals or clouds. It enables the characterization of convective updrafts from the surface to the top of the boundary layer (or the top of cumulus clouds), describing in particular the transition from the sub-cloud to the cloud layer, and retrieves plume characteristics, entrainment and detrainment rates, variances and fluxes. This sampling is used to analyze the contribution of boundary-layer thermals to vertical fluxes and variances.     

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.     

Validation of the “Lokal-Modell” over heterogeneous land surfaces using aircraft-based measurements of the REEEFA experiment and comparison with micro-scale simulations

Summary An aircraft-based experimental investigation of the atmospheric boundary layer (ABL) structure and of the energy exchange processes over heterogeneous land surfaces is presented. The measurements are used for the validation of the mesoscale atmospheric model “Lokal-Modell” (LM) of the German Weather Service with 2.8 km resolution. In addition, high-resolution simulations using the non-hydrostatic model FOOT3DK with 250 m resolution are performed in order to resolve detailed surface heterogeneities. Two special observation periods in May 1999 show comparable convective boundary layer (CBL) conditions. For one case study vertical profiles and area averages of meteorological quantities and energy fluxes are investigated in detail. The measured net radiation is highly dependent on surface albedo, and the latent heat flux exhibits a strong temporal variability in the investigation area. A reduction of this variability is possible by aggregation of multiple flight patterns. To calculate surface fluxes from aircraft measurements at low altitude, turbulent energy fluxes were extrapolated to the ground by the budget method, which turned out to be well applicable for the sensible heat flux, but not for the latent flux. The development of the ABL is well captured by the LM simulation. The comparison of spatiotemporal averages shows an underestimation of the observed net radiation, which is mainly caused by thin low-level clouds in the LM compared to observed scattered CBL clouds. The sensible heat flux is reproduced very well, while the latent flux is highly overestimated especially above forests. The realistic representation of surface heterogeneities in the investigation area in the FOOT3DK simulations leads to improvements for the energy fluxes, but an overestimation of the latent heat flux still persists. A study of upscaling effects yields more structures than the LM fields when averaged to the same scale, which are partly caused by the non-linear effects of parameter aggregation on the LM scale.     

A Case Study on the Effects of Heterogeneous Soil Moisture on Mesoscale Boundary-Layer Structure in the Southern Great Plains, U.S.A. Part I: Simple Prognostic Model

The atmospheric boundary-layer (ABL) depth was observed by airborne lidar and balloon soundings during the Southern Great Plains 1997 field study (SGP97). This paper is Part I of a two-part case study examining the relationship of surface heterogeneity to observed ABL structure. Part I focuses on observations. During two days (12–13 July 1997) following rain, midday convective ABL depth varied by as much as 1.5 km across 400 km, even with moderate winds. Variability in ABL depth was driven primarily by the spatial variation in surface buoyancy flux as measured from short towers and aircraft within the SGP97 domain. Strong correlation was found between time-integrated buoyancy flux and airborne remotely sensed surface soil moisture for the two case-study days, but only a weak correlation was found between surface energy fluxes and vegetation greenness as measured by satellite. A simple prognostic one-dimensional ABL model was applied to test to what extent the soil moisture spatial heterogeneity explained the variation in north–south ABL depth across the SGP97 domain. The model was able to better predict mean ABL depth and variations on horizontal scales of approximately 100 km using observed soil moisture instead of constant soil moisture. Subsidence, advection, convergence/divergence and spatial variability of temperature inversion strength also contributed to ABL depth variations. In Part II, assimilation of high-resolution soil moisture into a three-dimensional mesoscale model (MM5) is discussed and shown to improve predictions of ABL structure. These results have implications for ABL models and the influence of soil moisture on mesoscale meteorology     

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.     

Spatial and Temporal Characteristics of Horizontal Rolls and Cells in the Atmospheric Boundary Layer Based on Radar and in Situ Observations

This paper is the second of a series devoted to the observation and analysis of coherent structures in the cloudy Atmospheric Boundary Layer (ABL) such as horizontal rolls or thermal cells. In the first paper, the TRAC (Turbulence Radar Aircraft Cells) experiment which is the observational support of this investigation based on coupled radar-aircraft measurements, was presented along with an overview of the main results of this campaign held in June 1993 in France. Here the analysis is focused on the spatial characteristics (length-scale, orientation ... ) of the coherent structures, their temporal evolution and vertical distribution deduced from the radar reflectivity fields acquired in clear air at several levels for five different ABLs. For that, an original and efficient image processing method able to extract the major mode of the organisation was developed. These characteristics are examined in relation to the dynamic and thermodynamic state of the ABL using mean and turbulent information as observed by the in-situ aircraft. These experimental results are the basis of a future Large Eddy Simulation modeling of an organised ABL which is in progress and will be the concern of the third paper in the series.     

Land Breeze and Thermals： A Scale Threshold to Distinguish Their Effects

Land breeze is a type of mesoscale circulation developed due to thermal forcing over a heterogeneous landscape. It can contribute to atmospheric dynamic and hydrologic processes through affecting heat and water fluxes on the land-atmosphere interface and generating shallow convective precipitation. If the scale of the landscape heterogeneity is smaller than a certain size, however, the resulting land breeze becomes weak and becomes mixed up with other thermal convections like thermals. This study seeks to identify a scale threshold to distinguish the effects between land breeze and thermals. Two-dimensional simulations were performed with the Regional Atmospheric Modeling System （RAMS） to simulate thermals and land breeze. Their horizontal scale features were analyzed using the wavelet transform. The thermals developed over a homogeneous landscape under dry or wet conditions have an initial scale of 2-5 km during their early stage of development. The scale jumps to 10-15 km when condensation occurs. The solution of an analytical model indicates that the reduced degree of atmospheric instability due to the release of condensation potential heat could be one of the contributing factors for the increase in scale. The land breeze, on the other hand, has a major scale identical to the size of the landscape heterogeneity throughout various stages of development. The results suggest that the effects of land breeze can be clearly distinguished from those of thermals only if the size of the landscape heterogeneity is larger than the scale threshold of about 5 km for dry atmospheric processes or about 15 km for moist ones.     

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.     

A Case Study of Rayleigh-Bénard Convection with Clouds

A turbulence data set collected by the research aircraft Hercules and Falcon in the planetary boundary layer (PBL) over the North Sea during Rayleigh-Bénard convection (RBC) is analysed. Altogether nearly three hundred cell passages at different levels and in two different flight directions were sampled.The convective boundary-layer height (H) was about 1 km, and the RBC cells had a diameter D of roughly 2–3 km, resulting in an aspect ratio A = D/H 2–3. This value is also found in the case of RBC in laboratory-scale flows, whereas most of the recent PBL experimental work reports convection PBL rolls with A 3 and mesoscale cellular convection (MCC) with A 10–40 over the oceans.The large number of RBC cell passages made it possible to composite their average structure. Due to the more complex three-dimensional structure and the importance of thermals to the RBC dynamics, spectral, temporal and spatial decompositions and model calculations were necessary to illuminate structure, dynamics, energetics and organisation. The final impression is that the structure of RBC in the PBL is given by a honeycomb-like arrangement of short-lived mixed-layer thermals with more passive downward motions in between. The regularity of the Cu-cloud cover results partly from the more stationary flow in the cloud-free cell centres. On the other hand it is shown that active as well as inactive clouds contribute to the cloud cover. Thus, the PBL flow and the cloud cover are decoupled, at least temporarily and locally.Due to sparse observational and measured information about RBC occurrence and structure in the PBL, additional material was gathered, resulting in the impression that RBC is one additional realised mode of organised convection in the PBL, as has already been clarified for PBL rolls and MCC by recent investigations.     

The composite shape and structure of coherent eddies in the convective boundary layer

Conditional sampling is used to locate coherent structures in a large data set obtained from flights by an instrumented light aircraft in convective boundary layers over Eyre Peninsula, South Australia. The high resolution and excellent spatial coverage of the data enable a detailed study of the internal structure of surface-layer plumes and mixed-layer thermals. A compositing technique is used to construct averaged traverses through coherent structures located within aircraft data runs of different altitudes and directions. Groups of composites are combined to form horizontal and vertical cross-sections which describe the internal flow patterns and the distribution of physical variables associated with typical coherent structures and their environment. In addition to the well-known along-wind features of surface-layer plumes, a strong, consistent inflow/entrainment pattern is evident in the lateral direction. Air from the horizontal plane channels around the sides and then in behind the microfront present at the upstream edge. Forces set up by the driving instability in the along-wind and vertical directions are counter-balanced by organised flow in the across-wind direction. It is found that mixed-layer thermal towers have a relatively simple form, consisting primarily of large columns of warm, upward-moving turbulent air, which may occasionally be in a state of slow rotation. An analysis of possible geometrical distortions within the results is performed, leading to a comparison ofw/z estimated from the horizontal velocity convergence field inside plumes/thermals, with that computed from the slope of partitionedw profiles.     

Numerical Study Of The Impact Of Coherent Structures On Vertical Transfers In The Atmospheric Boundary Layer

In two preceding papers, coherent structures of theatmospheric boundary layer (ABL), such as rollvortices or cells, were investigated through radar andaircraft observations collected during the TRAC-93(Turbulence Radar Aircraft Cells) experiment held inFrance in June 1993. The analysis of this experimentaldata set provided information on the spatialcharacteristics of these organisations (length scale,orientation, type ... ), their temporal and verticalevolution, and their relation with the dynamic andthermodynamic conditions of the ABL. For the thirdpaper in this series, a large eddy simulation model is used to examine the impact of thecoherent structures on the ABL vertical fluxes. Theanalysis of the simulated horizontal fields is madewith two-dimensional auto and cross-correlationsapplied on different pertinent ABL variables. Theresults emphasise a directional anisotropy of theseorganised fields throughout the ABL, much morepronounced in the heat flux fields, not only at thelength scale of organisations but also at theturbulence scales. This finding has an importantconsequence for traditional ABL flux measurementsbased on the hypothesis of isotropic and homogeneousturbulence. It can explain part of the underestimationof the surface fluxes often mentioned in theliterature. This approach makes it possible tomodify the concept of diffusion time (in chemicalmodelling) and could also lead to revised ABLparameterisations in Range Scale models.     

Boundary-layer heat and moisture budgets from fife

Aircraft stacks were flown upwind and downwind of the First ISLSCP Field Experiment (FIFE) site in Kansas to measure the heat and moisture budgets of the boundary layer under fairly clear skies for four daytime periods. In this paper, we evaluate the terms in the conservation equation. The vertical flux divergence and advection do not account for the difference between surface and low-level aircraft flux estimates. Budget estimates of the surface fluxes using the aircraft data agree well with surface flux measurements, but extrapolation of the aircraft fluxes gives surface fluxes that are too low. With the 5 km cutoff filter used, the aircraft underestimate for sensible heat flux is about 40%, and for the latent heat flux about 30%. Part of the underestimation is attributable to long-wavelength contributions (longer than the 5 km filter), but more investigation is needed.     

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.     

Observations of planview flux patterns within convective structures of the marine atmospheric surface layer

Air/sea flux variability on horizontal scales from 50 m to several km results, in part, from the presence of coherent convective structures within the atmospheric boundary layer. The horizontal distribution of fluxes within these convective updrafts and downdrafts is, therefore, central to studies of air/sea interaction and remote sensing of sea surface wind and wave fields. This study derives these flux patterns from observations of the Marine Atmospheric Surface Layer (MASL).Research aircraft flights through the MASL provide an optimal means for sampling large numbers of the above-mentioned coherent structures. The NCAR Electra flew numerous legs through the MASL at a height of 50 m during the 1987 stratocumulus phase of Project FIRE (First ISSCP (International Satellite Cloud Climatology Program) Regional Experiment).In situ measurements from these legs serve as the dataset for this paper. The data are processed in such a way as to retain only the turbulence fluctuations. Conditional sampling, based on the vertical velocity field, results in the isolation of convective updrafts and downdrafts. Compositing of the data for these two classes of convective drafts results in horizontal planviews of the vertical fluxes of buoyancy, absolute humidity, along-meanwind component of momentum, and vertical velocity. To ensure dynamical similarity, these horizontal planviews are oriented in a coordinate system aligned with the mean wind.     

Aircraft Observations of Sea-Surface Turbulent Fluxes Near the California Coast

Aircraft turbulence data from the Autonomous Ocean Sampling Network project were analyzed and compared to the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk parametrization of turbulent fluxes in an ocean area near the coast of California characterized by complex atmospheric flow. Turbulent fluxes measured at about 35 m above the sea surface using the eddy-correlation method were lower than bulk estimates under unstable and stable atmospheric stratification for all but light winds. Neutral turbulent transfer coefficients were used in this comparison because they remove the effects of mean atmospheric conditions and atmospheric stability. Spectral analysis suggested that kilometre-scale longitudinal rolls affect significantly turbulence measurements even near the sea surface, depending on sampling direction. Cross-wind sampling tended to capture all the available turbulent energy. Vertical soundings showed low boundary-layer depths and high flux divergence near the sea surface in the case of sensible heat flux but minimal flux divergence for the momentum flux. Cross-wind sampling and flux divergence were found to explain most of the observed discrepancies between the measured and bulk flux estimates. At low wind speeds the drag coefficient determined with eddy correlation and an inertial dissipation method after corrections were applied still showed high values compared to bulk estimates. This discrepancy correlated with the dominance of sea swell, which was a usually observed condition under low wind speeds. Under stable atmospheric conditions measured sensible heat fluxes, which usually have low values over the ocean, were possibly affected by measurement errors and deviated significantly from bulk estimates.     

Boundary-layer structure near the cold front of a marine cyclone during “ERICA”

The boundary layer in the warm sector of a moderately deepening winter cyclone during the Experiment on Rapidly Intensifying Cyclones over the Atlantic (ERICA) is studied near the cold front. Data from the National Center for Atmospheric Research Electra research aircraft are used to examine mean and turbulence quantities. The aircraft data and supplemental data from ships, drifting buoys and moored buoys reveal an equivalent-barotropic pressure field. The area is found to be dominated by gradients in temperature and in turbulent fluxes, with changes occurring over 100 km horizontally being comparable to changes over 350 m vertically. The horizontal components of the gradients are found to be a maximum in a direction perpendicular to the front. Cross-sections perpendicular to the front are used to illustrate boundary-layer structure. Profiles of wind speed, stress, wind direction and stress direction are estimated from an Ekman model that is modified to take into account the equivalent-barotropic pressure field. Comparison of profiles from the model to the aircraft-measured data show reasonable agreement far from the front (100 km) when the model uses a constant eddy viscosity of approximately 6 kg m^–1 s^–1. Near the front there is less agreement with the model. Profiles of turbulent fluxes of momentum, heat and latent heat are divergent, with along-wind momentum flux negative and decreasing upward, cross-wind momentum flux positive and increasing upward, and heat flux and latent heat flux small, positive and decreasing upward. Far from the front, the turbulent kinetic energy budget shows that dissipation balances shear production. However, near-front behavior has an imbalance at low altitude, with shear production appearing as a TKE sink.     

Very-Large-Scale Motions in the Atmospheric Boundary Layer Educed by Snapshot Proper Orthogonal Decomposition

Large-eddy simulations of the atmospheric boundary layer (ABL) under a wide range of stabilities are conducted to educe very-large-scale motions and then to study their dynamics and how they are influenced by buoyancy. Preliminary flow visualizations suggest that smaller-scale motions that resemble hairpins are embedded in much larger scale streamwise meandering rolls. Using simulations that represent more than 150 h of physical time, many snapshots in the \(xy\) -, \(yz\) - and \(xz\) -planes are then collected to perform snapshot proper orthogonal decomposition and further investigate the large structures. These analyses confirm that large streamwise rolls that share several features with the very-large-scale motions observed in laboratory studies arise as the dominant modes under most stabilities, but the effect of the surface kinematic buoyancy flux on the energy content of these dominant modes is very significant. The first two modes in the \(yz\) -plane in the neutral case contain up to 3 % of the total turbulent kinetic energy; they also have a vertical tilt angle in the \(yz\) -plane of about 0 to 30 \(^\circ \) due to the turning effect associated with the Coriolis force. Unstable cases also feature streamwise rolls, but in the convective ABL they are strengthened by rising plumes in between them, with two to four rolls spanning the whole domain in the first few modes; the Coriolis effect is much weaker in the unstable ABL. These rolls are no longer the dominant modes under stable conditions where the first mode is observed to contain sheet-like motions with high turbulent kinetic energy. Using these proper orthogonal decomposition modes, we are also able to extract the vertical velocity fields corresponding to individual modes and then to correlate them with the horizontal velocity or temperature fields to obtain the momentum and heat flux carried by individual modes. Structurally, the fluxes are explained by the topology of their corresponding modes. However, the fraction of the fluxes produced by the modes is invariably smaller than the fraction of energy they contain, particularly under stable conditions where the first modes are found to perform weak counter-gradient fluxes.