Nonsystematic Rupture Directivity of Geothermal Energy Induced Microseismicity in Helsinki, Finland

The rupture behavior of microseismicity in fluid‐injection settings with low fault stresses is generally believed to be controlled by the pore pressure, including a tendency of the larger induced earthquakes to rupture into the perturbed volume toward the injection well, implying a degree of predictability. Here, we examine directivity patterns to identify fault planes and rupture directions of the 21 largest earthquakes (local magnitudes, ML ${M}_{L}$ 1.3–1.9) recorded during the 2018 St1 Deep Heat geothermal project near Helsinki, Finland. We use the Empirical Green's Function technique to retrieve per‐seismic‐station corner frequencies, earthquake durations, and directivity trends. After combining the directivity trends with focal mechanisms calculated using principle component analysis, we resolve rupture planes and rupture directions of 10 events. In contrast to studies of induced events at other sites, we find that one event rupture toward, two rupture away, and the remaining rupture parallel to the well. Furthermore, we find that the events prefer mode II failures rather than mode III. These observations provide new constraints for mechanical models of rupture growth in pore‐pressure dominated settings.

[1]  G. Hillers,et al.  Induced Earthquake Source Parameters, Attenuation, and Site Effects From Waveform Envelopes in the Fennoscandian Shield , 2022, Journal of Geophysical Research: Solid Earth.

[2]  Younghee Kim,et al.  Rupture Directivity of the 2021 ML 2.2 Gwangyang, Korea, Microearthquake: Toward Resolving High-Resolution Rupture Process of a Small Earthquake , 2022, The Seismic Record.

[3]  Zhigang Peng,et al.  Complex Source Behaviors and Spatiotemporal Evolution of Seismicity During the 2015–2016 Earthquake Sequence in Cushing, Oklahoma , 2021, Journal of Geophysical Research: Solid Earth.

[4]  D. Eaton,et al.  Large‐Scale Fracture Systems Are Permeable Pathways for Fault Activation During Hydraulic Fracturing , 2021, Journal of Geophysical Research: Solid Earth.

[5]  P. Heikkinen,et al.  Earthquake catalog of induced seismicity recorded during and after stimulation of Enhanced Geothermal System in Helsinki, Finland , 2021 .

[6]  Y. Ben‐Zion,et al.  Detection Limits and Near‐Field Ground Motions of Fast and Slow Earthquakes , 2020, Journal of Geophysical Research: Solid Earth.

[7]  G. Hillers,et al.  The 2018 Geothermal Reservoir Stimulation in Espoo/Helsinki, Southern Finland: Seismic Network Anatomy and Data Features , 2020, Seismological Research Letters.

[8]  G. Dresen,et al.  Seismicity during and after stimulation of a 6.1 km deep Enhanced Geothermal System in Helsinki, Finland , 2020 .

[9]  J. Boatwright,et al.  Directivity of M 3.1 Earthquake near Anza, California and the Effect on Peak Ground Motion , 2020 .

[10]  G. Atkinson,et al.  Stress Drops and Directivity of Induced Earthquakes in the Western Canada Sedimentary Basin , 2019, Bulletin of the Seismological Society of America.

[11]  S. Wiemer,et al.  Rupture Process of the Mw 3.3 Earthquake in the St. Gallen 2013 Geothermal Reservoir, Switzerland , 2019, Geophysical Research Letters.

[12]  M. Gutierrez,et al.  Patterns in complex hydraulic fractures observed by true-triaxial experiments and implications for proppant placement and stimulated reservoir volumes , 2019, Journal of Petroleum Exploration and Production Technology.

[13]  G. Dresen,et al.  Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland , 2019, Science Advances.

[14]  J. Verdon,et al.  Investigating the role of elastostatic stress transfer during hydraulic fracturing-induced fault activation , 2019, Geophysical Journal International.

[15]  Yihe Huang,et al.  Do Injection‐Induced Earthquakes Rupture Away from Injection Wells due to Fluid Pressure Change? , 2019, Bulletin of the Seismological Society of America.

[16]  Yajing Liu,et al.  Earthquake Stress Drop in the Charlevoix Seismic Zone, Eastern Canada , 2018, Geophysical Research Letters.

[17]  W. Fan,et al.  Investigating microearthquake finite source attributes with IRIS Community Wavefield Demonstration Experiment in Oklahoma , 2018 .

[18]  J. Renner,et al.  Sensitivity of Full Moment Tensors to Data Preprocessing and Inversion Parameters: A Case Study from the Salton Sea Geothermal Field , 2018 .

[19]  Chris Van Houtte,et al.  Improved Model Fitting for the Empirical Green's Function Approach Using Hierarchical Models , 2018 .

[20]  Rachel E. Abercrombie,et al.  Earthquake Directivity, Orientation, and Stress Drop Within the Subducting Plate at the Hikurangi Margin, New Zealand , 2017 .

[21]  P. Mai,et al.  Induced seismicity provides insight into why earthquake ruptures stop , 2017, Science Advances.

[22]  S. Wiemer,et al.  The induced earthquake sequence related to the St. Gallen deep geothermal project (Switzerland): Fault reactivation and fluid interactions imaged by microseismicity , 2017 .

[23]  V. Vavryčuk,et al.  Moment Tensor Inversion Based on the Principal Component Analysis of Waveforms: Method and Application to Microearthquakes in West Bohemia, Czech Republic , 2017 .

[24]  W. Ellsworth,et al.  Stress drops of induced and tectonic earthquakes in the central United States are indistinguishable , 2017, Science Advances.

[25]  G. Prieto,et al.  Earthquake rupture below the brittle-ductile transition in continental lithospheric mantle , 2017, Science Advances.

[26]  C. Ji,et al.  Moment rate scaling for earthquakes 3.3 ≤ M ≤ 5.3 with implications for stress drop , 2016 .

[27]  S. Shapiro,et al.  Rupture directivity of fluid‐induced microseismic events: Observations from an enhanced geothermal system , 2016 .

[28]  Rachel E. Abercrombie,et al.  Variability of earthquake stress drop in a subduction setting, the Hikurangi Margin, New Zealand , 2016, Geophysical Journal International.

[29]  Peter M. Shearer,et al.  New perspectives on self‐similarity for shallow thrust earthquakes , 2016 .

[30]  Grzegorz Kwiatek,et al.  HybridMT: A MATLAB/Shell Environment Package for Seismic Moment Tensor Inversion and Refinement , 2016 .

[31]  K. Imanishi,et al.  Small Earthquakes Deviate from the Omega‐Square Model as Revealed by Multiple Spectral Ratio Analysis , 2016 .

[32]  Valentin Gischig,et al.  The importance of earthquake interactions for injection‐induced seismicity: Retrospective modeling of the Basel Enhanced Geothermal System , 2016 .

[33]  M. Vallée,et al.  Stress-drop variability of shallow earthquakes extracted from a global database of source time functions , 2016 .

[34]  David Dempsey,et al.  Collective properties of injection‐induced earthquake sequences: 1. Model description and directivity bias , 2016 .

[35]  M. O'Sullivan,et al.  A damage mechanics approach to the simulation of hydraulic fracturing/shearing around a geothermal injection well , 2016 .

[36]  M. Ishii,et al.  Inversion for rupture properties based upon 3-D directivity effect and application to deep earthquakes in the Sea of Okhotsk region , 2015 .

[37]  V. Gischig Rupture propagation behavior and the largest possible earthquake induced by fluid injection into deep reservoirs , 2015 .

[38]  Agust Gudmundsson Elastic energy release in great earthquakes and eruptions , 2014, Front. Earth Sci..

[39]  Serge A. Shapiro,et al.  Magnitudes of induced earthquakes and geometric scales of fluid-stimulated rock volumes , 2011 .

[40]  G. Dresen,et al.  Source Parameters of Picoseismicity Recorded at Mponeng Deep Gold Mine, South Africa: Implications for Scaling Relations , 2011 .

[41]  R. Horne,et al.  Investigation of injection-induced seismicity using a coupled fluid flow and rate/state friction model , 2011 .

[42]  Won-Young Kim,et al.  The 2002 M5 Au Sable Forks, NY, earthquake sequence: Source scaling relationships and energy budget , 2010 .

[43]  R. Abercrombie,et al.  Source parameters and rupture velocity of small M ≤ 2.1 reservoir induced earthquakes , 2009 .

[44]  G. A. Prieto,et al.  A Fortran 90 library for multitaper spectrum analysis , 2009, Comput. Geosci..

[45]  G. Kwiatek Relative source time functions of seismic events at the Rudna copper mine, Poland: estimation of inversion uncertainties , 2008 .

[46]  S. Gibowicz,et al.  Comparison of source parameters estimated in the frequency and time domains for seismic events at the Rudna copper mine, Poland , 2008 .

[47]  B. Chiou,et al.  Directivity in NGA Earthquake Ground Motions: Analysis Using Isochrone Theory , 2008 .

[48]  Hiroo Kanamori,et al.  Static and Dynamic Scaling Relations for Earthquakes and Their Implications for Rupture Speed and Stress Drop , 2004 .

[49]  W. Ellsworth,et al.  Apparent break in earthquake scaling due to path and site effects on deep borehole recordings , 2003 .

[50]  B. Domański,et al.  Source Time Function of Seismic Events at Rudna Copper Mine, Poland , 2002 .

[51]  Ralph J. Archuleta,et al.  The Three-Dimensional Dynamics of Dipping Faults , 2000 .

[52]  Yoshiaki Hisada,et al.  A Theoretical Omega-Square Model Considering the Spatial Variation in Slip and Rupture Velocity , 2000 .

[53]  Dino Bindi,et al.  Source parameters of small events using constrained deconvolution with empirical Green’s functions , 1999 .

[54]  Hans-Peter Harjes,et al.  Source parameters of injection-induced microearthquakes at 9 km depth at the KTB Deep Drilling site, Germany , 1998, Bulletin of the Seismological Society of America.

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

[56]  C. Doll,et al.  Source characterization and fault plane determination for MbLg = 1.2 to 4.4 earthquakes in the Charlevoix Seismic Zone, Quebec, Canada , 1995, Bulletin of the Seismological Society of America.

[57]  J. Kakkuri,et al.  On horizontal crustal strain in Finland , 1992 .

[58]  Jim Mori,et al.  Source parameters for small events associated with the 1986 North Palm Springs, California, earthquake determined using empirical Green functions , 1990 .

[59]  Frank L. Vernon,et al.  Rupture characteristics and tomographic source imaging of ML ∼ 3 earthquakes near Anza, southern California , 1986 .

[60]  Charles S. Mueller,et al.  Source pulse enhancement by deconvolution of an empirical Green's function , 1985 .

[61]  J. Boatwright,et al.  A spectral theory for circular seismic sources; simple estimates of source dimension, dynamic stress drop, and radiated seismic energy , 1980 .

[62]  S. Hartzell Earthquake aftershocks as Green's functions , 1978 .

[63]  C. Bufe,et al.  Shear-wave attenuation along the San Andreas fault zone in central California , 1975, Bulletin of the Seismological Society of America.

[64]  J. Brune Tectonic stress and the spectra of seismic shear waves from earthquakes , 1970 .

[65]  W. E. Williams,et al.  Green's Functions , 1970, Nature.

[66]  N. A. Haskell Total energy and energy spectral density of elastic wave radiation from propagating faults , 1964 .