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 −1 . Plain Language SummaryWe describe the development and behavior of UKESM1, a novel climate model that includes improved representations of processes in the atmosphere, ocean, and on land. These processes are inter-related: For example, dust is produced on the land and blown up into the atmosphere where it affects the amount of sunlight falling on Earth. Dust can also be dissolved in the ocean, where it affects marine life. This in turn changes both the amount of carbon dioxide absorbed by the ocean and the material emitted from the surface into the atmosphere, which has an affect on the formation of clouds. UKESM1 includes many processes and interactions such as these, giving it a high level of complexity. Ensuring realistic process behavior is a major challenge in the development of our model, and we have carefully tested this. UKESM1 performs well, correctly exhibiting stable results from a continuous pre-industrial simulation (used to provide a reference for future experiments) and showing good agreement

Abstract. We describe Global Atmosphere 7.0 and Global Land 7.0 (GA7.0/GL7.0), the latest science configurations of the Met Office Unified Model (UM) and the Joint UK Land Environment Simulator (JULES) land surface model developed for use across weather and climate timescales. GA7.0 and GL7.0 include incremental developments and targeted improvements that, between them, address four critical errors identified in previous configurations: excessive precipitation biases over India, warm and moist biases in the tropical tropopause layer (TTL), a source of energy non-conservation in the advection scheme and excessive surface radiation biases over the Southern Ocean. They also include two new parametrisations, namely the UK Chemistry and Aerosol (UKCA) GLOMAP-mode (Global Model of Aerosol Processes) aerosol scheme and the JULES multi-layer snow scheme, which improve the fidelity of the simulation and were required for inclusion in the Global Atmosphere/Global Land configurations ahead of the 6th Coupled Model Intercomparison Project (CMIP6). In addition, we describe the GA7.1 branch configuration, which reduces an overly negative anthropogenic aerosol effective radiative forcing (ERF) in GA7.0 whilst maintaining the quality of simulations of the present-day climate. GA7.1/GL7.0 will form the physical atmosphere/land component in the HadGEM3–GC3.1 and UKESM1 climate model submissions to the CMIP6.

Aerosol processes and, in particular, aerosol‐cloud interactions cut across the traditional physical‐Earth system boundary of coupled Earth system models and remain one of the key uncertainties in estimating anthropogenic radiative forcing of climate. Here we calculate the historical aerosol effective radiative forcing (ERF) in the HadGEM3‐GA7 climate model in order to assess the suitability of this model for inclusion in the UK Earth system model, UKESM1. The aerosol ERF, calculated for the year 2000 relative to 1850, is large and negative in the standard GA7 model leading to an unrealistic negative total anthropogenic forcing over the twentieth century. We show how underlying assumptions and missing processes in both the physical model and aerosol parameterizations lead to this large aerosol ERF. A number of model improvements are investigated to assess their impact on the aerosol ERF. These include an improved representation of cloud droplet spectral dispersion, updates to the aerosol activation scheme, and black carbon optical properties. One of the largest contributors to the aerosol forcing uncertainty is insufficient knowledge of the preindustrial aerosol climate. We evaluate the contribution of uncertainties in the natural marine emissions of dimethyl sulfide and organic aerosol to the ERF. The combination of model improvements derived from these studies weakens the aerosol ERF by up to 50% of the original value and leads to a total anthropogenic historical forcing more in line with assessed values.

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

A new climate model, HadGEM3 N96ORCA1, is presented that is part of the GC3.1 configuration of HadGEM3. N96ORCA1 has a horizontal resolution of ~135 km in the atmosphere and 1° in the ocean and requires an order of magnitude less computing power than its medium‐resolution counterpart, N216ORCA025, while retaining a high degree of performance traceability. Scientific performance is compared to both observations and the N216ORCA025 model. N96ORCA1 reproduces observed climate mean and variability almost as well as N216ORCA025. Patterns of biases are similar across the two models. In the northwest Atlantic, N96ORCA1 shows a cold surface bias of up to 6 K, typical of ocean models of this resolution. The strength of the Atlantic meridional overturning circulation (16 to 17 Sv) matches observations. In the Southern Ocean, a warm surface bias (up to 2 K) is smaller than in N216ORCA025 and linked to improved ocean circulation. Model El Niño/Southern Oscillation and Atlantic Multidecadal Variability are close to observations. Both the cold bias in the Northern Hemisphere (N96ORCA1) and the warm bias in the Southern Hemisphere (N216ORCA025) develop in the first few decades of the simulations. As in many comparable climate models, simulated interhemispheric gradients of top‐of‐atmosphere radiation are larger than observations suggest, with contributions from both hemispheres. HadGEM3 GC3.1 N96ORCA1 constitutes the physical core of the UK Earth System Model (UKESM1) and will be used extensively in the Coupled Model Intercomparison Project 6 (CMIP6), both as part of the UK Earth System Model and as a stand‐alone coupled climate model.

[1] The relationship between the AOD and wind speed (U), obtained at the Mace Head atmospheric research station, is determined for clean marine (Northeast Atlantic) conditions under steady-state wind speed scenarios. AOD measurements were made under much higher wind speed conditions (up to 18 m s À1) than hitherto measured. The AOD is found to be approximately dependent on U 2 , with a high correlation (r 2 = 0.97) determined at all measured wavelengths of 862 nm, 500 nm, 412 nm and 368 nm. A corresponding anti-correlation was observed between the Å ngström exponent and wind speed (r 2 = 0.93). The results illustrate that under moderate-to-high wind conditions, sea-spray AOD can rival or even exceed that associated with pollution plumes over oceanic regions and can potentially contribute to radiative, dynamical and biogeochemical feedback processes associated with future climate change.

Abstract. We describe Global Atmosphere 7.0 and Global Land 7.0 (GA7.0/GL7.0): the latest science configurations of the Met Office Unified Model and JULES land surface model developed for use across weather and climate timescales. GA7.0 and GL7.0 include incremental developments and targeted improvements that between them address four critical errors identified in previous configurations: excessive precipitation biases over India, warm/moist biases in the tropical tropopause layer, a source of energy non-conservation in the advection scheme and excessive surface radiation biases over the Southern Ocean. They also 5 include two new parametrisations, namely the UKCA Glomap-mode aerosol scheme and the JULES multi-layer snow scheme, which improve the fidelity of the simulation and were required for the coupled climate models and Earth system models that will use these configurations in submissions to CMIP6.In addition, we describe the GA7.1 branch configuration, which reduces an overly negative anthropogenic aerosol effective radiative forcing in GA7.0, whilst maintaining the quality of simulations of the present-day climate. GA7.1/GL7.0 will form 10 the physical atmosphere/land component in the HadGEM3-GC3.1 and UKESM1 climate models.

We analyze the atmospheric processes that explain the large changes in radiative feedbacks between the two latest climate configurations of the Hadley Centre Global Environmental model. We use a large set of atmosphere‐only climate change simulations (amip and amip‐p4K) to separate the contributions to the differences in feedback parameter from all the atmospheric model developments between the two latest model configurations. We show that the differences are mostly driven by changes in the shortwave cloud radiative feedback in the midlatitudes, mainly over the Southern Ocean. Two new schemes explain most of the differences: the introduction of a new aerosol scheme and the development of a new mixed‐phase cloud scheme. Both schemes reduce the strength of the preexisting shortwave negative cloud feedback in the midlatitudes. The new aerosol scheme dampens a strong aerosol‐cloud interaction, and it also suppresses a negative clear‐sky shortwave feedback. The mixed‐phase scheme increases the amount of cloud liquid water path (LWP) in the present day and reduces the increase in LWP with warming. Both changes contribute to reducing the negative radiative feedback of the increase of LWP in the warmer climate. The mixed‐phase scheme also enhances a strong, preexisting, positive cloud fraction feedback. We assess the realism of the changes by comparing present‐day simulations against observations and discuss avenues that could help constrain the relevant processes.