Determination of shallow S-wave attenuation by down-hole waveform deconvolution: a case study in Istanbul (Turkey)

SUMMARY The estimation of shear wave velocity and attenuation in near-surface geology is of primary importance in engineering seismology. In fact, their knowledge is essential for site response studies when preparing improved seismic hazard scenarios. In this study, we propose two approaches for estimating the average shear wave quality factor Qs by using recordings of a vertical array of accelerometers. The methods are mainly based on the deconvolution of the wavefield recorded in a borehole with that recorded at the surface. The first method requires the Fourier transform of the deconvolved wavefield to be fitted with a theoretical transfer function valid for the vertical or nearly vertical (in the case at hand up to 30 ◦ incidence angle) propagation of S waves. The second method is based on the spectral fitting of the Fourier transform of only the acausal part of the deconvolved wavefield with a theoretical transfer function. Both methods can be applied without any prior knowledge of the subsoil structure (since they are based on empirical data analysis) and do not require a precise knowledge of the azimuthal orientation of the sensors in the boreholes (which is seldom available). First, we describe the theoretical framework of the proposed methodologies for Qs estimation, which are based on the assumption that the structure in the borehole is weakly heterogeneous in the vertical direction (i.e. no large impedance contrast exists between the borehole sensor and the surface). Second, by using synthetic accelerograms, we verify that in a realistic subsoil structure, the assumption of vertical homogeneity can hold and we investigate the robustness and the suitability of the proposed methods. Finally, only the method that was shown to provide the more stable results, based on fitting the borehole-to-surface spectral ratio with a theoretical function, is applied to earthquakes signals recorded by a vertical array of accelerometers installed in Atak¨ oy (western Istanbul). Results show that using borehole data provides a fair and robust estimate of an average Qs (of about 30, 46 and 99 for the 0–50, 0–70, 0–140 m depth ranges, respectively) that can be used for numerical simulations of ground motion.

[1]  Igor B. Morozov,et al.  Geometrical attenuation, frequency dependence of Q, and the absorption band problem , 2008 .

[2]  Stefano Parolai,et al.  Joint inversion of phase velocity dispersion and H/V ratio curves from seismic noise recordings using a genetic algorithm, considering higher modes , 2005 .

[3]  Erdal Safak,et al.  Models and Methods to Characterize Site Amplification from a Pair of Records , 1997 .

[4]  S. Pezeshk,et al.  Determination of Shallow Shear-Wave Attenuation in the Mississippi Embayment Using Vertical Seismic Profiling Data , 2009 .

[5]  Upmanu Lall,et al.  The Great Salt Lake: A Barometer of Low-Frequency Climatic Variability , 1995 .

[6]  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.

[7]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[8]  E. Woolery,et al.  Qs estimation for unconsolidated sediments using first-arrival SH wave critical refractions , 1994 .

[9]  H. Okada The Microtremor Survey Method , 2003 .

[10]  Peter M. Shearer,et al.  High-frequency borehole seismograms recorded in the San Jcinto Fault zone, Southern California Part 2. Attenuation and site effects , 1991, Bulletin of the Seismological Society of America.

[11]  Jeannot Trampert,et al.  SH Propagator Matrix and Qs Estimates From Borehole‐ and Surface‐Recorded Earthquake Data , 1993 .

[12]  Stefano Parolai,et al.  Joint inversion of H/V ratios and dispersion curves from seismic noise: Estimating the S‐wave velocity of bedrock , 2005 .

[13]  Roel Snieder,et al.  Equivalence of the virtual-source method and wave-field deconvolution in seismic interferometry. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  安芸 敬一 Space and time spectra of stationary stochastic waves, with special reference to microtremors , 1959 .

[15]  Peter M. Shearer,et al.  High-frequency borehole seismograms recorded in the San Jacinto Fault zone, Southern California. Part 1. Polarizations , 1991, Bulletin of the Seismological Society of America.

[16]  Marco Mucciarelli,et al.  Comparison of Empirical and Numerical Site Responses at the Tito Test Site, Southern Italy , 2007 .

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

[18]  Wei Li,et al.  Site Amplification and Attenuation via Downhole Array Seismogram Inversion: A Comparative Study of the 2003 Miyagi-Oki Aftershock Sequence , 2008 .

[19]  David Halliday,et al.  Seismic interferometry, surface waves and source distribution , 2008 .

[20]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

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

[22]  Atilla Ansal,et al.  The Ataköy vertical array (Turkey): Insights into seismic wave propagation in the shallow-most crustal layers by waveform deconvolution , 2009 .

[23]  Upmanu Lall,et al.  Decadal‐to‐centennial‐scale climate variability: Insights into the rise and fall of the Great Salt Lake , 1995 .

[24]  T. Fumal,et al.  Near-surface velocities and attenuation at two boreholes near Anza, California, from logging data , 1990, Bulletin of the Seismological Society of America.

[25]  Michel Campillo,et al.  Emergence of broadband Rayleigh waves from correlations of the ambient seismic noise , 2004 .

[26]  S. Kinoshita Deep-Borehole-Measured QP and QS Attenuation for Two Kanto Sediment Layer Sites , 2008 .

[27]  S. Parolai,et al.  S-wave Velocity Profiles for Earthquake Engineering Purposes for the Cologne Area (Germany) , 2006 .

[28]  Roel Snieder,et al.  Downhole Receiver Function: a Case Study , 2007 .

[29]  M. Ohori,et al.  A Comparison of ESAC and FK Methods of Estimating Phase Velocity Using Arbitrarily Shaped Microtremor Arrays , 2002 .

[30]  Richard L. 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.

[31]  Michael E. Mann,et al.  Spatial correlations of interdecadal variation in global surface temperatures , 1993 .

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

[33]  S. H. Seale,et al.  Site amplification and attenuation of strong ground motion , 1989 .

[34]  E. Hauksson,et al.  Results from a 1500 m deep, three-level downhole seismometer array: Site response, low Q values, and fmax , 1987 .

[35]  R. Tonn,et al.  THE DETERMINATION OF THE SEISMIC QUALITY FACTOR Q FROM VSP DATA: A COMPARISON OF DIFFERENT COMPUTATIONAL METHODS1 , 1991 .

[36]  Rongjiang Wang,et al.  A simple orthonormalization method for stable and efficient computation of Green's functions , 1999, Bulletin of the Seismological Society of America.

[37]  R. White,et al.  Measurements of Earth Attenuation from Downhole and Surface Seismic RECORDINGS , 1984 .

[38]  W. B. Joyner,et al.  The attenuation of seismic shear waves in quaternary alluvium in Santa Clara Valley, California , 1994, Bulletin of the Seismological Society of America.

[39]  T. Brocher Compressional and Shear-Wave Velocity versus Depth Relations for Common Rock Types in Northern California , 2008 .