Search

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm−3 to represent clean conditions, or as 2000 cm−3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 643 m3 and the turbulent kinetic energy dissipation rate of 10−3 m−2 s−3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s−1, are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Advanced microphysics schemes (such as Eulerian bin and Lagrangian superdroplet) are becoming standard tools for cloud physics research and parameterization development. This study compares a double-moment bin scheme and a Lagrangian superdroplet scheme via large-eddy simulations of nonprecipitating and precipitating cumulus congestus clouds. Cloud water mixing ratio in the bin simulations is reduced compared to the Lagrangian simulations in the upper part of the cloud, likely from numerical diffusion, which is absent in the Lagrangian approach. Greater diffusion in the bin simulations is compensated by more secondary droplet activation (activation above cloud base), leading to similar or somewhat higher droplet number concentrations and smaller mean droplet radius than the Lagrangian simulations for the nonprecipitating case. The bin scheme also produces a significantly larger standard deviation of droplet radius than the superdroplet method, likely due to diffusion associated with the vertical advection of bin variables. However, the spectral width in the bin simulations is insensitive to the grid spacing between 50 and 100 m, suggesting other mechanisms may be compensating for diffusion as the grid spacing is modified. For the precipitating case, larger spectral width in the bin simulations initiates rain earlier and enhances rain development in a positive feedback loop. However, with time, rain formation in the superdroplet simulations catches up to the bin simulations. Offline calculations using the same drop size distributions in both schemes show that the different numerical methods for treating collision–coalescence also contribute to differences in rain formation. The stochastic collision–coalescence in the superdroplet method introduces more variability in drop growth for a given rain mixing ratio.

Supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars—water vapor and temperature—leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

This paper examines the impact of cloud-base turbulence on activation of cloud condensation nuclei (CCN). Following our previous studies, we contrast activation within a nonturbulent adiabatic parcel and an adiabatic parcel filled with turbulence. The latter is simulated by applying a forced implicit large-eddy simulation within a triply periodic computational domain of 643 m3. We consider two monodisperse CCN. Small CCN have a dry radius of 0.01 μm and a corresponding activation (critical) radius and critical supersaturation of 0.6 μm and 1.3%, respectively. Large CCN have a dry radius of 0.2 μm and feature activation radius of 5.4 μm and critical supersaturation 0.15%. CCN are assumed in 200-cm−3 concentration in all cases. Mean cloud-base updraft velocities of 0.33, 1, and 3 m s−1 are considered. In the nonturbulent parcel, all CCN are activated and lead to a monodisperse droplet size distribution above the cloud base, with practically the same droplet size in all simulations. In contrast, turbulence can lead to activation of only a fraction of all CCN with a nonzero spectral width above the cloud base, of the order of 1 μm, especially in the case of small CCN and weak mean cloud-base ascent. We compare our results to studies of the turbulent single-size CCN activation in the Pi chamber. Sensitivity simulations that apply a smaller turbulence intensity, smaller computational domain, and modified initial conditions document the impact of specific modeling assumptions. The simulations call for a more realistic high-resolution modeling of turbulent cloud-base activation.

The impact of giant sea salt aerosols released from breaking waves on rain formation in marine boundary layer clouds is studied using large-eddy simulations (LES). We perform simulations of marine cumuli and stratocumuli for various concentrations of cloud condensation nuclei (CCN) and giant CCN (GCCN). Cloud microphysics are modeled with a Lagrangian method that provides key improvements in comparison to previous LES of GCCN that used Eulerian bin microphysics. We find that GCCN significantly increase precipitation in stratocumuli. This effect is strongest for low and moderate CCN concentrations. GCCN are found to have a smaller impact on precipitation formation in cumuli. These conclusions are in agreement with field measurements. We develop a simple parameterization of the effect of GCCN on precipitation, accretion, and autoconversion rates in marine stratocumuli. Significance Statement Breaking sea waves release salt particles into the atmosphere. Cloud droplets formed on these salt particles can grow larger than droplets formed on other smaller particles. Therefore, sea salt particles can be important for rain formation over oceans. To investigate this effect, we performed idealized computer simulations of stratocumulus and cumulus clouds. Sea salt particles were modeled with an unprecedented precision thanks to the use of an emerging modeling method. In our simulations sea salt particles significantly enhance rain formation in stratocumuli, but not in cumuli. Our study has implications for climate models, because stratocumuli are important for Earth’s energy budget and for rain enhancement experiments.

Motivated by recent discussions concerning differences of convective dynamics in polluted and pristine environments, the so-called convective invigoration in particular, this paper provides an analysis of factors affecting convective updraft buoyancy, such as the in-cloud supersaturation, condensate and precipitation loading, and entrainment. We use the deep convective period from simulations of daytime convection development over land discussed in our previous publications. An entraining parcel framework is used in the theoretical analysis. We show that for the specific case considered here, finite (positive) supersaturation noticeably reduces pseudo-adiabatic parcel buoyancy and cumulative convective available potential energy (cCAPE) in the lower troposphere. This comes from keeping a small fraction of the water vapor in a supersaturated state and thus reducing the latent heating. Such a lower-tropospheric impact is comparable to the effects of condensate loading and entrainment in the idealized parcel framework. For the entire tropospheric depth, loading and entrainment have a much more significant impact on the total CAPE. For the cloud model results, we compare ensemble simulations applying either a bulk microphysics scheme with saturation adjustment or a more comprehensive double-moment scheme with supersaturation prediction. We compare deep convective updraft velocities, buoyancies, and supersaturations from all ensembles. In agreement with the parcel analysis, the saturation-adjustment scheme provides noticeably stronger updrafts in the lower troposphere. For the simulations predicting supersaturation, there are small differences between pristine and polluted conditions below the freezing level that are difficult to explain by standard analysis of the in-cloud buoyancy components. By applying the piggybacking technique, we show that the lower-tropospheric buoyancy differences between pristine and polluted simulations come from a combination of temperature (i.e., latent heating) and condensate loading differences that work together to make polluted buoyancies and updraft velocities slightly larger when compared to their pristine analogues. Overall, the effects are rather small and contradict previous claims of a significant invigoration of deep convection in polluted environments.

This is a rebuttal of Fan and Khain’s comments (hereafter FK21) on a 2020 paper by Grabowski and Morrison (hereafter GM20) that questions the impact of ultrafine cloud condensation nuclei (CCN) on deep convection. GM20 argues that “cold invigoration,” an increase of the updraft speed from lofting and freezing of additional cloud water in polluted environments, is unlikely because the latent heating from freezing of this cloud water approximately recovers the negative impact on the buoyancy from the weight of this water. FK21 suggest a variety of processes that could invalidate our claim. We maintain that our argument is valid and invite the authors to compare their microphysics scheme with ours in the same simplified modeling framework. However, pollution does affect the partitioning of latent heating within the column and likely leads to convection changes beyond a single diurnal cycle through larger-scale circulation changes. This argument explains impacts seen in our idealized mesoscale simulations and in convective–radiative equilibrium simulations by others. We agree with FK21 on the existence of a “warm invigoration” mechanism but question its interpretation. Consistent with the simulations in GM20, we argue that changes in the buoyancy can be explained by the response of the supersaturation to droplet microphysical changes induced by pollution. The buoyancy change is determined by supersaturation differences between pristine and polluted conditions, while condensation rate responds to these supersaturation changes. Finally, we agree with FK21 that the piggybacking modeling technique cannot prove or disprove invigoration; rather, it is a diagnostic technique that can be used to understand mechanisms driving simulation differences.

A single nonprecipitating cumulus congestus setup is applied to compare droplet spectra grown by the diffusion of water vapor in Eulerian bin and particle-based Lagrangian microphysics schemes. Bin microphysics represent droplet spectral evolution applying the spectral density function. In the Lagrangian microphysics, computational particles referred to as superdroplets are followed in time and space with each superdroplet representing a multiplicity of natural cloud droplets. The same cloud condensation nuclei (CCN) activation and identical representation of the droplet diffusional growth allow the comparison. The piggybacking method is used with the two schemes operating in a single simulation, one scheme driving the dynamics and the other one piggybacking the simulated flow. Piggybacking allows point-by-point comparison of droplet spectra predicted by the two schemes. The results show the impact of inherent limitations of the two microphysics simulation methods, numerical diffusion in the Eulerian scheme and a limited number of superdroplets in the Lagrangian scheme. Numerical diffusion in the Eulerian scheme results in a more dilution of the cloud upper half and thus smaller cloud droplet mean radius. The Lagrangian scheme typically has larger spatial fluctuations of droplet spectral properties. A significantly larger mean spectral width in the bin microphysics across the entire cloud depth is the largest difference between the two schemes. A fourfold increase of the number of superdroplets per grid volume and a twofold increase of the spectral resolution and thus the number of bins have small impact on the results and provide only minor changes to the comparison between simulated cloud properties.

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10−6) s−1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Numerical simulations of the impact of ultrafine cloud condensation nuclei (CCN) on deep convection are analyzed to investigate the idea proposed by Fan et al. that addition of ultrafine CCN to an otherwise pristine environment leads to convective invigoration. The piggybacking methodology is applied, allowing rigorous separation of the impact of aerosols from different flow realizations that typically occur when even a small element of the model physics or modeling setup is changed. The setup follows the case of daytime convective development over land based on observations during the Large-Scale Biosphere–Atmosphere (LBA) experiment in Amazonia. Overall, the simulated impacts of ultrafine CCN are similar to the previous study by the authors on the impact of pollution on deep convection. There is no convective invigoration above the freezing level, but there is a small invigoration (increase in vertical velocities) below due to the supersaturation and buoyancy differences in conditions with additional ultrafine CCN compared to unperturbed pristine conditions. As in the previous study, the most significant impact is on the upper-tropospheric convective anvils that feature higher cloud fractions in conditions with ultrafine CCN. The increase comes from purely microphysical considerations as the increased cloud droplet concentrations from ultrafine CCN lead to increased ice crystal concentrations and, consequently, smaller sizes and fall velocities, and longer residence times. Mesoscale organization due to low-level shear has a small effect on the simulated ultrafine CCN impacts. Finally, an alternative explanation of increased lightning above oceanic shipping lines seen in satellite observations and argued to result from convective invigoration is provided.

Scaled-up DNS and implicit LES simulations are used to study turbulent cloud base CCN activation and early growth of cloud droplets. The simulation framework includes a triply periodic computational domain ~1,000 cubic meters filled with inertial-range homogeneous isotropic turbulence. The domain experiences decrease of the mean air temperature and reduction of the mean pressure, both mimicking the rise of an adiabatic air parcel through the cloud base. Results of turbulent simulations are compared to CCN activation and droplet growth within a classical nonturbulent rising parcel. The key difference is a blurriness of the separation between activated and nonactivated (haze) CCN, especially for weak mean ascent rates, when CCN activate and subsequently some deactivate instead of becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the adiabatic nonturbulent parcel once CCN activation is completed. Further increase of the spectral width in the turbulent parcel is similar to that for the initially-monodisperse droplets in the inertial-range homogeneous isotropic turbulence that we and others studied previously, with the standard deviation of the radius squared increasing approximately as the square root of time. This contrasts with the classical nonturbulent parcel framework for which the radius squared standard deviation above the cloud base remains constant because of the parabolic growth of cloud droplets once surface tension and dilute effects can be neglected.

In this study, the super droplet-method (SDM) is used in large-eddy simulations of an isolated cumulus congestus observed during the 1995 Small Cumulus Microphysics Study field project in order to ...

A hybrid bin microphysical scheme is developed in a parcel model framework to study how natural aerosol particles and different types of hygroscopic seeding materials affect the precipitation formation. A novel parameter is introduced to describe the impact of different seeding particles on the evolution of the drop size distribution. The results of more than 100 numerical experiments using the hybrid bin parcel model show that: (a) The Ostwald-ripening effect has a substantial contribution to the broadening of the drop size distribution near the cloud base. The efficiency of this effect increases as the updraft velocity decreases. (b) The efficiency of hygroscopic seeding is significant only if the size of the seeding particles is in the coarse particle size range. The presence of the water-soluble background coarse particles reduces the efficiency of the seeding. (c) The efficient broadening of the size distribution due to the seeding depends on the width of the size distribution of water drops in the control cases, but the relation is not as straightforward as in the case of the glaciogenic seeding.

Climate change affects the dynamics and thermodynamics of moist convection. Changes in the dynamics concern, for instance, an increase of convection strength due to increases of convective available potential energy (CAPE). Thermodynamics involve increases in water vapor that the warmer atmosphere can hold and convection can work with. Small-scale simulations are conducted to separate these two components for daytime development of unorganized convection over land. The simulations apply a novel modeling technique referred to as the piggybacking (or master–slave) approach and consider the global climate model (GCM)-predicted change of atmospheric temperature and moisture profiles in the Amazon region at the end of the century under a business-as-usual scenario. The simulations show that the dynamic impact dominates because changes in cloudiness and rainfall come from cloud dynamics considerations, such as the change in CAPE and convective inhibition (CIN) combined with the impact of environmental relative humidity (RH) on deep convection. The small RH reduction between the current and future climate significantly affects the mean surface rain accumulation as it changes from a small reduction to a small increase when the RH decrease is eliminated. The thermodynamic impact on cloudiness and precipitation is generally small, with the extreme rainfall intensifying much less than expected from an atmospheric moisture increase. These results are discussed in the context of previous studies concerning climate change–induced modifications of moist convection. Future research directions applying the piggybacking method are discussed.

Representation of cloud microphysics is a key aspect of simulating clouds. From the early days of cloud modeling, numerical models have relied on an Eulerian approach for all cloud and thermodynamic and microphysics variables. Over time the sophistication of microphysics schemes has steadily increased, from simple representations of bulk masses of cloud and rain in each grid cell, to including different ice particle types and bulk hydrometeor concentrations, to complex schemes referred to as bin or spectral schemes that explicitly evolve the hydrometeor size distributions within each model grid cell. As computational resources grow, there is a clear trend toward wider use of bin schemes, including their use as benchmarks to develop and test simplified bulk schemes. We argue that continuing on this path brings fundamental challenges difficult to overcome. The Lagrangian particle-based probabilistic approach is a practical alternative in which the myriad of cloud and precipitation particles present in a natural cloud is represented by a judiciously selected ensemble of point particles called superdroplets or superparticles. The advantages of the Lagrangian particle-based approach when compared to the Eulerian bin methodology are explained, and the prospects of applying the method to more comprehensive cloud simulations—for instance, targeting deep convection or frontal cloud systems—are discussed.

This paper discusses a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as super-droplets. Each super-droplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh-Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

Entrainment-mixing and turbulent fluctuations critically impact cloud droplet size distributions (DSDs) in cumulus clouds. This problem is investigated via a new sophisticated modeling framework using the CM1 LES model and a Lagrangian cloud microphysics scheme – the “super-droplet method” (SDM) – coupled with sub-grid-scale (SGS) schemes for particle transport and supersaturation fluctuations. This modeling framework is used to simulate a cumulus congestus cloud. Average DSDs in different cloud regions show broadening from entrainment and secondary cloud droplet activation (activation above the cloud base). DSD width increases with increasing entrainment-induced dilution as expected from past work, except in the most diluted cloud regions. The new modeling framework with SGS transport and supersaturation fluctuations allows a more sophisticated treatment of secondary activation compared to previous studies. In these simulations, it contributes about 25%of the cloud droplet population and impacts DSDs in two contrastingways: narrowing in extremely diluted regions and broadening in relatively less diluted. SGS supersaturation fluctuations contribute significantly to an increase in DSD width via condensation growth and evaporation. Mixing of super-droplets from SGS velocity fluctuations also broadens DSDs. The relative dispersion (ratio of DSD dispersion and mean radius) negatively correlates with grid-scale vertical velocity in updrafts, but is positively correlated in downdrafts. The latter is from droplet activation driven by positive SGS supersaturation fluctuations in grid-mean subsaturated conditions. Finally, the sensitivity to model grid length is evaluated. The SGS schemes have greater influence as the grid length is increased, and they partially compensate for the reduced model resolution.

Motivated by recent discussions concerning differences of convective dynamics in polluted and pristine environments, the so-called convective invigoration in particular, this paper provides an analysis of factors affecting convective updraft buoyancy, such as the in-cloud supersaturation, condensate and precipitation loading, and entrainment. We use the deep convective period from simulations of daytime convection development over land discussed in our previous publications. An entraining parcel framework in used in the theoretical analysis. We show that for the specific case considered here finite (positive) supersaturation noticeably reduces pseudo-adiabatic parcel buoyancy and cumulative CAPE in the lower troposphere. This comes from keeping a small fraction of the water vapor in a supersaturated state and thus reducing the latent heating. Such a lower-tropospheric impact is comparable to the effects of the condensate loading and entrainment in the idealized parcel framework. For the entire tropospheric depth, loading and entrainment have a much more significant impact on the total CAPE. For instance, an increase in the fractional entrainment rate from 0.05 km−1 to 0.3 km−1 reduces the theoretical level of neutral buoyancy from the upper to the middle troposphere and CAPE by a factor of 4. For the cloud model results, we compare ensemble simulations applying either a bulk microphysics scheme with saturation adjustment or a more comprehensive double-moment scheme with supersaturation prediction. The diagnosed bulk fractional entrainment rate, independent of the microphysics scheme applied in the simulations, is either 0.13 or 0.20 km−1 depending on whether we consider profiles of the upper end of the percentile range or of the mean in-cloud equivalent potential temperature. We compare deep convective updrafts, buoyancies, and supersaturations from all ensembles. In agreement with the parcel analysis, the saturation adjustment scheme provides noticeably stronger updrafts in the lower troposphere. For the simulations predicting supersaturation, there are small differences between pristine and polluted conditions below the freezing level that are difficult to explain by standard analysis of the in-cloud buoyancy components. By applying the piggybacking technique, we show that the lower-tropospheric buoyancy differences between pristine and polluted simulations come from a combination of temperature (i.e., latent heating) and condensate loading differences that work together to make polluted buoyancies and updraft velocities slightly larger when compared to their pristine analogues. Overall, the effects are rather small and contradict previous claims of a significant invigoration of deep convection in polluted environments.