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

Aerosol variations and trends over different land and ocean regions from 1980 to 2009 are analyzed with the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and observations from multiple satellite sensors and available ground-based networks. Excluding time periods with large volcanic influence, aerosol optical depth (AOD) and surface concentration over polluted land regions generally vary with anthropogenic emissions, but the magnitude of this association can be dampened by the presence of natural aerosols, especially dust. Over the 30-year period in this study, the largest reduction in aerosol levels occurs over Europe, where AOD has decreased by 40–60% on average and surface sulfate concentrations have declined by a factor of up to 3–4. In contrast, East Asia and South Asia show AOD increases, but the relatively high level of dust aerosols in Asia reduces the correlation between AOD and pollutant emission trends. Over major dust source regions, model analysis indicates that the change of dust emissions over the Sahara and Sahel has been predominantly driven by the change of near-surface wind speed, but over Central Asia it has been largely influenced by the change of the surface wetness. The decreasing dust trend in the North African dust outflow region of the tropical North Atlantic and the receptor sites of Barbados and Miami is closely associated with an increase of the sea surface temperature in the North Atlantic. This temperature increase may drive the decrease of the wind velocity over North Africa, which reduces the dust emission, and the increase of precipitation over the tropical North Atlantic, which enhances dust removal during transport. Despite significant trends over some major continental source regions, the model-calculated global annual average AOD shows little change over land and ocean in the past three decades, because opposite trends in different land regions cancel each other out in the global average, and changes over large open oceans are negligible. This highlights the necessity for regional-scale assessment of aerosols and their climate impacts, as global-scale average values can obscure important regional changes.

An assessment of global particulate nitrate and ammonium aerosol based on simulations from nine models participating in the Aerosol Comparisons between Observations and Models (AeroCom) phase III study is presented. A budget analysis was conducted to understand the typical magnitude, distribution, and diversity of the aerosols and their precursors among the models. To gain confidence regarding model performance, the model results were evaluated with various observations globally, including ground station measurements over North America, Europe, and east Asia for tracer concentrations and dry and wet depositions, as well as with aircraft measurements in the Northern Hemisphere mid-to-high latitudes for tracer vertical distributions. Given the unique chemical and physical features of the nitrate occurrence, we further investigated the similarity and differentiation among the models by examining (1) the pH-dependent NH ₃ wet deposition; (2) the nitrate formation via heterogeneous chemistry on the surface of dust and sea salt particles or thermodynamic equilibrium calculation including dust and sea salt ions; and (3) the nitrate coarse-mode fraction (i.e., coarse/total). It is found that HNO ₃ , which is simulated explicitly based on full O ₃ -HO _x -NO _x -aerosol chemistry by all models, differs by up to a factor of 9 among the models in its global tropospheric burden. This partially contributes to a large difference in NO ₃ ⁻ , whose atmospheric burden differs by up to a factor of 13. The atmospheric burdens of NH ₃ and NH ₄ ⁺ differ by 17 and 4, respectively. Analyses at the process level show that the large diversity in atmospheric burdens of NO ₃ ⁻ , NH ₃ , and NH ₄ ⁺ is also related to deposition processes. Wet deposition seems to be the dominant process in determining the diversity in NH ₃ and NH ₄ ⁺ lifetimes. It is critical to correctly account for contributions of heterogeneous chemical production of nitrate on dust and sea salt, because this process overwhelmingly controls atmospheric nitrate production (typically > 80 %) and determines the coarse- and fine-mode distribution of nitrate aerosol.

Simulated multi-model "diversity" in aerosol direct radiative forcing estimates is often perceived as a measure of aerosol uncertainty. However, current models used for aerosol radiative forcing calculations vary considerably in model components relevant for forcing calculations and the associated "host-model uncertainties" are generally convoluted with the actual aerosol uncertainty. In this AeroCom Prescribed intercomparison study we systematically isolate and quantify host model uncertainties on aerosol forcing experiments through prescription of identical aerosol radiative properties in twelve participating models. Even with prescribed aerosol radiative properties, simulated clear-sky and all-sky aerosol radiative forcings show significant diversity. For a purely scattering case with globally constant optical depth of 0.2, the global-mean all-sky top-of-atmosphere radiative forcing is −4.47 Wm⁻² and the inter-model standard deviation is 0.55 Wm⁻², corresponding to a relative standard deviation of 12%. For a case with partially absorbing aerosol with an aerosol optical depth of 0.2 and single scattering albedo of 0.8, the forcing changes to 1.04 Wm⁻², and the standard deviation increases to 1.01 W⁻², corresponding to a significant relative standard deviation of 97%. However, the top-of-atmosphere forcing variability owing to absorption (subtracting the scattering case from the case with scattering and absorption) is low, with absolute (relative) standard deviations of 0.45 Wm⁻² (8%) clear-sky and 0.62 Wm⁻² (11%) all-sky. Scaling the forcing standard deviation for a purely scattering case to match the sulfate radiative forcing in the AeroCom Direct Effect experiment demonstrates that host model uncertainties could explain about 36% of the overall sulfate forcing diversity of 0.11 Wm⁻² in the AeroCom Direct Radiative Effect experiment. Host model errors in aerosol radiative forcing are largest in regions of uncertain host model components, such as stratocumulus cloud decks or areas with poorly constrained surface albedos, such as sea ice. Our results demonstrate that host model uncertainties are an important component of aerosol forcing uncertainty that require further attention.

Atmospheric black carbon (BC) absorbs solar radiation, and exacerbates global warming through exerting positive radiative forcing (RF). However, the contribution of BC to ongoing changes in global climate is under debate. Anthropogenic BC emissions, and the resulting distribution of BC concentration, are highly uncertain. In particular, long-range transport and processes affecting BC atmospheric lifetime are poorly understood. Here we discuss whether recent assessments may have overestimated present-day BC radiative forcing in remote regions. We compare vertical profiles of BC concentration from four recent aircraft measurement campaigns to simulations by 13 aerosol models participating in the AeroCom Phase II intercomparison. An atmospheric lifetime of BC of less than 5 days is shown to be essential for reproducing observations in remote ocean regions, in line with other recent studies. Adjusting model results to measurements in remote regions, and at high altitudes, leads to a 25% reduction in AeroCom Phase II median direct BC forcing, from fossil fuel and biofuel burning, over the industrial era. The sensitivity of modelled forcing to BC vertical profile and lifetime highlights an urgent need for further flight campaigns, close to sources and in remote regions, to provide improved quantification of BC effects for use in climate policy.