The ECMWF OCEAN5 system is a global ocean and sea-ice ensemble of reanalysis and real-time analysis. This paper gives a full description of the OCEAN5 system, with the focus on upgrades of system components with respect to its predecessors, ORAS4 and ORAP5. An important novelty in OCEAN5 is the ensemble generation strategy that includes perturbation of initial conditions and a generic perturbation scheme for observations and forcing fields. Other upgrades include revisions to the a priori bias correction scheme, observation quality control and assimilation method for sea-level anomalies. The OCEAN5 historical reconstruction of the ocean and sea-ice state is the ORAS5 reanalysis, which includes five ensemble members and covers the period from 1979 onwards. Updated versions of observation data sets are used in ORAS5 production, with special attention devoted to the consistency of sea surface temperature (SST) and sea-ice observations. Assessment of ORAS5 through sensitivity experiments suggests that all system components contribute to an improved fit to observation in reanalyses, with the most prominent contribution from direct assimilation of ocean in situ observations. Results of observing system experiments further suggest that the Argo float is the most influential observation type in our data assimilation system. Assessment of ORAS5 has also been carried out for several key ocean state variables and verified against reference climate data sets from the ESA CCI (European Space Agency Climate Change Initiative) project. With respect to ORAS4, ORAS5 has improved ocean climate state and variability in terms of SST and sea level, mostly due to increased model resolution and updates in assimilated observation data sets. In spite of the improvements, ORAS5 still underestimates the temporal variance of sea level and continues exhibiting large SST biases in the Gulf Stream and its extension regions which are possibly associated with misrepresentation of front positions. Overall, the SST and sea-ice uncertainties estimated using five ORAS5 ensemble members have spatial patterns consistent with those of analysis error. The ensemble spread of sea ice is commensurable with the sea-ice analysis error. On the contrary, the ensemble spread is under-dispersive for SST.

We present an investigation of the performance of the ECMWF coupled atmosphere‐waves‐ocean model for different ocean and atmosphere resolutions on a series of tropical cyclones in the Western Pacific with the aim to better understand the coupled feedback mechanisms in these extreme conditions. For some of the test cases, we only find little impact of coupling the atmosphere to the ocean, while in others, we observe a very large impact. To further understand these differences, we have selected two tropical cyclones (TCs) as case studies: TC Haiyan (with small impact of coupling) and TC Neoguri (with large impact of coupling). The comparison between these two cases suggests that the upper ocean stratification is the key in determining the strength of the coupled feedback. A strong coupled feedback is found whenever the ocean heat content of the upper layer is low while a very weak coupled feedback is found whenever the ocean has a thick warm mixed layer. The oceanographic response to tropical cyclones for the two storms has been compared to sea surface temperature and derived surface currents from drifting buoys and to subsurface observations from Argo and ship launched XBT's. These comparisons show that we are able to realistically reproduce the atmospheric and oceanographic interaction during tropical cyclone conditions which gives us confidence that the coupled modeling system is physically sound.

We report on the first multi-year kilometre-scale global coupled simulations using ECMWF's Integrated Forecasting System (IFS) coupled to both the NEMO and FESOM ocean–sea ice models, as part of the H2020 Next Generation Earth Modelling Systems (nextGEMS) project. We focus mainly on an unprecedented IFS-FESOM coupled setup, with an atmospheric resolution of 4.4 km and a spatially varying ocean resolution that reaches locally below 5 km grid spacing. A shorter coupled IFS-FESOM simulation with an atmospheric resolution of 2.8 km has also been performed. A number of shortcomings in the original numerical weather prediction (NWP)-focused model configurations were identified and mitigated over several cycles collaboratively by the modelling centres, academia, and the wider nextGEMS community. The main improvements are (i) better conservation properties of the coupled model system in terms of water and energy budgets, which also benefit ECMWF's operational 9 km IFS-NEMO model; (ii) a realistic top-of-the-atmosphere (TOA) radiation balance throughout the year; (iii) improved intense precipitation characteristics; and (iv) eddy-resolving features in large parts of the mid- and high-latitude oceans (finer than 5 km grid spacing) to resolve mesoscale eddies and sea ice leads. New developments at ECMWF for a better representation of snow and land use, including a dedicated scheme for urban areas, were also tested on multi-year timescales. We provide first examples of significant advances in the realism and thus opportunities of these kilometre-scale simulations, such as a clear imprint of resolved Arctic sea ice leads on atmospheric temperature, impacts of kilometre-scale urban areas on the diurnal temperature cycle in cities, and better propagation and symmetry characteristics of the Madden–Julian Oscillation.

The Next Generation of Earth Modeling Systems (nextGEMS) project aimed to produce multidecadal climate simulations, for the first time, with resolved kilometer-scale (km-scale) processes in the ocean, land, and atmosphere. In only 3 years, nextGEMS achieved this milestone with the two km-scale Earth system models, ICOsahedral Non-hydrostatic model (ICON) and Integrated Forecasting System coupled to the Finite-volumE Sea ice-Ocean Model (IFS-FESOM). nextGEMS was based on three cornerstones: (1) developing km-scale Earth system models with small errors in the energy and water balance, (2) performing km-scale climate simulations with a throughput greater than 1 simulated year per day, and (3) facilitating new workflows for an efficient analysis of the large simulations with common data structures and output variables. These cornerstones shaped the timeline of nextGEMS, divided into four cycles. Each cycle marked the release of a new configuration of ICON and IFS-FESOM, which were evaluated at hackathons. The hackathon participants included experts from climate science, software engineering, and high-performance computing as well as users from the energy and agricultural sectors. The continuous efforts over the four cycles allowed us to produce 30-year simulations with ICON and IFS-FESOM, spanning the period 2020–2049 under the SSP3-7.0 scenario. The throughput was about 500 simulated days per day on the Levante supercomputer of the German Climate Computing Center (DKRZ). The simulations employed a horizontal grid of about 5 km resolution in the ocean and 10 km resolution in the atmosphere and land. Aside from this technical achievement, the simulations allowed us to gain new insights into the realism of ICON and IFS-FESOM. Beyond its time frame, nextGEMS builds the foundation of the Climate Change Adaptation Digital Twin developed in the Destination Earth initiative and paves the way for future European research on climate change.

The drive to develop environmental prediction systems that are seamless across both weather and climate timescales has culminated in the development and use of Earth system models, which include a coupled representation of the atmosphere, land, ocean and sea ice, for medium-range weather forecasts. One region where such a coupled Earth system approach has the potential to significantly influence the skill of weather forecasts is in the polar and sub-polar seas, where fluxes of heat, moisture and momentum are strongly influenced by the position of the sea ice edge. In this study we demonstrate that using a dynamically coupled ocean and sea ice model in the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System results in improved sea ice edge position forecasts in the Northern Hemisphere in the medium range. Further, this improves forecasts of boundary layer temperature and humidity downstream of the sea ice edge in some regions during periods of rapid change in the sea ice, compared to forecasts in which the sea surface temperature anomalies and sea ice concentration do not evolve throughout the forecasts. However, challenges remain, such as large errors in the position of the ice edge in the ocean analysis used to initialise the ocean component of the coupled system, which has an error of approximately 50 % of the total forecast error at day 9, suggesting there is much skill to be gained by improving the ocean analysis at and around the sea ice edge. The importance of the choice of sea ice analysis for verification is also highlighted, with a call for more guidance on the suitability of satellite sea ice products to verify forecasts on daily to weekly timescales and on meso-scales (< 500 km).

The ECMWF OCEAN5 system is a global ocean and sea-ice ensemble of reanalysis and real-time analysis. This manuscript gives a full description of the OCEAN5 system, with the focus on upgrades of system components with respect to its predecessors ORAS4 and ORAP5. An important novelty in OCEAN5 is the ensemble generation strategy that includes perturbation of initial conditions, and a generic perturbation scheme for observations and forcing fields. Other upgrades include revisions to the a-priori bias correction scheme, observation quality control and assimilation method for sea-level anomaly. The OCEAN5 historical reconstruction of the ocean and sea-ice state is the ORAS5 reanalysis, which includes 5 ensemble members and covers the period from 1979 onwards, and with a backward extension until 1958. Updated version of observation data sets are used in ORAS5 production, with special attention devoted to the consistency of sea surface temperature (SST) and sea-ice observations. Assessment of ORAS5 in the observation space suggests that assimilation of observations contribute to reducing the analysis error, with the most prominent contribution from direct assimilation of ocean in-situ observations. Results of observing system experiment further suggest that Argo float is the most influential observation type in our data assimilation system. Assessment of ORAS5 has also been carried out with several key ocean state variables and verified against independent observation data sets from ESA CCI project. With respect to ORAS4, ORAS5 has improved ocean climate state and variability in terms of SST and sea-level, mostly due to increased model resolution and updates in assimilated observation data sets. In spite of the improvements, ORAS5 still underestimates the temporal variance of sea level, and continue exhibiting large SST biases in the Gulf Stream and extension regions which is possibly associated with misrepresentation of front positions. Overall, the SST and sea-ice uncertainties estimated using five ORAS5 ensemble members have spatial patterns consistent with those of analysis error. The ensemble spread of sea-ice is commensurable with the sea-ice analysis error. On the contrary, the ensemble is under-dispersive for SST.