Comparison of generalized data assimilation method with an Ensemble Optimal Interpolation scheme

Author(s):  
Konstantin Belyaev ◽  
Andrey Kuleshov ◽  
Ilya Smirnov ◽  
Clemente A.S. Tanajura
2019 ◽  
Author(s):  
Konstantin Belyaev ◽  
Andrey Kuleshov ◽  
Ilya Smirnov ◽  
Clemente A. S. Tanajura

Abstract. An original hybrid data assimilation scheme recently developed is presented and tested. The scheme is based on the application of the theory of diffusion random processes. It is applied here in conjunction with the Hybrid-Coordinate Ocean Model (HYCOM) to assimilate altimetry data from the Archiving, Validating and Interpolating Satellite Oceanography Data (AVISO) in the Atlantic. Several numerical experiments were conducted and their results were analyzed. It is shown that the method is able to assimilate data and to produce analyses closer to observations. It also conserves the model balance. This method allows calculating the confidence range of the analyses by estimating their errors The presented method is compared with the Ensemble Optimal Interpolation scheme (EnOI) and it is shown that it has several advantages, in particular, it provides a better forecast and requires less computational cost.


2011 ◽  
Vol 28 (12) ◽  
pp. 1624-1640 ◽  
Author(s):  
Weiwei Fu ◽  
Jiang Zhu

Abstract Sea level anomalies (SLA) from the Ocean Topography Experiment (TOPEX)/Poseidon are assimilated with three-dimensional variational data assimilation (3DVAR) and ensemble optimal interpolation (EnOI) for the period of 1997–2001. When sea level data are assimilated, one major concern is how to project the surface information downward. In 3DVAR, downward projection is usually achieved by minimizing a cost function that computes the relations among temperature, salinity, and sea level. In EnOI, the surface information is propagated to other variables through a stationary ensemble. Their effects on the simulated variability are evaluated in a tropical Pacific Ocean model. When compared with different datasets, it is found that effects of 3DVAR and EnOI are different in several aspects. For sea level, the standard deviation is improved by both methods, but EnOI is more effective in the central/eastern Pacific. The SLA evolution is better reproduced with EnOI than with 3DVAR. For temperature, the model–reanalysis correlations are increased by 0.1–0.2 in the top 200 m with both methods, but EnOI is more effective, especially along the thermocline depth. When compared with the Tropical Atmosphere–Ocean array (TAO) profiles, evolution of the temperature reveals that 3DVAR tends to cause more errors during ENSO events. The correlations with TAO profile are increased by 0.1–0.3 with EnOI and are generally decreased by 0.1–0.3 with 3DVAR. For salinity, both methods have weak impact on the model–reanalysis correlations above the thermocline. Relative to 3DVAR, EnOI can increase the correlation by 0.2 below the thermocline. When compared with the TAO profiles, the differences are reduced to some extent with both methods, but 3DVAR is very negative on the simulated variability.


Mathematics ◽  
2021 ◽  
Vol 9 (19) ◽  
pp. 2371
Author(s):  
Konstantin Belyaev ◽  
Andrey Kuleshov ◽  
Ilya Smirnov ◽  
Clemente A. S. Tanajura

In this paper, we consider a recently developed data assimilation method, the Generalized Kalman Filter (GKF), which is a generalization of the widely-used Ensemble Optimal Interpolation (EnOI) method. Both methods are applied for modeling the Atlantic Ocean circulation using the known Hybrid Coordinate Ocean Model. The along-track altimetry data taken from the Archiving, Validating and Interpolating Satellite Oceanography Data (AVISO) were used for data assimilation and other data from independent archives of observations; particularly, the temperature and salinity data from the Pilot Research Array in the Tropical Atlantic were used for independent comparison. Several numerical experiments were performed with their results discussed and analyzed. It is shown that values of the ocean state variables obtained in the calculations using the GKF method are closer to the observations in terms of standard metrics in comparison with the calculations using the standard data assimilation method EnOI. Furthermore, the GKF method requires less computational effort compared to the EnOI method.


2021 ◽  
Vol 51 ◽  
pp. 101317
Author(s):  
Habib Toye ◽  
Peng Zhan ◽  
Furrukh Sana ◽  
Sivareddy Sanikommu ◽  
Naila Raboudi ◽  
...  

2005 ◽  
Vol 133 (8) ◽  
pp. 2310-2334 ◽  
Author(s):  
Anna Borovikov ◽  
Michele M. Rienecker ◽  
Christian L. Keppenne ◽  
Gregory C. Johnson

Abstract One of the most difficult aspects of ocean-state estimation is the prescription of the model forecast error covariances. The paucity of ocean observations limits our ability to estimate the covariance structures from model–observation differences. In most practical applications, simple covariances are usually prescribed. Rarely are cross covariances between different model variables used. Here a comparison is made between a univariate optimal interpolation (UOI) scheme and a multivariate OI algorithm (MvOI) in the assimilation of ocean temperature profiles. In the UOI case only temperature is updated using a Gaussian covariance function. In the MvOI, salinity, zonal, and meridional velocities as well as temperature are updated using an empirically estimated multivariate covariance matrix. Earlier studies have shown that a univariate OI has a detrimental effect on the salinity and velocity fields of the model. Apparently, in a sequential framework it is important to analyze temperature and salinity together. For the MvOI an estimate of the forecast error statistics is made by Monte Carlo techniques from an ensemble of model forecasts. An important advantage of using an ensemble of ocean states is that it provides a natural way to estimate cross covariances between the fields of different physical variables constituting the model-state vector, at the same time incorporating the model’s dynamical and thermodynamical constraints as well as the effects of physical boundaries. Only temperature observations from the Tropical Atmosphere–Ocean array have been assimilated in this study. To investigate the efficacy of the multivariate scheme, two data assimilation experiments are validated with a large independent set of recently published subsurface observations of salinity, zonal velocity, and temperature. For reference, a control run with no data assimilation is used to check how the data assimilation affects systematic model errors. While the performance of the UOI and MvOI is similar with respect to the temperature field, the salinity and velocity fields are greatly improved when the multivariate correction is used, as is evident from the analyses of the rms differences between these fields and independent observations. The MvOI assimilation is found to improve upon the control run in generating water masses with properties close to the observed, while the UOI fails to maintain the temperature and salinity structure.


Sign in / Sign up

Export Citation Format

Share Document