Ground penetrating radar (GPR) and electrical resistivity tomography (ERT) survey on the neotectonic activity of the Padul-Nigüelas Fault Zone (PNFZ), Southern Spain

The Padul-Nigüelas Fault Zone (PNFZ) is situated at the south-western mountain front of the Sierra Nevada (Spain) in an extensive regime and belongs to the internal zone of the Betic Cordilleras. The aim of this study is a collection of new evidence for neotectonic activity of the fault zone with classical geological field work and modern geophysical methods, such as ground penetrating radar (GPR). Among an apparently existing bed rock fault scarp with triangular facets, other evidences, such as deeply incised valleys and faults in the colluvial wedges, are present in the PNFZ. The preliminary results of our recent field work have shown that the synsedimentary faults within the colluvial sediments seem to propagate basinwards and the bed rock fault is only exhumed due to erosion for the studied segment (west of Marchena). We will use further GPR data and geomorphologic indices to gather further evidences of neotectonic activity of the PNFZ.

Soil moisture variability in a temperate deciduous forest: insights from electrical resistivity and throughfall data

Hydrogeophysics is a growing discipline that holds significant promise to help elucidate details of dynamic processes in the near surface, built on the ability of geophysical methods to measure properties from which hydrological and geochemical variables can be derived. For example, bulk electrical conductivity is governed by, amongst others, interstitial water content, fluid salinity, and temperature, and can be measured using a range of geophysical methods. In many cases, electrical resistivity tomography (ERT) is well suited to characterize these properties in multiple dimensions and to monitor dynamic processes, such as water infiltration and solute transport. In recent years, ERT has been used increasingly for ecosystem research in a wide range of settings; in particular to characterize vegetation-driven changes in root-zone and near-surface water dynamics. This increased popularity is due to operational factors (e.g., improved equipment, low site impact), data considerations (e.g., excellent repeatability), and the fact that ERT operates at scales significantly larger than traditional point sensors. Current limitations to a more widespread use of the approach include the high equipment costs, and the need for site-specific petrophysical relationships between properties of interest. In this presentation we will discuss recent equipment advances and theoretical and methodological aspects involved in the accurate estimation of soil moisture from ERT results. Examples will be presented from two studies in a temperate climate (Michigan, USA) and one from a humid tropical location (Tapajos, Brazil).

Relationship between electrical resistivity and basic geotechnical parameters for marine clays

Recently, considerable efforts have been made in the attempt to map quick clay areas using electrical resistivity measurements. However there is a lack of understanding regarding which soil parameters control the measured resistivity values. To address this issue, inverted resistivity values from 15 marine clay sites in Norway have been compared with basic geotechnical index properties. It was found that the resistivity value is strongly controlled by the salt content of the pore fluid. Resistivity decreases rapidly with increasing salt content. There is also a relatively clear trend of decreasing resistivity with increasing clay content and plasticity index. Resistivity values become very low (˜5 O·m) for high clay content (>50%), medium- to high-plasticity (Ip ˜ 20%) materials with salt content values greater than about 8 g/L (or corresponding remoulded shear strength values greater than 4 kPa). For the range of values studied, there is poor correlation between resistivity and bulk density and between resistivity and water content. The data studied suggest that the range of resistivity values corresponding to quick clay is 10 to 100 O·m, which is consistent with other published limits. A comparison is made between two-dimensional electrical resistivity tomography (ERT) and resistivity cone penetration test (RCPTU) data for two of the sites and the two sets of data show similar trends and values irrespective of scale effect.

Integrated Earth Resistivity Tomography (ERT) and Multilevel Gas Sampling: A Tool to Map Geogenic and Anthropogenic Methane Accumulation on Brownfield Sites

Soil gas emissions of methane and carbon dioxide on brownfield sites are usually attributed to anthropogenic activities; however geogenic sources of soil gas are often not considered during site investigation and risk management strategies. This paper presents a field study at a redeveloped brownfield site on a flood plain to identify accumulations of methane biogas trapped in underlying sediments. The investigation is based on a multidisciplinary approach using direct multi-level sampling measurements and Earth resistivity tomography . Resistivity imaging was applied to evaluate the feasibility of identifying the size and spatial continuity of soil gas accumulations in anthropogenic and naturally occurring deposits. As a result, biogas accumulations are described within both anthropogenic deposits and pristine organic sediments. This result is important to identify the correct approaches to identify and manage risks associated with soil gas emissions on brownfield and pristine sites. The organic-rich sediments in Quaternary fluvial environments of São Paulo Basin in particular the Tietê River, biogas reservoirs can be generated and trapped beneath geogenic and anthropogenic layers, potentially requiring the management of brownfield developments across this region.

Comparative analysis between biogas flow in landfill and electrical resistivity tomography in Rio Claro city, Brazil

The biogas originated from anaerobic degradation of organic matter in landfills consists basically in CH4, CO2, and H2O. The landfills represent an important depository of organic matter with high energetic potential in Brazil, although with inexpressive use in the present. The estimation of production of the productive rate of biogas represents one of the major difficulties of technical order to the planning of capture system for rational consumption of this resource. The applied geophysics consists in a set of methods and techniques with wide use in environmental and hydrogeological studies. The DC resistivity method is largely applied in environmental diagnosis of the contamination in soil and groundwater, due to the contrast of electrical properties frequent between contaminated areas and the natural environment. This paper aims to evaluate eventual relationships between biogas flows quantified in drains located in the landfill, with characteristic patterns of electrical resistivity in depth. The drain of higher flow (117 m3 /h) in depth was characterized for values between 8000 Ω⋅m and 100.000 Ω⋅m, in contrast with values below 2000 Ω⋅m, which characterize in subsurface the drain with less flow (37 m3 /h), besides intermediary flow and electrical resistivity values, attributed to the predominance of areas with accumulation or generation of biogas.

Geotechnical monitoring using a combination of LIDAR, real-time dGPS, and electrical tomography to assess geotechnical risk of a coastal landslide.

The Antrim Coast Road stretching from the seaport of Larne in the East of Northern Ireland has a well-deserved reputation for being one of the most spectacular roads in Europe (Day, 2006). However the problematic geology; Jurassic Lias Clay and Triassic Mudstone overlain by Cretaceous Limestone and Tertiary Basalt, and environmental variables result in frequent instances of slope instability manifested in both shallow debris flows and occasional massive rotational movements, creating a geotechnical risk to this highway. This paper describes how a variety of techniques are being used to both assess instability and monitor movement of these active slopes near one site at Straidkilly Point, Glenarm. An in-depth understanding of the geology was obtained via boreholes, resistivity surveys and laboratory testing. Environmental variables recorded by an on-site weather station were correlated with measured pore water pressure and soil moisture infiltration data. Terrestrial LiDAR (TLS), with surveys carried out on a bi-monthly basis allowed for the generation of Digital Elevation Models (DEMs) of difference, highlighting areas of recent movement, accumulation and depletion. Morphology parameters were generated from the DEMs and include slope, curvature and multiple measures of roughness. Changes in the structure of the slope coupled with morphological parameters were characterised and linked to progressive failures from the temporal monitoring. In addition to TLS monitoring, Aerial LiDAR datasets were used for the spatio-morphological characterisation of the slope on a macro scale. A Differential Global Positioning System (dGPS) was also deployed on site to provide a real-time warning system for gross movements, which were also correlated with environmental conditions. Frequent electrical resistivity tomography (ERT) surveys were also implemented to provide a better understanding of long-term changes in soil moisture and help to define the complex geology. The paper describes how the data obtained via a diverse range of methods has been combined to facilitate a more informed management regime of geotechnical risk by the Northern Ireland Roads Service.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

Appraisal of groundwater resource in Holocene soil deposits by resistivity, hydrochemical and granulomerial studies in the Gulf of Mannar Coast from Southern India

This article is aimed to delineate groundwater sources in Holocene deposits area in the Gulf of Mannar Coast from Southern India. For this purpose 2-D electrical resistivity tomography (ERT), hydrochemical and granulomerical studies were carried out and integrated to identify hydrogeological structures and portable groundwater resource in shallow depths which in general appears in the coastal tracts. The 2-D ERT was used to determine the two-dimensional subsurface geological formations by multicore cable with Wenner array. Low resistivity of 1-5 Omega m for saline water appeared due to calcite at the depth of about 5 m below the ground level (bgl). Sea water intrusion was observed around the maximum resistivity as 5 Omega m at the 8 m depth, bgl in the calcite environs, but the calcareous sandstone layer shows around 15-64 Omega m at the 6 m depth, bgl. The hydrochemical variation of TDS, HCO3-, Cl-, Na+, K+, Ca2+, and Mg2+ concentrations was observed for the saline and sea water intrusion in the groundwater system. The granulometic analysis shows that the study area was under the sea between 5400 and 3000 year ago. The events of ice melting an unnatural ice-stone rain/hail among 5000-4000 years ago resulted in the inundation of sea over the area and deposits of late Holocene marine transgression formation up to Puthukottai quartzite region for a stretch of around 17 km.

Application of Electrical Resistance Tomography System to Monitor Gas/Liquid Two-Phase Flow in a Horizontal Pipe

The Electrical Resistance Tomography (ERT) technique possesses great potential in monitoring widely exiting industrial two/multi-phase flow. For vertical pipe flow and inclined pipe flow, some application studies with exciting results have been reported, but there is rarely a paper regarding the application of ERT to horizontal gas/liquid pipe flow. This paper addresses this issue and proposes a smart method, Liquid Level Detection method, to conventional ERT system. The enhanced ERT system using the new method can monitor horizontal pipe flow effectively and its application is no longer restricted by the flow conditions. Some experimental results from monitoring an air/water slug pipe flow are presented.

MEASUREMENT OF SOLID SLURRY FLOW VIA CORRELATION OF ELECTROMAGNETIC FLOW METER, ELECTRICAL RESISTANCE TOMOGRAPHY AND MECHANISTIC MODELLING

The study presented here was carried out to obtain the actual solids flow rate by the combination of electrical resistance tomography and electromagnetic flow meter. A new in-situ measurement method based on measurements of the Electromagnetic Flow Meters (EFM) and Electrical Resistance Tomography (ERT) to study the flow rates of individual phases in a vertical flow was proposed. The study was based on laboratory experiments that were carried out with a 50 mm vertical flow rig for a number of sand concentrations and different mixture velocities. A range of sand slurries with median particle size from 212 mu m to 355 mu m was tested. The solid concentration by volume covered was 5% and 15%, and the corresponding density of 5% was 1078 kg/m(3) and of 15% was 1238 kg/m(3). The flow velocity was between 1.5 m/s and 3.0 m/s. A total of 6 experimental tests were conducted. The equivalent liquid model was adopted to validate in-situ volumetric solids fraction and calculate the slip velocity. The results show that the ERT technique can be used in conjunction with an electromagnetic flow meter as a way of measurement of slurry flow rate in a vertical pipe flow. However it should be emphasized that the EFM results must be treated with reservation when the flow pattern at the EFM mounting position is a non-homogenous flow. The flow rate obtained by the EFM should be corrected considering the slip velocity and the flow pattern.

Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers

The validation of variable-density flow models simulating seawater intrusion in coastal aquifers requires information about concentration distribution in groundwater. Electrical resistivity tomography (ERT) provides relevant data for this purpose. However, inverse modeling is not accurate because of the non-uniqueness of solutions. Such difficulties in evaluating seawater intrusion can be overcome by coupling geophysical data and groundwater modeling. First, the resistivity distribution obtained by inverse geo-electrical modeling is established. Second, a 3-D variable-density flow hydrogeological model is developed. Third, using Archie's Law, the electrical resistivity model deduced from salt concentration is compared to the formerly interpreted electrical model. Finally, aside from that usual comparison-validation, the theoretical geophysical response of concentrations simulated with the groundwater model can be compared to field-measured resistivity data. This constitutes a cross-validation of both the inverse geo-electrical model and the groundwater model.
[Comte, J.-C., and O. Banton (2007), Cross-validation of geo-electrical and hydrogeological models to evaluate seawater intrusion in coastal aquifers, Geophys. Res. Lett., 34, L10402, doi:10.1029/2007GL029981.]

Specific depth cone resistivity measurements to determine soil engineering properties

The University of British Columbia (UBC) began performing piezocone penetration tests (CPTU) with electrical resistivity measurements (RCPTU) in 1989. Since then, RCPTU research at UBC has focused on obtaining geo-environmental parameters such as fluid resistivity and soil engineering properties such as porosity and degree of saturation from measurements of bulk soil electrical resistivity using the empirical relationship proposed by Archie (1942). Within this framework, the paper illustrates and discusses important design and calibration issues for resistivity modules such as the use of isolated circuitry to achieve linear calibrations over large ranges of resistivity. The suitability of RCPTU measurements for determination of geo-environmental and geotechnical parameters are assessed using typical ranges of soil and groundwater properties and methods of isolating individual factors for study are discussed. Illustrative examples of RCPTU research efforts including the environmental characterization of mine tailings, delineation of saline water intrusions in fresh water aquifers and the quality control of geotechnical ground densification are presented throughout the text. It is shown that groundwater temperature and hence ion mobility is not significantly altered by frictional heat generated during piezocone penetration and that ratio-based approaches to monitoring soil porosity can be used to eliminate the requirement for extensive groundwater sampling programs. Lastly, it is shown that RCPTU measurements above the water table can only be made using resistivity modules that are stable over a large range of resistivities and that such measurements are the most difficult to interpret because of grain surface conduction effects and generally unknown fluid resistivities.

The effects of deep cracks on the rain-induced instability of slopes : A case study

Rainfall can disrupt the balance of natural soil slope. This imbalance will be accelerated by existence of cracks in soil slope, which lead to decreasing shear strength and increasing hydraulic conductivity of the soil slope. Some research works have been conducted on the effects of surface-cracks on slope stability. However, the influence of deep-cracks is yet to be investigated. Limited availability of deep crack data due to the lack of effective sub-soil investigation methods could be one of the obstacles. To emphasize the effects of deep cracks in soil slope on its rain-induced instability, a natural soil slope in Indonesia that failed in 31st October 2010 due to heavy rainfall was analyzed for stability with and without deep cracks in the slope. The slope stability analysis was conducted using SLOPE/W coupling with the results of transient seepage analysis (SEEP/W) that simulate the pore-water pressure development in the slope during the rainfall. The results of Electrical Resistivity Tomography (ERT) survey, bore-hole tests and geometrical survey conducted on the slope before its failure were used to identify the soil layers’ stratification including deep cracks, the properties of different soil layers, and geometrical parameters of the slope for the analysis. The results showed that it is vital to consider the existence of deep crack in soil slopes in analysing their instability induced by rainfalls.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.