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.