The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.

The concept of feedback is key in assessing whether a perturbation to a system is amplified or damped by mechanisms internal to the system. In polar regions, climate dynamics are controlled by both radiative and non-radiative interactions between the atmosphere, ocean, sea ice, ice sheets and land surfaces. Precisely quantifying polar feedbacks is required for a process-oriented evaluation of climate models, a clear understanding of the processes responsible for polar climate changes, and a reduction in uncertainty associated with model projections. This quantification can be performed using a simple and consistent approach that is valid for a wide range of feedbacks, offering the opportunity for more systematic feedback analyses and a better understanding of polar climate changes.

Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.

Climate tipping elements are large-scale subsystems of the Earth that may transgress critical thresholds (tipping points) under ongoing global warming, with substantial impacts on the biosphere and human societies. Frequently studied examples of such tipping elements include the Greenland Ice Sheet, the Atlantic Meridional Overturning Circulation (AMOC), permafrost, monsoon systems, and the Amazon rainforest. While recent scientific efforts have improved our knowledge about individual tipping elements, the interactions between them are less well understood. Also, the potential of individual tipping events to induce additional tipping elsewhere or stabilize other tipping elements is largely unknown. Here, we map out the current state of the literature on the interactions between climate tipping elements and review the influences between them. To do so, we gathered evidence from model simulations, observations, and conceptual understanding, as well as examples of paleoclimate reconstructions where multi-component or spatially propagating transitions were potentially at play. While uncertainties are large, we find indications that many of the interactions between tipping elements are destabilizing. Therefore, we conclude that tipping elements should not only be studied in isolation, but also more emphasis has to be put on potential interactions. This means that tipping cascades cannot be ruled out on centennial to millennial timescales at global warming levels between 1.5 and 2.0 ^∘C or on shorter timescales if global warming surpassed 2.0 ^∘C. At these higher levels of global warming, tipping cascades may then include fast tipping elements such as the AMOC or the Amazon rainforest. To address crucial knowledge gaps in tipping element interactions, we propose four strategies combining observation-based approaches, Earth system modeling expertise, computational advances, and expert knowledge.

A multimodel, multiresolution ensemble using Coupled Model Intercomparison Project Phase 6 (CMIP6) High Resolution Model Intercomparison Project (HighResMIP) coupled experiments is used to assess the performance of key aspects of the North Atlantic circulation. The Atlantic Meridional Overturning Circulation (AMOC), and related heat transport, tends to become stronger as ocean model resolution is enhanced, better agreeing with observations at 26.5°N. However, for most models the circulation remains too shallow compared to observations and has a smaller temperature contrast between the northward and southward limbs of the AMOC. These biases cause the northward heat transport to be systematically too low for a given overturning strength. The higher‐resolution models also tend to have too much deep mixing in the subpolar gyre. In the period 2015–2050 the overturning circulation tends to decline more rapidly in the higher‐resolution models, which is related to both the mean state and to the subpolar gyre contribution to deep water formation. The main part of the decline comes from the Florida Current component of the circulation. Such large declines in AMOC are not seen in the models with resolutions more typically used for climate studies, suggesting an enhanced risk for Northern Hemisphere climate change. However, only a small number of different ocean models are included in the study.

The Earth System Model Evaluation Tool (ESMValTool) is a community diagnostics and performance metrics tool designed to improve comprehensive and routine evaluation of Earth system models (ESMs) participating in the Coupled Model Intercomparison Project (CMIP). It has undergone rapid development since the first release in 2016 and is now a well-tested tool that provides end-to-end provenance tracking to ensure reproducibility. It consists of (1) an easy-to-install, well-documented Python package providing the core functionalities (ESMValCore) that performs common preprocessing operations and (2) a diagnostic part that includes tailored diagnostics and performance metrics for specific scientific applications. Here we describe large-scale diagnostics of the second major release of the tool that supports the evaluation of ESMs participating in CMIP Phase 6 (CMIP6). ESMValTool v2.0 includes a large collection of diagnostics and performance metrics for atmospheric, oceanic, and terrestrial variables for the mean state, trends, and variability. ESMValTool v2.0 also successfully reproduces figures from the evaluation and projections chapters of the Intergovernmental Panel on Climate Change (IPCC) Fifth Assessment Report (AR5) and incorporates updates from targeted analysis packages, such as the NCAR Climate Variability Diagnostics Package for the evaluation of modes of variability, the Thermodynamic Diagnostic Tool (TheDiaTo) to evaluate the energetics of the climate system, as well as parts of AutoAssess that contains a mix of top–down performance metrics. The tool has been fully integrated into the Earth System Grid Federation (ESGF) infrastructure at the Deutsches Klimarechenzentrum (DKRZ) to provide evaluation results from CMIP6 model simulations shortly after the output is published to the CMIP archive. A result browser has been implemented that enables advanced monitoring of the evaluation results by a broad user community at much faster timescales than what was possible in CMIP5.