Journal cover Journal topic
Solid Earth An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 2.921 IF 2.921
  • IF 5-year value: 3.087 IF 5-year
  • CiteScore value: 4.8 CiteScore
  • SNIP value: 1.314 SNIP 1.314
  • IPP value: 2.87 IPP 2.87
  • SJR value: 0.993 SJR 0.993
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 38 Scimago H
    index 38
  • h5-index value: 36 h5-index 36
Volume 4, issue 2
Solid Earth, 4, 357–371, 2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.
Solid Earth, 4, 357–371, 2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.

Method article 15 Oct 2013

Method article | 15 Oct 2013

An objective rationale for the choice of regularisation parameter with application to global multiple-frequency S-wave tomography

C. Zaroli1, M. Sambridge2, J.-J. Lévêque1, E. Debayle3, and G. Nolet4 C. Zaroli et al.
  • 1Institut de Physique du Globe de Strasbourg, UMR7516, Université de Strasbourg, EOST/CNRS, France
  • 2Research School of Earth Sciences, Australian National University, Canberra ACT 0200, Australia
  • 3Laboratoire de Géologie de Lyon, UMR5276, CNRS, Université Lyon 1 et ENS de Lyon, France
  • 4Géoazur, Université de Nice Sophia-Antipolis, UMR6526, CNRS, France

Abstract. In a linear ill-posed inverse problem, the regularisation parameter (damping) controls the balance between minimising both the residual data misfit and the model norm. Poor knowledge of data uncertainties often makes the selection of damping rather arbitrary. To go beyond that subjectivity, an objective rationale for the choice of damping is presented, which is based on the coherency of delay-time estimates in different frequency bands. Our method is tailored to the problem of global multiple-frequency tomography (MFT), using a data set of 287 078 S-wave delay times measured in five frequency bands (10, 15, 22, 34, and 51 s central periods). Whereas for each ray path the delay-time estimates should vary coherently from one period to the other, the noise most likely is not coherent. Thus, the lack of coherency of the information in different frequency bands is exploited, using an analogy with the cross-validation method, to identify models dominated by noise. In addition, a sharp change of behaviour of the model ℓ-norm, as the damping becomes lower than a threshold value, is interpreted as the signature of data noise starting to significantly pollute at least one model component. Models with damping larger than this threshold are diagnosed as being constructed with poor data exploitation. Finally, a preferred model is selected from the remaining range of permitted model solutions. This choice is quasi-objective in terms of model interpretation, as the selected model shows a high degree of similarity with almost all other permitted models (correlation superior to 98% up to spherical harmonic degree 80). The obtained tomographic model is displayed in the mid lower-mantle (660–1910 km depth), and is shown to be compatible with three other recent global shear-velocity models. A wider application of the presented rationale should permit us to converge towards more objective seismic imaging of Earth's mantle.

Publications Copernicus