Received: 19 Jan 2019 – Accepted for review: 28 Jan 2019 – Discussion started: 08 Feb 2019
Abstract. Statistical Data Assimilation (SDA) is the transfer of information from field or laboratory observations to a user selected model of the dynamical system producing those observations. The data is noisy and the model has errors; the information transfer addresses properties of the conditional probability distribution of the states of the model conditioned on the observations. The quantities of interest in SDA are the conditional expected values of functions of the model state, and these require the approximate evaluation of high dimensional integrals. We introduce a conditional probability distribution and use the Laplace method with annealing to identify the maxima of the conditional probability distribution. The annealing method slowly increases the precision term of the model as it enters the Laplace method. In this paper, we extend the idea of precision annealing (PA) to Monte Carlo calculations of conditional expected values using Metropolis-Hastings methods.
This preprint has been withdrawn.
How to cite. Wong, A. S., Hao, K., Fang, Z., and Abarbanel, H. D. I.: Precision Annealing Monte Carlo Methods for Statistical Data Assimilation: Metropolis-Hastings Procedures, Nonlin. Processes Geophys. Discuss. [preprint], https://doi.org/10.5194/npg-2019-1, 2019.
Our paper deals with data assimilation methods for chaotic systems and how one can make predictions from incomplete data of such systems. The method that we chose to explore in detail is a Monte Carlo method with an annealing heuristic. Our results show that Monte Carlo methods are a viable alternatives to the standard set of derivative-based methods. We verify the method using the Lorenz 96 system due to the simplicity of that system.
Our paper deals with data assimilation methods for chaotic systems and how one can make...