@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{DavidsenKwiatekCharalampidouetal.2017, author = {Davidsen, Joern and Kwiatek, Grzegorz and Charalampidou, Elli-Maria and Goebel, Thomas and Stanchits, Sergei and Rueck, Marc and Dresen, Georg}, title = {Triggering Processes in Rock Fracture}, series = {Physical review letters}, volume = {119}, journal = {Physical review letters}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.068501}, pages = {8}, year = {2017}, abstract = {We study triggering processes in triaxial compression experiments under a constant displacement rate on sandstone and granite samples using spatially located acoustic emission events and their focal mechanisms. We present strong evidence that event-event triggering plays an important role in the presence of large-scale or macrocopic imperfections, while such triggering is basically absent if no significant imperfections are present. In the former case, we recover all established empirical relations of aftershock seismicity including the Gutenberg-Richter relation, a modified version of the Omori-Utsu relation and the productivity relation-despite the fact that the activity is dominated by compaction-type events and triggering cascades have a swarmlike topology. For the Gutenberg-Richter relations, we find that the b value is smaller for triggered events compared to background events. Moreover, we show that triggered acoustic emission events have a focal mechanism much more similar to their associated trigger than expected by chance.}, language = {en} } @article{MartinezGarzonKwiatekSoneetal.2014, author = {Martinez-Garzon, Patricia and Kwiatek, Grzegorz and Sone, Hiroki and Bohnhoff, Marco and Dresen, Georg and Hartline, Craig}, title = {Spatiotemporal changes, faulting regimes, and source parameters of induced seismicity: A case study from the Geysers geothermal field}, series = {Journal of geophysical research : Solid earth}, volume = {119}, journal = {Journal of geophysical research : Solid earth}, number = {11}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1002/2014JB011385}, pages = {8378 -- 8396}, year = {2014}, abstract = {The spatiotemporal, kinematic, and source characteristics of induced seismicity occurring at different fluid injection rates are investigated to determine the predominant physical mechanisms responsible for induced seismicity at the northwestern part of The Geysers geothermal field, California. We analyze a relocated hypocenter catalog from a seismicity cluster where significant variations of the stress tensor orientation were previously observed to correlate with injection rates. We find that these stress tensor orientation changes may be related to increased pore pressure and the corresponding changes in poroelastic stresses at reservoir depth. Seismic events during peak injections tend to occur at greater distances from the injection well, preferentially trending parallel to the maximum horizontal stress direction. In contrast, at lower injection rates the seismicity tends to align in a different direction which suggests the presence of a local fault. During peak injection intervals, the relative contribution of strike-slip faulting mechanisms increases. Furthermore, increases in fluid injection rates also coincide with a decrease in b values. Our observations suggest that regardless of the injection stage, most of the induced seismicity results from thermal fracturing of the reservoir rock. However, during peak injection intervals, the increase in pore pressure may likewise be responsible for the induced seismicity. By estimating the thermal and hydraulic diffusivities of the reservoir, we confirm that the characteristic diffusion length for pore pressure is much greater than the corresponding length scale for temperature and also more consistent with the spatial extent of seismicity observed during different injection rates.}, 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{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{ZieglerReiterHeidbachetal.2015, author = {Ziegler, Moritz O. and Reiter, Karsten and Heidbach, Oliver and Zang, Arno and Kwiatek, Grzegorz and Stromeyer, Dietrich and Dahm, Torsten and Dresen, Georg and Hofmann, Gerhard}, title = {Mining-Induced Stress Transfer and Its Relation to a 1.9 Seismic Event in an Ultra-deep South African Gold Mine}, series = {Pure and applied geophysics}, volume = {172}, journal = {Pure and applied geophysics}, number = {10}, publisher = {Springer}, address = {Basel}, issn = {0033-4553}, doi = {10.1007/s00024-015-1033-x}, pages = {2557 -- 2570}, year = {2015}, abstract = {On 27 December 2007, a 1.9 seismic event occurred within a dyke in the deep-level Mponeng Gold Mine, South Africa. From the seismological network of the mine and the one from the Japanese-German Underground Acoustic Emission Research in South Africa (JAGUARS) group, the hypocentral depth (3,509 m), focal mechanism and aftershock location were estimated. Since no mining activity took place in the days before the event, dynamic triggering due to blasting can be ruled out as the cause. To investigate the hypothesis that stress transfer, due to excavation of the gold reef, induced the event, we set up a small-scale high-resolution three-dimensional (3D) geomechanical numerical model. The model consisted of the four different rock units present in the mine: quartzite (footwall), hard lava (hanging wall), conglomerate (gold reef) and diorite (dykes). The numerical solution was computed using a finite-element method with a discretised mesh of approximately elements. The initial stress state of the model is in agreement with in situ data from a neighbouring mine, and the step-wise excavation was simulated by mass removal from the gold reef. The resulting 3D stress tensor and its changes due to mining were analysed based on the Coulomb failure stress changes on the fault plane of the event. The results show that the seismic event was induced regardless of how the Coulomb failure stress changes were calculated and of the uncertainties in the fault plane solution. We also used the model to assess the seismic hazard due to the excavation towards the dyke. The resulting curve of stress changes shows a significant increase in the last in front of the dyke, indicating that small changes in the mining progress towards the dyke have a substantial impact on the stress transfer.}, 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.
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{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{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} }