Heat transport into active layer soils is important to understanding potential responses to changes in surface energy balance, particularly in the context of changing climate. Here we present results of a study to characterise soil thermal properties along a soil moisture gradient adjacent to Lake Fryxell in Taylor Valley, Antarctica. Our goals were to characterise the thermal characteristics of these relatively wet soils (compared to the rest of the McMurdo Dry Valleys landscape), and to assess the response of the active layer to possible increases in soil moisture. We measured subsurface temperatures at depths from 3 to 50 cm at four locations along a natural gradient of wet to dry soils adjacent to Lake Fryxell from January 2006 to January 2007. We used a numerical model to estimate apparent thermal diffusivity (ATD) and simulate observed temperature time series. Calculations of ATD at discrete locations yielded values ranging from 1.0 × 10−9 – 2.4 × 10−5 m2 s−1, and the corresponding range of bulk (i.e. depth averaged at a single surface location) ATD was 2.9 × 10−9–1.2 × 10−7 m2 s−1. Thawed soils had a range of bulk ATD during warming of 2.9 × 10−9–3.8 × 10−8 m2 s−1, and during cooling of 2.9 × 10−9–4.8 × 10−8 m2 s−1. When soils were frozen, however, the range of bulk ATD was 7.6 × 10−9–1.2 × 10−7 m2 s−1 during warming, and 7.8 × 10−9–1.1 × 10−7 m2 s−1 during cooling. Estimated bulk ATD values were consistently greater in locations of enhanced soil moisture, so lakeside soils were more likely to conduct energy into the subsurface. Increased soil moisture across the landscape would likely increase ATD, allowing for greater heat exchange between the atmosphere and the subsurface. Copyright © 2009 John Wiley & Sons, Ltd.

[1] A method is proposed to localize preferential fluid flow pathways in porous media on the basis of time-lapse self-potential measurements associated with salt tracer injection upstream. This method is first tested using laboratory data. A network of nonpolarizing electrodes located is connected to a highly sensitive voltmeter used to record the resulting electrical field fluctuations occurring over time at the surface of the tank. The transport of the conductive salt plume through the permeable porous materials changes the localized streaming potential coupling coefficient associated with the advective drag of the excess charge of the pore water and is also responsible for a diffusion current associated with the salinity gradient. Monitoring of the electrical potential distribution at the ground surface can be used to localize the pulse of saline water over time and to determine its velocity. This method applies in real time and can be used to track highly localized flow pathways characterized by high permeability. Our sandbox experiment demonstrates the applicability of this new method under well-controlled conditions with a coarse-sand channel embedded between fine-sand banks. A finite element model allows us to reproduce the time-lapse electrical potential distribution over the channel, but some discrepancies were observed on the banks. Finally, we performed a numerical simulation for a synthetic case study inspired by a recently published field case study. A Markov chain Monte Carlo sampler is used to determine the permeability and the porosity of the preferential fluid flow pathway of this synthetic case study.

Resistivity and self-potential tomography can be used to investigate anomalous seepage inside heterogeneous earthen dams. The self-potential (SP) signals provide a unique signature to groundwater flow because the source current density responsible for the SP signals is proportional to the Darcy velocity. The distribution of the SP signals is also influenced by the distribution of the resistivity; therefore, resistivity and SP need to be used in concert to elucidate groundwater flow pathways. In this study, a survey is conducted at a small earthen dam in Colorado where anomalous seepage is observed on the downstream face at the dam toe. The data reveal SP and direct current resistivity anomalies that are used to delineate three anomalous seepage zones within the dam and to estimate the source of the localized seepage discharge. The SP data are inverted in two dimensions using the resistivity distribution to determine the distribution of the Darcy velocity responsible for the observed seepage. The inverted Darcy velocity agrees with an estimation of the Darcy velocity from the hydraulic conductivity obtained from a slug test and the observed head gradient.

In south-central Texas the lower Guadalupe River has incised into the outcrop of the Carrizo-Wilcox aquifer. The river and the aquifer are hydraulically connected across the outcrop, although the connectivity is obscured at the surface by alluvium and surface-water and groundwater exchange dynamics are currently poorly understood. To investigate surface-water and groundwater exchange dynamics between the lower Guadalupe River and the Carrizo-Wilcox aquifer, a geophysical study was completed along a 14.86 km reach of the river by using water-borne gradient self-potential (SP) profiling and two-dimensional direct-current electric resistivity tomography. This paper explores the applicability of these water-borne geoelectric methods in delineating gaining and losing channel reaches, and demonstrates that geoelectric signals in the form of total electric field strength can be logged with an electric dipole and decomposed into component SP signals depicting regional and local groundwater flow patterns attributable to regional and localized hydraulic gradients. Localized SP anomalies of several tens of millivolts, indicative of hyporheic exchange flows, are observed and superimposed upon a 124 mV regional SP anomaly indicative of ambient groundwater exchange flows between the river and the aquifer. The observed SP signals are interpreted through two-dimensional finite-element modeling of streaming potentials attributable to ambient groundwater exchange and hyporheic exchange flow patterns. Variables of the channel environment such as temperature and concentration gradients, depth, and velocity are considered and subsequently eliminated as alternative sources of the SP signals that are presented.

A hydrogeophysical survey is performed at small earthen dam that overlies a confined aquifer. The structure of the dam has not shown evidence of anomalous seepage internally or through the foundation prior to the survey. However, the surface topography is mounded in a localized zone 150 m downstream, and groundwater discharges from this zone periodically when the reservoir storage is maximum. We use self-potential and electrical resistivity tomography surveys with seismic refraction tomography to (1) determine what underlying hydrogeologic factors, if any, have contributed to the successful long-term operation of the dam without apparent indicators of anomalous seepage through its core and foundation; and (2) investigate the hydraulic connection between the reservoir and the seepage zone to determine whether there exists a potential for this success to be undermined. Geophysical data are informed by hydraulic and geotechnical borehole data. Seismic refraction tomography is performed to determine the geometry of the phreatic surface. The hydro-stratigraphy is mapped with the resistivity data and groundwater flow patterns are determined with self-potential data. A self-potential model is constructed to represent a perpendicular profile extending out from the maximum cross-section of the dam, and self-potential data are inverted to recover the groundwater velocity field. The groundwater flow pattern through the aquifer is controlled by the bedrock topography and a preferential flow pathway exists beneath the dam. It corresponds to a sandy-gravel layer connecting the reservoir to the downstream seepage zone.