Downhole Receiver Function: a Case Study

Receiver function is defined as the spectral ratio of the radial component and the vertical component of the ground motion. It is used to characterize converted waves. We extend the use of the receiver function to downhole data using waves recorded in a borehole, excited by an earthquake of magnitude 4.0 near San Fran- cisco, California, on 26 June 1994. The focal depth of the event was 6.6 km and the epicenter was located at a distance of 12.6 km from the borehole array. Six three- component sensors were located at different depths in a borehole. To extract a co- herent response of the near-surface from the incoherent earthquake waves, we de- convolve the waves recorded by the sensors at different depths with the waves recorded by the sensor on the surface. Deconvolution applied to the waves in the S- time window recorded by the radial component result in an upgoing and adowngoing wave propagating with S-wave velocity. For the waves in the P-time window re- corded by the radial component, deconvolution also gives an upgoing and a down- going wave propagating with S-wave velocity. This interesting result suggests a P- to-S conversion at a depth below the deepest sensor. To diagnose this we compute the receiverfunction forthe borehole recording of the earthquakewaves.Thereceiver function shows an upgoing wave with an arrival close to time t ! 0 for the deepest sensor. The agreement of the upgoing wave in the receiver function with the travel- time curve for the P-to-S converted wave, calculated using the P- and the S-wave velocity profile, supports the hypothesis of a pronounced P-to-S conversion. We present a synthetic example to illustrate that the first arrival of the receiver function applied to borehole data gives the upward-propagating P-to-S converted wave. To corroborate the observation of the mode conversion, we apply receiver function to a different earthquake data recorded by the same borehole array in 1998. The focal depth of the event was 6.9 km and the epicenter was located at a distance of 13 km from the borehole array. The receiver function for these data also show an upgoing wave with a pulse close to time t ! 0 at the deepest sensor. The moveout of the upgoing wave agrees with the travel-time curve for the P-to-S converted wave, hence supporting our observation of the mode conversion.

[1]  Charles A. Langston,et al.  Corvallis, Oregon, crustal and upper mantle receiver structure from teleseismic P and S waves , 1977, Bulletin of the Seismological Society of America.

[2]  M. Fink,et al.  How to estimate the Green’s function of a heterogeneous medium between two passive sensors? Application to acoustic waves , 2003 .

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

[4]  Ahmed Elgamal,et al.  LOTUNG DOWNHOLE ARRAY. I: EVALUATION OF SITE DYNAMIC PROPERTIES , 1995 .

[5]  C. Langston,et al.  Effect of sinusoidal interfaces on teleseismic P-wave receiver functions , 1995 .

[6]  Fumio Yamazaki,et al.  Orientation error estimation of buried seismographs in array observation , 1992 .

[7]  Roel Snieder,et al.  Extracting the Building Response Using Seismic Interferometry: Theory and Application to the Millikan Library in Pasadena, California , 2006 .

[8]  Michel Campillo,et al.  High-Resolution Surface-Wave Tomography from Ambient Seismic Noise , 2005, Science.

[9]  J. Sheng,et al.  Interferometric/daylight seismic imaging , 2004 .

[10]  Ramesh Desikan,et al.  Receiver function analysis of the North American crust and upper mantle , 2002 .

[11]  Peter Gerstoft,et al.  Seismic interferometry-turning noise into signal , 2006 .

[12]  Vadim Levin,et al.  Receiver Functions from Multiple-Taper Spectral Correlation Estimates , 2000 .

[13]  Roel Snieder,et al.  Extraction of near-surface properties for a lossy layered medium using the propagator matrix , 2007 .

[14]  R. Snieder.  Extracting the Green's function from the correlation of coda waves: a derivation based on stationary phase. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  K. Wapenaar,et al.  Retrieving the Green’s function in an open system by cross correlation : A comparison of approaches (L) , 2005 .

[16]  Roel Snieder,et al.  Interferometric Imaging By Deconvolution: Theory And Numerical Examples , 2006 .

[17]  Hans-Peter Harjes,et al.  Imaging crustal discontinuities and the downgoing slab beneath western Crete , 2000 .

[18]  R. Weaver,et al.  On the emergence of the Green's function in the correlations of a diffuse field: pulse-echo using thermal phonons. , 2001, Ultrasonics.

[19]  S. Semken,et al.  Imaging the seismic structure of the crust and upper mantle beneath the Great Plains, Rio Grande Rift, and Colorado Plateau using receiver functions , 2005 .

[20]  Hideji Kawakami,et al.  Normalized input–output minimization analysis of wave propagation in buildings , 2003 .

[21]  Robert A. Phinney,et al.  Structure of the Earth's crust from spectral behavior of long‐period body waves , 1964 .

[22]  K. Wapenaar.  Retrieving the elastodynamic Green's function of an arbitrary inhomogeneous medium by cross correlation. , 2004, Physical review letters.

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

[24]  P. Molnar,et al.  Receiver functions in the western United States, with implications for upper mantle structure and dynamics , 2003 .

[25]  Kees Wapenaar,et al.  Retrieving the Green's Function by Cross-correlation: a Comparison of Approaches , 2004 .

[26]  Kenneth G. Dueker,et al.  Mantle discontinuity structure beneath the Colorado Rocky Mountains and High Plains , 1998 .