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

Journal metrics

Journal metrics

  • IF value: 1.129 IF 1.129
  • IF 5-year value: 1.519 IF 5-year 1.519
  • CiteScore value: 1.54 CiteScore 1.54
  • SNIP value: 0.798 SNIP 0.798
  • SJR value: 0.610 SJR 0.610
  • IPP value: 1.41 IPP 1.41
  • h5-index value: 21 h5-index 21
  • Scimago H index value: 48 Scimago H index 48
Volume 25, issue 1 | Copyright
Nonlin. Processes Geophys., 25, 99-123, 2018
https://doi.org/10.5194/npg-25-99-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 3.0 License.

Research article 16 Feb 2018

Research article | 16 Feb 2018

Ensemble Kalman filter for the reconstruction of the Earth's mantle circulation

Marie Bocher1,2, Alexandre Fournier3, and Nicolas Coltice2 Marie Bocher et al.
  • 1Computational Seismology, Institute of Geophysics, Department of Earth Sciences, ETH Zürich, Sonneggstrasse 5, 8092 Zürich, Switzerland
  • 2Laboratoire de Géologie de Lyon, Université Claude Bernard Lyon 1, Ecole Normale Supṕrieure de Lyon, Lyon, France
  • 3Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, CNRS (UMR 7154), Paris, France

Abstract. Recent advances in mantle convection modeling led to the release of a new generation of convection codes, able to self-consistently generate plate-like tectonics at their surface. Those models physically link mantle dynamics to surface tectonics. Combined with plate tectonic reconstructions, they have the potential to produce a new generation of mantle circulation models that use data assimilation methods and where uncertainties in plate tectonic reconstructions are taken into account. We provided a proof of this concept by applying a suboptimal Kalman filter to the reconstruction of mantle circulation (Bocher et al., 2016). Here, we propose to go one step further and apply the ensemble Kalman filter (EnKF) to this problem. The EnKF is a sequential Monte Carlo method particularly adapted to solve high-dimensional data assimilation problems with nonlinear dynamics. We tested the EnKF using synthetic observations consisting of surface velocity and heat flow measurements on a 2-D-spherical annulus model and compared it with the method developed previously. The EnKF performs on average better and is more stable than the former method. Less than 300 ensemble members are sufficient to reconstruct an evolution. We use covariance adaptive inflation and localization to correct for sampling errors. We show that the EnKF results are robust over a wide range of covariance localization parameters. The reconstruction is associated with an estimation of the error, and provides valuable information on where the reconstruction is to be trusted or not.

Publications Copernicus
Download
Short summary
We propose a new method to reconstruct the circulation in the Earth's mantle for the last 300 Myr. This method is based on the sequential assimilation of plate layouts obtained from plate tectonic reconstructions into mantle convection models. This method allows us to take into account uncertainties in plate tectonic reconstructions, and provides an estimation of the uncertainties in the final result. We test and validate this method in a controlled environment by using synthetic experiments.
We propose a new method to reconstruct the circulation in the Earth's mantle for the last...
Citation
Share