A mean field approximation in data assimilation for nonlinear dynamics

Gregory L. Eyink, Juan M. Restrepo, Francis J. Alexander

Research output: Contribution to journalArticlepeer-review

28 Scopus citations

Abstract

This paper considers the problem of data assimilation into nonlinear stochastic dynamic equations, from the point of view that the optimal solution is provided by the probabilities conditioned upon observations. An implementation of Bayes formula is described to calculate such probabilities. In the context of a simple model with multimodal statistics, it is shown that the conditional statistics succeed in tracking mode transitions where some standard suboptimal estimators fail. However, in complex models the exact conditional probabilities cannot be practically calculated. Instead, approximations to the conditional statistics must be sought. In this paper, attention is focused on approximations to the analysis step arising from the conditioning on observational data. A suboptimal mean-field conditional analysis is obtained from a statistical mechanics of time-histories. It is shown to have a variational formulation, reducing the approximate calculation of the conditional statistics to the minimization of the "effective action", a convex cost function. This mean-field analysis is compared with a standard linear analysis, based on a Kalman gain matrix. In the simple model problem, the mean-field conditional analysis is shown to approximate well the exact conditional statistics. Published by Elsevier B.V.

Original languageEnglish
Pages (from-to)347-368
Number of pages22
JournalPhysica D: Nonlinear Phenomena
Volume195
Issue number3-4
DOIs
StatePublished - Aug 15 2004
Externally publishedYes

Funding

We wish to thank C.D. Levermore and M. Ghil in particular for much useful advice on the subject of this paper and, also, the anonymous reviewers of a previous version of the paper. This work, LAUR 02-3058, was carried out in part at Los Alamos National Laboratory under the auspices of the Department of Energy and supported by LDRD-16 2000047. We also received support from NSF/ITR, grant DMS-0113649 (GLE, JMR), as well as from DOE DE-FG02-021625533 (JMR).

Keywords

  • Bayes formula
  • Data assimilation
  • Mean field

Fingerprint

Dive into the research topics of 'A mean field approximation in data assimilation for nonlinear dynamics'. Together they form a unique fingerprint.

Cite this