An Objective Rationale for the Choice of Regularisation Parameter with Application to Global Multiple-Frequency S-Wave Tomography

Thursday, 18 December 2014
Christophe Zaroli1, Malcolm Sambridge2, Jean-Jacques Leveque1, Eric Debayle3 and Guust Nolet4, (1)Institut de Physique du Globe Strasbourg, Strasbourg Cedex, France, (2)Australian National University, Canberra, ACT, Australia, (3)LGLTPE Laboratoire de Géologie de Lyon : Terre, Planètes et Environnement, Villeurbanne Cedex, France, (4)GeoAzur, Valbonne, France
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, using a data set of 287078 S-wave delay-times measured in five frequency bands (10, 15, 22, 34, 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 infinity-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.

The obtained tomographic model is displayed in mid lower-mantle (660-1910 km depth), and is shown to be mostly compatible with three other recent global shear-velocity models, while significant differences can be noticed. A wider application of the presented rationale should permit us to converge towards more objective seismic imaging of the Earth's mantle, using as much as possible of the relevant structural information in the data.

This work was recently published: Zaroli, C., Sambridge, M., Lévêque, J.-J., Debayle, E., and Nolet, G. (2013) - Solid Earth, 4, 357-371, doi:10.5194/se-4-357-2013