Search

Instability of slopes, embankments, and cuttings on the railway network is increasingly prevalent globally. Monitoring vulnerable infrastructure aids in geotechnical asset management, and improvements to transport safety and efficiency. Here, we examine the use of a novel, near-real-time Electrical Resistivity Tomography (ERT) monitoring system for assessing the stability of a railway cutting in Leicestershire, United Kingdom. In 2015, an ERT monitoring system was installed across a relict landslide (grassed) and an area of more stable ground on either side (wooded), to monitor changes in electrical resistivity through time and space, and to assess the influence of different types of vegetation on the stability of transportation infrastructure. Two years of 4-Dimensional ERT monitoring results are presented here, and petrophysical relationships developed in the laboratory are applied to calibrate the resistivity models in order to provide an insight into hydrogeological pathways within a railway cutting. The influence of vegetation type on subsurface moisture pathways and on slope stability is also assessed – here we find that seasonal subsurface changes in moisture content and soil suction are exacerbated by the presence of trees (wooded area). This results in shrink-swell behaviour of the clays comprising the railway cutting, resulting in fissuring and a reduction in shear strength, leading to instability. As such, it is proposed that on slopes comprised of expansive soils, grassed slopes are beneficial for stability. Insights into the use of 4-D ERT for monitoring railway infrastructure gained from this study may be applied to the monitoring of critical geotechnical assets elsewhere.

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.

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.

Understanding ecologically sensitive wetlands often requires non‐invasive methods to characterize their complex structure (e.g., deposit heterogeneity) and hydrogeological parameters (e.g., porosity and hydraulic conductivity). Here, electrical conductivities of a riparian wetland were obtained using frequency domain electromagnetic induction (EMI) methods. The wetland was previously characterized by extensive intrusive measurements and 3D electrical resistivity tomography (ERT) surveys and hence offers an ideal opportunity to objectively assess EMI methods. Firstly, approaches to obtain structural information (e.g., elevation and thickness of alluvium) from EMI data and inverted models were assessed. Regularized and sharp inversion algorithms were investigated for ERT calibrated EMI data. Moreover, the importance of EMI errors in inversion was investigated. The hydrological information content was assessed using correlations with piezometric data and petrophysical models. It was found that EMI data were dominated by the thickness of peaty alluvial soils and relatively insensitive to topography and total alluvial thickness. Furthermore, although error weighting in the inversion improved the accuracy of alluvial soil thickness predictions, the multi‐linear regression method performed the best. For instance, an iso‐conductivity method to estimate the alluvial soil thickness in the regularized models had a normalized mean absolute difference (NMAD) of 21.4%, and although this performed better than the sharp inversion algorithm (NMAD = 65.3%), the multi‐linear regression approach (using 100 intrusive observations) achieved a NMAD = 18.0%. In terms of hydrological information content, correlations between EMI results and piezometric data were poor, however robust relationships between petrophysically derived porosity and hydraulic conductivity were observed for the alluvial soils and gravels. Key Points: Raw ECa values are highly correlated with the thickness of alluvial soil in a riparian wetlandAlluvial soil thickness predictions from multi‐linear regressions were more accurate than from electromagnetic induction (EMI) inversion methodsRobust predictions of hydraulic conductivity across the field site require more extensive intrusive data [ABSTRACT FROM AUTHOR]

Electrical resistivity tomography (ERT) is widely used to image the Earth's subsurface and has proven to be an extremely useful tool in application to hydrological problems. Conventional smoothness-constrained inversion of ERT data is efficient and robust, and consequently very popular. However, it does not resolve well sharp interfaces of a resistivity field and tends to reduce and smooth resistivity variations. These issues can be problematic in a range of hydrological or near-surface studies, for example mapping regolith-bedrock interfaces. While fully Bayesian approaches, such as those using Markov chain Monte Carlo sampling, can address the above issues, their very high computation cost makes them impractical for many applications. Ensemble Kalman inversion (EKI) offers a computationally efficient alternative by approximating the Bayesian posterior distribution in a derivative-free manner, which means only a relatively small number of 'black-box' model runs are required. Although common limitations for ensemble Kalman filter-type methods apply to EKI, it is both efficient and generally captures uncertainty patterns correctly. We propose the use of a new EKI-based framework for ERT which estimates a resistivity model and its uncertainty at a modest computational cost. Our EKI framework uses a level-set parametrization of the unknown resistivity to allow efficient estimation of discontinuous resistivity fields. Instead of estimating level-set parameters directly, we introduce a second step to characterize the spatial variability of the resistivity field and infer length scale hyperparameters directly. We demonstrate these features by applying the method to a series of synthetic and field examples. We also benchmark our results by comparing them to those obtained from standard smoothness-constrained inversion. Resultant resistivity images from EKI successfully capture arbitrarily shaped interfaces between resistivity zones and the inverted resistivities are close to the true values in synthetic cases. We highlight its readiness and applicability to similar problems in geophysics. [ABSTRACT FROM AUTHOR]

The Ripley Landslide is a small (0.04 km2), slow‐moving landslide in the Thompson River Valley, British Columbia, that is threatening the serviceability of two national railway lines. Slope failures in this area are having negative impacts on railway infrastructure, terrestrial and aquatic ecosystems, public safety, communities, local heritage and the economy. This is driving the need for monitoring at the site, and in recent years there has been a shift from traditional geotechnical surveys and visual inspections for monitoring infrastructure assets toward less invasive, lower cost, and less time‐intensive methods, including geophysics. We describe the application of a novel electrical resistivity tomography system for monitoring the landslide. The system provides near‐real time geoelectrical imaging, with results delivered remotely via a modem, avoiding the need for costly repeat field visits, and enabling near‐real time interpretation of the four‐dimensional electrical resistivity tomography data. Here, we present the results of the electrical resistivity tomography monitoring alongside field sensor‐derived relationships between suction, resistivity, moisture content and continuous monitoring single‐frequency Global Navigation Satellite System stations. Four‐dimensional electrical resistivity tomography data allows us to monitor spatial and temporal changes in resistivity, and by extension, in moisture content and soil suction. The models reveal complex hydrogeological pathways, as well as considerable seasonal variation in the response of the subsurface to changing weather conditions, which cannot be predicted through interrogation of weather and sensor data alone, providing new insight into the subsurface processes active at the site of the Ripley Landslide. [ABSTRACT FROM AUTHOR]

Adaptive optimal experimental design methods use previous data and results to guide the choice and design of future experiments. This paper describes the formulation of an adaptive survey design technique to produce optimal resistivity imaging surveys for time-lapse geoelectrical monitoring experiments. These survey designs are time-dependent and, compared to dipole-dipole or static optimized surveys that do not change over time, focus a greater degree of the image resolution on regions of the subsurface that are actively changing. The adaptive optimization method is validated using a controlled laboratory monitoring experiment comprising a well-defined cylindrical target moving along a trajectory that changes its depth and lateral position. The algorithm is implemented on a standard PC in conjunction with a modified automated multichannel resistivity imaging system. Data acquisition using the adaptive survey designs requires no more time or power than with comparable standard surveys, and the algorithm processing takes place while the system batteries recharge. The results show that adaptively designed optimal surveys yield a quantitative increase in image quality over and above that produced by using standard dipole-dipole or static (time-independent) optimized surveys. [ABSTRACT FROM AUTHOR]

Landslides in the Thompson River Valley, British Columbia, Canada, threaten the serviceability of two railway lines that connect Vancouver to the rest of Canada and the US. To minimise the impact of slope instability on vital transport infrastructure, as well as on terrestrial and aquatic ecosystems, public safety, communities, local heritage, and the economy, and to better inform decision making, there is a need for monitoring. Since 2013, the Ripley Landslide – a small, slow-moving, translational landslide – has been the focus of monitoring efforts in the Thompson River Valley transportation corridor. In November 2017, a novel Electrical Resistivity Tomography (ERT) monitoring system was installed on the site, providing near-real-time data collection via a telemetric link. 4-Dimensional resistivity models are presented in the context of moisture content and soil suction, two parameters known to influence slope stability in the Thompson River Valley. Here, we discuss the development of laboratory-based petrophysical relationships that relate electrical resistivity to moisture content and soil suction directly, building on relationships developed in the field. The 4-D ERT models were calibrated using these petrophysical relationships to provide insights into the complex spatial and temporal variations in moisture content and soil suction. This study highlights the utility of geoelectrical monitoring for assessing slope stability in the context of moisture-driven landslides. • Direct relationships between resistivity, suction, and moisture content are developed in the laboratory. • Calibration of ERT images using petrophysical relationships showss complex spatial and temporal changes in moisture content through space and time. • 4-D ERT provides insight into subsurface changes that cannot be predicted through rainfall data alone. [ABSTRACT FROM AUTHOR]

If electrodes move during geoelectrical resistivity monitoring and their new positions are not incorporated in the inversion, then the resulting tomographic images exhibit artefacts that can obscure genuine time-lapse resistivity changes in the subsurface. The effects of electrode movements on time-lapse resistivity tomography are investigated using a simple analytical model and real data. The correspondence between the model and the data is sufficiently good to be able to predict the effects of electrode movements with reasonable accuracy. For the linear electrode arrays and 2-D inversions under consideration, the data are much more sensitive to longitudinal than transverse or vertical movements. Consequently the model can be used to invert the longitudinal offsets of the electrodes from their known baseline positions using only the time-lapse ratios of the apparent resistivity data. The example data sets are taken from a permanently installed electrode array on an active lobe of a landslide. Using two sets with different levels of noise and subsurface resistivity changes, it is found that the electrode positions can be recovered to an accuracy of 4 per cent of the baseline electrode spacing. This is sufficient to correct the artefacts in the resistivity images, and provides for the possibility of monitoring the movement of the landslide and its internal hydraulic processes simultaneously using electrical resistivity tomography only. [ABSTRACT FROM AUTHOR]

The effects of geometric errors on crosshole resistivity data are investigated using analytical methods. Geometric errors are systematic and can occur due to uncertainties in the individual electrode positions, the vertical spacing between electrodes in the same borehole, or the vertical offset between electrodes in opposite boreholes. An estimate of the sensitivity to geometric error is calculated for each of two generic types of four-electrode crosshole configuration: current flow and potential difference crosshole (XH) and in-hole (IH). It is found that XH configurations are not particularly sensitive to geometric error unless the boreholes are closely spaced on the scale of the vertical separation of the current and potential electrodes. However, extremely sensitive IH configurations are shown to exist for any borehole separation. Therefore, it is recommended that XH configurations be used in preference to IH schemes. The effects of geometric error are demonstrated using real XH data from a closely spaced line of boreholes designed to monitor bioremediation of chlorinated solvents at an industrial site. A small fraction of the data had physically unrealistic apparent resistivities, which were either negative or unexpectedly large. However by filtering out configurations with high sensitivities to geometric error, all of the suspect data were removed. This filtering also significantly improved the convergence between the predicted and the measured resistivities when the data were inverted. In addition to systematic geometric errors, the measured data also exhibit a high level of random noise. Despite this, the resulting inverted images correspond reasonably closely with the known geology and nearby cone penetrometer resistivity profiles. [ABSTRACT FROM AUTHOR]