Coulomb threshold rate-and-state model for fault reactivation: application to induced seismicity at Groningen

• Published in: Geophysical journal international Vol. 228; no. 3; pp. 2061 - 2072
• Main Authors: Heimisson, Elías R, Smith, Jonathan D, Avouac, Jean-Philippe, Bourne, Stephen J
• Format: Journal Article
• Language: English
• Published: Oxford University Press 01.03.2022
• Subjects: Instability analysis; Dynamics and mechanics of faulting; Earthquake interaction, forecasting, and prediction; Earthquake hazards; Statistical seismology; Europe
• ISSN: 0956-540X; 1365-246X
• DOI: 10.1093/gji/ggab467

SUMMARY A number of recent modelling studies of induced seismicity have used the 1994 rate-and-state friction model of Dieterich 1994 to account for the fact that earthquake nucleation is not instantaneous. Notably, the model assumes a population of seismic sources accelerating towards instability with a distribution of initial slip speeds such that they would produce earthquakes steadily in the absence of any perturbation to the system. This assumption may not be valid in typical intraplate settings where most examples of induced seismicity occur, since these regions have low stressing rates and initially low seismic activity. The goal of this paper is twofold. First, to derive a revised Coulomb rate-and-state model, which takes into account that seismic sources can be initially far from instability. Second, to apply and test this new model, called the Threshold rate-and-state model, on the induced seismicity of the Groningen gas field in the Netherlands. Stress changes are calculated based on a model of reservoir compaction since the onset of gas production. We next compare the seismicity predicted by our threshold model and Dieterich’s model with the observations. The two models yields comparable spatial distributions of earthquakes in good agreement with the observations. We find however that the Threshold model provides a better fit to the observed time-varying seismicity rate than Dieterich’s model, and reproduces better the onset, peak and decline of the observed seismicity rate. We compute the maximum magnitude expected for each model given the Gutenberg–Richter distribution and compare to the observations. We find that the Threshold model both shows better agreement with the observed maximum magnitude and provides result consistent with lack of observed seismicity prior to 1993. We carry out analysis of the model fit using a Chi-squared reduced statistics and find that the model fit is dramatically improved by smoothing the seismicity rate. We interpret this finding as possibly suggesting an influence of source interactions, or clustering, on a long timescale of about 3–5 yr.