TY - GEN A1 - Al-Halbouni, Djamil A1 - Holohan, Eoghan P. A1 - Taheri, Abbas A1 - Schöpfer, Martin P. J. A1 - Emam, Sacha A1 - Dahm, Torsten T1 - Geomechanical modelling of sinkhole development using distinct elements BT - model verification for a single void space and application to the Dead Sea area T2 - Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe N2 - Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 1061 KW - rock mass KW - karst KW - dissolution KW - reflection KW - subsidence KW - subrosion KW - collapse KW - simulation KW - scale KW - fault Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-468435 SN - 1866-8372 IS - 1061 ER - TY - JOUR A1 - Nooshiri, Nima A1 - Bean, Christopher J. A1 - Dahm, Torsten A1 - Grigoli, Francesco A1 - Kristjansdottir, Sigriour A1 - Obermann, Anne A1 - Wiemer, Stefan T1 - A multibranch, multitarget neural network for rapid point-source inversion in a microseismic environment BT - examples from the Hengill Geothermal Field, Iceland JF - Geophysical journal international N2 - Despite advanced seismological techniques, automatic source characterization for microseismic earthquakes remains difficult and challenging since current inversion and modelling of high-frequency signals are complex and time consuming. For real-time applications such as induced seismicity monitoring, the application of standard methods is often not fast enough for true complete real-time information on seismic sources. In this paper, we present an alternative approach based on recent advances in deep learning for rapid source-parameter estimation of microseismic earthquakes. The seismic inversion is represented in compact form by two convolutional neural networks, with individual feature extraction, and a fully connected neural network, for feature aggregation, to simultaneously obtain full moment tensor and spatial location of microseismic sources. Specifically, a multibranch neural network algorithm is trained to encapsulate the information about the relationship between seismic waveforms and underlying point-source mechanisms and locations. The learning-based model allows rapid inversion (within a fraction of second) once input data are available. A key advantage of the algorithm is that it can be trained using synthetic seismic data only, so it is directly applicable to scenarios where there are insufficient real data for training. Moreover, we find that the method is robust with respect to perturbations such as observational noise and data incompleteness (missing stations). We apply the new approach on synthesized and example recorded small magnitude (M <= 1.6) earthquakes at the Hellisheioi geothermal field in the Hengill area, Iceland. For the examined events, the model achieves excellent performance and shows very good agreement with the inverted solutions determined through standard methodology. In this study, we seek to demonstrate that this approach is viable for microseismicity real-time estimation of source parameters and can be integrated into advanced decision-support tools for controlling induced seismicity. KW - Neural networks KW - fuzzy logic KW - Computational seismology KW - Induced seismicity KW - Earthquake source observations Y1 - 2021 U6 - https://doi.org/10.1093/gji/ggab511 SN - 0956-540X SN - 1365-246X VL - 229 IS - 2 SP - 999 EP - 1016 PB - Oxford Univ. Press CY - Oxford ER - TY - GEN A1 - Polom, Ulrich A1 - Alrshdan, Hussam A1 - Al-Halbouni, Djamil A1 - Holohan, Eoghan P. A1 - Dahm, Torsten A1 - Sawarieh, Ali A1 - Atallah, Mohamad Y. A1 - Krawczyk, Charlotte M. T1 - Shear wave reflection seismic yields subsurface dissolution and subrosion patterns BT - application to the Ghor Al-Haditha sinkhole site, Dead Sea, Jordan T2 - Postprints der Universität Potsdam : Mathematisch Naturwissenschaftliche Reihe N2 - Near-surface geophysical imaging of alluvial fan settings is a challenging task but crucial for understating geological processes in such settings. The alluvial fan of Ghor Al-Haditha at the southeast shore of the Dead Sea is strongly affected by localized subsidence and destructive sinkhole collapses, with a significantly increasing sinkhole formation rate since ca. 1983. A similar increase is observed also on the western shore of the Dead Sea, in correlation with an ongoing decline in the Dead Sea level. Since different structural models of the upper 50 m of the alluvial fan and varying hypothetical sinkhole processes have been suggested for the Ghor Al-Haditha area in the past, this study aimed to clarify the subsurface characteristics responsible for sinkhole development. For this purpose, high-frequency shear wave reflection vibratory seismic surveys were carried out in the Ghor Al-Haditha area along several crossing and parallel profiles with a total length of 1.8 and 2.1 km in 2013 and 2014, respectively. The sedimentary architecture of the alluvial fan at Ghor Al-Haditha is resolved down to a depth of nearly 200 m at a high resolution and is calibrated with the stratigraphic profiles of two boreholes located inside the survey area. The most surprising result of the survey is the absence of evidence of a thick (> 2–10 m) compacted salt layer formerly suggested to lie at ca. 35–40 m depth. Instead, seismic reflection amplitudes and velocities image with good continuity a complex interlocking of alluvial fan deposits and lacustrine sediments of the Dead Sea between 0 and 200 m depth. Furthermore, the underground section of areas affected by sinkholes is characterized by highly scattering wave fields and reduced seismic interval velocities. We propose that the Dead Sea mud layers, which comprise distributed inclusions or lenses of evaporitic chloride, sulfate, and carbonate minerals as well as clay silicates, become increasingly exposed to unsaturated water as the sea level declines and are consequently destabilized and mobilized by both dissolution and physical erosion in the subsurface. This new interpretation of the underlying cause of sinkhole development is supported by surface observations in nearby channel systems. Overall, this study shows that shear wave seismic reflection technique is a promising method for enhanced near-surface imaging in such challenging alluvial fan settings. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 979 KW - salt dissolution KW - hazard KW - coast KW - area KW - subsidence KW - shoreline KW - karst KW - lake Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-459134 SN - 1866-8372 IS - 979 SP - 1079 EP - 1098 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 - Niemz, Peter A1 - Dahm, Torsten A1 - Milkereit, Claus A1 - Cesca, Simone A1 - Petersen, Gesa Maria A1 - Zang, Arno T1 - Insights into hydraulic fracture growth gained from a joint analysis of seismometer-derived tilt signals and scoustic emissions JF - Journal of geophysical research : Solid earth N2 - 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. KW - hydraulic fracturing KW - fracture growth KW - tilt KW - acoustic emissions KW - injections KW - broadband seismometer Y1 - 2021 U6 - https://doi.org/10.1029/2021JB023057 SN - 2169-9313 SN - 2169-9356 VL - 126 IS - 12 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Krüger, Frank A1 - Dahm, Torsten A1 - Hannemann, Katrin T1 - Mapping of Eastern North Atlantic Ocean seismicity from Po/So observations at a mid-aperture seismological broad-band deep sea array JF - Geophysical journal international N2 - 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. KW - body waves KW - earthquake source observations KW - seismicity and tectonics Y1 - 2020 U6 - https://doi.org/10.1093/gji/ggaa054 SN - 0956-540X SN - 1365-246X VL - 221 IS - 2 SP - 1055 EP - 1080 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Krüger, Frank A1 - Dahm, Torsten T1 - Higher degree moment inversin using far-field broad-band recordings : theory and evaluation of the method with application to the 1994 Bolivia deep earthauke Y1 - 1999 ER - TY - JOUR A1 - Eibl, Eva P. S. A1 - Hainzl, Sebastian A1 - Vesely, Nele I. K. A1 - Walter, Thomas R. A1 - Jousset, Philippe A1 - Hersir, Gylfi Pall A1 - Dahm, Torsten T1 - Eruption interval monitoring at strokkur Geyser, Iceland JF - Geophysical research letters N2 - 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.
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. Y1 - 2019 U6 - https://doi.org/10.1029/2019GL085266 SN - 0094-8276 SN - 1944-8007 VL - 47 IS - 1 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Nooshiri, Nima A1 - Saul, Joachim A1 - Heimann, Sebastian A1 - Tilmann, Frederik A1 - Dahm, Torsten T1 - Revision of earthquake hypocentre locations in global bulletin data sets using source-specific station terms JF - Geophysical journal international N2 - 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. KW - Seismicity and tectonics KW - Computational seismology KW - Subduction zone processes KW - Pacific Ocean KW - South America Y1 - 2016 U6 - https://doi.org/10.1093/gji/ggw405 SN - 0956-540X SN - 1365-246X VL - 208 IS - 2 SP - 589 EP - 602 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Woith, Heiko A1 - Petersen, Gesa Maria A1 - Hainzl, Sebastian A1 - Dahm, Torsten T1 - Review: Can Animals Predict Earthquakes? JF - Bulletin of the Seismological Society of America N2 - In public perception, abnormal animal behavior is widely assumed to be a potential earthquake precursor, in strong contrast to the viewpoint in natural sciences. Proponents of earthquake prediction via animals claim that animals feel and react abnormally to small changes in environmental and physico-chemical parameters related to the earthquake preparation process. In seismology, however, observational evidence for changes of physical parameters before earthquakes is very weak. In this study, we reviewed 180 publications regarding abnormal animal behavior before earthquakes and analyze and discuss them with respect to (1) magnitude-distance relations, (2) foreshock activity, and (3) the quality and length of the published observations. More than 700 records of claimed animal precursors related to 160 earthquakes are reviewed with unusual behavior of more than 130 species. The precursor time ranges from months to seconds prior to the earthquakes, and the distances from a few to hundreds of kilometers. However, only 14 time series were published, whereas all other records are single observations. The time series are often short (the longest is 1 yr), or only small excerpts of the full data set are shown. The probability density of foreshocks and the occurrence of animal precursors are strikingly similar, suggesting that at least parts of the reported animal precursors are in fact related to foreshocks. Another major difficulty for a systematic and statistical analysis is the high diversity of data, which are often only anecdotal and retrospective. The study clearly demonstrates strong weaknesses or even deficits in many of the published reports on possible abnormal animal behavior. To improve the research on precursors, we suggest a scheme of yes and no questions to be assessed to ensure the quality of such claims. Y1 - 2018 U6 - https://doi.org/10.1785/0120170313 SN - 0037-1106 SN - 1943-3573 VL - 108 IS - 3A SP - 1031 EP - 1045 PB - Seismological Society of America CY - Albany ER -