Search

The effects of aerosol perturbations on deep convective clouds have received considerable attention in the recent literature, especially from a modeling perspective. The published responses in precipitation amount and updraft strength vary in both sign and magnitude and may be the result of different models and parameterizations. Here, a simple numerical framework is employed to determine the potential effects of warming both below the freezing level (warm invigoration) and above the freezing level (mixed-phase invigoration) due to increased aerosol loading. The role of updraft width and slope in the same framework is also examined, highlighting the relative importance of each factor on the resulting updraft strength. The results show that the potential effects of warm invigoration are 2–3 times larger than for mixed-phase invigoration. However, a relatively small response in updraft velocity to warming is found, especially in comparison with the predicted changes in updraft velocity due to small differences in system slope and width, with 87.7% and 96.4% of the subadiabatic and adiabatic realizations, respectively, showing changes in updraft velocity of less than 15% for warmings of no more than 2°C. This result suggests that observations of the aerosol effect will be largely muddled by the natural variability of convective updraft width and slope (which are related to environmental wind shear).

Changes in the aerosol number concentration are reflected by changes in raindrop size and number concentration that ultimately affect the strength of cold pools via evaporation. Therefore, aerosol perturbations can potentially alter the balance between cold pool–induced and low-level wind shear–induced circulations. In the present work, simulations with increased aerosol loadings below approximately 3 km, between approximately 3 and 10 km, and at all vertical levels are performed to specifically address both the overall sensitivity of a squall line to the vertical distribution of aerosols and the extent to which low-level aerosols can affect the convective strength of the system. The results suggest that low-level aerosol perturbations have a negligible effect on the overall storm strength even though they act to enhance low-level latent heating rates. A tracer analysis shows that the low-level aerosols are either predominantly detrained at or below the freezing level or are rapidly lifted to the top of the troposphere or the lower stratosphere within the strongest convective cores. Moreover, it is shown that midlevel aerosol perturbations have nearly the same effect as perturbing the entire domain, increasing the convective updraft mass flux by more than 10%. These changes in strength are driven by a complex chain of events caused by smaller supercooled droplets, larger graupel, and larger raindrops. Combined, these changes tend to reduce the low-level bulk evaporation rate, thus weakening the cold pool and enhancing updraft strength. The results presented herein suggest that midlevel aerosol perturbations may exhibit a much larger effect on squall lines, at least in the context of this idealized framework.

The dynamical effects of increased aerosol loading on the strength and structure of numerically simulated squall lines are explored. Results are explained in the context of Rotunno–Klemp–Weisman (RKW) theory. Changes in aerosol loading lead to changes in raindrop size and number that ultimately affect the strength of the cold pool via changes in evaporation. Thus, the balance between cold pool and low-level wind shear–induced vorticities can be changed by an aerosol perturbation. Simulations covering a wide range of low-level wind shears are performed to study the sensitivity to aerosols in different environments and provide more general conclusions. Simulations with relatively weak low-level environmental wind shear (0.0024 s−1) have a relatively strong cold pool circulation compared to the environmental shear. An increase in aerosol loading leads to a weakening of the cold pool and, hence, a more optimal balance between the cold pool– and environmental shear–induced circulations according to RKW theory. Consequently, there is an increase in the convective mass flux of nearly 20% in polluted conditions relative to pristine. This strengthening coincides with more upright convective updrafts and a significant increase (nearly 20%) in cumulative precipitation. An increase in aerosol loading in a strong wind shear environment (0.0064 s−1) leads to less optimal storms and a suppression of the convective mass flux and precipitation. This occurs because the cold pool circulation is weak relative to the environmental shear when the shear is strong, and further weakening of the cold pool with high aerosol loading leads to an even less optimal storm structure (i.e., convective updrafts begin to tilt downshear).

A novel two-moment bulk aerosol parameterization is derived from a state-of-the-art 2D bin microphysics model using power-law relationships and a semi-analytical technique for activation. The activation scheme predicts both number and mass of a lognormal aerosol distribution and permits the evolution of the modal mass with time. The newly developed bulk aerosol scheme is formulated for use in traditional two-moment bulk microphysics models. The new explicit scheme is compared with the 2D bin scheme and a simple scaling aerosol parameterization, in which all the aerosol processes are scaled to the respective cloud process rates, in a kinematic model with a specified flow field. Hybrid simulations in which the explicit activation formulation is coupled to the scaling parameterization are also performed. Model results demonstrate the significance of including a physically realistic representation of aerosols contained in haze, cloud droplets, and rain. It is shown that the explicit aerosol parameterization and scaling method predict similar bulk aerosol quantities and match the results of the 2D bin model only if an explicit treatment of aerosol activation—that is, both aerosol number and mass transfer because of activation—is included in the microphysics model.

Three configurations of a bulk microphysics scheme in conjunction with a detailed bin scheme are implemented in the Weather Research and Forecasting (WRF) model to specifically address the role of the saturation adjustment assumption (i.e., condensing/evaporating the surplus/deficit water vapor relative to saturation in one time step) on aerosol-induced invigoration of deep convective clouds. The bulk model configurations are designed to treat cloud droplet condensation/evaporation using either saturation adjustment, as employed in most bulk models, or an explicit representation of supersaturation over a time step, as used in bin models. Results demonstrate that the use of saturation adjustment artificially enhances condensation and latent heating at low levels and limits the potential for an increase in aerosol concentration to increase buoyancy at mid to upper levels. This leads to a small weakening of the time- and domain-averaged convective mass flux (~ -3%) in polluted compared to clean conditions. In contrast, the bin model and bulk scheme with explicit prediction of supersaturation simulate an increase in latent heating aloft and the convective updraft mass flux is weakly invigorated (~5%). The bin model also produces a large increase in domain-mean cumulative surface precipitation in polluted conditions (~18%), while all of the bulk model configurations simulate little change in precipitation. Finally, it is shown that the cold pool weakens substantially with increased aerosol loading when saturation adjustment is applied, which acts to reduce the low-level convergence and weaken the convective dynamics. With an explicit treatment of supersaturation in the bulk and bin models there is little change in cold pool strength, so that the convective response to polluted conditions is influenced more by changes in latent heating aloft. It is concluded that the use of saturation adjustment can explain differences in the response of cold pool evolution and convective dynamics with aerosol loading simulated by the bulk and bin models, but cannot explain large differences in the response of surface precipitation between these models.

A two-dimensional (2-D) continuous spectral aerosol-droplet microphysics model is presented and implemented into the Weather Research and Forecasting (WRF) model for large-eddy simulations (LES) of warm clouds. Activation and regeneration of aerosols are treated explicitly in the calculation of condensation/evaporation. The model includes a 2-D spectrum that encompasses wet aerosol particles (i.e., haze droplets), cloud droplets, and drizzle droplets in a continuous and consistent manner and allows for the explicit tracking of aerosol size within cloud droplets due to collision-coalescence. The system of differential equations describing condensation/evaporation (i.e., mass conservation and energy conservation) is solved simultaneously within each grid cell. The model is demonstrated by simulating a marine stratocumulus deck for two different aerosol loadings (100 and 500 cm−3), and comparison with the more traditional microphysics modeling approaches (both 1-D bin and bulk schemes) is evaluated. The simulations suggest that in a 1-D bin microphysics scheme, without regeneration, too few particles are produced and hence the mode of the droplet size spectrum occurs at a larger size relative to the 2-D bin model results. Moreover, with regeneration, the 1-D scheme produces too many small droplets and thus shifts the mode toward smaller sizes. These large shifts in the droplet size distribution can potentially have significant effects on the efficiency of the collision-coalescence process, fall speeds, and ultimately precipitation.

Three-dimensional large-eddy simulations (LES) with detailed bin-resolved microphysics are performed to explore the diurnal variation of marine stratocumulus (MSc) clouds under clean and polluted conditions. The sensitivity of the aerosol-cloud-precipitation interactions to variation of sea surface temperature, free tropospheric humidity, large-scale divergence rate, and wind speed is assessed. The comprehensive set of simulations corroborates previous studies that (1) with moderate/heavy drizzle, an increase in aerosol leads to an increase in cloud thickness; and (2) with non/light drizzle, an increase in aerosol results in a thinner cloud, due to the pronounced effect on entrainment. It is shown that for higher SST, stronger large-scale divergence, drier free troposphere, or lower wind speed, the cloud thins and precipitation decreases. The sign and magnitude of the Twomey effect, droplet dispersion effect, cloud thickness effect, and cloud optical depth susceptibility to aerosol perturbations (i.e., change in cloud optical depth to change in aerosol number concentration) are evaluated by LES experiments and compared with analytical formulations. The Twomey effect emerges as dominant in total cloud optical depth susceptibility to aerosol perturbations. The dispersion effect, that of aerosol perturbations on the cloud droplet size spectrum, is positive (i.e., increase in aerosol leads to spectral narrowing) and accounts for 3% to 10% of the total cloud optical depth susceptibility at nighttime, with greater influence in heavier drizzling clouds. The cloud thickness effect is negative (i.e., increase in aerosol leads to thinner cloud) for non/light drizzling cloud and positive for a moderate/heavy drizzling clouds; the cloud thickness effect contributes 5% to 22% of the nighttime total cloud susceptibility. Overall, the total cloud optical depth susceptibility ranges from ~0.28 to 0.53 at night; an increase in aerosol concentration enhances cloud optical depth, especially with heavier precipitation and in a more pristine environment. During the daytime, the range of magnitude for each effect is more variable owing to cloud thinning and decoupling. The good agreement between LES experiments and analytical formulations suggests that the latter may be useful in evaluations of the total cloud susceptibility. The ratio of the magnitude of the cloud thickness effect to that of the Twomey effect depends on cloud base height and cloud thickness in unperturbed (clean) clouds.

The potential effects of increased aerosol loading on the development of deep convective clouds and resulting precipitation amounts are studied by employing the Weather Research and Forecasting (WRF) model as a detailed high-resolution cloud resolving model (CRM) with both detailed bulk and bin microphysics schemes. The bulk microphysics scheme incorporates a physically based parameterization of cloud droplet activation as well as homogeneous and heterogeneous freezing in order to explicitly resolve the possible aerosol-induced effects on the cloud microphysics. These parameterizations allow one to segregate the effects of an increase in the aerosol number concentration into enhanced cloud condensation nuclei (CCN) and/or ice nuclei (IN) concentrations using bulk microphysics. The bin microphysics scheme, with its explicit calculations of cloud particle collisions, is shown to better predict cumulative precipitation. Increases in the CCN number concentration may not have a monotonic influence on the cumulative precipitation resulting from deep convective clouds. We demonstrate that the aerosol-induced effect is controlled by the balance between latent heating and the increase in condensed water aloft, each having opposing effects on buoyancy. It is also shown that under polluted conditions and in relatively dry environments, increases in the CCN number concentration reduce the cumulative precipitation due to the competition between the sedimentation and evaporation/sublimation timescales. The effect of an increase in the IN number concentration on the dynamics of deep convective clouds is small, but may act to suppress precipitation. A comparison of the predictions using the bin and bulk microphysics schemes demonstrate a significant difference between the predicted precipitation and the influence of aerosol perturbations on updraft velocity within the convective core. The bulk microphysics scheme is shown to be unable to capture the changes in latent heating that occur as a result of changes in the CCN number concentration, while the bin microphysics scheme demonstrates significant increases in the latent heating aloft with increasing CCN number concentration. This suggests that a detailed two-bulk microphysics scheme, which is more computationally efficient than bin microphysics schemes, may not be sufficient, even when coupled to a detailed activation scheme, to predict small changes that result from perturbations in aerosol loading.