A Multi-Model Ensemble Kalman Filter for Data Assimilation and Forecasting
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
© 2023 The Authors. Journal of Advances in Modeling Earth Systems published by Wiley Periodicals LLC on behalf of American Geophysical Union. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
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
Name | Size | Download all |
---|---|---|
md5:f0c08d5a41cdbf36f5cf5727d788f6bf
|
1.6 MB | Preview Download |
Additional details
- Ministère de l'Europe et des Affaires étrangères
- MOPGA-977406H
- European Commission
- 820970
- European Institute of Innovation and Technology
- Available
-
2023-01-19Issue Online
- Accepted
-
2022-12-28Manuscript Accepted
- Caltech groups
- Division of Geological and Planetary Sciences
- Publication Status
- Published