S1: Diagnosis and Metrics in Climate Models

Principal investigators: Prof. Hans Burchard (Leibniz Institute for Baltic Sea Research Warnemünde), Prof. Veronika Eyring (MARUM/University of Bremen), Prof. Thomas Jung (Alfred Wegener Institute for Polar and Marine Research/MARUM/University of Bremen), Dr. Martin Losch (Alfred-Wegener-Institut, Helmholtz Zentrum für Polar- und Meeresforschung), Prof. Jens Rademacher (Universität Hamburg), Prof. Bjorn Stevens (Max Planck Institut für Meteorologie)

In a continuation and extension of the S1 subproject from the first phase of the CRC, S1 will assess the impact of model developments – proposed by the other project areas – on the quality of two climate models developed and used by the German scientific community, namely ICONa/ICON-o and OpenIFS/FESOM2. Major parts of the work will be organized around performing numerical simulations with new parameterizations, and development and application of new metrics, diagnostics and tools for model evaluation in terms of model biases and response to forcing.

Numerical simulations will be performed for a hierarchy of model setups, ranging from simplified channel models to realistic global coupled high-resolution simulations following CMIP6 protocols. These simulations will provide the basis for further model evaluation. Scalable open source tools for processing and evaluation of this data will be further developed and shared with the community. We will continue using diagnostics developed in the first phase, and implementing additional diagnostics already published in the literature. Moreover, we will develop new diagnostics that will help us to evaluate models in the context of their energetic consistency.

The new diagnostics will extend the validation of the scale-dependent spatial variability performed by the MODES software to the temporal domain. The extension will include a novel method that quantifies the effects of model biases on the simulated spatio-temporal variability. Remote effects of model biases – here called bias teleconnections – will be initially investigated in a GCM of intermediate complexity, then applied to ICON-a/ICON-o and OpenIFS/FESOM2 including new parameterizations.

Numerical mixing and dissipation are potential sources of large errors in ocean components of climate models. New numerical schemes developed in the CRC to tackle these shortcomings need to be properly evaluated. Therefore, methods to quantify numerical errors in water mass transformations will be further developed in M5 during the first two years and implemented and evaluated in FESOM2 during the last two years of the second phase of S1.

Results of ICON-a/ICON-o and OpenIFS/FESOM2 model simulations will be put into a broader climate modelling context (with emphasis on CMIP6 simulations). The Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for routine evaluation of climate models, will be enhanced with the offline diagnostics developed within this project. By making additional use of a wide range of diagnostics already included in the new ESMValTool version 2.0 (v2.0), the impact of model developments will be thoroughly evaluated, also in the context of CMIP6 models.

FESOM model

One of the two ocean models used in the TRR181 project is FESOM (Finite-Element/volumE Sea ice-Ocean Model), which is part of the AWI-CM climate model of the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI).

FESOM is the only ocean model participating in CMIP6, that is formulated on an unstructured mesh. This allows scientists to flexibly increase horizontal resolution of the model in more energetically active areas of the ocean like the Gulf stream or Agulhas current. The resulting model resolves important details of ocean circulation, but is still computationally efficient.

FESOM ocean currents in Indian and Pacific oceans (100m)

The Earth’s energy budget and other funny aspects of the thermodynamics of the climate system

State-of-the-art climate models still struggle to reproduce a reasonably energetically consistent system, even though outstanding improvements have been achieved in the recent past.

Valerio Lembo, Postdoc in S1

The idea of this subproject is assessing the impact of introducing new numerical schemes and physical parametrizations developed in the TRR181 for the energy closure of state-of-the-art climate models. We provide diagnostic tools that allow for evaluation and intercomparison of climate models, starting from their outputted datasets.

It might sound trivial, expecting that the climate system, if in steady state, is also in thermodynamic equilibrium. This is at least what our studies of classical thermodynamics suggest. The problem is that the system constantly exchanges energy with its exterior, i.e. the outer space, and within its interior. In steady state conditions, the net exchange of energy with the exterior has to be null. In other words, the climate system is in thermodynamic equilibrium, once we averaged out the modulation of the solar energy input to an appropriately long timescale and all the energy exchanges occurring in its interior, shaping the solar “reflection” and the thermal energy output. This is a clear example of what is called a “non-equilibrium dissipative steady state thermodynamical system”.

State-of-the-art climate models still struggle to reproduce a reasonably energetically consistent system, even though outstanding improvements have been achieved in the recent past. This points to the very basic reasons for climate modeling, on one hand reflecting the lack of understanding of some processes involving energy exchanges and the limits of the discretization/truncation of the real world in finite dimension models, on the other hand preventing us from correctly evaluating the impact of the various forcings for reconstructed and projected climate change.

As TRR181, we are participating to the community effort called “ESMValTool”, whose aim is providing a set of standardized diagnostics for the evaluation of state-of-the-art and forthcoming multi-model ensembles. In our diagnostics, we try to address specifically the Earth’s energy budget and its atmospheric and oceanic components, and the atmospheric energy exchanges, including the Lorenz Energy Cycle, which describes the energy exchanges in the extratropical synoptic eddies. We also provide an estimate of the atmospheric material entropy production, i.e. the entropy production through irreversible processes, and the water mass budget, which is known to be one of the main sources of uncertainty for the modeled energy budget.

The diagnostic tool is currently being ported from version 1 to version 2 of ESMValTool, and will be hopefully soon publicly released. A report for the ESMValTool version 2 is being written, with contributions by all groups in the community, and another paper, focused on potential applications of the tool in various fields of climate science, will be submitted.

Metrics and Diagnostics for model improvements

The proper protocol and experiment setup for numerical experiments is crucial.

Nikolay Koldunov, Postdoc in S1

I am Nikolay Koldunov, Post Doc at MARUM and Alfred Wegener Institute. Since October I begin to work at Research Project S1: Diagnosis and Metrics in Climate Models. The main aim of the project is to integrate and synthesize work done in other parts of the TRR181. In particular we will provide metrics and diagnostics to help access the impact of model improvements suggested by TRR181 on quality of the climate models. One of the main challenges is to create model diagnostics that would not only quantify improvements, but also allow to clearly identify the cause of changes in model behavior. In this respect the proper protocol and experiment setup for numerical experiments is crucial and its development will be important part of my work. The resulting diagnostics will become available for the wider research community through the ESMValTool, that is going to be one of the main instruments of model analysis for the CMIP6 project.