@article{DahmHeimannMetzetal.2021,
  author    = {Dahm, Torsten and Heimann, Sebastian and Metz, Malte and Isken, Marius Paul},
  title     = {A self-similar dynamic rupture model based on the simplified wave-rupture analogy},
  series = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society},
  volume    = {225},
  journal   = {Geophysical journal international / the Royal Astronomical Society, the Deutsche Geophysikalische Gesellschaft and the European Geophysical Society},
  number    = {3},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggab045},
  pages     = {1586 -- 1604},
  year      = {2021},
  abstract  = {The investigation of stresses, faults, structure and seismic hazards requires a good understanding and mapping of earthquake rupture and slip. Constraining the finite source of earthquakes from seismic and geodetic waveforms is challenging because the directional effects of the rupture itself are small and dynamic numerical solutions often include a large number of free parameters. The computational effort is large and therefore difficult to use in an exploratory forward modelling or inversion approach. Here, we use a simplified self-similar fracture model with only a few parameters, where the propagation of the fracture front is decoupled from the calculation of the slip. The approximative method is flexible and computationally efficient. We discuss the strengths and limitations of the model with real-case examples of well-studied earthquakes. These include the M-w 8.3 2015 Illapel, Chile, megathrust earthquake at the plate interface of a subduction zone and examples of continental intraplate strike-slip earthquakes like the M-w 7.1 2016 Kumamoto, Japan, multisegment variable slip event or the M-w 7.5 2018 Palu, Indonesia, supershear earthquake. Despite the simplicity of the model, a large number of observational features ranging from different rupture-front isochrones and slip distributions to directional waveform effects or high slip patches are easy to model. The temporal evolution of slip rate and rise time are derived from the incremental growth of the rupture and the stress drop without imposing other constraints. The new model is fast and implemented in the open-source Python seismology toolbox Pyrocko, ready to study the physics of rupture and to be used in finite source inversions.},
  language  = {en}
}
@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{DahmKuehnOhrnbergeretal.2010,
  author    = {Dahm, Torsten and Kuehn, Daniela and Ohrnberger, Matthias and Kroeger, Jens and Wiederhold, Helga and Reuther, Claus-Dieter and Dehghani, Ali and Scherbaum, Frank},
  title     = {Combining geophysical data sets to study the dynamics of shallow evaporites in urban environments : application to Hamburg, Germany},
  issn      = {0956-540X},
  doi       = {10.1111/j.1365-246X.2010.04521.x},
  year      = {2010},
  abstract  = {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.},
  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. <br /> 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{KruegerDahmHannemann2020,
  author    = {Kr{\"u}ger, Frank and Dahm, Torsten and Hannemann, Katrin},
  title     = {Mapping of Eastern North Atlantic Ocean seismicity from Po/So observations at a mid-aperture seismological broad-band deep sea array},
  series = {Geophysical journal international},
  volume    = {221},
  journal   = {Geophysical journal international},
  number    = {2},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggaa054},
  pages     = {1055 -- 1080},
  year      = {2020},
  abstract  = {A mid-aperture broad-band test array (OBS array DOCTAR) was deployed from June 2011 to April 2012 about 100 km north of the Gloria fault in the Eastern North Atlantic in about 5000 m water depth. In addition arrays were installed on Madeira Island and in western Portugal mainland. For the first time in the Eastern North Atlantic, we recorded a large number of high frequency Po and So waves from local and regional small and moderate earthquakes (M-L < 4). An incoherent beamforming method was adapted to scan continuous data for such Po and So arrivals applying a sliding window waveform migration and frequency-wavenumber technique. We identify about 320 Po and 1550 So arrivals and compare the phase onsets with the ISC catalogue (ISC 2015) for the same time span. Up to a distance of 6 degrees to the DOCTAR stations all events listed in the ISC catalogue could be associated to Po and So phases. Arrivals from events in more than 10 degrees distance could be identified only in some cases. Only few Po and/or So arrivals were detected for earthquakes from the European and African continental area, the continental shelf regions and for earthquakes within or northwest of the Azores plateau. Unexpectedly, earthquake clusters are detected within the oceanic plates north and south of the Gloria fault and far from plate boundaries, indicating active intraplate structures. We also observe and locate numerous small magnitude earthquakes on the segment of the Gloria fault directly south of DOCTAR, which likely coincides with the rupture of the 25 November 1941 event. Local small magnitude earthquakes located beneath DOCTAR show hypocentres up to 30 km depth and strike-slip focal mechanisms. A comparison with detections at temporary mid-aperture arrays on Madeira and in western Portugal shows that the deep ocean array performs much better than the island and the continental array regarding the detection threshold for events in the oceanic plates. We conclude that sparsely distributed mid-aperture seismic arrays in the deep ocean could decrease the detection and location threshold for seismicity with M-L < 4 in the oceanic plate and might constitute a valuable tool to monitor oceanic plate seismicity.},
  language  = {en}
}
@article{KaramzadehToularoudHeimannDahmetal.2020,
  author    = {Karamzadeh Toularoud, Nasim and Heimann, Sebastian and Dahm, Torsten and Kr{\"u}ger, Frank},
  title     = {Earthquake source arrays},
  series = {Geophysical journal international},
  volume    = {221},
  journal   = {Geophysical journal international},
  number    = {1},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggaa002},
  pages     = {352 -- 370},
  year      = {2020},
  abstract  = {A collection of earthquake sources recorded at a single station, under specific conditions, are considered as a source array (SA), that is interpreted as if earthquake sources originate at the station location and are recorded at the source location. Then, array processing methods, that is array beamforming, are applicable to analyse the recorded signals. A possible application is to use source array multiple event techniques to locate and characterize near-source scatterers and structural interfaces. In this work the aim is to facilitate the use of earthquake source arrays by presenting an automatic search algorithm to configure the source array elements. We developed a procedure to search for an optimal source array element distribution given an earthquake catalogue including accurate origin time and hypocentre locations. The objective function of the optimization process can be flexibly defined for each application to ensure the prerequisites (criteria) of making a source array. We formulated four quantitative criteria as subfunctions and used the weighted sum technique to combine them in one single scalar function. The criteria are: (1) to control the accuracy of the slowness vector estimation using the time domain beamforming method, (2) to measure the waveform coherency of the array elements, (3) to select events with lower location error and (4) to select traces with high energy of specific phases, that is, sp- or ps-phases. The proposed procedure is verified using synthetic data as well as real examples for the Vogtland region in Northwest Bohemia. We discussed the possible application of the optimized source arrays to identify the location of scatterers in the velocity model by presenting a synthetic test and an example using real waveforms.},
  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{CescaGrigoliHeimannetal.2014,
  author    = {Cesca, Simone and Grigoli, Francesco and Heimann, Sebastian and Gonzalez, Alvaro and Buforn, Elisa and Maghsoudi, Samira and Blanch, Estefania and Dahm, Torsten},
  title     = {The 2013 September-October seismic sequence offshore Spain: a case of seismicity triggered by gas injection?},
  series = {Geophysical journal international},
  volume    = {198},
  journal   = {Geophysical journal international},
  number    = {2},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggu172},
  pages     = {941 -- 953},
  year      = {2014},
  abstract  = {A spatially localized seismic sequence originated few tens of kilometres offshore the Mediterranean coast of Spain, close to the Ebro river delta, starting on 2013 September 5, and lasting at least until 2013 October. The sequence culminated in a maximal moment magnitude M-w 4.3 earthquake, on 2013 October 1. The most relevant seismogenic feature in the area is the Fosa de Amposta fault system, which includes different strands mapped at different distances to the coast, with a general NE-SW orientation, roughly parallel to the coastline. However, no significant known historical seismicity has involved this fault system in the past. The epicentral region is also located near the offshore platform of the Castor project, where gas is conducted through a pipeline from mainland and where it was recently injected in a depleted oil reservoir, at about 2 km depth. We analyse the temporal evolution of the seismic sequence and use full waveform techniques to derive absolute and relative locations, estimate depths and focal mechanisms for the largest events in the sequence (with magnitude mbLg larger than 3), and compare them to a previous event (2012 April 8, mbLg 3.3) taking place in the same region prior to the gas injection. Moment tensor inversion results show that the overall seismicity in this sequence is characterized by oblique mechanisms with a normal fault component, with a 30A degrees low-dip angle plane oriented NNE-SSW and a subvertical plane oriented NW-SE. The combined analysis of hypocentral location and focal mechanisms could indicate that the seismic sequence corresponds to rupture processes along shallow low-dip surfaces, which could have been triggered by the gas injection in the reservoir, and excludes the activation of the Amposta fault, as its known orientation is inconsistent with focal mechanism results. An alternative scenario includes the iterated triggering of a system of steep faults oriented NW-SE, which were identified by prior marine seismics investigations.},
  language  = {en}
}
@article{GrigoliCescaAmorosoetal.2014,
  author    = {Grigoli, Francesco and Cesca, Simone and Amoroso, Ortensia and Emolo, Antonio and Zollo, Aldo and Dahm, Torsten},
  title     = {Automated seismic event location by waveform coherence analysis},
  series = {Geophysical journal international},
  volume    = {196},
  journal   = {Geophysical journal international},
  number    = {3},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggt477},
  pages     = {1742 -- 1753},
  year      = {2014},
  abstract  = {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.},
  language  = {en}
}
@article{KruegerDahm1999,
  author    = {Kr{\"u}ger, Frank and Dahm, Torsten},
  title     = {Higher degree moment inversin using far-field broad-band recordings : theory and evaluation of the method with application to the 1994 Bolivia deep earthauke},
  year      = {1999},
  language  = {en}
}
@article{JamalreyhaniRezapourCescaetal.2022,
  author    = {Jamalreyhani, Mohammadreza and Rezapour, Mehdi and Cesca, Simone and Dahm, Torsten and Heimann, Sebastian and Sudhaus, Henriette and Isken, Marius Paul},
  title     = {Insight into the 2017-2019 Lurestan arc seismic sequence (Zagros, Iran); complex earthquake interaction in the basement and sediments},
  series = {Geophysical journal international},
  volume    = {230},
  journal   = {Geophysical journal international},
  number    = {1},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggac057},
  pages     = {114 -- 130},
  year      = {2022},
  abstract  = {Despite its high-seismogenic potential, the details of the seismogenic processes of Zagros Simply Folded Belt (SFB) remains debated. Three large earthquakes (M-w 7.3, 5.9 and 6.3) struck in the Lurestan arc of the Zagros SFB in 2017 and 2018. The sequence was recorded by seismic stations at regional, and teleseismic distances. Coseismic surface displacements, measured by Sentinel-1A/B satellites, provide additional data and a unique opportunity to study these earthquakes in detail. Here, we complement previous studies of the coseismic slip distribution of the 12 November 2017 M-w 7.3 Ezgeleh earthquake by a detailed analysis of its aftershocks, and we analysed the rupture process of the two interrelated earthquakes (25 August 2018 M-w 5.9 Tazehabad and the 25 November 2018 M-w 6.3 Sarpol-e Zahab earthquakes). We model the surface displacements obtained from Interferometric Synthetic Aperture Radar (InSAR) measurements and seismic records. We conduct non-linear probabilistic optimizations based on joint InSAR and seismic data to obtain finite-fault rupture of these earthquakes. The Lurestan arc earthquakes were followed by a sustained aftershock activity, with 133 aftershocks exceeding M-n 4.0 until 30 December 2019. We rely on the permanent seismic networks of Iran and Iraq to relocate similar to 700 M-n 3 + events and estimate moment tensor solutions for 85 aftershocks down to M-w 4.0. The 2017 Ezgeleh earthquake has been considered to activate a low-angle (similar to 17 degrees) dextral-thrust fault at the depth of 10-20 km. However, most of its aftershocks have shallow centroid depths (8-12 km). The joint interpretation of finite source models, moment tensor and hypocentral location indicate that the 2018 Tazehabad and Sarpol-e Zahab earthquakes ruptured different strike-slip structures, providing evidence for the activation of the sinistral and dextral strike-slip faults, respectively. The deformation in the Lurestan arc is seismically accommodated by a complex fault system involving both thrust and strike-slip faults. Knowledge about the deformation characteristics is important for the understanding of crustal shortening, faulting and hazard and risk assessment in this region.},
  language  = {en}
}
@article{EiblHainzlVeselyetal.2019,
  author    = {Eibl, Eva P. S. and Hainzl, Sebastian and Vesely, Nele I. K. and Walter, Thomas R. and Jousset, Philippe and Hersir, Gylfi Pall and Dahm, Torsten},
  title     = {Eruption interval monitoring at strokkur Geyser, Iceland},
  series = {Geophysical research letters},
  volume    = {47},
  journal   = {Geophysical research letters},
  number    = {1},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1029/2019GL085266},
  pages     = {10},
  year      = {2019},
  abstract  = {Geysers are hot springs whose frequency of water eruptions remain poorly understood. We set up a local broadband seismic network for 1 year at Strokkur geyser, Iceland, and developed an unprecedented catalog of 73,466 eruptions. We detected 50,135 single eruptions but find that the geyser is also characterized by sets of up to six eruptions in quick succession. The number of single to sextuple eruptions exponentially decreased, while the mean waiting time after an eruption linearly increased (3.7 to 16.4 min). While secondary eruptions within double to sextuple eruptions have a smaller mean seismic amplitude, the amplitude of the first eruption is comparable for all eruption types. We statistically model the eruption frequency assuming discharges proportional to the eruption multiplicity and a constant probability for subsequent events within a multituple eruption. The waiting time after an eruption is predictable but not the type or amplitude of the next one. <br /> Plain Language Summary Geysers are springs that often erupt in hot water fountains. They erupt more often than volcanoes but are quite similar. Nevertheless, it is poorly understood how often volcanoes and also geysers erupt. We created a list of 73,466 eruption times of Strokkur geyser, Iceland, from 1 year of seismic data. The geyser erupted one to six times in quick succession. We found 50,135 single eruptions but only 1 sextuple eruption, while the mean waiting time increased from 3.7 min after single eruptions to 16.4 min after sextuple eruptions. Mean amplitudes of each eruption type were higher for single eruptions, but all first eruptions in a succession were similar in height. Assuming a constant heat inflow at depth, we can predict the waiting time after an eruption but not the type or amplitude of the next one.},
  language  = {en}
}
@article{ValenzuelaMalebranCescaLopezCominoetal.2022,
  author    = {Valenzuela-Malebran, Carla and Cesca, Simone and Lopez-Comino, Jos{\´e} {\´A}ngel and Zeckra, Martin and Kr{\"u}ger, F. and Dahm, Torsten},
  title     = {Source mechanisms and rupture processes of the Jujuy seismic nest, Chile-Argentina border},
  series = {Journal of South American earth sciences},
  volume    = {117},
  journal   = {Journal of South American earth sciences},
  publisher = {Elsevier},
  address   = {Oxford},
  issn      = {0895-9811},
  doi       = {10.1016/j.jsames.2022.103887},
  pages     = {13},
  year      = {2022},
  abstract  = {The Altiplano-Puna plateau, in Central Andes, is the second-largest continental plateau on Earth, extending between 22 degrees and 27 degrees S at an average altitude of 4400 m. The Puna plateau has been formed in consequence of the subduction of the oceanic Nazca Plate beneath the continental South American plate, which has an average crustal thickness of 50 km at this location. A large seismicity cluster, the Jujuy cluster, is observed at depth of 150-250 km beneath the central region of the Puna plateau. The cluster is seismically very active, with hundreds of earthquakes reported and a peak magnitude MW 6.6 on 25th August 2006. The cluster is situated in one of three band of intermediate-depth focus seismicity, which extend parallel to the trench roughly North to South. It has been hypothesized that the Jujuy cluster could be a seismic nest, a compact seismogenic region characterized by a high stationary activity relative to its surroundings. In this study, we collected more than 40 years of data from different catalogs and proof that the cluster meets the three conditions of a seismic nest. Compared to other known intermediate depth nests at Hindu Kush (Afganisthan) or Bucaramanga (Colombia), the Jujuy nest presents an outstanding seismicity rate, with more than 100 M4+ earthquakes per year. We additionally performed a detailed analysis of the rupture process of some of the largest earthquakes in the nest, by means of moment tensor inversion and directivity analysis. We focused on the time period 2017-2018, where the seismic monitoring was the most extended. Our results show that earthquakes in the nest take place within the eastward subducting oceanic plate, but rupture along sub-horizontal planes dipping westward. We suggest that seismicity at Jujuy nest is controlled by dehydration processes, which are also responsible for the generation of fluids ascending to the crust beneath the Puna volcanic region. We use the rupture plane and nest geometry to provide a constraint to maximal expected magnitude, which we estimate as MW -6.7.},
  language  = {en}
}
@misc{DahmBeckerBischoffetal.2013,
  author    = {Dahm, Torsten and Becker, Dirk and Bischoff, Monika and Cesca, Simone and Dost, B. and Fritschen, R. and Hainzl, Sebastian and Klose, C. D. and Kuhn, D. and Lasocki, S. and Meier, Thomas and Ohrnberger, Matthias and Rivalta, Eleonora and Wegler, Ulrich and Husen, Stephan},
  title     = {Recommendation for the discrimination of human-related and natural seismicity},
  series = {Journal of seismology},
  volume    = {17},
  journal   = {Journal of seismology},
  number    = {1},
  publisher = {Springer},
  address   = {Dordrecht},
  issn      = {1383-4649},
  doi       = {10.1007/s10950-012-9295-6},
  pages     = {197 -- 202},
  year      = {2013},
  abstract  = {Various techniques are utilized by the seismological community, extractive industries, energy and geoengineering companies to identify earthquake nucleation processes in close proximity to engineering operation points. These operations may comprise fluid extraction or injections, artificial water reservoir impoundments, open pit and deep mining, deep geothermal power generations or carbon sequestration. In this letter to the editor, we outline several lines of investigation that we suggest to follow to address the discrimination problem between natural seismicity and seismic events induced or triggered by geoengineering activities. These suggestions have been developed by a group of experts during several meetings and workshops, and we feel that their publication as a summary report is helpful for the geoscientific community. Specific investigation procedures and discrimination approaches, on which our recommendations are based, are also published in this Special Issue (SI) of Journal of Seismology.},
  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{SenCescaBischoffetal.2013,
  author    = {Sen, Ali Tolga and Cesca, Simone and Bischoff, Monika and Meier, Thomas and Dahm, Torsten},
  title     = {Automated full moment tensor inversion of coal mining-induced seismicity},
  series = {Geophysical journal international},
  volume    = {195},
  journal   = {Geophysical journal international},
  number    = {2},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggt300},
  pages     = {1267 -- 1281},
  year      = {2013},
  abstract  = {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.},
  language  = {en}
}
@article{HeimannGonzalezWangetal.2013,
  author    = {Heimann, Sebastian and Gonzalez, Alvaro and Wang, Rongjiang and Cesca, Simone and Dahm, Torsten},
  title     = {Seismic characterization of the Chelyabinsk Meteor's terminal explosion},
  series = {Seismological research letters},
  volume    = {84},
  journal   = {Seismological research letters},
  number    = {6},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0895-0695},
  doi       = {10.1785/0220130042},
  pages     = {1021 -- 1025},
  year      = {2013},
  language  = {en}
}
@article{GrigoliCescaVassalloetal.2013,
  author    = {Grigoli, Francesco and Cesca, Simone and Vassallo, Maurizio and Dahm, Torsten},
  title     = {Automated seismic event location by travel-time stacking an application to mining induced seismicity},
  series = {Seismological research letters},
  volume    = {84},
  journal   = {Seismological research letters},
  number    = {4},
  publisher = {Seismological Society of America},
  address   = {Albany},
  issn      = {0895-0695},
  doi       = {10.1785/0220120191},
  pages     = {666 -- 677},
  year      = {2013},
  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{HannemannEulenfeldKruegeretal.2021,
  author    = {Hannemann, Katrin and Eulenfeld, Tom and Kr{\"u}ger, Frank and Dahm, Torsten},
  title     = {Seismic scattering and absorption of oceanic lithospheric S waves in the Eastern North Atlantic},
  series = {Geophysical journal international},
  volume    = {229},
  journal   = {Geophysical journal international},
  number    = {2},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggab493},
  pages     = {948 -- 961},
  year      = {2021},
  abstract  = {The scattering and absorption of high-frequency seismic waves in the oceanic lithosphere is to date only poorly constrained by observations. Such estimates would not only improve our understanding of the propagation of seismic waves, but also unravel the small-scale nature of the lithosphere and its variability. Our study benefits from two exceptional situations: (1) we deployed over 10 months a mid-aperture seismological array in the central part of the Eastern North Atlantic in 5 km water depth and (2) we could observe in total 340 high-frequency (up to 30 Hz) Po and So arrivals with tens to hundreds of seconds long seismic coda from local and regional earthquakes in a wide range of backazimuths and epicentral distances up to 850 km with a travel path in the oceanic lithosphere. Moreover, the array was located about 100 km north of the Gloria fault, defining the plate boundary between the Eurasian and African plates at this location which also allows an investigation of the influence of an abrupt change in lithospheric age (20 Ma in this case) on seismic waves. The waves travel with velocities indicating upper-mantle material. We use So waves and their coda of pre-selected earthquakes to estimate frequency-dependent seismic scattering and intrinsic attenuation parameters. The estimated scattering attenuation coefficients are between 10(-4) and 4 x 10(-5) m(-1) and are typical for the lithosphere or the upper mantle. Furthermore, the total quality factors for So waves below 5 Hz are between 20 and 500 and are well below estimates from previous modelling for observations in the Pacific Ocean. This implies that the Atlantic Ocean is more attenuative for So waves compared to the Pacific Ocean, which is inline with the expected behaviour for the lithospheric structures resulting from the slower spreading rates in the Atlantic Ocean. The results for the analysed events indicate that for frequencies above 3 Hz, intrinsic attenuation is equal to or slightly stronger than scattering attenuation and that the So-wave coda is weakly influenced by the oceanic crust. Both observations are in agreement with the proposed propagation mechanism of scattering in the oceanic mantle lithosphere. Furthermore, we observe an age dependence which shows that an increase in lithospheric age is associated with a decrease in attenuation. However, we also observe a trade-off of this age-dependent effect with either a change in lithospheric thickness or thermal variations, for example due to small-scale upwellings in the upper mantle in the southeast close to Madeira and the Canaries. Moreover, the influence of the nearby Gloria fault is visible in a reduction of the intrinsic attenuation below 3 Hz for estimates across the fault. This is the first study to estimate seismic scattering and absorption parameters of So waves for an area with several hundreds of kilometres radius centred in the Eastern North Atlantic and using them to characterize the nature of the oceanic lithosphere.},
  language  = {en}
}