Entrainment, Detrainment, and Dilution of Dry and Moist Atmospheric Thermals
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
There is evidence that most convective detrainment occurs below 355-K potential temperature (14 km). On the other hand, it is also generally accepted that convection does occasionally reach the altitude of the tropopause, typically approx. 16 km, or higher. The question is whether the relatively infrequent penetrations above 14 km are sufficient to influence the composition of the atmosphere above 14 km. Using measurements of O3 and CO as tracers of convective activity, I quantify the height distribution of mass detraining from convection near the tropopause by means of a simple mixing model similar to those used to calculate mixing into the stratospheric "tropical pipe". I calculate that significant detrainment from convection is occurring as high as approx. 380 K (16 km), an altitude typical of the tropopause. Further, I calculate that as much as 40% of the mass crossing the 380-K surface at this time detrained above 370 K. I conclude that the "mixing barrier" at 14 km is more porous than previously thought.
Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Further, comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.
Partitioning deep convective cloud condensates into components that sediment and detrain, known to be a challenge for global climate models, is important for cloud vertical distribution and anvil cloud formation. In this study, we address this issue by improving the convective microphysics scheme in the National Center for Atmospheric Research Community Atmosphere Model version 5.3 (CAM5.3). The improvements include: (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. Results from the default and modified CAM5.3 models were evaluated against observations from the U.S. Department of Energy Tropical Warm Pool-International Cloud Experiment (TWP-ICE) field campaign. The default model overestimates ice amount, which is largely attributed to the underestimation of convective ice particle sedimentation. By considering cloud ice sedimentation and rimed ice particles and applying a new convective snow terminal velocity parameterization, the vertical distribution of ice amount is much improved in the midtroposphere and upper troposphere when compared to observations. The vertical distribution of ice condensate also agrees well with observational best estimates upon considering snow detrainment. Comparison with observed convective updrafts reveals that current bulk model fails to reproduce the observed updraft magnitude and occurrence frequency, suggesting spectral distributions be required to simulate the subgrid updraft heterogeneity.
The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.
Abstract Vertical wind shear has long been known to tilt convective towers and reduce thermal ascent rates. The purpose of this study is to better understand the physical mechanisms responsible for reduced ascent rates in shallow convection. In particular, the study focuses on cloud-edge mass flux to assess how shear impacts mass-flux profiles of both the ensemble and individual clouds of various depths. A compositing algorithm is used to distill large-eddy simulation (LES) output to focus on up- and down-shear cloud edges that are not influenced by complex cloud geometry or nearby clouds. A direct entrainment algorithm is used to estimate the mass flux through the cloud surface. We find that the dynamics on the up- and down-shear sides are fundamentally different, with the entrainment of environmental momentum and dilution of buoyancy being primarily responsible for the reduced down-shear ascent rates. Direct estimates of fluid flow through the cloud interface indicate a counter-shear organized flow pattern that entrains on the down-shear side and detrains on the up-shear side, resulting from the subcloud shear being lifted into the cloud layer by the updraft. In spite of organized regions of entrainment and detrainment, the overall net lateral mass flux remains unchanged with respect to the no shear run, with weak detrainment present throughout cloud depth.
Climate models often ignore the influence of ice-phase physics (IPP) of hydrometeors as a second order effect. This has also been true for McRAS (Microphysics of clouds with Relaxed Arakawa Schubert Scheme) developed by the authors. Recognizing that the temperature sounding is critical for moist-convection, and, that IPP would modify it, we investigated the influence of introducing IPP into McRAS coupled to FvGCM (finite volume General Circulation Model with NCAR physics). We analyzed three 3-yr long simulations; the first called Control Case, CC and had no IPP; the other two called Experiments El and E2 had IPP introduced with two different in-cloud freezing assumptions. Simulation El assumed that all hydrometeors remain liquid in the updraft and freeze upon detrainment. Simulation E2 invoked the in-cloud freezing of new condensate generated at subfreezing temperatures in the updraft while old cloud water continued to ascend as liquid. Upon detrainment, this cloud water also froze like in E1. With these assumptions, about 50% of hydrometeors froze in the tower and the rest froze in the anvil. However, in both El and E2, the frozen hydrometeors melted during fall at the first encounter of above freezing ambient temperature. Comparative analysis revealed that El simulated far more mid-level and far less deep clouds while E2 had modified deep and more mid-level clouds as compared to CC along with some major changes around the melt-level. We infer that IPP produced a more realistic response in E2. At the basic level, the results show that ice-phase processes influence convective detrainment at mid- and deep levels in accord with TOGAGOARE observations. The results suggest that IPP can help to mitigate less-than-observed mid-level and over-abundance of deep convective clouds in McRAS.
The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.
Here, we introduce a physics-based aerosol wet removal scheme with unified treatments of aerosol transport and removal by convective clouds into the Community Atmosphere Model version 6. Since several important physical processes are still neglected or poorly represented in this new physics-based scheme, we develop secondary improvements to the parameterizations of aerosol activation, resuspension, and cloud-borne aerosol detrainment in this new scheme. Changes in the aerosol wet removal scheme cause tropospheric aerosol concentrations to decrease to different extents: compared to the control run, the physics-based scheme significantly decreases aerosol burdens by up to 60% over the southern Pacific Ocean, whereas the secondary improvements mitigate the decreasing tendency. The burden changes also depend on aerosol chemical components: the sulfate mass decrease is compensated by secondary production, black carbon (BC) is effectively removed via increasing the hygroscopicity of particulate organic matter from 0 to 0.2, and dust shows the most spatially heterogeneous changes. Simulated aerosol profiles are evaluated against aircraft-based observations over the Pacific and Atlantic Oceans. The secondary-improved scheme reduces the overestimations of upper tropospheric BC and sea salt concentrations by a factor of 10 and 1,000, respectively, and reproduces the dependence of BC mass decrease rates on cloud types. Consideration of convective cloud-borne aerosol detrainment plays the most important role in enhancing the aerosol wet removal and decreasing the positive biases of tropospheric BC and sea salt concentrations. We also summarize unresolved issues related to convective cloud genesis and microphysics, cloud-borne aerosol evolution, and BC and dust emissions.
This article presents a new deep convective parameterization that determines cloud characteristics based on a specified cloud size distribution. The vertical profiles of cloud properties are determined by analytic equations, which formulate entrainment with an inverse relationship to cloud width. In line with recent studies of large eddy simulations (LES), cloud widths are assumed to be constant with height and vertical mass flux (M) characteristics of the clouds are therefore regulated by the vertical velocity profile. The parameterization is configured to work with existing cloud base M closure formulations, with the closure predicting the total cloud area rather than the cloud base M directly. Analytic formula are also used to connect the vertical wind shear magnitude to the cloud size distribution, wherein larger shear magnitudes result in more numerous large updrafts than weaker shear magnitudes, which is in line with recent research results. The parameterization is compared against 10 deep convective LES with varying thermodynamic and vertical wind shear profiles. Results show dramatic improvements in the prediction of normalized M, detrainment, and the properties of detrained air over the existing Zhang and McFarlane (1995) scheme. In particular, the new model is able to correctly portray the transition from a bottom-heavy M profile in weakly sheared environments, to a top-heavy M profile in strongly sheared environments.
Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.
Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .
Abstract Interfaces at the edge of an idealized, nonprecipitating, warm cloud are studied using direct numerical simulation (DNS) complemented with a Lagrangian particle tracking routine. Once a shell has formed, four zones can be distinguished: the cloud core, visible shell, invisible shell, and the environment. The union of the visible and invisible regions is the shell commonly referred to in literature. The boundary between the invisible shell and the environment is the turbulent–nonturbulent interface (TNTI), which is typically not considered in cloud studies. Three million particles were seeded homogeneously across the domain and properties were recorded along individual trajectories. The results demonstrate that the traditional cloud boundary (separating cloudy and noncloudy regions using thresholds applied on liquid condensate or updraft velocity) are some distance away from the TNTI. Furthermore, there is no dynamic difference between the traditional liquid-condensate boundary and the region extending to the TNTI. However, particles crossing the TNTI exhibit a sharp jump in enstrophy and a smooth increase in buoyancy. The traditional cloud boundary coincides with the location of minimum buoyancy in the shell. The shell premixes the entraining and detraining air and analysis reveals a highly skewed picture of entrainment and detrainment at the traditional cloud boundary. A preferential entrainment of particles with velocity and specific humidity higher than the mean values in the shell is observed. Large-eddy simulation of a more realistic setup detects an interface with similar properties using the same thresholds as in the DNS, indicating that the DNS results extrapolate beyond their idealized conditions.
A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.
A global climatology of stratospheric aerosol is created by combining nearly a decade (1979-1981 and 1984-1990) of contemporaneous observations from the Stratospheric Aerosol and Gas Experiment (SAGE I and II) and Stratospheric Aerosol Measurement (SAM II) instruments. One goal of this work is to provide a representative distribution of the aerosol layer for use in radiative and chemical modeling. A table of decadal average 1 micron extinction values is included, extending from the tropopause to 35 km and 80 deg S to 85 deg N, which allows estimation of surface area density. We find that the aerosol layer is distinctly volcanic in nature and suggest that the decadal average is a more useful estimate of future aerosol loading than a 'background' loading, which is never clearly achieved during the data record. This climatology lends insight into the general circulation of the stratosphere. Latitude - altitude sections of extinction radio at 1 micron are shown, average by decade, season, and phase of the quasi-biennial oscillation (QBO). A tropical reservoir region is diagnosed, with an 'upper' and a 'lower' transport regime. In the tropics above 22 km (upper regime), enhanced lofting occurs in the summer, with suppressed lofting or eddy dilution in the winter. In the extratropics within two scale heights of the tropopause (lower regime), poleward and downward transport is most robust during winter, especially in the northern hemisphere. The transport patterns persist into the subsequent equinoctial season. Ascent associated with QBO easterly shear favors detrainment in the upper regime, while relative descent and poleward spreading during QBO westerly shear favors detrainment in the lower regime. Extinction radio differences between the winter-spring and summer-fall hemispheres, and differences between the two phases of the QBO, are typically 20-50%. Dynamical implications of the aerosol distributions are explored, with focus on interhemispheric differences, strong subtropical gradients, and the pronounced annual cycle.
The design of a prognostic cloud scheme named McRAS (Microphysics of clouds with Relaxed Arakawa-Schubert Scheme) for general circulation models (GCMs) will be discussed. McRAS distinguishes three types of clouds: (1) convective, (2) stratiform, and (3) boundary-layer types. The convective clouds transform and merge into stratiform clouds on an hourly time-scale, while the boundary-layer clouds merge into the stratiform clouds instantly. The cloud condensate converts into precipitation following the auto-conversion equations of Sundqvist that contain a parametric adaptation for the Bergeron-Findeisen process of ice crystal growth and collection of cloud condensate by precipitation. All clouds convect, advect, as well as diffuse both horizontally and vertically with a fully interactive cloud-microphysics throughout the life-cycle of the cloud, while the optical properties of clouds are derived from the statistical distribution of hydrometeors and idealized cloud geometry. An evaluation of McRAS in a single column model (SCM) with the GATE Phase III and 5-ARN CART datasets has shown that together with the rest of the model physics, McRAS can simulate the observed temperature, humidity, and precipitation without many systematic errors. The time history and time mean incloud water and ice distribution, fractional cloudiness, cloud optical thickness, origin of precipitation in the convective anvil and towers, and the convective updraft and downdraft velocities and mass fluxes all show a realistic behavior. Performance of McRAS in GEOS 11 GCM shows several satisfactory features but some of the remaining deficiencies suggest need for additional research involving convective triggers and inhibitors, provision for continuously detraining updraft, a realistic scheme for cumulus gravity wave drag, and refinements to physical conditions for ascertaining cloud detrainment level.
Compared to satellite-derived heating profiles, the Goddard Institute for Space Studies general circulation model (GCM) convective heating is too deep and its stratiform upper-level heating is too weak. This deficiency highlights the need for GCMs to parameterize the mesoscale organization of convection. Cloud-resolving model simulations of convection near Darwin, Australia, in weak wind shear environments of different humidities are used to characterize mesoscale organization processes and to provide parameterization guidance. Downdraft cold pools appear to stimulate further deep convection both through their effect on eddy size and vertical velocity. Anomalously humid air surrounds updrafts, reducing the efficacy of entrainment. Recovery of cold pool properties to ambient conditions over 5-6 h proceeds differently over land and ocean. Over ocean increased surface fluxes restore the cold pool to prestorm conditions. Over land surface fluxes are suppressed in the cold pool region; temperature decreases and humidity increases, and both then remain nearly constant, while the undisturbed environment cools diurnally. The upper-troposphere stratiform rain region area lags convection by 5-6 h under humid active monsoon conditions but by only 1-2 h during drier break periods, suggesting that mesoscale organization is more readily sustained in a humid environment. Stratiform region hydrometeor mixing ratio lags convection by 0-2 h, suggesting that it is strongly influenced by detrainment from convective updrafts. Small stratiform region temperature anomalies suggest that a mesoscale updraft parameterization initialized with properties of buoyant detrained air and evolving to a balance between diabatic heating and adiabatic cooling might be a plausible approach for GCMs.
It is well known that stratospheric humidity is primarily controlled by freeze drying (ice crystal growth and sedimentation) of air ascending across the cold tropical tropopause. However, the suggestion of an important source of water vapor from deep convection that extends above the tropical tropopause has persisted. There exists much anecdotal evidence of direct convective hydration of the lower stratosphere based on measurements from high altitude aircraft campaigns, but quantifying the impact of deep convection on the overall budget of stratospheric water vapor has proven challenging. The role of convection on the humidity of the upper troposphere and lower stratosphere (UTLS) is investigated in simulations of cirrus clouds along trajectories launched from given potential temperature level surfaces. The one-dimensional (vertical) cloud model tracks individual ice crystals through their life cycles, beginning with nucleation or detrainment from convection, followed by deposition growth, sedimentation and sublimation. Convective influence of the parcels is diagnosed by tracing the trajectories through time-dependent fields of convective cloud-top height adjusted to match the CloudSAT and CALIPSO statistics. Model simulations of UTLS water vapor and cloud fields are evaluated and constrained by comparison with MLS and CALIPSO measurements. The simulation results indicate that the overall impact of convection on water vapor near the tropical tropopause is 10-15%, while the impact on lower stratospheric humidity is no more than a few percent. Ice crystals detrained from deep convection have relatively small effect. The general implications for the importance of deep convection on UTLS humidity and cirrus cloud fraction will be discussed.