An optimal data assimilation method and its application to the numerical simulation of the ocean dynamics

• Published in: Mathematical and computer modelling of dynamical systems Vol. 24; no. 1; pp. 12 - 25
• Main Authors: Belyaev, Konstantin, Kuleshov, Andrey, Tuchkova, Natalia, Tanajura, Clemente A.S.
• Format: Journal Article
• Language: English
• Published: Abingdon Taylor & Francis 02.01.2018; Taylor & Francis Ltd
• Subjects: Archiving; Computational fluid dynamics; Computer simulation; Dynamic structural analysis; Gas dynamics; Kalman filters; Mathematical analysis; Mathematical models; Numerical modelling; ocean circulation; Ocean dynamics; Oceanography; Sea level; stochastic filtering
• ISSN: 1387-3954; 1744-5051
• DOI: 10.1080/13873954.2017.1338300

An original data assimilation (DA) scheme with a general dynamics model is considered. It is shown that this scheme can be approximated by the stochastic diffusion process. The sufficient conditions to provide this approximation are formulated. Based on this algorithm a new DA method is developed. The method combines variational and statistical approaches commonly used in DA theory and minimizes the variance of the trajectory of a diffusion process in conjunction with a dynamics numerical model. In this sense the method is optimal in contrast to other DA approaches. The proposed scheme takes the model dynamics into account and in this way it differs from the well-known Kalman filter. Furthermore, the derived DA method can be applied to a very wide field of dynamical systems, for example, gas dynamics, fluid dynamics and other disciplines. However, the current study deals with oceanography and DA in oceanography specifically. Then the method is applied to the HYbrid Coordinate Ocean Model and assimilates satellite sea level anomaly data from the Archiving, Validating and Interpolating Satellite Oceanography Data over the Atlantic Ocean to correct the model state. Several numerical experiments have been performed. The experiments show that the method substantially changes the synoptic and mesoscale structure of ocean dynamics. Also, the distribution of the obtained result is estimated through the solution of the Fokker-Planck-Kolmogorov equation.