2017, Hübner, Rico, Günther, Thomas, Heller, Katja, Noell, Ursula, Kleber, Arno

Identifying principles of water movement in the shallow subsurface is crucial for adequate process-based hydrological models. Hillslopes are the essential interface for water movement in catchments. The shallow subsurface on slopes typically consists of different layers with varying characteristics. The aim of this study was to draw conclusions about the infiltration behaviour, to identify water flow pathways and derive some general interpretations for the validity of the water movement on a hillslope with periglacial slope deposits (cover beds), where the layers differ in their sedimentological and hydrological properties. Especially the described varying influence of the basal layer (LB) as an impeding layer on the one hand and as a remarkable pathway for rapid subsurface stormflow on the other. We used a time lapse 3-D electrical resistivity tomography (ERT) approach combined with punctual hydrometric data to trace the spreading and the progression of an irrigation plume in layered slope deposits during two irrigation experiments. This multi-technical approach enables us to connect the high spatial resolution of the 3-D ERT with the high temporal resolution of the hydrometric devices. Infiltration through the uppermost layer was dominated by preferential flow, whereas the water flow in the deeper layers was mainly matrix flow. Subsurface stormflow due to impeding characteristic of the underlying layer occurs in form of "organic layer interflow" and at the interface to the first basal layer (LB1). However, the main driving factor for subsurface stormflow is the formation of a capillary barrier at the interface to the second basal layer (LB2). The capillary barrier prevents water from entering the deeper layer under unsaturated conditions and diverts the seepage water according to the slope inclination. With higher saturation, the capillary barrier breaks down and water reaches the highly conductive deeper layer. This highlights the importance of the capillary barrier effect for the prevention or activation of different flow pathways under variable hydrological conditions.

2017, Rücker, Carsten, Günther, Thomas, Wagner, Florian M.

Many tasks in applied geosciences cannot be solved by single measurements, but require the integration of geophysical, geotechnical and hydrological methods. Numerical simulation techniques are essential both for planning and interpretation, as well as for the process understanding of modern geophysical methods. These trends encourage open, simple, and modern software architectures aiming at a uniform interface for interdisciplinary and flexible modelling and inversion approaches. We present pyGIMLi (Python Library for Inversion and Modelling in Geophysics), an open-source framework that provides tools for modelling and inversion of various geophysical but also hydrological methods. The modelling component supplies discretization management and the numerical basis for finite-element and finite-volume solvers in 1D, 2D and 3D on arbitrarily structured meshes. The generalized inversion framework solves the minimization problem with a Gauss-Newton algorithm for any physical forward operator and provides opportunities for uncertainty and resolution analyses. More general requirements, such as flexible regularization strategies, time-lapse processing and different sorts of coupling individual methods are provided independently of the actual methods used. The usage of pyGIMLi is first demonstrated by solving the steady-state heat equation, followed by a demonstration of more complex capabilities for the combination of different geophysical data sets. A fully coupled hydrogeophysical inversion of electrical resistivity tomography (ERT) data of a simulated tracer experiment is presented that allows to directly reconstruct the underlying hydraulic conductivity distribution of the aquifer. Another example demonstrates the improvement of jointly inverting ERT and ultrasonic data with respect to saturation by a new approach that incorporates petrophysical relations in the inversion. Potential applications of the presented framework are manifold and include time-lapse, constrained, joint, and coupled inversions of various geophysical and hydrological data sets.

2017, Ronczka, Mathias, Hellman, Kristofer, Günther, Thomas, Wisén, Roger, Dahlin, Torleif

Tunnelling below water passages is a challenging task in terms of planning, pre-investigation and construction. Fracture zones in the underlying bedrock lead to low rock quality and thus reduced stability. For natural reasons, they tend to be more frequent at water passages. Ground investigations that provide information on the subsurface are necessary prior to the construction phase, but these can be logistically difficult. Geophysics can help close the gaps between local point information by producing subsurface images. An approach that combines seismic refraction tomography and electrical resistivity tomography has been tested at the Äspö Hard Rock Laboratory (HRL). The aim was to detect fracture zones in a well-known but logistically challenging area from a measuring perspective. The presented surveys cover a water passage along part of a tunnel that connects surface facilities with an underground test laboratory. The tunnel is approximately 100 m below and 20 m east of the survey line and gives evidence for one major and several minor fracture zones. The geological and general test site conditions, e.g. with strong power line noise from the nearby nuclear power plant, are challenging for geophysical measurements. Co-located positions for seismic and ERT sensors and source positions are used on the 450 m underwater section of the 700 m profile. Because of a large transition zone that appeared in the ERT result and the missing coverage of the seismic data, fracture zones at the southern and northern parts of the underwater passage cannot be detected by separated inversion. Synthetic studies show that significant three-dimensional (3-D) artefacts occur in the ERT model that even exceed the positioning errors of underwater electrodes. The model coverage is closely connected to the resolution and can be used to display the model uncertainty by introducing thresholds to fade-out regions of medium and low resolution. A structural coupling cooperative inversion approach is able to image the northern fracture zone successfully. In addition, previously unknown sedimentary deposits with a significantly large thickness are detected in the otherwise unusually well-documented geological environment. The results significantly improve the imaging of some geologic features, which would have been undetected or misinterpreted otherwise, and combines the images by means of cluster analysis into a conceptual subsurface model.