Abstract. The Copernicus Atmosphere Monitoring Service (CAMS) reanalysis is the latest global reanalysis dataset of atmospheric composition produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), consisting of three-dimensional time-consistent atmospheric composition fields, including aerosols and chemical species. The dataset currently covers the period 2003–2016 and will be extended in the future by adding 1 year each year. A reanalysis for greenhouse gases is being produced separately. The CAMS reanalysis builds on the experience gained during the production of the earlier Monitoring Atmospheric Composition and Climate (MACC) reanalysis and CAMS interim reanalysis. Satellite retrievals of total column CO; tropospheric column NO2; aerosol optical depth (AOD); and total column, partial column and profile ozone retrievals were assimilated for the CAMS reanalysis with ECMWF's Integrated Forecasting System. The new reanalysis has an increased horizontal resolution of about 80 km and provides more chemical species at a better temporal resolution (3-hourly analysis fields, 3-hourly forecast fields and hourly surface forecast fields) than the previously produced CAMS interim reanalysis. The CAMS reanalysis has smaller biases compared with most of the independent ozone, carbon monoxide, nitrogen dioxide and aerosol optical depth observations used for validation in this paper than the previous two reanalyses and is much improved and more consistent in time, especially compared to the MACC reanalysis. The CAMS reanalysis is a dataset that can be used to compute climatologies, study trends, evaluate models, benchmark other reanalyses or serve as boundary conditions for regional models for past periods.

A large number of processes are involved in the chain from emissions of aerosol precursor gases and primary particles to impacts on cloud radiative forcing. Those processes are manifest in a number of relationships that can be expressed as factors dlnX/dlnY driving aerosol effects on cloud radiative forcing. These factors include the relationships between cloud condensation nuclei (CCN) concentration and emissions, droplet number and CCN concentration, cloud fraction and droplet number, cloud optical depth and droplet number, and cloud radiative forcing and cloud optical depth. The relationship between cloud optical depth and droplet number can be further decomposed into the sum of two terms involving the relationship of droplet effective radius and cloud liquid water path with droplet number. These relationships can be constrained using observations of recent spatial and temporal variability of these quantities. However, we are most interested in the radiative forcing since the preindustrial era. Because few relevant measurements are available from that era, relationships from recent variability have been assumed to be applicable to the preindustrial to present-day change. Our analysis of Aerosol Comparisons between Observations and Models (AeroCom) model simulations suggests that estimates of relationships from recent variability are poor constraints on relationships from anthropogenic change for some terms, with even the sign of some relationships differing in many regions. Proxies connecting recent spatial/temporal variability to anthropogenic change, or sustained measurements in regions where emissions have changed, are needed to constrain estimates of anthropogenic aerosol impacts on cloud radiative forcing.

Abstract. 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 (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) 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 AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc 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 AODc. 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 SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry 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 %.

Abstract. This article describes the IFS-AER aerosol module used operationally in the Integrated Forecasting System (IFS) cycle 45R1, operated by the European Centre for Medium-Range Weather Forecasts (ECMWF) in the framework of the Copernicus Atmospheric Monitoring Services (CAMS). We describe the different parameterizations for aerosol sources, sinks, and its chemical production in IFS-AER, as well as how the aerosols are integrated in the larger atmospheric composition forecasting system. The focus is on the entire 45R1 code base, including some components that are not used operationally, in which case this will be clearly specified. This paper is an update to the Morcrette et al. (2009) article that described aerosol forecasts at the ECMWF using cycle 32R2 of the IFS. Between cycles 32R2 and 45R1, a number of source and sink processes have been reviewed and/or added, notably increasing the complexity of IFS-AER. A greater integration with the tropospheric chemistry scheme of the IFS has been achieved for the sulfur cycle and for nitrate production. Two new species, nitrate and ammonium, have also been included in the forecasting system. Global budgets and aerosol optical depth (AOD) fields are shown, as is an evaluation of the simulated particulate matter (PM) and AOD against observations, showing an increase in skill from cycle 40R2, used in the CAMS interim ReAnalysis (CAMSiRA), to cycle 45R1.

Abstract. The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration.This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs.This study focuses on the effect of the variance (σ 2 w ) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σ w (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σ w depends on turbulent kinetic energy are tested. Results from the mid-range fixed σ w and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations.The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m −2 with σ w = 0.1 m s −1 , −2.1 W m −2 with σ w derived from TKE, −2.25 W m −2 with σ w = 0.4 m s −1 , and −2.3 W m −2 with σ w = 0.7 m s −1 . The breadth of this range is 0.4 W m −2 , which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σ w would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects.Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Abstract. We document and evaluate the aerosol schemes as implemented in the physical and Earth system models, the Global Coupled 3.1 configuration of the Hadley Centre Global Environment Model version 3 (HadGEM3-GC3.1) and the United Kingdom Earth System Model (UKESM1), which are contributing to the sixth Coupled Model Intercomparison Project (CMIP6). The simulation of aerosols in the present-day period of the historical ensemble of these models is evaluated against a range of observations. Updates to the aerosol microphysics scheme are documented as well as differences in the aerosol representation between the physical and Earth system configurations. The additional Earth system interactions included in UKESM1 lead to differences in the emissions of natural aerosol sources such as dimethyl sulfide, mineral dust and organic aerosol and subsequent evolution of these species in the model. UKESM1 also includes a stratospheric–tropospheric chemistry scheme which is fully coupled to the aerosol scheme, while GC3.1 employs a simplified aerosol chemistry mechanism driven by prescribed monthly climatologies of the relevant oxidants. Overall, the simulated speciated aerosol mass concentrations compare reasonably well with observations. Both models capture the negative trend in sulfate aerosol concentrations over Europe and the eastern United States of America (US) although the models tend to underestimate sulfate concentrations in both regions. Interactive emissions of biogenic volatile organic compounds in UKESM1 lead to an improved agreement of organic aerosol over the US. Simulated dust burdens are similar in both models despite a 2-fold difference in dust emissions. Aerosol optical depth is biased low in dust source and outflow regions but performs well in other regions compared to a number of satellite and ground-based retrievals of aerosol optical depth. Simulated aerosol number concentrations are generally within a factor of 2 of the observations, with both models tending to overestimate number concentrations over remote ocean regions, apart from at high latitudes, and underestimate over Northern Hemisphere continents. Finally, a new primary marine organic aerosol source is implemented in UKESM1 for the first time. The impact of this new aerosol source is evaluated. Over the pristine Southern Ocean, it is found to improve the seasonal cycle of organic aerosol mass and cloud droplet number concentrations relative to GC3.1 although underestimations in cloud droplet number concentrations remain. This paper provides a useful characterisation of the aerosol climatology in both models and will facilitate understanding in the numerous aerosol–climate interaction studies that will be conducted as part of CMIP6 and beyond.

Evaluation of the aerosol schemes in current climate models is dependent upon the available observational data. In-situ observations from flight campaigns can provide valuable data about the vertical distribution of aerosol that is difficult to obtain from satellite or ground-based platforms, although they are localised in space and time. Using single-particle soot-photometer (SP2) measurements from the HIAPER Pole-to-Pole Observations (HIPPO) campaign, which consists of many vertical profiles over a large region of the Pacific, we evaluate the meridional and vertical distribution of black carbon (BC) aerosol simulated by the HadGEM3–UKCA and ECHAM5–HAM2 models. Both models show a similar pattern of overestimating the BC column burden compared to that derived from the observations, in many areas by an order of magnitude. However, by sampling the simulated BC mass mixing ratio along the flight track and comparing to the observations, we show that this discrepancy has a rather different vertical structure in the two models: in HadGEM3–UKCA the discrepancy is dominated by excess aerosol in the tropical upper troposphere, while in ECHAM5–HAM2 areas of discrepancy are spread across many different latitudes and altitudes.

Using this methodology, we conduct sensitivity tests on two specific elements of the models: biomass-burning emissions and scavenging by convective precipitation. We show that, by coupling the convective scavenging more tightly with convective transport, both the column burden and vertical distribution of BC in HadGEM3–UKCA are much improved with respect to the observations, with a substantial and statistically significant increase in correlation – this demonstrates the importance of a realistic representation of this process. In contrast, updating from GFED2 to GFED3.1 biomass-burning emissions makes a more modest improvement in both models, which is not statistically significant. By comparing our results with a more traditional approach using regional- and monthly-mean vertical profile curves, we show that the point-by-point analysis allows the model improvements to be demonstrated more clearly.

We also demonstrate the important role that nudged simulations (where the large-scale model dynamics are continuously relaxed towards a reanalysis) can play in this type of evaluation, allowing statistically significant differences between configurations of the aerosol scheme to be seen where the differences between the corresponding free-running simulations would not be significant

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