Inverting interpolated receiver functions with surface wave dispersion and gravity: Application to the western U.S. and adjacent Canada and Mexico

We use P wave receiver functions from the western U.S. and adjacent regions to construct a receiver function wavefield interpolation scheme that helps to equalize the lateral sampling of the receiver functions and the surface wave dispersion and to greatly simplify the receiver functions. Spatial interpolation and smoothing suppress poorly sampled and difficult to interpret back azimuthal variations and allow the extraction of the first-order features in the receiver function wavefield, including observations from several ray parameter ranges. We combine the interpolated receiver functions with Rayleigh wave dispersion estimates and surface gravity observations to estimate the 3-D shear wave speed beneath the region. Speed variations in the 3-D model correlate strongly with expected geologic variations and illuminate broad-scale features of the western U.S. crust and upper mantle. The model is smooth, self-consistent, and demonstrates the compatibility of the interpolated receiver functions and dispersion observations.

[1]  Charles J. Ammon,et al.  Iterative deconvolution and receiver-function estimation , 1999 .

[2]  R. Herrmann,et al.  Lithospheric structure of the Arabian Shield from the joint inversion of receiver functions and surface-wave group velocities , 2003 .

[3]  N. Biswas.  Earth-flattening procedure for the propagation of Rayleigh wave , 1972 .

[4]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[5]  Donald Plouff,et al.  GRAVITY AND MAGNETIC FIELDS OF POLYGONAL PRISMS AND APPLICATION TO MAGNETIC TERRAIN CORRECTIONS , 1976 .

[6]  A. Levander,et al.  Continuing Colorado plateau uplift by delamination-style convective lithospheric downwelling , 2011, Nature.

[7]  Charles J. Ammon,et al.  Lithospheric Structure of the Arabian Shield from the Joint Inversion of Receiver Function and Surface-Wave Dispersion Observations , 2000 .

[8]  J. Brune,et al.  Shear-wave velocity structure in the northern Basin and Range province from the combined analysis of receiver functions and surface waves , 1997, Bulletin of the Seismological Society of America.

[9]  R. Allen,et al.  Seismic imaging east of the Rocky Mountains with USArray , 2014 .

[10]  Charles A. Langston,et al.  Structure under Mount Rainier, Washington, inferred from teleseismic body waves , 1979 .

[11]  Remko Scharroo,et al.  Generic Mapping Tools: Improved Version Released , 2013 .

[12]  Walter R. Roest,et al.  Age, spreading rates, and spreading asymmetry of the world's ocean crust , 2008 .

[13]  Charles J. Ammon,et al.  The isolation of receiver effects from teleseismic P waveforms , 1991, Bulletin of the Seismological Society of America.

[14]  B. Zurek,et al.  Thick-structured Proterozoic lithosphere of the Rocky Mountain region , 2001 .

[15]  Robert B. Herrmann,et al.  Computer Programs in Seismology: An Evolving Tool for Instruction and Research , 2013 .

[16]  Carl Tape,et al.  Adjoint Tomography of the Southern California Crust , 2009, Science.

[17]  Gary L. Pavlis,et al.  Model Update March 2011: Upper Mantle Heterogeneity beneath North America from Traveltime Tomography with Global and USArray Transportable Array Data , 2012 .

[18]  T. Zheng,et al.  Receiver Function Imaging with Reconstructed Wavefields from Sparsely Scattered Stations , 2015 .

[19]  K. Fischer,et al.  The lithosphere–asthenosphere boundary and the tectonic and magmatic history of the northwestern United States , 2014 .

[20]  Heather Bedle,et al.  S velocity variations beneath North America , 2009 .

[21]  Michael H. Ritzwoller,et al.  Joint inversion of surface wave dispersion and receiver functions: a Bayesian Monte-Carlo approach , 2013 .

[22]  R. Allen,et al.  Lithosphere-asthenosphere interaction beneath the western United States from the joint inversion of body-wave traveltimes and surface-wave phase velocities , 2011 .

[23]  C. Amante,et al.  ETOPO1 arc-minute global relief model : procedures, data sources and analysis , 2009 .

[24]  K. Fischer,et al.  Contrasting lithospheric signatures across the western United States revealed by Sp receiver functions , 2014 .

[25]  E. R. Engdahl,et al.  Constraints on seismic velocities in the Earth from traveltimes , 1995 .

[26]  Jeroen Tromp,et al.  Imaging lateral heterogeneity in the northern Apennines from time reversal of reflected surface waves , 2007 .

[27]  John E. Nafe,et al.  PHYSICAL PROPERTIES OF MARINE SEDIMENTS , 1961 .

[28]  M. Ritzwoller,et al.  A 3‐D model of the crust and uppermost mantle beneath the Central and Western US by joint inversion of receiver functions and surface wave dispersion , 2013 .

[29]  C. Ammon,et al.  Evidence for mafic lower crust in Tanzania, East Africa, from joint inversion of receiver functions and Rayleigh wave dispersion velocities , 2005 .

[30]  George E. Randall,et al.  Efficient calculation of complete differential seismograms for laterally homogeneous earth models , 1994 .

[31]  G. Masters,et al.  Update on CRUST1.0 - A 1-degree Global Model of Earth's Crust , 2013 .

[32]  F. Birch.  The velocity of compressional waves in rocks to 10 kilobars: 1. , 1960 .

[33]  G. Ekström.  A global model of Love and Rayleigh surface wave dispersion and anisotropy, 25–250 s , 2011 .

[34]  Lion Krischer,et al.  ObsPy: A Python Toolbox for Seismology , 2010 .

[35]  B. Romanowicz,et al.  Tectonic regionalization without a priori information: A cluster analysis of upper mantle tomography , 2011 .

[36]  J. Pederson,et al.  Colorado Plateau magmatism and uplift by warming of heterogeneous lithosphere , 2009, Nature.

[37]  Monica Maceira,et al.  Joint inversion of surface wave velocity and gravity observations and its application to central Asian basins shear velocity structure , 2009 .

[38]  Georges Balmino,et al.  Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies , 2012, Journal of Geodesy.

[39]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[40]  Gary L. Pavlis,et al.  Three‐dimensional, prestack, plane wave migration of teleseismic P‐to‐S converted phases: 1. Theory , 2003 .

[41]  G. Randall,et al.  On the nonuniqueness of receiver function inversions , 1990 .