Implementing positivity constraints in 4-D resistivity time-lapse inversion

Over the last 25 years 2-D and 3-D resistivity surveys have been used for a wide range of engineering, environmental, hydrological and mineral exploration surveys (Loke et al. 2013). In some surveys, the purpose includes the monitoring of subsurface changes with time (Chambers et al. 2014). The 4-D smoothness-constrained inversion method (Loke et al. 2014) has proved to be a stable and robust method for the inversion of time-lapse data sets. This method inverts the data sets measured at different times simultaneously and it includes a temporal smoothness constraint to ensure that the resistivity changes in a smooth manner with time. In some surveys, such as infiltration experiments (Kuras et al. 2016), it is known that the subsurface resistivity should only decrease (or increase) with time. As the standard 4-D inversion method does not explicitly constrain the direction of the changes with time, this could result in artefacts where an increase in the resistivity is obtained in the inverse model while it is only expected to decrease (or vice versa). In this paper we describe a modification of the 4-D smoothness-constrained inversion method to remove such temporal artefacts

Landslide characterization using P- and S-wave seismic refraction tomography — The importance of elastic moduli

© 2016 In the broad spectrum of natural hazards, landslides in particular are capable of changing the landscape and causing significant human and economic losses. Detailed site investigations form an important component in the landslide risk mitigation and disaster risk reduction process. These investigations usually rely on surface observations, discrete sampling of the subsurface, and laboratory testing to examine properties that are deemed representative of entire slopes. Often this requires extensive interpolations and results in large uncertainties. To compliment and extend these approaches, we present a study from an active landslide in a Lias Group clay slope, North Yorkshire, UK, examining combined P- and S-wave seismic refraction tomography (SRT) as a means of providing subsurface volumetric imaging of geotechnical proxies. The distributions of seismic wave velocities determined from SRT at the study site indicated zones with higher porosity and fissure density that are interpreted to represent the extent and depth of mass movements and weathered bedrock zones. Distinguishing the lithological units was facilitated by deriving the Poisson's ratio from the SRT data as saturated clay and partially saturated sandy silts showed distinctively different Poisson's ratios. Shear and Young's moduli derived from the SRT data revealed the weak nature of the materials in active parts of the landslide (i.e. 25 kPa and 100 kPa respectively). The SRT results are consistent with intrusive (i.e. cone penetration tests), laboratory, and additional geoelectrical data from this site. This study shows that SRT forms a cost-effective method that can significantly reduce uncertainties in the conceptual ground model of geotechnical and hydrological conditions that govern landslide dynamics

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 2-D sections, which are insufficient to capture the full degree of spatial heterogeneity. This study is the first to employ 4-D (i.e., 3-D time lapse) resistivity imaging on an active landslide, providing long-term data (3 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 2 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

Assessment of ground-based monitoring techniques applied to landslide investigations

A landslide complex in the Whitby Mudstone Formation at Hollin Hill, North Yorkshire, UK is periodically re-activated in response to rainfall-induced pore-water pressure fluctuations. This paper compares long-term measurements (i.e., 2009 – 2014) obtained from a combination of monitoring techniques that have been employed together for the first time on an active landslide. The results highlight the relative performance of the different techniques, and can provide guidance for researchers and practitioners for selecting and installing appropriate monitoring techniques to assess unstable slopes. Particular attention is given to the spatial and temporal resolution offered by the different approaches that include: Real Time Kinematic-GPS (RTK-GPS) monitoring of a ground surface marker array, conventional inclinometers, Shape Acceleration Arrays (SAA), tilt meters, active waveguides with Acoustic Emission (AE) monitoring, and piezometers. High spatial resolution information has allowed locating areas of stability and instability across a large slope. This has enabled identification of areas where further monitoring efforts should be focused. High temporal resolution information allowed the capture of S’-shaped slope displacement-time behaviour (i.e. phases of slope acceleration, deceleration and stability) in response to elevations in pore-water pressures. This study shows that a well-balanced suite of monitoring techniques that provides high temporal and spatial resolution on both measurement and slope scale is necessary to fully understand failure and movement mechanisms of slopes. In the case of the Hollin Hill landslide it enabled detailed interpretation of the geomorphological processes governing landslide activity. It highlights the benefit of regularly surveying a network of GPS markers to determine areas for installation of movement monitoring techniques that offer higher resolution both temporally and spatially. The small sensitivity of tilt meter measurements to translational movements limited the ability to record characteristic ‘S’-shaped landslide movements at Hollin Hill, which were identified using SAA and AE measurements. This high sensitivity to landslide movements indicates the applicability of SAA and AE monitoring to be used in early warning systems, through detecting and quantifying accelerations of slope movement

Time-lapse 4-D resistivity imaging inversion with positivity constraints

© 2019 24th European Meeting of Environmental and Engineering Geophysics. All rights reserved. Time-lapse resistivity surveys are used to monitor changes in the subsurface. In some situations, it is known that the resistivity will only decrease (or increase) with time. The 4-D ERT smoothness-constrained inversion method, that includes temporal smoothness constraint, has proved to be a robust method that reduces artefacts due to noise. However, in some cases, the time-lapse inverse models might show an increase in the resistivity with time where it is only expected to decrease. We modify the 4-D ERT inverse method to include a constraint that removes this artefact. The standard 4-D ERT inversion algorithm is first used to generate an initial model. If the resistivity is expected to decrease with time, for the model cells that show a resistivity increase with time, a truncation procedure is used where the resistivities of the different time models are reset to the mean value (corresponding to zero change with time). We then use the method of transformations in the inversion method that ensures the resistivities of the later time models are always less than the first model. The constraints can be modified so that they are only applied to selected regions in the model in cases where additional information is available