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

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

Abstract. The core version of the Norwegian Climate Center's Earth System Model, named NorESM1-M, is presented. The NorESM family of models are based on the Community Climate System Model version 4 (CCSM4) of the University Corporation for Atmospheric Research, but differs from the latter by, in particular, an isopycnic coordinate ocean model and advanced chemistry-aerosol-cloud-radiation interaction schemes. NorESM1-M has a horizontal resolution of approximately 2 • for the atmosphere and land components and 1 • for the ocean and ice components. NorESM is also available in a lower resolution version (NorESM1-L) and a version that includes prognostic biogeochemical cycling (NorESM1-ME). The latter two model configurations are not part of this paper. Here, a first-order assessment of the model stability, the mean model state and the internal variability based on the model experiments made available to CMIP5 are presented. Further analysis of the model performance is provided in an accompanying paper , presenting the corresponding climate response and scenario projections made with NorESM1-M.

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.

Published by Copernicus Publications on behalf of the European Geosciences Union. 8698 J. Quaas et al.: Aerosol indirect effects -general circulation model intercomparisonbetween τ a and f cld . 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 N d . 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 clearand cloudy-sky forcings with estimates of anthropogenic τ a and satellite-retrieved N d -τ 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 .

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 agi...

Precipitation is expected to respond differently to various drivers of anthropogenic climate change. We present the first results from the Precipitation Driver and Response Model Intercomparison Project (PDRMIP), where nine global climate models have perturbed CO2, CH4, black carbon, sulfate, and solar insolation. We divide the resulting changes to global mean and regional precipitation into fast responses that scale with changes in atmospheric absorption and slow responses scaling with surface temperature change. While the overall features are broadly similar between models, we find significant regional intermodel variability, especially over land. Black carbon stands out as a component that may cause significant model diversity in predicted precipitation change. Processes linked to atmospheric absorption are less consistently modeled than those linked to top‐of‐atmosphere radiative forcing. We identify a number of land regions where the model ensemble consistently predicts that fast precipitation responses to climate perturbations dominate over the slow, temperature‐driven responses.