Search

Despite a large number of studies, out of all drivers of radiative forcing, the effect of aerosols has the largest uncertainty in global climate model radiative forcing estimates. There have been studies of aerosol optical properties in climate models, but the effects of particle number size distribution need a more thorough inspection. We investigated the trends and seasonality of particle number concentrations in nucleation, Aitken, and accumulation modes at 21 measurement sites in Europe and the Arctic. For 13 of those sites, with longer measurement time series, we compared the field observations with the results from five climate models, namely EC-Earth3, ECHAM-M7, ECHAM-SALSA, NorESM1.2, and UKESM1. This is the first extensive comparison of detailed aerosol size distribution trends between in situ observations from Europe and five earth system models (ESMs). We found that the trends of particle number concentrations were mostly consistent and decreasing in both measurements and models. However, for many sites, climate models showed weaker decreasing trends than the measurements. Seasonal variability in measured number concentrations, quantified by the ratio between maximum and minimum monthly number concentration, was typically stronger at northern measurement sites compared to other locations. Models had large differences in their seasonal representation, and they can be roughly divided into two categories: for EC-Earth and NorESM, the seasonal cycle was relatively similar for all sites, and for other models the pattern of seasonality varied between northern and southern sites. In addition, the variability in concentrations across sites varied between models, some having relatively similar concentrations for all sites, whereas others showed clear differences in concentrations between remote and urban sites. To conclude, although all of the model simulations had identical input data to describe anthropogenic mass emissions, trends in differently sized particles vary among the models due to assumptions in emission sizes and differences in how models treat size-dependent aerosol processes. The inter-model variability was largest in the accumulation mode, i.e. sizes which have implications for aerosol–cloud interactions. Our analysis also indicates that between models there is a large variation in efficiency of long-range transportation of aerosols to remote locations. The differences in model results are most likely due to the more complex effect of different processes instead of one specific feature (e.g. the representation of aerosol or emission size distributions). Hence, a more detailed characterization of microphysical processes and deposition processes affecting the long-range transport is needed to understand the model variability.

A mathematical model called BUOYANT has previously been developed for the evaluation of the dispersion of positively buoyant plumes originating from major warehouse fires. The model addresses the variations of the cross-plume integrated properties (i.e., the average properties along a trajectory) of a rising plume in a vertically varying atmosphere and the atmospheric dispersion after the plume rise regime. We have described in this article an extension of the BUOYANT model to include a detailed treatment of the early evolution of the fire plumes before the plume rise and atmospheric dispersion regimes. The input and output of the new module consist of selected characteristics of forest or pool fires and the properties of a source term for the plume rise module, respectively. The main structure of this source term module is based on the differential equations for low-momentum releases of buoyant material, which govern the evolution of the plume radius, as well as velocity and density differences. The source term module is also partially based on various experimental results on fire plumes. We have evaluated the refined BUOYANT model by comparing the model predictions against the experimental field-scale data from the Prescribed Fire Combustion and Atmospheric Dynamics Research Experiment, RxCADRE. The predicted concentrations of CO2 agreed fairly well with the aircraft measurements conducted in the RxCADRE campaign. We have also compiled an operational version of the model. The operational model can be used for emergency contingency planning and the training of emergency personnel in case of major forest and pool fires.

The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.

We have implemented and evaluated a secondary organic aerosol scheme within the chemistry transport model TM5-MP in this work. In earlier versions of TM5-MP the secondary organic aerosol (SOA) was emitted as Aitken-sized particle mass emulating the condensation. In the current scheme we simulate the formation of secondary organic aerosol from oxidation of isoprene and monoterpenes by ozone and hydroxyl radicals, which produce semi-volatile organic compounds (SVOCs) and extremely low-volatility compounds (EVOCs). Subsequently, SVOCs and ELVOCs can condense on particles. Furthermore, we have introduced a new particle formation mechanism depending on the concentration of ELVOCs. For evaluation purposes, we have simulated the year 2010 with the old and new scheme; we see an increase in simulated production of SOA from 39.9 Tg yr−1 with the old scheme to 52.5 Tg yr−1 with the new scheme. For more detailed analysis, the particle mass and number concentrations and their influence on the simulated aerosol optical depth are compared to observations. Phenomenologically, the new particle formation scheme implemented here is able to reproduce the occurrence of observed particle formation events. However, the modelled concentrations of formed particles are clearly lower than in observations, as is the subsequent growth to larger sizes. Compared to the old scheme, the new scheme increases the number concentrations across the observation stations while still underestimating the observations. The organic aerosol mass concentrations in the US show a much better seasonal cycle and no clear overestimation of mass concentrations anymore. In Europe the mass concentrations are lowered, leading to a larger underestimation of observations. Aerosol optical depth (AOD) is generally slightly increased except in the northern high latitudes. This brings the simulated annual global mean AOD closer to the observational estimate. However, as the increase is rather uniform, biases tend to be reduced only in regions where the model underestimates the AOD. Furthermore, the correlations with satellite retrievals and ground-based sun-photometer observations of AOD are improved. Although the process-based approach to SOA formation causes a reduction in model performance in some areas, overall the new scheme improves the simulated aerosol fields.

This paper documents the global climate model EC-Earth3-AerChem, one of the members of the EC-Earth3 family of models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6). EC-Earth3-AerChem has interactive aerosols and atmospheric chemistry and contributes to the Aerosols and Chemistry Model Intercomparison Project (AerChemMIP). In this paper, we give an overview of the model, describe in detail how it differs from the other EC-Earth3 configurations, and outline the new features compared with the previously documented version of the model (EC-Earth 2.4). We explain how the model was tuned and spun up under preindustrial conditions and characterize the model's general performance on the basis of a selection of coupled simulations conducted for CMIP6. The net energy imbalance at the top of the atmosphere in the preindustrial control simulation is on average −0.09 W m−2 with a standard deviation due to interannual variability of 0.25 W m−2, showing no significant drift. The global surface air temperature in the simulation is on average 14.08 ∘C with an interannual standard deviation of 0.17 ∘C, exhibiting a small drift of 0.015 ± 0.005 ∘C per century. The model's effective equilibrium climate sensitivity is estimated at 3.9 ∘C, and its transient climate response is estimated at 2.1 ∘C. The CMIP6 historical simulation displays spurious interdecadal variability in Northern Hemisphere temperatures, resulting in a large spread across ensemble members and a tendency to underestimate observed annual surface temperature anomalies from the early 20th century onwards. The observed warming of the Southern Hemisphere is well reproduced by the model. Compared with the ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis version 5 (ERA5), the surface air temperature climatology for 1995–2014 has an average bias of −0.86 ± 0.05 ∘C with a standard deviation across ensemble members of 0.35 ∘C in the Northern Hemisphere and 1.29 ± 0.02 ∘C with a corresponding standard deviation of 0.05 ∘C in the Southern Hemisphere. The Southern Hemisphere warm bias is largely caused by errors in shortwave cloud radiative effects over the Southern Ocean, a deficiency of many climate models. Changes in the emissions of near-term climate forcers (NTCFs) have significant effects on the global climate from the second half of the 20th century onwards. For the SSP3-7.0 Shared Socioeconomic Pathway, the model gives a global warming at the end of the 21st century (2091–2100) of 4.9 ∘C above the preindustrial mean. A 0.5 ∘C stronger warming is obtained for the AerChemMIP scenario with reduced emissions of NTCFs. With concurrent reductions of future methane concentrations, the warming is projected to be reduced by 0.5 ∘C.

This paper presents an analysis of the mineral dust aerosol modelled by five Earth system models (ESMs) within the project entitled Coordinated Research in Earth Systems and Climate: Experiments, kNowledge, Dissemination and Outreach (CRESCENDO). We quantify the global dust cycle described by each model in terms of global emissions, together with dry and wet deposition, reporting large differences in the ratio of dry over wet deposition across the models not directly correlated with the range of particle sizes emitted. The multi-model mean dust emissions with five ESMs is 2836 Tg yr−1 but with a large uncertainty due mainly to the difference in the maximum dust particle size emitted. The multi-model mean of the subset of four ESMs without particle diameters larger than 10 µ m is 1664 (σ=651) Tg yr−1. Total dust emissions in the simulations with identical nudged winds from reanalysis give us better consistency between models; i.e. the multi-model mean global emissions with three ESMs are 1613 (σ=278) Tg yr−1, but 1834 (σ=666) Tg yr−1 without nudged winds and the same models. Significant discrepancies in the globally averaged dust mass extinction efficiency explain why even models with relatively similar global dust load budgets can display strong differences in dust optical depth. The comparison against observations has been done in terms of dust optical depths based on MODIS (Moderate Resolution Imaging Spectroradiometer) satellite products, showing global consistency in terms of preferential dust sources and transport across the Atlantic. The global localisation of source regions is consistent with MODIS, but we found regional and seasonal differences between models and observations when we quantified the cross-correlation of time series over dust-emitting regions. To faithfully compare local emissions between models we introduce a re-gridded normalisation method that can also be compared with satellite products derived from dust event frequencies. Dust total deposition is compared with an instrumental network to assess global and regional differences. We find that models agree with observations within a factor of 10 for data stations distant from dust sources, but the approximations of dust particle size distribution at emission contributed to a misrepresentation of the actual range of deposition values when instruments are close to dust-emitting regions. The observed dust surface concentrations also are reproduced to within a factor of 10. The comparison of total aerosol optical depth with AERONET (AErosol RObotic NETwork) stations where dust is dominant shows large differences between models, although with an increase in the inter-model consistency when the simulations are conducted with nudged winds. The increase in the model ensemble consistency also means better agreement with observations, which we have ascertained for dust total deposition, surface concentrations and optical depths (against both AERONET and MODIS retrievals). We introduce a method to ascertain the contributions per mode consistent with the multi-modal direct radiative effects, which we apply to study the direct radiative effects of a multi-modal representation of the dust particle size distribution that includes the largest particles.

This work documents and evaluates the tropospheric gas-phase chemical mechanism MOGUNTIA in the three-dimensional chemistry transport model TM5-MP. Compared to the modified CB05 (mCB05) chemical mechanism previously used in the model, MOGUNTIA includes a detailed representation of the light hydrocarbons (C1–C4) and isoprene, along with a simplified chemistry representation of terpenes and aromatics. Another feature implemented in TM5-MP for this work is the use of the Rosenbrock solver in the chemistry code, which can replace the classical Euler backward integration method of the model. Global budgets of ozone (O3), carbon monoxide (CO), hydroxyl radicals (OH), nitrogen oxides (NOx), and volatile organic compounds (VOCs) are analyzed, and their mixing ratios are compared with a series of surface, aircraft, and satellite observations for the year 2006. Both mechanisms appear to be able to satisfactorily represent observed mixing ratios of important trace gases, with the MOGUNTIA chemistry configuration yielding lower biases than mCB05 compared to measurements in most of the cases. However, the two chemical mechanisms fail to reproduce the observed mixing ratios of light VOCs, indicating insufficient primary emission source strengths, oxidation that is too fast, and/or a low bias in the secondary contribution to C2–C3 organics via VOC atmospheric oxidation. Relative computational memory and time requirements of the different model configurations are also compared and discussed. Overall, the MOGUNTIA scheme simulates a large suite of oxygenated VOCs that are observed in the atmosphere at significant levels. This significantly expands the possible applications of TM5-MP.

It is important to understand how future environmental policies will impact both climate change and air pollution. Although targeting near-term climate forcers (NTCFs), defined here as aerosols, tropospheric ozone, and precursor gases, should improve air quality, NTCF reductions will also impact climate. Prior assessments of the impact of NTCF mitigation on air quality and climate have been limited. This is related to the idealized nature of some prior studies, simplified treatment of aerosols and chemically reactive gases, as well as a lack of a sufficiently large number of models to quantify model diversity and robust responses. Here, we quantify the 2015–2055 climate and air quality effects of non-methane NTCFs using nine state-of-the-art chemistry–climate model simulations conducted for the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP). Simulations are driven by two future scenarios featuring similar increases in greenhouse gases (GHGs) but with “weak” (SSP3-7.0) versus “strong” (SSP3-7.0-lowNTCF) levels of air quality control measures. As SSP3-7.0 lacks climate policy and has the highest levels of NTCFs, our results (e.g., surface warming) represent an upper bound. Unsurprisingly, we find significant improvements in air quality under NTCF mitigation (strong versus weak air quality controls). Surface fine particulate matter (PM2.5) and ozone (O3) decrease by -2.2±0.32 µg m−3 and -4.6±0.88 ppb, respectively (changes quoted here are for the entire 2015–2055 time period; uncertainty represents the 95 % confidence interval), over global land surfaces, with larger reductions in some regions including south and southeast Asia. Non-methane NTCF mitigation, however, leads to additional climate change due to the removal of aerosol which causes a net warming effect, including global mean surface temperature and precipitation increases of 0.25±0.12 K and 0.03±0.012 mm d−1, respectively. Similarly, increases in extreme weather indices, including the hottest and wettest days, also occur. Regionally, the largest warming and wetting occurs over Asia, including central and north Asia (0.66±0.20 K and 0.03±0.02 mm d−1), south Asia (0.47±0.16 K and 0.17±0.09 mm d−1), and east Asia (0.46±0.20 K and 0.15±0.06 mm d−1). Relatively large warming and wetting of the Arctic also occur at 0.59±0.36 K and 0.04±0.02 mm d−1, respectively. Similar surface warming occurs in model simulations with aerosol-only mitigation, implying weak cooling due to ozone reductions. Our findings suggest that future policies that aggressively target non-methane NTCF reductions will improve air quality but will lead to additional surface warming, particularly in Asia and the Arctic. Policies that address other NTCFs including methane, as well as carbon dioxide emissions, must also be adopted to meet climate mitigation goals.

Biogenic volatile organic compounds (BVOCs) emitted from vegetation are oxidised in the atmosphere and can form aerosol particles either by contributing to new particle formation or by condensing onto existing aerosol particles. As the understanding of the importance of BVOCs for aerosol formation has increased over the years, these processes have made their way into Earth system models (ESMs). In this study, sensitivity experiments are run with three different ESMs (the Norwegian Earth System Model (NorESM), EC-Earth and ECHAM) to investigate how the direct and indirect aerosol radiative effects are affected by changes in the formation of secondary organic aerosol (SOA) from BVOCs. In the first two sensitivity model experiments, the yields of SOA precursors from oxidation of BVOCs are changed by ±50 %. For the third sensitivity test, the formed oxidation products do not participate in the formation of new particles but are only allowed to condense onto existing aerosols. In the last two sensitivity experiments, the emissions of BVOC compounds (isoprene and monoterpenes) are turned off, one at a time. The goal of the study is to investigate whether it is of importance to treat SOA formation processes correctly in the models rather than to evaluate the correctness of the current treatment in the models. The results show that the impact on the direct radiative effect (DRE) is linked to the changes in the SOA production in the models, where more SOA leads to a stronger DRE and vice versa. However, the magnitude by which the DRE changes (maximally 0.15 W m−2 globally averaged) in response to the SOA changes varies between the models, with EC-Earth displaying the largest changes. The results for the cloud radiative effects (CREs) are more complicated than for the DRE. The changes in CRE differ more among the ESMs, and for some sensitivity experiments they even have different signs. The most sensitive models are NorESM and EC-Earth, which have CRE changes of up to 0.82 W m−2. The varying responses in the different models are connected to where in the aerosol size distributions the changes in mass and number due to SOA formation occur, in combination with the aerosol number concentration levels in the models. We also find that interactive gas-phase chemistry as well as the new particle formation parameterisation has important implications for the DRE and CRE in some of the sensitivity experiments. The results from this study indicate that BVOC-SOA treatment in ESMs can have a substantial impact on the modelled climate but that the sensitivity varies greatly between the models. Since BVOC emissions have changed historically and will continue to change in the future, the spread in model results found in this study implies uncertainty into ESM estimates of aerosol forcing from land-use change and BVOC feedback strengths.

A total of 16 global chemistry transport models and general circulation models have participated in this study; 14 models have been evaluated with regard to their ability to reproduce the near-surface observed number concentration of aerosol particles and cloud condensation nuclei (CCN), as well as derived cloud droplet number concentration (CDNC). Model results for the period 2011–2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of −24 % and −35 % for particles with dry diameters >50 and >120 nm, as well as −36 % and −34 % for CCN at supersaturations of 0.2 % and 1.0 %, respectively. However, they seem to behave differently for particles activating at very low supersaturations ( %) than at higher ones. A total of 15 models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. A global mean reduction of a factor of about 2 is found in the model diversity for CCN at a supersaturation of 0.2 % (CCN0.2) compared to that for N3, maximizing over regions where new particle formation is important. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of the seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120 nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40 % during winter and 20 % in summer. In contrast to the large spread in simulated aerosol particle and CCN number concentrations, the CDNC derived from simulated CCN spectra is less diverse and in better agreement with CDNC estimates consistently derived from the observations (average NMB −13 % and −22 % for updraft velocities 0.3 and 0.6 m s−1, respectively). In addition, simulated CDNC is in slightly better agreement with observationally derived values at lower than at higher updraft velocities (index of agreement 0.64 vs. 0.65). The reduced spread of CDNC compared to that of CCN is attributed to the sublinear response of CDNC to aerosol particle number variations and the negative correlation between the sensitivities of CDNC to aerosol particle number concentration (∂Nd/∂Na) and to updraft velocity (∂Nd/∂w). Overall, we find that while CCN is controlled by both aerosol particle number and composition, CDNC is sensitive to CCN at low and moderate CCN concentrations and to the updraft velocity when CCN levels are high. Discrepancies are found in sensitivities ∂Nd/∂Na and ∂Nd/∂w; models may be predisposed to be too “aerosol sensitive” or “aerosol insensitive” in aerosol–cloud–climate interaction studies, even if they may capture average droplet numbers well. This is a subtle but profound finding that only the sensitivities can clearly reveal and may explain inter-model biases on the aerosol indirect effect.

In this paper, we present the implementation and evaluation of the aerosol microphysics module SALSA2.0 in the framework of the aerosol–chemistry–climate model ECHAM-HAMMOZ. It is an alternative microphysics module to the default modal microphysics scheme M7 in ECHAM-HAMMOZ. The SALSA2.0 implementation within ECHAM-HAMMOZ is evaluated against observations of aerosol optical properties, aerosol mass, and size distributions, comparing also to the skill of the M7 implementation. The largest differences between the implementation of SALSA2.0 and M7 are in the methods used for calculating microphysical processes, i.e., nucleation, condensation, coagulation, and hydration. These differences in the microphysics are reflected in the results so that the largest differences between SALSA2.0 and M7 are evident over regions where the aerosol size distribution is heavily modified by the microphysical processing of aerosol particles. Such regions are, for example, highly polluted regions and regions strongly affected by biomass burning. In addition, in a simulation of the 1991 Mt. Pinatubo eruption in which a stratospheric sulfate plume was formed, the global burden and the effective radii of the stratospheric aerosol are very different in SALSA2.0 and M7. While SALSA2.0 was able to reproduce the observed time evolution of the global burden of sulfate and the effective radii of stratospheric aerosol, M7 strongly overestimates the removal of coarse stratospheric particles and thus underestimates the effective radius of stratospheric aerosol. As the mode widths of M7 have been optimized for the troposphere and were not designed to represent stratospheric aerosol, the ability of M7 to simulate the volcano plume was improved by modifying the mode widths, decreasing the standard deviations of the accumulation and coarse modes from 1.59 and 2.0, respectively, to 1.2 similar to what was observed after the Mt. Pinatubo eruption. Overall, SALSA2.0 shows promise in improving the aerosol description of ECHAM-HAMMOZ and can be further improved by implementing methods for aerosol processes that are more suitable for the sectional method, e.g., size-dependent emissions for aerosol species and size-resolved wet deposition.

Aerosols have important impacts on air quality and climate, but the processes affecting their removal from the atmosphere are not fully understood and are poorly constrained by observations. This makes modelled aerosol lifetimes uncertain. In this study, we make use of an observational constraint on aerosol lifetimes provided by radionuclide measurements and investigate the causes of differences within a set of global models. During the Fukushima Dai-Ichi nuclear power plant accident of March 2011, the radioactive isotopes cesium-137 (137Cs) and xenon-133 (133Xe) were released in large quantities. Cesium attached to particles in the ambient air, approximately according to their available aerosol surface area. 137Cs size distribution measurements taken close to the power plant suggested that accumulation-mode (AM) sulfate aerosols were the main carriers of cesium. Hence, 137Cs can be used as a proxy tracer for the AM sulfate aerosol's fate in the atmosphere. In contrast, the noble gas 133Xe behaves almost like a passive transport tracer. Global surface measurements of the two radioactive isotopes taken over several months after the release allow the derivation of a lifetime of the carrier aerosol. We compare this to the lifetimes simulated by 19 different atmospheric transport models initialized with identical emissions of 137Cs that were assigned to an aerosol tracer with each model's default properties of AM sulfate, and 133Xe emissions that were assigned to a passive tracer. We investigate to what extent the modelled sulfate tracer can reproduce the measurements, especially with respect to the observed loss of aerosol mass with time. Modelled 137Cs and 133Xe concentrations sampled at the same location and times as station measurements allow a direct comparison between measured and modelled aerosol lifetime. The e-folding lifetime τe, calculated from station measurement data taken between 2 and 9 weeks after the start of the emissions, is 14.3 days (95 % confidence interval 13.1–15.7 days). The equivalent modelled τe lifetimes have a large spread, varying between 4.8 and 26.7 days with a model median of 9.4 ± 2.3 days, indicating too fast a removal in most models. Because sufficient measurement data were only available from about 2 weeks after the release, the estimated lifetimes apply to aerosols that have undergone long-range transport, i.e. not for freshly emitted aerosol. However, modelled instantaneous lifetimes show that the initial removal in the first 2 weeks was quicker (lifetimes between 1 and 5 days) due to the emissions occurring at low altitudes and co-location of the fresh plume with strong precipitation. Deviations between measured and modelled aerosol lifetimes are largest for the northernmost stations and at later time periods, suggesting that models do not transport enough of the aerosol towards the Arctic. The models underestimate passive tracer (133Xe) concentrations in the Arctic as well but to a smaller extent than for the aerosol (137Cs) tracer. This indicates that in addition to too fast an aerosol removal in the models, errors in simulated atmospheric transport towards the Arctic in most models also contribute to the underestimation of the Arctic aerosol concentrations.

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3–UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3–UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN > 3 nm), while the profiles of larger particles (e.g. CN > 100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.

Recently developed parameterizations for the sea spray aerosol source flux, encapsulating wave state, and its organic fraction were incorporated into the aerosol–climate model ECHAM-HAMMOZ to investigate the direct and indirect radiative effects of sea spray aerosol particles. Our simulated global sea salt emission of 805 Tg yr−1 (uncertainty range 378–1233 Tg yr−1) was much lower than typically found in previous studies. Modelled sea salt and sodium ion concentrations agreed relatively well with measurements in the smaller size ranges at Mace Head (annual normalized mean model bias −13% for particles with vacuum aerodynamic diameter Dva < 1 μm), Point Reyes (−29% for particles with aerodynamic diameter Da < 2.5 μm) and Amsterdam Island (−52% for particles with Da < 1 μm) but the larger sizes were overestimated (899% for particles with 2.5 μm < Da < 10 μm) at Amsterdam Island. This suggests that at least the high end of the previous estimates of sea spray mass emissions is unrealistic. On the other hand, the model clearly underestimated the observed concentrations of organic or total carbonaceous aerosol at Mace Head (−82%) and Amsterdam Island (−68%). The large overestimation (212%) of organic matter at Point Reyes was due to the contribution of continental sources. At the remote Amsterdam Island site, the organic concentration was underestimated especially in the biologically active months, suggesting a need to improve the parameterization of the organic sea spray fraction. Globally, the satellite-retrieved AOD over the oceans, using PARASOL data, was underestimated by the model (means over ocean 0.16 and 0.10, respectively); however, in the pristine region around Amsterdam Island the measured AOD fell well within the simulated uncertainty range. The simulated sea spray aerosol contribution to the indirect radiative effect was positive (0.3 W m−2), in contrast to previous studies. This positive effect was ascribed to the tendency of sea salt aerosol to suppress both the in-cloud supersaturation and the formation of cloud condensation nuclei from sulfate. These effects can be accounted for only in models with sufficiently detailed aerosol microphysics and physics-based parameterizations of cloud activation. However, due to a strong negative direct effect, the simulated effective radiative forcing (total radiative) effect was −0.2 W m−2. The simulated radiative effects of the primary marine organic emissions were small, with a direct effect of 0.03 W m−2 and an indirect effect of −0.07 W m−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 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.

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 12 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 data set 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.

We present the implementation and evaluation of a sectional aerosol microphysics module SALSA within the aerosol-climate model ECHAM5-HAM. This aerosol microphysics module has been designed to be flexible and computationally efficient so that it can be implemented in regional or global scale models. The computational efficiency has been achieved by minimising the number of variables needed to describe the size and composition distribution. The aerosol size distribution is described using 10 size classes with parallel sections which can have different chemical compositions. Thus in total, the module tracks 20 size sections which cover diameters ranging from 3 nm to 10 μm and are divided into three subranges, each with an optimised selection of processes and compounds. The implementation of SALSA into ECHAM5-HAM includes the main aerosol processes in the atmosphere: emissions, removal, radiative effects, liquid and gas phase sulphate chemistry, and the aerosol microphysics. The aerosol compounds treated in the module are sulphate, organic carbon, sea salt, black carbon, and mineral dust. In its default configuration, ECHAM5-HAM treats aerosol size distribution using the modal method. In this implementation, the aerosol processes were converted to be used in a sectional model framework. The ability of the module to describe the global aerosol properties was evaluated by comparing against (1) measured continental and marine size distributions, (2) observed variability of continental number concentrations, (3) measured sulphate, organic carbon, black carbon and sea-salt mass concentrations, (4) observations of aerosol optical depth (AOD) and other aerosol optical properties from satellites and AERONET network, (5) global aerosol budgets and concentrations from previous model studies, and (6) model results using M7, which is the default aerosol microphysics module in ECHAM5-HAM. The evaluation shows that the global aerosol properties can be reproduced reasonably well using a coarse resolution of 10 sections in size space. The simulated global aerosol budgets are within the range of previous studies. Surface concentrations of sulphate and carbonaceous species have an annual mean within a factor of two of the observations. The simulated sea-salt concentrations reproduce the observations within a factor of two, apart from the Southern Ocean over which the concentrations are within a factor of five. Regionally, AOD is in a relatively good agreement with the observations (within a factor of two). At mid-latitudes the observed AOD is captured well, while at high-latitudes as well as in some polluted and dust regions the modelled AOD is significantly lower than observed. Regarding most of the investigated aerosol properties, the SALSA and the modal aerosol module M7 perform comparably well against observations. However, SALSA reproduces the observed number concentrations and the size distribution of CCN sized particles much more accurately than M7, and is therefore a good choice for aerosol-cloud interaction studies in global models. Our study also shows that when activation type nucleation in the boundary layer is included, the observed concentration of particles under 50 nm in diameter are reproduced much better compared to when only binary nucleation in the free troposphere is assumed.

Biogenic aerosol formation is likely to contribute significantly to the global aerosol load. In recent years, new-particle formation has been observed in various ecosystems around the world but hardly any measurements have taken place in the terrestrial Southern Hemisphere. Here, we report the first results of atmospheric ion and charged particle concentrations as well as of new-particle formation in a Eucalypt forest in Tumbarumba, South-East Australia, from July 2005 to October 2006. The measurements were carried out with an Air Ion Spectrometer (AIS) with a size range from 0.34 to 40 nm. The Eucalypt forest was a very strong source of new aerosol particles. Daytime aerosol formation took place on 52% of days with acceptable data, which is 2–3 times as often as in the Nordic boreal zone. Average growth rates for negative/positive 1.5–3 nm particles during these formation events were 2.89/2.68 nmh−1, respectively; for 3-7 nm particles 4.26/4.03, and for 7–20 nm particles 8.90/7.58 nmh−1, respectively. The growth rates for large ions were highest when the air was coming from the native forest which suggests that the Eucalypts were a strong source of condensable vapours. Average concentrations of cluster ions (0.34–1.8 nm) were 2400/1700 cm−3 for negative/positive ions, very high compared to most other measurements around the world. One reason behind these high concentrations could be the strong radon efflux from the soils around the Tumbarumba field site. Furthermore, comparison between night-time and daytime concentrations supported the view that cluster ions are produced close to the surface within the boundary layer also at night but that large ions are mostly produced in daytime. Finally, a previously unreported phenomenon, nocturnal aerosol formation, appeared in 32% of the analysed nights but was clustered almost entirely within six months from summer to autumn in 2006. From January to May, nocturnal formation was 2.5 times as frequent as daytime formation. Therefore, it appears that in summer and autumn, nocturnal production was the major mechanism for aerosol formation in Tumbarumba.

We have measured the size distributions of air ions (0.42–7.5 nm in diameter) with the Balanced Scanning Mobility Analyzer in boreal forest, in Southern Finland since spring 2003. The size range covers the size range of cluster ions (approximately 0.42–1.6 nm) and naturally charged nanometre aerosol particles (1.6–7.5 nm) or intermediate air ions. Based on the measurements from April 2003 to March 2006 we studied the characteristics of charged aerosol particle formation by classifying each day either as a particle formation event, undefined or non-event day. The principal of the classification, as well as the statistical description of the charged aerosol particle formation events are given. We found in total 270 (26% of the analysed days) and 226 (22% of the analysed days) particle formation days for negative and positive intermediate ions, respectively. For negatively charged particles we classified 411 (40% of the analysed days) undefined and 348 (34% of the analysed days) non-event days whereas for positively charged particles 343 (33% of the analysed days) undefined and 460 (45% of the analysed days) non-event days. The results were compared with the ordinary classification based on the Differential Mobility Particle Sizer (DMPS) measurements carried out at the same place. The above-presented values differed slightly from that found from the DMPS data, with a lower particle diameter of 3 nm. In addition, we have found the rain-induced intermediate ion bursts frequently. The rain effect was detected on 163 days by means of negative ions and on 105 days by positive ones. Another interesting phenomenon among the charged aerosol particles was the appearance and existence of intermediate ions during the snowfall. We observed this phenomenon 24 times with negatively charged particles and 21 times with positively charged ones during winter months (October–April). These intermediate air ions were seen during the snowfall and may be caused by ice crystals, although the origin of these intermediate ions is unclear at the moment.