John W. Lane

Time-lapse imaging of saline-tracer transport in fractured rock using difference-attenuation radar tomography

[1] Accurate characterization of fractured-rock aquifer heterogeneity remains one of the most challenging and important problems in groundwater hydrology. We demonstrate a promising strategy to identify preferential flow paths in fractured rock using a combination of geophysical monitoring and conventional hydrogeologic tests. Cross-well differenceattenuation ground-penetrating radar was used to monitor saline-tracer migration in an experiment at the U.S. Geological Survey Fractured Rock Hydrology Research Site in Grafton County, New Hampshire. Radar data sets were collected every 10 min in three adjoining planes for 5 hours during each of 12 tracer tests. An innovative inversion method accounts for data acquisition times and temporal changes in attenuation during data collection. The inverse algorithm minimizes a combination of two functions. The first is the sum of weighted squared data residuals. Second is a measure of solution complexity based on an a priori space-time covariance function, subject to constraints that limit radarattenuation changes to regions of the tomograms traversed by high difference-attenuation ray paths. The time series of tomograms indicate relative tracer concentrations and tracer arrival times in the image planes; from these we infer the presence and location of a preferential flow path within a previously identified zone of transmissive fractures. These results provide new insights into solute channeling and the nature of aquifer heterogeneity at the site. INDEX TERMS: 0910 Exploration Geophysics: Data processing; 0915 Exploration Geophysics: Downhole methods; 1829 Hydrology: Groundwater hydrology; 1832 Hydrology: Groundwater transport; 1894 Hydrology: Instruments and techniques; KEYWORDS: radar tomography, fractured rock, ground-penetrating radar, geophysics, hydrogeophysics

Delineating a shallow fault zone and dipping bedrock strata using multichannal analysis of surface waves with a land streamer

The multichannel analysis of surface waves (MASW) seismic method was used to delineate a fault zone and gently dipping sedimentary bedrock at a site overlain by several meters of regolith. Seismic data were collected rapidly and inexpensively using a towed 30-channel land streamer and a rubberband-accelerated weight-drop seismic source. Data processed using the MASW method imaged the subsurface to a depth of about 20 m and allowed detection of the overburden, gross bedding features, and fault zone. The fault zone was characterized by a lower shear-wave velocity ( Vs ) than the competent bedrock, consistent with a large-scale fault, secondary fractures, and in-situ weathering. The MASW 2D Vs section was further interpreted to identify dipping beds consistent with local geologic mapping. Mapping of shallow-fault zones and dipping sedimentary rock substantially extends the applications of the MASW method.

Radar attenuation tomography using the centroid frequency downshift method

A method for tomographically estimating electromagnetic (EM) wave attenuation based on analysis of centroid frequency downshift (CFDS) of impulse radar signals is described and applied to cross-hole radar data. The method is based on a constant-Q model, which assumes a linear frequency dependence of attenuation for EM wave propagation above the transition frequency. The method uses the CFDS to construct the projection function. In comparison with other methods for estimating attenuation, the CFDS method is relatively insensitive to the effects of geometric spreading, instrument response, and antenna coupling and radiation pattern, but requires the data to be broadband so that the frequency shift and variance can be easily measured. The method is well-suited for difference tomography experiments using electrically conductive tracers. The CFDS method was tested using cross-hole radar data collected at the U.S. Geological Survey Fractured Rock Research Site at Mirror Lake, New Hampshire (NH) during a saline-tracer injection experiment. The attenuation-difference tomogram created with the CFDS method outlines the spatial distribution of saline tracer within the tomography plane.

Continuous resistivity profiling to delineate submarine groundwater discharge—examples and limitations

Aquifer-ocean interaction, saline intrusion, and submarine groundwater discharge (SGD) are emerging topics in hydrology and oceanography with important implications for water-resource management and estuarine ecology. Although the threat of saltwater intrusion has long been recognized in coastal areas, SGD has, until recently, received much less attention. It is clear that SGD constitutes a major nutrient flux to coastal waters, with implications for estuarine ecology, eutrophication, and loss of coral reefs; however, fundamental questions regarding SGD remain unanswered: What are the spatial and temporal distributions of SGD offshore? How do seasonal and storm-related variations in aquifer recharge affect SGD flux and nutrient loading? What controls do aquifer structure and heterogeneity impose? How are SGD and saline recirculation related? Geophysical methods can provide insights to help answer these questions and improve the understanding of this intriguing and environmentally relevant hydrologic pheno...

Origin and Extent of Fresh Paleowaters on the Atlantic Continental Shelf, USA

While the existence of relatively fresh groundwater sequestered within permeable, porous sediments beneath the Atlantic continental shelf of North and South America has been known for some time, these waters have never been assessed as a potential resource. This fresh water was likely emplaced during Pleistocene sea-level low stands when the shelf was exposed to meteoric recharge and by elevated recharge in areas overrun by the Laurentide ice sheet at high latitudes. To test this hypothesis, we present results from a high-resolution paleohydrologic model of groundwater flow, heat and solute transport, ice sheet loading, and sea level fluctuations for the continental shelf from New Jersey to Maine over the last 2 million years. Our analysis suggests that the presence of fresh to brackish water within shallow Miocene sands more than 100 km offshore of New Jersey was facilitated by discharge of submarine springs along Baltimore and Hudson Canyons where these shallow aquifers crop out. Recharge rates four times modern levels were computed for portions of New Englands continental shelf that were overrun by the Laurentide ice sheet during the last glacial maximum. We estimate the volume of emplaced Pleistocene continental shelf fresh water (less than 1 ppt) to be 1300 km(3) in New England. We also present estimates of continental shelf fresh water resources for the U.S. Atlantic eastern seaboard (10(4) km(3)) and passive margins globally (3 x 10(5) km(3)). The simulation results support the hypothesis that offshore fresh water is a potentially valuable, albeit nonrenewable resource for coastal megacities faced with growing water shortages.

Geoelectrical evidence of bicontinuum transport in groundwater

Bicontinuum models and rate-limited mass transfer (RLMT) explain complex transport behavior (e.g., long tailing and rebound) in heterogeneous geologic media, but experimental verification is problematic because geochemical samples represent the mobile component of the pore space. Here, we present geophysical evidence of RLMT at the field scale during an aquifer-storage and recovery experiment in a fractured limestone aquifer in Charleston, South Carolina. We observe a hysteretic relation between measurements of pore-fluid conductivity and bulk electrical conductivity; this hysteresis contradicts advective-dispersive transport and the standard petrophysical model relating porefluid and bulk conductivity, but can be explained by considering bicontinuum transport models that include first-order RLMT. Using a simple numerical model, we demonstrate that geoelectrical measurements are sensitive to bicontinuum transport and RLMT parameters, which are otherwise difficult to infer from direct, hydrologic measurements.

A Small-Diameter NMR Logging Tool for Groundwater Investigations

A small-diameter nuclear magnetic resonance (NMR) logging tool has been developed and field tested at various sites in the United States and Australia. A novel design approach has produced relatively inexpensive, small-diameter probes that can be run in open or PVC-cased boreholes as small as 2 inches in diameter. The complete system, including surface electronics and various downhole probes, has been successfully tested in small-diameter monitoring wells in a range of hydrogeological settings. A variant of the probe that can be deployed by a direct-push machine has also been developed and tested in the field. The new NMR logging tool provides reliable, direct, and high-resolution information that is of importance for groundwater studies. Specifically, the technology provides direct measurement of total water content (total porosity in the saturated zone or moisture content in the unsaturated zone), and estimates of relative pore-size distribution (bound vs. mobile water content) and hydraulic conductivity. The NMR measurements show good agreement with ancillary data from lithologic logs, geophysical logs, and hydrogeologic measurements, and provide valuable information for groundwater investigations.

Understanding Water Column and Streambed Thermal Refugia for Endangered Mussels in the Delaware River

Groundwater discharge locations along the upper Delaware River, both discrete bank seeps and diffuse streambed upwelling, may create thermal niche environments that benefit the endangered dwarf wedgemussel (Alasmidonta heterodon). We seek to identify whether discrete or diffuse groundwater inflow is the dominant control on refugia. Numerous springs and seeps were identified at all locations where dwarf wedgemussels still can be found. Infrared imagery and custom high spatial resolution fiber-optic distributed temperature sensors reveal complex thermal dynamics at one of the seeps with a relatively stable, cold groundwater plume extending along the streambed/water-column interface during midsummer. This plume, primarily fed by a discrete bank seep, was shown through analytical and numerical heat-transport modeling to dominate temperature dynamics in the region of potential habitation by the adult dwarf wedgemussel.

Attenuation-difference radar tomography: results of a multiplane experiment at the U.S. Geological Survey Fractured-Rock Research Site, Mirror Lake, New Hampshire

Attenuation-difference, borehole-radar tomography was used to monitor a series of sodium chloride tracer injection tests conducted within the FSE wellfield at the U.S. Geological Survey Fractured-Rock Hydrology Research Site in Grafton County, New Hampshire, USA. Borehole-radar tomography surveys were conducted by using the sequential-scanning and injection method in three boreholes that form a triangular prism of adjoining tomographic image planes. Attenuation-difference data were inverted by using a weighted damped least-squares (WDLS) inversion method and several different solution simplicity schemes to suppress tomogram artifacts induced by the acquisition geometry, the location and magnitude of the anomaly, and noisy data. Qualitatively, flat tomograms generated by minimizing the norm of the first spatial derivative were most effective at suppressing artifacts. Although artifact suppression measures can be somewhat effective, negative consequences or artifact suppression include (1) reduction in the magnitude of the attenuation differences in the vicinity of the target anomaly and (2) blurring of the anomalies. These effects distort estimates of the location and magnitude of attenuation anomalies. In order to estimate robustly the location and magnitude of attenuation differences, areal constraints were imposed on the WDLS inversions to confine changes in attenuation to regions that are intersected by rays with large attenuation differences and bounded by rays with insignificant attenuation differences. Results from forward modeling support the application of these constraints. The resolution matrix was used to model the effects of acquisition geometry, target anomaly shape, location, and magnitude. For this study, the method of forward modeling indicates that estimates of pixel attenuation are improved by applying areal constraints to the WDLS inversions.

Dual‐domain mass‐transfer parameters from electrical hysteresis: Theory and analytical approach applied to laboratory, synthetic streambed, and groundwater experiments

Models of dual-domain mass transfer (DDMT) are used to explain anomalous aquifer transport behavior such as the slow release of contamination and solute tracer tailing. Traditional tracer experiments to characterize DDMT are performed at the flow path scale (meters), which inherently incorporates heterogeneous exchange processes; hence, estimated “effective” parameters are sensitive to experimental design (i.e., duration and injection velocity). Recently, electrical geophysical methods have been used to aid in the inference of DDMT parameters because, unlike traditional fluid sampling, electrical methods can directly sense less-mobile solute dynamics and can target specific points along subsurface flow paths. Here we propose an analytical framework for graphical parameter inference based on a simple petrophysical model explaining the hysteretic relation between measurements of bulk and fluid conductivity arising in the presence of DDMT at the local scale. Analysis is graphical and involves visual inspection of hysteresis patterns to (1) determine the size of paired mobile and less-mobile porosities and (2) identify the exchange rate coefficient through simple curve fitting. We demonstrate the approach using laboratory column experimental data, synthetic streambed experimental data, and field tracer-test data. Results from the analytical approach compare favorably with results from calibration of numerical models and also independent measurements of mobile and less-mobile porosity. We show that localized electrical hysteresis patterns resulting from diffusive exchange are independent of injection velocity, indicating that repeatable parameters can be extracted under varied experimental designs, and these parameters represent the true intrinsic properties of specific volumes of porous media of aquifers and hyporheic zones.