Efficient methods for computing observation impact in 4D-Var data assimilation

Files
TR Number
Date
2013-12-01
Journal Title
Journal ISSN
Volume Title
Publisher
Springer
Abstract

This paper presents a practical computational approach to quantify the effect of individual observations in estimating the state of a system. Such an analysis can be used for pruning redundant measurements, and for designing future sensor networks.

The mathematical approach is based on computing the sensitivity of the reanalysis (unconstrained optimization solution) with respect to the data. The computational cost is dominated by the solution of a linear system, whose matrix is the Hessian of the cost function, and is only available in operator form. The right hand side is the gradient of a scalar cost function that quantities the forecast error of the numerical model. The use of adjoint models to obtain the necessary first and second order derivatives is discussed. We study various strategies to accelerate the computation, including matrix-free iterative solvers, preconditioners, and an in-house multigrid solver. Experiments are conducted on both a small-size shallow-water equations model, and on a large-scale numerical weather prediction model, in order to illustrate the capabilities of the new methodology.

Description
Keywords
Technology, Computer Science, Interdisciplinary Applications, Geosciences, Multidisciplinary, Computer Science, Geology, Variational data assimilation, Sensitivity analysis, Observation impact, Matrix-free solvers, Preconditioners, VARIATIONAL DATA ASSIMILATION, ENSEMBLE DATA ASSIMILATION, ADJOINT SENSITIVITY, ADAPTIVE OBSERVATIONS, TARGETED OBSERVATIONS, KALMAN FILTER, SYSTEM, MODELS, FORECASTS, METRICS
Citation