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: Among the numerous atmospheric constituents, aerosols play a unique role on climate, due to their scattering and absorbing capabilities, visibility degradation and their effect on incoming and outgoing radiation. The most important optical properties are the aerosol optical depth (AOD), the asymmetry parameter (g) and the single scattering albedo (SSA). Uncertainties in aerosol microphysics in global models, which in turn affect their optical properties, propagate to uncertainties on the effect of aerosols on climate. This study aims to estimate the uncertainty of AOD, g and SSA attributable to the aerosol representation in models, namely mixing state, aerosol size and aerosol associated water. As a reference, the monthly mean output of the general circulation model LMDz-INCA from the international comparison exercise AEROCOM B was used. For the optical properties calculations, aerosols were considered either externally mixed, homogeneously internally mixed or coated spheres. The radius was allowed to vary by ±20% (with 2% intervals) and the aerosol water content by ±50% (with 5% intervals) with respect to the reference model output. All of these possible combinations were assumed to be equally likely and the optical properties were calculated for each one of them. A probability density function (PDF) was constructed at each model grid point for AOD, g and SSA. From this PDF, the 1σ and 2σ uncertainties of the AOD, g and SSA were calculated and are available as global maps for each month. For the range of the cases studied, we derive a maximum 2σ uncertainty range in AOD of 70%, while for g and SSA the maxima reach 18% and 28% respectively. The mixing state was calculated to be important, with the aerosol absorption and SSA being the most affected properties when absorbing aerosols are present.

Description: We examine the response of Arctic gas and aerosol concentrations to perturbations in pollutant emissions from Europe, East and South Asia, and North America using results from a coordinated model intercomparison. These sensitivities to regional emissions (mixing ratio change per unit emission) vary widely across models and species. Intermodel differences are systematic, however, so that the relative importance of different regions is robust. North America contributes the most to Arctic ozone pollution. For aerosols and CO, European emissions dominate at the Arctic surface but Asian emissions become progressively more important with altitude, and are dominant in the upper troposphere. Sensitivities show strong seasonality: surface sensitivities typically maximize during boreal winter for European and during spring for East Asian and North American emissions. Mid-tropospheric sensitivities, however, nearly always maximize during spring or summer for all regions. Deposition of black carbon (BC) onto Greenland is most sensitive to North American emissions. North America and Europe each contribute ~40% of total BC deposition to Greenland, with ~20% from East Asia. Elsewhere in the Arctic, both sensitivity and total BC deposition are dominated by European emissions. Model diversity for aerosols is especially large, resulting primarily from differences in aerosol physics and removal. Comparison of aerosols with observations indicates problems in either the models or interpretation of the measurements. For gas phase pollutants such as CO and O3, which are relatively well-simulated, the processes contributing most to uncertainties depend on the source region. Uncertainties in the Arctic surface CO response to emissions perturbations are dominated by emissions for East Asian sources, while uncertainties in transport, emissions, and oxidation are comparable for European and North American sources. At higher levels, model-to-model variations in transport and oxidation are most important. Differences in photochemistry appear to play the largest role in the intermodel variations in Arctic ozone sensitivity.

Description: This is the first article of a series presenting a detailed analysis of bromine chemistry simulated with the atmospheric chemistry general circulation model ECHAM5/MESSy. Release from sea salt is an important bromine source, hence the model explicitly calculates aerosol chemistry and phase partitioning for coarse mode aerosol particles. Many processes including chemical reaction rates are influenced by the particle size distribution, and aerosol associated water strongly affects the aerosol pH. Knowledge of the aerosol pH is important as it determines the aerosol chemistry, e.g., the efficiency of sulphur oxidation and bromine release. Here, we focus on the simulated sea salt aerosol size distribution and the coarse mode aerosol pH. A comparison with available field data shows that the simulated aerosol distributions agree reasonably well within the range of measurements. In spite of the small number of aerosol pH measurements and the uncertainty in its experimental determination, the simulated aerosol pH compares well with the observations. The aerosol pH ranges from alkaline aerosol in areas of strong production down to pH values of 1 over regions of medium sea salt production and high levels of gas phase acids, mostly polluted regions over the oceans in the northern hemisphere.

Description: The effects of unified aerosol sources on global aerosol fields simulated by different models are examined in this paper. We compare results from two AeroCom experiments, one with different (ExpA) and one with unified emissions, injection heights, and particle sizes at the source (ExpB). Surprisingly, harmonization of aerosol sources has only a small impact on the simulated diversity for aerosol burden, and consequently optical properties, as the results are largely controlled by model-specific transport, removal, chemistry (leading to the formation of secondary aerosols) and parameterizations of aerosol microphysics (e.g. the split between deposition pathways) and to a lesser extent on the spatial and temporal distributions of the (precursor) emissions. The burdens of black carbon and especially sea salt become more coherent in ExpB only, because the large ExpA diversity for these two species was caused by few outliers. The experiment also indicated that despite prescribing emission fluxes and size distributions, ambiguities in the implementation in individual models can lead to substantial differences. These results indicate the need for a better understanding of aerosol life cycles at process level (including spatial dispersal and interaction with meteorological parameters) in order to obtain more reliable results from global aerosol simulations. This is particularly important as such model results are used to assess the consequences of specific air pollution abatement strategies.

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: The CALIOPE air quality modelling system is developed and applied to Europe with high spatial resolution (12 km × 12 km). The modelled daily-to-seasonal aerosol variability over Europe in 2004 is evaluated and analysed. Aerosols are estimated from two models, CMAQv4.5 (AERO4) and BSC-DREAM8b. CMAQv4.5 calculates biogenic, anthropogenic and sea salt aerosol and BSC-DREAM8b provides the natural mineral dust contribution from North African deserts. For the evaluation, we use daily PM10, PM2.5 and aerosol components data from 55 stations of the EMEP/CREATE network and total, coarse and fine aerosol optical depth (AOD) data from 35 stations of the AERONET sun photometer network. Annual correlations between modelled and observed values for PM10 and PM2.5 are 0.55 and 0.47, respectively. Correlations for total, coarse and fine AOD are 0.51, 0.63, and 0.53, respectively. The higher correlations of the PM10 and the coarse mode AOD are largely due to the accurate representation of the African dust influence in the forecasting system. Overall PM and AOD levels are underestimated. The evaluation of the aerosol components highlights underestimations in the fine fraction of carbonaceous matter (EC and OC) and secondary inorganic aerosols (SIA; i.e. nitrate, sulphate and ammonium). The scores of the bulk parameters are significantly improved after applying a simple model bias correction based on the observed aerosol composition. The simulated PM10 and AOD present maximum values over the industrialized and populated Po Valley and Benelux regions. SIA are dominant in the fine fraction representing up to 80% of the aerosol budget in latitudes north of 40° N. In southern Europe, high PM10 and AOD are linked to the desert dust transport from the Sahara which contributes up to 40% of the aerosol budget. Maximum seasonal ground-level concentrations (PM10 〉 30 μg m−3) are found between spring and early autumn. We estimate that desert dust causes daily exceedances of the PM10 European air quality limit value (50 μg m−3) in large areas south of 45° N with more than 75 exceedances per year in the southernmost regions.

Description: The AeroCom exercise diagnoses multi-component aerosol modules in global modeling. In an initial assessment simulated global distributions for mass and mid-visible aerosol optical thickness (aot) were compared among 20 different modules. Model diversity was also explored in the context of previous comparisons. For the component combined aot general agreement has improved for the annual global mean. At 0.11 to 0.14, simulated aot values are at the lower end of global averages suggested by remote sensing from ground (AERONET ca. 0.135) and space (satellite composite ca. 0.15). More detailed comparisons, however, reveal that larger differences in regional distribution and significant differences in compositional mixture remain. Of particular concern are large model diversities for contributions by dust and carbonaceous aerosol, because they lead to significant uncertainty in aerosol absorption (aab). Since aot and aab, both, influence the aerosol impact on the radiative energy-balance, the aerosol (direct) forcing uncertainty in modeling is larger than differences in aot might suggest. New diagnostic approaches are proposed to trace model differences in terms of aerosol processing and transport: These include the prescription of common input (e.g. amount, size and injection of aerosol component emissions) and the use of observational capabilities from ground (e.g. measurements networks) or space (e.g. correlations between aerosol and clouds).

Description: Simulation results of global aerosol models have been assembled in the framework of the AeroCom intercomparison exercise. In this paper, we analyze the life cycles of dust, sea salt, sulfate, black carbon and particulate organic matter as simulated by sixteen global aerosol models. The differences among the results (model diversities) for sources and sinks, burdens, particle sizes, water uptakes, and spatial dispersals have been established. These diversities have large consequences for the calculated radiative forcing and the aerosol concentrations at the surface. Processes and parameters are identified which deserve further research. The AeroCom all-models-average emissions are dominated by the mass of sea salt (SS), followed by dust (DU), sulfate (SO4), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead generally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO4, is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO4-sources is mainly caused by the disagreement on depositional loss of precursor gases and on chemical production. The diversities of the emissions are passed on to the burdens, but the latter are also strongly affected by the model-specific treatments of transport and aerosol processes. The burdens of dry masses decrease from largest to smallest: DU, SS, SO4, POM, and BC. The all-models-average residence time is shortest for SS with about half a day, followed by SO4 and DU with four days, and POM and BC with six and seven days, respectively. The wet deposition rate is controlled by the solubility and increases from DU, BC, POM to SO4 and SS. It is the dominant sink for SO4, BC, and POM, and contributes about one third to the total removal of SS and DU species. For SS and DU we find high diversities for the removal rate coefficients and deposition pathways. Models do neither agree on the split between wet and dry deposition, nor on that between sedimentation and other dry deposition processes. We diagnose an extremely high diversity for the uptake of ambient water vapor that influences the particle size and thus the sink rate coefficients. Furthermore, we find little agreement among the model results for the partitioning of wet removal into scavenging by convective and stratiform rain. Large differences exist for aerosol dispersal both in the vertical and in the horizontal direction. In some models, a minimum of total aerosol concentration is simulated at the surface. Aerosol dispersal is most pronounced for SO4 and BC and lowest for SS. Diversities are higher for meridional than for vertical dispersal, they are similar for the individual species and highest for SS and DU. For these two components we do not find a correlation between vertical and meridional aerosol dispersal. In addition the degree of dispersals of SS and DU is not related to their residence times. SO4, BC, and POM, however, show increased meridional dispersal in models with larger vertical dispersal, and dispersal is larger for longer simulated residence times.