2015, Hübner, R., Heller, K., Günther, T., Kleber, A.

Besides floodplains, hillslopes are basic units that mainly control water movement and flow pathways within catchments of subdued mountain ranges. The structure of their shallow subsurface affects water balance, e.g. infiltration, retention, and runoff. Nevertheless, there is still a gap in the knowledge of the hydrological dynamics on hillslopes, notably due to the lack of generalization and transferability. This study presents a robust multi-method framework of electrical resistivity tomography (ERT) in addition to hydrometric point measurements, transferring hydrometric data into higher spatial scales to obtain additional patterns of distribution and dynamics of soil moisture on a hillslope. A geoelectrical monitoring in a small catchment in the eastern Ore Mountains was carried out at weekly intervals from May to December 2008 to image seasonal moisture dynamics on the hillslope scale. To link water content and electrical resistivity, the parameters of Archie's law were determined using different core samples. To optimize inversion parameters and methods, the derived spatial and temporal water content distribution was compared to tensiometer data. The results from ERT measurements show a strong correlation with the hydrometric data. The response is congruent to the soil tension data. Water content calculated from the ERT profile shows similar variations as that of water content from soil moisture sensors. Consequently, soil moisture dynamics on the hillslope scale may be determined not only by expensive invasive punctual hydrometric measurements, but also by minimally invasive time-lapse ERT, provided that pedo-/petrophysical relationships are known. Since ERT integrates larger spatial scales, a combination with hydrometric point measurements improves the understanding of the ongoing hydrological processes and better suits identification of heterogeneities.

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, 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, 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.

2013, Beff, L., Günther, T., Vandoorne, B., Couvreur, V., Javaux, M.

A good understanding of the soil water content (SWC) distribution at the field scale is essential to improve the management of water, soil and crops. Recent studies proved that Electrical Resistivity Tomography (ERT) opens interesting perspectives in the determination of the SWC distribution in 3 dimensions (3-D). This study was conducted (i) to check and validate how ERT is able to monitor SWC distribution in a maize field during the late growing season; and (ii) to investigate how maize plants and rainfall affect the dynamics of SWC distribution. Time Domain Reflectometry (TDR) measurements were used to validate ERT-inverted SWC values. Evolution of water mass balance was also calculated to check whether ERT was capable of giving a reliable estimate of soil water stock evolution. It is observed that ERT was able to give the same average SWC as TDR (R2 = 0.98). In addition, ERT gives better estimates of the water stock than TDR thanks to its higher spatial resolution. The high resolution of ERT measurements also allows for the discrimination of SWC heterogeneities. The SWC distribution showed that alternation of maize rows and inter-rows was the main influencing factor of the SWC distribution. The drying patterns were linked to the root profiles, with drier zones under the maize rows. During short periods, with negligible rainfall, the SWC decrease took place mainly in the two upper soil horizons and in the inter-row area. In contrast, rainfall increased the SWC mostly under the maize rows and in the upper soil layer. Nevertheless, the total amount of rainfall during the growing season was not sufficient to modify the SWC patterns induced by the maize rows. During the experimental time, there was hardly any SWC redistribution from maize rows to inter-rows. Yet, lateral redistribution from inter-rows to maize rows induced by potential gradient generates SWC decrease in the inter-row area and in the deeper soil horizons.

2012, Garré, S., Günther, T., Diels, J., Vanderborght, J.

Contour hedgerow intercropping systems have been proposed as an alternative to traditional agricultural practice with a single crop, as they are effective in reducing run-off and soil erosion. However, competition for water and nutrients between crops and associated hedgerows may reduce the overall performance of these systems. To get a more detailed understanding of the competition for water, spatially resolved monitoring of soil water contents in the soil-plant-atmosphere system is necessary. Electrical resistivity tomography (ERT) is potentially a valuable technique to monitor changes in soil moisture in space and time. In this study, the performance of different ERT electrode arrays to detect the soil moisture dynamics in a mono- and an intercropping system was tested. Their performance was analyzed based on a synthetic study using geophysical measures, such as data recovery and resolution, and using spatial statistics of retrieved water content, such as an adjusted coefficient of variation and semivariances. The synthetic ERT measurements detected differences between the cropping systems and retrieved spatial structure of the soil moisture distribution, but the variance and semivariance were underestimated. Sharp water content contrasts between horizons or in the neighborhood of a root water uptake bulb were smoothened. The addition of electrodes deeper in the soil improved the performance, but sometimes only marginally. ERT is therefore a valuable tool for soil moisture monitoring in the field under different cropping systems if an electrode array is used which can resolve the patterns expected to be present in the medium. The use of spatial statistics allowed to not only identify the overall model recovery, but also to quantify the recovery of spatial structures.