Title :
An error estimate of Gaussian recursive filter in 3Dvar problem
Author :
Cuomo, Salvatore ; Farina, Raffaele ; Galletti, Ardelio ; Marcellino, Livia
Author_Institution :
Dept. of Math. & Applic. “R. Caccioppoli”, Univ. of Naples Federico II, Naples, Italy
Abstract :
Computational kernel of the three-dimensional variational data assimilation (3D-Var) problem is a linear system, generally solved by means of an iterative method. The most costly part of each iterative step is a matrix-vector product with a very large covariance matrix having Gaussian correlation structure. This operation may be interpreted as a Gaussian convolution, that is a very expensive numerical kernel. Recursive Filters (RFs) are a well known way to approximate the Gaussian convolution and are intensively applied in the meteorology, in the oceanography and in forecast models. In this paper, we deal with an oceanographic 3D-Var data assimilation scheme, named OceanVar, where the linear system is solved by using the Conjugate Gradient (GC) method by replacing, at each step, the Gaussian convolution with RFs. Here we give theoretical issues on the discrete convolution approximation with a first order (1st-RF) and a third order (3rd-RF) recursive filters. Numerical experiments confirm given error bounds and show the benefits, in terms of accuracy and performance, of the 3-rd RF.
Keywords :
Gaussian processes; conjugate gradient methods; covariance matrices; data assimilation; geophysical signal processing; iterative methods; oceanographic techniques; recursive filters; 1st-RF; 3D variational data assimilation problem; 3Dvar problem; 3rd-RF; Gaussian convolution; Gaussian correlation structure; Gaussian recursive filter; OceanVar; conjugate gradient method; covariance matrix; discrete convolution approximation; first order recursive filters; iterative method; linear system; matrix-vector product; numerical kernel; oceanographic 3D-Var data assimilation scheme; third order recursive filters; Approximation methods; Convolution; Covariance matrices; Data assimilation; Linear systems; Radio frequency; Vectors;
Conference_Titel :
Computer Science and Information Systems (FedCSIS), 2014 Federated Conference on
Conference_Location :
Warsaw