Abstract
A Bayesian data assimilation scheme is formulated for advection-dominated or hyperbolic evolutionary problems, and observations. It uses the physics to dynamically update the likelihood in order to extend the impact of the likelihood on the posterior, a strategy that would be particularly useful when the observation network is sparse in space and time and the associated measurement uncertainties are low. The filter is applied to a problem with linear dynamics and Gaussian statistics, and compared to the exact estimate, a model outcome, and the Kalman filter estimate. By comparing to the exact estimate the dynamic likelihood filter is shown to be superior to model outcomes and to the Kalman estimate, when the observation system is sparse. The added computational expense of the method is linear in the number of observations and thus computationally efficient, suggesting that the method is practical even if the space dimensions of the physical problem are large.
Original language | English |
---|---|
Pages (from-to) | 2915-2924 |
Number of pages | 10 |
Journal | Quarterly Journal of the Royal Meteorological Society |
Volume | 143 |
Issue number | 708 |
DOIs | |
State | Published - Oct 2017 |
Externally published | Yes |
Funding
This work was supported by Pacific Earthquake Engineering Research Center (PEER) research grant no. 1123-NCTRYH, and National Science Foundation NSF/OCE no. 1434198. I wish to thank Stockholm University Rossby Fellowship Program, for their hospitality. The referees and the handling editor provided very useful comments that were used to improve the article.
Funders | Funder number |
---|---|
National Science Foundation NSF/OCE | 1434198 |
Pacific Earthquake Engineering Research Center | 1123-NCTRYH |
Keywords
- Kalman filter
- data assimilation
- dynamic likelihood
- forecasting
- hyperbolic equation
- wave equation