Fault roughness controls injection-induced seismicity.

Surface roughness ubiquitously prevails in natural faults across various length scales. Despite extensive studies highlighting the important role of fault geometry in the dynamics of tectonic earthquakes, whether and how fault roughness affects fluid-induced seismicity remains elusive. Here, we investigate the effects of fault geometry and stress heterogeneity on fluid-induced fault slip and associated seismicity characteristics using laboratory experiments and numerical modeling. We perform fluid injection experiments on quartz-rich sandstone samples containing either a smooth or a rough fault. We find that geometrical roughness slows down injection-induced fault slip and reduces macroscopic slip velocities and fault slip-weakening rates. Stress heterogeneity and roughness control hypocenter distribution, frequency-magnitude characteristics, and source mechanisms of injection-induced acoustic emissions (AEs) (analogous to natural seismicity). In contrast to smooth faults where injection-induced AEs are uniformly distributed, slip on rough faults produces spatially localized AEs with pronounced non-double-couple source mechanisms. We demonstrate that these clustered AEs occur around highly stressed asperities where induced local slip rates are higher, accompanied by lower Gutenberg-Richter b-values. Our findings suggest that real-time monitoring of induced microseismicity during fluid injection may allow identifying progressive localization of seismic activity and improve forecasting of runaway events.

Triggering Processes in Rock Fracture.

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.

A laboratory perspective on accelerating preparatory processes before earthquakes and implications for foreshock detectability.

Dynamic failure in the laboratory is commonly preceded by many foreshocks which accompany premonitory aseismic slip. Aseismic slip is also thought to govern earthquake nucleation in nature, yet, foreshocks are rare. Here, we examine how heterogeneity due to different roughness, damage and pore pressures affects premonitory slip and acoustic emission characteristics. High fluid pressures increase stiffness and reduce heterogeneity which promotes more rapid slip acceleration and shorter precursory periods, similar to the effect of low geometric heterogeneity on smooth faults. The associated acoustic emission activity in low-heterogeneity samples becomes increasingly dominated by earthquake-like double-couple focal mechanisms. The similarity of fluid pressure increase and roughness reduction suggests that increased stress and geometric homogeneity may substantially shorten the duration of foreshock activity. Gradual fault activation and extended foreshock activity is more likely observable on immature faults at shallow depth.

Earthquake energy dissipation in a fracture mechanics framework.

Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake's potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics.

No evidence of magnitude clustering in an aftershock sequence of nano- and picoseismicity.

One of the hallmarks of our current understanding of seismicity as highlighted by the epidemic-type-aftershock sequence model is that the magnitudes of earthquakes are independent of one another and can be considered as randomly drawn from the Gutenberg-Richter distribution. This assumption forms the basis of many approaches for forecasting seismicity rates and hazard assessment. Recently, it has been suggested that the assumption of independent magnitudes is not valid. It was subsequently argued that this conclusion was not supported by the original earthquake data from California. One of the main challenges is the lack of completeness of earthquake catalogs. Here, we study an aftershock sequence of nano- and picoseismicity as observed at the Mponeng mine, for which the issue of incompleteness is much less pronounced. We show that this sequence does not exhibit any significant evidence of magnitude correlations.

Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland.

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 m3 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.

Microearthquakes preceding a M4.2 Earthquake Offshore Istanbul.

A primary hurdle in observing small foreshocks is the detection-limit of most seismic networks, which is typically about magnitude M1-1.5. We show that a start-up test of a borehole-based seismic network with a much lower detection limit overcame this problem for an Mw4.2 earthquake. This earthquake occurred offshore of Istanbul, Turkey, on a fault system that is likely to rupture in an M > 7 event in the coming decades. In the three days before and two after, a total of 62 or more earthquakes, including at least 18 foreshocks, came from the mainshock source area. The signal similarity of the foreshocks shows a clear increase during the hours before the Mw4.2 mainshock. Similar foreshock sequences have recently been reported for a few well monitored M > 7 plate-boundary earthquakes. The sequence surrounding the Mw4.2 gives the impression of stochastic failures that ended up interactively unloading stress concentrations. The Mw4.2 mainshock then resulted from the accumulated release of significantly smaller events, as suggested by other field and laboratory studies.

Universal Quake Statistics: From Compressed Nanocrystals to Earthquakes.

Slowly-compressed single crystals, bulk metallic glasses (BMGs), rocks, granular materials, and the earth all deform via intermittent slips or "quakes". We find that although these systems span 12 decades in length scale, they all show the same scaling behavior for their slip size distributions and other statistical properties. Remarkably, the size distributions follow the same power law multiplied with the same exponential cutoff. The cutoff grows with applied force for materials spanning length scales from nanometers to kilometers. The tuneability of the cutoff with stress reflects "tuned critical" behavior, rather than self-organized criticality (SOC), which would imply stress-independence. A simple mean field model for avalanches of slipping weak spots explains the agreement across scales. It predicts the observed slip-size distributions and the observed stress-dependent cutoff function. The results enable extrapolations from one scale to another, and from one force to another, across different materials and structures, from nanocrystals to earthquakes.

An earthquake gap south of Istanbul.

Over the last century the North Anatolian Fault Zone in Turkey has produced a remarkable sequence of large earthquakes. These events have now left an earthquake gap south of Istanbul and beneath the Marmara Sea, a gap that has not been filled for 250 years. Here we investigate the nature of the eastern end of this gap using microearthquakes recorded by seismographs primarily on the Princes Islands offshore Istanbul. This segment lies at the western terminus of the 1999 Mw7.4 Izmit earthquake. Starting from there, we identify a 30-km-long fault patch that is entirely aseismic down to a depth of 10 km. Our evidence indicates that this patch is locked and is therefore a potential nucleation point for another Marmara segment earthquake-a potential that has significant natural hazards implications for the roughly 13 million Istanbul residents immediately to its north.

Scaling and universality in rock fracture.

We present a detailed statistical analysis of acoustic emission time series from laboratory rock fracture obtained from different experiments on different materials including acoustic emission controlled triaxial fracture and punch-through tests. In all considered cases, the waiting time distribution can be described by a unique scaling function indicating its universality. This scaling function is even indistinguishable from that for earthquakes suggesting its general validity for fracture processes independent of time, space, and magnitude scales.