Sebastian Uhlemann

Derivation of lowland riparian wetland deposit architecture using geophysical image analysis and interface detection

For groundwater-surface water interactions to be understood in complex wetland settings, the architecture of the underlying deposits requires investigation at a spatial resolution sufficient to characterize significant hydraulic pathways. Discrete intrusive sampling using conventional approaches provides insufficient sample density and can be difficult to deploy on soft ground. Here a noninvasive geophysical imaging approach combining three-dimensional electrical resistivity tomography (ERT) and the novel application of gradient and isosurface-based edge detectors is considered as a means of illuminating wetland deposit architecture. The performance of three edge detectors were compared and evaluated against ground truth data, using a lowland riparian wetland demonstration site. Isosurface-based methods correlated well with intrusive data and were useful for defining the geometries of key geological interfaces (i.e., peat/gravels and gravels/Chalk). The use of gradient detectors approach was unsuccessful, indicating that the assumption that the steepest resistivity gradient coincides with the associated geological interface can be incorrect. These findings are relevant to the application of this approach in settings with a broadly layered geology with strata of contrasting resistivities. In addition, ERT revealed substantial structures in the gravels related to the depositional environment (i.e., braided fluvial system) and a complex distribution of low-permeability putty Chalk at the bedrock surface—with implications for preferential flow and variable exchange between river and groundwater systems. These results demonstrate that a combined approach using ERT and edge detectors can provide valuable information to support targeted monitoring and inform hydrological modeling of wetlands.

Integrated time-lapse geoelectrical imaging of wetland hydrological processes

Wetlands provide crucial habitats, are critical in the global carbon cycle, and act as key biogeochemical and hydrological buffers. The effectiveness of these services is mainly controlled by hydrological processes, which can be highly variable both spatially and temporally due to structural complexity and seasonality. Spatial analysis of 2D geoelectrical monitoring data integrated into the interpretation of conventional hydrological data has been implemented to provide a detailed understanding of hydrological processes in a riparian wetland. This study shows that a combination of processes can define the resistivity signature of the shallow subsurface, highlighting the seasonality of these processes and its corresponding effect on biogeochemical processesthe wetland hydrology. Groundwater exchange between peat and the underlying river terrace deposits, spatially and temporally defined by geoelectrical imaging and verified by point sensor data, highlighted the groundwater dependent nature of the wetland. A 30 % increase in peat resistivity was shown to be caused by a nearly entire exchange of the saturating groundwater. For the first time, we showed that automated interpretation of geoelectrical data can be used to quantify shrink-swell of expandable soils, affecting hydrological parameters, such as, porosity, water storage capacity, and permeability. This study shows that an integrated interpretation of hydrological and geophysical data can significantly improve the understanding of wetland hydrological processes. Potentially, this approach can provide the basis for the evaluation of ecosystem services and may aid in the optimization of wetland management strategies.

Reconstruction of landslide movements by inversion of 4‐D electrical resistivity tomography monitoring data

Reliable tomographic inversion of geoelectrical monitoring data from unstable slopes relies critically on knowing the electrode positions, which may move over time. We develop and present an innovative inverse method to recover movements in both surface directions from geoelectrical measurements made on a grid of monitoring electrodes. For the first time, we demonstrate this method using field data from an active landslide to recover sequences of movement over timescales of days to years. Comparison with GPS measurements demonstrated an accuracy of within 10 % of the electrode spacing, sufficient to correct the majority of artefacts that would occur in subsequent image reconstructions if incorrect positions are used. Over short timescales where the corresponding subsurface resistivity changes were smaller, the constraints could be relaxed and an order-of-magnitude better accuracy was achievable. This enabled the onset and acceleration of landslide activity to be detected with a temporal resolution of a few days.

Four-dimensional imaging of moisture dynamics during landslide reactivation

Landslides pose significant risks to communities and infrastructure, and mitigating these risks relies on understanding landslide causes and triggering processes. It has been shown that geophysical surveys can significantly contribute to the characterization of unstable slopes. However, hydrological processes can be temporally and spatially heterogeneous, requiring their related properties to be monitored over time. Geoelectrical monitoring can provide temporal and volumetric distributions of electrical resistivity, which are directly related to moisture content. To date, studies demonstrating this capability have been restricted to 2D sections, which are insufficient to capture the full degree of spatial heterogeneity. This study is the first to employ 4D (i.e., 3D time-lapse) resistivity imaging on an active landslide, providing long-term data (three years) highlighting the evolution of moisture content prior to landslide reactivation and showing its decline post reactivation. Crucially the time-lapse inversion methodology employed here incorporates movements of the electrodes on the unstable surface. Although seasonal characteristics dominate the shallow moisture dynamics during the first two years with surficial drying in summer and wetting in winter, in the months preceding reactivation, moisture content increased by more than 45 % throughout the slope. This is in agreement with independent data showing a significant rise in piezometric heads and shallow soil moisture contents as a result of prolonged and intense rainfall. Based on these results, remediation measures could be designed and early-warning systems implemented. Thus, resistivity monitoring that can allow for moving electrodes provides a new means for the effective mitigation of landslide risk.

Geoelectrical monitoring of simulated subsurface leakage to support high-hazard nuclear decommissioning at the Sellafield Site, UK

A full-scale field experiment applying 4D (3D time-lapse) cross-borehole Electrical Resistivity Tomography (ERT) to the monitoring of simulated subsurface leakage was undertaken at a legacy nuclear waste silo at the Sellafield Site, UK. The experiment constituted the first application of geoelectrical monitoring in support of decommissioning work at a UK nuclear licensed site. Images of resistivity changes occurring since a baseline date prior to the simulated leaks revealed likely preferential pathways of silo liquor simulant flow in the vadose zone and upper groundwater system. Geophysical evidence was found to be compatible with historic contamination detected in permeable facies in sediment cores retrieved from the ERT boreholes. Results indicate that laterally discontinuous till units forming localized hydraulic barriers substantially affect flow patterns and contaminant transport in the shallow subsurface at Sellafield. We conclude that only geophysical imaging of the kind presented here has the potential to provide the detailed spatial and temporal information at the (sub-)meter scale needed to reduce the uncertainty in models of subsurface processes at nuclear sites.

Monitoring rock freezing and thawing by novel geoelectrical and acoustic techniques

Automated monitoring of freeze-thaw cycles and fracture propagation in mountain rockwalls is needed to provide early warning about rockfall hazards. Conventional geoelectrical methods such as electrical resistivity tomography (ERT) are limited by large and variable ohmic contact resistances, requiring galvanic coupling with metal electrodes inserted into holes drilled into rock, and which can be loosened by rock weathering. We report a novel experimental methodology that combined capacitive resistivity imaging (CRI), ERT, and microseismic event recording to monitor freeze-thaw of six blocks of hard and soft limestones under conditions simulating an active layer above permafrost and seasonally frozen rock in a nonpermafrost environment. Our results demonstrate that the CRI method is highly sensitive to freeze-thaw processes; it yields property information equivalent to that obtained with conventional ERT and offers a viable route for nongalvanic long-term geoelectrical monitoring, extending the benefits of the methodology to soft/hard rock environments. Contact impedances achieved with CRI are less affected by seasonal temperature changes, the aggregate state of the pore water (liquid or frozen), and the presence of low-porosity rock with high matrix resistivities than those achieved with ERT. Microseismic monitoring has the advantage over acoustic emissions that events were recorded in relevant field distances of meters to decameters from cracking events. For the first time we recorded about 1000 microcracking events and clustered them in four groups according to frequency and waveform. Compared to previous studies, mainly on ice-cracking in glaciers, the groups are attributed to single- or multiple-stage cracking events such as crack coalescence.

High resolution profile of inorganic aqueous geochemistry and key redox zones in an arsenic bearing aquifer in Cambodia

Arsenic contamination of groundwaters in South and Southeast Asia is a major threat to public health. In order to better understand the geochemical controls on the mobility of arsenic in a heavily arsenic-affected aquifer in northern Kandal Province, Cambodia, key changes in inorganic aqueous geochemistry have been monitored at high vertical and lateral resolution along dominant groundwater flow paths along two distinct transects. The two transects are characterized by differing geochemical, hydrological and lithological conditions. Arsenic concentrations in groundwater are highly heterogenous, and are broadly positively associated with iron and negatively associated with sulfate and dissolved oxygen. The observed correlations are generally consistent with arsenic mobilization by reductive-dissolution of iron (hydr)oxides. Key redox zones, as identified using groupings of the PHREEQC model equilibrium electron activity of major redox couples (notably ammonium/nitrite; ammonium/nitrate; nitrite/nitrate; dissolved oxygen/water) have been identified and vary with depth, site and season. Mineral saturation is also characterized. Seasonal changes in groundwater chemistry were observed in areas which were (i) sandy and of high permeability; (ii) in close proximity to rivers; and/or (iii) in close proximity to ponds. Such changes are attributed to monsoonal-driven surface-groundwater interactions and are consistent with the separate provenance of recharge sources as identified using stable isotope mixing models.

Optimised sequential experimental design for Geoelectrical Resistivity Monitoring Surveys

Sequential experimental design methods use previous data and results to guide the choice and design of future experiments. This paper describes the application of a sequential design technique to produce optimal resistivity imaging surveys for time-lapse geoelectrical monitoring experiments. These survey designs are time-dependent, and are optimised to focus a greater degree of the image resolution on the regions of the subsurface that are actively changing than static optimised surveys that do not change over time. The sequential design method is applied to a synthetic 2.5D monitoring experiment comprising a well-defined cylindrical target moving along a trajectory that changes its depth and lateral position. The data are simulated to be as realistic as possible, incorporating survey design constraints for a real resistivity monitoring system and realistic levels and distributions of random noise, in order to match a forthcoming experimental test of the method. The results of the simulations indicate that sequentially designed optimal surveys yield an increase in image quality over and above that produced by using a static (time-independent) optimised survey.

Current and future role of instrumentation and monitoring in the performance of transport infrastructure slopes

Instrumentation is often used to monitor the performance of engineered infrastructure slopes. This paper looks at the current role of instrumentation and monitoring, including the reasons for monitoring infrastructure slopes, the instrumentation typically installed and parameters measured. The paper then investigates recent developments in technology and considers how these may change the way that monitoring is used in the future, and tries to summarize the barriers and challenges to greater use of instrumentation in slope engineering. The challenges relate to economics of instrumentation within a wider risk management system, a better understanding of the way in which slopes perform and/or lose performance, and the complexities of managing and making decisions from greater quantities of data.

Numerical Simulations of Capacitive Resistivity Imaging (CRI) Measurements

Electrical resistivity tomography (ERT) is a well-developed geophysical technique that is used to study a variety of geoscientific problems. In recent years it has been applied to the study of permafrost processes at both field and laboratory scale. However, highly resistive surface conditions limit its applicability due to high and variable contact resistances. The use of capacitively coupled sensors is expected to overcome this problem by providing a steady contact impedance regime. Although the theory of capacitive resistivity imaging (CRI) is well understood, a point-pole approximation of the sensors is typically used to show the similarity between CRI and ERT. Due to their nature, capacitive sensors cannot be designed as point-poles as they require a finite extent. This paper assesses the effects the finite size of sensors has on the applicability of CRI theory and aims to provide an improved understanding of the measured data. We employ finite-element numerical modelling to simulate CRI measurements over a homogeneous halfspace and on a finite rock sample. The results of a parameter study over a homogeneous halfspace are compared to an analytical solution. Observed discrepancies between the two solutions clearly indicate that large sensor sizes and small sensor separations violate the point-pole assumption of the analytical solution. In terms of data interpretation, this dictates that sensor separations smaller than twice the sensor size have to be avoided in order to remain below a generic error threshold of 5%. We show that sensor elevation, halfspace resistivity, halfspace permittivity, and measurement frequency have only minor effects on the discrepancy between simulation and analytical solution. The simulation of sequential CRI measurements on a finite rock sample suggests that, in line with expectations, the measured signals lie mainly in the 4th quadrant of the complex plane. However, we can also observe data with negative geometric factors, which are related to uncommon array. A comparison between simulated and measured data showed very good agreement; it validated the simulations and explained the measured data acquired using a prototype multisensor CRI system. We show that a comparison of simulated and measured imaginary parts of the transfer impedance can be used to assess CRI measurement errors. Our work demonstrates that finite-element numerical modelling of CRI measurements is a valuable tool with which to define limitations on array design and to assess data quality.