Abstract

Abstract. Data assimilation integrates information from observational measurements with numerical models. When used with coupled models of Earth system compartments, e.g., the atmosphere and the ocean, consistent joint states can be estimated. A common approach for data assimilation is ensemble-based methods which utilize an ensemble of state realizations to estimate the state and its uncertainty. These methods are far more costly to compute than a single coupled model because of the required integration of the ensemble. However, with uncoupled models, the ensemble methods also have been shown to exhibit a particularly good scaling behavior. This study discusses an approach to augment a coupled model with data assimilation functionality provided by the Parallel Data Assimilation Framework (PDAF). Using only minimal changes in the codes of the different compartment models, a particularly efficient data assimilation system is generated that utilizes parallelization and in-memory data transfers between the models and the data assimilation functions and hence avoids most of the file reading and writing, as well as model restarts during the data assimilation process. This study explains the required modifications to the programs with the example of the coupled atmosphere–sea-ice–ocean model AWI-CM (AWI Climate Model). Using the case of the assimilation of oceanic observations shows that the data assimilation leads only to small overheads in computing time of about 15 % compared to the model without data assimilation and a very good parallel scalability. The model-agnostic structure of the assimilation software ensures a separation of concerns in which the development of data assimilation methods can be separated from the model application.

Highlights

  • Data assimilation (DA) methods are used to combine observational information with models

  • Ensemble DA (EnDA) methods use an ensemble of model state realizations to represent the state estimate and the uncertainty of this estimate given by the ensemble spread

  • The error-subspace transform Kalman filter (ESTKF) is an efficient formulation of the ensemble-based Kalman filters (EnKFs) that has been applied in different studies to assimilate satellite data into sea-ice–ocean models (e.g., Kirchgessner et al, 2017; Mu et al, 2018; Androsov et al, 2019) and biogeochemical ocean models (e.g., Pradhan et al, 2019; Goodliff et al, 2019)

Read more

Summary

Introduction

Data assimilation (DA) methods are used to combine observational information with models. The core part of the filter, which computes the corrected state vector (the socalled “analysis state”) taking into account the observational information, does not need to know how the state vector is constructed This property is important for coupled DA, where the state vector will be distributed over different compartments, such as the atmosphere and the ocean. Karspeck et al (2018) have discussed a coupled atmosphere–ocean DA system They apply the DART software and perform weakly coupled DA using two separate ensemble-based filters for the ocean and atmosphere, which produce restart files for each model compartment. These are used to initialize the ensemble integration of the coupled model.

Ensemble filters
Filter algorithms
The ESTKF
Weakly coupled ensemble filtering
Strongly coupled ensemble filtering
Setup of data assimilation program
Augmenting a coupled model for ensemble data assimilation
Parallelization for coupled ensemble data assimilation
Call-back routines for handling of model fields and observations
Scalability
Performance tuning
Application example
Findings
Discussion
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call