Receiver Functions from Multiple-Taper Spectral Correlation Estimates

Teleseismic P waves are followed by a series of scattered waves, particularly P -to- S converted phases, that form a coda. The sequence of scattered waves on the horizontal components can be represented by the receiver function (RF) for the station and may vary with the approach angle and azimuth of the incoming P wave. We have developed a frequency-domain RF inversion algorithm using multiple-taper correlation (MTC) estimates, instead of spectral division, using the pre-event noise spectrum for frequency-dependent damping. The multitaper spectrum estimates are leakage resistant, so low-amplitude portions of the P -wave spectrum can contribute usefully to the RF estimate. The coherence between vertical and horizontal components can be used to obtain a frequency-dependent uncertainty for the RF. We compare the MTC method with two popular methods for RF estimation, time-domain deconvolution (TDD), and spectral division (SPD), both with damping to avoid numerical instabilities. Deconvolution operators are often biased toward the frequencies where signal is strongest. Spectral-division schemes with constant water-level damping can suffer from the same problem in the presence of strong signal-generated noise. Estimates of uncertainty are scarce for TDD and SPD, which impedes developing a weighted average of RF estimates from multiple events. Multiple-taper correlation RFs are more resistant to signal-generated noise in test cases, though a “coherent” scattering effect, like a strong near-surface organ-pipe resonance in soft sediments, will overprint the Ps conversions from deeper interfaces. The MTC RF analysis confirms the broad features of an earlier RF study for the Urals foredeep by Levin and Park (1997a) using station ARU of the Global Seismographic Network (GSN), but adds considerable detail, resolving P -to- S converted energy up to f = 4.0 Hz.

[1]  M. Bostock.  Anisotropic upper-mantle stratigraphy and architecture of the Slave craton , 1997, Nature.

[2]  A. Sheehan,et al.  Mantle discontinuity structure from midpoint stacks of converted P to S waves across the Yellowstone hotspot track , 1997 .

[3]  Vadim Levin,et al.  P-SH conversions in a flat-layered medium with anisotropy of arbitrary orientation , 1997 .

[4]  P-SV resonances in irregular low-velocity surface layers , 1985 .

[5]  D. Thomson,et al.  Spectrum estimation and harmonic analysis , 1982, Proceedings of the IEEE.

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

[7]  J. Lees,et al.  Multiple-taper spectral analysis: a stand-alone C-subroutine , 1995 .

[8]  L. J. Burdick,et al.  Modeling crustal structure through the use of converted phases in teleseismic body-wave forms , 1977, Bulletin of the Seismological Society of America.

[9]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

[10]  M. Régnier.  Lateral variation of upper mantle structure beneath New Caledonia determined from P-wave receiver function : evidence for a fossil subduction zone , 1988 .

[11]  Frank L. Vernon,et al.  Multitaper spectral analysis of high-frequency seismograms , 1987 .

[12]  David J. Thomson,et al.  Coherence established between atmospheric carbon dioxide and global temperature , 1990, Nature.

[13]  Arthur L. Lerner-Lam,et al.  Crustal thickness variations across the Colorado Rocky Mountains from teleseismic receiver functions , 1995 .

[14]  J. Minster,et al.  Constraints on crustal structure and complex Moho topography beneath Piñon Flat, California, from teleseismic receiver functions , 1996, Bulletin of the Seismological Society of America.

[15]  T. J. Owens,et al.  Shallow structure effects on broadband teleseismic P waveforms , 1988 .

[16]  W. Menke,et al.  Shear wave splitting in the Appalachians and the Urals: A case for multilayered anisotropy , 1999 .

[17]  J. Cassidy,et al.  Numerical experiments in broadband receiver function analysis , 1992, Bulletin of the Seismological Society of America.

[18]  C. Langston.  Evidence for the subducting lithosphere under southern Vancouver Island and western Oregon from teleseismic P wave conversions , 1981 .

[19]  J. Montagner,et al.  Shear wave splitting in the mantle Ps phases , 1996 .

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

[21]  E. Gordeev,et al.  Seismicity and structure of the Kamchatka subduction zone , 1997 .

[22]  G. Zandt,et al.  Velocity structure beneath College Station Alaska from receiver functions , 1996 .

[23]  G. Abers,et al.  Source scaling of earthquakes in the shumagin region, Alaska: time-domain inversions of regional waveforms , 1995 .

[24]  J. Minster,et al.  Thickness estimates of the upper-mantle transition zone from bootstrapped velocity spectrum stacks of receiver functions , 1998 .

[25]  Mauricio D. Sacchi,et al.  Deconvolution of teleseismic recordings for mantle structure , 1997 .

[26]  M. Stiller,et al.  Preserved Collisional Crustal Structure of the Southern Urals Revealed by Vibroseis Profiling , 1996, Science.

[27]  G. Abers.  Array measurements of phases used in receiver-function calculations: Importance of scattering , 1998, Bulletin of the Seismological Society of America.

[28]  Dogan Seber,et al.  Grid search modeling of receiver functions: Implications for crustal structure in the Middle East and North Africa , 1998 .

[29]  Vadim Levin,et al.  P-SH Conversions in Layered Media with Hexagonally Symmetric Anisotropy: A CookBook , 1998 .

[30]  Vadim Levin,et al.  Crustal anisotropy in the Ural Mountains Foredeep from teleseismic receiver functions , 1997 .

[31]  Jeffrey Park.  Surface waves in layered anisotropic structures , 1996 .

[32]  E. Sembera,et al.  Coherence of seismic body waves from local events as measured by a small‐aperture array , 1991 .

[33]  Sobolev,et al.  Seismic Evidence for a Detached Indian Lithospheric Mantle Beneath Tibet. , 1999, Science.

[34]  L. Vinnik,et al.  Anisotropy of the mantle inferred from observations of P to S converted waves , 1984 .

[35]  Robert W. Clayton,et al.  Source shape estimation and deconvolution of teleseismic bodywaves , 1976 .

[36]  S. Hough.  Constraining sediment thickness in the San Francisco Bay Area using observed resonances and P‐to‐S conversions , 1990 .

[37]  J. Bernard Minster,et al.  Simultaneous time-domain deconvolution with application to the computation of receiver functions , 1995 .

[38]  M. Savage.  Lower crustal anisotropy or dipping boundaries? Effects on receiver functions and a case study in New Zealand , 1998 .

[39]  Jeffrey Park,et al.  Shear zones in the Proterozoic lithosphere of the Arabian shield and the nature of the Hales discontinuity , 2000 .

[40]  J. Cassidy,et al.  A relic proterozoic subduction zone in western Canada: New evidence from seismic reflection and receiver function data , 1996 .