Search

Aerosol radiative forcing uncertainty affects estimates of climate sensitivity and limits model skill at making climate projections. Efforts to improve the representations of physical processes in climate models, including extensive comparisons with observations, have not significantly constrained the range of possible aerosol forcing values. A far stronger constraint, in particular for the lower (most-negative) bound, can be achieved using global mean energy-balance arguments based on observed changes in historical temperature. Here, we show that structural deficiencies in a climate model, revealed as inconsistencies among observationally constrained cloud properties in the model, limit the effectiveness of observational constraint of the uncertain physical processes. We sample uncertainty in 37 model parameters related to aerosols, clouds and radiation in a perturbed parameter ensemble of the UK Earth System Model and evaluate 1 million model variants (different parameter settings from Gaussian Process emulators) against satellite-derived observations over several cloudy regions. We show that it is possible to reduce the parametric uncertainty in global mean aerosol forcing by more than 50 %, constraining it to a range in close agreement with energy-balance constraints (around −1.3 to −0.1 W m−2). However, our analysis of a very large set of model variants exposes model internal inconsistencies that would not be apparent in a small set of model simulations. Incorporating observations associated with these inconsistencies weakens the forcing constraint because they require a wider range of parameter values to accommodate conflicting information. Our estimated aerosol forcing range is the maximum feasible constraint using our structurally imperfect model and the chosen observations. Structural model developments targeted at the identified inconsistencies would enable a larger set of observations to be used for constraint, which would then narrow the uncertainty further. Such an approach provides a rigorous pathway to improved model realism and reduced uncertainty that has so far not been achieved through the normal model development approach.

In this study we evaluate simulated surface SO2 and sulfate (SO42-) concentrations from the United Kingdom Earth System Model (UKESM1) against observations from ground-based measurement networks in the USA and Europe for the period 1987–2014. We find that UKESM1 captures the historical trend for decreasing concentrations of atmospheric SO2 and SO42- in both Europe and the USA over the period 1987–2014. However, in the polluted regions of the eastern USA and Europe, UKESM1 over-predicts surface SO2 concentrations by a factor of 3 while under-predicting surface SO42- concentrations by 25 %–35 %. In the cleaner western USA, the model over-predicts both surface SO2 and SO42- concentrations by factors of 12 and 1.5 respectively. We find that UKESM1’s bias in surface SO2 and SO42- concentrations is variable according to region and season. We also evaluate UKESM1 against total column SO2 from the Ozone Monitoring Instrument (OMI) using an updated data product. This comparison provides information about the model's global performance, finding that UKESM1 over-predicts total column SO2 over much of the globe, including the large source regions of India, China, the USA, and Europe as well as over outflow regions. Finally, we assess the impact of a more realistic treatment of the model's SO2 dry deposition parameterization. This change increases SO2 dry deposition to the land and ocean surfaces, thus reducing the atmospheric loading of SO2 and SO42-. In comparison with the ground-based and satellite observations, we find that the modified parameterization reduces the model's over-prediction of surface SO2 concentrations and total column SO2. Relative to the ground-based observations, the simulated surface SO42- concentrations are also reduced, while the simulated SO2 dry deposition fluxes increase.

Despite a large number of studies, 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 different sizes in total of for 21 measurement sites in Europe and Arctic. For 13 of those, 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 (ESM). 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, were typically stronger in 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, for others, 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 showing clear differences in concentrations between remote and urban sites. To conclude, although anthropogenic mass emissions are harmonized in models, trends in different sized particles vary among the model 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.

Quantifying forcings from anthropogenic perturbations to the Earth System (ES) is important for understanding changes in climate since the pre-industrial period. In this paper, we quantify and analyse a wide range of present-day (PD) anthropogenic climate forcings with the UK's Earth System Model (ESM), UKESM1, following the protocols defined by the Radiative Forcing Model Intercomparison Project (RFMIP) and the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP). In particular, by quantifying effective radiative forcings (ERFs) that include rapid adjustments within a full ESM, it enables the role of various climate-chemistry-aerosol-cloud feedbacks to be quantified. Global mean ERFs are 1.83, 0.13, −0.33, and 0.93 W m−2 at the PD (Year 2014) relative to the pre-industrial (PI; Year 1850) for carbon dioxide, nitrous oxide, ozone-depleting substances, and methane, respectively. The PD total greenhouse gas ERF is 2.89 W m−2, larger than the sum of the individual GHG ERFs. UKESM1 has an aerosol forcing of −1.13 W m−2. A relatively strong negative forcing from aerosol-cloud interactions and a small negative instantaneous forcing from aerosol-radiation interactions are partially offset by a substantial forcing from black carbon absorption. Internal mixing and chemical interactions mean that neither the forcing from aerosol-radiation interactions nor aerosol-cloud interactions are linear, making the total aerosol ERF less than the sum of the individual speciated aerosol ERFs. Tropospheric ozone precursors, in addition to exerting a positive forcing due to ozone, lead to oxidant changes which in turn cause an indirect aerosol ERF, altering the sign of the net ERF from nitrogen oxide emissions. Together, aerosol and tropospheric ozone precursors (near-term climate forcers, NTCFs) exert a global mean ERF of −1.12 W m−2, mainly due to changes in the cloud radiative effect. There is also a negative PD ERF from land use (−0.32 W m−2). It is outside the range of previous estimates, and is most likely due to too strong an albedo response. In combination, the net anthropogenic ERF is potentially biased low (1.61 W m−2) relative to other estimates, due to the inclusion of non-linear feedbacks and ES interactions. By including feedbacks between greenhouse gases, stratospheric and tropospheric ozone, aerosols, and clouds, some of which act non-linearly, this work demonstrates the importance of ES interactions when quantifying climate forcing. It also suggests that rapid adjustments need to include chemical as well as physical adjustments to fully account for complex ES interactions.