Role of the Observability Gramian in Parameter Estimation: Application to Nonchaotic and Chaotic Systems via the Forward Sensitivity Method

Data assimilation in chaotic regimes is challenging, and among the challenging aspects is placement of observations to induce convexity of the cost function in the space of control. This problem is examined by using Saltzman’s spectral model of convection that admits both chaotic and nonchaotic regi...

Full description

Bibliographic Details
Main Authors: John M. Lewis, Sivaramakrishnan Lakshmivarahan
Format: Article
Language:English
Published: MDPI AG 2022-10-01
Series:Atmosphere
Subjects:
Online Access:https://www.mdpi.com/2073-4433/13/10/1647
Description
Summary:Data assimilation in chaotic regimes is challenging, and among the challenging aspects is placement of observations to induce convexity of the cost function in the space of control. This problem is examined by using Saltzman’s spectral model of convection that admits both chaotic and nonchaotic regimes and is controlled by two parameters—Rayleigh and Prandtl numbers. The problem is simplified by stripping the seven-variable constraint to a three-variable constraint. Since emphasis is placed on observation positioning to avoid cost-function flatness, forecast sensitivity to controls is needed. Four-dimensional variational assimilation (4D-Var) is silent on this issue of observation placement while Forecast Sensitivity Method (FSM) delivers sensitivities used in placement. With knowledge of the temporal forecast sensitivity matrix <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mi>V</mi></semantics></math></inline-formula>, derivatives of the forecast variables to controls, the cost function can be expressed as a function of the observability Gramian <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msup><mi>V</mi><mi>T</mi></msup><mi>V</mi></mrow></semantics></math></inline-formula> using first-order Taylor series expansion. The goal is to locate observations at places that force the Gramian positive definite. Further, locations are chosen such that the condition number of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><mrow><msup><mi>V</mi><mi>T</mi></msup><mi>V</mi></mrow></semantics></math></inline-formula> is small and this guarantees convexity in the vicinity of the cost function minimum. Four numerical experiments are executed, and results are compared with the structure of the cost function independently determined though arduous computation over a wide range of the two nondimensional numbers. The results are especially good based on reduction in cost function value and comparison with cost function structure.
ISSN:2073-4433