Joint inversion of receiver function and surface wave dispersion observations |
| |
| Authors: | J Julià C J Ammon R B Herrmann A M Correig |
| |
| Institution: | ;Departament d'Astronomia i Meteorologia, Facultat de F´sica, Universitat de Barcelona, Mart´i Franqués, 1, E08028 Barcelona, Spain. Email: ;Department of Earth and Atmospheric Sciences, Saint Louis University, 3507 Laclede Avenue, St. Louis, MO 63103, USA. Emails:, , ;Laboratori d'Estudis Geof´sics 'Eduard Fontserè', Institut d'Estudis Catalans, Carme, 47, E08001 Barcelona, Spain |
| |
| Abstract: | We implement a method to invert jointly teleseismic P wave receiver functions and surface wave group and phase velocities for a mutually consistent estimate of earth structure. Receiver functions are primarily sensitive to shear wave velocity contrasts and vertical traveltimes, and surface wave dispersion measurements are sensitive to vertical shear wave velocity averages. Their combination may bridge resolution gaps associated with each individual data set. We formulate a linearized shear velocity inversion that is solved using a damped leastsquares scheme that incorporates a priori smoothness constraints for velocities in adjacent layers. The data sets are equalized for the number of data points and physical units in the inversion process. The combination of information produces a relatively simple model with a minimal number of sharp velocity contrasts. We illustrate the approach using noisefree and realistic noise simulations and conclude with an inversion of observations from the Saudi Arabian Shield. Inversion results for station SODA, located in the Arabian Shield, include a crust with a sharp gradient near the surface (shear velocity changing from 1.8 to 3.5 km s1 in 3 km) underlain by a 5kmthick layer with a shear velocity of 3.5 km s1 and a 27kmthick layer with a shear velocity of 3.8 km s1, and an upper mantle with an average shear velocity of 4.7 km s1. The crustmantle transition has a significant gradient, with velocity values varying from 3.8 to 4.7 km s1 between 35 and 40 km depth. Our results are compatible with independent inversions for crustal structure using refraction data. |
| |
| Keywords: | body waves crustal structure inverse problem inversion layered media surface waves |
|
|
