Abstract. 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 (SO 4 ), particulate organic matter (POM), and finally black carbon (BC). Interactive parameterizations of the emissions and contrasting particles sizes of SS and DU lead genCorrespondence to: C. Textor (christiane.textor@cea.fr) erally to higher diversities of these species, and for total aerosol. The lower diversity of the emissions of the fine aerosols, BC, POM, and SO 4 , is due to the use of similar emission inventories, and does therefore not necessarily indicate a better understanding of their sources. The diversity of SO 4 -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, SO 4 , POM, and BC.The all-models-average residence time is shortest for SS with about half a day, followed by SO 4 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 SO 4 and SS. It is the dominant sink for SO 4 , BC, and POM, and contributes about one third to the total removal of SS and DU species. For SS Published by Copernicus GmbH on behalf of the European Geosciences Union. C. Textor et al.: Diversities of aerosol life cycles within AeroComand 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 SO 4 and BC and lowest for SS. Diversities are higher for meridional than for verti...

Abstract. Inventories for global aerosol and aerosol precursor emissions have been collected (based on published inventories and published simulations), assessed and prepared for the year 2000 (present-day conditions) and for the year 1750 (pre-industrial conditions). These global datasets establish a comprehensive source for emission input to global modeling, when simulating the aerosol impact on climate with state-of-the-art aerosol component modules. As these modules stratify aerosol into dust, sea-salt, sulfate, organic matter and soot, for all these aerosol types global fields on emission strength and recommendations for injection altitude and particulate size are provided. Temporal resolution varies between daily (dust and sea-salt), monthly (wild-land fires) and annual (all other emissions). These datasets benchmark aerosol emissions according to the knowledge in the year 2004. They are intended to serve as systematic constraints in sensitivity studies of the AeroCom initiative, which seeks to quantify (actual) uncertainties in aerosol global modeling.

This study presents the results of a broad intercomparison of a total of 15 global aerosol models within the AeroCom project. Each model is compared to observations related to desert dust aerosols, their direct radiative effect, and their impact on the biogeochemical cycle, i.e., aerosol optical depth (AOD) and dust deposition. Additional com parisons to Angstrom exponent (AE), coarse mode AOD and dust surface concentrations are included to extend the assessment of model performance and to identify common biases present in models. These data comprise a benchmark dataset that is proposed for model inspection and future dust model development. There are large differences among the global models that simulate the dust cycle and its impact on climate. In general, models simulate the climatology of vertically integrated parameters (AOD and AE) within a factor of two whereas the total deposition and surface concentration are reproduced within a factor of 10. In addition, smaller m! ean normalized bias and root mean square errors are obtained for the climatology of AOD and AE than for total deposition and surface concentration. Characteristics of the datasets used and their uncertainties may influence these differences. Large uncertainties still exist with respect to the deposition fluxes in the southern oceans. Further measurements and model studies are necessary to assess the general model performance to reproduce dust deposition in ocean regions sensible to iron contributions. Models overestimate the wet deposition in regions dominated by dry deposition. They generally simulate more realistic surface concentration at stations downwind of the main sources than at remote ones. Most models simulate the gradient in AOD and AE between the different dusty regions. However the seasonality and magnitude of both variables is better simulated at African stations than Middle East ones. The models simulate the offshore transport of West Africa throughout the year but they overestimate the AOD and they transport too fine particles . The models also reproduce the dust transport across the Atlantic in the summer in terms of both AOD and AE but not so well in winter-spring nor the southward displacement of the dust cloud that is responsible of the dust transport into South America. Based on the dependency of AOD on aerosol burden and size distribution we use model bias with respect to AOD and AE to infer the bias of the dust emissions in Africa and the Middle East. According to this analysis we suggest that a range of possible emissions for North Africa is 400 to 2200 Tg yr(-1) and in the Middle East 26 to 526 Tg yr(-1

We have developed detailed emission inventories for the amount of both black and organic carbon particles from biomass burning sources (wood fuel, charcoal burning, dung, charcoal production, agricultural, savanna and forest fires). We have also estimated an inventory for organic carbon particles from fossil fuel burning and urban activities from an existing inventory for fossil fuel sources of black carbon. We also provide an estimate for the natural source of organic matter. These emissions have been used together with our global aerosol model to study the global distribution of carbonaceous aerosols. The accuracy of the inventories and the model formulation has been tested by comparing the model simulations of carbonaceous aerosols in the atmosphere and in precipitation with observations reported in the literature. For most locations and seasons, the predicted concentrations are in reasonable agreement with the observations, although the model underpredicts black carbon concentrations in polar regions. The predicted concentrations in remote areas are extremely sensitive to both the rate of removal by wet deposition and the height of injection of the aerosols. Finally, a global map of the aerosol single scattering albedo was developed from the simulated carbonaceous particle distribution and a previously developed model for aerosol sulfates. The computed aerosol single scattering albedos compare well with observations, suggesting that most of the important aerosol species have been included in the model. For most locations and seasons, the single scattering albedo is larger than 0.85, indicating that these aerosols, in general, lead to a net cooling.

For CG energies we adopt a production rate of 10 x 1016 molecules NO/J based on the current literature. Using a method to simulate global lightning frequencies from satelliteobserved cloud data, we have calculated the LNOx on various spatial (regional, zonal, meridional, and global) and temporal scales (daily, monthly, seasonal, and interannual). Regionally, the production of LNO• is concentrated over tropical continental regions, predominantly in the summer hemisphere. The annual mean production rate is calculated to be 12.2 Tg N/yr, and we believe it extremely unlikely that this number is less than 5 or more than 20 Tg N/yr. Although most of LNO• is produced in the lowest 5 km by CG lightning, convective mixing in the thunderstorms is likely to deposit large amounts of NO• in the upper troposphere where it is important in ozone production. On an annual basis, 64% of the LNO• is produced in the northern hemisphere, implying that the northern hemisphere should have natural ozone levels as much as 2 times greater than the southern hemisphere, even before anthropogenic influences. The amount of 03 produced from this NOx is expected to exceed the stratospheric source by a factor of 1.5, and thus the hemispheric asymmetry in LNOx would lead to a significant excess of northern hemisphere 03 even in the preindustrial troposphere. (The monthly climatologies for LNO• on a 1 ø x 1 ø latitude-longitude grid can be obtained by e-mail to cprice@flash.tau.ac.il.)

Abstract. Nine different global models with detailed aerosol modules have independently produced instantaneous direct radiative forcing due to anthropogenic aerosols. The anthropogenic impact is derived from the difference of two model simulations with prescribed aerosol emissions, one for present-day and one for pre-industrial conditions. The difference in the solar energy budget at the top of the atmosphere (ToA) yields a new harmonized estimate for the aerosol direct radiative forcing (RF) under all-sky conditions. On a global annual basis RF is −0.22 Wm −2 , ranging from +0.04 to −0.41 Wm −2 , with a standard deviation of ±0.16 Wm −2 . Anthropogenic nitrate and dust are not included in this estimate. No model shows a significant positive all-sky RF. The corresponding clear-sky RF is −0.68 Wm −2 . The cloud-sky RF was derived based on all-sky and clear-sky RF and modelled cloud cover. It was significantly different from zero and ranged between −0.16 and +0.34 Wm −2 . A sensitivity analysis shows that the total aerosol RF is influenced by considerable diversity in simulated residence times, mass extinction coefficients and most importantly forcing efficiencies (forcing per unit optical depth). The clear-sky forcing efficiency (forcing per unit optical depth) has diversity comparable to that for the all-sky/ clear-sky forcing ratio. While the diversity in clear-sky forcing efficiency is impacted by factors Correspondence to: M. Schulz (michael.schulz@cea.fr) such as aerosol absorption, size, and surface albedo, we can show that the all-sky/clear-sky forcing ratio is important because all-sky forcing estimates require proper representation of cloud fields and the correct relative altitude placement between absorbing aerosol and clouds. The analysis of the sulphate RF shows that long sulphate residence times are compensated by low mass extinction coefficients and vice versa. This is explained by more sulphate particle humidity growth and thus higher extinction in those models where short-lived sulphate is present at lower altitude and vice versa. Solar atmospheric forcing within the atmospheric column is estimated at +0.82±0.17 Wm −2 . The local annual average maxima of atmospheric forcing exceed +5 Wm −2 confirming the regional character of aerosol impacts on climate. The annual average surface forcing is −1.02±0.23 Wm −2 . With the current uncertainties in the modelling of the radiative forcing due to the direct aerosol effect we show here that an estimate from one model is not sufficient but a combination of several model estimates is necessary to provide a mean and to explore the uncertainty.

Abstract. We evaluate black carbon (BC) model predictions from the AeroCom model intercomparison project by considering the diversity among year 2000 model simulations and comparing model predictions with available measurements. These model-measurement intercomparisons include BC surface and aircraft concentrations, aerosol absorption optical depth (AAOD) retrievals from AERONET and Ozone Monitoring Instrument (OMI) and BC column estimations based on AERONET. In regions other than Asia, most models are biased high compared to surface concentration measurements. However compared with (column) AAOD or BC burden retreivals, the models are generally biased low. The average ratio of model to retrieved AAOD is less than 0.7 in South American and 0.6 in African biomass burning regions; both of these regions lack surface concentration measurements. In Asia the average model to observed ratio is 0.7 for AAOD and 0.5 for BC surface concentrations. Compared with aircraft measurements over the Americas at latitudes between 0 and 50N, the average model is a factor of 8 larger than observed, and most models exceed the measured BC standard deviation in the mid to upper troposphere. At higher latitudes the average model to aircraft BC ratio is 0.4 and models underestimate the observed BC loading in the lower and middle troposphere associated with springtime Arctic haze. Low model bias for AAOD but overestimation of surface and upper atmospheric BC concentrations at lower latitudes suggests that most models are underestimating BC absorption and should improve estimates for refractive index, particle size, and optical effects of BC coating. Retrieval uncertainties and/or differences with model diagnostic treatment may also contribute to the model-measurement disparity. Largest AeroCom model diversity occurred in northern Eurasia and the remote Arctic, regions influenced by anthropogenic sources. Changing emissions, aging, removal, or optical properties within a single model generated a smaller change in model predictions than the range represented by the full set of AeroCom models. Upper tropospheric concentrations of BC mass from the aircraft measurements are suggested to provide a unique new benchmark to test scavenging and vertical dispersion of BC in global models.

Global aviation operations contribute to anthropogenic climate change via a complex set of processes that lead to a net surface warming. Of importance are aviation emissions of carbon dioxide (CO 2 ), nitrogen oxides (NO x ), water vapor, soot and sulfate aerosols, and increased cloudiness due to contrail formation. Aviation grew strongly over the past decades (1960–2018) in terms of activity, with revenue passenger kilometers increasing from 109 to 8269 billion km yr −1 , and in terms of climate change impacts, with CO 2 emissions increasing by a factor of 6.8–1034 Tg CO 2 yr −1 . Over the period 2013–2018, the growth rates in both terms show a marked increase. Here, we present a new comprehensive and quantitative approach for evaluating aviation climate forcing terms. Both radiative forcing (RF) and effective radiative forcing (ERF) terms and their sums are calculated for the years 2000–2018. Contrail cirrus, consisting of linear contrails and the cirrus cloudiness arising from them, yields the largest positive net (warming) ERF term followed by CO 2 and NO x emissions. The formation and emission of sulfate aerosol yields a negative (cooling) term. The mean contrail cirrus ERF/RF ratio of 0.42 indicates that contrail cirrus is less effective in surface warming than other terms. For 2018 the net aviation ERF is +100.9 mW (mW) m −2 (5–95% likelihood range of (55, 145)) with major contributions from contrail cirrus (57.4 mW m −2 ), CO 2 (34.3 mW m −2 ), and NO x (17.5 mW m −2 ). Non-CO 2 terms sum to yield a net positive (warming) ERF that accounts for more than half (66%) of the aviation net ERF in 2018. Using normalization to aviation fuel use, the contribution of global aviation in 2011 was calculated to be 3.5 (4.0, 3.4) % of the net anthropogenic ERF of 2290 (1130, 3330) mW m −2 . Uncertainty distributions (5%, 95%) show that non-CO 2 forcing terms contribute about 8 times more than CO 2 to the uncertainty in the aviation net ERF in 2018. The best estimates of the ERFs from aviation aerosol-cloud interactions for soot and sulfate remain undetermined. CO 2 -warming-equivalent emissions based on global warming potentials (GWP* method) indicate that aviation emissions are currently warming the climate at approximately three times the rate of that associated with aviation CO 2 emissions alone. CO 2 and NO x aviation emissions and cloud effects remain a continued focus of anthropogenic climate change research and policy discussions.