Filter accuracy for the Lorenz 96 model: Fixed versus adaptive observation operators

In the context of filtering chaotic dynamical systems it is well-known that partial observations, if sufficiently informative, can be used to control the inherent uncertainty due to chaos. The purpose of this paper is to investigate, both theoretically and numerically, conditions on the observations of chaotic systems under which they can be accurately filtered. In particular, we highlight the advantage of adaptive observation operators over fixed ones. The Lorenz ’96 model is used to exemplify our findings. We consider discrete-time and continuous-time observations in our theoretical developments. We prove that, for fixed observation operator, the 3DVAR filter can recover the system state within a neighbourhood determined by the size of the observational noise. It is required that a sufficiently large proportion of the state vector is observed, and an explicit form for such sufficient fixed observation operator is given. Numerical experiments, where the data is incorporated by use of the 3DVAR and extended Kalman filters, suggest that less informative fixed operators than given by our theory can still lead to accurate signal reconstruction. Adaptive observation operators are then studied numerically; we show that, for carefully chosen adaptive observation operators, the proportion of the state vector that needs to be observed is drastically smaller than with a fixed observation operator. Indeed, we show that the number of state coordinates that need to be observed may even be significantly smaller than the total number of positive Lyapunov exponents of the underlying system.

[1]  Richard E. Mortensen,et al.  Infinite-Dimensional Dynamical Systems in Mechanics and Physics (Roger Temam) , 1991, SIAM Rev..

[2]  G. Benettin,et al.  Kolmogorov Entropy and Numerical Experiments , 1976 .

[3]  D. Menemenlis Inverse Modeling of the Ocean and Atmosphere , 2002 .

[4]  Mark Kostuk,et al.  Synchronization and statistical methods for the data assimilation of HVc neuron models , 2012 .

[5]  Henry D. I. Abarbanel,et al.  Predicting the Future: Completing Models of Observed Complex Systems , 2013 .

[6]  Edriss S. Titi,et al.  Continuous Data Assimilation Using General Interpolant Observables , 2013, J. Nonlinear Sci..

[7]  A. M. Stuart,et al.  Accuracy and stability of the continuous-time 3DVAR filter for the Navier–Stokes equation , 2012, 1210.1594.

[8]  Edriss S. Titi,et al.  Determining Modes for Continuous Data Assimilation in 2D Turbulence , 2003 .

[9]  D. S. McCormick,et al.  Accuracy and stability of filters for dissipative PDEs , 2012, 1203.5845.

[10]  T. Tarn,et al.  Observers for nonlinear stochastic systems , 1975, 1975 IEEE Conference on Decision and Control including the 14th Symposium on Adaptive Processes.

[11]  Edriss S. Titi,et al.  Discrete data assimilation in the Lorenz and 2D Navier–Stokes equations , 2010, 1010.6105.

[12]  Ning Liu,et al.  Inverse Theory for Petroleum Reservoir Characterization and History Matching , 2008 .

[13]  Anna Trevisan,et al.  Assimilation of Standard and Targeted Observations within the Unstable Subspace of the Observation–Analysis–Forecast Cycle System , 2004 .

[14]  Andrew J. Majda,et al.  Filtering Complex Turbulent Systems , 2012 .

[15]  Istvan Szunyogh,et al.  A Local Ensemble Kalman Filter for Atmospheric Data Assimilation , 2002 .

[16]  Eugenia Kalnay,et al.  Atmospheric Modeling, Data Assimilation and Predictability , 2002 .

[17]  A. Stuart,et al.  Data Assimilation: A Mathematical Introduction , 2015, 1506.07825.

[18]  Luigi Palatella,et al.  Nonlinear Processes in Geophysics On the Kalman Filter error covariance collapse into the unstable subspace , 2011 .

[19]  K. Emanuel,et al.  Optimal Sites for Supplementary Weather Observations: Simulation with a Small Model , 1998 .

[20]  A. Stuart,et al.  Analysis of the 3DVAR filter for the partially observed Lorenz'63 model , 2012, 1212.4923.