ALBERT

Description: Much of the large uncertainty in estimates of anthropogenic aerosol effects on climate arises from the multi-scale nature of the interactions between aerosols, clouds and dynamics, which are difficult to represent in conventional general circulation models (GCMs). In this study, we use a multi-scale aerosol-climate model that treats aerosols and clouds across multiple scales to study aerosol indirect effects. This multi-scale aerosol-climate model is an extension of a multi-scale modeling framework (MMF) model that embeds a cloud-resolving model (CRM) within each vertical column of a GCM grid. The extension allows a more physically-based treatment of aerosol-cloud interactions in both stratiform and convective clouds on the global scale in a computationally feasible way. Simulated model fields, including liquid water path (LWP), ice water path, cloud fraction, shortwave and longwave cloud forcing, precipitation, water vapor, and cloud droplet number concentration are in reasonable agreement with observations. The new model performs quantitatively similar to the previous version of the MMF model in terms of simulated cloud fraction and precipitation. The simulated change in shortwave cloud forcing from anthropogenic aerosols is −0.77 W m−2, which is less than half of that (−1.79 W m−2) calculated by the host GCM (NCAR CAM5) with traditional cloud parameterizations and is also at the low end of the estimates of other conventional global aerosol-climate models. The smaller forcing in the MMF model is attributed to a smaller (3.9 %) increase in LWP from preindustrial conditions (PI) to present day (PD) compared with 15.6 % increase in LWP in stratiform clouds in CAM5. The difference is caused by a much smaller response in LWP to a given perturbation in cloud condensation nuclei (CCN) concentrations from PI to PD in the MMF (about one-third of that in CAM5), and, to a lesser extent, by a smaller relative increase in CCN concentrations from PI to PD in the MMF (about 26 % smaller than that in CAM5). The smaller relative increase in CCN concentrations in the MMF is caused in part by a smaller increase in aerosol lifetime from PI to PD in the MMF, a positive feedback in aerosol indirect effects induced by cloud lifetime effects from aerosols. The smaller response in LWP to anthropogenic aerosols in the MMF model is consistent with observations and with high resolution model studies, which may indicate that aerosol indirect effects simulated in conventional global climate models are overestimated and point to the need to use global high resolution models, such as MMF models or global CRMs, to study aerosol indirect effects. The simulated total anthropogenic aerosol effect in the MMF is −1.05 W m−2, which is close to the Murphy et al. (2009) inverse estimate of −1.1±0.4 W m−2 (1σ) based on the examination of the Earth's energy balance. Further improvements in the representation of ice nucleation and low clouds in MMF are needed to refine the aerosol indirect effect estimate.

Description: Detailed measurements of ice crystals in cirrus clouds were used to compare with results from the Community Atmospheric Model Version 5 (CAM5) global climate model. The observations are from two different field campaigns with contrasting conditions: Atmospheric Radiation Measurements Spring Cloud Intensive Operational Period in 2000 (ARM-IOP), which was characterized primarily by midlatitude frontal clouds and cirrus, and Tropical Composition, Cloud and Climate Coupling (TC4), which was dominated by anvil cirrus. Results show that the model typically overestimates the slope parameter of the exponential size distributions of cloud ice and snow, while the variation with temperature (height) is comparable. The model also overestimates the ice/snow number concentration (0th moment of the size distribution) and underestimates higher moments (2nd through 5th), but compares well with observations for the 1st moment. Overall the model shows better agreement with observations for TC4 than for ARM-IOP in regards to the moments. The mass-weighted terminal fall speed is lower in the model compared to observations for both ARM-IOP and TC4, which is partly due to the overestimation of the size distribution slope parameter. Sensitivity tests with modification of the threshold size for cloud ice to snow autoconversion (Dcs) do not show noticeable improvement in modeled moments, slope parameter and mass weighed fall speed compared to observations. Further, there is considerable sensitivity of the cloud radiative forcing to Dcs, consistent with previous studies, but no value of Dcs improves modeled cloud radiative forcing compared to measurements. Since the autoconversion of cloud ice to snow using the threshold size Dcs has little physical basis, future improvement to combine cloud ice and snow into a single category, eliminating the need for autoconversion, is suggested.

Description: A new stratiform cloud scheme including a two-moment bulk microphysics module, a cloud cover parameterization allowing ice supersaturation, and an ice nucleation parameterization has been implemented into the recently developed GFDL AM3 general circulation model (GCM) as part of an effort to treat aerosol-cloud-radiation interactions more realistically. Unlike the original scheme, the new scheme facilitates the study of cloud-ice-aerosol interactions via influences of dust and sulfate on ice nucleation. While liquid and cloud ice water path associated with stratiform clouds are similar for the new and the original scheme, column integrated droplet numbers and global frequency distributions (PDFs) of droplet effective radii differ significantly. This difference is in part due to a difference in the implementation of the Wegener-Bergeron-Findeisen (WBF) mechanism, which leads to a larger contribution from super-cooled droplets in the original scheme. Clouds are more likely to be either completely glaciated or liquid due to the WBF mechanism in the new scheme. Super-saturations over ice simulated with the new scheme are in qualitative agreement with observations, and PDFs of ice numbers and effective radii appear reasonable in the light of observations. Especially, the temperature dependence of ice numbers qualitatively agrees with in-situ observations. The global average long-wave cloud forcing decreases in comparison to the original scheme as expected when super-saturation over ice is allowed. Anthropogenic aerosols lead to a larger decrease in short-wave absorption (SWABS) in the new model setup, but outgoing long-wave radiation (OLR) decreases as well, so that the net effect of including anthropogenic aerosols on the net radiation at the top of the atmosphere (netradTOA = SWABS-OLR) is of similar magnitude for the new and the original scheme.

Description: Cloud microphysical process rates control the amount of condensed water in clouds and impact the susceptibility of precipitation to cloud-drop number and aerosols. The relative importance of different microphysical processes in a climate model is analyzed, and the autoconversion and accretion processes are found to be critical to the condensate budget in most regions. A simple steady-state model of warm rain formation is used to illustrate that the diagnostic rain formulations typical of climate models may result in excessive contributions from autoconversion, compared to observations and large eddy simulation models with explicit bin-resolved microphysics and rain formation processes. The behavior does not appear to be caused by the bulk process rate formulations themselves, because the steady-state model with the same bulk accretion and autoconversion has reduced contributions from autoconversion. Sensitivity tests are conducted to analyze how perturbations to the precipitation microphysics for stratiform clouds impact process rates, precipitation susceptibility and aerosol–cloud interactions (ACI). With similar liquid water path, corrections for the diagnostic rain assumptions in the GCM based on the steady-state model to boost accretion indicate that the radiative effects of ACI may decrease by 20% in the GCM. Links between process rates, susceptibility and ACI are not always clear in the GCM. Better representation of the precipitation process, for example by prognosticating precipitation mass and number, may help better constrain these effects in global models with bulk microphysics schemes.

Description: Aerosol indirect effects continue to constitute one of the most important uncertainties for anthropogenic climate perturbations. Within the international AEROCOM initiative, the representation of aerosol-cloud-radiation interactions in ten different general circulation models (GCMs) is evaluated using three satellite datasets. The focus is on stratiform liquid water clouds since most GCMs do not include ice nucleation effects, and none of the model explicitly parameterises aerosol effects on convective clouds. We compute statistical relationships between aerosol optical depth (τa) and various cloud and radiation quantities in a manner that is consistent between the models and the satellite data. It is found that the model-simulated influence of aerosols on cloud droplet number concentration (Nd) compares relatively well to the satellite data at least over the ocean. The relationship between τa and liquid water path is simulated much too strongly by the models. This suggests that the implementation of the second aerosol indirect effect mainly in terms of an autoconversion parameterisation has to be revisited in the GCMs. A positive relationship between total cloud fraction (fcld) and τa as found in the satellite data is simulated by the majority of the models, albeit less strongly than that in the satellite data in most of them. In a discussion of the hypotheses proposed in the literature to explain the satellite-derived strong fcld–τa relationship, our results indicate that none can be identified as a unique explanation. Relationships similar to the ones found in satellite data between τa and cloud top temperature or outgoing long-wave radiation (OLR) are simulated by only a few GCMs. The GCMs that simulate a negative OLR–τa relationship show a strong positive correlation between τa and fcld. The short-wave total aerosol radiative forcing as simulated by the GCMs is strongly influenced by the simulated anthropogenic fraction of τa, and parameterisation assumptions such as a lower bound on Nd. Nevertheless, the strengths of the statistical relationships are good predictors for the aerosol forcings in the models. An estimate of the total short-wave aerosol forcing inferred from the combination of these predictors for the modelled forcings with the satellite-derived statistical relationships yields a global annual mean value of −1.5±0.5 Wm−2. In an alternative approach, the radiative flux perturbation due to anthropogenic aerosols can be broken down into a component over the cloud-free portion of the globe (approximately the aerosol direct effect) and a component over the cloudy portion of the globe (approximately the aerosol indirect effect). An estimate obtained by scaling these simulated clear- and cloudy-sky forcings with estimates of anthropogenic τa and satellite-retrieved Nd–τa regression slopes, respectively, yields a global, annual-mean aerosol direct effect estimate of −0.4±0.2 Wm−2 and a cloudy-sky (aerosol indirect effect) estimate of −0.7±0.5 Wm−2, with a total estimate of −1.2±0.4 Wm−2.

Description: 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.

Description: Cloud-system resolving simulations with the chemistry version of the Weather Research and Forecasting (WRF-Chem) model are used to quantify the relative impacts of regional anthropogenic and oceanic emissions on changes in aerosol properties, cloud macro- and microphysics, and cloud radiative forcing over the Southeast Pacific (SEP) during the VAMOS Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx) (15 October–16 November 2008). Two distinct regions are identified. The near-coast polluted region is characterized by low surface precipitation rates, the strong suppression of non-sea-salt particle activation due to sea-salt particles, a predominant albedo effect in aerosol indirect effects, and limited impact of aerosols associated with anthropogenic emissions on clouds. Opposite sensitivities to natural marine and anthropogenic aerosol perturbations are seen in cloud properties (e.g., cloud optical depth and cloud-top and cloud-base heights), precipitation, and the top-of-atmosphere and surface shortwave fluxes over this region. The relatively clean remote region is characterized by large contributions of aerosols from non-regional sources (lateral boundaries) and much stronger drizzle at the surface. Under a scenario of five-fold increase in regional anthropogenic emissions, this relatively clean region shows large cloud responses, for example, a 13% increase in cloud-top height and a 9% increase in albedo in response to a moderate increase (25% of the reference case) in cloud condensation nuclei (CCN) concentration. The reduction of precipitation due to this increase in anthropogenic aerosols more than doubles the aerosol lifetime in the clean marine boundary layer. Therefore, the aerosol impacts on precipitation are amplified by the positive feedback of precipitation on aerosol, which ultimately alters the cloud micro- and macro-physical properties, leading to strong aerosol-cloud-precipitation interactions. The high sensitivity is also related to an increase in cloud-top entrainment rate (by 16% at night) due to the increased anthropogenic aerosols. The simulated aerosol-cloud-precipitation interactions due to the increased anthropogenic aerosols have a stronger diurnal cycle over the clean region compared to the near-coast region with stronger interactions at night. During the day, solar heating results in more frequent decoupling of the cloud and sub-cloud layers, thinner clouds, reduced precipitation, and reduced sensitivity to the increase in anthropogenic emissions. This study shows the importance of natural aerosols in accurately quantifying anthropogenic forcing within a regional modeling framework. The results of this study also imply that the energy balance perturbations from increased anthropogenic emissions are larger in the more susceptible clean environment than in already polluted environment and are larger than possible from the first indirect effect alone.

Description: A cloud system-resolving model (the Weather Research and Forecasting model) with 1 km horizontal grid spacing is used to investigate the response of an idealized supercell storm to increased cloud droplet concentrations associated with polluted conditions. The primary focus is on exploring robustness of simulated aerosol effects in the face of complex process interactions and feedbacks between the cloud microphysics and dynamics. Simulations are run using sixteen different model configurations with various microphysical or thermodynamic processes modified or turned off. Robustness of the storm response to polluted conditions is also explored for each configuration by performing additional simulations with small perturbations to the initial conditions. Differences in the domain-mean accumulated surface precipitation and convective mass flux between polluted and pristine conditions are small for almost all model configurations, with relative differences in each quantity generally less than 15%. Configurations that produce a decrease (increase) in cold pool strength in polluted conditions also tend to simulate a decrease (increase) in surface precipitation and convective mass flux. Combined with an analysis of the dynamical and thermodynamic fields, these results indicate the importance of interactions between microphysics, cold pool evolution, and dynamics along outflow boundaries in explaining the system response. Several model configurations, including the baseline, produce an overall similar storm response (weakening) in polluted conditions despite having different microphysical or thermodynamic processes turned off. With hail initiation turned off or the hail fallspeed-size relation set to that of snow, the model produces an invigoration instead of weakening of the storm in polluted conditions. These results highlight the difficulty of foreseeing impacts of changes to model parameterizations and isolating process interactions that drive the system response to aerosols. Overall, these findings are robust, in a qualitative sense, to small perturbations in the initial conditions. However, there is sensitivity in the magnitude, and in some cases sign, of the storm response to polluted conditions with small perturbations in the temperature of the thermal used to initiate convection (less than ±0.5 K) or the vertical shear of the environmental wind (±5%). It is concluded that reducing uncertainty in simulations of aerosol effects on individual deep convective storms will likely require ensemble methods in addition to continued improvement of model parameterizations.