2017-12-19, Wadas, Sonja H., Tanner, David C., Polom, Ulrich, Krawczyk, Charlotte M.

In November 2010, a large sinkhole opened up in the urban area of Schmalkalden, Germany. To determine the key factors which benefited the development of this collapse structure and therefore the dissolution, we carried out several shear-wave reflection-seismic profiles around the sinkhole. In the seismic sections we see evidence of the Mesozoic tectonic movement in the form of a NW–SE striking, dextral strike-slip fault, known as the Heßleser Fault, which faulted and fractured the subsurface below the town. The strike-slip faulting created a zone of small blocks ( < 100 m in size), around which steep-dipping normal faults, reverse faults and a dense fracture network serve as fluid pathways for the artesian-confined groundwater. The faults also acted as barriers for horizontal groundwater flow perpendicular to the fault planes. Instead groundwater flows along the faults which serve as conduits and forms cavities in the Permian deposits below ca. 60 m depth. Mass movements and the resulting cavities lead to the formation of sinkholes and dissolution-induced depressions. Since the processes are still ongoing, the occurrence of a new sinkhole cannot be ruled out. This case study demonstrates how S-wave seismics can characterize a sinkhole and, together with geological information, can be used to study the processes that result in sinkhole formation, such as a near-surface fault zone located in soluble rocks. The more complex the fault geometry and interaction between faults, the more prone an area is to sinkhole occurrence.

2016, Wadas, Sonja H., Polom, Ulrich, Krawczyk, Charlotte M.

Subrosion is the subsurface leaching of soluble rocks that results in the formation of depression and collapse structures. This global phenomenon is a geohazard in urban areas. To study near-surface subrosion structures, four shear-wave seismic reflection profiles, with a total length of ca. 332 m, were carried out around the famous leaning church tower of Bad Frankenhausen in northern Thuringia, Germany, which shows an inclination of 4.93° from the vertical. Most of the geological underground of Thuringia is characterized by soluble Permian deposits, and the Kyffhäuser Southern Margin Fault is assumed to be a main pathway for water to leach the evaporite. The seismic profiles were acquired with the horizontal micro-vibrator ELVIS, developed at Leibniz Institute for Applied Geophysics (LIAG), and a 72 m long landstreamer equipped with 72 horizontal geophones. The high-resolution seismic sections show subrosion-induced structures to a depth of ca. 100 m and reveal five features associated with the leaching of Permian deposits: (1) lateral and vertical varying reflection patterns caused by strongly heterogeneous strata, (2) discontinuous reflectors, small offsets, and faults, which show the underground is heavily fractured, (3) formation of depression structures in the near-surface, (4) diffractions in the unmigrated seismic sections that indicate increased scattering of the seismic waves, and (5) varying seismic velocities and low-velocity zones that are presumably caused by fractures and upward-migrating cavities. A previously undiscovered southward-dipping listric normal fault was also found, to the north of the church. It probably serves as a pathway for water to leach the Permian formations below the church and causes the tilting of the church tower. This case study shows the potential of horizontal shear-wave seismic reflection to image near-surface subrosion structures in an urban environment with a horizontal resolution of less than 1 m in the uppermost 10–15 m.

2012, Burschil, T., Scheer, W., Kirsch, R., Wiederhold, H.

Within the scope of climatic change and associated sea level rise, coastal aquifers are endangered and are becoming more a focus of research to ensure the future water supply in coastal areas. For groundwater modelling a good understanding of the geological/hydrogeological situation and the aquifer behavior is necessary. In preparation of groundwater modelling and assessment of climate change impacts on coastal water resources, we setup a geological/hydrogeological model for the North Sea Island of Föhr. Data from different geophysical methods applied from the air, the surface and in boreholes contribute to the 3-D model, e.g. airborne electromagnetics (SkyTEM) for spatial mapping the resistivity of the entire island, seismic reflections for detailed cross-sections in the groundwater catchment area, and geophysical borehole logging for calibration of these measurements. An iterative and integrated evaluation of the results from the different geophysical methods contributes to reliable data as input for the 3-D model covering the whole island and not just the well fields. The complex subsurface structure of the island is revealed. The local waterworks use a freshwater body embedded in saline groundwater. Several glaciations reordered the youngest Tertiary and Quaternary sediments by glaciotectonic thrust faulting, as well as incision and refill of glacial valleys. Both subsurface structures have a strong impact on the distribution of freshwater-bearing aquifers. A digital geological 3-D model reproduces the hydrogeological structure of the island as a base for a groundwater model. In the course of the data interpretation, we deliver a basis for rock identification. We demonstrate that geophysical investigation provide petrophysical parameters and improve the understanding of the subsurface and the groundwater system. The main benefit of our work is that the successful combination of electromagnetic, seismic and borehole data reveals the complex geology of a glacially-affected island. A sound understanding of the subsurface structure and the compilation of a 3-D model is imperative and the basis for a groundwater flow model to predict climate change effects on future water resources.

2013, Krawczyk, C.M., Buddensiek, M.-L., Oncken, O., Kukowski, N.

With the study and technical development introduced here, we combine analogue sandbox simulation techniques with seismic physical modelling of sandbox models. For that purpose, we designed and developed a new mini-seismic facility for laboratory use, comprising a seismic tank, a PC-driven control unit, a positioning system, and piezoelectric transducers used here for the first time in an array mode. To assess the possibilities and limits of seismic imaging of small-scale structures in sandbox models, different geometry setups were tested in the first 2-D experiments that also tested the proper functioning of the device and studied the seismo-elastic properties of the granular media used. Simple two-layer models of different materials and layer thicknesses as well as a more complex model comprising channels and shear zones were tested using different acquisition geometries and signal properties. We suggest using well sorted and well rounded grains with little surface roughness (glass beads). Source receiver-offsets less than 14 cm for imaging structures as small as 2.0–1.5 mm size have proven feasible. This is the best compromise between wide beam and high energy output, and is applicable with a consistent waveform. Resolution of the interfaces of layers of granular materials depends on the interface preparation rather than on the material itself. Flat grading of interfaces and powder coverage yields the clearest interface reflections. Finally, sandbox seismic sections provide images of high quality showing constant thickness layers as well as predefined channel structures and indications of the fault traces from shear zones. Since these were artificially introduced in our test models, they can be regarded as zones of disturbance rather than tectonic shear zones characterized by decompaction. The multiple-offset surveying introduced here, improves the quality with respect to S / N ratio and source signature even more; the maximum depth penetration in glass-bead layers thereby amounts to 5 cm. Thus, the presented mini-seismic device is already able to resolve structures within simple models of saturated porous media, so that multiple-offset seismic imaging of shallow sandbox models, that are structurally evolving, is generally feasible.

2015, Burschil, T., Beilecke, T., Krawczyk, C.M.

High-resolution reflection seismic methods are an established non-destructive tool for engineering tasks. In the near surface, shear-wave reflection seismic measurements usually offer a higher spatial resolution in the same effective signal frequency spectrum than P-wave data, but data quality varies more strongly. To discuss the causes of these differences, we investigated a P-wave and a SH-wave seismic reflection profile measured at the same location on the island of Föhr, Germany and applied seismic reflection processing to the field data as well as finite-difference modelling of the seismic wave field. The simulations calculated were adapted to the acquisition field geometry, comprising 2 m receiver distance (1 m for SH wave) and 4 m shot distance along the 1.5 km long P-wave and 800 m long SH-wave profiles. A Ricker wavelet and the use of absorbing frames were first-order model parameters. The petrophysical parameters to populate the structural models down to 400 m depth were taken from borehole data, VSP (vertical seismic profile) measurements and cross-plot relations. The simulation of the P-wave wave-field was based on interpretation of the P-wave depth section that included a priori information from boreholes and airborne electromagnetics. Velocities for 14 layers in the model were derived from the analysis of five nearby VSPs (vP =1600–2300 m s-1). Synthetic shot data were compared with the field data and seismic sections were created. Major features like direct wave and reflections are imaged. We reproduce the mayor reflectors in the depth section of the field data, e.g. a prominent till layer and several deep reflectors. The SH-wave model was adapted accordingly but only led to minor correlation with the field data and produced a higher signal-to-noise ratio. Therefore, we suggest to consider for future simulations additional features like intrinsic damping, thin layering, or a near-surface weathering layer. These may lead to a better understanding of key parameters determining the data quality of near-surface shear-wave seismic measurements.