Electrical resistivity tomography (ERT) is increasingly being used to investigate unstable slopes and monitor the hydrogeological processes within. But movement of electrodes or incorrect placement of electrodes with respect to an assumed model can introduce significant resistivity artefacts into the reconstruction. In this work, we demonstrate a joint resistivity and electrode movement reconstruction algorithm within an iterative Gauss-Newton framework. We apply this to ERT monitoring data from an active slow-moving landslide in the UK. Results show fewer resistivity artefacts and suggest that electrode movement and resistivity can be reconstructed at the same time under certain conditions. A new 2.5-D formulation for the electrode position Jacobian is developed and is shown to give accurate numerical solutions when compared to the adjoint method on 3-D models. On large finite element meshes, the calculation time of the newly developed approach was also proven to be orders of magnitude faster than the 3-D adjoint method and addressed modelling errors in the 2-D perturbation and adjoint electrode position Jacobian.

Within geoelectrical imaging, the choice of measurement configurations and electrode locations is known to control the image resolution. Previous work has shown that optimized survey designs can provide a model resolution that is superior to standard survey designs. This paper demonstrates a methodology to optimize resolution within a target area, while limiting the number of required electrodes, thereby selecting optimal electrode locations. This is achieved by extending previous work on the 'Compare-R' algorithm, which by calculating updates to the resolution matrix optimizes the model resolution in a target area. Here, an additional weighting factor is introduced that allows to preferentially adding measurement configurations that can be acquired on a given set of electrodes. The performance of the optimization is tested on two synthetic examples and verified with a laboratory study. The effect of the weighting factor is investigated using an acquisition layout comprising a single line of electrodes. The results show that an increasing weight decreases the area of improved resolution, but leads to a smaller number of electrode positions. Imaging results superior to a standard survey design were achieved using 56 per cent fewer electrodes. The performance was also tested on a 3-D acquisition grid, where superior resolution within a target at the base of an embankment was achieved using 22 per cent fewer electrodes than a comparable standard survey. The effect of the underlying resistivity distribution on the performance of the optimization was investigated and it was shown that even strong resistivity contrasts only have minor impact. The synthetic results were verified in a laboratory tank experiment, where notable image improvements were achieved. This work shows that optimized surveys can be designed that have a resolution superior to standard survey designs, while requiring significantly fewer electrodes. This methodology thereby provides a means for improving the efficiency of geoelectrical imaging.

Measurement sensors permanently installed on landslides will inevitably change their position over time due to mass movements. To interpret and correct the recorded data, these movements have to be determined. This is especially important in the case of geoelectrical monitoring, where incorrect sensor positions produce strong artefacts in the resulting resistivity models. They may obscure real changes, which could indicate triggering mechanisms for landslide failure or reactivation. In this paper we introduce a methodology to interpolate movements from a small set of sparsely distributed reference points to a larger set of electrode locations. Within this methodology we compare three interpolation techniques, i.e., a piecewise planar, bi-linear spline, and a kriging based interpolation scheme. The performance of these techniques is tested on a synthetic and a real-data example, showing a recovery rate of true movements to about 1% and 10% of the electrode spacing, respectively. The significance for applying the proposed methodology is demonstrated by inverse modelling of 4D electrical resistivity tomography data, where it is shown that by correcting for sensor movements corresponding artefacts can virtually be removed and true resistivity changes be imaged.

Measurement errors can play a pivotal role in geophysical inversion. Most inverse models require users to prescribe or assume a statistical model of data errors before inversion. Wrongly prescribed errors can lead to over- or under-fitting of data; however, the derivation of models of data errors is often neglected. With the heightening interest in uncertainty estimation within hydrogeophysics, better characterisation and treatment of measurement errors is needed to provide improved image appraisal. Here we focus on the role of measurement errors in electrical resistivity tomography (ERT). We have analysed two time-lapse ERT datasets: one contains 96 sets of direct and reciprocal data collected from a surface ERT line within a 24 h timeframe; the other is a two-year-long cross-borehole survey at a UK nuclear site with 246 sets of over 50,000 measurements. Our study includes the characterisation of the spatial and temporal behaviour of measurement errors using autocorrelation and correlation coefficient analysis. We find that, in addition to well-known proportionality effects, ERT measurements can also be sensitive to the combination of electrodes used, i.e. errors may not be uncorrelated as often assumed. Based on these findings, we develop a new error model that allows grouping based on electrode number in addition to fitting a linear model to transfer resistance. The new model explains the observed measurement errors better and shows superior inversion results and uncertainty estimates in synthetic examples. It is robust, because it groups errors together based on the electrodes used to make the measurements. The new model can be readily applied to the diagonal data weighting matrix widely used in common inversion methods, as well as to the data covariance matrix in a Bayesian inversion framework. We demonstrate its application using extensive ERT monitoring datasets from the two aforementioned sites.

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.

This paper describes the Late Pleistocene to Holocene stratigraphy of the River Lambourn; a minor headwater of the River Thames in the Berkshire Downs. The Quaternary valley-fill comprises around 5-8. m of Late Pleistocene gravels overlain by Holocene peats and chalky clays. Quaternary deposits overlie an irregular rockhead erosion surface with deep scouring particularly evident on prominent bends in the valley. The gravels subdivide into a lower unit of chalky gravels overlain by coarse flint gravels. Ground penetrating radar suggests that gravels at depth are relatively structureless, but at the top show well-developed point-bar accretion surfaces which occur in association with peat-filled sinuous channels. These probably date from around the Pleistocene-Holocene boundary and may have formed in response to climate change and increased groundwater outflow as stream hydrology changed from the short-duration, high-magnitude flows of the Lower Dryas to the uniform, low-magnitude flows of the Holocene. Holocene peats initially infilled abandoned floodplain channels at around 10. kyr BP but later encroached over much of the Lambourn floodplain. A progressive upward decrease in organic material and an increase in the proportion of chalky clays from around 4. kyr BP probably occurred in response to floodplain accretion coupled with increased erosion of the chalk catchment related to agricultural clearance and a wetter climate.