@article{PirliHainzlSchweitzeretal.2018, author = {Pirli, Myrto and Hainzl, Sebastian and Schweitzer, Johannes and K{\"o}hler, Andreas and Dahm, Torsten}, title = {Localised thickening and grounding of an Antarctic ice shelf from tidal triggering and sizing of cryoseismicity}, series = {Earth \& planetary science letters}, volume = {503}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2018.09.024}, pages = {78 -- 87}, year = {2018}, abstract = {We observe remarkably periodic patterns of seismicity rates and magnitudes at the Fimbul Ice Shelf, East Antarctica, correlating with the cycles of the ocean tide. Our analysis covers 19 years of continuous seismic recordings from Antarctic broadband stations. Seismicity commences abruptly during austral summer 2011 at a location near the ocean front in a shallow water region. Dozens of highly repetitive events occur in semi-diurnal cycles, with magnitudes and rates fluctuating steadily with the tide. In contrast to the common unpredictability of earthquake magnitudes, the event magnitudes show deterministic trends within single cycles and strong correlations with spring tides and tide height. The events occur quasi-periodically and the highly constrained event sources migrate landwards during rising tide. We show that a simple, mechanical model can explain most of the observations. Our model assumes stick-slip motion on a patch of grounded ice shelf, which is forced by the variations of the ocean-tide height and ice flow. The well fitted observations give new insights into the general process of frictional triggering of earthquakes, while providing independent evidence of variations in ice shelf thickness and grounding.}, language = {en} } @article{PassarelliHainzlCescaetal.2015, author = {Passarelli, Luigi and Hainzl, Sebastian and Cesca, Simone and Maccaferri, Francesco and Mucciarelli, Marco and R{\"o}ßler, Dirk and Corbi, Fabio and Dahm, Torsten and Rivalta, Eleonora}, title = {Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy}, series = {Geophysical journal international}, volume = {201}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggv111}, pages = {1553 -- 1567}, year = {2015}, abstract = {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.}, language = {en} } @misc{PassarelliHainzlCescaetal.2016, author = {Passarelli, Luigi and Hainzl, Sebastian and Cesca, Simone and Maccaferri, Francesco and Mucciarelli, Marco and Roessler, Dirk and Corbi, Fabio and Dahm, Torsten and Rivalta, Eleonora}, title = {Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy (vol 201, pg 1553, 2015)}, series = {Geophysical journal international}, volume = {204}, journal = {Geophysical journal international}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggv425}, pages = {365 -- 365}, year = {2016}, language = {en} } @article{NooshiriSaulHeimannetal.2017, author = {Nooshiri, Nima and Saul, Joachim and Heimann, Sebastian and Tilmann, Frederik and Dahm, Torsten}, title = {Revision of earthquake hypocentre locations in global bulletin data sets using source-specific station terms}, series = {Geophysical journal international}, volume = {208}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggw405}, pages = {589 -- 602}, year = {2017}, abstract = {Global earthquake locations are often associated with very large systematic travel-time residuals even for clear arrivals, especially for regional and near-regional stations in subduction zones because of their strongly heterogeneous velocity structure. Travel-time corrections can drastically reduce travel-time residuals at regional stations and, in consequence, improve the relative location accuracy. We have extended the shrinking-box source-specific station terms technique to regional and teleseismic distances and adopted the algorithm for probabilistic, nonlinear, global-search location. We evaluated the potential of the method to compute precise relative hypocentre locations on a global scale. The method has been applied to two specific test regions using existing P- and pP-phase picks. The first data set consists of 3103 events along the Chilean margin and the second one comprises 1680 earthquakes in the Tonga-Fiji subduction zone. Pick data were obtained from the GEOFON earthquake bulletin, produced using data from all available, global station networks. A set of timing corrections varying as a function of source position was calculated for each seismic station. In this way, we could correct the systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure without explicitly solving for a velocity model. Residual statistics show that the median absolute deviation of the travel-time residuals is reduced by 40-60 per cent at regional distances, where the velocity anomalies are strong. Moreover, the spread of the travel-time residuals decreased by similar to 20 per cent at teleseismic distances (>28 degrees). Furthermore, strong variations in initial residuals as a function of recording distance are smoothed out in the final residuals. The relocated catalogues exhibit less scattered locations in depth and sharper images of the seismicity associated with the subducting slabs. Comparison with a high-resolution local catalogue reveals that our relocation process significantly improves the hypocentre locations compared to standard locations.}, language = {en} } @article{NiemzDahmMilkereitetal.2021, author = {Niemz, Peter and Dahm, Torsten and Milkereit, Claus and Cesca, Simone and Petersen, Gesa Maria and Zang, Arno}, title = {Insights into hydraulic fracture growth gained from a joint analysis of seismometer-derived tilt signals and scoustic emissions}, series = {Journal of geophysical research : Solid earth}, volume = {126}, journal = {Journal of geophysical research : Solid earth}, number = {12}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2021JB023057}, pages = {14}, year = {2021}, abstract = {Hydraulic fracturing is performed to enhance rock permeability, for example, in the frame of geothermal energy production or shale gas exploitation, and can potentially trigger induced seismicity. The tracking of increased permeabilities and the fracturing extent is often based on the microseismic event distribution within the stimulated rock volume, but it is debated whether the microseismic activity adequately depicts the fracture formation. We are able to record tilt signals that appear as long-period transients (<180 s) on two broadband seismometers installed close (17-72 m) to newly formed, meter-scale hydraulic fractures. With this observation, we can overcome the limitations of the microseismic monitoring alone and verify the fracture mapping. Our analysis for the first time combines a catalog of previously analyzed acoustic emissions ([AEs] durations of 20 ms), indirectly mapping the fractures, with unique tilt signals, that provide independent, direct insights into the deformation of the rock. The analysis allows to identify different phases of the fracturing process including the (re)opening, growth, and aftergrowth of fractures. Further, it helps to differentiate between the formation of complex fracture networks and single macrofractures, and it validates the AE fracture mapping. Our findings contribute to a better understanding of the fracturing processes, which may help to reduce fluid-injection-induced seismicity and validate efficient fracture formation.
Plain Language Summary Hydraulic fracturing (HF) describes the opening of fractures in rocks by injecting fluids under high pressure. The new fractures not only can facilitate the extraction of shale gas but can also be used to heat up water in the subsurface in enhanced geothermal systems, a corner stone of renewable energy production. The fracture formation is inherently accompanied by small, nonfelt earthquakes (microseismic events). Occasionally, larger events felt by the population can be induced by the subsurface operations. Avoiding such events is important for the acceptance of HF operations and requires a detailed knowledge about the fracture formation. We jointly analyze two very different data sets recorded during mine-scale HF experiments: (a) the tilting of the ground caused by the opening of the fractures, as recorded by broadband seismometers-usually deployed for earthquake monitoring-installed close to the experiments and (b) a catalog of acoustic emissions, seismic signals of few milliseconds emitted by tiny cracks around the forming hydraulic fracture. The novel joint analysis allows to characterize the fracturing processes in greater detail, contributing to the understanding of the physical processes, which may help to understand fluid-injection-induced seismicity and validate the formation of hydraulic fractures.}, language = {en} } @article{NiemzCescaHeimannetal.2020, author = {Niemz, Peter and Cesca, Simone and Heimann, Sebastian and Grigoli, Francesco and von Specht, Sebastian and Hammer, Conny and Zang, Arno and Dahm, Torsten}, title = {Full-waveform-based characterization of acoustic emission activity in a mine-scale experiment}, series = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, volume = {222}, journal = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0955-419X}, doi = {10.1093/gji/ggaa127}, pages = {189 -- 206}, year = {2020}, abstract = {Understanding fracturing processes and the hydromechanical relation to induced seismicity is a key question for enhanced geothermal systems (EGS). Commonly massive fluid injection, predominately causing hydroshearing, are used in large-scale EGS but also hydraulic fracturing approaches were discussed. To evaluate the applicability of hydraulic fracturing techniques in EGS, six in situ, multistage hydraulic fracturing experiments with three different injection schemes were performed under controlled conditions in crystalline rock at the Aspo Hard Rock Laboratory (Sweden). During the experiments the near-field ground motion was continuously recorded by 11 piezoelectric borehole sensors with a sampling rate of 1 MHz. The sensor network covered a volume of 30x30x30 m around a horizontal, 28-m-long injection borehole at a depth of 410 m. To extract and characterize massive, induced, high-frequency acoustic emission (AE) activity from continuous recordings, a semi-automated workflow was developed relying on full waveform based detection, classification and location procedures. The approach extended the AE catalogue from 196 triggered events in previous studies to more than 19600 located AEs. The enhanced catalogue, for the first time, allows a detailed analysis of induced seismicity during single hydraulic fracturing experiments, including the individual fracturing stages and the comparison between injection schemes. Beside the detailed study of the spatio-temporal patterns, event clusters and the growth of seismic clouds, we estimate relative magnitudes and b-values of AEs for conventional, cyclic progressive and dynamic pulse injection schemes, the latter two being fatigue hydraulic fracturing techniques. While the conventional fracturing leads to AE patterns clustered in planar regions, indicating the generation of a single main fracture plane, the cyclic progressive injection scheme results in a more diffuse, cloud-like AE distribution, indicating the activation of a more complex fracture network. For a given amount of hydraulic energy (pressure multiplied by injected volume) pumped into the system, the cyclic progressive scheme is characterized by a lower rate of seismicity, lower maximum magnitudes and significantly larger b-values, implying an increased number of small events relative to the large ones. To our knowledge, this is the first direct comparison of high resolution seismicity in a mine-scale experiment induced by different hydraulic fracturing schemes.}, language = {en} } @article{MaghsoudiHainzlCescaetal.2014, author = {Maghsoudi, Samira and Hainzl, Sebastian and Cesca, Simone and Dahm, Torsten and Kaiser, Diethelm}, title = {Identification and characterization of growing large-scale en-echelon fractures in a salt mine}, series = {Geophysical journal international}, volume = {196}, journal = {Geophysical journal international}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggt443}, pages = {1092 -- 1105}, year = {2014}, abstract = {The spatiotemporal seismicity of acoustic emission (AE) events recorded in the Morsleben salt mine is investigated. Almost a year after backfilling of the cavities from 2003, microevents are distributed with distinctive stripe shapes above cavities at different depth levels. The physical forces driving the creation of these stripes are still unknown. This study aims to find the active stripes and track fracture developments over time by combining two different temporal and spatial clustering techniques into a single methodological approach. Anomalous seismicity parameters values like sharp b-value changes for two active stripes are good indicators to explain possible stress accumulation at the stripe tips. We identify the formation of two new seismicity stripes and show that the AE activities in active clusters are migrated mostly unidirectional to eastward and upward. This indicates that the growth of underlying macrofractures is controlled by the gradient of extensional stress. Studying size distribution characteristic in terms of frequency-magnitude distribution and b-value in active phase and phase with constant seismicity rate show that deviations from the Gutenberg-Richter power law can be explained by the inclusion of different activity phases: (1) the inactive period before the formation of macrofractures, which is characterized by a deficit of larger events (higher b-values) and (2) the period of fracture growth characterized by the occurrence of larger events (smaller b-values).}, language = {en} } @article{MaghsoudiCescaHainzletal.2013, author = {Maghsoudi, Samira and Cesca, Simone and Hainzl, Sebastian and Kaiser, Diethelm and Becker, Dirk and Dahm, Torsten}, title = {Improving the estimation of detection probability and magnitude of completeness in strongly heterogeneous media, an application to acoustic emission (AE)}, series = {Geophysical journal international}, volume = {193}, journal = {Geophysical journal international}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggt049}, pages = {1556 -- 1569}, year = {2013}, abstract = {Reliable estimations of magnitude of completeness (M-c) are essential for a correct interpretation of seismic catalogues. The spatial distribution of M-c may be strongly variable and difficult to assess in mining environments, owing to the presence of galleries, cavities, fractured regions, porous media and different mineralogical bodies, as well as in consequence of inhomogeneous spatial distribution of the seismicity. We apply a 3-D modification of the probabilistic magnitude of completeness (PMC) method, which relies on the analysis of network detection capabilities. In our approach, the probability to detect an event depends on its magnitude, source receiver Euclidian distance and source receiver direction. The suggested method is proposed for study of the spatial distribution of the magnitude of completeness in a mining environment and here is applied to a 2-months acoustic emission (AE) data set recorded at the Morsleben salt mine, Germany. The dense seismic network and the large data set, which includes more than one million events, enable a detailed testing of the method. This method is proposed specifically for strongly heterogeneous media. Besides, it can also be used for specific network installations, with sensors with a sensitivity, dependent on the direction of the incoming wave (e.g. some piezoelectric sensors). In absence of strong heterogeneities, the standards PMC approach should be used. We show that the PMC estimations in mines strongly depend on the source receiver direction, and cannot be correctly accounted using a standard PMC approach. However, results can be improved, when adopting the proposed 3-D modification of the PMC method. Our analysis of one central horizontal and vertical section yields a magnitude of completeness of about M-c approximate to 1 (AE magnitude) at the centre of the network, which increases up to M-c approximate to 4 at further distances outside the network; the best detection performance is estimated for a NNE-SSE elongated region, which corresponds to the strike direction of the low-attenuating salt body. Our approach provides us with small-scale details about the capability of sensors to detect an earthquake, which can be linked to the presence of heterogeneities in specific directions. Reduced detection performance in presence of strong structural heterogeneities (cavities) is confirmed by synthetic waveform modelling in heterogeneous media.}, language = {en} } @article{MaghsoudiCescaHainzletal.2015, author = {Maghsoudi, Samira and Cesca, Simone and Hainzl, Sebastian and Dahm, Torsten and Z{\"o}ller, Gert and Kaiser, Diethelm}, title = {Maximum Magnitude of Completeness in a Salt Mine}, series = {Bulletin of the Seismological Society of America}, volume = {105}, journal = {Bulletin of the Seismological Society of America}, number = {3}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120140039}, pages = {1491 -- 1501}, year = {2015}, abstract = {In this study, we analyze acoustic emission (AE) data recorded at the Morsleben salt mine, Germany, to assess the catalog completeness, which plays an important role in any seismicity analysis. We introduce the new concept of a magnitude completeness interval consisting of a maximum magnitude of completeness (M-c(max)) in addition to the well-known minimum magnitude of completeness. This is required to describe the completeness of the catalog, both for the smallest events (for which the detection performance may be low) and for the largest ones (which may be missed because of sensors saturation). We suggest a method to compute the maximum magnitude of completeness and calculate it for a spatial grid based on (1) the prior estimation of saturation magnitude at each sensor, (2) the correction of the detection probability function at each sensor, including a drop in the detection performance when it saturates, and (3) the combination of detection probabilities of all sensors to obtain the network detection performance. The method is tested using about 130,000 AE events recorded in a period of five weeks, with sources confined within a small depth interval, and an example of the spatial distribution of M-c(max) is derived. The comparison between the spatial distribution of M-c(max) and of the maximum possible magnitude (M-max), which is here derived using a recently introduced Bayesian approach, indicates that M-max exceeds M-c(max) in some parts of the mine. This suggests that some large and important events may be missed in the catalog, which could lead to a bias in the hazard evaluation.}, language = {en} } @article{KuehnHainzlDahmetal.2022, author = {K{\"u}hn, Daniela and Hainzl, Sebastian and Dahm, Torsten and Richter, Gudrun and Vera Rodriguez, Ismael}, title = {A review of source models to further the understanding of the seismicity of the Groningen field}, series = {Netherlands journal of geosciences : NJG}, volume = {101}, journal = {Netherlands journal of geosciences : NJG}, publisher = {Cambridge Univ. Press}, address = {Cambridge}, issn = {0016-7746}, doi = {10.1017/njg.2022.7}, pages = {12}, year = {2022}, abstract = {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.}, language = {en} }