TY - JOUR 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 JF - Solid earth 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. Y1 - 2018 U6 - https://doi.org/10.5194/se-9-1079-2018 SN - 1869-9510 SN - 1869-9529 VL - 9 IS - 5 SP - 1079 EP - 1098 PB - Copernicus CY - Göttingen ER - TY - JOUR 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 JF - Solid earth 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. Y1 - 2018 U6 - https://doi.org/10.5194/se-9-1341-2018 SN - 1869-9510 SN - 1869-9529 VL - 9 IS - 6 SP - 1341 EP - 1373 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Hensch, Martin A1 - Dahm, Torsten A1 - Ritter, Joachim A1 - Heimann, Sebastian A1 - Schmidt, Bernd A1 - Stange, Stefan A1 - Lehmann, Klaus T1 - Deep low-frequency earthquakes reveal ongoing magmatic recharge beneath Laacher See Volcano (Eifel, Germany) JF - Geophysical journal international N2 - The occurrence of deep low-frequency (DLF) microearthquakes beneath volcanoes is commonly attributed to mass transport in the volcanic plumbing system and used to infer feeding channels from and into magma reservoirs. The key question is how magmas migrate from depth to the shallow crust and whether magma reservoirs are currently being recharged. For the first time since the improvement of the local seismic networks in the East Eifel region (Rhineland-Palatinate, Germany), we detect and locate recurrent DLF earthquakes in the lower crust and upper mantle beneath the Laacher See Volcano (LSV), using a joint data set of permanent sensors and a temporary deployment. So far, eight DLF earthquake sequences were observed in four distinct clusters between 10 and 40 km depth. These clusters of weak events (M-L< 2) align along an approximately 80. southeast dipping line south of the LSV. Moment tensor solutions of these events have large shear components, and the irregular dispersion and long coda of body waves indicate interaction processes between shear cracks and fluids. We find a rotation of P-axes orientation for shallow tectonic earthquakes compared to DLF events, indicating that the stress field in the depth interval of DLF events might favour a vertical migration of magma or magmatic fluids. The caldera of the LSV was formed by the last major eruption of the East Eifel Volcanic Field only 12.9 kyr ago, fed by a shallow magma chamber at 5-8 km depth and erupting a total magma volume of 6.7 km(3). The observed DLF earthquake activity and continuous volcanic gas emissions around the LSV indicate an active magmatic system, possibly connected with an upper mantle melt zone. KW - Waveform inversion KW - Volcano seismology KW - Magma migration and fragmentation KW - Volcano monitoring Y1 - 2019 U6 - https://doi.org/10.1093/gji/ggy532 SN - 0956-540X SN - 1365-246X VL - 216 IS - 3 SP - 2025 EP - 2036 PB - Oxford Univ. Press CY - Oxford ER - TY - JOUR A1 - Cesca, Simone A1 - Heimann, Sebastian A1 - Kriegerowski, Marius A1 - Saul, Joachim A1 - Dahm, Torsten T1 - Moment tensor inversion for nuclear explosions BT - what can we learn from the 6 January and 9 September 2016 Nuclear Tests, North Korea? JF - Seismological research letters N2 - Two nuclear explosions were carried out by the Democratic People’s Republic of North Korea in January and September 2016. Epicenters were located close to those of the 2006, 2009, and 2013 previous explosions. We perform a seismological analysis of the 2016 events combining the analysis of full waveforms at regional distances and seismic array beams at teleseismic distances. We estimate the most relevant source parameters, such as source depth, moment release, and full moment tensor (MT). The best MT solution can be decomposed into an isotropic source, directly related with the explosion and an additional deviatoric term, likely due to near‐source interactions with topographic and/or underground facilities features. We additionally perform an accurate resolution test to assess source parameters uncertainties and trade‐offs. This analysis sheds light on source parameters inconsistencies among studies on previous shallow explosive sources. The resolution of the true MT is hindered by strong source parameters trade‐offs, so that a broad range of well‐fitting MT solutions can be found, spanning from a dominant positive isotropic term to a dominant negative vertical compensated linear vector dipole. The true mechanism can be discriminated by additionally modeling first‐motion polarities at seismic arrays at teleseismic distances. A comparative assessment of the 2016 explosion with earlier nuclear tests documents similar vertical waveforms but a significant increase of amplitude for the 2016 explosions, which proves that the 9 September 2016 was the largest nuclear explosion ever performed in North Korea with a magnitude Mw 4.9 and a shallow depth of less than 2 km, although there are no proofs of a fusion explosion. Modeling transversal component waveforms suggests variable size and orientation of the double‐couple components of the 2009, 2013, and 2016 sources. Y1 - 2017 U6 - https://doi.org/10.1785/0220160139 SN - 0895-0695 SN - 1938-2057 VL - 88 IS - 2A SP - 300 EP - 310 PB - Seismological Society of America CY - Albany 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 - GEN A1 - Dahm, Torsten A1 - Heimann, Sebastian A1 - Funke, Sigward A1 - Wendt, Siegfried A1 - Rappsilber, Ivo A1 - Bindi, Dino A1 - Plenefisch, Thomas A1 - Cotton, Fabrice Pierre T1 - Correction to: Seismicity in the block mountains between Halle and Leipzig, Central Germany: centroid moment tensors, ground motion simulation, and felt intensities of two M approximate to 3 earthquakes in 2015 and 2017 (vol 22, pg 985, 2018) T2 - Journal of seismology Y1 - 2018 U6 - https://doi.org/10.1007/s10950-018-9773-6 SN - 1383-4649 SN - 1573-157X VL - 22 IS - 6 SP - 1669 EP - 1671 PB - Springer CY - Dordrecht ER - TY - JOUR A1 - Al-Halbouni, Djamil A1 - Holohan, Eoghan P. A1 - Taheri, Abbas A1 - Watson, Robert A. A1 - Polom, Ulrich A1 - Schoepfer, Martin P. J. A1 - Emam, Sacha A1 - Dahm, Torsten T1 - Distinct element geomechanical modelling of the formation of sinkhole clusters within large-scale karstic depressions JF - Solid earth N2 - The 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvalalike depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea. Y1 - 2019 U6 - https://doi.org/10.5194/se-10-1219-2019 SN - 1869-9510 SN - 1869-9529 VL - 10 IS - 4 SP - 1219 EP - 1241 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Karamzadeh, Nasim Toularoud A1 - Kühn, Daniela A1 - Kriegerowski, Marius A1 - López-Comino, José Ángel A1 - Cesca, Simone A1 - Dahm, Torsten T1 - Small-aperture array as a tool to monitor fluid injection- and extraction-induced microseismicity BT - applications and recommendations JF - Acta Geophysica N2 - The monitoring of microseismicity during temporary human activities such as fluid injections for hydrofracturing, hydrothermal stimulations or wastewater disposal is a difficult task. The seismic stations often cannot be installed on hard rock, and at quiet places, noise is strongly increased during the operation itself and the installation of sensors in deep wells is costly and often not feasible. The combination of small-aperture seismic arrays with shallow borehole sensors offers a solution. We tested this monitoring approach at two different sites: (1) accompanying a fracking experiment in sedimentary shale at 4km depth and (2) above a gas field under depletion. The small-aperture arrays were planned according to theoretical wavenumber studies combined with simulations considering the local noise conditions. We compared array recordings with recordings available from shallow borehole sensors and give examples of detection and location performance. Although the high-frequency noise on the 50-m-deep borehole sensors was smaller compared to the surface noise before the injection experiment, the signals were highly contaminated during injection by the pumping activities. Therefore, a set of three small-aperture arrays at different azimuths was more suited to detect small events, since noise recorded on these arrays is uncorrelated with each other. Further, we developed recommendations for the adaptation of the monitoring concept to other sites experiencing induced seismicity. KW - Microseismic monitoring KW - Induced seismicity KW - Array seismology KW - Shallow borehole sensors Y1 - 2019 U6 - https://doi.org/10.1007/s11600-018-0231-1 SN - 1895-6572 SN - 1895-7455 VL - 67 IS - 1 SP - 311 EP - 326 PB - Springer CY - Cham 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 - Maghsoudi, Samira A1 - Cesca, Simone A1 - Hainzl, Sebastian A1 - Kaiser, Diethelm A1 - Becker, Dirk A1 - Dahm, Torsten T1 - Improving the estimation of detection probability and magnitude of completeness in strongly heterogeneous media, an application to acoustic emission (AE) JF - Geophysical journal international N2 - 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. KW - Seismic attenuation KW - Statistical seismology Y1 - 2013 U6 - https://doi.org/10.1093/gji/ggt049 SN - 0956-540X VL - 193 IS - 3 SP - 1556 EP - 1569 PB - Oxford Univ. Press CY - Oxford ER -