Application of Recursive Estimation to Heat Tracing for Groundwater/Surface-Water Exchange.

We present and demonstrate a recursive-estimation framework to infer groundwater/surface-water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat-transport problem as a state-space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time-varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch-Tung-Striebel Smoother (ERTSS). Whereas the EKF is suited to real-time ("online") applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near-real time or batch-processing ("offline") applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field-experimental data and are shown to be efficient and rapid for the estimation of time-varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time-varying groundwater/surface-water exchange.

An overview of geophysical technologies appropriate for characterization and monitoring at fractured-rock sites.

Geophysical methods are used increasingly for characterization and monitoring at remediation sites in fractured-rock aquifers. The complex heterogeneity of fractured rock poses enormous challenges to groundwater remediation professionals, and new methods are needed to cost-effectively infer fracture and fracture-zone locations, orientations and properties, and to develop conceptual site models for flow and transport. Despite the potential of geophysical methods to "see" between boreholes, two issues have impeded the adoption of geophysical methods by remediation professionals. First, geophysical results are commonly only indirectly related to the properties of interest (e.g., permeability) to remediation professionals, and qualitative or quantitative interpretation is required to convert geophysical results to hydrogeologic information. Additional demonstration/evaluation projects are needed in the site remediation literature to fully transfer geophysical methods from research to practice. Second, geophysical methods are commonly viewed as inherently risky by remediation professionals. Although it is widely understood that a given method may or may not work at a particular site, the reasons are not always clear to end users of geophysical products. Synthetic modeling tools are used in research to assess the potential of a particular method to successfully image a target, but these tools are not widely used in industry. Here, we seek to advance the application of geophysical methods to solve problems facing remediation professionals with respect to fractured-rock aquifers. To this end, we (1) provide an overview of geophysical methods applied to characterization and monitoring of fractured-rock aquifers; (2) review case studies showcasing different geophysical methods; and (3) discuss best practices for method selection and rejection based on synthetic modeling and decision support tools.

Geophysical signatures of soil AFFF contamination from spectral induced polarization and low field nuclear magnetic resonance methods.

Few field methods are available for characterizing source zones impacted with aqueous film forming foam (AFFF). Non-invasive geophysical characterization of AFFF source zone contamination in situ could assist with the delineation and characterization of these sites, allowing for more informed sampling regimes aimed at quantifying subsurface poly- and perfluoroalkyl substances (PFAS) contamination. We present initial results from the investigation of the sensitivity of two existing surface and borehole-deployable geophysical technologies, spectral induced polarization (SIP), and low field nuclear magnetic resonance (NMR), to soils impacted with AFFF. To investigate the sensitivity of these methods to AFFF-impacted soil, bench-scale column experiments were conducted on samples consisting of natural and synthetic soils and groundwater. While our findings do not show strong evidence of NMR sensitivity to soil PFAS contamination, we do find evidence that SIP has sufficient sensitivity to detect sorption of AFFF constituents (including PFAS) to soils. This finding is based on evidence that AFFF constituents associated with the pore surface produce a measurable polarization response in both freshly impacted synthetic soils and in soils historically impacted with AFFF. Our findings encourage further exploration of the SIP method as a technology for characterizing contaminant concentrations across AFFF source zones.

Field evaluation of semi-automated moisture estimation from geophysics using machine learning.

Geophysical methods can provide three-dimensional (3D), spatially continuous estimates of soil moisture. However, point-to-point comparisons of geophysical properties to measure soil moisture data are frequently unsatisfactory, resulting in geophysics being used for qualitative purposes only. This is because (1) geophysics requires models that relate geophysical signals to soil moisture, (2) geophysical methods have potential uncertainties resulting from smoothing and artifacts introduced from processing and inversion, and (3) results from multiple geophysical methods are not easily combined within a single soil moisture estimation framework. To investigate these potential limitations, an irrigation experiment was performed wherein soil moisture was monitored through time, and several surface geophysical datasets indirectly sensitive to soil moisture were collected before and after irrigation: ground penetrating radar, electrical resistivity tomography (ERT), and frequency domain electromagnetics (FDEM). Data were exported in both raw and processed form, and then snapped to a common 3D grid to facilitate moisture prediction by standard calibration techniques, multivariate regression, and machine learning. A combination of inverted ERT data, raw FDEM, and inverted FDEM data was most informative for predicting soil moisture using a random regression forest model (one-thousand 60/40 training/test cross-validation folds produced root mean squared errors ranging from 0.025-0.046 cm3/cm3). This cross-validated model was further supported by a separate evaluation using a test set from a physically separate portion of the study area. Machine learning was conducive to a semi-automated model-selection process that could be used for other sites and datasets to locally improve accuracy.

GW/SW-MST: A Groundwater/Surface-Water Method Selection Tool.

Groundwater/surface-water (GW/SW) exchange and hyporheic processes are topics receiving increasing attention from the hydrologic community. Hydraulic, chemical, temperature, geophysical, and remote sensing methods are used to achieve various goals (e.g., inference of GW/SW exchange, mapping of bed materials, etc.), but the application of these methods is constrained by site conditions such as water depth, specific conductance, bed material, and other factors. Researchers and environmental professionals working on GW/SW problems come from diverse fields and rarely have expertise in all available field methods; hence there is a need for guidance to design field campaigns and select methods that both contribute to study goals and are likely to work under site-specific conditions. Here, we present the spreadsheet-based GW/SW-Method Selection Tool (GW/SW-MST) to help practitioners identify methods for use in GW/SW and hyporheic studies. The GW/SW-MST is a Microsoft Excel-based decision support tool in which the user selects answers to questions about GW/SW-related study goals and site parameters and characteristics. Based on user input, the tool indicates which methods from a toolbox of 32 methods could potentially contribute to achieving the specified goals at the site described.

DTSGUI: A Python Program to Process and Visualize Fiber-Optic Distributed Temperature Sensing Data.

Fiber-optic distributed temperature sensing (FO-DTS) has proven to be a transformative technology for the hydrologic sciences, with application to diverse problems including hyporheic exchange, groundwater/surface-water interaction, fractured-rock characterization, and cold regions hydrology. FO-DTS produces large, complex, and information-rich datasets. Despite the potential of FO-DTS, adoption of the technology has been impeded by lack of tools for data processing, analysis, and visualization. New tools are needed to efficiently and fully capitalize on the information content of FO-DTS datasets. To this end, we present DTSGUI, a public-domain Python-based software package for editing, parsing, processing, statistical analysis, georeferencing, and visualization of FO-DTS data.

Geophysical Assessment of a Proposed Landfill Site in Fredericktown, Missouri.

In cooperation with the U.S. Environmental Protection Agency (EPA), the U.S. Geological Survey (USGS) collected numerous types of geophysical data at a site in Fredericktown, Missouri, in June 2018. Various geophysical surveys were collectively used to help evaluate the overall suitability of the site for use as a mine waste-soil and sediment repository, and to evaluate the suite of geophysical methods for rapid collection and preliminary assessment of sites with shallow sediments. Land-based geophysical methods, which included frequency-domain electromagnetic induction (FDEM), electrical resistivity tomography (ERT), horizontal-to-vertical spectral ratio passive seismic (HVSR), and shear-wave refraction, were used to determine the depths to crystalline bedrock and characterize the overlying unconsolidated sediments (or regolith). Water-borne FDEM profiles and forward-looking infrared (FLIR) thermal image surveys were conducted along the Fredericktown City Lake shoreline to identify locations of potential interactions between groundwater and surface water. Sediment temperature profilers were installed at two locations along the shoreline to characterize shallow unconsolidated sediment thermal properties and support the interpretation of the other geophysical surveys. Geophysical reconnaissance methods including the FDEM and HVSR methods, were used to rapidly evaluate the vertical and lateral extent of overburden, or unconsolidated sediments, overlying the bedrock at the site. The results of these methods were compared to reference geophysical methods of ERT and shear-wave refraction surveys that have greater accuracy and are more labor intensive and time-consuming. A goal of the project was the evaluation of the validity and reliability of this suite of reconnaissance geophysical methods as a means by which shallow (less than 3 meters (m)) sediments can be rapidly assessed. Two orthogonal ERT survey profiles, which used 28 electrodes spaced 1 m apart in dipole-dipole and combined Wenner-Schlumberger configurations, were collected to determine the subsurface resistivity. The results were inverted to produce electrical resistivity profiles that were compared to the FDEM and HVSR survey results. The FDEM data were collected along cleared paths through the proposed disposal cell locations. The data were inverted to generate depth-dependent estimates of electrical conductivity along the transects. An analysis of the depth of investigation (DOI) indicated the FDEM imaged to depths of about 3 m below land surface. The ERT, FDEM, and HVSR indicated the depth to crystalline bedrock was approximately 1.5 m below land surface with shallower and deeper areas. Results from this investigation indicate this suite of methods will likely perform well at sites with shallow depths to bedrock and strong conductivity and acoustic impedance contrasts, where the FDEM and HVSR methods can provide estimates of the depth to bedrock, and ERT and shear-wave refraction surveys might not be worth the added time and expense.

Case Histories of GPR for Animal Burrows Mapping and Geometry.

Ground-penetrating radar (GPR) has a wide range of applications, from geologic mapping to concrete inspection. A recently emerging GPR application is deployment in biological investigations as a non-invasive technique. Geophysical mapping of features such as tree roots and turtle burrows has proved valuable for the understanding of these subsurface systems for ecological, environmental, or engineering purposes. Four case histories of GPR investigations pertaining to animal burrows are described: cutter ants in Brazil, groundhogs in Michigan, and groundhogs, and burrowing bees in Pennsylvania. Cutter ants (Atta spp.) in Amazonian Brazil are known to construct burrows of nearly the same dimensions as groundhogs as they excavate galleries up to 7 m deep for leaf storage. Cutter ant burrows are hazardous to heavy equipment and may also cause loss of mud circulation during rotary drilling. Groundhogs (Marmota monax), found throughout the United States, cause unseen hazards, particularly for equestrian facilities where a sudden collapse can cause severe injuries to both horse and rider. Burrowing bees (Colletes inaequalis) are common in the northeastern United States. The size of the bee burrows is significantly smaller than that of the cutter ants and the groundhogs. The data for these surveys were collected over a twenty-year span, crossing several generations of survey equipment and processing techniques. Together, these four case histories highlight the historic and current capabilities of GPR systems applied to mapping subsurface burrow systems. These examples demonstrate the important impact near surface heterogeneities have in altering ecological, environmental, or engineering systems and the utility of GPR for mapping such heterogeneities.

The Dual-Domain Porosity Apparatus: Characterizing Dual Porosity at the Sediment/Water Interface.

The characterization of pore-space connectivity in porous media at the sediment/water interface is critical in understanding contaminant transport and reactive biogeochemical processes in zones of groundwater and surface-water exchange. Previous in situ studies of dual-domain (i.e., mobile/less-mobile porosity) systems have been limited to solute tracer injections at scales of meters to hundreds of meters and subsequent numerical model parameterization using fluid concentration histories. Pairing fine-scale (e.g., sub-meter) geoelectrical measurements with fluid tracer data over time alleviates dependence on flowpath-scale experiments, enabling spatially targeted characterization of shallow sediment/water interface media where biogeochemical reactivity is often high. The Dual-Domain Porosity Apparatus is a field-tested device capable of variable rate-controlled downward flow experiments. The Dual-Domain Porosity Apparatus facilitates inference of dual-domain parameters, i.e., mobile/less-mobile exchange rate coefficient and the ratio of less mobile to mobile porosity. The Dual-Domain Porosity Apparatus experimental procedure uses water electrical conductivity as a conservative tracer of differential loading and flushing of pore spaces within the region of measurement. Variable injection rates permit the direct quantification of the flow-dependence of dual-domain parameters, which has been theorized for decades but remains challenging to assess using existing experimental methodologies.

MoisturEC: A New R Program for Moisture Content Estimation from Electrical Conductivity Data.

Noninvasive geophysical estimation of soil moisture has potential to improve understanding of flow in the unsaturated zone for problems involving agricultural management, aquifer recharge, and optimization of landfill design and operations. In principle, several geophysical techniques (e.g., electrical resistivity, electromagnetic induction, and nuclear magnetic resonance) offer insight into soil moisture, but data-analysis tools are needed to "translate" geophysical results into estimates of soil moisture, consistent with (1) the uncertainty of this translation and (2) direct measurements of moisture. Although geostatistical frameworks exist for this purpose, straightforward and user-friendly tools are required to fully capitalize on the potential of geophysical information for soil-moisture estimation. Here, we present MoisturEC, a simple R program with a graphical user interface to convert measurements or images of electrical conductivity (EC) to soil moisture. Input includes EC values, point moisture estimates, and definition of either Archie parameters (based on experimental or literature values) or empirical data of moisture vs. EC. The program produces two- and three-dimensional images of moisture based on available EC and direct measurements of moisture, interpolating between measurement locations using a Tikhonov regularization approach.

Hydrogeochemical controls on brook trout spawning habitats in a coastal stream.

Brook trout (Salvelinus fontinalis) spawn in fall and overwintering egg development can benefit from stable, relatively warm temperatures in groundwater-seepage zones. However, eggs are also sensitive to dissolved oxygen concentration, which may be reduced in discharging groundwater (i.e., seepage). We investigated a 2 km reach of the coastal Quashnet River in Cape Cod, Massachusetts, USA, to relate preferred fish spawning habitats to geology, geomorphology, and discharging groundwater geochemistry. Thermal reconnaissance methods were used to locate zones of rapid groundwater discharge, which were predominantly found along the central channel of a wider stream valley section. Pore-water chemistry and temporal vertical groundwater flux were measured at a subset of these zones during field campaigns over several seasons. Seepage zones in open-valley sub-reaches generally showed suboxic conditions and higher dissolved solutes compared to the underlying glacial outwash aquifer. These discharge zones were cross-referenced with preferred brook trout redds and evaluated during 10 years of observation, all of which were associated with discrete alcove features in steep cutbanks, where stream meander bends intersect the glacial valley walls. Seepage in these repeat spawning zones was generally stronger and more variable than in open-valley sites, with higher dissolved oxygen and reduced solute concentrations. The combined evidence indicates that regional groundwater discharge along the broader valley bottom is predominantly suboxic due to the influence of near-stream organic deposits; trout show no obvious preference for these zones when spawning. However, the meander bends that cut into sandy deposits near the valley walls generate strong oxic seepage zones that are utilized routinely for redd construction and the overwintering of trout eggs. Stable water isotopic data support the conclusion that repeat spawning zones are located directly on preferential discharges of more localized groundwater. In similar coastal systems with extensive valley peat deposits, the specific use of groundwater-discharge points by brook trout may be limited to morphologies such as cutbanks, where groundwater flow paths do not encounter substantial buried organic material and remain oxygen-rich.

Specific polarizability of sand-clay mixtures with varying ethanol concentration.

We utilise a concept of specific polarizability (c s ), represented as the ratio of mineral-fluid interface polarization per pore-normalised surface area S p , to demonstrate the influence of clay-organic interaction on complex conductivity measurements. Complex conductivity measurements were performed on kaolinite- and illite-sand mixtures as a function of varying ethanol (EtOH) concentration (10% and 20% v/v). The specific surface area of each clay type and Ottawa sand was determined by nitrogen-gas-adsorption Brunauer-Emmett-Teller method. We also calculated the porosity and saturation of each mixture based on weight loss of dried samples. Debye decomposition, a phenomenological model, was applied to the complex conductivity data to determine normalised chargeability (m n ). Specific polarizability estimates from previous complex conductivity measurements for bentonite-sand mixtures were compared with our dataset. The c s for all sand-clay mixtures decreased as the EtOH concentration increased from 0% to 10% to 20% v/v. We observe similar c s responses to EtOH concentration for all sand-clay mixtures. Analysis of variance with a level of significance α = 0.05 suggests that the suppression in c s responses with increasing EtOH concentration was statistically significant for all sand-clay mixtures. On the other hand, real conductivity showed only 10% to 20% v/v changes with increasing EtOH concentration. The c s estimates reflect the sensitivity of complex conductivity measurements to alteration in surface chemistry at available surface adsorption sites for different clay types, likely resulting from ion exchange at the clay surface and associated with kinetic reactions in the electrical double layer of the clay-water-EtOH media. Our results indicate a much larger influence of specific surface area and ethanol concentration on clay-driven polarization relative to changes in clay mineralogy.

Heat as a groundwater tracer in shallow and deep heterogeneous media: Analytical solution, spreadsheet tool, and field applications.

Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction-dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson (2004), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater-surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux-LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi-decadal surface signals have disturbed subsurface thermal regimes.

Scenario Evaluator for Electrical Resistivity Survey Pre-modeling Tool.

Geophysical tools have much to offer users in environmental, water resource, and geotechnical fields; however, techniques such as electrical resistivity imaging (ERI) are often oversold and/or overinterpreted due to a lack of understanding of the limitations of the techniques, such as the appropriate depth intervals or resolution of the methods. The relationship between ERI data and resistivity is nonlinear; therefore, these limitations depend on site conditions and survey design and are best assessed through forward and inverse modeling exercises prior to field investigations. In this approach, proposed field surveys are first numerically simulated given the expected electrical properties of the site, and the resulting hypothetical data are then analyzed using inverse models. Performing ERI forward/inverse modeling, however, requires substantial expertise and can take many hours to implement. We present a new spreadsheet-based tool, the Scenario Evaluator for Electrical Resistivity (SEER), which features a graphical user interface that allows users to manipulate a resistivity model and instantly view how that model would likely be interpreted by an ERI survey. The SEER tool is intended for use by those who wish to determine the value of including ERI to achieve project goals, and is designed to have broad utility in industry, teaching, and research.

1DTempPro V2: New Features for Inferring Groundwater/Surface-Water Exchange.

A new version of the computer program 1DTempPro extends the original code to include new capabilities for (1) automated parameter estimation, (2) layer heterogeneity, and (3) time-varying specific discharge. The code serves as an interface to the U.S. Geological Survey model VS2DH and supports analysis of vertical one-dimensional temperature profiles under saturated flow conditions to assess groundwater/surface-water exchange and estimate hydraulic conductivity for cases where hydraulic head is known.

1DTempPro: analyzing temperature profiles for groundwater/surface-water exchange.

A new computer program, 1DTempPro, is presented for the analysis of vertical one-dimensional (1D) temperature profiles under saturated flow conditions. 1DTempPro is a graphical user interface to the U.S. Geological Survey code Variably Saturated 2-Dimensional Heat Transport (VS2DH), which numerically solves the flow and heat-transport equations. Pre- and postprocessor features allow the user to calibrate VS2DH models to estimate vertical groundwater/surface-water exchange and also hydraulic conductivity for cases where hydraulic head is known.

Electrical signatures of ethanol-liquid mixtures: implications for monitoring biofuels migration in the subsurface.

Ethanol (EtOH), an emerging contaminant with potential direct and indirect environmental effects, poses threats to water supplies when spilled in large volumes. A series of experiments was directed at understanding the electrical geophysical signatures arising from groundwater contamination by ethanol. Conductivity measurements were performed at the laboratory scale on EtOH-water mixtures (0 to 0.97 v/v EtOH) and EtOH-salt solution mixtures (0 to 0.99 v/v EtOH) with and without a sand matrix using a conductivity probe and a four-electrode electrical measurement over the low frequency range (1-1000 Hz). A Lichtenecker-Rother (L-R) type mixing model was used to simulate electrical conductivity as a function of EtOH concentration in the mixture. For all three experimental treatments increasing EtOH concentration resulted in a decrease in measured conductivity magnitude (|σ|). The applied L-R model fitted the experimental data at concentration ≤0.4v/v EtOH, presumably due to predominant and symmetric intermolecular (EtOH-water) interaction in the mixture. The deviation of the experimental |σ| data from the model prediction at higher EtOH concentrations may be associated with hydrophobic effects of EtOH-EtOH interactions in the mixture. The |σ| data presumably reflected changes in relative strength of the three types of interactions (water-water, EtOH-water, and EtOH-EtOH) occurring simultaneously in EtOH-water mixtures as the ratio of EtOH to water changed. No evidence of measurable polarization effects at the EtOH-water and EtOH-water-mineral interfaces over the investigated frequency range was found. Our results indicate the potential for using electrical measurements to characterize and monitor EtOH spills in the subsurface.

Complex resistivity signatures of ethanol biodegradation in porous media.

Numerous adverse effects are associated with the accidental release of ethanol (EtOH) and its persistence in the subsurface. Geophysical techniques may permit non-invasive, real time monitoring of microbial degradation of hydrocarbon. We performed complex resistivity (CR) measurements in conjunction with geochemical data analysis on three microbial-stimulated and two control columns to investigate changes in electrical properties during EtOH biodegradation processes in porous media. A Debye Decomposition approach was applied to determine the chargeability (m), normalized chargeability (m(n)) and time constant (τ) of the polarization magnitude and relaxation length scale as a function of time. The CR responses showed a clear distinction between the bioaugmented and control columns in terms of real (σ') and imaginary (σ″) conductivity, phase (ϕ) and apparent formation factor (F(app)). Unlike the control columns, a substantial decrease in σ' and increase in F(app) occurred at an early time (within 4 days) of the experiment for all three bioaugmented columns. The observed decrease in σ' is opposite to previous studies on hydrocarbon biodegradation. These columns also exhibited increases in ϕ (up to ~9 mrad) and σ″ (up to two order of magnitude higher) 5 weeks after microbial inoculation. Variations in m and m(n) were consistent with temporal changes in ϕ and σ″ responses, respectively. Temporal geochemical changes and high resolution scanning electron microscopy imaging corroborated the CR findings, thus indicating the sensitivity of CR measurements to EtOH biodegradation processes. Our results offer insight into the potential application of CR measurements for long-term monitoring of biogeochemical and mineralogical changes during intrinsic and induced EtOH biodegradation in the subsurface.

Complex resistivity signatures of ethanol in sand-clay mixtures.

We performed complex resistivity (CR) measurements on laboratory columns to investigate changes in electrical properties as a result of varying ethanol (EtOH) concentration (0% to 30% v/v) in a sand-clay (bentonite) matrix. We applied Debye decomposition, a phenomenological model commonly used to fit CR data, to determine model parameters (time constant: τ, chargeability: m, and normalized chargeability: mn). The CR data showed a significant (P≤0.001) time-dependent variation in the clay driven polarization response (~12 mrad) for 0% EtOH concentration. This temporal variation probably results from the clay-water reaction kinetics trending towards equilibrium in the sand-clay-water system. The clay polarization is significantly suppressed (P≤0.001) for both measured phase (ϕ) and imaginary conductivity (σ″) with increasing EtOH concentration. Normalized chargeability consistently decreases (by up to a factor of ~2) as EtOH concentration increases from 0% to 10% and 10 to 20%, respectively. We propose that such suppression effects are associated with alterations in the electrical double layer (EDL) at the clay-fluid interface due to (a) strong EtOH adsorption on clay, and (b) complex intermolecular EtOH-water interactions and subsequent changes in ionic mobility on the surface in the EDL. Changes in the CR data following a change of the saturating fluid from EtOH 20% to plain water indicate strong hysteresis effects in the electrical response, which we attribute to persistent EtOH adsorption on clay. Our results demonstrate high sensitivity of CR measurements to clay-EtOH interactions in porous media, indicating the potential application of this technique for characterization and monitoring of ethanol contamination in sediments containing clays.