English

A probabilistic scheme for joint parameter estimation and state prediction in complex dynamical systems

Computation 2018-12-12 v2 Probability Chaotic Dynamics

Abstract

Many problems in the geophysical sciences demand the ability to calibrate the parameters and predict the time evolution of complex dynamical models using sequentially-collected data. Here we introduce a general methodology for the joint estimation of the static parameters and the forecasting of the state variables of nonlinear, and possibly chaotic, dynamical models. The proposed scheme is essentially probabilistic. It aims at recursively computing the sequence of joint posterior probability distributions of the unknown model parameters and its (time varying) state variables conditional on the available observations. The latter are possibly partial and contaminated by noise. The new framework combines a Monte Carlo scheme to approximate the posterior distribution of the fixed parameters with filtering (or {\em data assimilation}) techniques to track and predict the distribution of the state variables. For this reason, we refer to the proposed methodology as {\em nested filtering}. In this paper we specifically explore the use of Gaussian filtering methods, but other approaches fit naturally within the new framework. As an illustrative example, we apply three different implementations of the methodology to the tracking of the state, and the estimation of the fixed parameters, of a stochastic two-scale Lorenz 96 system. This model is commonly used to assess data assimilation procedures in meteorology. For this example, we compare different nested filters and show estimation and forecasting results for a 4,000-dimensional system.

Keywords

Cite

@article{arxiv.1708.03730,
  title  = {A probabilistic scheme for joint parameter estimation and state prediction in complex dynamical systems},
  author = {Sara Pérez-Vieites and Inés P. Mariño and Joaquín Míguez},
  journal= {arXiv preprint arXiv:1708.03730},
  year   = {2018}
}

Comments

In this version we have tried to clarify the presentation of the methodology and the analysis. The body of the manuscript has been simplified (in terms of the notation as well) by moving some algorithms and the proofs to appendices. We have added a new variant of the algorithm based on a sequential quasi-Monte Carlo scheme. All numerical results have been updated

R2 v1 2026-06-22T21:13:00.673Z