Geology topics

Frederick D. Day-Lewis

Publications and source records attributed to Frederick D. Day-Lewis.

At least 19 recordsLinked to original sources

USGS2021

Nuclear magnetic resonanance logs of fractured bedrock at the Hidden Lane Landfill site, Culpeper Basin, Virginia

In May 2018, the U.S. Geological Survey (USGS) in cooperation with the U.S. Environmental Protection Agency (EPA) collected borehole nuclear magnetic resonance (bNMR) logs in three boreholes completed in sandstone and siltstone of the Balls Bluff Member of the Bull Run Formation at a Superfund Site in Culpeper Basin, Virginia. The bNMR logs were used to aid in the evaluation of the aquifer by measuring the porosity, determining the mobile and immobile fractions of water, and estimating the hydraulic conductivity, to evaluate the potential storage and transport properties at the site. The bNMR method measures the transverse (T 2 ) decay in nuclear magnetism in response to radio-frequency pulses. The relaxation decay is related to water content and the size of the pores where the water resides. In addition, the relaxation decay parameters are used to estimate hydraulic conductivity. The results were compared to other borehole logs collected at the site and to regional groundwater investigations in similar rock formations.

Virginia↗

USGS2020

Formation criteria for hyporheic anoxic microzones: Assessing interactions of hydraulics, nutrients and biofilms

Recent experimental studies have detected the presence of anoxic microzones in hyporheic sediments. These microzones are small‐scale anoxic pores, embedded within oxygen‐rich porous media and can act as anaerobic reaction sites producing reduction compounds such as nitrous oxide, a potent greenhouse gas. Microbes are a key control on nutrient transformation in hyporheic sediment, but their associated biomass growth is also capable of altering hydraulic flux, leading to potential bioclogging. Here, we developed one of the first computational modeling approaches that combined hydraulics and microbial conditions to explore the continuous evolution of microzones in stream sediments. The model assessed stream and sediment conditions with different hydraulic flux (0.1–1.0 m/day Darcy flux), nutrient concentrations (O 2 = 8 mg/L, OrgC = 20 mg/L, NO − 3 = 1.5–3 mg/L, and NH 3 = 0.5–1 mg/L), and biomass scenarios (with and without). The model domain is a pore network model with random sized pore‐throat radii creating heterogeneous and anisotropic flow that is representative of a natural streambed composed of medium sand with a hydraulic conductivity of 0.8 m/day. Results from 30 day simulations indicate that hyporheic microzone formation will occur and microzone distributions are not simply controlled by residence time alone, rather by the complex interactions of hydraulic flux, nutrient concentrations, and biomass, with bioclogging having strong feedbacks on both hydraulics and nutrients. Under all conditions with biomass growth, anoxic microzones were unstable, perishing a few days after formation, because bioclogging primarily occurs near the influent (downwelling) area of the hyporheic zone. In turn, this bioclogging shifts transport conditions from advection‐dominated to diffusion‐dominated transport, removing all oxic regions in the hyporheic zone. Overall, results from the modeling show that anoxic microzones are likely to form under many hyporheic zone conditions, and be dynamic through space and time as they are dependent on both hydraulic flux and nutrient transport.

Water Resources Research↗

USGS2020

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.

Groundwater↗

USGS2019

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

Journal of Geophysical Research: Biogeosciences↗

USGS2019

Multi-scale preferential flow processes in an urban streambed under variable hydraulic conditions

Spatially preferential flow processes occur at nested scales at the sediment-water interface (SWI), due in part to sediment heterogeneities, which may be enhanced in flashy urban streams with heavy road sand influence. However, several factors, including the flow-rate dependence of preferential hyporheic flow and discrete groundwater discharge zones are commonly overlooked in reach-scale models of groundwater/surface water exchange. Using a series of controlled-head tracer-injection experiments coupled with cm-scale geophysics within the highly reactive upper 30 cm of the hyporheic zone of an urban stream, we quantified the flow dependence of local less-mobile porosity volume, mass-transfer rate coefficient, and the resulting local residence time in the less-mobile pore space at three controlled downward fluid fluxes (0.8, 2, and 3 m/d). Experiments were performed in two adjacent streambed locations, representing different sediment bulk vertical permeability. Less-mobile porosity parameters were generally substantial and similar between the two streambed locations; though a more competent, thin, organic layer at ∼15 cm depth in one location strongly impacted tracer loading, flushing dynamics, and local residence times. Increased downward flux led to (1) a decrease in less-mobile porosity residence time in all experiments, and (2) an increase in less-mobile porosity fraction for most experiments. Additionally, at the larger stream reach-scale, surface electrodes for electrical resistivity measurement were installed along 22 m of the wetted stream channel. These surface electrode measurements were collected during a natural storm flow event, which revealed widespread, short-term, flushing (e.g. <3 h) of the hyporheic zone with stream water, followed by longer-term (e.g. >60 h) flushing of the SWI with riparian zone groundwater. Flow dependence of preferential hyporheic zone flowpaths, like in the controlled tracer experiments, was also observed in these reach-scale electrical resistivity tomography measurements. Our findings reveal that the spatial and temporal dependence of preferential flow processes create highly dynamic SWI conditions that will affect the physical and coupled biogeochemical functions of the SWI in urbanized, sand-impacted streams.

Journal of Hydrology↗

USGS2019

Geophysical mapping of plume discharge to surface water at a crude oil spill site: Inversion versus machine learning

The interpretation of geophysical survey results to answer hydrologic, engineering, and geologic questions is critical to diverse problems for management of water, energy, and mineral resources. Although geophysical images provide valuable qualitative insight into subsurface architecture and conditions, translating geophysical images into quantitative information (e.g., saturation, concentration, and hydraulic properties) often involves substantial nonuniqueness and uncertainty owing to the limited resolution of geophysical imaging and uncertainty in petrophysical relations. We have developed a machine-learning approach to address these challenges in the context of a field-based investigation to map zones where a hydrocarbon plume was discharging to surface water at the National Crude Oil Spill Fate and Natural Attenuation Research Site in Bemidji, Minnesota, USA. The two-step approach combines multiple types of geophysical and direct information and effectively bypasses inversion and its associated assumptions. Integrating multifrequency electromagnetic induction, ground-penetrating radar, and fluid-sampling data, we first identify discharge zones and second estimate specific conductance versus depth. Compared with conventional inversion results, the machine-learning results (1) directly address the study objectives (delineating the discharge zones); (2) better extract depth-dependent information from the data, for which sensitivity diminishes rapidly with depth; and (3) quantify the uncertainty of the predictions (i.e., discharge versus nondischarge zones), rather than the uncertainty of the geophysical estimates (i.e., the standard error of estimation for the logarithm of electrical conductivity).

Minnesota↗

USGS2018

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 to 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) studies have been limited to solute tracer injections at scales of meters to 100s 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 meter-scale 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.

Groundwater↗

USGS2018

Direct observations of hydrologic exchange occurring with less‐mobile porosity and the development of anoxic microzones in sandy lakebed sediments

Quantifying coupled mobile/less‐mobile porosity dynamics is critical to the prediction of biogeochemical storage, release, and transformation processes in the zone where groundwater and surface water exchange. The recent development of fine‐scale geoelectrical monitoring paired with pore‐water sampling in groundwater systems enables direct characterization of hydrologic exchange between more‐ and less‐mobile porosity during tracer tests. We adapt this technique to sandy interface sediments at a groundwater flow‐through kettle lake. Tracer experiments were conducted within controlled‐head permeameters over a range of specified downward flow conditions over several days. Although the bed was predominantly composed of highly permeable sands and gravels, cobble inclusions created less‐mobile flow zones at the centimeter scale. Less‐mobile porosity fractions, residence times, and rates of exchange were inferred from paired bulk and fluid electrical conductivity data, without the need for inverse model calibration. The conservative solute experiments were paired with 15 NO 3 − and other reactive amendments, revealing anaerobic processes occurring at shallow sediment depths where pore‐water sampling indicated bulk‐oxic conditions. The average less‐mobile porosity residence times as evaluated with the geoelectrical method were on 1‐hr timescales, which appear to be biogeochemically important in the context of creating anoxic microzones within less‐mobile porosity of sandy interface sediments.

Water Resources Research↗

USGS2018

Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media

Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.

Hydrological Processes↗

USGS2018

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.

Groundwater↗

USGS2018

An on-campus well field for hydrogeophysics education and undergraduate research

The emerging subdiscipline of hydrogeophysics is underdeveloped in undergraduate curricula relative to its importance in professional engineering/environmental practice. In 2001, the Bucknell Department of Geology initiated efforts to refocus an undergraduate geophysics course on near-surface geophysical methods for hydrologic, environmental, and engineering problems. In addition to offering students practical experience, treatment of hydrogeophysics provides important pedagogical opportunities. Field-based hydrogeophysics labs challenge students to integrate concepts from other geology courses, as well as from physics, math, and chemistry. We faced two challenges in revising our geophysics course: (1) access to wells for field exercises on borehole geophysics; and (2) the costs of acquiring and maintaining equipment. We pursued two strategies to solve these problems. First, we established an on-campus well field, which serves as a field laboratory for downhole and cross-hole experiments. Second, we incorporated field demonstrations and lectures by professional geoscientists, including alumni, into our courses. By adding field exercises to our syllabi and promoting undergraduate research, we are building a cutting-edge dataset that includes televiewer and standard wellbore logging, cross-hole tomography, and hydraulic-test data. Student-led analysis of these data has already provided valuable insights into the control of fractures on aquifer properties, and these observations are being combined with outcrop studies to place our results in a more regional context.

Pennsylvania↗

USGS2018

Evaluating long-term patterns of decreasing groundwater discharge through a lake-bottom permeable reactive barrier

Identifying and quantifying groundwater exchange is critical when considering contaminant fate and transport at the groundwater/surface-water interface. In this paper, areally distributed temperature and point seepage measurements are used to efficiently assess spatial and temporal groundwater discharge patterns through a glacial-kettle lakebed area containing a zero-valent iron permeable reactive barrier (PRB). Concern was that the PRB was becoming less permeable with time owing to biogeochemical processes within the PRB. Patterns of groundwater discharge over an 8-year period were examined using fiber-optic distributed temperature sensing (FO-DTS) and snapshot-in-time point measurements of temperature. The resulting thermal maps show complex and uneven distributions of temperatures across the lakebed and highlight zones of rapid seepage near the shoreline and along the outer boundaries of the PRB. Repeated thermal mapping indicates an increase in lakebed temperatures over time at periods of similar stage and surface-water temperature. Flux rates in six seepage meters permanently installed on the lakebed in the PRB area decreased on average by 0.021 md -1 (or about 4.5 percent) annually between 2004 and 2015. Modeling of diurnal temperature signals from shallow vertical profiles yielded mean flux values ranging from 0.39 to 1.15 md -1 , with stronger fluxes generally related to colder lakebed temperatures. The combination of an increase in lakebed temperatures, declines in direct seepage, and observations of increased cementation of the lakebed surface provide in situ evidence that the permeability of the PRB is declining. The presence of temporally persistent rapid seepage zones is also discussed.

Journal of Environmental Management↗

USGS2017

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.

Groundwater↗

USGS2016

Surface geophysical methods for characterising frozen ground in transitional permafrost landscapes

The distribution of shallow frozen ground is paramount to research in cold regions, and is subject to temporal and spatial changes influenced by climate, landscape disturbance and ecosystem succession. Remote sensing from airborne and satellite platforms is increasing our understanding of landscape-scale permafrost distribution, but typically lacks the resolution to characterise finer-scale processes and phenomena, which are better captured by integrated surface geophysical methods. Here, we demonstrate the use of electrical resistivity imaging (ERI), electromagnetic induction (EMI), ground penetrating radar (GPR) and infrared imaging over multiple summer field seasons around the highly dynamic Twelvemile Lake, Yukon Flats, central Alaska, USA. Twelvemile Lake has generally receded in the past 30 yr, allowing permafrost aggradation in the receded margins, resulting in a mosaic of transient frozen ground adjacent to thick, older permafrost outside the original lakebed. ERI and EMI best evaluated the thickness of shallow, thin permafrost aggradation, which was not clear from frost probing or GPR surveys. GPR most precisely estimated the depth of the active layer, which forward electrical resistivity modelling indicated to be a difficult target for electrical methods, but could be more tractable in time-lapse mode. Infrared imaging of freshly dug soil pit walls captured active-layer thermal gradients at unprecedented resolution, which may be useful in calibrating emerging numerical models. GPR and EMI were able to cover landscape scales (several kilometres) efficiently, and new analysis software showcased here yields calibrated EMI data that reveal the complicated distribution of shallow permafrost in a transitional landscape.

Alaska↗

USGS2016

Application of frequency- and time-domain electromagnetic surveys to characterize hydrostratigraphy and landfill construction at the Amargosa Desert Research Site, Beatty, Nevada

In 2014 and 2015, the U.S. Geological Survey (USGS), conducted frequency-domain electromagnetic (FDEM) surveys at the USGS Amargosa Desert Research Site (ADRS), approximately 17 kilometers (km) south of Beatty, Nevada. The FDEM surveys were conducted within and adjacent to a closed low-level radioactive waste disposal site located at the ADRS. FDEM surveys were conducted on a grid of north-south and east-west profiles to assess the locations and boundaries of historically recorded waste-disposal trenches. In 2015, the USGS conducted time-domain (TDEM) soundings along a profile adjacent to the disposal site (landfill) in cooperation with the U.S. Environmental Protection Agency (USEPA), to assess the thickness and characteristics of the underlying deep unsaturated zone, and the hydrostratigraphy of the underlying saturated zone. FDEM survey results indicate the general location and extent of the waste-disposal trenches and reveal potential differences in material properties and the type and concentration of waste in several areas of the landfill. The TDEM surveys provide information on the underlying hydrostratigraphy and characteristics of the unsaturated zone that inform the site conceptual model and support an improved understanding of the hydrostratigraphic framework. Additional work is needed to interpret the TDEM results in the context of the local and regional structural geology.

Nevada↗

USGS2015

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.

Groundwater Monitoring & Remediation↗

USGS2015

Imaging pathways in fractured rock using three-dimensional electrical resistivity tomography

Major challenges exist in delineating bedrock fracture zones because these cause abrupt changes in geological and hydrogeological properties over small distances. Borehole observations cannot sufficiently capture heterogeneity in these systems. Geophysical techniques offer the potential to image properties and processes in between boreholes. We used three-dimensional cross borehole electrical resistivity tomography (ERT) in a 9 m (diameter) × 15 m well field to capture high-resolution flow and transport processes in a fractured mudstone contaminated by chlorinated solvents, primarily trichloroethylene. Conductive (sodium bromide) and resistive (deionized water) injections were monitored in seven boreholes. Electrode arrays with isolation packers and fluid sampling ports were designed to enable acquisition of ERT measurements during pulsed tracer injections. Fracture zone locations and hydraulic pathways inferred from hydraulic head drawdown data were compared with electrical conductivity distributions from ERT measurements. Static ERT imaging has limited resolution to decipher individual fractures; however, these images showed alternating conductive and resistive zones, consistent with alternating laminated and massive mudstone units at the site. Tracer evolution and migration was clearly revealed in time-lapse ERT images and supported by in situ borehole vertical apparent conductivity profiles collected during the pulsed tracer test. While water samples provided important local information at the extraction borehole, ERT delineated tracer migration over spatial scales capturing the primary hydrogeological heterogeneity controlling flow and transport. The fate of these tracer injections at this scale could not have been quantified using borehole logging and/or borehole sampling methods alone.

New Jersey↗

USGS2015

Anomalous solute transport in saturated porous media: Relating transport model parameters to electrical and nuclear magnetic resonance properties

The advection-dispersion equation (ADE) fails to describe commonly observed non-Fickian solute transport in saturated porous media, necessitating the use of other models such as the dual-domain mass-transfer (DDMT) model. DDMT model parameters are commonly calibrated via curve fitting, providing little insight into the relation between effective parameters and physical properties of the medium. There is a clear need for material characterization techniques that can provide insight into the geometry and connectedness of pore spaces related to transport model parameters. Here, we consider proton nuclear magnetic resonance (NMR), direct-current (DC) resistivity, and complex conductivity (CC) measurements for this purpose, and assess these methods using glass beads as a control and two different samples of the zeolite clinoptilolite, a material that demonstrates non-Fickian transport due to intragranular porosity. We estimate DDMT parameters via calibration of a transport model to column-scale solute tracer tests, and compare NMR, DC resistivity, CC results, which reveal that grain size alone does not control transport properties and measured geophysical parameters; rather, volume and arrangement of the pore space play important roles. NMR cannot provide estimates of more-mobile and less-mobile pore volumes in the absence of tracer tests because these estimates depend critically on the selection of a material-dependent and flow-dependent cutoff time. Increased electrical connectedness from DC resistivity measurements are associated with greater mobile pore space determined from transport model calibration. CC was hypothesized to be related to length scales of mass transfer, but the CC response is unrelated to DDMT.

Water Resources Research↗