3D Bayesian Variational Surface Wave Tomography and Application to the Southwest China

TL;DR

This paper extends variational Bayesian methods to 3D surface wave tomography, demonstrating efficient and detailed velocity structure imaging with uncertainty quantification using synthetic and real data from Southwest China.

Contribution

It introduces three variational inference techniques for 3D surface wave tomography, improving computational efficiency and uncertainty estimation over traditional methods.

Findings

All methods yield accurate velocity estimates.

sSVGD provides more reliable uncertainty estimates.

Variational methods produce more detailed structures than traditional approaches.

Abstract

Seismic surface wave tomography uses surface wave information to obtain velocity structures in the subsurface. Due to data noise and nonlinearity of the problem, surface wave tomography often has non-unique solutions. It is therefore required to quantify uncertainty of the results in order to better interpret the resulting images. Bayesian inference is the most widely-used method for this purpose. However, the commonly-used Monte Carlo methods require huge computational cost and remains intractable in high-dimensional problems. Variational inference uses optimization to solve Bayesian inverse problems, and therefore can be more efficient in the case of large datasets and high-dimensional parameter spaces. Variational inference has been applied to 2-D surface wave tomographic problems. In this study, we extend the method to 3-D surface wave tomography by directly inverting for 3-D spatial structures from frequency-dependent travel time measurements. Specifically, we apply three variational methods, mean-field automatic differential variational inference (mean-field ADVI), physically structured variational inference (PSVI) and stochastic Stein varational gradient descent (sSVGD) to surface wave tomographic problems using both synthetic data and real data in the Southwest China. The results show that all methods can provide accurate velocity estimates, while sSVGD produces more reasonable uncertainty estimates than mean-field ADVI and PSVI because of Gaussian assumption used in the later methods. In the real data case, the variational methods provide more detailed velocity structures than those obtained using traditional methods, along with reliable uncertainty estimates. We therefore conclude that variational surface wave tomograph can be applied fruitfully to many realistic problems.