The earthquake size distribution follows, in most instances, a power law, with the slope of this power law, the 'b value', commonly used to describe the relative occurrence of large and small events ...(a high b value indicates a larger proportion of small earthquakes, and vice versa). Statistically significant variations of b values have been measured in laboratory experiments, mines and various tectonic regimes such as subducting slabs, near magma chambers, along fault zones and in aftershock zones. However, it has remained uncertain whether these differences are due to differing stress regimes, as it was questionable that samples in small volumes (such as in laboratory specimens, mines and the shallow Earth's crust) are representative of earthquakes in general. Given the lack of physical understanding of these differences, the observation that b values approach the constant 1 if large volumes are sampled was interpreted to indicate that b = 1 is a universal constant for earthquakes in general. Here we show that the b value varies systematically for different styles of faulting. We find that normal faulting events have the highest b values, thrust events the lowest and strike-slip events intermediate values. Given that thrust faults tend to be under higher stress than normal faults we infer that the b value acts as a stress meter that depends inversely on differential stress.
Due to the deep socioeconomic implications, induced seismicity is a timely and increasingly relevant topic of interest for the general public. Cases of induced seismicity have a global distribution ...and involve a large number of industrial operations, with many documented cases from as far back to the beginning of the twentieth century. However, the sparse and fragmented documentation available makes it difficult to have a clear picture on our understanding of the physical phenomenon and consequently in our ability to mitigate the risk associated with induced seismicity. This review presents a unified and concise summary of the still open questions related to monitoring, discrimination, and management of induced seismicity in the European context and, when possible, provides potential answers. We further discuss selected critical European cases of induced seismicity, which led to the suspension or reduction of the related industrial activities.
Key Points
We provide a unified and concise summary about the still open questions on monitoring, discrimination, and management of induced seismicity
We review critical cases of induced seismicity in Europe which led to the suspension of the related industrial activities
This study outlines the scientific and societal challenges posed by the induced seismicity in a European perspective
Laboratory experiments highlight a systematic b value decrease during the stress increase period before failure, and some large natural events are known to show a precursory decrease in the b value. ...However, short‐term forecast models currently consider only the generic probability that an event can trigger subsequent seismicity in the near field. While the probability increase over a stationary Poissonian background is substantial, selected case studies have shown through cost‐benefit analysis that the absolute main shock probability remains too low to warrant significant mitigation actions. We analyze the probabilities considering both changes in the seismicity rates and temporal changes in the b value. The precursory b value decrease in the 2009 L'Aquila case results in an additional fiftyfold probability increase for a M6.3 event. Translated into time‐varying hazard and risk, these changes surpass the cost‐benefit threshold for short‐term evacuation.
Key Points
L'Aquila's foreshocks show significantly lower b values than background seismicity
Using time‐dependent b values increases forecasted main shock probabilities
Such probabilities surpass levels to retrospectively warrant short‐term evacuation
Immediately after a large earthquake, the main question asked by the public and decision-makers is whether it was the mainshock or a foreshock to an even stronger event yet to come. So far, ...scientists can only offer empirical evidence from statistical compilations of past sequences, arguing that normally the aftershock sequence will decay gradually whereas the occurrence of a forthcoming larger event has a probability of a few per cent. Here we analyse the average size distribution of aftershocks of the recent Amatrice-Norcia and Kumamoto earthquake sequences, and we suggest that in many cases it may be possible to discriminate whether an ongoing sequence represents a decaying aftershock sequence or foreshocks to an upcoming large event. We propose a simple traffic light classification to assess in real time the level of concern about a subsequent larger event and test it against 58 sequences, achieving a classification accuracy of 95 per cent.
We estimate Coulomb stress variations (ΔCFS) by considering interactions among 163 earthquakes (ML ≤ 3.4) that occurred during the hydraulic stimulation of the enhanced geothermal system in Basel, ...Switzerland, in 2006. We observe that overall 75% of event locations are consistent with positive ΔCFS. The performance of the model increases with time and distance from injection, accordingly with the presumed less dominant role of the pore pressure further from the injection well and after shut‐in. We analyze the sensitivity of results to location and focal mechanism uncertainties, the fault‐plane ambiguity, and the friction coefficient. We find that ΔCFS are highly sensitive to location accuracy but robust with regard to uncertainties of the other parameters. Our results suggest that (i) the Coulomb model may complement conventional pore‐pressure triggering models and (ii) implementing it for time‐dependent seismic hazard assessment during fluid injection may improve the forecasting ability but will require highly accurate hypocenter estimates.
Key Points
Static stress changes are estimated for an injection‐induced sequence
Earthquake interactions alone play a role in the triggering process
A sensitivity analysis of the Coulomb model is performed
During deep reservoir engineering projects, in which permeability is enhanced by high-pressure fluid injection, seismicity is invariably induced, posing nuisance to the local population and a ...potential hazard for structures. Hazard and risk assessment tools that can operate in real-time during reservoir stimulation depend on the ability to efficiently model induced seismicity. We here propose a novel modelling approach based on a combination of physical considerations and stochastic elements. It can model a large number of synthetic event catalogues, and at the same time is constrained by observations of hydraulic behaviour in the injection well. We model fluid flow using non-linear pressure diffusion equations, in which permeability increases irreversibly above a prescribed pressure threshold. The transient pressure field is used to trigger events at so-called 'seed points' that are distributed randomly in space and represent potential earthquake hypocentres. We assign to each seed point a differential stress based on the mean estimates of the in situ stress field and add a normal distributed random value. Assuming a fault orientation with respect to the stress field and a Mohr-Coulomb failure criterion, we evaluate at each time step, if a seed point is triggered through a pressure increase. A negative proportional relationship between differential stress and b values is further assumed as observed from tectonic earthquakes and in laboratory experiments. As soon as an event is triggered, we draw a random magnitude from a power-law distribution with a b value corresponding to the differential stress at the triggered seed point. We thus obtain time-dependent catalogues of seismic events including magnitude. The strategy of modelling flow and seismicity in a decoupled manner ensures efficiency and flexibility of the model. The model parameters are calibrated using observations from the Basel deep geothermal experiment in 2006. We are able to reproduce the hydraulic behaviour, the space-time evolution of the seismicity and its frequency-magnitude distribution. A large number of simulations of the calibrated model are then used to capture the variability of the process, an important input to compute probabilistic seismic hazard. We also use the calibrated model to explore alternative injection scenarios by varying injection volume, pressure as well as depth, and show the possible effect of those parameters on seismic hazard.
We explore the role of earthquake interactions during an injection‐induced seismic sequence. We propose a model, which considers both a transient pressure and static stress redistribution due to ...event interactions as triggering mechanisms. By calibrating the model against observations at the Enhanced Geothermal System of Basel, Switzerland, we are able to reproduce the time behavior of the seismicity rate. We observe that considering earthquake interactions in the modeling leads to a larger number of expected seismic events (24% more) if compared to a pressure‐induced seismicity only. The increase of the number of events is particularly evident after the end of the injection. We conclude that implementing a model for estimating the static stress changes due to mutual event interactions increases significantly the understanding of the process and the behavior of induced seismicity.
Key Points
We model synthetic catalogues for induced seismicity accounting for earthquake interactions in terms of static stress transfer
Static stress interactions affect the number of seismic events, which increase (24% more) with respect to a pressure‐induced only case
Static stress transfer has a significant role in understanding the rate and total number of events, as well as their spatial distribution
We determine stress drops from P‐wave spectra of about 1000 earthquakes induced by hydraulic stimulation in crystalline rock for a deep heat mining project in Basel, Switzerland. We observe an ...increase in stress drop by about a factor of five with radial distance from 10 m to 300 m, which suggests that stress drop correlates with pore pressure perturbations due to the injection. We test this hypothesis by calculating the injection‐related pore pressure perturbation based on a simple linear pore pressure diffusion model and find a good correlation of the expected pore pressure perturbation with the estimated stress drops.
Key Points
We compute stress drop from the source spectra of induced microearthquakes
Stress drop increases with distance from the injection point
We seek to explain this with pore pressure diffusion