Bayesian Inversion of Receiver Functions and Surface Wave Dispersion Data in the Brazilian Northeast

2017 ◽  
Author(s):  
Thuany Patrícia Costa de Lima ◽  
Hrvoje Tkalčić ◽  
Seongryong Kim ◽  
Jordi Julià
2019 ◽  
Vol 24 (1) ◽  
pp. 101-120
Author(s):  
Kajetan Chrapkiewicz ◽  
Monika Wilde-Piórko ◽  
Marcin Polkowski ◽  
Marek Grad

AbstractNon-linear inverse problems arising in seismology are usually addressed either by linearization or by Monte Carlo methods. Neither approach is flawless. The former needs an accurate starting model; the latter is computationally intensive. Both require careful tuning of inversion parameters. An additional challenge is posed by joint inversion of data of different sensitivities and noise levels such as receiver functions and surface wave dispersion curves. We propose a generic workflow that combines advantages of both methods by endowing the linearized approach with an ensemble of homogeneous starting models. It successfully addresses several fundamental issues inherent in a wide range of inverse problems, such as trapping by local minima, exploitation of a priori knowledge, choice of a model depth, proper weighting of data sets characterized by different uncertainties, and credibility of final models. Some of them are tackled with the aid of novel 1D checkerboard tests—an intuitive and feasible addition to the resolution matrix. We applied our workflow to study the south-western margin of the East European Craton. Rayleigh wave phase velocity dispersion and P-wave receiver function data were gathered in the passive seismic experiment “13 BB Star” (2013–2016) in the area of the crust recognized by previous borehole and refraction surveys. Final models of S-wave velocity down to 300 km depth beneath the array are characterized by proximity in the parameter space and very good data fit. The maximum value in the mantle is higher by 0.1–0.2 km/s than reported for other cratons.


1992 ◽  
Vol 82 (2) ◽  
pp. 962-979
Author(s):  
Paul C. Yao ◽  
James Dorman

Abstract Group velocity dispersion of explosion-generated seismic surface waves with periods ranging from 0.2 to 1.5 sec is used to investigate shallow crustal structure of eastern and central Tennessee. Several modes of both Rayleigh and Love waves can be identified and separated on the seismograms of seven SARSN regional network stations by zero-phase digital filtering. Dispersion data for sinusoidal wave motion were based on digitized zero-crossing times. By forward modeling, we find that a wave guide of at least two layers over a half-space can adequately represent our particular multi-mode, narrow-band observations. In a layered section about 3 km thick, lower velocities characterize outcropping clastic rocks of the Cumberland plateau, and higher velocities correspond to shallow carbonate rocks of the Nashville Dome. Half-space shear velocities of about 3.9 km/sec appear to represent lower Paleozoic carbonate lithology deeper than 2 to 4 km on most of the observed paths. Our best data, interpreted jointly with earlier data of Oliver and Ewing (1958) and of Chen et al. (1989), have a composite period range of 0.2 to 40 sec, but they represent different Appalachian paths. Group velocities over this broad spectrum are satisfied by a complex model with two low-velocity layers. The uniqueness of this model cannot be demonstrated, but it represents important hypotheses concerning regional geologic features that can be tested more rigorously by improved surface-wave dispersion data.


2020 ◽  
Vol 221 (2) ◽  
pp. 938-950
Author(s):  
Pingping Wu ◽  
Handong Tan ◽  
Changhong Lin ◽  
Miao Peng ◽  
Huan Ma ◽  
...  

SUMMARY Multiphysics imaging for data inversion is of growing importance in many branches of science and engineering. Cross-gradient constraint has been considered as a feasible way to reduce the non-uniqueness problem inherent in inversion process by finding geometrically consistent images from multigeophysical data. Based on OCCAM inversion algorithm, a direct inversion method of 2-D profile velocity structure with surface wave dispersion data is proposed. Then we jointly invert the profiles of magnetotelluric and surface wave dispersion data with cross-gradient constraints. Three synthetic models, including block homogeneous or heterogeneous models with consistent or inconsistent discontinuities in velocity and resistivity, are presented to gauge the performance of the joint inversion scheme. We find that owning to the complementary advantages of the two geophysical data sets, the models recovered with structure coupling constraints exhibit higher resolution in the classification of complex geologic units and settle some imaging problems caused by the separate inversion methods. Finally, a realistic velocity model from the NE Tibetan Plateau and its corresponding resistivity model calculated by empirical law are used to test the effectiveness of the joint inversion scheme in the real geological environment.


2017 ◽  
Vol 211 (1) ◽  
pp. 528-540 ◽  
Author(s):  
Jeremy M. Gosselin ◽  
Stan E. Dosso ◽  
John F. Cassidy ◽  
Jorge E. Quijano ◽  
Sheri Molnar ◽  
...  

2014 ◽  
Vol 119 (2) ◽  
pp. 1079-1093 ◽  
Author(s):  
G. Burgos ◽  
J.-P. Montagner ◽  
E. Beucler ◽  
Y. Capdeville ◽  
A. Mocquet ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document