Search Results

Now showing 1 - 10 of 11
  • Item
    Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise
    (Munich : EGU, 2012) Jørgensen, F.; Scheer, W.; Thomsen, S.; Sonnenborg, T.O.; Hinsby, K.; Wiederhold, H.; Schamper, C.; Burschil, T.; Roth, B.; Kirsch, R.; Auken, E.
    Geophysical techniques are increasingly being used as tools for characterising the subsurface, and they are generally required to develop subsurface models that properly delineate the distribution of aquifers and aquitards, salt/freshwater interfaces, and geological structures that affect groundwater flow. In a study area covering 730 km2 across the border between Germany and Denmark, a combination of an airborne electromagnetic survey (performed with the SkyTEM system), a high-resolution seismic survey and borehole logging has been used in an integrated mapping of important geological, physical and chemical features of the subsurface. The spacing between flight lines is 200–250 m which gives a total of about 3200 line km. About 38 km of seismic lines have been collected. Faults bordering a graben structure, buried tunnel valleys, glaciotectonic thrust complexes, marine clay units, and sand aquifers are all examples of geological structures mapped by the geophysical data that control groundwater flow and to some extent hydrochemistry. Additionally, the data provide an excellent picture of the salinity distribution in the area and thus provide important information on the salt/freshwater boundary and the chemical status of groundwater. Although the westernmost part of the study area along the North Sea coast is saturated with saline water and the TEM data therefore are strongly influenced by the increased electrical conductivity there, buried valleys and other geological elements are still revealed. The mapped salinity distribution indicates preferential flow paths through and along specific geological structures within the area. The effects of a future sea level rise on the groundwater system and groundwater chemistry are discussed with special emphasis on the importance of knowing the existence, distribution and geometry of the mapped geological elements, and their control on the groundwater salinity distribution is assessed.
  • Item
    Zero-Offset VSP Monitoring of CO2Storage: Impedance Inversion and Wedge Modelling at the Ketzin Pilot Site
    (New York, NY : Hindawi, 2014) Götz, Julia; Lüth, Stefan; Krawczyk, Charlotte M.; Cosma, Calin
    At the CO2 storage pilot site near the town of Ketzin (35 km west of Berlin, Germany) the sandstone reservoir at 630 m–650 m depth is thin and heterogeneous. The time-lapse analysis of zero-offset VSP measurements shows that CO2-induced amplitude changes can be observed on near-well corridor stacks. Further, we investigate whether CO2-induced amplitude changes in the monitoring data can be used to derive geometrical and petrophysical parameters governing the migration of CO2 within a brine saturated sandstone aquifer. 2D seismic-elastic modelling is done to test the processing workflow and to perform a wedge modelling study for estimation of the vertical expansion of the CO2 plume. When using the NRMS error as a measure for the similarity between the modelled and recorded repeat traces, the best match is achieved for a plume thickness of 6-7 m within the reservoir sandstone of 8 m thickness. With band limited impedance inversion a velocity reduction at the top of the reservoir of 30%, influenced by casing reverberations as well as CO2 injection, is found. The relation of seismic amplitude to CO2 saturated layer thickness and CO2-induced changes in P-wave velocities are important parameters for the quantification of the injected CO2 volume.
  • Item
    Numerical modelling of climate change impacts on freshwater lenses on the North Sea Island of Borkum using hydrological and geophysical methods
    (Munich : EGU, 2012) Sulzbacher, H.; Wiederhold, H.; Siemon, B.; Grinat, M.; Igel, J.; Burschil, T.; Günther, T.; Hinsby, K.
    A numerical, density dependent groundwater model is set up for the North Sea Island of Borkum to estimate climate change impacts on coastal aquifers and especially the situation of barrier islands in the Wadden Sea. The database includes information from boreholes, a seismic survey, a helicopter-borne electromagnetic (HEM) survey, monitoring of the freshwater-saltwater boundary by vertical electrode chains in two boreholes, measurements of groundwater table, pumping and slug tests, as well as water samples. Based on a statistical analysis of borehole columns, seismic sections and HEM, a hydrogeological model is set up. The groundwater model is developed using the finite-element programme FEFLOW. The density dependent groundwater model is calibrated on the basis of hydraulic, hydrological and geophysical data, in particular spatial HEM and local monitoring data. Verification runs with the calibrated model show good agreement between measured and computed hydraulic heads. A good agreement is also obtained between measured and computed density or total dissolved solids data for both the entire freshwater lens on a large scale and in the area of the well fields on a small scale. For simulating future changes in this coastal groundwater system until the end of the current century, we use the climate scenario A2, specified by the Intergovernmental Panel on Climate Change and, in particular, the data for the German North Sea coast. Simulation runs show proceeding salinisation with time beneath the well fields of the two waterworks Waterdelle and Ostland. The modelling study shows that the spreading of well fields is an appropriate protection measure against excessive salinisation of the water supply until the end of the current century.
  • Item
    Spectral induced polarization measurements for predicting the hydraulic conductivity in sandy aquifers
    (Munich : EGU, 2013) Attwa, M.; Günther, T.
    Field and laboratory spectral induced polarization (SIP) measurements are integrated to characterize the hydrogeological conditions at the Schillerslage test site in Germany. The phase images are capable of monitoring thin peat layers within the sandy aquifers. However, the field results show limitations of decreasing resolution with depth. In comparison with the field inversion results, the SIP laboratory measurements show a certain shift in SIP response due to different compaction and sorting of the samples. The SIP data are analyzed to derive an empirical relationship for predicting the hydraulic conductivity (K). In particular, two significant but weak correlations between individual real resistivities (ρ') and relaxation times (τ), based on a Debye decomposition (DD) model, with measured K are found for the upper groundwater aquifer. The maximum relaxation time (τmax) and logarithmically weighted average relaxation time (τlw) show a better relation with K values than the median value τ50. A combined power law relation between individual ρ' and τ with K is developed with an expression of A · (ρ')B · (τlw)C, where A, B and C are determined using a least-squares fit between the measured and predicted K. The suggested approach with the calculated coefficients of the first aquifer is applied for the second. Results show good correlation with the measured K indicating that the derived relationship is superior to single phase angle models as Börner or Slater models.
  • Item
    Ground-penetrating radar insight into a coastal aquifer: the freshwater lens of Borkum Island
    (Munich : EGU, 2013) Igel, J.; Günther, T.; Kuntzer, M.
    Freshwater lenses, as important resource for drinking water, are sensitive to climate changes and sea level rise. To simulate this impact on the groundwater systems, hydraulic subsurface models have to be designed. Geophysical techniques can provide information for generating realistic models. The aim of our work is to show how ground-penetrating radar (GPR) investigations can contribute to such hydrological simulations. In the pilot area, Borkum island, GPR was used to map the shape of the groundwater table (GWT) and to characterise the aquifer. In total, 20 km of constant offset (CO) profiles were measured with centre frequencies of 80 and 200 MHz. Wave velocities were determined by common midpoint (CMP) measurements and vertical radar profiling (VRP) in a monitoring well. The 80 MHz CO data show a clear reflection at the groundwater table, whereas the reflection is weaker for the 200 MHz data. After correcting the GPR water tables for the capillary rise, they are in good accordance with the pressure heads of the observation wells in the area. In the centre of the island, the groundwater table is found up to 3.5 m above sea level, however it is lower towards the coastline and marshland. Some local depressions are observed in the region of dune valleys and around pumping stations of the local water supplier. GPR also reveals details within the sediments and highly-permeable aeolian sands can be distinguished from less-permeable marine sediments. Further, a silt loam layer below the water table could be mapped on a large area. The reflection characteristics indicates scattered erosion channels in this layer that cause it to be an aquitard with some leakage. GPR provides a high resolution map of the groundwater table and insight into the stratigraphy of the sediments and their hydraulic properties. This is valuable complementary information to the observation of sparsely distributed monitoring wells as input to hydraulic simulation.
  • Item
    Compiling geophysical and geological information into a 3-D model of the glacially-affected island of Föhr
    (Munich : EGU, 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.
  • Item
    Hydraulic properties at the North Sea island of Borkum derived from joint inversion of magnetic resonance and electrical resistivity soundings
    (Munich : EGU, 2012) Günther, T.; Müller-Petke, M.
    For reliably predicting the impact of climate changes on salt/freshwater systems below barrier islands, a long-term hydraulic modelling is inevitable. As input we need the parameters porosity, salinity and hydraulic conductivity at the catchment scale, preferably non-invasively acquired with geophysical methods. We present a methodology to retrieve the searched parameters and a lithological interpretation by the joint analysis of magnetic resonance soundings (MRS) and vertical electric soundings (VES). Both data sets are jointly inverted for resistivity, water content and decay time using a joint inversion scheme. Coupling is accomplished by common layer thicknesses. We show the results of three soundings measured on the eastern part of the North Sea island of Borkum. Pumping test data is used to calibrate the petrophysical relationship for the local conditions in order to estimate permeability from nuclear magnetic resonance (NMR) data. Salinity is retrieved from water content and resistivity using a modified Archie equation calibrated by local samples. As a result we are able to predict porosity, salinity and hydraulic conductivities of the aquifers, including their uncertainties. The joint inversion significantly improves the reliability of the results. Verification is given by comparison with a borehole. A sounding in the flooding area demonstrates that only the combined inversion provides a correct subsurface model. Thanks to the joint application, we are able to distinguish fluid conductivity from lithology and provide reliable hydraulic parameters as shown by uncertainty analysis. These findings can finally be used to build groundwater flow models for simulating climate changes. This includes the improved geometry and lithological attribution, and also the parameters and their uncertainties. © Author(s) 2012.
  • Item
    Evaluation of deep crustal earthquakes in northern Germany - Possible tectonic causes
    (Oxford [u.a.] : Wiley-Blackwell, 2019) Brandes, Christian; Plenefisch, Thomas; Tanner, David C.; Gestermann, Nicolai; Steffen, Holger
    We present new evidence for seven deep crustal, intraplate earthquakes in northern Germany, a region regarded as an area of low seismicity. From 2000 to 2018, seven earthquakes with magnitudes of ML 1.3–3.1, were detected at depths of 17.0–31.4 km. By placing the earthquake hypocentres in a geological three-dimensional model, we can correlate two of the earthquakes with the Thor Suture, a major fault zone in this area. Five of the earthquakes group in the lower crust near the Moho, which implies that parts of the lower crust and the crust/mantle boundary in northern Germany act as a structural discontinuity on which deformation localizes. Numerical simulation implies that stress changes due to glacial isostatic adjustment most likely triggered these deep crustal earthquakes.
  • Item
    Finite-difference modelling to evaluate seismic P-wave and shear-wave field data
    (Göttingen : Copernicus Publ., 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.
  • Item
    High-resolution shear-wave seismic reflection as a tool to image near-surface subrosion structures – a case study in Bad Frankenhausen, Germany
    (Göttingen : Copernicus Publ., 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.