Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union

Journal metrics

  • IF value: 1.329 IF 1.329
  • IF 5-year<br/> value: 1.394 IF 5-year
  • CiteScore<br/> value: 1.27 CiteScore
  • SNIP value: 0.903 SNIP 0.903
  • SJR value: 0.709 SJR 0.709
  • IPP value: 1.455 IPP 1.455
  • h5-index value: 20 h5-index 20
Nonlin. Processes Geophys., 21, 869-885, 2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
Research article
25 Aug 2014
A non-Gaussian analysis scheme using rank histograms for ensemble data assimilation
S. Metref1, E. Cosme1, C. Snyder2, and P. Brasseur1 1CNRS – Univ. Grenoble Alpes, LGGE (UMR5183), 38041 Grenoble, France
2National Center for Atmospheric Research, Boulder, Colorado, USA
Abstract. One challenge of geophysical data assimilation is to address the issue of non-Gaussianities in the distributions of the physical variables ensuing, in many cases, from nonlinear dynamical models. Non-Gaussian ensemble analysis methods fall into two categories, those remapping the ensemble particles by approximating the best linear unbiased estimate, for example, the ensemble Kalman filter (EnKF), and those resampling the particles by directly applying Bayes' rule, like particle filters. In this article, it is suggested that the most common remapping methods can only handle weakly non-Gaussian distributions, while the others suffer from sampling issues. In between those two categories, a new remapping method directly applying Bayes' rule, the multivariate rank histogram filter (MRHF), is introduced as an extension of the rank histogram filter (RHF) first introduced by Anderson (2010). Its performance is evaluated and compared with several data assimilation methods, on different levels of non-Gaussianity with the Lorenz 63 model. The method's behavior is then illustrated on a simple density estimation problem using ensemble simulations from a coupled physical–biogeochemical model of the North Atlantic ocean. The MRHF performs well with low-dimensional systems in strongly non-Gaussian regimes.

Citation: Metref, S., Cosme, E., Snyder, C., and Brasseur, P.: A non-Gaussian analysis scheme using rank histograms for ensemble data assimilation, Nonlin. Processes Geophys., 21, 869-885,, 2014.
Publications Copernicus