Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model

Citation
P. Brasseur et al., Assimilation of altimetric data in the mid-latitude oceans using the Singular Evolutive Extended Kalman filter with an eddy-resolving, primitive equation model, J MAR SYST, 22(4), 1999, pp. 269-294
Citations number
43
Categorie Soggetti
Aquatic Sciences","Earth Sciences
Journal title
JOURNAL OF MARINE SYSTEMS
ISSN journal
09247963 → ACNP
Volume
22
Issue
4
Year of publication
1999
Pages
269 - 294
Database
ISI
SICI code
0924-7963(199911)22:4<269:AOADIT>2.0.ZU;2-Z
Abstract
A new data assimilation scheme has been elaborated for ocean circulation mo dels based on the concept of an evolutive, reduced-order Kalman filter. The dimension of the assimilation problem is reduced by expressing the initial error covariance matrix as a truncated series of orthogonal perturbations. This error sub-space evolves during the assimilation so as to capture the growing modes of the estimation error. The algorithm has been formulated in quite a general fashion to make it tractable with a large variety of ocean models and measurement types. In the present paper, we have examined three possible strategies to compute the evolution of the error subspace in the so-called Singular Evolutive Extended Kalman (SEEK) filter: the steady filt er considers a time-independent error sub-space, the apprentice filter prog ressively enriches the error sub-space with the information learned from th e innovation vector after each analysis step, and the dynamical filter upda tes the error sub-space according to the model dynamics. The SEEK filter ha s been implemented to assimilate synthetic observations of the surface topo graphy in a non-linear, primitive equation model that uses density as verti cal coordinate. A simplified box configuration has been adopted to simulate a Gulf Stream-like current and its associated eddies and gyres with a reso lution of 20 km in the horizontal, and 4 levels in the vertical. The concep t of twin experiments is used to demonstrate that the conventional SEEK fil ter must be complemented by a learning mechanism in order to model the misr epresented tail of the error covariance matrix. An approach based on the ve rtical physics of the isopycnal model, is shown particularly robust to cont rol the velocity field in deep layers with surface observations only. The c ost of the method makes it a suitable candidate for large-size assimilation problems and operational applications. (C) 1999 Elsevier Science B.V. All rights reserved.