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⁻¹.

This paper reviews recent literature concerning a wide range of processes through which climate change could potentially impact global-scale agricultural productivity, and presents projections of changes in relevant meteorological, hydrological and plant physiological quantities from a climate model ensemble to illustrate key areas of uncertainty. Few global-scale assessments have been carried out, and these are limited in their ability to capture the uncertainty in climate projections, and omit potentially important aspects such as extreme events and changes in pests and diseases. There is a lack of clarity on how climate change impacts on drought are best quantified from an agricultural perspective, with different metrics giving very different impressions of future risk. The dependence of some regional agriculture on remote rainfall, snowmelt and glaciers adds to the complexity. Indirect impacts via sea-level rise, storms and diseases have not been quantified. Perhaps most seriously, there is high uncertainty in the extent to which the direct effects of CO2 rise on plant physiology will interact with climate change in affecting productivity. At present, the aggregate impacts of climate change on global-scale agricultural productivity cannot be reliably quantified.

A decline in the carbon fertilization effect: One source of uncertainty in climate science is how the carbon fertilization effect (CFE) will contribute to mitigation of anthropogenic climate change. Wang et al. explored the temporal dynamics of CFE on vegetation photosynthesis at the global scale. There has been a decline over recent decades in the contribution of CFE to vegetation photosynthesis, perhaps owing to the limiting effects of plant nutrients such as nitrogen and phosphorus. This declining trend has not been adequately accounted for in carbon cycle models. CFE thus has limitations for long-term mitigation of climate change, and future warming might currently be underestimated. Science , this issue p. 1295

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.

Understanding future changes in the terrestrial carbon cycle is important for reliable projections of climate change and impacts on ecosystems. It is well known that nitrogen (N) could limit plants' response to increased atmospheric carbon dioxide and it is therefore important to include a representation of the N cycle in Earth system models. Here we present the implementation of the terrestrial nitrogen cycle in the Joint UK Land Environment Simulator (JULES) – the land surface scheme of the UK Earth System Model (UKESM). Two configurations are discussed – the first one (JULES-CN) has a bulk soil biogeochemical model and the second one is a development configuration that resolves the soil biogeochemistry with depth (JULES-CN_layer). In JULES the nitrogen (N) cycle is based on the existing carbon (C) cycle and represents all the key terrestrial N processes in a parsimonious way. Biological N fixation is dependent on net primary productivity, and N deposition is specified as an external input. Nitrogen leaves the vegetation and soil system via leaching and a bulk gas loss term. Nutrient limitation reduces carbon-use efficiency (CUE – ratio of net to gross primary productivity) and can slow soil decomposition. We show that ecosystem level N limitation of net primary productivity (quantified in the model by the ratio of the potential amount of C that can be allocated to growth and spreading of the vegetation compared with the actual amount achieved in its natural state) falls at the lower end of the observational estimates in forests (approximately 1.0 in the model compared with 1.01 to 1.38 in the observations). The model shows more N limitation in the tropical savanna and tundra biomes, consistent with the available observations. Simulated C and N pools and fluxes are comparable to the limited available observations and model-derived estimates. The introduction of an N cycle improves the representation of interannual variability of global net ecosystem exchange, which was more pronounced in the C-cycle-only versions of JULES (JULES-C) than shown in estimates from the Global Carbon Project. It also reduces the present-day CUE from a global mean value of 0.45 for JULES-C to 0.41 for JULES-CN and 0.40 for JULES-CN_layer, all of which fall within the observational range. The N cycle also alters the response of the C fluxes over the 20th century and limits the CO₂ fertilisation effect, such that the simulated current-day land C sink is reduced by about 0.5 Pg C yr⁻¹ compared to the version with no N limitation. JULES-CN_layer additionally improves the representation of soil biogeochemistry, including turnover times in the northern high latitudes. The inclusion of a prognostic land N scheme marks a step forward in functionality and realism for the JULES and UKESM models.

Variability in climate exerts a strong influence on vegetation productivity (gross primary productivity; GPP), and therefore has a large impact on the land carbon sink. However, no direct observations of global GPP exist, and estimates rely on models that are constrained by observations at various spatial and temporal scales. Here, we assess the consistency in GPP from global products which extend for more than three decades; two observation‐based approaches, the upscaling of FLUXNET site observations (FLUXCOM) and a remote sensing derived light use efficiency model (RS‐LUE), and from a suite of terrestrial biosphere models (TRENDYv6). At local scales, we find high correlations in annual GPP among the products, with exceptions in tropical and high northern latitudes. On longer time scales, the products agree on the direction of trends over 58% of the land, with large increases across northern latitudes driven by warming trends. Further, tropical regions exhibit the largest interannual variability in GPP, with both rainforests and savannas contributing substantially. Variability in savanna GPP is likely predominantly driven by water availability, although temperature could play a role via soil moisture‐atmosphere feedbacks. There is, however, no consensus on the magnitude and driver of variability of tropical forests, which suggest uncertainties in process representations and underlying observations remain. These results emphasize the need for more direct long‐term observations of GPP along with an extension of in situ networks in underrepresented regions (e.g., tropical forests). Such capabilities would support efforts to better validate relevant processes in models, to more accurately estimate GPP.