English

Bayesian eikonal tomography using Gaussian processes

Geophysics 2023-02-21 v2

Abstract

Eikonal tomography has become a popular methodology for deriving phase velocity maps from surface wave phase delay measurements. Its high efficiency makes it popular for handling datasets deriving from large-N arrays, in particular in the ambient-noise tomography setting. However, the results of eikonal tomography are crucially dependent on the way in which phase delay measurements are interpolated, a point which has not been thoroughly investigated. In this work, I provide a rigorous formulation for eikonal tomography using Gaussian processes (GPs) to interpolate phase delay measurements, including uncertainties. GPs allow the posterior phase delay gradient to be analytically derived. From the phase delay gradient, an excellent approximate solution for phase velocities can be obtained using the saddlepoint method. The result is a fully Bayesian result for phase velocities of surface waves, incorporating the nonlinear wavefront bending inherent in eikonal tomography, with no sampling required. The results of this analysis imply that the uncertainties reported for eikonal tomography are often underestimated.

Keywords

Cite

@article{arxiv.2302.05267,
  title  = {Bayesian eikonal tomography using Gaussian processes},
  author = {Jack B. Muir},
  journal= {arXiv preprint arXiv:2302.05267},
  year   = {2023}
}

Comments

Updated figures, text and acknowledgements

R2 v1 2026-06-28T08:37:04.756Z