Time-lapse Electrical Resistivity Tomography Research Articles

Electrical resistivity tomography: A reliable tool to monitor the efficiency of different irrigation systems in horticulture field

Water management in agricultural systems is essential for optimal crop yields without incurring excessive water costs and wastage. The choice of irrigation method is crucial for better water management and distribution. The drip system appears to be among the best methods in the field of precision agriculture. In addition to the irrigation system, mulching with ridge plastic film to drain excess water is widely used to increase crop yields in terms of plant water availability. In this study, the time-lapse Electrical Resistivity Tomography (ERT), a not-invasive geophysical technique, is proposed as a simple and reliable method to evaluate the effectiveness of the irrigation systems and to monitor the changes in water content over time and over a volume of soil. ERTs data were compared to moisture ones retrieved from sensors that record continuously over time, but punctually. The ERT investigations were conducted in melon-growing lands in southern Tuscany (Italy). Measurements were carried out on two different fields in two periods: spring and summer. The aim of the work was to evaluate, by means of volumetric measures of the soil conductivity, the effectiveness of three different drip systems and of the mulch ridge. In both the monitored fields the ridge was created in a half portion of the field itself, while the other part of the land was left plat. Geoelectrical investigations associated with humidity sensors have shown that in the summer a too high mulch ridge quickly drains the irrigation water, bringing the root zone into a water deficit. The ERTs also provided good results relating to the irrigation system, demonstrating that a three-lines drip irrigation system, compared to a two-lines one, manages to distribute the irrigation water homogeneously, guaranteeing a constant water content for the plants over time.

Employing automated electrical resistivity tomography for detecting short- and long-term changes in permafrost and active-layer dynamics in the maritime Antarctic

Abstract. Repeated electrical resistivity tomography (ERT) surveys can substantially advance the understanding of spatial and temporal freeze–thaw dynamics in remote regions, such as Antarctica, where the evolution of permafrost has been poorly investigated. To enable time-lapse ERT surveys in Antarctica, an automated ERT (A-ERT) system is required, as regular site visits are not feasible. In this context, we developed a robust A-ERT prototype and installed it at the Crater Lake CALM-S site on Deception Island, Antarctica, to collect quasi-continuous ERT measurements. We developed an automated data processing workflow to efficiently filter and invert the A-ERT datasets and extract the key information required for a detailed investigation of permafrost and active-layer dynamics. In this paper, we report on the results of two complete year-round A-ERT datasets collected in 2010 and 2019 at the Crater Lake CALM-S site and compare them with available climate and borehole data. The A-ERT profile has a length of 9.5 m with an electrode spacing of 0.5 m, enabling a maximum investigation depth of approximately 2 m. Our detailed investigation of the A-ERT data and inverted results shows that the A-ERT system can detect the active-layer freezing and thawing events with high temporal resolution. The resistivity of the permafrost zone in 2019 is very similar to the values found in 2010, suggesting the stability of the permafrost over almost 1 decade at this site. The evolution of thaw depth exhibits a similar pattern in both years, with the active-layer thickness fluctuating between 0.20–0.35 m. However, a slight thinning of the active layer is evident in early 2019, compared to the equivalent period in 2010. These findings show that A-ERT datasets, combined with the new processing workflow that we developed, are an effective tool for studying permafrost and active-layer dynamics with very high resolution and minimal environmental disturbance. The ability of the A-ERT setup to monitor the spatiotemporal progression of thaw depth in two dimensions, and potentially in three dimensions, and to detect brief surficial refreezing and thawing of the active layer reveals the significance of the automatic ERT monitoring system to record continuous resistivity changes. An A-ERT monitoring setup with a longer profile length can investigate greater depths, offering effective monitoring at sites where boreholes are costly and invasive techniques are unsuitable. This shows that the A-ERT setup described in this paper can be a significant addition to the Global Terrestrial Network for Permafrost (GTN-P) and the Circumpolar Active Layer Monitoring (CALM) networks to further investigate the impact of fast-changing climate and extreme meteorological events on the upper soil horizons and to work towards establishing an early warning system for the consequences of climate change.

Mapping and monitoring dense non-aqueous phase liquid source zone by fused surface and cross-borehole electrical resistivity tomography

Effective characterization of dense non-aqueous phase liquid (DNAPL) source zones is crucial for remediating polluted sites. DNAPL often reside as residuals or pools within high-permeability lenses and above impermeable layers due to soil heterogeneity, gravity, and capillary barriers. Given the high cost of drilling, electrical resistivity tomography (ERT) techniques—including surface ERT and cross-borehole ERT, are commonly used for DNAPL source zone mapping and monitoring. However, the low spatial resolution of ERT increases uncertainty in source zone investigations. This study proposes a method for improving DNAPL mapping and monitoring by fusing surface and cross-borehole ERT data. Sandbox experiments were conducted to simulate a heterogeneous DNAPL source zone, employing both ERT methods for static mapping and dynamic monitoring. Reflective light imaging (RLM) was used to visualize DNAPL migration and provide saturation data, allowing for the quantification of ERT's effectiveness in characterizing DNAPL distribution. The results indicate that individual ERT methods face significant challenges in DNAPL source zone mapping due to background interference. Surface ERT alone tends to underestimate the extent of deeper DNAPL source zones. However, fusing surface and cross-borehole ERT results in a complementary enhancement of vertical spatial resolution, thereby improving the characterization of DNAPL source zones. The fusion of static and time-lapse ERT data substantially enhances DNAPL source zone mapping and monitoring capabilities. By calculating the ratio of the ERT-monitored area to the actual area using resistivity change contours (5 %, 10 %, 15 %), it was found that fusing surface and cross-borehole ERT data improved monitoring resolution by 50.48 % compared to surface ERT alone and by 22.95 % compared to cross-borehole ERT. Principal component analysis (PCA) was effective in fusing time-lapse data, while the weighted average method (WAM) outperformed PCA for static resistivity data fusion.

Organic soil amendment effects on soil hydrology in an almond orchard evaluated using time-lapse electrical resistivity tomography

Adding organic amendments to soils in orchards has been suggested as a climate-smart agricultural practice that can increase resilience to extremes such as drought. The benefits of adding almond hulls and shells to soil include releasing potassium into the soil, improving water infiltration, reducing soil water evaporation, and enhancing water-holding capacity. In this study, soil infiltration and root water uptake (RWU) patterns of almond trees in amended and control soil treatments were investigated. An almond orchard was mulched with a mixture of almond hulls and shells used as surface-applied organic matter amendments. The combined use of time-lapse electrical resistivity tomography (ERT), stem water potential (SWP) and leaf water potential (LWP) measurements were used to evaluate RWU, soil infiltration patterns, and tree water status at different times during the study. The results of the ERT showed that the RWU patterns for the almond trees are distinct between the amended and the control treatments after the irrigation was applied. Compared to the control treatment, the amended treatment allowed the soil to store more water. Regardless of treatment, the majority of RWU patterns were observed in the top 0.5–1 m root zone depth. Almond trees began to recover from stress, as indicated by the SWP and LWP values, four hours after the start of irrigation. In addition, this was demonstrated by the ERT measurements, which revealed the RWU activity four hours after the irrigation had been applied. This research reveals that time-lapse ERT surveys combined with soil and tree water status data can infer patterns of RWU in almond trees grown with organic soil amendments and under controlled conditions. Furthermore, it was concluded that using almond hulls and shells as organic matter amendments can help almond growers improve infiltration, making almond production more resilient to climate change-related extremes such as droughts.

The contraction of freshwater lenses in barrier island: A combined geophysical and numerical analysis

Coastal aquifers, vital for supplying freshwater to over one billion people worldwide, often face saltwater intrusion, with barrier island aquifers playing a crucial role in this coastal system. Despite the proliferation of studies exploring the impacts of sea-level rise, storm surges, and over-pumping, the effect of droughts on barrier island aquifers salinization remains largely under-investigated. This lack of attention is particularly concerning given the heightened vulnerability of barrier island aquifers; most ultimately rely solely on aerial recharge compared to continental coastal aquifers. An in-depth understanding of recharge and salinization processes is important for sustainably managing these islands' marginal freshwater resources, especially considering climate change impacts. This study presents an evaluation approach for the response of a freshwater lens (FWL) to drought conditions, incorporating in-situ observations, geophysical measurements, and numerical modeling. The study examines the response of a shallow unconfined aquifer on an Atlantic coast barrier island to the 2020 Northeast extreme hydrological drought. It encompasses a density-driven flow model informed by time-lapse electrical resistivity imaging and in-situ measurements, including groundwater head and salinity collected from monitoring wells. Results from our approach indicate that FWL volume is reduced by 11% during the drought, but it returns to its previous volume over the following spring season, indicating that the net volume of the FWL remains unchanged on an annual scale. The mean groundwater residence time on the island is approximately a year, indicating that the barrier island aquifer is a hydrodynamically active coastal zone. These findings highlight the vulnerability and resilience of shallow unconfined barrier island aquifers to droughts and climate change.

Monitoring moisture dynamics in multi-layer cover systems for mine tailings reclamation using autonomous and remote time-lapse electrical resistivity tomography

The dynamics of moisture content in cover systems constructed on mining wastes were monitored at the pilot scale using 2D autonomous, remote, and non-invasive time-lapse electrical resistivity tomography combined with conventional point sensors. A methodology was proposed to process the daily hydrogeophysical datasets from 23 m-long instrumented sections of covers with capillary barrier effects (CCBEs) designed to act as oxygen barriers, and covers with low saturated hydraulic conductivity layers (LSHCCs) designed to limit the water infiltration rate. Hydrogeophysical monitoring suggested that CCBEs were able to maintain high degrees of saturation in the moisture-retaining layer throughout the one-year monitoring period, which would make it an efficient oxygen barrier. Larger spatio-temporal changes in moisture content were observed in LSHCCs and most of the low hydraulic conductivity layers remained below 85% saturation, which was attributed to the combined effect of low precipitation, rapid vegetation development, and water percolation through the cover. The methodology proposed in this pilot-scale "proof-of-concept" study allowed the hydrogeological behavior of the cover systems to be monitored in the 23 m-long instrumented sections using continuous geoelectrical data, which demonstrated that this innovative monitoring technique could be useful for geochemical and geotechnical monitoring programs in large-scale mining waste storage facilities.

Non-invasive geophysical methods for monitoring the shallow aquifer based on time-lapse electrical resistivity tomography, magnetic resonance sounding, and spontaneous potential methods

The Ningdong coalfield has played a pivotal role in advancing local economic development and meeting national energy. Nevertheless, mining operations have engendered ecological challenges encompassing subterranean water depletion, land desertification, and ground subsidence, primarily stemming from the disruption of coal seam roof strata. Consequently, the local ecosystem has incurred substantial harm. Water-preserved coal mining presently constitutes the pivotal technology in mitigating this problem. The primary challenge of this technique lies in identifying critical aquifer layers and understanding the heights of water-conducting fracture zones. To obtain a precise comprehension of the seepage patterns within the upper coal seam aquifer during mining, delineate the extent of water-conducting fracture zones, non-invasive geophysical techniques such as time-lapse electrical resistivity tomography (TL-ERT), magnetic resonance sounding (MRS), and spontaneous potential (SP) have been employed to monitor alterations within the shallow coalfield’s aquifer throughout the mining process in the Ningdong coalfield. By conducting meticulous examinations of fluctuations in resistivity, moisture content, and self-potential within the superjacent strata during coal seam extraction, the predominant underground water infiltration strata were ascertained, concurrently enabling the estimation of the development elevation of water-conducting fracture zones. This outcome furnishes a geophysical underpinning for endeavors concerning local water-preserved coal mining and ecological rehabilitation.

Assessing the efficiency of the irrigation system in a horticulture field through time-lapse electrical resistivity tomography

AbstractThe characterization by means of geophysical techniques of agricultural soils subjected to continuous irrigation cycles makes it possible to study the heterogeneity of a soil and the preferential pathways of water flows without disturbing soil and plants. A better knowledge of soil heterogeneity enables optimal water resource management in terms of crop, yield, and sustainability. In this study, time-lapse monitoring using electrical resistivity tomographies (ERT) is proposed as a reliable and non-invasive technique to quantify the movement of water flows and thus the variation of soil water content during the irrigation process. ERT surveys have been conducted in melon-growing soils in southern Tuscany (Italy). Five survey campaigns have been carried out between June and August 2022, in which ERT data have been collected by taking measurements before (T0), during (T1), and after (T2) the irrigation phase. The interpretation of the ERT results provided information on the spatial and temporal distribution of water fluxes in the soil and root zone of melons during the irrigation phases. The investigation made it possible to identify the preferential pathways of infiltration of irrigation water, the points where water is absorbed by the roots, and the points where water follows a preferential pathway instead distributing itself entirely below the root growth zone. Thus, this research suggests that the ERT technique can be used to evaluate the efficiency of the irrigation system in order to achieve optimal management of the water resource, avoiding preferential flow paths that lead to less water availability for the plant.

Time-Lapse Geophysical Measurements for Monitoring Coastal Groundwater Dynamics in an Unconfined Aquifer.

The coastal zone, which is the interface between land and sea, is hydrodynamically very active due to the complex interactions of various hydrological controls and variable-density fluids. These forces vary over time, resulting in a state of dynamic equilibrium in the system. The major hydrological processes in coastal aquifer systems are salt water intrusion and submarine groundwater discharge, which are interdependent. Monitoring these complex processes is crucial for sustainable coastal zone management but poses a significant research challenge. In this study, we demonstrate the effectiveness of non-invasive geophysical techniques, specifically the time-lapse electrical resistivity imaging method, in conjunction with groundwater monitoring, for monitoring coastal groundwater dynamics in an unconfined aquifer at varying time scales and hydrogeological settings present at formerly glaciated sites worldwide. We generated two-dimensional baseline salt water intrusion maps for the test site, located on the coast of Rhode Island, USA. The time-lapse electrical resistivity survey method enables the rapid estimation of fresh groundwater discharge. Our approach offers insight into the mechanisms and seasonably variable salt water-freshwater interactions in unconfined heterogeneous aquifers. Although the results are site-specific, their implications are broad and may stimulate other studies related to sea to land pollution (sea water intrusion) and land to sea pollution (groundwater discharge) in heterogeneous coastal aquifer settings.

Detecting soil water redistribution in subsurface drip irrigated processing tomatoes using electrical resistivity tomography, proximal sensing and hydrological modelling

In this study, multiple soil-plant-atmosphere continuum (SPAC) monitoring methodologies, including electrical resistivity tomography (ERT), proximal thermal sensing techniques, and micrometeorological data, were combined with two-dimensional (2-D) soil hydrological modelling using HYDRUS 2-D to explore the soil water redistribution, and infer the relative crop water status in a subsurface drip irrigated (SDI) processing tomato field located in California (Yolo County, USA). Specifically, time-lapse ERT surveys were performed at two transects distributed parallel and perpendicular, respectively, to the SDI line, during an irrigation event. The ERT results were compared to HYDRUS 2-D outputs and the relative differences were explained in the form of local heterogeneities in electrical resistivity (ER) changes, as a proxy for soil water content (SWC) variations. Concurrent simultaneous soil wetting and root water uptake during the last irrigation event of the season caused negligible changes in ER in the active root zone. Slight differences in ER were observed in the top 20 cm along the dripline, confirming that the emitter spacing is small enough to create a wetted strip along the processing tomato bed. These changes were also compared to SWC values measured with time domain reflectometry soil moisture sensors. A comparison between HYDRUS 2-D and ERT confirmed negligible changes in ER during irrigation due to simultaneous wetting and root water uptake processes. In addition, a good correlation was observed between the proximal sensed and the ERT results. Finally, the findings of this study underscore the necessity of using multiple methods for improving our knowledge of the SPAC system under real field conditions.

A multiscale accuracy assessment of moisture content predictions using time-lapse electrical resistivity tomography in mine tailings

Accurate and large-scale assessment of volumetric water content (VWC) plays a critical role in mining waste monitoring to mitigate potential geotechnical and environmental risks. In recent years, time-lapse electrical resistivity tomography (TL-ERT) has emerged as a promising monitoring approach that can be used in combination with traditional invasive and point-measurements techniques to estimate VWC in mine tailings. Moreover, the bulk electrical conductivity (EC) imaged using TL-ERT can be converted into VWC in the field using petrophysical relationships calibrated in the laboratory. This study is the first to assess the scale effect on the accuracy of ERT-predicted VWC in tailings. Simultaneous and co-located monitoring of bulk EC and VWC are carried out in tailings at five different scales, in the laboratory and in the field. The hydrogeophysical datasets are used to calibrate a petrophysical model used to predict VWC from TL-ERT data. Overall, the accuracy of ERT-predicted VWC is pm 0.03~textrm{m}^3/textrm{m}^3, and the petrophysical models determined at sample-scale in the laboratory remain valid at larger scales. Notably, the impact of temperature and pore water EC evolution plays a major role in VWC predictions at the field scale (tenfold reduction of accuracy) and, therefore, must be properly taken into account during the TL-ERT data processing using complementary hydrogeological sensors. Based on these results, we suggest that future studies using TL-ERT to predict VWC in mine tailings could use sample-scale laboratory apparatus similar to the electrical resistivity Tempe cell presented here to calibrate petrophysical models and carefully upscale them to field applications.

LNAPL migration processes based on time-lapse electrical resistivity tomography

Contamination from light non-aqueous phase liquids (LNAPLs) and their derivatives, arising from exploration, production, and transportation, has become a prevalent pollution source. This poses direct threats to human health. However, conventional investigative methods face limitations when applied to studying the extent and migration process of LNAPL contamination, as well as the redistribution of LNAPL during groundwater level fluctuations. Conventional methods lack the ability to rapidly, efficiently, and in real-time acquire information about contaminated areas. Therefore, this study utilizes time-lapse electrical resistivity tomography to investigate the migration mechanism of LNAPL under unsaturated conditions, constant groundwater levels, and groundwater level reductions. A relationship between resistivity and water and oil contents was established and used for inverse calculation of LNAPL content via resistivity inversion. Time-lapse electrical resistivity tomography revealed LNAPL migration in a “concave” shape across three conditions. Groundwater presence notably slowed migration, hindering downward movement and leading to a floating oil band. A robust mathematical model was established to derive the relationship between resistivity and water and oil contents. Finally, LNAPL distribution under unsaturated conditions was inversely obtained from resistivity data, showing highest content at the top leak point, obstructed area, and bottom of soil column. Consequently, time-lapse electrical resistivity tomography demonstrates a notable capacity to characterize the LNAPL migration process. This technique constitutes an effective geophysical method for monitoring and describing the characteristics of LNAPL migration. Its significance lies in enhancing our understanding of remediation for LNAPL-induced groundwater and land contamination.

Data-model comparisons of isotopic and hydrological variabilities of the karstic vadose zone above Villars Cave, SW-France based on 20 years' monitoring record

Because of the great heterogeneity of the limestone fissures, conduits, and reservoirs in karst areas, the understanding of the infiltration in caves and karst terrains is complicated. For example, the time delays observed between rainfall and the various flows into caves, such as shaft flow and seepage flow, exhibit significant variability. To better understand the infiltration above a typical shallow cave, Villars Cave, South-west France, we used both a very long monitoring (>20 years) isotopic data set of rainfall and dripping rates at two different cave gallery levels (14 m and 36 m below the surface, respectively) and a conceptual model (KarstFor) that mix water flows and isotopes. In order to determine the most appropriate reservoir and flow settings, we used one time-lapse electrical resistivity tomography (ERT) image made above the cave, which were incorporated into the KarstFor model. The modeled variations in storage water and infiltration flows are in a good agreement with the observed dripping discharge rates, indicating seasonal variability and various time delays between water excess ((Rainfall-Evapotranspiration) > 0) and drip rates in the different cave levels. The fact that the model strikingly simulates the small, but constant, isotopic difference of drip water δ18O (δ18Od) values between the upper and the lower galleries of Villars Cave, not only reinforces reliability in the model itself but also reveals the impact of multiple infiltration routes and karst reservoir dynamics on the δ18Od. On an inter-annual scale, the observed δ18Od values follow the variations of the water excess and of the cave temperature, which, in the period analyzed, show a marked trend, possibly induced by the climate variability patterns above the cave. Overall, these results are not only important to better understand the infiltration in karstic zones, but also to better interpret speleothem isotopic records widely used in paleo-climate studies.

Imaging hydrological dynamics in karst unsaturated zones by time-lapse electrical resistivity tomography

The hydrodynamics of karst terrain are highly complex due to the diverse fractures and reservoirs within limestone formations. The time delay between rainfall events and subsequent flow into reservoirs exhibits significant variability. However, these hydrological processes are not easily visualized in karst topography. Subsurface geophysics, specifically 2D time-lapse electrical resistivity tomography (ERT), provides an effective method for studying the relationships between hydrological and geophysical features. In our research, we adopted ERT in the Karst Critical Zone (KCZ) to visualize specific karstic zones, including cave galleries, water storage reservoirs, wetting fronts, soil layers, and potential preferential flow paths down to a depth of 20 m. To capture spatial and seasonal variations in resistivity, we presented a comprehensive approach by combining sixteen inversion models obtained between February 2020 and September 2022 above the Villars Cave in SW-France—a well-known prehistoric cave. We used a multi-dimensional statistical technique called Hierarchical Agglomerative Clustering (HAC) to create a composite model that divided the synthetic ERT image into eight clusters representing different karst critical zones. The ERT image clearly visualized the cave gallery with high resistivity values that remained consistent throughout the seasons. Our analysis revealed a close seasonal relationship between water excess and resistivity variations in most infiltration zones, with time delays increasing with depth. The karst reservoirs, located at significant depths compared to other clusters, displayed sensitivity to changes in water excess but were primarily affected by fluctuations in water conductivity, particularly during summer or dry periods. These findings have significant implications for predicting rainwater infiltration pathways into caves, thereby assisting in the conservation and preservation of prehistoric caves and their cultural heritage.

Contaminant source and aquifer characterization: An application of ES-MDA demonstrating the assimilation of geophysical data

Contaminant source and aquifer characterization (CSAC) is critical in groundwater pollution evaluation and remediation. The ensemble smoother with multiple data assimilation (ES-MDA) is utilized to jointly identify contaminant source information and hydraulic conductivities by assimilating time-lapse electrical resistivity tomography (ERT) data. In a synthetic profile with a time-varying release history in a heterogeneous aquifer, we verify the performance of the proposed data assimilation framework. The results show that the CSAC problem could be handled by the proposed approach. The time-varying release history and the high permeability area can be identified with adequate time-lapse ERT measurements. Further reproduction of the evolution of the plume after CSAC also shows consistency with the reference plume. The poor conditioning inversion caused by the filter inbreeding is analyzed by comparing four scenarios with different apparent resistivity measurements. Furthermore, we also evaluate the impact of uncertainties in the petrophysical properties and geophysical observations on our data assimilation framework. The results show that the proposed ES-MDA data assimilation framework could provide a convincing inversion of the time-varying release history and hydraulic conductivities.

Assessing the risk of slope failure to highway infrastructure using automated time-lapse electrical resistivity tomography monitoring

Electrical resistivity tomography (ERT) monitoring provides time-lapse images of the subsurface. These images can be used to assess spatiotemporal variation in moisture content, which is a key driver of slope failure, making ERT monitoring an effective tool to evaluate precursory conditions of failure. This work presents the results of ERT monitoring on a slope above a major highway located on the border between England and Wales. During highway construction in the 1960s the slope was subject to several large landslide events which resulted in the re-design of the carriageway and installation of engineered mitigation measures. A section of the slope known as the ‘partially slipped area’ exhibited partial displacement during this time but did not progress to full slope failure, and therefore presents an ongoing risk to the highway, even though it does not experience ongoing displacement. An ERT monitoring system was installed across this area to monitor subsurface variations in moisture content. The results show a complex pattern of subsurface moisture dynamics within the partially slipped area when compared to the adjacent area of stable slope. This is most likely a result of the uneven and hummocky terrain in the partially slipped area and its effects on rainfall infiltration, storage and drainage, combined with the displacement-induced jointing present in the underlying sandstone units. The ERT results are used to assess the volume of unstable ground, placing the volume at the upper end of estimates from previous studies. Furthermore, analysis of the ERT dataset for surface displacements shows no movement at the site, which is confirmed by analysis of differential LiDAR plots and ground motion data derived from InSAR. This study demonstrates the application of ERT monitoring on a low activity, high risk slope, highlighting the need to understand subsurface processes at the slope-scale to inform long-term slope management.

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

Groundwater monitoring and specific yield estimation using time-lapse electrical resistivity imaging and machine learning

This paper presents an alternative method for monitoring groundwater levels and estimating specific yields of an unconfined aquifer under different seasonal conditions. The approach employs the Time-Lapse Electrical Resistivity Imaging (TL-ERI) method and machine learning-based time series clustering. A TL-ERI survey was conducted at ten sites (WS01-WS10 sites) throughout the dry and wet seasons, with five-time measurements collected for each site, in the Taichung-Nantou Basin along the Wu River, Central Taiwan. The obtained resistivity raw data was inverted and converted into normalized water content values using Archie’s law, followed by applying the Van Genuchten (VG) model for the Soil Water Characteristic Curve to estimate the Groundwater Level (GWL), and estimated the theoretical specific yield (Sy) by computing the difference between the saturated and residual water contents of the fitted VG model. In addition, the specific yield capacity (Sc), representing the nature of the storage capacity in the aquifer, was also calculated. The results showed that this approach was able to estimate those hydrogeological parameters. The spatial distribution of the GWL reveals that during the dry-wet seasons from February to July, there was a high GWL that extended from southeast to northwest. Conversely, during the wet-dry seasons from July to October, the high GWL shrank, which can be attributed to recharge variations from rainfall events. The determined spatial distribution of Sy and Sc fall within the range of 0.03–0.24 and 0.14–0.25, respectively. To quantitatively establish areas of similar groundwater level changes along with the VG model parameter variations during the study period, a Time series Clustering analysis (TSC) was performed by utilizing Hierarchical Agglomerative Clustering (HAC). The findings suggest that the WS03 site is a promising area for further investigation due to its highest Sc value with a slight change in groundwater levels during the dry and wet seasons. This study brings an advanced development of the geoelectrical method to estimate regional hydrogeological parameters in an area with limited available groundwater observation wells, in different seasonal conditions for groundwater management purposes.

Monitoring and characterization of water infiltration in soil unsaturated zone through an integrated geophysical approach

This paper presents an integrated geophysical method to characterize the vadose zone by conducting lab-scale infiltration experiments, during which time-lapse electrical resistivity tomography (ERT) sections and ground penetrating radar (GPR) profiles were simultaneously and continuously captured together with soil moisture and temperature probing devices installed at various depths. Time series water content data obtained from the moisture sensors, coupled with discreetly calibrated petrophysical relationships and time-lapse datasets from ERT and GPR, were combined to monitor the spatial–temporal evolution of soil water content during wetting and drying and to track the downward progression of wetting front. Reduction in electrical resistivity as well as increment in dielectric constant were converted to rising soil moisture, and then further scaled into cumulative infiltration through integration with corresponding wetting depth as observed through the experiment. Robust estimates regarding controlling vadose zone parameters such as infiltration rate and hydraulic conductivity were subsequently achieved by fitting the results from complementary investigations with an unsaturated infiltration model. With calibrated saturated hydraulic conductivity fitted well with reported and experimentally derived ranges, it’s suggested that the model is capable of reflecting average changes in infiltration rate at the lab scale. However, the model intrinsically lacks the ability to provide spatial variations of hydraulic parameters, thereby further investigation is essential to establish a suitable model for characterizing the more complex vadose zone at field scale.

4D interpretation of time-lapse electrical resistivity monitoring data to identify preferential flow path in a landfill, South Korea.

Monitoring the leakage of leachate from a landfill is critical in preventing possible contamination in the surrounding area. Time-lapse (TL) electrical resistivity tomography (ERT) has been performed along eleven survey lines at four different time points in a landfill in Korea. The TL data sets were interpreted using an in-house 4D inversion algorithm. Changes in 4D inversion results were analyzed in order to interpret a leachate-contaminant region. Since the rainy season started during obtaining TL ERT data sets, the effects of precipitation on TL ERT data are also analyzed. Changes in electrical resistivity (ER) showed that precipitation increases ER of contaminant zones. As hydrogeochemical data offer contamination information in some areas where boreholes are located, these are helpful to interpret and compare with ERT inversion results to evaluate the extent of the contaminated plume. We also classified soil textures from particle size analysis on soil samples and analyzed electrical conductivity (EC) and dissolved oxygen (DO) using groundwater samples obtained from observation wells in the survey site. The information on soil structure as well as the results of 4D inversion provided insight into the location of a preferential flow path.