Abstract. Many of the next generation of global climate models will include aerosol schemes which explicitly simulate the microphysical processes that determine the particle size distribution. These models enable aerosol optical properties and cloud condensation nuclei (CCN) concentrations to be determined by fundamental aerosol processes, which should lead to a more physically based simulation of aerosol direct and indirect radiative forcings. This study examines the global variation in particle size distribution simulated by twelve global aerosol microphysics models to quantify model diversity and to identify any common biases against observations. Evaluation against size distribution measurements from a new European network of aerosol supersites shows that the mean model agrees quite well with the observations at many sites on the annual mean, but there are some seasonal biases common to many sites. In particular, at many of these European sites, the accumulation mode number concentration is biased low during winter and Aitken mode concentrations tend to be overestimated in winter and underestimated in summer. At high northern latitudes, the models strongly underpredict Aitken and accumulation particle concentrations compared to the measurements, consistent with previous studies that have highlighted the poor performance of global aerosol models in the Arctic. In the marine boundary layer, the models capture the observed meridional variation in the size distribution, which is dominated by the Aitken mode at high latitudes, with an increasing concentration of accumulation particles with decreasing latitude. Considering vertical profiles, the models reproduce the observed peak in total particle concentrations in the upper troposphere due to new particle formation, although modelled peak concentrations tend to be biased high over Europe. Overall, the multi-model-mean dataset simulates the global variation of the particle size distribution with a good degree of skill, suggesting that most of the individual global aerosol microphysics models are performing well, although the large model diversity indicates that some models are in poor agreement with the observations. Further work is required to better constrain size-resolved primary and secondary particle number sources, and an improved understanding of nucleation and growth (e.g. the role of nitrate and secondary organics) will improve the fidelity of simulated particle size distributions.

The Hadley Centre Global Environmental Model (HadGEM) includes two aerosol schemes: the Coupled Large-scale Aerosol Simulator for Studies in Climate (CLASSIC), and the new Global Model of Aerosol Processes (GLOMAP-mode). GLOMAP-mode is a modal aerosol microphysics scheme that simulates not only aerosol mass but also aerosol number, represents internally-mixed particles, and includes aerosol microphysical processes such as nucleation. In this study, both schemes provide hindcast simulations of natural and anthropogenic aerosol species for the period 2000–2006. HadGEM simulations using GLOMAP-mode compare better than CLASSIC against a data-assimilated aerosol re-analysis and aerosol ground-based observations. GLOMAP-mode sulphate aerosol residence time is two days longer than CLASSIC sulphate aerosols, whereas black carbon residence time is much shorter. As a result, CLASSIC underestimates aerosol optical depths in continental regions of the Northern Hemisphere and likely overestimates absorption in remote regions. Aerosol direct and first indirect radiative forcings are computed from simulations of aerosols with emissions for the year 1850 and 2000. In 1850, GLOMAP-mode predicts lower aerosol optical depths and higher cloud droplet number concentrations than CLASSIC. Consequently, simulated clouds are much less susceptible to natural and anthropogenic aerosol changes when the microphysical scheme is used. In particular, the response of cloud condensation nuclei to an increase in dimethyl sulphide emissions becomes a factor of four smaller. The combined effect of different 1850 baselines, residence times, and cloud susceptibilities, leads to substantial differences in the aerosol forcings simulated by the two schemes. GLOMAP-mode finds a present-day direct aerosol forcing of −0.49 W m^−2 on a global average, 72% stronger than the corresponding forcing from CLASSIC. This difference is compensated by changes in first indirect aerosol forcing: the forcing of −1.17 W m^−2 obtained with GLOMAP-mode is 20% weaker than with CLASSIC. Results suggest that mass-based schemes such as CLASSIC lack the necessary sophistication to provide realistic input to aerosol-cloud interaction schemes. Furthermore, the importance of the 1850 baseline highlights how model skill in predicting present-day aerosol does not guarantee reliable forcing estimates. Those findings suggest that the more complex representation of aerosol processes in microphysical schemes improves the fidelity of simulated aerosol forcings

Regional patterns of aerosol radiative forcing are important for understanding climate change on decadal time scales. Uncertainty in aerosol forcing is likely to vary regionally and seasonally because of the short aerosol lifetime and heterogeneous emissions. Here the sensitivity of regional aerosol cloud albedo effect (CAE) forcing to 31 aerosol process parameters and emission fluxes is quantified between 1978 and 2008. The effects of parametric uncertainties on calculations of the balance of incoming and outgoing radiation are found to be spatially and temporally dependent. Regional uncertainty contributions of opposite sign cancel in global-mean forcing calculations, masking the regional importance of some parameters. Parameters that contribute little to uncertainty in Earth's global energy balance during recent decades make significant contributions to regional forcing variance. Aerosol forcing sensitivities are quantified within 11 climatically important regions, where surface temperatures are thought to influence large-scale climate effects. Substantial simulated uncertainty in CAE forcing in the eastern Pacific leaves open the possibility that apparent shifts in the mean ENSO state may result from a forced aerosol signal on multidecadal time scales. A likely negative aerosol CAE forcing in the tropical North Atlantic calls into question the relationship between Northern Hemisphere aerosol emission reductions and CAE forcing of sea surface temperatures in the main Atlantic hurricane development region on decadal time scales. Simulated CAE forcing uncertainty is large in the North Pacific, suggesting that the role of the CAE in altering Pacific tropical storm frequency and intensity is also highly uncertain.

A global modal aerosol microphysics module (GLOMAP-mode) is evaluated and improved by comparing against a sectional version (GLOMAP-bin) and observations in the same 3-D global offline chemistry transport model. With both schemes, the model captures the main features of the global particle size distribution, with sub-micron aerosol approximately unimodal in continental regions and bi-modal in marine regions. Initial bin-mode comparisons showed that various size distribution parameter settings (mode widths and inter-modal separation sizes) resulted in clear biases compared to the sectional scheme. By adjusting these parameters in the modal scheme, much better agreement is achieved against the bin scheme and observations. Surface mass of sulphate, sea-salt, black carbon (BC) and organic carbon (OC) are, on the annual mean, within 25 % in the two schemes in nearly all regions. On the annual mean, surface level concentrations of condensation nuclei (CN), cloud condensation nuclei (CCN), surface area density and condensation sink also compare within 25 % in most regions. However, marine CCN concentrations between 30° N and 30° S are systematically higher in the modal scheme, by 25–60 %, which we attribute to differences in size-resolved particle growth or cloud-processing. Larger differences also exist in regions or seasons dominated by biomass burning and in free-troposphere and high-latitude regions. Indeed, in the free-troposphere, GLOMAP-mode BC is a factor 2–4 higher than GLOMAP-bin, likely due to differences in size-resolved scavenging. Nevertheless, in most parts of the atmosphere, we conclude that bin-mode differences are much less than model-observation differences, although some processes are missing in these runs which may pose a bigger challenge to modal schemes (e.g. boundary layer nucleation, ultra-fine sea-spray). The findings here underline the need for a spectrum of complexity in global models, with size-resolved aerosol properties predicted by modal schemes needing to be continually benchmarked and improved against freely evolving sectional schemes and observations

Abstract. The emission of dimethylsulphide (DMS) gas by phytoplankton and the subsequent formation of aerosol has long been suggested as an important climate regulation mechanism. The key aerosol quantity is the number concentration of cloud condensation nuclei (CCN), but until recently global models did not include the necessary aerosol physics to quantify CCN. Here we use a global aerosol microphysics model to calculate the sensitivity of CCN to changes in DMS emission using multiple present-day and future sea-surface DMS climatologies. Calculated annual fluxes of DMS to the atmosphere for the five model-derived and one observations based present day climatologies are in the range 15.1 to 32.3 Tg a−1 sulphur. The impact of DMS climatology on surface level CCN concentrations was calculated in terms of summer and winter hemispheric mean values of ΔCCN/ΔFluxDMS, which varied between −51 and +147 cm−3/(mg m−2 day−1 sulphur), with a mean of 56 cm−3/(mg m−2 day−1 sulphur). The range is due to CCN production in the atmosphere being strongly dependent on the spatial distribution of the emitted DMS. The DMS flux from a future globally warmed climatology was 0.2 Tg a−1 sulphur higher than present day with a mean CCN response of 95 cm−3/(mg m−2 day−1 sulphur) relative to present day. The largest CCN response was seen in the southern Ocean, contributing to a Southern Hemisphere mean annual increase of less than 0.2%. We show that the changes in DMS flux and CCN concentration between the present day and global warming scenario are similar to interannual differences due to variability in windspeed. In summary, although DMS makes a significant contribution to global marine CCN concentrations, the sensitivity of CCN to potential future changes in DMS flux is very low. This finding, together with the predicted small changes in future seawater DMS concentrations, suggests that the role of DMS in climate regulation is very weak.