Seismic wave attenuation from borehole and surface records in the top 2.5km beneath the city of Basel, Switzerland

SUMMARY We investigate attenuation (Q−1) of sediments of 2.5–3.5 km thickness underneath the city of Basel, Switzerland. We use recordings of 195 induced events that were obtained during and after the stimulation of a reservoir for a Deep Heat Mining Project in 2006 and 2007. The data set is ideally suited to estimate Q as all events are confined to a small source volume and were recorded by a dense surface network as well as six borehole sensors at various depths. The deepest borehole sensor is positioned at a depth of 2.7 km inside the crystalline basement at a mean hypocentral distance of 1.8 km. This allows us to measure Q for frequencies between 10 and 130 Hz. We apply two different methods to estimate Q. First, we use a standard spectral ratio technique to obtain Q, and as a second measure we estimate Q in the time domain, by convolving signals recorded by the deepest sensor with a Q operator and then comparing the convolved signals to recordings at the shallower stations. Both methods deliver comparable values for Q. We also observe similar attenuation for P- and S- waves (QP∼QS). As expected, Q increases with depth, but with values around 30–50, it is low even for the consolidated Permian and Mesozoic sediments between 500 and 2700 m.

[1]  K. Seitz,et al.  Determination of displacements in the upper Rhine graben Area from GPS and leveling data , 2005 .

[2]  Donat Fäh,et al.  The 1356 Basel earthquake: an interdisciplinary revision , 2009 .

[3]  Rachel E. Abercrombie,et al.  Near-surface attenuation and site effects from comparison of surface and deep borehole recordings , 1997, Bulletin of the Seismological Society of America.

[4]  Donat Fäh,et al.  Ground motion and macroseismic intensities of a seismic event related to geothermal reservoir stimulation below the city of Basel—observations and modelling , 2009 .

[5]  G. M. Jackson,et al.  Principal component transforms of triaxial recordings by singular value decomposition , 1991 .

[6]  R. Parker,et al.  Reducing the bias of multitaper spectrum estimates , 2007 .

[7]  N. Deichmann Far-field pulse shapes from circular sources with variable rupture velocities , 1997, Bulletin of the Seismological Society of America.

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

[9]  Hiroaki Niitsuma,et al.  Analysis of Microseismic Events from a Stimulation at Basel, Switzerland , 2007 .

[10]  Donat Fäh,et al.  Earthquakes in Switzerland and surrounding regions during 2006 , 2007 .

[11]  Gregory C. Beroza,et al.  Does apparent stress vary with earthquake size? , 2001 .

[12]  Peter Huggenberger,et al.  The southeastern border of the Upper Rhine Graben: a 3D geological model and its importance for tectonics and groundwater flow , 2005 .

[13]  Tamao Sato,et al.  Seismic radiation from circular cracks growing at variable rupture velocity , 1994, Bulletin of the Seismological Society of America.

[14]  Anna Sommaruga Décollement tectonics in the Jura forelandfold-and-thrust belt , 1999 .

[15]  Benjamin Edwards,et al.  Attenuation of seismic shear wave energy in Switzerland , 2011 .

[16]  P. Blümling,et al.  Recent stress field and neotectonics in the Eastern Jura Mountains, Switzerland , 1987 .

[17]  Stefan M. Schmid,et al.  Latest Pliocene to recent thick-skinned tectonics at the Upper Rhine Graben – Jura Mountains junction , 2007 .

[18]  Kurt S. Riedel,et al.  Minimum bias multiple taper spectral estimation , 2018, IEEE Trans. Signal Process..

[19]  Thomas P. Bronez,et al.  On the performance advantage of multitaper spectral analysis , 1992, IEEE Trans. Signal Process..

[20]  Donat Fäh,et al.  Deterministic earthquake scenario for the Basel area: Simulating strong motions and site effects for Basel, Switzerland , 2005 .

[21]  D. Jongmans,et al.  Microearthquake S-wave observations from 0 to 1 km in the Varian well at Parkfield, California , 1995, Bulletin of the Seismological Society of America.

[22]  R. Abercrombie A Summary of Attenuation Measurements from Borehole Recordings of Earthquakes: The 10 Hz Transition Problem , 1998 .

[23]  W. Walter,et al.  Evidence for non‐constant energy/moment scaling from coda‐derived source spectra , 2005 .

[24]  K. Aki,et al.  Separation of intrinsic and scattering attenuation in southern California using TERRAscope data , 1994 .

[25]  T. Spillmann,et al.  Application of microseismic multiplet analysis to the Basel geothermal reservoir stimulation events , 2010 .

[26]  B. Dyer,et al.  Microseismic imaging of a geothermal reservoir stimulation , 2008 .

[27]  John A. Orcutt,et al.  Surface and near-surface effects on seismic waves—theory and borehole seismometer results , 1987 .

[28]  Donat Fäh,et al.  The 1356 earthquake: an interdisciplinary revision , 2009 .

[29]  Domenico Giardini,et al.  Earthquakes Induced By the Stimulation of an Enhanced Geothermal System Below Basel (Switzerland) , 2009 .

[30]  Mitsuyuki Hoshiba,et al.  Separation of scattering attenuation and intrinsic absorption in Japan using the multiple lapse time window analysis of full seismogram envelope , 1993 .

[31]  T. Urbancic,et al.  Rock-mass characterization using intrinsic and scattering attenuation estimates at frequencies from 400 to 1600 Hz , 1996 .

[32]  A Summary of Attenuation Measurements from Borehole Recordings of Earthquakes: The 10 Hz Transition Problem , 1998 .

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

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

[35]  A Seismic-Ray Tracing Method for Two-Dimensional Inhomogeneous Media , 1976 .

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

[37]  H. Kanamori,et al.  Determination of earthquake energy release and ML using TERRAscope , 1993 .

[38]  George L. Choy,et al.  Global patterns of radiated seismic energy and apparent stress , 1995 .

[39]  Peter M. Shearer,et al.  Earthquake source scaling and self-similarity estimation by stacking P and S spectra , 2004 .

[40]  W. Ellsworth,et al.  Observations of Earthquake Source Parameters at 2 km Depth in the Long Valley Caldera, Eastern California , 2001 .

[41]  S. K. Singh,et al.  Seismic energy release in Mexican subduction zone earthquakes , 1994 .

[42]  K. Evans,et al.  Stress orientation to 5 km depth in the basement below Basel (Switzerland) from borehole failure analysis , 2009 .

[43]  Ru-Shan Wu,et al.  Multiple scattering and energy transfer of seismic waves - separation of scattering effect from intrinsic attenuation - I. Theoretical modelling , 1985 .

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

[45]  A. Becker The Jura Mountains — an active foreland fold-and-thrust belt? , 2000 .

[46]  Rachel E. Abercrombie,et al.  Earthquake source scaling relationships from −1 to 5 ML using seismograms recorded at 2.5‐km depth , 1995 .

[47]  Markus Häring,et al.  Characterisation of the Basel 1 enhanced geothermal system , 2008 .

[48]  Jacques R. Ernst,et al.  Earthquake focal mechanisms of the induced seismicity in 2006 and 2007 below Basel (Switzerland) , 2009 .