TY - JOUR A1 - Heimann, Sebastian A1 - Vasyura-Bathke, Hannes A1 - Sudhaus, Henriette A1 - Isken, Marius Paul A1 - Kriegerowski, Marius A1 - Steinberg, Andreas A1 - Dahm, Torsten T1 - A Python framework for efficient use of pre-computed Green's functions in seismological and other physical forward and inverse source problems JF - Solid earth N2 - The computation of such synthetic GFs is computationally and operationally demanding. As a consequence, the onthe-fly recalculation of synthetic GFs in each iteration of an optimisation is time-consuming and impractical. Therefore, the pre-calculation and efficient storage of synthetic GFs on a dense grid of source to receiver combinations enables the efficient lookup and utilisation of GFs in time-critical scenarios. We present a Python-based framework and toolkit - Pyrocko-GF - that enables the pre-calculation of synthetic GF stores, which are independent of their numerical calculation method and GF transfer function. The framework aids in the creation of such GF stores by interfacing a suite of established numerical forward modelling codes in seismology (computational back ends). So far, interfaces to back ends for layered Earth model cases have been provided; however, the architecture of Pyrocko-GF is designed to cover back ends for other geometries (e.g. full 3-D heterogeneous media) and other physical quantities (e.g. gravity, pressure, tilt). Therefore, Pyrocko-GF defines an extensible GF storage format suitable for a wide range of GF types, especially handling elasticity and wave propagation problems. The framework assists with visualisations, quality control, and the exchange of GF stores, which is supported through an online platform that provides many pre-calculated GF stores for local, regional, and global studies. The Pyrocko-GF toolkit comes with a well-documented application programming interface (API) for the Python programming language to efficiently facilitate forward modelling of geophysical processes, e.g. synthetic waveforms or static displacements for a wide range of source models. Y1 - 2019 U6 - https://doi.org/10.5194/se-10-1921-2019 SN - 1869-9510 SN - 1869-9529 VL - 10 IS - 6 SP - 1921 EP - 1935 PB - Copernicus CY - Göttingen 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 - Karamzadeh Toularoud, Nasim A1 - Heimann, Sebastian A1 - Dahm, Torsten A1 - Krüger, Frank T1 - Application based seismological array design by seismicity scenario modelling JF - Geophysical journal international N2 - The design of an array configuration is an important task in array seismology during experiment planning. Often the array response function (ARF), which depends on the relative position of array stations and frequency content of the incoming signals, is used as the array design criterion. In practice, additional constraints and parameters have to be taken into account, for example, land ownership, site-specific noise levels or characteristics of the seismic sources under investigation. In this study, a flexible array design framework is introduced that implements a customizable scenario modelling and optimization scheme by making use of synthetic seismograms. Using synthetic seismograms to evaluate array performance makes it possible to consider additional constraints. We suggest to use synthetic array beamforming as an array design criterion instead of the ARF. The objective function of the optimization scheme is defined according to the monitoring goals, and may consist of a number of subfunctions. The array design framework is exemplified by designing a seven-station small-scale array to monitor earthquake swarm activity in Northwest Bohemia/Vogtland in central Europe. Two subfunctions are introduced to verify the accuracy of horizontal slowness estimation; one to suppress aliasing effects due to possible secondary lobes of synthetic array beamforming calculated in horizontal slowness space and the other to reduce the event’s mislocation caused by miscalculation of the horizontal slowness vector. Subsequently, a weighting technique is applied to combine the subfunctions into one single scalar objective function to use in the optimization process. KW - Array Seismology KW - Array design KW - Seismicity modelling Y1 - 2018 SN - 0956-540X SN - 1365-246X VL - 216 IS - 3 SP - 1711 EP - 1727 PB - Oxford Univ. Press CY - Oxford 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 - Sen, Ali Tolga A1 - Cesca, Simone A1 - Bischoff, Monika A1 - Meier, Thomas A1 - Dahm, Torsten T1 - Automated full moment tensor inversion of coal mining-induced seismicity JF - Geophysical journal international N2 - Seismicity induced by coal mining in the Ruhr region, Germany, has been monitored continuously over the last 25 yr. In 2006, a dense temporary network (HAMNET) was deployed to locally monitor seismicity induced by longwall mining close to the town of Hamm. Between 2006 July and 2007 July, more than 7000 events with magnitudes M-L from -1.7 to 2.0 were detected. The spatiotemporal distribution of seismicity shows high correlation with the mining activity. In order to monitor rupture processes, we set up an automated source inversion routine and successfully perform double couple and full moment tensor (MT) inversions for more than 1000 events with magnitudes above M-L -0.5. The source inversion is based on a full waveform approach, both in the frequency and in the time domain, providing information about the centroid location, focal mechanism, scalar moment and full MT. Inversion results indicate a strong dominance of normal faulting focal mechanisms, with a steeper plane and a subhorizontal one. Fault planes are oriented parallel to the mining stopes. We classify the focal mechanisms based on their orientation and observe different frequency-magnitude distributions for families of events with different focal mechanisms; the overall frequency-magnitude distribution is not fitting the Gutenberg-Richter relation. Full MTs indicate that non-negligible opening tensile components accompanied normal faulting source mechanisms. Finally, extended source models are investigated for largest events. Results suggest that the rupture processes mostly occurred along the subvertical planes. KW - Geomechanics KW - Fracture and flow KW - Earthquake source observations KW - Seismicity and tectonics Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt300 SN - 0956-540X SN - 1365-246X VL - 195 IS - 2 SP - 1267 EP - 1281 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Vassallo, Maurizio A1 - Dahm, Torsten T1 - Automated seismic event location by travel-time stacking an application to mining induced seismicity JF - Seismological research letters Y1 - 2013 U6 - https://doi.org/10.1785/0220120191 SN - 0895-0695 VL - 84 IS - 4 SP - 666 EP - 677 PB - Seismological Society of America CY - Albany ER - TY - JOUR A1 - Grigoli, Francesco A1 - Cesca, Simone A1 - Amoroso, Ortensia A1 - Emolo, Antonio A1 - Zollo, Aldo A1 - Dahm, Torsten T1 - Automated seismic event location by waveform coherence analysis JF - Geophysical journal international N2 - Automated location of seismic events is a very important task in microseismic monitoring operations as well for local and regional seismic monitoring. Since microseismic records are generally characterized by low signal-to-noise ratio, automated location methods are requested to be noise robust and sufficiently accurate. Most of the standard automated location routines are based on the automated picking, identification and association of the first arrivals of P and S waves and on the minimization of the residuals between theoretical and observed arrival times of the considered seismic phases. Although current methods can accurately pick P onsets, the automatic picking of the S onset is still problematic, especially when the P coda overlaps the S wave onset. In this paper, we propose a picking free earthquake location method based on the use of the short-term-average/long-term-average (STA/LTA) traces at different stations as observed data. For the P phases, we use the STA/LTA traces of the vertical energy function, whereas for the S phases, we use the STA/LTA traces of a second characteristic function, which is obtained using the principal component analysis technique. In order to locate the seismic event, we scan the space of possible hypocentral locations and origin times, and stack the STA/LTA traces along the theoretical arrival time surface for both P and S phases. Iterating this procedure on a 3-D grid, we retrieve a multidimensional matrix whose absolute maximum corresponds to the spatial coordinates of the seismic event. A pilot application was performed in the Campania-Lucania region (southern Italy) using a seismic network (Irpinia Seismic Network) with an aperture of about 150 km. We located 196 crustal earthquakes (depth < 20 km) with magnitude range 1.1 < M-L < 2.7. A subset of these locations were compared with accurate manual locations refined by using a double-difference technique. Our results indicate a good agreement with manual locations. Moreover, our method is noise robust and performs better than classical location methods based on the automatic picking of the P and S waves first arrivals. KW - Time-series analysis KW - Inverse theory KW - Earthquake source observations KW - Seismicity and tectonics KW - Early warning Y1 - 2014 U6 - https://doi.org/10.1093/gji/ggt477 SN - 0956-540X SN - 1365-246X VL - 196 IS - 3 SP - 1742 EP - 1753 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Dahm, Torsten A1 - Kuehn, Daniela A1 - Ohrnberger, Matthias A1 - Kroeger, Jens A1 - Wiederhold, Helga A1 - Reuther, Claus-Dieter A1 - Dehghani, Ali A1 - Scherbaum, Frank T1 - Combining geophysical data sets to study the dynamics of shallow evaporites in urban environments : application to Hamburg, Germany N2 - Shallowly situated evaporites in built-up areas are of relevance for urban and cultural development and hydrological regulation. The hazard of sinkholes, subrosion depressions and gypsum karst is often difficult to evaluate and may quickly change with anthropogenic influence. The geophysical exploration of evaporites in metropolitan areas is often not feasible with active industrial techniques. We collect and combine different passive geophysical data as microgravity, ambient vibrations, deformation and hydrological information to study the roof morphology of shallow evaporites beneath Hamburg, Northern Germany. The application of a novel gravity inversion technique leads to a 3-D depth model of the salt diapir under study. We compare the gravity-based depth model to pseudo-depths from H/V measurements and depth estimates from small-scale seismological array data. While the general range and trend of the diapir roof is consistent, a few anomalous regions are identified where H/V pseudo-depths indicate shallower structures not observed in gravity or array data. These are interpreted by shallow residual caprock floaters and zones of increased porosity. The shallow salt structure clearly correlates with a relative subsidence in the order of 2 mm yr(-1). The combined interpretation of roof morphology, yearly subsidence rates, chemical analyses of groundwater and of hydraulic head in aquifers indicates that the salt diapir beneath Hamburg is subject to significant ongoing dissolution that may possibly affect subrosion depressions, sinkhole distribution and land usage. The combined analysis of passive geophysical data may be exemplary for the study of shallow evaporites beneath other urban areas. Y1 - 2010 UR - http://onlinelibrary.wiley.com/journal/10.1111/(ISSN)1365-246X U6 - https://doi.org/10.1111/j.1365-246X.2010.04521.x SN - 0956-540X 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 -