scholarly journals Numerical Modeling of Marine Circulation with 4D Variational Data Assimilation

2020 ◽  
Vol 8 (7) ◽  
pp. 503
Author(s):  
Vladimir Zalesny ◽  
Valeriy Agoshkov ◽  
Victor Shutyaev ◽  
Eugene Parmuzin ◽  
Natalia Zakharova

The technology is presented for modeling and prediction of marine hydrophysical fields based on the 4D variational data assimilation technique developed at the Marchuk Institute of Numerical Mathematics, Russian Academy of Sciences (INM RAS). The technology is based on solving equations of marine hydrodynamics using multicomponent splitting, thereby solving an optimality system that includes adjoint equations and covariance matrices of observation errors. The hydrodynamic model is described by primitive equations in the sigma-coordinate system, which is solved by finite-difference methods. The technology includes original algorithms for solving the problems of variational data assimilation using modern iterative processes with a special choice of iterative parameters. The methods and technology are illustrated by the example of solving the problem of circulation of the Baltic Sea with 4D variational data assimilation of sea surface temperature information.

2018 ◽  
Vol 25 (2) ◽  
pp. 429-439 ◽  
Author(s):  
Victor Shutyaev ◽  
Francois-Xavier Le Dimet ◽  
Eugene Parmuzin

Abstract. The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find unknown parameters of the model. The observation data, and hence the optimal solution, may contain uncertainties. A response function is considered as a functional of the optimal solution after assimilation. Based on the second-order adjoint techniques, the sensitivity of the response function to the observation data is studied. The gradient of the response function is related to the solution of a nonstandard problem involving the coupled system of direct and adjoint equations. The nonstandard problem is studied, based on the Hessian of the original cost function. An algorithm to compute the gradient of the response function with respect to observations is presented. A numerical example is given for the variational data assimilation problem related to sea surface temperature for the Baltic Sea thermodynamics model.


Author(s):  
François-Xavier Le Dimet ◽  
Victor Shutyaev ◽  
Eugene Parmuzin

AbstractThe problem of variational data assimilation for a model of ocean thermodynamics is formulated as an optimal control problem to find the boundary heat flux. The sensitivity of functionals of the optimal solution with respect to observations is studied. Computing the gradient of the functionals is reduced to the solution of a non-standard problem which is a coupled system involving direct and adjoint equations with mutually dependent boundary values. Solvability of the non-standard problem is studied based on the Hessian of the original cost function. An algorithm for computing the gradient of the response function related to the mean surface temperature is developed and justified. Numerical examples are presented for the Baltic Sea thermodynamics model.


2018 ◽  
Author(s):  
Victor Shutyaev ◽  
Francois-Xavier Le Dimet ◽  
Eugene Parmuzin

Abstract. The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find unknown parameters of the model. The observation data, and hence the optimal solution, may contain uncertainties. A response function is considered as a functional of the optimal solution after assimilation. Based on the second-order adjoint techniques, the sensitivity of the response function to the observation data is studied. The gradient of the response function is related to the solution of a non-standard problem involving the coupled system of direct and adjoint equations. The non-standard problem is studied, based on the Hessian of the original cost function. An algorithm to compute the gradient of the response function with respect to observations is presented. Numerical example is given for variational data assimilation problem related to sea surface temperature for the Baltic Sea thermodynamics model.


2018 ◽  
Vol 33 (2) ◽  
pp. 137-147 ◽  
Author(s):  
Victor P. Shutyaev

Abstract The problem of variational data assimilation for an evolution model is considered with the aim to identify the initial condition. The solvability of the optimality system is studied. Based on the adjoint equations, iterative algorithms for solving the problem are developed and justified.


Author(s):  
Victor Shutyaev ◽  
Francois-Xavier Le Dimet ◽  
Elena Shubina

AbstractThe problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find the initial condition function (analysis). The observation data, and hence the optimal solution, may contain uncertainties. A response function is considered as a functional of the optimal solution after assimilation. Based on the second-order adjoint techniques, the sensitivity of the response function to the observation data is studied. The gradient of the response function is related to the solution of a non-standard problem involving the coupled system of direct and adjoint equations. The solvability of the non-standard problem is studied, based on the Hessian of the original cost function. An algorithm to compute the gradient of the response function with respect to observations is developed and justified.


2019 ◽  
Vol 486 (4) ◽  
pp. 421-425
Author(s):  
V. P. Shutyaev ◽  
F.-X. Le Dimet

The problem of variational data assimilation for a nonlinear evolutionary model is formulated as an optimal control problem to find simultaneously unknown parameters and the initial state of the model. The response function is considered as a functional of the optimal solution found as a result of assimilation. The sensitivity of the functional to observational data is studied. The gradient of the functional with respect to observations is associated with the solution of a nonstandard problem involving a system of direct and adjoint equations. On the basis of the Hessian of the original cost function, the solvability of the nonstandard problem is studied. An algorithm for calculating the gradient of the response function with respect to observational data is formulated and justified.


2019 ◽  
Vol 7 (6) ◽  
pp. 195 ◽  
Author(s):  
Valery Agoshkov ◽  
Natalia Lezina ◽  
Tatiana Sheloput

One of the modern fields in mathematical modelling of water areas is developing hybrid coastal ocean models based on domain decomposition. In coastal ocean modelling a problem to be solved is setting open boundary conditions. One of the methods dealing with open boundaries is variational data assimilation. The purpose of this work is to apply the domain decomposition method to the variational data assimilation problem. The method to solve the problem of restoring boundary functions at the liquid boundaries for a system of linearized shallow water equations is studied. The problem of determining additional unknowns is considered as an inverse problem and solved using well-known approaches. The methodology based on the theory of optimal control and adjoint equations is used. In the paper the theoretical study of the problem is carried out, unique and dense solvability of the problem is proved, an iterative algorithm is proposed and its convergence is studied. The results of the numerical experiments are presented and discussed.


Sign in / Sign up

Export Citation Format

Share Document