Spectral estimates for saddle point matrices arising in weak constraint four-dimensional variational data assimilation

TL;DR

This paper analyzes the spectral properties of saddle point matrices in weak constraint 4D-Var data assimilation, providing bounds on eigenvalues and examining their sensitivity to observation counts, with numerical validation.

Contribution

It offers a detailed spectral analysis and eigenvalue bounds for saddle point matrices in 4D-Var data assimilation, enhancing understanding of their numerical behavior.

Findings

Eigenvalue bounds for saddle point matrices

Spectral sensitivity to number of observations

Numerical experiments confirm theoretical bounds

Abstract

We consider the large-sparse symmetric linear systems of equations that arise in the solution of weak constraint four-dimensional variational data assimilation, a method of high interest for numerical weather prediction. These systems can be written as saddle point systems with a 3x3 block structure but block eliminations can be performed to reduce them to saddle point systems with a 2x2 block structure, or further to symmetric positive definite systems. In this paper, we analyse how sensitive the spectra of these matrices are to the number of observations of the underlying dynamical system. We also obtain bounds on the eigenvalues of the matrices. Numerical experiments are used to confirm the theoretical analysis and bounds.