ALBERT

Description: Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137Cs and 133Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τe, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τe lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

Description: The local and regional influence of elevated point sources on summertime aerosol forcing and cloud-aerosol interactions in northeastern North America was investigated using the WRF-Chem community model. The direct effects of aerosols on incoming solar radiation were simulated using existing modules to relate aerosol sizes and chemical composition to aerosol optical properties. Indirect effects were simulated by adding a prognostic treatment of cloud droplet number and adding modules that activate aerosol particles to form cloud droplets, simulate aqueous-phase chemistry, and tie a two-moment treatment of cloud water (cloud water mass and cloud droplet number) to an existing radiation scheme. Fully interactive feedbacks thus were created within the modified model, with aerosols affecting cloud droplet number and cloud radiative properties, and clouds altering aerosol size and composition via aqueous processes, wet scavenging, and gas-phase-related photolytic processes. Comparisons of a baseline simulation with observations show that the model captured the general temporal cycle of aerosol optical depths (AODs) and produced clouds of comparable thickness to observations at approximately the proper times and places. The model overpredicted SO2 mixing ratios and PM2.5 mass, but reproduced the range of observed SO2 to sulfate aerosol ratios, suggesting that atmospheric oxidation processes leading to aerosol sulfate formation are captured in the model. The baseline simulation was compared to a sensitivity simulation in which all emissions at model levels above the surface layer were set to zero, thus removing stack emissions. Instantaneous, site-specific differences for aerosol and cloud related properties between the two simulations could be quite large, as removing above-surface emission sources influenced when and where clouds formed within the modeling domain. When summed spatially over the finest resolution model domain (the extent of which corresponds to the typical size of a single GCM grid cell) and temporally over a three day analysis period, total rainfall in the sensitivity simulation increased by 31% over that in the baseline simulation. Fewer optically thin clouds, arbitrarily defined as a cloud exhibiting an optical depth less than 1, formed in the sensitivity simulation. Domain-averaged AODs dropped from 0.46 in the baseline simulation to 0.38 in the sensitivity simulation. The overall net effect of additional aerosols attributable to primary particulates and aerosol precursors from point source emissions above the surface was a domain-averaged reduction of 5 W m−2 in mean daytime downwelling shortwave radiation.

Description: The local and regional influence of elevated point sources on summertime aerosol forcing and cloud-aerosol interactions in northeastern North America was investigated using the WRF-Chem community model. The direct effects of aerosols on incoming solar radiation were simulated using existing modules to relate aerosol sizes and chemical composition to aerosol optical properties. Indirect effects were simulated by adding a prognostic treatment of cloud droplet number and adding modules that activate aerosol particles to form cloud droplets, simulate aqueous-phase chemistry, and tie a two-moment treatment of cloud water (cloud water mass and cloud droplet number) to precipitation and an existing radiation scheme. Fully interactive feedbacks thus were created within the modified model, with aerosols affecting cloud droplet number and cloud radiative properties, and clouds altering aerosol size and composition via aqueous processes, wet scavenging, and gas-phase-related photolytic processes. Comparisons of a baseline simulation with observations show that the model captured the general temporal cycle of aerosol optical depths (AODs) and produced clouds of comparable thickness to observations at approximately the proper times and places. The model overpredicted SO2 mixing ratios and PM2.5 mass, but reproduced the range of observed SO2 to sulfate aerosol ratios, suggesting that atmospheric oxidation processes leading to aerosol sulfate formation are captured in the model. The baseline simulation was compared to a sensitivity simulation in which all emissions at model levels above the surface layer were set to zero, thus removing stack emissions. Instantaneous, site-specific differences for aerosol and cloud related properties between the two simulations could be quite large, as removing above-surface emission sources influenced when and where clouds formed within the modeling domain. When summed spatially over the finest resolution model domain (the extent of which corresponds to the typical size of a single global climate model grid cell) and temporally over a three day analysis period, total rainfall in the sensitivity simulation increased by 31% over that in the baseline simulation. Fewer optically thin clouds, arbitrarily defined as a cloud exhibiting an optical depth less than 1, formed in the sensitivity simulation. Domain-averaged AODs dropped from 0.46 in the baseline simulation to 0.38 in the sensitivity simulation. The overall net effect of additional aerosols attributable to primary particulates and aerosol precursors from point source emissions above the surface was a domain-averaged reduction of 5 W m−2 in mean daytime downwelling shortwave radiation.

Description: We use global models to explore the microphysical effects of carbonaceous aerosols on liquid clouds. Although absorption of solar radiation by soot warms the atmosphere, soot may cause climate cooling due to its contribution to cloud condensation nuclei (CCN) and therefore cloud brightness. Six global models conducted three soot experiments; four of the models had detailed aerosol microphysical schemes. The average cloud radiative response to biofuel soot (black and organic carbon), including both indirect and semi-direct effects, is −0.11 Wm−2, comparable in size but opposite in sign to the respective direct effect. In a more idealized fossil fuel black carbon experiment, some models calculated a positive cloud response because soot provides a deposition sink for sulfuric and nitric acids and secondary organics, decreasing nucleation and evolution of viable CCN. Biofuel soot particles were also typically assumed to be larger and more hygroscopic than for fossil fuel soot and therefore caused more negative forcing, as also found in previous studies. Diesel soot (black and organic carbon) experiments had relatively smaller cloud impacts with five of the models

Description: A clear understanding of particle formation mechanisms is critical for assessing aerosol indirect radiative forcing and associated climate feedback processes. Recent studies reveal the importance of ion-mediated nucleation (IMN) in generating new particles and cloud condensation nuclei (CCN) in the atmosphere. Here we implement the IMN scheme into the Community Atmosphere Model version 5 (CAM5). Our simulations show that, compared to globally averaged results based on H2SO4-H2O binary homogeneous nucleation (BHN), the presence of ionization (i.e., IMN) halves H2SO4 column burden, but increases the column integrated nucleation rate by around one order of magnitude, total particle number burden by a factor of ~3, CCN burden by ~10% (at 0.2% supersaturation) to 65% (at 1.0% supersaturation), and cloud droplet number burden by ~18%. Compared to BHN, IMN increases cloud liquid water path by 7.5%, decreases precipitation by 1.1%, and increases total cloud cover by 1.9%. This leads to an increase of total shortwave cloud radiative forcing (SWCF) by 3.67 W m−2 (more negative) and longwave cloud forcing by 1.78 W m−2 (more positive), with large spatial variations. The effect of ionization on SWCF derived from this study (3.67 W m−2) is a factor of ~3 higher that of a previous study (1.15 W m−2) based on a different ion nucleation scheme and climate model. Based on the present CAM5 simulation, the 5-yr mean impacts of solar cycle induced changes in ionization rates on CCN and cloud forcing are small (~−0.02 W m−2) but have larger inter-annual (from −0.18 to 0.17 W m−2) and spatial variations.

Description: This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

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: The impact of black carbon (BC) aerosols on the global radiation balance is not well constrained. Here twelve global aerosol models are used to show that at least 20% of the present uncertainty in modeled BC direct radiative forcing (RF) is due to diversity in the simulated vertical profile of BC mass. Results are from phases 1 and 2 of the global aerosol model intercomparison project (AeroCom). Additionally, a significant fraction of the variability is shown to come from high altitudes, as, globally, more than 40% of the total BC RF is exerted above 5 km. BC emission regions and areas with transported BC are found to have differing characteristics. These insights into the importance of the vertical profile of BC lead us to suggest that observational studies are needed to better characterize the global distribution of BC, including in the upper troposphere.

Description: We report on the AeroCom Phase II direct aerosol effect (DAE) experiment where 16 detailed global aerosol models have been used to simulate the changes in the aerosol distribution over the industrial era. All 16 models have estimated the radiative forcing (RF) of the anthropogenic DAE, and have taken into account anthropogenic sulphate, black carbon (BC) and organic aerosols (OA) from fossil fuel, biofuel, and biomass burning emissions. In addition several models have simulated the DAE of anthropogenic nitrate and anthropogenic influenced secondary organic aerosols (SOA). The model simulated all-sky RF of the DAE from total anthropogenic aerosols has a range from −0.58 to −0.02 Wm−2, with a mean of −0.27 Wm−2 for the 16 models. Several models did not include nitrate or SOA and modifying the estimate by accounting for this with information from the other AeroCom models reduces the range and slightly strengthens the mean. Modifying the model estimates for missing aerosol components and for the time period 1750 to 2010 results in a mean RF for the DAE of −0.35 Wm−2. Compared to AeroCom Phase I (Schulz et al., 2006) we find very similar spreads in both total DAE and aerosol component RF. However, the RF of the total DAE is stronger negative and RF from BC from fossil fuel and biofuel emissions are stronger positive in the present study than in the previous AeroCom study. We find a tendency for models having a strong (positive) BC RF to also have strong (negative) sulphate or OA RF. This relationship leads to smaller uncertainty in the total RF of the DAE compared to the RF of the sum of the individual aerosol components. The spread in results for the individual aerosol components is substantial, and can be divided into diversities in burden, mass extinction coefficient (MEC), and normalized RF with respect to AOD. We find that these three factors give similar contributions to the spread in results.