
Multimodal Layered Transdimensional Inversion of Seismic Dispersion Curves With Depth Constraints
Author(s) -
Killingbeck S. F.,
Livermore P. W.,
Booth A. D.,
West L. J.
Publication year - 2018
Publication title -
geochemistry, geophysics, geosystems
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.928
H-Index - 136
ISSN - 1525-2027
DOI - 10.1029/2018gc008000
Subject(s) - geology , inversion (geology) , ground penetrating radar , markov chain monte carlo , bayesian probability , radar , seismology , algorithm , mathematics , statistics , computer science , telecommunications , tectonics
MuLTI (Multimodal Layered Transdimensional Inversion) is a Markov chain Monte Carlo implementation of Bayesian inversion for the probability distribution of shear wave velocity (Vs) as a function of depth. Based on Multichannel Analysis of Surface Wave methods, it requires as data (i) a Rayleigh‐wave dispersion curve and (ii) additional layer depth constraints, the latter we show significantly improve resolution compared to conventional unconstrained inversions. Such depth constraints may be drawn from any source (e.g., boreholes, complementary geophysical data) provided they also represent a seismic interface. We apply MuLTI to a Norwegian glacier, Midtdalsbreen, in which ground‐penetrating radar was used to constrain internal layers of snow, ice, and subglacial sediments, with transitions at 2 and 25.5 m, and whose Vs is assumed to be in the range 500–1,700, 1,700–1,950, and 200–2,800 m/s, respectively. Synthetic modeling demonstrates that MuLTI recovers the true model of Vs variation with depth. Significantly, compared to inversions without depth constraints, in this synthetic case MuLTI not only reduces by about a factor of 10 the error between the true and the best fitting model, but also reduces by a factor of 2 the vertically averaged spread of the distribution of Vs based on the 95% credible intervals. We further show that using frequencies above about 100 Hz lead to unconverged solutions due to mode ambiguities associated with fine spatial structures. For our acquired data on Midtdalsbreen, we use 14‐100 Hz data for which MuLTI produces a large‐scale converged inversion.