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.

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.

This paper quantifies the pre-industrial (1850) to present-day (2014) effective radiative forcing (ERF) of anthropogenic emissions of NO_X, volatile organic compounds (VOCs; including CO), SO₂, NH₃, black carbon, organic carbon, and concentrations of methane, N₂O and ozone-depleting halocarbons, using CMIP6 models. Concentration and emission changes of reactive species can cause multiple changes in the composition of radiatively active species: tropospheric ozone, stratospheric ozone, stratospheric water vapour, secondary inorganic and organic aerosol, and methane. Where possible we break down the ERFs from each emitted species into the contributions from the composition changes. The ERFs are calculated for each of the models that participated in the AerChemMIP experiments as part of the CMIP6 project, where the relevant model output was available. The 1850 to 2014 multi-model mean ERFs (± standard deviations) are −1.03 ± 0.37 W m⁻² for SO₂ emissions, −0.25 ± 0.09 W m⁻² for organic carbon (OC), 0.15 ± 0.17 W m⁻² for black carbon (BC) and −0.07 ± 0.01 W m⁻² for NH₃. For the combined aerosols (in the piClim-aer experiment) it is −1.01 ± 0.25 W m⁻². The multi-model means for the reactive well-mixed greenhouse gases (including any effects on ozone and aerosol chemistry) are 0.67 ± 0.17 W m⁻² for methane (CH₄), 0.26 ± 0.07 W m⁻² for nitrous oxide (N₂O) and 0.12 ± 0.2 W m⁻² for ozone-depleting halocarbons (HC). Emissions of the ozone precursors nitrogen oxides (NO_x), volatile organic compounds and both together (O₃) lead to ERFs of 0.14 ± 0.13, 0.09 ± 0.14 and 0.20 ± 0.07 W m⁻² respectively. The differences in ERFs calculated for the different models reflect differences in the complexity of their aerosol and chemistry schemes, especially in the case of methane where tropospheric chemistry captures increased forcing from ozone production.

Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼ 300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer–Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions.

Hydrogen is expected to play a key role in the global energy transition to net zero emissions in many scenarios. However, fugitive emissions of hydrogen into the atmosphere during its production, storage, distribution and use could reduce the climate benefit and also have implications for air quality. Here, we explore the atmospheric composition and climate impacts of increases in atmospheric hydrogen abundance using the UK Earth System Model (UKESM1) chemistry–climate model. Increases in hydrogen result in increases in methane, tropospheric ozone and stratospheric water vapour, resulting in a positive radiative forcing. However, some of the impacts of hydrogen leakage are partially offset by potential reductions in emissions of methane, carbon monoxide, nitrogen oxides and volatile organic compounds from the consumption of fossil fuels. We derive a refined methodology for determining indirect global warming potentials (GWPs) from parameters derived from steady-state simulations, which is applicable to both shorter-lived species and those with intermediate and longer lifetimes, such as hydrogen. Using this methodology, we determine a 100-year global warming potential for hydrogen of 12 ± 6. Based on this GWP and hydrogen leakage rates of 1 % and 10 %, we find that hydrogen leakage offsets approximately 0.4 % and 4 % respectively of total equivalent CO₂ emission reductions in our global hydrogen economy scenario. To maximise the benefit of hydrogen as an energy source, emissions associated with hydrogen leakage and emissions of the ozone precursor gases need to be minimised.