Sensitivity with respect to observations in variational data assimilation

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.

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.


2021 ◽  
Vol 2099 (1) ◽  
pp. 012031
Author(s):  
V P Shutyaev ◽  
E I Parmuzin ◽  
I Yu Gejadze

Abstract The sensitivity of functionals of the optimal solution to a variational data assimilation problem for the sea thermodynamics model is studied. The variational data assimilation problem is formulated as an optimal control problem to find the initial state and the boundary condition. The sensitivity of the response functions as functionals of the optimal solution with respect to the observation data is determined by the gradient of the response function and reduces to the solution of a non-standard problem being a coupled system of direct and adjoint equations with mutually dependent initial and boundary values. The algorithm to compute the gradient of the response function is presented, based on the Hessian of the original cost functional. The sensitivity analysis of the response function with respect to errors of observation data is carried out. Numerical examples are presented for the Black Sea thermodynamics model.


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.


2021 ◽  
Vol 36 (6) ◽  
pp. 347-357
Author(s):  
Victor Shutyaev ◽  
Eugene Parmuzin ◽  
Igor Gejadze

Abstract The problem of stability and sensitivity of functionals of the optimal solution of the variational data assimilation of sea surface temperature for the model of sea thermodynamics is considered. The variational data assimilation problem is formulated as an optimal control problem to find the initial state and the boundary heat flux. The sensitivity of the response functions as functionals of the optimal solution with respect to the observation data is studied. Computing the gradient of the response function reduces to the solution of a non-standard problem being a coupled system of direct and adjoint equations with mutually dependent initial and boundary values. The algorithm to compute the gradient of the response function is presented, based on the Hessian of the original cost functional. Stability analysis of the response function with respect to uncertainties of input data is given. Numerical examples are presented for the Black and Azov seas 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.


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.


Author(s):  
Victor Shutyaev ◽  
Arthur Vidard ◽  
François-Xavier Le Dimet ◽  
Igor Gejadze

AbstractThe problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find the initial condition. The optimal solution (analysis) error arises due to the errors in the input data (background and observation errors). Under the Gaussian assumption the optimal solution error covariance can be constructed using the Hessian of the auxiliary data assimilation problem. The aim of this paper is to study the evolution of model errors via data assimilation. The optimal solution error covariances are derived in the case of imperfect model and for the weak constraint formulation, when the model euations determine the cost functional.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document