Search Results

Now showing 1 - 6 of 6
  • Item
    Impacts of enhanced weathering on biomass production for negative emission technologies and soil hydrology
    (Katlenburg-Lindau [u.a.] : Copernicus, 2020) De Oliveira Garcia, Wagner; Amann, Thorben; Hartmann, Jens; Karstens, Kristine; Popp, Alexander; Boysen, Lena R.; Smith, Pete; Goll, Daniel
    Limiting global mean temperature changes to well below 2 °C likely requires a rapid and large-scale deployment of negative emission technologies (NETs). Assessments so far have shown a high potential of biomass-based terrestrial NETs, but only a few assessments have included effects of the commonly found nutrient-deficient soils on biomass production. Here, we investigate the deployment of enhanced weathering (EW) to supply nutrients to areas of afforestation-reforestation and naturally growing forests (AR) and bioenergy grasses (BG) that are deficient in phosphorus (P), besides the impacts on soil hydrology. Using stoichiometric ratios and biomass estimates from two established vegetation models, we calculated the nutrient demand of AR and BG. Insufficient geogenic P supply limits C storage in biomass. For a mean P demand by AR and a lowgeogenic-P-supply scenario, AR would sequester 119 Gt C in biomass; for a high-geogenic-P-supply and low-AR-Pdemand scenario, 187 Gt C would be sequestered in biomass; and for a low geogenic P supply and high AR P demand, only 92 GtC would be accumulated by biomass. An average amount of ∼ 150 Gt basalt powder applied for EW would be needed to close global P gaps and completely sequester projected amounts of 190 Gt C during the years 2006-2099 for the mean AR P demand scenario (2-362 Gt basalt powder for the low-AR-P-demand and for the high-AR-P-demand scenarios would be necessary, respectively). The average potential of carbon sequestration by EW until 2099 is ∼ 12 GtC (∼ 0:2-∼ 27 Gt C) for the specified scenarios (excluding additional carbon sequestration via alkalinity production). For BG, 8 kg basaltm2 a1 might, on average, replenish the exported potassium (K) and P by harvest. Using pedotransfer functions, we show that the impacts of basalt powder application on soil hydraulic conductivity and plant-Available water, to close predicted P gaps, would depend on basalt and soil texture, but in general the impacts are marginal. We show that EW could potentially close the projected P gaps of an AR scenario and nutrients exported by BG harvest, which would decrease or replace the use of industrial fertilizers. Besides that, EW ameliorates the soil's capacity to retain nutrients and soil pH and replenish soil nutrient pools. Lastly, EW application could improve plant-Available-water capacity depending on deployed amounts of rock powder - adding a new dimension to the coupling of land-based biomass NETs with EW. © 2020 Royal Society of Chemistry. All rights reserved.
  • Item
    Combined rock slope stability and shallow landslide susceptibility assessment of the Jasmund cliff area (Rügen Island, Germany)
    (Katlenburg-Lindau : European Geophysical Society, 2009-5-8) Günther, A.; Thiel, C.
    In this contribution we evaluated both the structurally-controlled failure susceptibility of the fractured Cretaceous chalk rocks and the topographically-controlled shallow landslide susceptibility of the overlying glacial sediments for the Jasmund cliff area on Rügen Island, Germany. We employed a combined methodology involving spatially distributed kinematical rock slope failure testing with tectonic fabric data, and both physically- and inventory-based shallow landslide susceptibility analysis. The rock slope failure susceptibility model identifies areas of recent cliff collapses, confirming its value in predicting the locations of future failures. The model reveals that toppling is the most important failure type in the Cretaceous chalk rocks of the area. The shallow landslide susceptibility analysis involves a physically-based slope stability evaluation which utilizes material strength and hydraulic conductivity data, and a bivariate landslide susceptibility analysis exploiting landslide inventory data and thematic information on ground conditioning factors. Both models show reasonable success rates when evaluated with the available inventory data, and an attempt was made to combine the individual models to prepare a map displaying both terrain instability and landslide susceptibility. This combination highlights unstable cliff portions lacking discrete landslide areas as well as cliff sections highly affected by past landslide events. Through a spatial integration of the rock slope failure susceptibility model with the combined shallow landslide assessment we produced a comprehensive landslide susceptibility map for the Jasmund cliff area.
  • Item
    Hydraulic characterisation of iron-oxide-coated sand and gravel based on nuclear magnetic resonance relaxation mode analyses
    (Munich : EGU, 2018) Costabel, Stephan; Weidner, Christoph; Müller-Petke, Mike; Houben, Georg
    The capability of nuclear magnetic resonance (NMR) relaxometry to characterise hydraulic properties of iron-oxide-coated sand and gravel was evaluated in a laboratory study. Past studies have shown that the presence of paramagnetic iron oxides and large pores in coarse sand and gravel disturbs the otherwise linear relationship between relaxation time and pore size. Consequently, the commonly applied empirical approaches fail when deriving hydraulic quantities from NMR parameters. Recent research demonstrates that higher relaxation modes must be taken into account to relate the size of a large pore to its NMR relaxation behaviour in the presence of significant paramagnetic impurities at its pore wall. We performed NMR relaxation experiments with water-saturated natural and reworked sands and gravels, coated with natural and synthetic ferric oxides (goethite, ferrihydrite), and show that the impact of the higher relaxation modes increases significantly with increasing iron content. Since the investigated materials exhibit narrow pore size distributions, and can thus be described by a virtual bundle of capillaries with identical apparent pore radius, recently presented inversion approaches allow for estimation of a unique solution yielding the apparent capillary radius from the NMR data. We found the NMR-based apparent radii to correspond well to the effective hydraulic radii estimated from the grain size distributions of the samples for the entire range of observed iron contents. Consequently, they can be used to estimate the hydraulic conductivity using the well-known Kozeny–Carman equation without any calibration that is otherwise necessary when predicting hydraulic conductivities from NMR data. Our future research will focus on the development of relaxation time models that consider pore size distributions. Furthermore, we plan to establish a measurement system based on borehole NMR for localising iron clogging and controlling its remediation in the gravel pack of groundwater wells.
  • 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
    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
    Two-dimensional numerical investigations on the termination of bilinear flow in fractures
    (Göttingen : Copernicus Publ., 2013) Ortiz R., A.E.; Jung, R.; Renner, J.
    Bilinear flow occurs when fluid is drained from a permeable matrix by producing it through an enclosed fracture of finite conductivity intersecting a well along its axis. The terminology reflects the combination of two approximately linear flow regimes: one in the matrix with flow essentially perpendicular to the fracture, and one along the fracture itself associated with the non-negligible pressure drop in it. We investigated the characteristics, in particular the termination, of bilinear flow by numerical modeling allowing for an examination of the entire flow field without prescribing the flow geometry in the matrix. Fracture storage capacity was neglected relying on previous findings that bilinear flow is associated with a quasi-steady flow in the fracture. Numerical results were generalized by dimensionless presentation. Definition of a dimensionless time that, other than in previous approaches, does not use geometrical parameters of the fracture permitted identifying the dimensionless well pressure for the infinitely long fracture as the master curve for type curves of all fractures with finite length from the beginning of bilinear flow up to fully developed radial flow. In log–log scale the master curve's logarithmic derivative initially follows a 1/4-slope straight line (characteristic for bilinear flow) and gradually bends into a horizontal line (characteristic for radial flow) for long times. During the bilinear flow period, isobars normalized to well pressure propagate with the fourth and second root of time in fracture and matrix, respectively. The width-to-length ratio of the pressure field increases proportional to the fourth root of time during the bilinear period, and starts to deviate from this relation close to the deviation of well pressure and its derivative from their fourth-root-of-time relations. At this time, isobars are already significantly inclined with respect to the fracture. The type curves of finite fractures all deviate counterclockwise from the master curve instead of clockwise or counterclockwise from the 1/4-slope straight line as previously proposed. The counterclockwise deviation from the master curve was identified as the arrival of a normalized isobar reflected at the fracture tip 16 times earlier. Nevertheless, two distinct regimes were found in regard to pressure at the fracture tip when bilinear flow ends. For dimensionless fracture conductivities TD < 1, a significant pressure increase is not observed at the fracture tip until bilinear flow is succeeded by radial flow at a fixed dimensionless time. For TD > 10, the pressure at the fracture tip has reached substantial fractions of the associated change in well pressure when the flow field transforms towards intermittent formation linear flow at times that scale inversely with the fourth power of dimensionless fracture conductivity. Our results suggest that semi-log plots of normalized well pressure provide a means for the determination of hydraulic parameters of fracture and matrix after shorter test duration than for conventional analysis.