TY - JOUR A1 - Weber, Michael H. A1 - Abu-Ayyash, Khalil A1 - Abueladas, Abdel-Rahman A1 - Agnon, Amotz A1 - Alasonati-Tašárová, Zuzana A1 - Al-Zubi, Hashim A1 - Babeyko, Andrey A1 - Bartov, Yuval A1 - Bauer, Klaus A1 - Becken, Michael A1 - Bedrosian, Paul A. A1 - Ben-Avraham, Zvi A1 - Bock, Günter A1 - Bohnhoff, Marco A1 - Bribach, Jens A1 - Dulski, Peter A1 - Ebbing, Joerg A1 - El-Kelani, Radwan J. A1 - Foerster, Andrea A1 - Förster, Hans-Jürgen A1 - Frieslander, Uri A1 - Garfunkel, Zvi A1 - Götze, Hans-Jürgen A1 - Haak, Volker A1 - Haberland, Christian A1 - Hassouneh, Mohammed A1 - Helwig, Stefan L. A1 - Hofstetter, Alfons A1 - Hoffmann-Rothe, Arne A1 - Jaeckel, Karl-Heinz A1 - Janssen, Christoph A1 - Jaser, Darweesh A1 - Kesten, Dagmar A1 - Khatib, Mohammed Ghiath A1 - Kind, Rainer A1 - Koch, Olaf A1 - Koulakov, Ivan A1 - Laske, Maria Gabi A1 - Maercklin, Nils T1 - Anatomy of the Dead Sea transform from lithospheric to microscopic scale N2 - Fault zones are the locations where motion of tectonic plates, often associated with earthquakes, is accommodated. Despite a rapid increase in the understanding of faults in the last decades, our knowledge of their geometry, petrophysical properties, and controlling processes remains incomplete. The central questions addressed here in our study of the Dead Sea Transform (DST) in the Middle East are as follows: (1) What are the structure and kinematics of a large fault zone? (2) What controls its structure and kinematics? (3) How does the DST compare to other plate boundary fault zones? The DST has accommodated a total of 105 km of left-lateral transform motion between the African and Arabian plates since early Miocene (similar to 20 Ma). The DST segment between the Dead Sea and the Red Sea, called the Arava/Araba Fault (AF), is studied here using a multidisciplinary and multiscale approach from the mu m to the plate tectonic scale. We observe that under the DST a narrow, subvertical zone cuts through crust and lithosphere. First, from west to east the crustal thickness increases smoothly from 26 to 39 km, and a subhorizontal lower crustal reflector is detected east of the AF. Second, several faults exist in the upper crust in a 40 km wide zone centered on the AF, but none have kilometer-size zones of decreased seismic velocities or zones of high electrical conductivities in the upper crust expected for large damage zones. Third, the AF is the main branch of the DST system, even though it has accommodated only a part (up to 60 km) of the overall 105 km of sinistral plate motion. Fourth, the AF acts as a barrier to fluids to a depth of 4 km, and the lithology changes abruptly across it. Fifth, in the top few hundred meters of the AF a locally transpressional regime is observed in a 100-300 m wide zone of deformed and displaced material, bordered by subparallel faults forming a positive flower structure. Other segments of the AF have a transtensional character with small pull-aparts along them. The damage zones of the individual faults are only 5-20 m wide at this depth range. Sixth, two areas on the AF show mesoscale to microscale faulting and veining in limestone sequences with faulting depths between 2 and 5 km. Seventh, fluids in the AF are carried downward into the fault zone. Only a minor fraction of fluids is derived from ascending hydrothermal fluids. However, we found that on the kilometer scale the AF does not act as an important fluid conduit. Most of these findings are corroborated using thermomechanical modeling where shear deformation in the upper crust is localized in one or two major faults; at larger depth, shear deformation occurs in a 20-40 km wide zone with a mechanically weak decoupling zone extending subvertically through the entire lithosphere. Y1 - 2009 UR - http://www.agu.org/journals/rg/ U6 - https://doi.org/10.1029/2008rg000264 SN - 8755-1209 ER - TY - JOUR A1 - Martinez-Garzon, Patricia A1 - Kwiatek, Grzegorz A1 - Sone, Hiroki A1 - Bohnhoff, Marco A1 - Dresen, Georg A1 - Hartline, Craig T1 - Spatiotemporal changes, faulting regimes, and source parameters of induced seismicity: A case study from the Geysers geothermal field JF - Journal of geophysical research : Solid earth N2 - 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. KW - thermal effect KW - focal mechanisms KW - geothermal KW - pore pressure KW - fluid-induced seismicity KW - reservoir characterization Y1 - 2014 U6 - https://doi.org/10.1002/2014JB011385 SN - 2169-9313 SN - 2169-9356 VL - 119 IS - 11 SP - 8378 EP - 8396 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Dresen, Georg A1 - Bohnhoff, Marco A1 - Sone, Hiroki A1 - Hartline, Craig T1 - Effects of long-term fluid injection on induced seismicity parameters and maximum magnitude in northwestern part of The Geysers geothermal field JF - Journal of geophysical research : Solid earth N2 - The long-term temporal and spatial changes in statistical, source, and stress characteristics of one cluster of induced seismicity recorded at The Geysers geothermal field (U.S.) are analyzed in relation to the field operations, fluid migration, and constraints on the maximum likely magnitude. Two injection wells, Prati-9 and Prati-29, located in the northwestern part of the field and their associated seismicity composed of 1776 events recorded throughout a 7year period were analyzed. The seismicity catalog was relocated, and the source characteristics including focal mechanisms and static source parameters were refined using first-motion polarity, spectral fitting, and mesh spectral ratio analysis techniques. The source characteristics together with statistical parameters (b value) and cluster dynamics were used to investigate and understand the details of fluid migration scheme in the vicinity of injection wells. The observed temporal, spatial, and source characteristics were clearly attributed to fluid injection and fluid migration toward greater depths, involving increasing pore pressure in the reservoir. The seasonal changes of injection rates were found to directly impact the shape and spatial extent of the seismic cloud. A tendency of larger seismic events to occur closer to injection wells and a correlation between the spatial extent of the seismic cloud and source sizes of the largest events was observed suggesting geometrical constraints on the maximum likely magnitude and its correlation to the average injection rate and volume of fluids present in the reservoir. KW - fluid-induced seismicity KW - maximum magnitude KW - reservoir characterization KW - source parameters KW - passive seismic monitoring Y1 - 2015 U6 - https://doi.org/10.1002/2015JB012362 SN - 2169-9313 SN - 2169-9356 VL - 120 IS - 10 SP - 7085 EP - 7101 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Martinez-Garzon, Patricia A1 - Kwiatek, Grzegorz A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Impact of fluid injection on fracture reactivation at The Geysers geothermal field JF - Journal of geophysical research : Solid earth N2 - We analyze the spatiotemporal distribution of fault geometries from seismicity induced by fluid injection at The Geysers geothermal field. The consistency of these faults with the local stress field is investigated using (1) the fault instability coefficient I comparing the orientation of a fault with the optimal orientation for failure in the assumed stress field and (2) the misfit angle beta between slip vectors observed from focal mechanisms and predicted from stress tensor. A statistical approach is applied to calculate the most likely fault instabilities considering the uncertainties from focal mechanisms and stress inversion. We find that faults activated by fluid injection may display a broad range in orientations. About 72% of the analyzed seismicity occurs on faults with favorable orientation for failure with respect to the stress field. However, a number of events are observed either to occur on severely misoriented faults or to slip in a different orientation than predicted from stress field. These events mostly occur during periods of high injection rates and are located in proximity to the injection wells. From the stress inversion, the friction coefficient providing the largest overall instability is mu = 0.5. About 91% of the events are activated with an estimated excess pore pressure <10 MPa, in agreement with previous models considering the combined effect of thermal and poroelastic stress changes from fluid injection. Furthermore, high seismic activity and largest magnitudes occur on favorably oriented faults with large instability coefficients and low slip misfit angles. Y1 - 2016 U6 - https://doi.org/10.1002/2016JB013137 SN - 2169-9313 SN - 2169-9356 VL - 121 SP - 7432 EP - 7449 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Ziegler, Moritz O. A1 - Heidbach, Oliver A1 - Zang, Arno A1 - Martinez-Garzon, Patricia A1 - Bohnhoff, Marco T1 - Estimation of the differential stress from the stress rotation angle in low permeable rock JF - Geophysical research letters N2 - Rotations of the principal stress axes are observed as a result of fluid injection into reservoirs. We use a generic, fully coupled 3-D thermo-hydro-mechanical model to investigate systematically the dependence of this stress rotation on different reservoir properties and injection scenarios. We find that permeability, injection rate, and initial differential stress are the key factors, while other reservoir properties only play a negligible role. In particular, we find that thermal effects do not significantly contribute to stress rotations. For reservoir types with usual differential stress and reservoir treatment the occurrence of significant stress rotations is limited to reservoirs with a permeability of less than approximately 10(-12)m(2). Higher permeability effectively prevents stress rotations to occur. Thus, according to these general findings, the observed principal stress axes rotation can be used as a proxy of the initial differential stress provided that rock permeability and fluid injection rate are known a priori. Y1 - 2017 U6 - https://doi.org/10.1002/2017GL073598 SN - 0094-8276 SN - 1944-8007 VL - 44 SP - 6761 EP - 6770 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Martínez-Garzón, Patricia A1 - Kwiatek, Grzegorz A1 - Bohnhoff, Marco A1 - Dresen, Georg T1 - Volumetric components in the earthquake source related to fluid injection and stress state JF - Geophysical research letters N2 - 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. KW - non-double-couple components KW - induced seismicity KW - geothermal KW - stress state KW - tensile opening KW - pore pressure Y1 - 2017 U6 - https://doi.org/10.1002/2016GL071963 SN - 0094-8276 SN - 1944-8007 VL - 44 IS - 2 SP - 800 EP - 809 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Plenkers, K. A1 - Leonhardt, Maria A1 - Zang, Arno A1 - von Specht, Sebastian A1 - Dresen, Georg A1 - Bohnhoff, Marco T1 - Insights into complex subdecimeter fracturing processes occurring during a water injection experiment at depth in Aspo Hard Rock Laboratory, Sweden JF - Journal of geophysical research : Solid earth N2 - 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. KW - induced seismicity KW - fracking KW - picoseismicity KW - seismomechanics KW - source parameters KW - maximum magnitude Y1 - 2018 U6 - https://doi.org/10.1029/2017JB014715 SN - 2169-9313 SN - 2169-9356 VL - 123 IS - 8 SP - 6616 EP - 6635 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Hofmann, Hannes A1 - Zimmermann, Günter A1 - Farkas, Márton Pál A1 - Huenges, Ernst A1 - Zang, Arno A1 - Leonhardt, Maria A1 - Kwiatek, Grzegorz A1 - Martinez-Garzon, Patricia A1 - Bohnhoff, Marco A1 - Min, Ki-Bok A1 - Fokker, Peter A1 - Westaway, Rob A1 - Bethmann, Falko A1 - Meier, Peter A1 - Yoon, Kern Shin A1 - Choi, JaiWon A1 - Lee, Tae Jong A1 - Kim, Kwang Yeom T1 - First field application of cyclic soft stimulation at the Pohang Enhanced Geothermal System site in Korea JF - Geophysical journal international N2 - Large-magnitude fluid-injection induced seismic events are a potential risk for geothermal energy developments worldwide. One potential risk mitigation measure is the application of cyclic injection schemes. After validation at small (laboratory) and meso (mine) scale, the concept has now been applied for the first time at field scale at the Pohang Enhanced Geothermal System (EGS) site in Korea. From 7 August until 14 August 2017 a total of 1756 m(3) of surface water was injected into Pohang well PX-1 at flow rates between 1 and 10 l s(-1), with a maximum wellhead pressure (WHP) of 22.8 MPa, according to a site-specific cyclic soft stimulation schedule and traffic light system. A total of 52 induced microearthquakes were detected in real-time during and shortly after the injection, the largest of M-w 1.9. After that event a total of 1771 m(3) of water was produced back from the well over roughly 1 month, during which time no larger-magnitude seismic event was observed. The hydraulic data set exhibits pressure-dependent injectivity increase with fracture opening between 15 and 17 MPa WHP, but no significant permanent transmissivity increase was observed. The maximum magnitude of the induced seismicity during the stimulation period was below the target threshold of M-w 2.0 and additional knowledge about the stimulated reservoir was gained. Additionally, the technical feasibility of cyclic injection at field scale was evaluated. The major factors that limited the maximum earthquake magnitude are believed to be: limiting the injected net fluid volume, flowback after the occurrence of the largest induced seismic event, using a cyclic injection scheme, the application of a traffic light system, and including a priori information from previous investigations and operations in the treatment design. KW - Cyclic soft stimulation (CSS) KW - induced seismicity KW - risk mitigation KW - enhanced geothermal systems (EGS) KW - granite KW - Pohang (Korea) Y1 - 2019 U6 - https://doi.org/10.1093/gji/ggz058 SN - 0956-540X SN - 1365-246X VL - 217 IS - 2 SP - 926 EP - 949 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Bentz, Stephan A1 - Martinez-Garzon, Patricia A1 - Kwiatek, Grzegorz A1 - Dresen, Georg A1 - Bohnhoff, Marco T1 - Analysis of Microseismicity Framing M-L > 2.5 Earthquakes at The Geysers Geothermal Field, California JF - Journal of geophysical research : Solid earth N2 - Preparatory mechanisms accompanying or leading to nucleation of larger earthquakes have been observed at both laboratory and field scales, but conditions favoring the occurrence of observable preparatory processes are still largely unknown. In particular, it remains a matter of debate why some earthquakes occur spontaneously without noticeable precursors as opposed to events that are preceded by an extended failure process. In this study, we have generated new high-resolution seismicity catalogs framing the occurrence of 20 M-L > 2.5 earthquakes at The Geysers geothermal field in California. To this end, a seismicity catalog of the 11 days framing each large event was created. We selected 20 sequences sampling different hypocentral depths and hydraulic conditions within the field. Seismic activity and magnitude frequency distributions displayed by the different earthquake sequences are correlated with their location within the reservoir. Sequences located in the northwestern part of the reservoir show overall increased seismic activity and low b values, while the southeastern part is dominated by decreased seismic activity and higher b values. Periods of high injection coincide with high b values and vice versa. These observations potentially reflect varying differential and mean stresses and damage of the reservoir rocks across the field. About 50% of analyzed sequences exhibit no change in seismicity rate in response to the large main event. However, we find complex waveforms at the onset of the main earthquake, suggesting that small ruptures spontaneously grow into or trigger larger events. KW - induced seismicity KW - earthquake nucleation KW - The Geysers KW - earthquake sequences Y1 - 2019 U6 - https://doi.org/10.1029/2019JB017716 SN - 2169-9313 SN - 2169-9356 VL - 124 IS - 8 SP - 8823 EP - 8843 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Kwiatek, Grzegorz A1 - Saarno, Tero A1 - Ader, Thomas A1 - Blümle, Felix A1 - Bohnhoff, Marco A1 - Chendorain, Michael A1 - Dresen, Georg A1 - Heikkinen, Pekka A1 - Kukkonen, Ilmo A1 - Leary, Peter A1 - Leonhardt, Maria A1 - Malin, Peter A1 - Martinez-Garzon, Patricia A1 - Passmore, Kevin A1 - Passmore, Paul A1 - Valenzuela, Sergio A1 - Wollin, Christopher T1 - Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland JF - Science Advances N2 - 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. Y1 - 2019 U6 - https://doi.org/10.1126/sciadv.aav7224 SN - 2375-2548 VL - 5 IS - 5 PB - American Association for the Advancement of Science CY - Washington ER -