@article{DurandBentzKwiateketal.2020,
  author    = {Durand, Virginie and Bentz, Stephan and Kwiatek, Grzegorz and Dresen, Georg and Wollin, Christopher and Heidbach, Oliver and Martinez-Garzon, Patricia and Cotton, Fabrice and Nurlu, Murat and Bohnhoff, Marco},
  title     = {A two-scale preparation phase preceded an M-w 5.8 earthquake in the sea of marmara offshore Istanbul, Turkey},
  series = {Seismological research letters},
  volume    = {91},
  journal   = {Seismological research letters},
  number    = {6},
  address   = {Boulder},
  issn      = {0895-0695},
  doi       = {10.1785/0220200110},
  pages     = {3139 -- 3147},
  year      = {2020},
  abstract  = {We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.},
  language  = {en}
}
@article{WangKwiatekRybackietal.2020,
  author    = {Wang, Lei and Kwiatek, Grzegorz and Rybacki, Erik and Bonnelye, Audrey and Bohnhoff, Marco and Dresen, Georg},
  title     = {Laboratory study on fluid-induced fault slip behavior: the role of fluid pressurization rate},
  series = {Geophysical research letters : GRL},
  volume    = {47},
  journal   = {Geophysical research letters : GRL},
  number    = {6},
  publisher = {Wiley},
  address   = {Hoboken, NJ},
  issn      = {0094-8276},
  doi       = {10.1029/2019GL086627},
  pages     = {12},
  year      = {2020},
  abstract  = {Understanding the physical mechanisms governing fluid-induced fault slip is important for improved mitigation of seismic risks associated with large-scale fluid injection. We conducted fluid-induced fault slip experiments in the laboratory on critically stressed saw-cut sandstone samples with high permeability using different fluid pressurization rates. Our experimental results demonstrate that fault slip behavior is governed by fluid pressurization rate rather than injection pressure. Slow stick-slip episodes (peak slip velocity < 4 mu m/s) are induced by fast fluid injection rate, whereas fault creep with slip velocity < 0.4 mu m/s mainly occurs in response to slow fluid injection rate. Fluid-induced fault slip may remain mechanically stable for loading stiffness larger than fault stiffness. Independent of fault slip mode, we observed dynamic frictional weakening of the artificial fault at elevated pore pressure. Our observations highlight that varying fluid injection rates may assist in reducing potential seismic hazards of field-scale fluid injection projects. <br /> Plain Language Summary Human-induced earthquakes from field-scale fluid injection projects including enhanced geothermal system and deep wastewater injection have been documented worldwide. Although it is clear that fluid pressure plays a crucial role in triggering fault slip, the physical mechanism behind induced seismicity still remains poorly understood. We performed laboratory tests, and here we present two fluid-induced slip experiments conducted on permeable Bentheim sandstone samples crosscut by a fault that is critically stressed. Fault slip is then triggered by pumping the water from the bottom end of the sample at different fluid injection rates. Our results show that fault slip is controlled by fluid pressure increase rate rather than by the absolute magnitude of fluid pressure. In contrast to episodes of relatively rapid but stable sliding events caused by a fast fluid injection rate, fault creep is observed during slow fluid injection. Strong weakening of the dynamic friction coefficient of the experimental fault is observed at elevated pore pressure, independent of fault slip mode. These results may provide a better understanding of the complex behavior of fluid-induced fault slip on the field scale.},
  language  = {en}
}
@article{BentzKwiatekMartinezGarzonetal.2020,
  author    = {Bentz, Stephan and Kwiatek, Grzegorz and Martinez-Garzon, Patricia and Bohnhoff, Marco and Dresen, Georg},
  title     = {Seismic moment evolution during hydraulic stimulations},
  series = {Geophysical research letters},
  volume    = {47},
  journal   = {Geophysical research letters},
  number    = {5},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1029/2019GL086185},
  pages     = {9},
  year      = {2020},
  abstract  = {Analysis of past and present stimulation projects reveals that the temporal evolution and growth of maximum observed moment magnitudes may be linked directly to the injected fluid volume and hydraulic energy. Overall evolution of seismic moment seems independent of the tectonic stress regime and is most likely governed by reservoir specific parameters, such as the preexisting structural inventory. Data suggest that magnitudes can grow either in a stable way, indicating the constant propagation of self-arrested ruptures, or unbound, for which the maximum magnitude is only limited by the size of tectonic faults and fault connectivity. Transition between the two states may occur at any time during injection or not at all. Monitoring and traffic light systems used during stimulations need to account for the possibility of unstable rupture propagation from the very beginning of injection by observing the entire seismicity evolution in near-real time and at high resolution for an immediate reaction in injection strategy. Plain Language Summary Predicting and controlling the size of earthquakes caused by fluid injection is currently the major concern of many projects associated with geothermal energy production. Here, we analyze the magnitude and seismic moment evolution with injection parameters for prominent geothermal and scientific projects to date. Evolution of seismicity seems to be largely independent of the tectonic stress background and seemingly depends on reservoir specific characteristics. We find that the maximum observed magnitudes relate linearly to the injected volume or hydraulic energy. A linear relation suggests stable growth of induced ruptures, as predicted by current models, or rupture growth may no longer depend on the stimulated volume but on tectonics. A system may change between the two states during the course of fluid injection. Close-by and high-resolution monitoring of seismic and hydraulic parameters in near-real time may help identify these fundamental changes in ample time to change injection strategy and manage maximum magnitudes.},
  language  = {en}
}
@article{DresenKwiatekGoebeletal.2020,
  author    = {Dresen, Georg and Kwiatek, Grzegorz and Goebel, Thomas and Ben-Zion, Yehuda},
  title     = {Seismic and aseismic preparatory processes before large stick-slip failure},
  series = {Pure and applied geophysics},
  volume    = {177},
  journal   = {Pure and applied geophysics},
  number    = {12},
  publisher = {Springer},
  address   = {Basel},
  issn      = {0033-4553},
  doi       = {10.1007/s00024-020-02605-x},
  pages     = {5741 -- 5760},
  year      = {2020},
  abstract  = {Natural earthquakes often have very few observable foreshocks which significantly complicates tracking potential preparatory processes. To better characterize expected preparatory processes before failures, we study stick-slip events in a series of triaxial compression tests on faulted Westerly granite samples. We focus on the influence of fault roughness on the duration and magnitude of recordable precursors before large stick-slip failure. Rupture preparation in the experiments is detectable over long time scales and involves acoustic emission (AE) and aseismic deformation events. Preparatory fault slip is found to be accelerating during the entire pre-failure loading period, and is accompanied by increasing AE rates punctuated by distinct activity spikes associated with large slip events. Damage evolution across the fault zones and surrounding wall rocks is manifested by precursory decrease of seismic b-values and spatial correlation dimensions. Peaks in spatial event correlation suggest that large slip initiation occurs by failure of multiple asperities. Shear strain estimated from AE data represents only a small fraction (< 1\%) of total shear strain accumulated during the preparation phase, implying that most precursory deformation is aseismic. The relative contribution of aseismic deformation is amplified by larger fault roughness. Similarly, seismic coupling is larger for smooth saw-cut faults compared to rough faults. The laboratory observations point towards a long-lasting and continuous preparation process leading to failure and large seismic events. The strain partitioning between aseismic and observable seismic signatures depends on fault structure and instrument resolution.},
  language  = {en}
}
@article{WangKwiatekRybackietal.2020,
  author    = {Wang, Lei and Kwiatek, Grzegorz and Rybacki, Erik and Bohnhoff, Marco and Dresen, Georg},
  title     = {Injection-induced seismic moment release and laboratory fault slip},
  series = {Geophysical research letters},
  volume    = {47},
  journal   = {Geophysical research letters},
  number    = {22},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1029/2020GL089576},
  pages     = {11},
  year      = {2020},
  abstract  = {Understanding the relation between injection-induced seismic moment release and operational parameters is crucial for early identification of possible seismic hazards associated with fluid-injection projects. We conducted laboratory fluid-injection experiments on permeable sandstone samples containing a critically stressed fault at different fluid pressurization rates. The observed fluid-induced fault deformation is dominantly aseismic. Fluid-induced stick-slip and fault creep reveal that total seismic moment release of acoustic emission (AE) events is related to total injected volume, independent of respective fault slip behavior. Seismic moment release rate of AE scales with measured fault slip velocity. For injection-induced fault slip in a homogeneous pressurized region, released moment shows a linear scaling with injected volume for stable slip (steady slip and fault creep), while we find a cubic relation for dynamic slip. Our results highlight that monitoring evolution of seismic moment release with injected volume in some cases may assist in discriminating between stable slip and unstable runaway ruptures.},
  language  = {en}
}
@article{MartinezGarzonKwiatekBohnhoffetal.2017,
  author    = {Mart{\´i}nez-Garz{\´o}n, Patricia and Kwiatek, Grzegorz and Bohnhoff, Marco and Dresen, Georg},
  title     = {Volumetric components in the earthquake source related to fluid injection and stress state},
  series = {Geophysical research letters},
  volume    = {44},
  journal   = {Geophysical research letters},
  number    = {2},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1002/2016GL071963},
  pages     = {800 -- 809},
  year      = {2017},
  abstract  = {We investigate source processes of fluid-induced seismicity from The Geysers geothermal reservoir in California to determine their relation with hydraulic operations and improve the corresponding seismic hazard estimates. Analysis of 869 well-constrained full moment tensors (M-w 0.8-3.5) reveals significant non-double-couple components (>25\%) for about 65\% of the events. Volumetric deformation is governed by cumulative injection rates with larger non-double-couple components observed near the wells and during high injection periods. Source mechanisms are magnitude dependent and vary significantly between faulting regimes. Normal faulting events (M-w<2) reveal substantial volumetric components indicating dilatancy in contrast to strike-slip events that have a dominant double-couple source. Volumetric components indicating closure of cracks in the source region are mostly found for reverse faulting events with M-w>2.5. Our results imply that source processes and magnitudes of fluid-induced seismic events are strongly affected by the hydraulic operations, the reservoir stress state, and the faulting regime.},
  language  = {en}
}
@article{ZangStephanssonStenbergetal.2017,
  author    = {Zang, Arno and Stephansson, Ove and Stenberg, Leif and Plenkers, Katrin and von Specht, Sebastian and Milkereit, Claus and Schill, Eva and Kwiatek, Grzegorz and Dresen, Georg and Zimmermann, G{\"u}nter and Dahm, Torsten and Weber, Michael},
  title     = {Hydraulic fracture monitoring in hard rock at 410 m depth with an advanced fluid-injection protocol and extensive sensor array},
  series = {Geophysical journal international},
  volume    = {208},
  journal   = {Geophysical journal international},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  pages     = {790 -- 813},
  year      = {2017},
  abstract  = {In this paper, an underground experiment at the Aspo Hard Rock Laboratory (HRL) is described. Main goal is optimizing geothermal heat exchange in crystalline rock mass at depth by multistage hydraulic fracturing with minimal impact on the environment, that is, seismic events. For this, three arrays with acoustic emission, microseismicity and electromagnetic sensors are installed mapping hydraulic fracture initiation and growth. Fractures are driven by three different water injection schemes (continuous, progressive and pulse pressurization). After a brief review of hydraulic fracture operations in crystalline rock mass at mine scale, the site geology and the stress conditions at Aspo HRL are described. Then, the continuous, single-flow rate and alternative, multiple-flow rate fracture breakdown tests in a horizontal borehole at depth level 410 m are described together with the monitoring networks and sensitivity. Monitoring results include the primary catalogue of acoustic emission hypocentres obtained from four hydraulic fractures with the in situ trigger and localizing network. The continuous versus alternative water injection schemes are discussed in terms of the fracture breakdown pressure, the fracture pattern from impression packer result and the monitoring at the arrays. An example of multistage hydraulic fracturing with several phases of opening and closing of fracture walls is evaluated using data from acoustic emissions, seismic broad-band recordings and electromagnetic signal response. Based on our limited amount of in situ tests (six) and evaluation of three tests in Avro granodiorite, in the multiple-flow rate test with progressively increasing target pressure, the acoustic emission activity starts at a later stage in the fracturing process compared to the conventional fracturing case with continuous water injection. In tendency, also the total number and magnitude of acoustic events are found to be smaller in the progressive treatment with frequent phases of depressurization.},
  language  = {en}
}
@article{KwiatekMartinezGarzonPlenkersetal.2018,
  author    = {Kwiatek, Grzegorz and Martinez-Garzon, Patricia and Plenkers, K. and Leonhardt, Maria and Zang, Arno and von Specht, Sebastian and Dresen, Georg and Bohnhoff, Marco},
  title     = {Insights into complex subdecimeter fracturing processes occurring during a water injection experiment at depth in Aspo Hard Rock Laboratory, Sweden},
  series = {Journal of geophysical research : Solid earth},
  volume    = {123},
  journal   = {Journal of geophysical research : Solid earth},
  number    = {8},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1029/2017JB014715},
  pages     = {6616 -- 6635},
  year      = {2018},
  abstract  = {We investigate the source characteristics of picoseismicity (M-w < -2) recorded during a hydraulic fracturing in situ experiment performed in the underground Aspo Hard Rock Laboratory, Sweden. The experiment consisted of six stimulations driven by three different water injection schemes and was performed inside a 28-m-long, horizontal borehole located at 410-m depth. The fracturing processes were monitored with a variety of seismic networks including broadband seismometers, geophones, high-frequency accelerometers, and acoustic emission sensors thereby covering a wide frequency band between 0.01 and 100,000Hz. Here we study the high-frequency signals with dominant frequencies exceeding 1000 Hz. The combined seismic network allowed for detection and detailed analysis of 196 small-scale seismic events with moment magnitudes M-W < -3.5 (source sizes of decimeter scale) that occurred solely during the stimulations and shortly after. The double-difference relocated hypocenter catalog as well as source parameters were used to study the physical characteristics of the induced seismicity and then compared to the stimulation parameters. We observe a spatiotemporal migration of the picoseismic events away and toward the injection intervals in direct correlation with changes in the hydraulic energy (product of fluid injection pressure and injection rate). We find that the total radiated seismic energy is extremely low with respect to the product of injected fluid volume and pressure (hydraulic energy). The radiated seismic energy correlates well with the hydraulic energy rate. The obtained fault plane solutions for particularly well-characterized events signify the reactivation of preexisting rock defects under influence of increased pore fluid pressure on fault plane orientations in good correspondence with the local stress field orientation.},
  language  = {en}
}
@article{MartinezGarvonBohnhoffMencinetal.2019,
  author    = {Martinez-Garvon, Patricia and Bohnhoff, Marco and Mencin, David and Kwiatek, Grzegorz and Dresen, Georg and Hodgkinson, Kathleen and Nurlu, Murat and Kadirioglu, Filiz Tuba and Kartal, Recai Feyiz},
  title     = {Slow strain release along the eastern Marmara region offshore Istanbul in conjunction with enhanced local seismic moment release},
  series = {Earth \& planetary science letters},
  volume    = {510},
  journal   = {Earth \& planetary science letters},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2019.01.001},
  pages     = {209 -- 218},
  year      = {2019},
  abstract  = {We analyze a large transient strainmeter signal recorded at 62.5 m depth along the southern shore of the eastern Sea of Marmara region in northwestern Turkey. This region represents a passage of stress transfer from the Izmit rupture to the Marmara seismic gap. The strain signal was recorded at the Esenkoy site by one of the ICDP-GONAF (International Continental Drilling Programme - Geophysical Observatory at the North Anatolian Fault) strainmeters on the Armutlu peninsula with a maximum amplitude of 5 microstrain and lasting about 50 days. The onset of the strain signal coincided with the origin time of a M-w 4.4 earthquake offshore Yalova, which occurred as part of a seismic sequence including eight M-w >= 3.5 earthquakes. The Mw 4.4 event occurred at a distance of about 30 km from Esenkoy on June 25th 2016 representing the largest earthquake in this region since 2008. Before the event, the maximum horizontal strain was subparallel to the regional maximum horizontal stress derived from stress inversion of local seismicity. During the strain transient, we observe a clockwise rotation in the local horizontal strain field of about 20 degrees. The strain signal does not correlate with known environmental parameters such as annual changes of sea level, rainfall or temperature. The strain signal could indicate local slow slip on the Cinarcik fault and thus a transfer of stress to the eastern Marmara seismic gap.},
  language  = {en}
}
@article{KwiatekSaarnoAderetal.2019,
  author    = {Kwiatek, Grzegorz and Saarno, Tero and Ader, Thomas and Bl{\"u}mle, Felix and Bohnhoff, Marco and Chendorain, Michael and Dresen, Georg and Heikkinen, Pekka and Kukkonen, Ilmo and Leary, Peter and Leonhardt, Maria and Malin, Peter and Martinez-Garzon, Patricia and Passmore, Kevin and Passmore, Paul and Valenzuela, Sergio and Wollin, Christopher},
  title     = {Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland},
  series = {Science Advances},
  volume    = {5},
  journal   = {Science Advances},
  number    = {5},
  publisher = {American Association for the Advancement of Science},
  address   = {Washington},
  issn      = {2375-2548},
  doi       = {10.1126/sciadv.aav7224},
  pages     = {11},
  year      = {2019},
  abstract  = {We show that near-real-time seismic monitoring of fluid injection allowed control of induced earthquakes during the stimulation of a 6.1-km-deep geothermal well near Helsinki, Finland. A total of 18,160 m(3) of fresh water was pumped into crystalline rocks over 49 days in June to July 2018. Seismic monitoring was performed with a 24-station borehole seismometer network. Using near-real-time information on induced-earthquake rates, locations, magnitudes, and evolution of seismic and hydraulic energy, pumping was either stopped or varied-in the latter case, between well-head pressures of 60 and 90 MPa and flow rates of 400 and 800 liters/min. This procedure avoided the nucleation of a project-stopping magnitude M-W 2.0 induced earthquake, a limit set by local authorities. Our results suggest a possible physics-based approach to controlling stimulation-induced seismicity in geothermal projects.},
  language  = {en}
}