We document the development of the first version of the U.K. Earth System Model UKESM1. The model represents a major advance on its predecessor HadGEM2‐ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well‐resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric‐stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two‐moment, five‐species, modal aerosol; and ocean biogeochemistry with two‐way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre‐industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger‐than‐observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC⁻¹.

Here we present a description of the UKCA StratTrop chemical mechanism, which is used in the UKESM1 Earth system model for CMIP6. The StratTrop chemical mechanism is a merger of previously well-evaluated tropospheric and stratospheric mechanisms, and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained for tropospheric ozone and its budget terms from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address.

We present an overview of state-of-the-art chemistry–climate and chemistry transport models that are used within phase 1 of the Chemistry–Climate Model Initiative (CCMI-1). The CCMI aims to conduct a detailed evaluation of participating models using process-oriented diagnostics derived from observations in order to gain confidence in the models' projections of the stratospheric ozone layer, tropospheric composition, air quality, where applicable global climate change, and the interactions between them. Interpretation of these diagnostics requires detailed knowledge of the radiative, chemical, dynamical, and physical processes incorporated in the models. Also an understanding of the degree to which CCMI-1 recommendations for simulations have been followed is necessary to understand model responses to anthropogenic and natural forcing and also to explain inter-model differences. This becomes even more important given the ongoing development and the ever-growing complexity of these models. This paper also provides an overview of the available CCMI-1 simulations with the aim of informing CCMI data users.

We describe the scientific and technical implementation of two models for a core set of experiments contributing to the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The models used are the physical atmosphere‐land‐ocean‐sea ice model HadGEM3‐GC3.1 and the Earth system model UKESM1 which adds a carbon‐nitrogen cycle and atmospheric chemistry to HadGEM3‐GC3.1. The model results are constrained by the external boundary conditions (forcing data) and initial conditions. We outline the scientific rationale and assumptions made in specifying these. Notable details of the implementation include an ozone redistribution scheme for prescribed ozone simulations (HadGEM3‐GC3.1) to avoid inconsistencies with the model's thermal tropopause, and land use change in dynamic vegetation simulations (UKESM1) whose influence will be subject to potential biases in the simulation of background natural vegetation. We discuss the implications of these decisions for interpretation of the simulation results. These simulations are expensive in terms of human and CPU resources and will underpin many further experiments; we describe some of the technical steps taken to ensure their scientific robustness and reproducibility.

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.

The evolution of tropospheric ozone from 1850 to 2100 has been studied using data from Phase 6 of the Coupled Model Intercomparison Project (CMIP6). We evaluate long-term changes using coupled atmosphere–ocean chemistry–climate models, focusing on the CMIP Historical and ScenarioMIP ssp370 experiments, for which detailed tropospheric-ozone diagnostics were archived. The model ensemble has been evaluated against a suite of surface, sonde and satellite observations of the past several decades and found to reproduce well the salient spatial, seasonal and decadal variability and trends. The multi-model mean tropospheric-ozone burden increases from 247 ± 36 Tg in 1850 to a mean value of 356 ± 31 Tg for the period 2005–2014, an increase of 44 %. Modelled present-day values agree well with previous determinations (ACCENT: 336 ± 27 Tg; Atmospheric Chemistry and Climate Model Intercomparison Project, ACCMIP: 337 ± 23 Tg; Tropospheric Ozone Assessment Report, TOAR: 340 ± 34 Tg). In the ssp370 experiments, the ozone burden increases to 416 ± 35 Tg by 2100. The ozone budget has been examined over the same period using lumped ozone production ($P_{{\mathrm{O}}_{\mathrm{3}}}$ ) and loss ($L_{{\mathrm{O}}_{\mathrm{3}}}$ ) diagnostics. Both ozone production and chemical loss terms increase steadily over the period 1850 to 2100, with net chemical production ($P_{{\mathrm{O}}_{\mathrm{3}}}$ -$L_{{\mathrm{O}}_{\mathrm{3}}}$ ) reaching a maximum around the year 2000. The residual term, which contains contributions from stratosphere–troposphere transport reaches a minimum around the same time before recovering in the 21st century, while dry deposition increases steadily over the period 1850–2100. Differences between the model residual terms are explained in terms of variation in tropopause height and stratospheric ozone burden.