TY - JOUR A1 - Bayona Viveros, Jose Antonio A1 - von Specht, Sebastian A1 - Strader, Anne A1 - Hainzl, Sebastian A1 - Cotton, Fabrice Pierre A1 - Schorlemmer, Danijel T1 - A Regionalized Seismicity Model for Subduction Zones Based on Geodetic Strain Rates, Geomechanical Parameters, and Earthquake-Catalog Data JF - Bulletin of the Seismological Society of America N2 - The Seismic Hazard Inferred from Tectonics based on the Global Strain Rate Map (SHIFT_GSRM) earthquake forecast was designed to provide high-resolution estimates of global shallow seismicity to be used in seismic hazard assessment. This model combines geodetic strain rates with global earthquake parameters to characterize long-term rates of seismic moment and earthquake activity. Although SHIFT_GSRM properly computes seismicity rates in seismically active continental regions, it underestimates earthquake rates in subduction zones by an average factor of approximately 3. We present a complementary method to SHIFT_GSRM to more accurately forecast earthquake rates in 37 subduction segments, based on the conservation of moment principle and the use of regional interface seismicity parameters, such as subduction dip angles, corner magnitudes, and coupled seismogenic thicknesses. In seven progressive steps, we find that SHIFT_GSRM earthquake-rate underpredictions are mainly due to the utilization of a global probability function of seismic moment release that poorly captures the great variability among subduction megathrust interfaces. Retrospective test results show that the forecast is consistent with the observations during the 1 January 1977 to 31 December 2014 period. Moreover, successful pseudoprospective evaluations for the 1 January 2015 to 31 December 2018 period demonstrate the power of the regionalized earthquake model to properly estimate subduction-zone seismicity. Y1 - 2019 U6 - https://doi.org/10.1785/0120190034 SN - 0037-1106 SN - 1943-3573 VL - 109 IS - 5 SP - 2036 EP - 2049 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Kühn, Daniela A1 - Hainzl, Sebastian A1 - Dahm, Torsten A1 - Richter, Gudrun A1 - Vera Rodriguez, Ismael T1 - A review of source models to further the understanding of the seismicity of the Groningen field JF - Netherlands journal of geosciences : NJG N2 - The occurrence of felt earthquakes due to gas production in Groningen has initiated numerous studies and model attempts to understand and quantify induced seismicity in this region. The whole bandwidth of available models spans the range from fully deterministic models to purely empirical and stochastic models. In this article, we summarise the most important model approaches, describing their main achievements and limitations. In addition, we discuss remaining open questions and potential future directions of development. KW - deterministic KW - empirical KW - hybrid KW - machine learning KW - seismicity model Y1 - 2022 U6 - https://doi.org/10.1017/njg.2022.7 SN - 0016-7746 SN - 1573-9708 VL - 101 PB - Cambridge Univ. Press CY - Cambridge ER - TY - JOUR A1 - Zoller, Gert A1 - Hainzl, Sebastian A1 - Holschneider, Matthias A1 - Ben-Zion, Yehuda T1 - Aftershocks resulting from creeping sections in a heterogeneous fault N2 - We show that realistic aftershock sequences with space-time characteristics compatible with observations are generated by a model consisting of brittle fault segments separated by creeping zones. The dynamics of the brittle regions is governed by static/kinetic friction, 3D elastic stress transfer and small creep deformation. The creeping parts are characterized by high ongoing creep velocities. These regions store stress during earthquake failures and then release it in the interseismic periods. The resulting postseismic deformation leads to aftershock sequences following the modified Omori law. The ratio of creep coefficients in the brittle and creeping sections determines the duration of the postseismic transients and the exponent p of the modified Omori law Y1 - 2005 SN - 0094-8276 ER - TY - JOUR A1 - Passarelli, Luigi A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Maccaferri, Francesco A1 - Mucciarelli, Marco A1 - Rößler, Dirk A1 - Corbi, Fabio A1 - Dahm, Torsten A1 - Rivalta, Eleonora T1 - Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy JF - Geophysical journal international N2 - Tectonic earthquake swarms challenge our understanding of earthquake processes since it is difficult to link observations to the underlying physical mechanisms and to assess the hazard they pose. Transient forcing is thought to initiate and drive the spatio-temporal release of energy during swarms. The nature of the transient forcing may vary across sequences and range from aseismic creeping or transient slip to diffusion of pore pressure pulses to fluid redistribution and migration within the seismogenic crust. Distinguishing between such forcing mechanisms may be critical to reduce epistemic uncertainties in the assessment of hazard due to seismic swarms, because it can provide information on the frequency-magnitude distribution of the earthquakes (often deviating from the assumed Gutenberg-Richter relation) and on the expected source parameters influencing the ground motion (for example the stress drop). Here we study the ongoing Pollino range (Southern Italy) seismic swarm, a long-lasting seismic sequence with more than five thousand events recorded and located since October 2010. The two largest shocks (magnitude M-w = 4.2 and M-w = 5.1) are among the largest earthquakes ever recorded in an area which represents a seismic gap in the Italian historical earthquake catalogue. We investigate the geometrical, mechanical and statistical characteristics of the largest earthquakes and of the entire swarm. We calculate the focal mechanisms of the M-l > 3 events in the sequence and the transfer of Coulomb stress on nearby known faults and analyse the statistics of the earthquake catalogue. We find that only 25 per cent of the earthquakes in the sequence can be explained as aftershocks, and the remaining 75 per cent may be attributed to a transient forcing. The b-values change in time throughout the sequence, with low b-values correlated with the period of highest rate of activity and with the occurrence of the largest shock. In the light of recent studies on the palaeoseismic and historical activity in the Pollino area, we identify two scenarios consistent with the observations and our analysis: This and past seismic swarms may have been 'passive' features, with small fault patches failing on largely locked faults, or may have been accompanied by an 'active', largely aseismic, release of a large portion of the accumulated tectonic strain. Those scenarios have very different implications for the seismic hazard of the area. KW - Seismicity and tectonics KW - Statistical seismology KW - Dynamics: seismotectonics Y1 - 2015 U6 - https://doi.org/10.1093/gji/ggv111 SN - 0956-540X SN - 1365-246X VL - 201 IS - 3 SP - 1553 EP - 1567 PB - Oxford Univ. Press CY - Oxford ER - TY - GEN A1 - Passarelli, Luigi A1 - Hainzl, Sebastian A1 - Cesca, Simone A1 - Maccaferri, Francesco A1 - Mucciarelli, Marco A1 - Roessler, Dirk A1 - Corbi, Fabio A1 - Dahm, Torsten A1 - Rivalta, Eleonora T1 - Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy (vol 201, pg 1553, 2015) T2 - Geophysical journal international KW - Seismicity and tectonics KW - Statistical seismology KW - Dynamics: seismotectonics Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggv425 SN - 0956-540X SN - 1365-246X VL - 204 SP - 365 EP - 365 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Wang, Lifeng A1 - Hainzl, Sebastian A1 - Zöller, Gert T1 - Assessment of stress coupling among the inter-, co- and post-seismic phases related to the 2004 M6 Parkfield earthquake JF - Geophysical journal international N2 - Due to large uncertainties and non-uniqueness in fault slip inversion, the investigation of stress coupling based on the direct comparison of independent slip inversions, for example, between the coseismic slip distribution and the interseismic slip deficit, may lead to ambiguous conclusions. In this study, we therefore adopt the stress-constrained joint inversion in the Bayesian approach of Wang et al., and implement the physical hypothesis of stress coupling as a prior. We test the hypothesis that interseismic locking is coupled with the coseismic rupture, and the early post-seismic deformation is a stress relaxation process in response to the coseismic stress perturbation. We characterize the role of stress coupling in the seismic cycle by evaluating the efficiency of the model to explain the available data. Taking the 2004 M6 Parkfield earthquake as a study case, we find that the stress coupling hypothesis is in agreement with the data. The coseismic rupture zone is found to be strongly locked during the interseismic phase and the post-seismic slip zone is indicated to be weakly creeping. The post-seismic deformation plays an important role to rebuild stress in the coseismic rupture zone. Based on our results for the stress accumulation during both inter- and post-seismic phase in the coseismic rupture zone, together with the coseismic stress drop, we estimate a recurrence time of M6 earthquake in Parkfield around 23-41 yr, suggesting that the duration of 38 yr between the two recent M6 events in Parkfield is not a surprise. KW - Inverse theory KW - Seismic cycle KW - Dynamics and mechanics of faulting Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggu102 SN - 0956-540X SN - 1365-246X VL - 197 IS - 3 SP - 1858 EP - 1868 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Daniel, G. A1 - Prono, E. A1 - Renard, F. A1 - Thouvenot, F. A1 - Hainzl, Sebastian A1 - Marsan, D. A1 - Helmstetter, A. A1 - Traversa, P. A1 - Got, J. L. A1 - Jenatton, L. A1 - Guiguet, R. T1 - Changes in effective stress during the 2003-2004 Ubaye seismic swarm, France JF - Journal of geophysical research : Solid earth N2 - We study changes in effective stress (normal stress minus pore pressure) that occurred in the French Alps during the 2003-2004 Ubaye earthquake swarm. Two complementary data sets are used. First, a set of 974 relocated events allows us to finely characterize the shape of the seismogenic area and the spatial migration of seismicity during the crisis. Relocations are performed by a double-difference algorithm. We compute differences in travel times at stations both from absolute picking times and from cross-correlation delays of multiplets. The resulting catalog reveals a swarm alignment along a single planar structure striking N130 degrees E and dipping 80 degrees W. This relocated activity displays migration properties consistent with a triggering by a diffusive fluid overpressure front. This observation argues in favor of a deep-seated fluid circulation responsible for a significant part of the seismic activity in Ubaye. Second, we analyze time series of earthquake detections at a single seismological station located just above the swarm. This time series forms a dense chronicle of +16,000 events. We use it to estimate the history of effective stress changes during this sequence. For this purpose we model the rate of events by a stochastic epidemic-type aftershock sequence model with a nonstationary background seismic rate lambda(0)(t). This background rate is estimated in discrete time windows. Window lengths are determined optimally according to a new change-point method on the basis of the interevent times distribution. We propose that background events are triggered directly by a transient fluid circulation at depth. Then, using rate-and-state constitutive friction laws, we estimate changes in effective stress for the observed rate of background events. We assume that changes in effective stress occurred under constant shear stressing rate conditions. We finally obtain a maximum change in effective stress close to -8 MPa, which corresponds to a maximum fluid overpressure of about 8 MPa under constant normal stress conditions. This estimate is in good agreement with values obtained from numerical modeling of fluid flow at depth, or with direct measurements reported from fluid injection experiments. Y1 - 2011 U6 - https://doi.org/10.1029/2010JB007551 SN - 2169-9313 SN - 2169-9356 VL - 116 IS - 4 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Zöller, Gert A1 - Hainzl, Sebastian A1 - Tilmann, Frederik A1 - Woith, Heiko A1 - Dahm, Torsten T1 - Comment on: Wikelski, Martin; Müller, Uschi; Scocco, Paola; Catorci, Andrea; Desinov, Lev V.; Belyaev, Mikhail Y.; Keim, Daniel A.; Pohlmeier, Winfried; Fechteler, Gerhard; Mai, Martin P. : Potential short-term earthquake forecasting by farm animal monitoring. - Ethology. - 126 (2020), 9. - S. 931 - 941. -ISSN 0179-1613. - eISSN 1439-0310. - doi 10.1111/eth.13078 JF - Ethology N2 - Based on an analysis of continuous monitoring of farm animal behavior in the region of the 2016 M6.6 Norcia earthquake in Italy, Wikelski et al., 2020; (Seismol Res Lett, 89, 2020, 1238) conclude that animal activity can be anticipated with subsequent seismic activity and that this finding might help to design a "short-term earthquake forecasting method." We show that this result is based on an incomplete analysis and misleading interpretations. Applying state-of-the-art methods of statistics, we demonstrate that the proposed anticipatory patterns cannot be distinguished from random patterns, and consequently, the observed anomalies in animal activity do not have any forecasting power. KW - animal behavior KW - earthquake precursor KW - error diagram KW - prediction KW - randomness KW - statistics Y1 - 2020 U6 - https://doi.org/10.1111/eth.13105 SN - 0179-1613 SN - 1439-0310 VL - 127 IS - 3 SP - 302 EP - 306 PB - Wiley CY - Hoboken ER - TY - JOUR A1 - Hainzl, Sebastian A1 - Zöller, Gert A1 - Brietzke, Gilbert B. A1 - Hinzen, Klaus-G. T1 - Comparison of deterministic and stochastic earthquake simulators for fault interactions in the Lower Rhine Embayment, Germany JF - Geophysical journal international N2 - Time-dependent probabilistic seismic hazard assessment requires a stochastic description of earthquake occurrences. While short-term seismicity models are well-constrained by observations, the recurrences of characteristic on-fault earthquakes are only derived from theoretical considerations, uncertain palaeo-events or proxy data. Despite the involved uncertainties and complexity, simple statistical models for a quasi-period recurrence of on-fault events are implemented in seismic hazard assessments. To test the applicability of statistical models, such as the Brownian relaxation oscillator or the stress release model, we perform a systematic comparison with deterministic simulations based on rate- and state-dependent friction, high-resolution representations of fault systems and quasi-dynamic rupture propagation. For the specific fault network of the Lower Rhine Embayment, Germany, we run both stochastic and deterministic model simulations based on the same fault geometries and stress interactions. Our results indicate that the stochastic simulators are able to reproduce the first-order characteristics of the major earthquakes on isolated faults as well as for coupled faults with moderate stress interactions. However, we find that all tested statistical models fail to reproduce the characteristics of strongly coupled faults, because multisegment rupturing resulting from a spatiotemporally correlated stress field is underestimated in the stochastic simulators. Our results suggest that stochastic models have to be extended by multirupture probability distributions to provide more reliable results. KW - Earthquake interaction KW - forecasting KW - and prediction KW - Seismicity and tectonics KW - Statistical seismology Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt271 SN - 0956-540X SN - 1365-246X VL - 195 IS - 1 SP - 684 EP - 694 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Hainzl, Sebastian A1 - Ogata, Y. T1 - Detecting fluid signals in seismicity data through statistical earthquake modeling N2 - [1] According to the well-known Coulomb failure criterion the variation of either stress or pore pressure can result in earthquake rupture. Aftershock sequences characterized by the Omori law are often assumed to be the consequence of varying stress, whereas earthquake swarms are thought to be triggered by fluid intrusions. The role of stress triggering can be analyzed by modeling solely three-dimensional (3-D) elastic stress changes in the crust, but fluid flows which initiate seismicity cannot be investigated without considering complex seismicity patterns resulting from both pore pressure variations and earthquake-connected stress field changes. We show that the epidemic-type aftershock sequence (ETAS) model is an appropriate tool to extract the primary fluid signal from such complex seismicity patterns. We analyze a large earthquake swarm that occurred in 2000 in Vogtland/NW Bohemia, central Europe. By fitting the stochastic ETAS model, we find that stress triggering is dominant in creating the observed seismicity patterns and explains the observed fractal interevent time distribution. External forcing, identified with pore pressure changes due to fluid intrusion, is found to directly trigger only a few percent of the total activity. However, temporal deconvolution indicates that a pronounced fluid signal initiated the swarm. These results are confirmed by our analogous investigation of model simulations in which earthquakes are triggered by fluid intrusion as well as stress transfers on a fault plane embedded in a 3-D elastic half-space. The deconvolution procedure based on the ETAS model is able to reveal the underlying pore pressure variations Y1 - 2005 SN - 0148-0227 ER -