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

Present day tropospheric ozone and its changes between 1850 and 2100 are considered, analysing 15 global models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean compares well against present day observations. The seasonal cycle correlates well, except for some locations in the tropical upper troposphere. Most (75 %) of the models are encompassed with a range of global mean tropospheric ozone column estimates from satellite data, but there is a suggestion of a high bias in the Northern Hemisphere and a low bias in the Southern Hemisphere, which could indicate deficiencies with the ozone precursor emissions. Compared to the present day ensemble mean tropospheric ozone burden of 337 ± 23 Tg, the ensemble mean burden for 1850 time slice is ~30% lower. Future changes were modelled using emissions and climate projections from four Representative Concentration Pathways (RCPs). Compared to 2000, the relative changes in the ensemble mean tropospheric ozone burden in 2030 (2100) for the different RCPs are: −4% (−16%) for RCP2.6, 2% (−7%) for RCP4.5, 1% (−9%) for RCP6.0, and 7% (18%) for RCP8.5. Model agreement on the magnitude of the change is greatest for larger changes. Reductions in most precursor emissions are common across the RCPs and drive ozone decreases in all but RCP8.5, where doubled methane and a 40–150% greater stratospheric influx (estimated from a subset of models) increase ozone. While models with a high ozone burden for the present day also have high ozone burdens for the other time slices, no model consistently predicts large or small ozone changes; i.e. the magnitudes of the burdens and burden changes do not appear to be related simply, and the models are sensitive to emissions and climate changes in different ways. Spatial patterns of ozone changes are well correlated across most models, but are notably different for models without time evolving stratospheric ozone concentrations. A unified approach to ozone budget specifications and a rigorous investigation of the factors that drive tropospheric ozone is recommended to help future studies attribute ozone changes and inter-model differences more clearly

Abstract. The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) examined the short-lived drivers of climate change in current climate models. Here we evaluate the 10 ACCMIP models that included aerosols, 8 of which also participated in the Coupled Model Intercomparison Project phase 5 (CMIP5). The models reproduce present-day total aerosol optical depth (AOD) relatively well, though many are biased low. Contributions from individual aerosol components are quite different, however, and most models underestimate east Asian AOD. The models capture most 1980–2000 AOD trends well, but underpredict increases over the Yellow/Eastern Sea. They strongly underestimate absorbing AOD in many regions. We examine both the direct radiative forcing (RF) and the forcing including rapid adjustments (effective radiative forcing; ERF, including direct and indirect effects). The models' all-sky 1850 to 2000 global mean annual average total aerosol RF is (mean; range) −0.26 W m−2; −0.06 to −0.49 W m−2. Screening based on model skill in capturing observed AOD yields a best estimate of −0.42 W m−2; −0.33 to −0.50 W m−2, including adjustment for missing aerosol components in some models. Many ACCMIP and CMIP5 models appear to produce substantially smaller aerosol RF than this best estimate. Climate feedbacks contribute substantially (35 to −58%) to modeled historical aerosol RF. The 1850 to 2000 aerosol ERF is −1.17 W m−2; −0.71 to −1.44 W m−2. Thus adjustments, including clouds, typically cause greater forcing than direct RF. Despite this, the multi-model spread relative to the mean is typically the same for ERF as it is for RF, or even smaller, over areas with substantial forcing. The largest 1850 to 2000 negative aerosol RF and ERF values are over and near Europe, south and east Asia and North America. ERF, however, is positive over the Sahara, the Karakoram, high Southern latitudes and especially the Arctic. Global aerosol RF peaks in most models around 1980, declining thereafter with only weak sensitivity to the Representative Concentration Pathway (RCP). One model, however, projects approximately stable RF levels, while two show increasingly negative RF due to nitrate (not included in most models). Aerosol ERF, in contrast, becomes more negative during 1980 to 2000. During this period, increased Asian emissions appear to have a larger impact on aerosol ERF than European and North American decreases due to their being upwind of the large, relatively pristine Pacific Ocean. There is no clear relationship between historical aerosol ERF and climate sensitivity in the CMIP5 subset of ACCMIP models. In the ACCMIP/CMIP5 models, historical aerosol ERF of about −0.8 to −1.5 W m−2 is most consistent with observed historical warming. Aerosol ERF masks a large portion of greenhouse forcing during the late 20th and early 21st century at the global scale. Regionally, aerosol ERF is so large that net forcing is negative over most industrialized and biomass burning regions through 1980, but remains strongly negative only over east and southeast Asia by 2000. Net forcing is strongly positive by 1980 over most deserts, the Arctic, Australia, and most tropical oceans. Both the magnitude of and area covered by positive forcing expand steadily thereafter.

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

Ozone (O₃) from 17 atmospheric chemistry models taking part in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) has been used to calculate tropospheric ozone radiative forcings (RFs). All models applied a common set of anthropogenic emissions, which are better constrained for the present-day than the past. Future anthropogenic emissions follow the four Representative Concentration Pathway (RCP) scenarios, which define a relatively narrow range of possible air pollution emissions. We calculate a value for the pre-industrial (1750) to present-day (2010) tropospheric ozone RF of 410 mW m⁻². The model range of pre-industrial to present-day changes in O₃ produces a spread (±1 standard deviation) in RFs of ±17%. Three different radiation schemes were used – we find differences in RFs between schemes (for the same ozone fields) of ±10%. Applying two different tropopause definitions gives differences in RFs of ±3%. Given additional (unquantified) uncertainties associated with emissions, climate-chemistry interactions and land-use change, we estimate an overall uncertainty of ±30% for the tropospheric ozone RF. Experiments carried out by a subset of six models attribute tropospheric ozone RF to increased emissions of methane (44±12%), nitrogen oxides (31 ± 9%), carbon monoxide (15 ± 3%) and non-methane volatile organic compounds (9 ± 2%); earlier studies attributed more of the tropospheric ozone RF to methane and less to nitrogen oxides. Normalising RFs to changes in tropospheric column ozone, we find a global mean normalised RF of 42 mW m⁻² DU⁻¹, a value similar to previous work. Using normalised RFs and future tropospheric column ozone projections we calculate future tropospheric ozone RFs (mW m⁻²; relative to 1750) for the four future scenarios (RCP2.6, RCP4.5, RCP6.0 and RCP8.5) of 350, 420, 370 and 460 (in 2030), and 200, 300, 280 and 600 (in 2100). Models show some coherent responses of ozone to climate change: decreases in the tropical lower troposphere, associated with increases in water vapour; and increases in the sub-tropical to mid-latitude upper troposphere, associated with increases in lightning and stratosphere-to-troposphere transport. Climate change has relatively small impacts on global mean tropospheric ozone RF