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⁻¹ (range 34–144 Tg a⁻¹) and the median SOA source strength (natural and anthropogenic) is 19 Tg a⁻¹ (range 13–121 Tg a⁻¹). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a⁻¹ (range 16–121 Tg a⁻¹), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a⁻¹; range 13–20 Tg a⁻¹, with one model at 37 Tg a⁻¹). 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⁻¹ (range 28–209 Tg a⁻¹), 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 aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately.

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AOD_f, AOD_c), Ångström exponent (AE), dry surface scattering (SC_dry), and absorption (AC_dry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AOD_c (∼ −45 % ± 25 %) than in fine AOD_f (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AOD_c bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AOD_c. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SC_dry and AC_dry, respectively. The larger underestimate of SC_dry than AC_dry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SC_dry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

We introduce and evaluate aerosol simulations with the global aerosol–climate model ECHAM6.3–HAM2.3, which is the aerosol component of the fully coupled aerosol–chemistry–climate model ECHAM–HAMMOZ. Both the host atmospheric climate model ECHAM6.3 and the aerosol model HAM2.3 were updated from previous versions. The updated version of the HAM aerosol model contains improved parameterizations of aerosol processes such as cloud activation, as well as updated emission fields for anthropogenic aerosol species and modifications in the online computation of sea salt and mineral dust aerosol emissions. Aerosol results from nudged and free-running simulations for the 10-year period 2003 to 2012 are compared to various measurements of aerosol properties. While there are regional deviations between the model and observations, the model performs well overall in terms of aerosol optical thickness, but may underestimate coarse-mode aerosol concentrations to some extent so that the modeled particles are smaller than indicated by the observations. Sulfate aerosol measurements in the US and Europe are reproduced well by the model, while carbonaceous aerosol species are biased low. Both mineral dust and sea salt aerosol concentrations are improved compared to previous versions of ECHAM–HAM. The evaluation of the simulated aerosol distributions serves as a basis for the suitability of the model for simulating aerosol–climate interactions in a changing climate.

Uncertainty in the representation of biomass burning (BB) aerosol composition and optical properties in climate models contributes to a range in modeled aerosol effects on incoming solar radiation. Depending on the model, the top-of-the-atmosphere BB aerosol effect can range from cooling to warming. By relating aerosol absorption relative to extinction and carbonaceous aerosol composition from 12 observational datasets to nine state-of-the-art Earth system models/chemical transport models, we identify varying degrees of overestimation in BB aerosol absorptivity by these models. Modifications to BB aerosol refractive index, size, and mixing state improve the Community Atmosphere Model version 5 (CAM5) agreement with observations, leading to a global change in BB direct radiative effect of −0.07 W m⁻², and regional changes of −2 W m⁻² (Africa) and −0.5 W m⁻² (South America/Temperate). Our findings suggest that current modeled BB contributes less to warming than previously thought, largely due to treatments of aerosol mixing state.

We analyse the development of a radiation fog event and its gradual transition from optically thin fog in a stable boundary layer to well-mixed optically thick fog. A comparison of observations and a detailed large-eddy simulation demonstrate that aerosol growth and activation is the key process in determining the onset of adiabatic fog. Weak turbulence and low supersaturations lead to the growth of aerosol particles which can significantly affect the visibility but do not significantly interact with the long-wave radiation, allowing the atmosphere to remain stable. Only when a substantial fraction of the aerosol activates into cloud droplets can the fog interact with the radiation, becoming optically thick and well mixed. Modifications to the parameterisation of cloud droplet numbers in fog, resulting in lower and more realistic concentrations, are shown to give significant improvements to an NWP model, which initially struggled to accurately simulate the transition. Finally, the consequences of this work for common aerosol activation parameterisations used in climate models are discussed, demonstrating that many schemes are reliant on an artificial minimum value when activating aerosol in fog, and adjustment of this minimum can significantly affect the sensitivity of the climate system to aerosol radiative forcing.

The global aerosol–climate model ECHAM6.3–HAM2.3 (E63H23) as well as the previous model versions ECHAM5.5–HAM2.0 (E55H20) and ECHAM6.1–HAM2.2 (E61H22) are evaluated using global observational datasets for clouds and precipitation. In E63H23, the amount of low clouds, the liquid and ice water path, and cloud radiative effects are more realistic than in previous model versions. E63H23 has a more physically based aerosol activation scheme, improvements in the cloud cover scheme, changes in the detrainment of convective clouds, changes in the sticking efficiency for the accretion of ice crystals by snow, consistent ice crystal shapes throughout the model, and changes in mixed-phase freezing; an inconsistency in ice crystal number concentration (ICNC) in cirrus clouds was also removed. Common biases in ECHAM and in E63H23 (and in previous ECHAM–HAM versions) are a cloud amount in stratocumulus regions that is too low and deep convective clouds over the Atlantic and Pacific oceans that form too close to the continents (while tropical land precipitation is underestimated). There are indications that ICNCs are overestimated in E63H23. Since clouds are important for effective radiative forcing due to aerosol–radiation and aerosol–cloud interactions (ERF_ari+aci) and equilibrium climate sensitivity (ECS), differences in ERF_ari+aci and ECS between the model versions were also analyzed. ERF_ari+aci is weaker in E63H23 (−1.0 W m⁻²) than in E61H22 (−1.2 W m⁻²) (or E55H20; −1.1 W m⁻²). This is caused by the weaker shortwave ERF_ari+aci (a new aerosol activation scheme and sea salt emission parameterization in E63H23, more realistic simulation of cloud water) overcompensating for the weaker longwave ERF_ari+aci (removal of an inconsistency in ICNC in cirrus clouds in E61H22). The decrease in ECS in E63H23 (2.5 K) compared to E61H22 (2.8 K) is due to changes in the entrainment rate for shallow convection (affecting the cloud amount feedback) and a stronger cloud phase feedback. Experiments with minimum cloud droplet number concentrations (CDNCmin) of 40 cm⁻³ or 10 cm⁻³ show that a higher value of CDNCmin reduces ERF_ari+aci as well as ECS in E63H23.