Randomised preconditioning for the forcing formulation of weak constraint 4D‐Var

[thumbnail of Open Access]
Preview
Text (Open Access) - Published Version
· Available under License Creative Commons Attribution.
· Please see our End User Agreement before downloading.
| Preview
Available under license: Creative Commons Attribution
[thumbnail of Dauzickate2021b.pdf]
Text - Accepted Version
· Restricted to Repository staff only
Restricted to Repository staff only

Please see our End User Agreement.

It is advisable to refer to the publisher's version if you intend to cite from this work. See Guidance on citing.

Add to AnyAdd to TwitterAdd to FacebookAdd to LinkedinAdd to PinterestAdd to Email

Daužickaitė, I. orcid id iconORCID: https://orcid.org/0000-0002-1285-1764, Lawless, A. S. orcid id iconORCID: https://orcid.org/0000-0002-3016-6568, Scott, J. A. orcid id iconORCID: https://orcid.org/0000-0003-2130-1091 and Van Leeuwen, P. J. orcid id iconORCID: https://orcid.org/0000-0003-2325-5340 (2021) Randomised preconditioning for the forcing formulation of weak constraint 4D‐Var. Quarterly Journal of the Royal Meteorological Society, 147 (740). pp. 3719-3734. ISSN 1477-870X doi: 10.1002/qj.4151

Abstract/Summary

There is growing awareness that errors in the model equations cannot be ignored in data assimilation methods such as four-dimensional variational assimilation (4D-Var). If allowed for, more information can be extracted from observations, longer time windows are possible, and the minimisation process is easier, at least in principle. Weak constraint 4D-Var estimates the model error and minimises a series of linear least-squares cost functions, which can be achieved using the conjugate gradient (CG) method; minimising each cost function is called an inner loop. CG needs preconditioning to improve its performance. In previous work, limited memory preconditioners (LMPs) have been constructed using approximations of the eigenvalues and eigenvectors of the Hessian in the previous inner loop. If the Hessian changes significantly in consecutive inner loops, the LMP may be of limited usefulness. To circumvent this, we propose using randomised methods for low rank eigenvalue decomposition and use these approximations to cheaply construct LMPs using information from the current inner loop. Three randomised methods are compared. Numerical experiments in idealized systems show that the resulting LMPs perform better than the existing LMPs. Using these methods may allow more efficient and robust implementations of incremental weak constraint 4D-Var.

Altmetric Badge

Item Type Article
URI https://reading-clone.eprints-hosting.org/id/eprint/99965
Identification Number/DOI 10.1002/qj.4151
Refereed Yes
Divisions Science > School of Mathematical, Physical and Computational Sciences > Department of Mathematics and Statistics
Publisher Royal Meteorological Society
Download/View statistics View download statistics for this item

Downloads

Downloads per month over past year

University Staff: Request a correction | Centaur Editors: Update this record

Search Google Scholar