Saved in:
| Main Authors: | , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2503.00253 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915551941492736 |
|---|---|
| author | Armstrong, Robin Grooms, Ian |
| author_facet | Armstrong, Robin Grooms, Ian |
| contents | Geoscientific applications of ensemble Kalman filters face several computational challenges arising from the high dimensionality of the forecast covariance matrix, particularly when this matrix incorporates localization. For square-root filters, updating the perturbations of the ensemble members from their mean is an especially challenging step, one which generally requires approximations that introduce a trade-off between accuracy and computational cost. This paper describes an ensemble square-root filter which achieves a favorable trade-off between these factors by discretizing an integral representation of the Kalman filter update equations, and in doing so, avoids a direct evaluation of the matrix square-root in the perturbation update stage. This algorithm, which we call InFo-ESRF ("Integral-Form Ensemble Square-Root Filter"), is parallelizable and uses a preconditioned Krylov method to update perturbations to a high degree of accuracy. Through numerical experiments with both a Gaussian forecast model and a multi-layer Lorenz-type system, we demonstrate that InFo-ESRF is competitive or superior to several existing localized square-root filters in terms of accuracy and cost. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2503_00253 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Data Assimilation With An Integral-Form Ensemble Square-Root Filter Armstrong, Robin Grooms, Ian Computational Physics Geoscientific applications of ensemble Kalman filters face several computational challenges arising from the high dimensionality of the forecast covariance matrix, particularly when this matrix incorporates localization. For square-root filters, updating the perturbations of the ensemble members from their mean is an especially challenging step, one which generally requires approximations that introduce a trade-off between accuracy and computational cost. This paper describes an ensemble square-root filter which achieves a favorable trade-off between these factors by discretizing an integral representation of the Kalman filter update equations, and in doing so, avoids a direct evaluation of the matrix square-root in the perturbation update stage. This algorithm, which we call InFo-ESRF ("Integral-Form Ensemble Square-Root Filter"), is parallelizable and uses a preconditioned Krylov method to update perturbations to a high degree of accuracy. Through numerical experiments with both a Gaussian forecast model and a multi-layer Lorenz-type system, we demonstrate that InFo-ESRF is competitive or superior to several existing localized square-root filters in terms of accuracy and cost. |
| title | Data Assimilation With An Integral-Form Ensemble Square-Root Filter |
| topic | Computational Physics |
| url | https://arxiv.org/abs/2503.00253 |