Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published January 2023 | Published
Journal Article Open

A Multi-Model Ensemble Kalman Filter for Data Assimilation and Forecasting

  • 1. ROR icon California Institute of Technology
  • 2. ROR icon École Normale Supérieure
  • 3. ROR icon University of California, Los Angeles

Abstract

Data assimilation (DA) aims to optimally combine model forecasts and observations that are both partial and noisy. Multi‐model DA generalizes the variational or Bayesian formulation of the Kalman filter, and we prove that it is also the minimum variance linear unbiased estimator. Here, we formulate and implement a multi‐model ensemble Kalman filter (MM‐EnKF) based on this framework. The MM‐EnKF can combine multiple model ensembles for both DA and forecasting in a flow‐dependent manner; it uses adaptive model error estimation to provide matrix‐valued weights for the separate models and the observations. We apply this methodology to various situations using the Lorenz96 model for illustration purposes. Our numerical experiments include multiple models with parametric error, different resolved scales, and different fidelities. The MM‐EnKF results in significant error reductions compared to the best model, as well as to an unweighted multi‐model ensemble, with respect to both probabilistic and deterministic error metrics.

Copyright and License

Acknowledgement

We thank Marc Bocquet for several helpful suggestions, V. Balaji for discussions on correlated model error, Safa Mote for helpful discussions regarding hybrid methods, Tapio Schneider for discussions on smoothing applications, and four anonymous referees for additional suggestions. E.B. was funded by the Make Our Planet Great Again (MOPGA) postdoctoral program of the French Ministry for Europe and Foreign Affairs. The present work is TiPES contribution #142; the TiPES (Tipping Points in the Earth System) project has received funding from the European Union's Horizon 2020 research and innovation program under Grant Agreement No. 820970. M.G. acknowledges support by the EIT Climate-KIC; EIT Climate-KIC is supported by the European Institute of Innovation & Technology (EIT), a body of the European Union.

Funding

E.B. was funded by the Make Our Planet Great Again (MOPGA) postdoctoral program of the French Ministry for Europe and Foreign Affairs. The present work is TiPES contribution #142; the TiPES (Tipping Points in the Earth System) project has received funding from the European Union's Horizon 2020 research and innovation program under Grant Agreement No. 820970. M.G. acknowledges support by the EIT Climate-KIC; EIT Climate-KIC is supported by the European Institute of Innovation & Technology (EIT), a body of the European Union.

Code Availability

Version 2022-12 of the Julia code implementing the MM-EnKF used in this manuscript is preserved at Bach (2022), available via the MIT License and developed openly at https://github.com/eviatarbach/mmda. No data was used in this study. Scripts for numerical experiments are available in the MM-EnKF repository.

Errata

In the originally published version of this article, a few symbols in the first column of Table 1 were not boldfaced as intended. The boldface has been added to bm, B, E{f,a}, K, ρ, P{f,a}, Q, R, x{t,f,a} and X{f,a}, and this may be considered the authoritative version of record.

Files

J Adv Model Earth Syst - 2023 - Bach - A Multi‐Model Ensemble Kalman Filter for Data Assimilation and Forecasting.pdf

Additional details

Created:
September 5, 2024
Modified:
October 25, 2024