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 March 15, 2023 | public
Journal Article

Deep learning-enhanced ensemble-based data assimilation for high-dimensional nonlinear dynamical systems

Abstract

Data assimilation (DA) is a key component of many forecasting models in science and engineering. DA allows one to estimate better initial conditions using an imperfect dynamical model of the system and noisy/sparse observations available from the system. Ensemble Kalman filter (EnKF) is a DA algorithm that is widely used in applications involving high-dimensional nonlinear dynamical systems. However, EnKF requires evolving large ensembles of forecasts using the dynamical model of the system. This often becomes computationally intractable, especially when the number of states of the system is very large, e.g., for weather prediction. With small ensembles, the estimated background error covariance matrix in the EnKF algorithm suffers from sampling error, leading to an erroneous estimate of the analysis state (initial condition for the next forecast cycle). In this work, we propose hybrid ensemble Kalman filter (H-EnKF), which is applied to a two-layer quasi-geostrophic turbulent flow as a test case. This framework utilizes a pre-trained deep learning-based data-driven surrogate that inexpensively generates and evolves a large data-driven ensemble of the states to accurately compute the background error covariance matrix with smaller sampling errors. The H-EnKF framework outperforms EnKF with only dynamical model or only the data-driven surrogate, and estimates a better initial condition without the need for any ad-hoc localization strategies. H-EnKF can be extended to any ensemble-based DA algorithm, e.g., particle filters, which are currently too expensive to use for high-dimensional systems.

Additional Information

© 2023 Elsevier. We thank Matti Morzfeld and Yonquiang Sun for insightful comments and discussions. This work was supported by an award from the ONR Young Investigator Program (N00014-20-1-2722), a grant from the NSF CSSI program (OAC-2005123), and NASA grant 80NSSC17K0266 to P.H. Computational resources were provided by NSF XSEDE (allocation ATM170020) to use Bridges GPU and the Rice University Center for Research Computing. The codes for H-EnKF are publicly available at https://github.com/ashesh6810/Hybrid-Ensemble-Kalman-Filter. CRediT authorship contribution statement. AC and PH formulated the problem. AC and EN wrote the codes and performed the analysis. All the authors discussed and analyzed the results. All authors contributed to the writing of the manuscript. Data availability. We have already shared all our codes and data in the manuscript. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional details

Created:
August 22, 2023
Modified:
October 25, 2023