Description of the resolution hierarchy of the global coupled HadGEM3-GC3.1 model as used in CMIP6 HighResMIP experiments

CMIP6 HighResMIP is a new experimental design for global climate model simulations that aims to assess the impact of model horizontal resolution on climate simulation fidelity. We describe a hierarchy of global coupled model resolutions based on the HadGEM3-GC3.1 model that range from an atmosphere-ocean resolution of 130 km-1° to 25 km-1/12°, all using the same forcings and initial conditions. In order to make such high resolution simulations possible, the experiments have a short 30 year spinup, followed by at least century-long simulations with both constant forcing (to assess drift and the focus of this work), and historic forcing. We assess the change in model biases as a function of both atmosphere and ocean resolution, together with the effectiveness and robustness of this new experimental design. We find reductions in the biases in top of atmosphere radiation components and cloud forcing. There are significant reductions in some common surface climate model biases as resolution is increased, particularly in the Atlantic for sea surface temperature and precipitation, primarily driven by increased ocean resolution. There is also a reduction in drift from the initial conditions both at the surface and in the deeper ocean at higher resolution. Using an eddy-present and eddy-rich ocean resolution enhances the strength of the North Atlantic ocean circulation (boundary currents, overturning circulation and heat transports), while an eddy-present ocean resolution has a considerably reduced Antarctic Circumpolar Current strength. All models have a reasonable representation of El Nino – Southern Oscillation. In general the biases present after 30 years of simulations do not change character markedly over longer timescales, justifying the experimental design

The low-resolution version of HadGEM3 GC3.1: development and evaluation for global climate

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 both to 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 north-west Atlantic, N96ORCA1 shows a cold surface bias of up to 6K, 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 2K) 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 UKESM1 and as a stand-alone coupled climate model

Resolving and parameterising the ocean mesoscale in earth system models

Purpose of Review. Assessment of the impact of ocean resolution in Earth System models on the mean state, variability, and future projections and discussion of prospects for improved parameterisations to represent the ocean mesoscale. Recent Findings. The majority of centres participating in CMIP6 employ ocean components with resolutions of about 1 degree in their full Earth Systemmodels (eddy-parameterising models). In contrast, there are alsomodels submitted toCMIP6 (both DECK and HighResMIP) that employ ocean components of approximately 1/4 degree and 1/10 degree (eddy-present and eddy-rich models). Evidence to date suggests that whether the ocean mesoscale is explicitly represented or parameterised affects not only the mean state of the ocean but also the climate variability and the future climate response, particularly in terms of the Atlantic meridional overturning circulation (AMOC) and the Southern Ocean. Recent developments in scale-aware parameterisations of the mesoscale are being developed and will be included in future Earth System models. Summary. Although the choice of ocean resolution in Earth System models will always be limited by computational considerations, for the foreseeable future, this choice is likely to affect projections of climate variability and change as well as other aspects of the Earth System. Future Earth System models will be able to choose increased ocean resolution and/or improved parameterisation of processes to capture physical processes with greater fidelity

Predicting September Arctic Sea Ice: A Multi-Model Seasonal Skill Comparison

Abstract This study quantifies the state-of-the-art in the rapidly growing field of seasonal Arctic sea ice prediction. A novel multi-model dataset of retrospective seasonal predictions of September Arctic sea ice is created and analyzed, consisting of community contributions from 17 statistical models and 17 dynamical models. Prediction skill is compared over the period 2001–2020 for predictions of Pan-Arctic sea ice extent (SIE), regional SIE, and local sea ice concentration (SIC) initialized on June 1, July 1, August 1, and September 1. This diverse set of statistical and dynamical models can individually predict linearly detrended Pan-Arctic SIE anomalies with skill, and a multi-model median prediction has correlation coefficients of 0.79, 0.86, 0.92, and 0.99 at these respective initialization times. Regional SIE predictions have similar skill to Pan-Arctic predictions in the Alaskan and Siberian regions, whereas regional skill is lower in the Canadian, Atlantic, and Central Arctic sectors. The skill of dynamical and statistical models is generally comparable for Pan-Arctic SIE, whereas dynamical models outperform their statistical counterparts for regional and local predictions. The prediction systems are found to provide the most value added relative to basic reference forecasts in the extreme SIE years of 1996, 2007, and 2012. SIE prediction errors do not show clear trends over time, suggesting that there has been minimal change in inherent sea ice predictability over the satellite era. Overall, this study demonstrates that there are bright prospects for skillful operational predictions of September sea ice at least three months in advance.</jats:p

UKESM1: description and evaluation of the UK Earth System Model

We document the development of the first version of the United Kingdom 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 idealised simulations show a high climate sensitivity relative to previous generations of models: equilibrium climate sensitivity (ECS) is 5.4 K, transient climate response (TCR) ranges from 2.68 K to 2.85 K, and transient climate response to cumulative emissions (TCRE) is 2.49 K/TtC to 2.66 K/TtC

Predicting September Arctic Sea Ice: A Multi-Model Seasonal Skill Comparison

This study quantifies the state-of-the-art in the rapidly growing field of seasonal Arctic sea ice prediction. A novel multi-model dataset of retrospective seasonal predictions of September Arctic sea ice is created and analyzed, consisting of community contributions from 17 statistical models and 17 dynamical models. Prediction skill is compared over the period 2001–2020 for predictions of Pan-Arctic sea ice extent (SIE), regional SIE, and local sea ice concentration (SIC) initialized on June 1, July 1, August 1, and September 1. This diverse set of statistical and dynamical models can individually predict linearly detrended Pan-Arctic SIE anomalies with skill, and a multi-model median prediction has correlation coefficients of 0.79, 0.86, 0.92, and 0.99 at these respective initialization times. Regional SIE predictions have similar skill to Pan-Arctic predictions in the Alaskan and Siberian regions, whereas regional skill is lower in the Canadian, Atlantic, and Central Arctic sectors. The skill of dynamical and statistical models is generally comparable for Pan-Arctic SIE, whereas dynamical models outperform their statistical counterparts for regional and local predictions. The prediction systems are found to provide the most value added relative to basic reference forecasts in the extreme SIE years of 1996, 2007, and 2012. SIE prediction errors do not show clear trends over time, suggesting that there has been minimal change in inherent sea ice predictability over the satellite era. Overall, this study demonstrates that there are bright prospects for skillful operational predictions of September sea ice at least three months in advance

Bringing it all together: science priorities for improved understanding of Earth system change and to support international climate policy

We review how the international modelling community, encompassing integrated assessment models, global and regional Earth system and climate models, and impact models, has worked together over the past few decades to advance understanding of Earth system change and its impacts on society and the environment and thereby support international climate policy. We go on to recommend a number of priority research areas for the coming decade, a timescale that encompasses a number of newly starting international modelling activities, as well as the IPCC Seventh Assessment Report (AR7) and the second UNFCCC Global Stocktake. Progress in these priority areas will significantly advance our understanding of Earth system change and its impacts, increasing the quality and utility of science support to climate policy. [...

The potential of numerical prediction systems to support the design of Arctic observing systems: insights from the APPLICATE and YOPP projects

Numerical systems used for weather and climate predictions have substantially improved over past decades. We argue that despite a continued need for further addressing remaining limitations of their key components, numerical prediction systems have reached a sufficient level of maturity to examine and critically assess the suitability of Earth's current observing systems – remote and in situ, for prediction purposes; and that they can provide evidence-based support for the deployment of future observational networks. We illustrate this point by presenting recent, co-ordinated international efforts focused on Arctic observing systems, led in the framework of the Year of Polar Prediction and the H2020 project APPLICATE. The Arctic, one of the world's most rapidly changing regions, is relatively poorly covered in terms of in situ data but richly covered in terms of satellite data. In this study, we demonstrate that existing state-of-the-art datasets and targeted sensitivity experiments produced with numerical prediction systems can inform us of the added value of existing or even hypothetical Arctic observations, in the context of predictions from hourly to interannual time-scales. Furthermore, we argue that these datasets and experiments can also inform us how the uptake of Arctic observations in numerical prediction systems can be enhanced to maximise predictive skill. Based on these efforts we suggest that (a) conventional in situ observations in the Arctic play a particularly important role in initializing numerical weather forecasts during the winter season, (b) observations from satellite microwave sounders play a particularly important role during the summer season, and their enhanced usage over snow and sea ice is expected to further improve their impact on predictive skill in the Arctic region and beyond, (c) the deployment of a small number of in situ sea-ice thickness monitoring devices at strategic sampling sites in the Arctic could be sufficient to monitor most of the large-scale sea-ice volume variability, and (d) sea-ice thickness observations can improve the simulation of both the sea ice and near-surface air temperatures on seasonal time-scales in the Arctic and beyond

An inter-comparison of the mass budget of the Arctic sea ice in CMIP6 models

We compare the mass budget of the Arctic sea ice for 15 models submitted to the latest Coupled Model Intercomparison Project (CMIP6), using new diagnostics that have not been available for previous model intercomparisons. These diagnostics allow us to look beyond the standard metrics of ice cover and thickness to compare the processes of sea ice growth and loss in climate models in a more detailed way than has previously been possible. For the 1960–1989 multi-model mean, the dominant processes causing annual ice growth are basal growth and frazil ice formation, which both occur during the winter. The main processes by which ice is lost are basal melting, top melting and advection of ice out of the Arctic. The first two processes occur in summer, while the latter process is present all year. The sea ice budgets for individual models are strikingly similar overall in terms of the major processes causing ice growth and loss and in terms of the time of year during which each process is important. However, there are also some key differences between the models, and we have found a number of relationships between model formulation and components of the ice budget that hold for all or most of the CMIP6 models considered here. The relative amounts of frazil and basal ice formation vary between the models, and the amount of frazil ice formation is strongly dependent on the value chosen for the minimum frazil ice thickness. There are also differences in the relative amounts of top and basal melting, potentially dependent on how much shortwave radiation can penetrate through the sea ice into the ocean. For models with prognostic melt ponds, the choice of scheme may affect the amount of basal growth, basal melt and top melt, and the choice of thermodynamic scheme is important in determining the amount of basal growth and top melt. As the ice cover and mass decline during the 21st century, we see a shift in the timing of the top and basal melting in the multi-model mean, with more melt occurring earlier in the year and less melt later in the summer. The amountof basal growth reduces in the autumn, but it increases in the winter due to thinner sea ice over the course of the 21st century. Overall, extra ice loss in May–June and reduced ice growth in October–November are partially offset by reduced ice melt in August and increased ice growth in January– February. For the individual models, changes in the budget components vary considerably in terms of magnitude and timing of change. However, when the evolving budget terms are considered as a function of the changing ice state itself, behaviours common to all the models emerge, suggesting that the sea ice components of the models are fundamentally responding in a broadly consistent way to the warming climate. It is possible that this similarity in the model budgets may represent a lack of diversity in the model physics of the CMIP6 models considered here. The development of new observational datasets for validating the budget terms would help to clarify this