Time-parallel incremental 4D-Var for two-dimensional rotating shallow water system

[thumbnail of Open Access]
Preview
Text (Open Access)
- Published Version
· Available under License Creative Commons Attribution.
[thumbnail of jcp_pint_da.pdf]
Text
- Accepted Version
· 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

Bhatt, R. ORCID: https://orcid.org/0000-0001-7061-1817, Debreu, L. and Vidard, A. (2026) Time-parallel incremental 4D-Var for two-dimensional rotating shallow water system. Journal of Computational Physics, 564. 115164. ISSN 1090-2716 doi: 10.1016/j.jcp.2026.115164

Abstract/Summary

Incremental four-dimensional variational data assimilation (4D-Var) relies on repeated serial integrations of tangent-linear and adjoint models within its inner loop to minimize a cost function. While classical domain decomposition offers some spatial parallelism, it faces scalability issues as model resolution and complexity increase, motivating the need for effective time-parallel strategies. This study presents a practical approach for implementing a time-parallel incremental 4D-Var based on the Parareal algorithm applied to a two-dimensional rotating shallow water system. A Krylov-enhanced Parareal variant is introduced to improve convergence for relatively more hyperbolic problems and is coupled with an inexact conjugate gradient (CG) method that adaptively adjusts Parareal precision during minimization. Numerical experiments show that the proposed implementation attains the same minimization accuracy as standard CG while achieving a theoretical time-parallel speedup of about 6.8 when communication and implementation-specific overheads are neglected. In addition, the adaptive Parareal tolerance reduces the total Parareal iterations by roughly half relative to when standard CG is used with a fixed Parareal tolerance. Further inexact CG experiments under realistic observational settings demonstrate measurable wall-clock speedups and indicate that communication overhead and sequential adjoint integrations remain important limitations on achievable parallel efficiency. Diagnostics based on geostrophic balance confirm the dynamical consistency of the retrieved initial state, demonstrating the approach’s feasibility for moderately sized PDE-constrained optimization problems.

Altmetric Badge

Dimensions Badge

Item Type Article
URI https://centaur.reading.ac.uk/id/eprint/130875
Identification Number/DOI 10.1016/j.jcp.2026.115164
Refereed Yes
Divisions Science > School of Mathematical, Physical and Computational Sciences > Department of Meteorology
Publisher Elsevier
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