Exploring Stratocumulus Cloud-Top Entrainment Processes and Parameterizations by Using Doppler Cloud Radar Observations

Abstract Observations made at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site during uniform nonprecipitating stratocumulus cloud conditions for a 14-h period are used to examine cloud-top entrainment processes and parameterizations. The observations from a vertically pointing Doppler cloud radar provide estimates of vertical velocity variance and energy dissipation rate (EDR) terms in the parameterized turbulent kinetic energy (TKE) budget of the entrainment zone. Hourly averages of the vertical velocity variance term in the TKE entrainment formulation correlated strongly (r = 0.72) with the dissipation rate term in the entrainment zone, with an increased correlation (r = 0.92) when accounting for the nighttime decoupling of the boundary layer. Independent estimates of entrainment rates were obtained from an inversion-height budget using the local time derivative and horizontal advection of cloud-top height together with large-scale vertical velocity at the boundary layer inversion from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis model. The mean entrainment rate from the inversion-height budget during the 14-h period was 0.74 ± 0.15 cm s−1 and was used to calculate bulk coefficients for entrainment parameterizations based on convective velocity scale w* and TKE budgets of the entrainment zone. The hourly values of entrainment rates calculated using these coefficients exhibited good agreement with those calculated from the inversion-height budget associated with substantial changes in surface buoyancy production and cloud-top radiative cooling. The results indicate a strong potential for making entrainment rate estimates directly from radar vertical velocity variance and the EDR measurements.

Download Full-text

Buoyancy effects on large-scale motions in convective atmospheric boundary layers: implications for modulation of near-wall processes

Journal of Fluid Mechanics ◽

10.1017/jfm.2018.711 ◽

2018 ◽

Vol 856 ◽

pp. 135-168 ◽

Cited By ~ 19

Author(s):

S. T. Salesky ◽

W. Anderson

Keyword(s):

Boundary Layer ◽

Velocity Field ◽

Boundary Layers ◽

Vertical Velocity ◽

Amplitude Modulation ◽

Large Scale ◽

Streamwise Velocity ◽

Small Scale ◽

Convective Conditions ◽

Wide Range

A number of recent studies have demonstrated the existence of so-called large- and very-large-scale motions (LSM, VLSM) that occur in the logarithmic region of inertia-dominated wall-bounded turbulent flows. These regions exhibit significant streamwise coherence, and have been shown to modulate the amplitude and frequency of small-scale inner-layer fluctuations in smooth-wall turbulent boundary layers. In contrast, the extent to which analogous modulation occurs in inertia-dominated flows subjected to convective thermal stratification (low Richardson number) and Coriolis forcing (low Rossby number), has not been considered. And yet, these parameter values encompass a wide range of important environmental flows. In this article, we present evidence of amplitude modulation (AM) phenomena in the unstably stratified (i.e. convective) atmospheric boundary layer, and link changes in AM to changes in the topology of coherent structures with increasing instability. We perform a suite of large eddy simulations spanning weakly ($-z_{i}/L=3.1$) to highly convective ($-z_{i}/L=1082$) conditions (where$-z_{i}/L$is the bulk stability parameter formed from the boundary-layer depth$z_{i}$and the Obukhov length $L$) to investigate how AM is affected by buoyancy. Results demonstrate that as unstable stratification increases, the inclination angle of surface layer structures (as determined from the two-point correlation of streamwise velocity) increases from$\unicode[STIX]{x1D6FE}\approx 15^{\circ }$for weakly convective conditions to nearly vertical for highly convective conditions. As$-z_{i}/L$increases, LSMs in the streamwise velocity field transition from long, linear updrafts (or horizontal convective rolls) to open cellular patterns, analogous to turbulent Rayleigh–Bénard convection. These changes in the instantaneous velocity field are accompanied by a shift in the outer peak in the streamwise and vertical velocity spectra to smaller dimensionless wavelengths until the energy is concentrated at a single peak. The decoupling procedure proposed by Mathiset al.(J. Fluid Mech., vol. 628, 2009a, pp. 311–337) is used to investigate the extent to which amplitude modulation of small-scale turbulence occurs due to large-scale streamwise and vertical velocity fluctuations. As the spatial attributes of flow structures change from streamwise to vertically dominated, modulation by the large-scale streamwise velocity decreases monotonically. However, the modulating influence of the large-scale vertical velocity remains significant across the stability range considered. We report, finally, that amplitude modulation correlations are insensitive to the computational mesh resolution for flows forced by shear, buoyancy and Coriolis accelerations.

Download Full-text

Observed Boundary Layer Controls on Shallow Cumulus at the ARM Southern Great Plains Site

Journal of the Atmospheric Sciences ◽

10.1175/jas-d-17-0244.1 ◽

2018 ◽

Vol 75 (7) ◽

pp. 2235-2255 ◽

Cited By ~ 15

Author(s):

Neil P. Lareau ◽

Yunyan Zhang ◽

Stephen A. Klein

Keyword(s):

Boundary Layer ◽

Great Plains ◽

Vertical Velocity ◽

Mass Flux ◽

Cloud Base ◽

Liquid Water Path ◽

Southern Great Plains ◽

Clear Sky ◽

Shallow Cumulus ◽

In Situ Sensors

Abstract The boundary layer controls on shallow cumulus (ShCu) convection are examined using a suite of remote and in situ sensors at ARM Southern Great Plains (SGP). A key instrument in the study is a Doppler lidar that measures vertical velocity in the CBL and along cloud base. Using a sample of 138 ShCu days, the composite structure of the ShCu CBL is examined, revealing increased vertical velocity (VV) variance during periods of medium cloud cover and higher VV skewness on ShCu days than on clear-sky days. The subcloud circulations of 1791 individual cumuli are also examined. From these data, we show that cloud-base updrafts, normalized by convective velocity, vary as a function of updraft width normalized by CBL depth. It is also found that 63% of clouds have positive cloud-base mass flux and are linked to coherent updrafts extending over the depth of the CBL. In contrast, negative mass flux clouds lack coherent subcloud updrafts. Both sets of clouds possess narrow downdrafts extending from the cloud edges into the subcloud layer. These downdrafts are also present adjacent to cloud-free updrafts, suggesting they are mechanical in origin. The cloud-base updraft data are subsequently combined with observations of convective inhibition to form dimensionless “cloud inhibition” (CI) parameters. Updraft fraction and liquid water path are shown to vary inversely with CI, a finding consistent with CIN-based closures used in convective parameterizations. However, we also demonstrate a limited link between CBL vertical velocity variance and cloud-base updrafts, suggesting that additional factors, including updraft width, are necessary predictors for cloud-base updrafts.

Download Full-text

How large-scale subsidence affects stratocumulus transitions

Atmospheric Chemistry and Physics Discussions ◽

10.5194/acpd-15-17229-2015 ◽

2015 ◽

Vol 15 (12) ◽

pp. 17229-17250 ◽

Cited By ~ 1

Author(s):

J. J. van der Dussen ◽

S. R. de Roode ◽

A. P. Siebesma

Keyword(s):

Boundary Layer ◽

Large Scale ◽

Climate Modeling ◽

Hadley Circulation ◽

Cloud Base ◽

Entrainment Rate ◽

Liquid Water Path ◽

Trade Winds ◽

Enhanced Absorption ◽

Subsidence Rate

Abstract. Some climate modeling results suggest that the Hadley circulation might weaken in a future climate, causing a subsequent reduction in the large-scale subsidence velocity in the subtropics. In this study we analyze the cloud liquid water path (LWP) budget from large-eddy simulation (LES) results of three idealized stratocumulus transition cases each with a different subsidence rate. As shown in previous studies a reduced subsidence is found to lead to a deeper stratocumulus-topped boundary layer, an enhanced cloud-top entrainment rate and a delay in the transition of stratocumulus clouds into shallow cumulus clouds during its equatorwards advection by the prevailing trade winds. The effect of a reduction of the subsidence rate can be summarized as follows. The initial deepening of the stratocumulus layer is partly counteracted by an enhanced absorption of solar radiation. After some hours the deepening of the boundary layer is accelerated by an enhancement of the entrainment rate. Because this is accompanied by a change in the cloud-base turbulent fluxes of moisture and heat, the net change in the LWP due to changes in the turbulent flux profiles is negligibly small.

Download Full-text

How is the marine atmospheric boundary layer turbulence organized in the trades ?

10.5194/egusphere-egu21-12339 ◽

2021 ◽

Author(s):

Pierre-Etienne Brilouet ◽

Marie Lothon ◽

Sandrine Bony

Keyword(s):

Boundary Layer ◽

Atmospheric Boundary Layer ◽

Vertical Velocity ◽

Large Scale ◽

Layer Structure ◽

Potential Temperature ◽

Surface Wind ◽

Cloud Amount ◽

Velocity Spectrum ◽

Marine Atmospheric Boundary Layer

Tradewind clouds can exhibit a wide diversity of mesoscale organizations, and the turbulence of marine atmospheric boundary layer (MABL) can exhibit coherent structures and mesoscale circulations. One of the objectives of the EUREC4A (Elucidating the role of cloud-circulation coupling in climate) field experiment was to better understand the tight interplay between the mesoscale organization of clouds, boundary-layer processes, and the large-scale environment.During the experiment, that took place East of Barbados over the Western Tropical Atlantic Ocean in Jan-Feb 2020, the French ATR-42 research aircraft was devoted to the characterization of the cloud amount and of the subcoud layer structure. During its 17 research flights, it sampled a large diversity of large scale conditions and cloud patterns. Multiple sensors onboard the aircraft measured high-frequency fluctuations of potential temperature, water vapour mixing ratio and wind , allowing for an extensive characterization of the turbulence within the subcloud layer. A quality-controled and calibrated turbulence dataset was produced on the basis of these measurements, which is now available on the EUREC4A AERIS data portal.The MABL turbulent structure is studied using this dataset, through a spectral analysis of the vertical velocity. Vertical profiles of characteristic length scales reveal a non-isotropic structure with a stretching of the eddies along the mean wind. The organization strength of the turbulent field is also explored by defining a diagnostic based on the shape of the vertical velocity spectrum. The structure and the degree of organization of the subcloud layer are characterized for different types of mesoscale convective pattern and as a function of the large-scale environment, including near-surface wind and lower-tropospheric stability conditions.&#160;

Download Full-text

Vertical velocity and turbulence aspects during Mistral events as observed by UHF wind profilers

Annales Geophysicae ◽

10.5194/angeo-22-3927-2004 ◽

2004 ◽

Vol 22 (11) ◽

pp. 3927-3936 ◽

Cited By ~ 12

Author(s):

J.-L. Caccia ◽

V. Guénard ◽

B. Benech ◽

B. Campistron ◽

P. Drobinski

Keyword(s):

Boundary Layer ◽

Planetary Boundary Layer ◽

Vertical Velocity ◽

Convective Boundary Layer ◽

Dissipation Rate ◽

Convective Boundary ◽

Turbulent Dissipation ◽

The Alps ◽

Rhone Valley ◽

Wind Profilers

Abstract. The general purpose of this paper is to experimentally study mesoscale dynamical aspects of the Mistral in the coastal area located at the exit of the Rhône-valley. The Mistral is a northerly low-level flow blowing in southern France along the Rhône-valley axis, located between the French Alps and the Massif Central, towards the Mediterranean Sea. The experimental data are obtained by UHF wind profilers deployed during two major field campaigns, MAP (Mesoscale Alpine Program) in autumn 1999, and ESCOMPTE (Expérience sur Site pour COntraindre les Modèles de Pollution atmosphériques et de Transports d'Emission) in summer 2001. Thanks to the use of the time evolution of the vertical profile of the horizontal wind vector, recent works have shown that the dynamics of the Mistral is highly dependent on the season because of the occurrence of specific synoptic patterns. In addition, during summer, thermal forcing leads to a combination of sea breeze with Mistral and weaker Mistral due to the enhanced friction while, during autumn, absence of convective turbulence leads to substantial acceleration as low-level jets are generated in the stably stratified planetary boundary layer. At the exit of the Rhône valley, the gap flow dynamics dominates, whereas at the lee of the Alps, the dynamics is driven by the relative contribution of "flow around" and "flow over" mechanisms, upstream of the Alps. This paper analyses vertical velocity and turbulence, i.e. turbulent dissipation rate, with data obtained by the same UHF wind profilers during the same Mistral events. In autumn, the motions are found to be globally and significantly subsident, which is coherent for a dry, cold and stable flow approaching the sea, and the turbulence is found to be of pure dynamical origin (wind shears and mountain/lee wave breaking), which is coherent with non-convective situations. In summer, due to the ground heating and to the interactions with thermal circulation, the vertical motions are less pronounced and no longer have systematic subsident charateristics. In addition, those vertical motions are found to be much less developed during the nighttimes because of the stabilization of the nocturnal planetary boundary layer due to a ground cooling. The enhanced turbulent dissipation-rate values found at lower levels during the afternoons of weak Mistral cases are consistent with the installation of the summer convective boundary layer and show that, as expected in weaker Mistral events, the convection is the preponderant factor for the turbulence generation. On the other hand, for stronger cases, such a convective boundary layer installation is perturbed by the Mistral.

Download Full-text

Estimating the Shallow Convective Mass Flux from the Subcloud-Layer Mass Budget

Journal of the Atmospheric Sciences ◽

10.1175/jas-d-19-0135.1 ◽

2020 ◽

Vol 77 (5) ◽

pp. 1559-1574 ◽

Cited By ~ 1

Author(s):

Raphaela Vogel ◽

Sandrine Bony ◽

Bjorn Stevens

Keyword(s):

Vertical Velocity ◽

Diurnal Cycle ◽

Mass Flux ◽

Large Scale ◽

Cloud Base ◽

Entrainment Rate ◽

Field Campaign ◽

Mass Budget ◽

Convective Mixing ◽

Convective Mass Flux

Abstract This paper develops a method to estimate the shallow-convective mass flux M at the top of the subcloud layer as a residual of the subcloud-layer mass budget. The ability of the mass-budget estimate to reproduce the mass flux diagnosed directly from the cloud-core area fraction and vertical velocity is tested using real-case large-eddy simulations over the tropical Atlantic. We find that M reproduces well the magnitude, diurnal cycle, and day-to-day variability of the core-sampled mass flux, with an average root-mean-square error of less than 30% of the mean. The average M across the four winter days analyzed is 12 mm s−1, where the entrainment rate E contributes on average 14 mm s−1 and the large-scale vertical velocity W contributes −2 mm s−1. We find that day-to-day variations in M are mostly explained by variations in W, whereas E is very similar among the different days analyzed. Instead E exhibits a pronounced diurnal cycle, with a minimum of about 10 mm s−1 around sunset and a maximum of about 18 mm s−1 around sunrise. Application of the method to dropsonde data from an airborne field campaign in August 2016 yields the first measurements of the mass flux derived from the mass budget, and supports the result that the variability in M is mostly due to the variability in W. Our analyses thus suggest a strong coupling between the day-to-day variability in shallow convective mixing (as measured by M) and the large-scale circulation (as measured by W). Application of the method to the EUREC4A field campaign will help evaluate this coupling, and assess its implications for cloud-base cloudiness.

Download Full-text

Bulk Models of the Sheared Convective Boundary Layer: Evaluation through Large Eddy Simulations

Journal of the Atmospheric Sciences ◽

10.1175/jas3870.1 ◽

2007 ◽

Vol 64 (3) ◽

pp. 786-807 ◽

Cited By ~ 37

Author(s):

Robert Conzemius ◽

Evgeni Fedorovich

Keyword(s):

Boundary Layer ◽

Convective Boundary Layer ◽

Entrainment Rate ◽

Convective Boundary ◽

Surface Shear ◽

Large Eddy ◽

Entrainment Zone ◽

Sheared Convective Boundary Layer ◽

Bulk Model ◽

Approach Zero

Abstract A set of first-order model (FOM) equations, describing the sheared convective boundary layer (CBL) evolution, is derived. The model output is compared with predictions of the zero-order bulk model (ZOM) for the same CBL type. Large eddy simulation (LES) data are employed to test both models. The results show an advantage of the FOM over the ZOM in the prediction of entrainment, but in many CBL cases, the predictions by the two models are fairly close. Despite its relative simplicity, the ZOM is able to quantify the effects of shear production and dissipation in an integral sense—as long as the constants describing the integral dissipation of shear- and buoyancy-produced turbulence kinetic energy (TKE) are prescribed appropriately and the shear is weak enough that the denominator of the ZOM entrainment equation does not approach zero, causing a numerical instability in the solutions. Overall, the FOM better predicts the entrainment rate due to its ability to avoid this instability. Also, the FOM in a more physically consistent manner reproduces the sheared CBL entrainment zone, whose depth is controlled by a balance among shear generation, buoyancy consumption, and dissipation of TKE. Such balance is manifested by nearly constant values of Richardson numbers observed in the entrainment zone of simulated sheared CBLs. Conducted model tests support the conclusion that the surface shear generation of TKE and its corresponding dissipation, as well as the nonstationary terms, can be omitted from the integral TKE balance equation.

Download Full-text

Turbulence in the stratified boundary layer under ice: observations from Lake Baikal and a new similarity model

Hydrology and Earth System Sciences ◽

10.5194/hess-24-1691-2020 ◽

2020 ◽

Vol 24 (4) ◽

pp. 1691-1708 ◽

Cited By ~ 2

Author(s):

Georgiy Kirillin ◽

Ilya Aslamov ◽

Vladimir Kozlov ◽

Roman Zdorovennov ◽

Nikolai Granin

Keyword(s):

Boundary Layer ◽

Heat Flux ◽

Lake Baikal ◽

Dissipation Rate ◽

Large Scale ◽

Shear Velocity ◽

Ice Cover ◽

Solid Boundary ◽

Buoyancy Frequency ◽

Ice Water

Abstract. Seasonal ice cover on lakes and polar seas creates seasonally developing boundary layer at the ice base with specific features: fixed temperature at the solid boundary and stable density stratification beneath. Turbulent transport in the boundary layer determines the ice growth and melting conditions at the ice–water interface, especially in large lakes and marginal seas, where large-scale water circulation can produce highly variable mixing conditions. Since the boundary mixing under ice is difficult to measure, existing models of ice cover dynamics usually neglect or parameterize it in a very simplistic form. We present the first detailed observations on mixing under ice of Lake Baikal, obtained with the help of advanced acoustic methods. The dissipation rate of the turbulent kinetic energy (TKE) was derived from correlations (structure functions) of current velocities within the boundary layer. The range of the dissipation rate variability covered 2 orders of magnitude, demonstrating strongly turbulent conditions. Intensity of mixing was closely connected to the mean speeds of the large-scale under-ice currents. Mixing developed on the background of stable density (temperature) stratification, which affected the vertical structure of the boundary layer. To account for stratification effects, we propose a model of the turbulent energy budget based on the length scale incorporating the dissipation rate and the buoyancy frequency (Dougherty–Ozmidov scaling). The model agrees well with the observations and yields a scaling relationship for the ice–water heat flux as a function of the shear velocity squared. The ice–water heat fluxes in the field were the largest among all reported in lakes (up to 40 W m−2) and scaled well against the proposed relationship. The ultimate finding is that of a strong dependence of the water–ice heat flux on the shear velocity under ice. The result suggests large errors in the heat flux estimations when the traditional “bulk” approach is applied to stratified boundary layers. It also implies that under-ice currents may have much stronger effect on the ice melt than estimated by traditional models.

Download Full-text

Turbulence Process Domination under the Combined Forcings of Wind Stress, the Langmuir Vortex Force, and Surface Cooling

Journal of Physical Oceanography ◽

10.1175/jpo-d-13-021.1 ◽

2014 ◽

Vol 44 (1) ◽

pp. 44-67 ◽

Cited By ~ 21

Author(s):

A. E. Gargett ◽

C. E. Grosch

Keyword(s):

Boundary Layer ◽

Wind Stress ◽

Vertical Velocity ◽

Ocean Surface ◽

The Other ◽

Surface Boundary ◽

Langmuir Circulation ◽

Surface Cooling ◽

Wide Range ◽

Velocity Variance

Abstract Turbulence in the ocean surface layer is generated by time-varying combinations of destabilizing surface buoyancy flux, wind stress forcing, and wave forcing through a vortex force associated with the surface wave field. Observations of time- and depth-averaged vertical velocity variance of full-depth turbulence in shallow unstratified water columns under destabilizing buoyancy forcing are used to determine when process domination can be assigned over a wide range of mixed forcings. The properties of two turbulence archetypes, one representing full-depth Langmuir circulations and the other representing full-depth convection, are described in detail. It is demonstrated that these archetypes lie in distinct regions of the plane of , where and are Langmuir and Rayleigh numbers, respectively, derived from scaling with surface stress velocity and a time scale characteristic of the growth of Langmuir circulation , where and are mean and Stokes velocities, respectively. Situations in which neither process dominates lie between the two end members, with relative dominance given by proximity to one or the other. Cases dominated by direct stress forcing are conspicuous by their absence. In cases of Langmuir domination, surface Stokes velocity is linearly related to , making it impossible to differentiate between scaling depth-averaged vertical velocity variance with , and any other scaling involving both and . A third nondimensional parameter is introduced and used to assess the importance of bottom boundary layer turbulence in a depth-limited system. Questions of time dependence and applicability of results to the open ocean surface boundary layer are considered.

Download Full-text

Vertical Velocity Variance Measurements from Wind Profiling Radars

10.5194/amt-2016-299 ◽

2016 ◽

Cited By ~ 1

Author(s):

Katherine McCaffrey ◽

Laura Bianco ◽

Paul Johnston ◽

James M. Wilczak

Keyword(s):

Boundary Layer ◽

Time Series ◽

Planetary Boundary Layer ◽

Vertical Velocity ◽

Prediction Models ◽

Weather Prediction ◽

Small Scale ◽

Stable Conditions ◽

Sonic Anemometers ◽

Velocity Variance

Abstract. Observations of turbulence in the planetary boundary layer are critical for developing and evaluating boundary layer parameterizations in mesoscale numerical weather prediction models. These observations, however, are expensive, and rarely profile the entire boundary layer. Using optimized configurations for 449 MHz and 915 MHz wind profiling radars during the eXperimental Planetary boundary layer Instrumentation Assessment, improvements have been made to the historical methods of measuring vertical velocity variance through the time series of vertical velocity, as well as the Doppler spectral width. Using six heights of sonic anemometers mounted on a 300-m tower, correlations of up to R2 = 0.74 are seen in measurements of the large-scale variances from the radar time series, and R2 = 0.79 in measurements of small-scale variance from radar spectral widths. The total variance, measured as the sum of the small- and large-scales agrees well with sonic anemometers, with R2 = 0.79. Correlation is higher in daytime, convective boundary layers than nighttime, stable conditions when turbulence levels are smaller. With the good agreement with the in situ measurements, highly-resolved profiles up to 2 km can be accurately observed from the 449 MHz radar, and 1 km from the 915 MHz radar. This optimized configuration will provide unique observations for the verification and improvement to boundary layer parameterizations in mesoscale models.

Download Full-text