Geology topics

W. E. Sanford

Publications and source records attributed to W. E. Sanford.

At least 19 recordsLinked to original sources

USGS2013

Characterization and conceptualization of groundwater flow systems

This chapter discusses some of the fundamental concepts, data needs and approaches that aid in developing a general understanding of a groundwater system. Principles of the hydrological cycle are reviewed; the processes of recharge and discharge in aquifer systems; types of geological, hydrological and hydraulic data needed to describe the hydrogeological framework of an aquifer system; factors affecting the distribution of recharge to aquifers; and uses of groundwater chemistry, geochemical modelling, environmental tracers and age interpretations in groundwater studies. Together, these concepts and observations aid in developing a conceptualization of groundwater flow systems and provide input to the development of numerical models of a flow system. Conceptualization of the geology, hydrology, geochemistry, and hydrogeological and hydrochemical framework can be quite useful in planning, study design, guiding sampling campaigns, acquisition of new data and, ultimately, developing numerical models capable of assessing a wide variety of societal issues — for example, sustainability of groundwater resources in response to real or planned withdrawals from the system, CO 2 sequestration or other waste isolation issues (such as nuclear waste disposal).

Book chapter↗

USGS2013

Defining groundwater age

This book investigates applications of selected chemical and isotopic substances that can be used to recognize and interpret age information pertaining to ‘old’ groundwater (defined as water that was recharged on a timescale from approximately 1000 to more than 1 000 000 a). However, as discussed below, only estimates of the ‘age’ of water extracted from wells can be inferred. These groundwater age estimates are interpreted from measured concentrations of chemical and isotopic substances in the groundwater. Even then, there are many complicating factors, as discussed in this book. In spite of these limitations, much can be learned about the physics of groundwater flow and about the temporal aspects of groundwater systems from age interpretations of measured concentrations of environmental tracers in groundwater systems. This chapter puts the concept of ‘age’ into context, including its meaning and interpretation, and attempts to provide a unifying usage for the rest of the book.

Book chapter↗

USGS2012

Assessment of a ground water flow model of the Bangkok Basin, Thailand, using carbon-14-based ages and paleohydrology

A study was undertaken to understand the groundwater flow conditions in the Bangkok Basin, Thailand, by comparing 14 C-based and simulated groundwater ages. 14 C measurements were made on about 50 water samples taken from wells throughout the basin. Simulated ages were obtained using 1) backward-pathline tracking based on the well locations, and 2) results from a three-dimensional groundwater flow model. Comparisons of ages at these locations reveal a large difference between 14 C-based ages and ages predicted by the steady-state groundwater flow model. Mainly, 14 C and 13 C analyses indicate that groundwater in the Bangkok area is about 20,000 years old, whereas steady-state flow and transport simulations imply that groundwater in the Bangkok area is 50,000–100,000 years old. One potential reason for the discrepancy between simulated and 14 C-based ages is the assumption in the model of steady-state flow. Groundwater velocities were probably greater in the region before about 10,000 years ago, during the last glacial maximum, because of the lower position of sea level and the absence of the surficial Bangkok Clay. Paleoflow conditions were estimated and then incorporated into a second set of simulations. The new assumption was that current steady-state flow conditions existed for the last 8,000 years but were preceded by steady-state conditions representative of flow during the last glacial maximum. This “transient” paleohydrologic simulation yielded a mean simulated age that more closely agrees with the mean 14 C-based age, especially if the 14 C-based age corrected for diffusion into clay layers. Although the uncertainties in both the simulated and 14 C-based ages are nontrivial, the magnitude of the improved match in the mean age using a paleohydrologic simulation instead of a steady-state simulation suggests that flow conditions in the basin have changed significantly over the last 10,000–20,000 years. Given that the valid age range of 14 C-dating methods and the timing of the last glacial maximum are of similar magnitude, adjustments for paleohydrologic conditions may be required for many such studies.

Hydrogeology Journal↗

USGS2010

Inferring shallow groundwater flow in saprolite and fractured rock using environmental tracers

The Ridge and Valley Province of eastern Tennessee is characterized by (1) substantial topographic relief, (2) folded and highly fractured rocks of various lithologies that have low primary permeability and porosity, and (3) a shallow residuum of medium permeability and high total porosity. Conceptual models of shallow groundwater flow and solute transport in this system have been developed but are difficult to evaluate using physical characterization or short‐term tracer methods due to extreme spatial variability in hydraulic properties. In this paper we describe how chlorofluorocarbon 12, 3 H, and 3 He were used to infer groundwater flow and solute transport in saprolite and fractured rock near Oak Ridge, Tennessee. In the shallow residuum, fracture spacings are <0.05 m, suggesting that concentrations of these tracers in fractures and in the matrix have time to diffusionally equilibrate. The relatively smooth nature of tracer concentrations with depth in the residuum is consistent with this model and quantitatively suggests recharge fluxes of 0.2 to 0.4 m yr −1 . In contrast, groundwater flow within the unweathered rock appears to be controlled by fractures with spacings of the order of 2 to 5 m, and diffusional equilibration of fractures and matrix has not occurred. For this reason, vertical fluid fluxes in the unweathered rock cannot be estimated from the tracer data.

Water Resources Research↗

USGS2009

Pore-water chemistry from the ICDP-USGS coer hole in the Chesapeake Bay impact structure--Implications for paleohydrology, microbial habitat, and water resources

We investigated the groundwater system of the Chesapeake Bay impact structure by analyzing the pore-water chemistry in cores taken from a 1766-m-deep drill hole 10 km north of Cape Charles, Virginia. Pore water was extracted using high-speed centrifuges from over 100 cores sampled from a 1300 m section of the drill hole. The pore-water samples were analyzed for major cations and anions, stable isotopes of water and sulfate, dissolved and total carbon, and bioavailable iron. The results reveal a broad transition between fresh and saline water from 100 to 500 m depth in the post-impact sediment section, and an underlying syn-impact section that is almost entirely filled with brine. The presence of brine in the lowermost post-impact section and the trend in the dissolved chloride with depth suggest a transport process dominated by molecular diffusion and slow, compaction-driven, upward flow. Major ion results indicate residual effects of diagenesis from heating, and a pre-impact origin for the brine. High levels of dissolved organic carbon (6-95 mg/L) and the distribution of electron acceptors indicate an environment that may be favorable for microbial activity throughout the drilled section. The concentration and extent of the brine is much greater than had previously been observed, suggesting its occurrence may be common in the inner crater. However, groundwater flow conditions in the structure may reduce the salt-water-intrusion hazard associated with the brine.

Special Paper of the Geological Society of America↗

USGS2008

Temporal response of hydraulic head, temperature, and chloride concentrations to sea-level changes, Floridan aquifer system, USA

Three-dimensional density-dependent flow and transport modeling of the Floridan aquifer system, USA shows that current chloride concentrations are not in equilibrium with current sea level and, second, that the geometric configuration of the aquifer has a significant effect on system responses. The modeling shows that hydraulic head equilibrates first, followed by temperatures, and then by chloride concentrations. The model was constructed using a modified version of SUTRA capable of simulating multi-species heat and solute transport, and was compared to pre-development conditions using hydraulic heads, chloride concentrations, and temperatures from 315 observation wells. Three hypothetical, sinusoidal sea-level changes occurring over 100,000 years were used to evaluate how the simulated aquifer responds to sea-level changes. Model results show that hydraulic head responses lag behind sea-level changes only where the Miocene Hawthorn confining unit is thick and represents a significant restriction to flow. Temperatures equilibrate quickly except where the Hawthorn confining unit is thick and the duration of the sea-level event is long (exceeding 30,000 years). Response times for chloride concentrations to equilibrate are shortest near the coastline and where the aquifer is unconfined; in contrast, chloride concentrations do not change significantly over the 100,000-year simulation period where the Hawthorn confining unit is thick.

Hydrogeology Journal↗

USGS2008

Deep drilling into the Chesapeake Bay impact structure

Samples from a 1.76-kilometer-deep corehole drilled near the center of the late Eocene Chesapeake Bay impact structure (Virginia, USA) reveal its geologic, hydrologic, and biologic history. We conducted stratigraphic and petrologic analyses of the cores to elucidate the timing and results of impact-melt creation and distribution, transient-cavity collapse, and ocean-water resurge. Comparison of post-impact sedimentary sequences inside and outside the structure indicates that compaction of the crater fill influenced long-term sedimentation patterns in the mid-Atlantic region. Salty connate water of the target remains in the crater fill today, where it poses a potential threat to the regional groundwater resource. Observed depth variations in microbial abundance indicate a complex history of impact-related thermal sterilization and habitat modification, and subsequent post-impact repopulation.

Science↗

USGS2008

Subsurface microbial diversity in deep-granitic-fracture water in Colorado

A microbial community analysis using 16S rRNA gene sequencing was performed on borehole water and a granite rock core from Henderson Mine, a >1,000-meter-deep molybdenum mine near Empire, CO. Chemical analysis of borehole water at two separate depths (1,044 m and 1,004 m below the mine entrance) suggests that a sharp chemical gradient exists, likely from the mixing of two distinct subsurface fluids, one metal rich and one relatively dilute; this has created unique niches for microorganisms. The microbial community analyzed from filtered, oxic borehole water indicated an abundance of sequences from iron-oxidizing bacteria ( Gallionella spp.) and was compared to the community from the same borehole after 2 weeks of being plugged with an expandable packer. Statistical analyses with UniFrac revealed a significant shift in community structure following the addition of the packer. Phospholipid fatty acid (PLFA) analysis suggested that Nitrosomonadales dominated the oxic borehole, while PLFAs indicative of anaerobic bacteria were most abundant in the samples from the plugged borehole. Microbial sequences were represented primarily by Firmicutes , Proteobacteria , and a lineage of sequences which did not group with any identified bacterial division; phylogenetic analyses confirmed the presence of a novel candidate division. This “Henderson candidate division” dominated the clone libraries from the dilute anoxic fluids. Sequences obtained from the granitic rock core (1,740 m below the surface) were represented by the divisions Proteobacteria (primarily the family Ralstoniaceae ) and Firmicutes . Sequences grouping within Ralstoniaceae were also found in the clone libraries from metal-rich fluids yet were absent in more dilute fluids. Lineage-specific comparisons, combined with phylogenetic statistical analyses, show that geochemical variance has an important effect on microbial community structure in deep, subsurface systems.

Applied and Environmental Microbiology↗

USGS2008

Rapid response of a hydrologic system to volcanic activity: Masaya volcano, Nicaragua

Hydrologic systems change in response to volcanic activity, and in turn may be sensitive indicators of volcanic activity. Here we investigate the coupled nature of magmatic and hydrologic systems using continuous multichannel time series of soil temperature collected on the flanks of Masaya volcano, Nicaragua, one of the most active volcanoes in Central America. The soil temperatures were measured in a low-temperature fumarole field located 3.5 km down the flanks of the volcano. Analysis of these time series reveals that they respond extremely rapidly, on a time scale of minutes, to changes in volcanic activity also manifested at the summit vent. These rapid temperature changes are caused by increased flow of water vapor through flank fumaroles during volcanism. The soil temperature response, ~5 °C, is repetitive and complex, with as many as 13 pulses during a single volcanic episode. Analysis of the frequency spectrum of these temperature time series shows that these anomalies are characterized by broad frequency content during volcanic activity. They are thus easily distinguished from seasonal trends, diurnal variations, or individual rainfall events, which triggered rapid transient increases in temperature during 5% of events. We suggest that the mechanism responsible for the distinctive temperature signals is rapid change in pore pressure in response to magmatism, a response that can be enhanced by meteoric water infiltration. Monitoring of distal fumaroles can therefore provide insight into coupled volcanic-hydrologic-meteorologic systems, and has potential as an inexpensive monitoring tool.

Masaya Volcano↗

USGS2007

Atmospheric bromine flux from the coastal Abu Dhabi sabkhat: A ground-water mass-balance investigation

A solute mass-balance study of ground water of the 3000 km2 coastal sabkhat (salt flats) of the Emirate of Abu Dhabi, United Arab Emirates, documents an annual bromide loss of approximately 255 metric tons (0.0032 Gmoles), or 85 kg/km2. This value is an order of magnitude greater than previously published direct measurements from the atmosphere over an evaporative environment of a salar in Bolivia. Laboratory evidence, consistent with published reports, suggests that this loss is by vapor transport to the atmosphere. If this bromine flux to the atmosphere is representative of the total earth area of active salt flats then it is a significant, and generally under recognized, input to the global atmospheric bromide flux.

Geophysical Research Letters↗

USGS2007

Investigation of the groundwater system at Masaya Caldera, Nicaragua, using transient electromagnetics and numerical simulation

The distribution of groundwater beneath Masaya Volcano, in Nicaragua, and its surrounding caldera was characterized using the transient electromagnetic method (TEM). Multiple soundings were conducted at 30 sites. Models of the TEM data consistently indicate a resistive layer that is underlain by one or more conductive layers. These two layers represent the unsaturated and saturated zones, respectively, with the boundary between them indicating the water-table elevation. A map of the TEM data shows that the water table in the caldera is a subdued replica of the topography, with higher elevations beneath the edifice in the south-central caldera and lower elevations in the eastern caldera, coinciding with the elevation of Laguna de Masaya. These TEM data, combined with regional hydrologic data, indicate that the caldera in hydrologically isolated from the surrounding region, with as much as 60??m of difference in elevation of the groundwater table across caldera-bounding faults. The water-table information and estimates of fluxes of water through the system were used to constrain a numerical simulation of groundwater flow. The simulation results indicate that basalt flows in the outer parts of the caldera have a relatively high transmissivity, whereas the central edifice has a substantially lower transmissivity. A layer of relatively high transmissivity must be present at depth within the edifice in order to deliver the observed flux of water and steam to the active vent. This hydrologic information about the caldera provides a baseline for assessing the response of this isolated groundwater system to future changes in magmatic activity. ?? 2007.

Journal of Volcanology and Geothermal Research↗

USGS2007

A simulation of groundwater discharge and nitrate delivery to chesapeake bay from the lowermost delmarva peninsula, USA

A groundwater model has been developed for the lowermost Delmarva Peninsula, USA, that simulates saltwater intrusion into local confined aquifers and nitrate delivery to the Chesapeake Bay from the surficial aquifer. A flow path and groundwater-age analysis was performed using the model to estimate the timing of nitrate delivery to the bay over the next several decades. The simulated mean and median residence times of groundwater in the lowermost peninsula are 30 and 15 years, respectively. Current and future nitrate concentrations in coastal groundwater discharge were simulated based on local well data that include nitrate concentrations and groundwater age. A simulated future-trends analysis indicates that nitrate that has been applied to agricultural regions over the last few decades will continue to discharge into the bay for several decades to come. This study highlights the importance of considering the groundwater lag time that affects the mean transport time from diffuse contamination sources.

Conference Paper↗

USGS2006

Research opportunities in interdisciplinary ground-water science in the U.S. Geological Survey

This report is written for the scientifically literate reader but is not limited to those who are involved in ground-water science. The report is intended to encourage U.S. Geological Survey scientists to develop a sense of excitement about ground-water science in the agency, to inform scientists about existing and potential ground-water science opportunities, and to engage scientists and managers in interdisciplinary discussions and collaboration. The report is intended for use by U.S. Geological Survey and Department of the Interior management to formulate long-term ground-water science programs and to continue sustained support of ground-water monitoring and research, some of which may not have an immediate impact. Finally, the report can be used to communicate the U.S. Geological Survey's vision of ground-water science to Congress, partners, other agencies, and the research community at large with the goals of enhancing collaborative opportunities, sharing information, and maintaining dialogue regarding the directions of U.S. Geological Survey ground-water science.

Circular↗

USGS2005

Correcting for diffusion in carbon-14 dating of ground water

It has generally been recognized that molecular diffusion can be a significant process affecting the transport of carbon-14 in the subsurface when occurring either from a permeable aquifer into a confining layer or from a fracture into a rock matrix. An analytical solution that is valid for steady-state radionuclide transport through fractured rock is shown to be applicable to many multilayered aquifer systems. By plotting the ratio of the rate of diffusion to the rate of decay of carbon-14 over the length scales representative of several common hydrogeologic settings, it is demonstrated that diffusion of carbon-14 should often be not only a significant process, but a dominant one relative to decay. An age-correction formula is developed and applied to the Bangkok Basin of Thailand, where a mean carbon-14-based age of 21,000 years was adjusted to 11,000 years to account for diffusion. This formula and its graphical representation should prove useful for many studies, for they can be used first to estimate the potential role of diffusion and then to make a simple first-order age correction if necessary.

Groundwater↗

USGS2005

A simulation of the hydrothermal response to the Chesapeake Bay bolide impact

Groundwater more saline than seawater has been discovered in the tsunami breccia of the Chesapeake Bay impact Crater. One hypothesis for the origin of this brine is that it may be a liquid residual following steam separation in a hydrothermal system that evolved following the impact. Initial scoping calculations have demonstrated that it is feasible such a residual brine could have remained in the crater for the 35 million years since impact. Numerical simulations have been conducted using the code HYDROTHERM to test whether or not conditions were suitable in the millennia following the impact for the development of a steam phase in the hydrothermal system. Hydraulic and thermal parameters were estimated for the bedrock underlying the crater and the tsunami breccia that fills the crater. Simulations at three different breccia permeabilities suggest that the type of hydrothermal system that might have developed would have been very sensitive to the permeability. A relatively low breccia permeability (1 ?? 10-16 m2) results in a system partitioned into a shallow water phase and a deeper superheated steam phase. A moderate breccia permeability (1 ?? 10-15 m2 ) results in a system with regionally extensive multiphase conditions. A relatively high breccia permeability (1 ?? 10-14 m2 ) results in a system dominated by warm-water convection cells. The permeability of the crater breccia could have had any of these values at given depths and times during the hydrothermal system evolution as the sediments compacted. The simulations were not able to take into account transient permeability conditions, or equations of state that account for the salt content of seawater. Results suggest, however, that it is likely that steam conditions existed at some time in the system following impact, providing additional evidence that is consistent with a hydrothermal origin for the crater brine. ?? Blackwell Publishing Ltd.

Geofluids↗

USGS2005

Chemical openness and potential for misinterpretation of the solute environment of coastal sabkhat

Sabkha deposits in the geologic record are commonly used to interpret the environmental conditions of deposition. Implicit in this use is the assumption that the solute system is chemically closed, that is, the authigenic minerals represent the composition of the fluids in their environment of origin. Thermodynamic and mass-balance calculations based on measurements of water and solute flux of contemporary Abu Dhabi coastal sabkha system, however, demonstrate that the system is open for sodium and chloride, where nearly half of the input is lost, but closed for sulfur, where nearly 100% is retained. Sulfur and chloride isotopes were consistent with this observation. If these sabkha deposits were preserved in the geologic record, they would suggest a solute environment rich in sulfate and poor in chloride; yet the reverse is true. In most coastal-sabkha environments, capillary forces bring solutes and water to the surface, where the water evaporates and halite, carnallite, sylvite, and other soluble minerals are precipitated. Retrograde minerals, such as anhydrite, calcite, dolomite, and gypsum, however, precipitate and accumulate in the capillary zone beneath the surface of the coastal sabkha. Because they possess relatively low solubility and are below the surface, these retrograde minerals are protected from dissolution and physical erosion occurring from infrequent but intense rainfall events. Thus, they are more likely to be preserved in the geological record than highly soluble minerals formed on the surface. ?? 2004 Elsevier B.V. All rights reserved.

Chemical Geology↗

USGS2004

Hydrochemical tracers in the middle Rio Grande Basin, USA: 1. Conceptualization of groundwater flow

Chemical and isotopic data for groundwater from throughout the Middle Rio Grande Basin, central New Mexico, USA, were used to identify and map groundwater flow from 12 sources of water to the basin, evaluate radiocarbon ages, and refine the conceptual model of the Santa Fe Group aquifer system. Hydrochemical zones, representing groundwater flow over thousands to tens of thousands of years, can be traced over large distances through the primarily siliciclastic aquifer system. The locations of the hydrochemical zones mostly reflect the “modern” predevelopment hydraulic-head distribution, but are inconsistent with a trough in predevelopment water levels in the west-central part of the basin, indicating that this trough is a transient rather than a long-term feature of the aquifer system. Radiocarbon ages adjusted for geochemical reactions, mixing, and evapotranspiration/dilution processes in the aquifer system were nearly identical to the unadjusted radiocarbon ages, and ranged from modern to more than 30 ka. Age gradients from piezometer nests ranged from 0.1 to 2 year cm –1 and indicate a recharge rate of about 3 cm year –1 for recharge along the eastern mountain front and infiltration from the Rio Grande near Albuquerque. There has been appreciably less recharge along the eastern mountain front north and south of Albuquerque.

Rio Grande Basin↗

USGS2004

Hydrochemical tracers in the middle Rio Grande Basin, USA: 2. Calibration of a groundwater-flow model

The calibration of a groundwater model with the aid of hydrochemical data has demonstrated that low recharge rates in the Middle Rio Grande Basin may be responsible for a groundwater trough in the center of the basin and for a substantial amount of Rio Grande water in the regional flow system. Earlier models of the basin had difficulty reproducing these features without any hydrochemical data to constrain the rates and distribution of recharge. The objective of this study was to use the large quantity of available hydrochemical data to help calibrate the model parameters, including the recharge rates. The model was constructed using the US Geological Survey's software MODFLOW, MODPATH, and UCODE, and calibrated using 14C activities and the positions of certain flow zones defined by the hydrochemical data. Parameter estimation was performed using a combination of nonlinear regression techniques and a manual search for the minimum difference between field and simulated observations. The calibrated recharge values were substantially smaller than those used in previous models. Results from a 30,000-year transient simulation suggest that recharge was at a maximum about 20,000 years ago and at a minimum about 10,000 years ago.

Rio Grande Basin↗