USGS ScienceSearch

Geology topics

C.I. Voss

Publications and source records attributed to C.I. Voss.

17 recordsLinked to original sources

A simulation-based approach for designing effective field-sampling programs to evaluate contamination risk of groundwater supplies

An approach to model discrimination and network design for evaluation of groundwater contamination risk is proposed and demonstrated by application to a site in a glaciofluvial aquifer in Sweden. The approach consists of first hypothesizing alternative conceptual models of hydrogeology at the site on the basis of both quantitative data and qualitative information. The conceptual models are then expressed as two-dimensional numerical models of groundwater flow and solute transport, and model attributes controlling risk to the water supply are determined by simulation. Model predictions of response to a specific field test are made with each model that affects risk. Regions for effective measurement networks are then identified. Effective networks are those that capture sufficient information to determine which of the hypothesized models best describes the system with a minimum of measurement points. For the example site in Sweden, the network is designed such that important system parameters may be accurately estimated at the same time as model discrimination is carried out.

Hydrogeology Journal

Sensitivity analysis of lake mass balance in discontinuous permafrost: the example of disappearing Twelvemile Lake, Yukon Flats, Alaska (USA)

Many lakes in northern high latitudes have undergone substantial changes in surface area over the last four decades, possibly as a result of climate warming. In the discontinuous permafrost of Yukon Flats, interior Alaska (USA), these changes have been non-uniform across adjacent watersheds, suggesting local controls on lake water budgets. Mechanisms that could explain the decreasing mass of one lake in Yukon Flats since the early 1980s, Twelvemile Lake, are identified via a scoping analysis that considers plausible changes in snowmelt mass and infiltration, permafrost distribution, and climate warming. Because predicted changes in evaporation (2 cmyr -1 ) are inadequate to explain the observed 17.5 cmyr -1 reduction in mass balance, other mechanisms are required. The most important potential mechanisms are found to involve: (1) changes in shallow, lateral groundwater flow to the lake possibly facilitated by vertical freeze-thaw migration of the permafrost table in gravel; (2) increased loss of lake water as downward groundwater flow through an open talik to a permeable subpermafrost flowpath; and (3) reduced snow meltwater inputs due to decreased snowpack mass and increased infiltration of snowmelt into, and subsequent evaporation from, fine-grained sediment mantling the permafrost-free lake basin.

Alaska

Glaciation and regional groundwater flow in the Fennoscandian shield

Regional-scale groundwater flow modeling of the Fennoscandian shield suggests that groundwater flow can be strongly affected by future climate change and glaciation. We considered variable-density groundwater flow in a 1500-km-long and approximately 10-km-deep cross-section through southern Sweden. Groundwater flow and shield brine transport in the cross-sectional model were analyzed under projected surface conditions for the next 140 ka. Simulations suggest that blockage of recharge and discharge by low-permeability permafrost or cold-based ice causes sinking of brine and consequent freshening of near-surface water in areas of natural discharge. Although recharge of basal meltwater is limited by the requirement that water pressure at the base of the ice sheet not exceed the pressure exerted by the weight of the ice, warm-based ice with basal melting creates a potential for groundwater recharge rates much larger than those of present, ice-free conditions. In the simulations, regional-scale redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (approximately 10 ka). During glacial retreat, significant upward flow of groundwater may occur below the ice sheet owing to pressure release. If the mechanical loading efficiency of the rocks is high, both subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are reduced because of loading-induced pressure changes. The maximum rate of groundwater discharge in the simulations occurs at the receding ice margin, and some discharge occurs below incursive postglacial seas. Recharge of basal meltwater could decrease the concentration of dissolved solids significantly below present-day levels at depths of up to several kilometers and may bring oxygenated conditions to an otherwise reducing chemical environment for periods exceeding 10 ka.

Fennoscandian Shield

Laboratory investigation of longitudinal dispersion in anisotropic porous media

In this study, laboratory experiments were used to investigate mechanisms that may cause anisotropy in the dispersion coefficient and to investigate the relation between anisotropy in hydraulic conductivity and anisotropy in longitudinal dispersion. Measurements of sodium chloride concentration (used as a tracer) were made at 105 in situ sampling locations in a new type of sand box designed to allow flow in either of two perpendicular directions. Two types of hydraulic anisotropy were examined. The first consisted of structured zones of increased hydraulic conductivity within a lower-conductivity medium. The second type involved low-conductivity platelike inclusions within a homogeneous, isotropic medium. The plates were aligned such that the tortuosity was increased only in one principal direction of permeability. Results using two examples of the first type of media showed that the apparent longitudinal dispersivities for flow parallel to the high-conductivity direction were greater than those perpendicular to this direction. Two examples of the second type of media produced smaller apparent longitudinal dispersivities for flow parallel to the high-conductivity direction. The results suggest that the mechanisms causing dispersive anisotropy can be related, conceptually, to the factors causing hydraulic anisotropy.

Water Resources Research

Vulnerability of deep groundwater in the Bengal Aquifer System to contamination by arsenic

Shallow groundwater, the primary water source in the Bengal Basin, contains up to 100 times the World Health Organization (WHO) drinking-water guideline of 10g l 1 arsenic (As), threatening the health of 70 million people. Groundwater from a depth greater than 150m, which almost uniformly meets the WHO guideline, has become the preferred alternative source. The vulnerability of deep wells to contamination by As is governed by the geometry of induced groundwater flow paths and the geochemical conditions encountered between the shallow and deep regions of the aquifer. Stratification of flow separates deep groundwater from shallow sources of As in some areas. Oxidized sediments also protect deep groundwater through the ability of ferric oxyhydroxides to adsorb As. Basin-scale groundwater flow modelling suggests that, over large regions, deep hand-pumped wells for domestic supply may be secure against As invasion for hundreds of years. By contrast, widespread deep irrigation pumping might effectively eliminate deep groundwater as an As-free resource within decades. Finer-scale models, incorporating spatial heterogeneity, are needed to investigate the security of deep municipal abstraction at specific urban locations. ?? 2010 Macmillan Publishers Limited. All rights reserved.

Nature Geoscience

Estimation of regional-scale groundwater flow properties in the Bengal Basin of India and Bangladesh

Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, 14C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on 14C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 104. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system. ?? US Government 2009.

Hydrogeology Journal

Finite-difference grid for a doublet well in an anisotropic aquifer

The U.S. Geological Survey is modeling hydraulic flow and thermal-energy transport at a two-well injection/ withdrawal system in St. Paul, Minnesota. The design of the finite-difference model grid for the doublet-well system is complicated because the aquifer is anisotropic and the principal axes of transmissivity are not aligned with the axis between the two wells. An analytical solution for flow in a doublet-well system in an infinite anisotropic aquifer was employed in the design of a grid with artificial boundaries placed in the midst of the flow field. Flow-net analysis was used to determine water flux across an cquipotential boundary and to assign approximate flux values at model boundaries. This enabled the simulation of the effects of the entire flow field, although only a small part was modeled. The validity of the flux values at the model boundaries for the isothermal case was tested by simulation of an eight-day injection test of ambient-temperature water. Model-computed pressures compared very favorably with field-observed pressures. The validity of boundary-flux values also was tested for nonisothermal conditions by simulation of injection of 300o F water at 300 gallons per minute for eight days.

Ground Water

Inverse modeling for seawater intrusion in coastal aquifers: Insights about parameter sensitivities, variances, correlations and estimation procedures derived from the Henry problem

Inverse modeling studies employing data collected from the classic Henry seawater intrusion problem give insight into several important aspects of inverse modeling of seawater intrusion problems and effective measurement strategies for estimation of parameters for seawater intrusion. Despite the simplicity of the Henry problem, it embodies the behavior of a typical seawater intrusion situation in a single aquifer. Data collected from the numerical problem solution are employed without added noise in order to focus on the aspects of inverse modeling strategies dictated by the physics of variable-density flow and solute transport during seawater intrusion. Covariances of model parameters that can be estimated are strongly dependent on the physics. The insights gained from this type of analysis may be directly applied to field problems in the presence of data errors, using standard inverse modeling approaches to deal with uncertainty in data. Covariance analysis of the Henry problem indicates that in order to generally reduce variance of parameter estimates, the ideal places to measure pressure are as far away from the coast as possible, at any depth, and the ideal places to measure concentration are near the bottom of the aquifer between the center of the transition zone and its inland fringe. These observations are located in and near high-sensitivity regions of system parameters, which may be identified in a sensitivity analysis with respect to several parameters. However, both the form of error distribution in the observations and the observation weights impact the spatial sensitivity distributions, and different choices for error distributions or weights can result in significantly different regions of high sensitivity. Thus, in order to design effective sampling networks, the error form and weights must be carefully considered. For the Henry problem, permeability and freshwater inflow can be estimated with low estimation variance from only pressure or only concentration observations. Permeability, freshwater inflow, solute molecular diffusivity, and porosity can be estimated with roughly equivalent confidence using observations of only the logarithm of concentration. Furthermore, covariance analysis allows a logical reduction of the number of estimated parameters for ill-posed inverse seawater intrusion problems. Ill-posed problems may exhibit poor estimation convergence, have a non-unique solution, have multiple minima, or require excessive computational effort, and the condition often occurs when estimating too many or co-dependent parameters. For the Henry problem, such analysis allows selection of the two parameters that control system physics from among all possible system parameters.

Advances in Water Resources

Three-dimensional variable-density flow simulation of a coastal aquifer in southern Oahu, Hawaii, USA

Three-dimensional modeling of groundwater flow and solute transport in the Pearl Harbor aquifer, southern Oahu, Hawaii, shows that the readjustment of the freshwater–saltwater transition zone takes a long time following changes in pumping, irrigation, or recharge in the aquifer system. It takes about 50 years for the transition zone to move 90% of the distance to its new steady position. Further, the Ghyben–Herzberg estimate of the freshwater/saltwater interface depth occurred between the 10 and 50% simulated seawater concentration contours in a complex manner during 100 years of the pumping history of the aquifer. Thus, it is not a good predictor of the depth of potable water. Pre-development recharge was used to simulate the 1880 freshwater-lens configuration. Historical pumpage and recharge distributions were used and the resulting freshwater-lens size and position were simulated through 1980. Simulations show that the transition zone moved upward and landward during the period simulated. Previous groundwater flow models for Oahu have been limited to areal models that simulate a sharp interface between freshwater and saltwater or solute-transport models that simulate a vertical aquifer section. The present model is based on the US Geological Survey’s three-dimensional solute transport (3D SUTRA) computer code. Using several new tools for pre- and post-processing of model input and results have allowed easy model construction and unprecedented visualization of the freshwater lens and underlying transition zone in Hawaii’s most developed aquifer.

Hawaii

Conceptual uncertainty in crystalline bedrock: Is simple evaluation the only practical approach?

A simple evaluation can be used to characterize the capacity of crystalline bedrock to act as a barrier to release radionuclides from a nuclear waste repository. Physically plausible bounds on groundwater flow and an effective transport-resistance parameter are estimated based on fundamental principles and idealized models of pore geometry. Application to an intensively characterized site in Sweden shows that, due to high spatial variability and uncertainty regarding properties of transport paths, the uncertainty associated with the geological barrier is too high to allow meaningful discrimination between good and poor performance. Application of more complex (stochastic-continuum and discrete-fracture-network) models does not yield a significant improvement in the resolution of geological barrier performance. Comparison with seven other less intensively characterized crystalline study sites in Sweden leads to similar results, raising a question as to what extent the geological barrier function can be characterized by state-of-the art site investigation methods prior to repository construction. A simple evaluation provides a simple and robust practical approach for inclusion in performance assessment.

IAHS-AISH Publication

Conceptual uncertainty in crystalline bedrock: Is simple evaluation the only practical approach?

A simple evaluation can be used to characterise the capacity of crystalline bedrock to act as a barrier to releases of radionuclides from a nuclear waste repository. Physically plausible bounds on groundwater flow and an effective transport-resistance parameter are estimated based on fundamental principles and idealised models of pore geometry. Application to an intensively characterised site in Sweden shows that, due to high spatial variability and uncertainty regarding properties of transport paths, the uncertainty associated with the geological barrier is too high to allow meaningful discrimination between good and poor performance. Application of more complex (stochastic-continuum and discrete-fracture-network) models does not yield a significant improvement in the resolution of geologic-barrier performance. Comparison with seven other less intensively characterised crystalline study sites in Sweden leads to similar results, raising a question as to what extent the geological barrier function can be characterised by state-of-the art site investigation methods prior to repository construction. A simple evaluation provides a simple and robust practical approach for inclusion in performance assessment.

Acta Universitatis Carolinae, Geologica

Factors controlling the configuration of the fresh-saline water interface in the Dead Sea coastal aquifers: Synthesis of TDEM surveys and numerical groundwater modeling

TDEM (time domain electromagnetic) traverses in the Dead Sea (DS) coastal aquifer help to delineate the configuration of the interrelated fresh-water and brine bodies and the interface in between. A good linear correlation exists between the logarithm of TDEM resistivity and the chloride concentration of groundwater, mostly in the higher salinity range, close to that of the DS brine. In this range, salinity is the most important factor controlling resistivity. The configuration of the fresh-saline water interface is dictated by the hydraulic gradient, which is controlled by a number of hydrological factors. Three types of irregularities in the configuration of fresh-water and saline-water bodies were observed in the study area: 1. Fresh-water aquifers underlying more saline ones ("Reversal") in a multi-aquifer system. 2. "Reversal" and irregular residual saline-water bodies related to historical, frequently fluctuating DS base level and respective interfaces, which have not undergone complete flushing. A rough estimate of flushing rates may be obtained based on knowledge of the above fluctuations. The occurrence of salt beds is also a factor affecting the interface configuration. 3. The interface steepens towards and adjacent to the DS Rift fault zone. Simulation analysis with a numerical, variable-density flow model, using the US Geological Survey's SUTRA code, indicates that interface steep- ening may result from a steep water-level gradient across the zone, possibly due to a low hydraulic conductivity in the immediate vicinity of the fault.

Hydrogeology Journal

Effects of clay dispersion on aquifer storage and recovery in coastal aquifers

Cyclic injection, storage, and withdrawal of freshwater in brackish aquifers is a form of aquifer storage and recovery (ASR) that can beneficially supplement water supplies in coastal areas. A 1970s field experiment in Norfolk, Virginia, showed that clay dispersion in the unconsolidated sedimentary aquifer occurred because of cation exchange on clay minerals as freshwater displaced brackish formation water. Migration of interstitial clay particles clogged pores, reduced permeability, and decreased recovery efficiency, but a calcium preflush was found to reduce clay dispersion and lead to a higher recovery efficiency. Column experiments were performed in this study to quantify the relations between permeability changes and clay mineralogy, clay content, and initial water salinity. The results of these experiments indicate that dispersion of montmorillonite clay is a primary contributor to formation damage. The reduction in permeability by clay dispersion may be expressed as a linear function of chloride content. Incorporating these simple functions into a radial, cross-sectional, variable-density, ground-water flow and transport model yielded a satisfactory simulation of the Norfolk field test - and represented an improvement over the model that ignored changes in permeability. This type of model offers a useful planning and design tool for ASR operations in coastal clastic aquifer systems.

Transport in Porous Media

Analysis of an anisotropic coastal aquifer system using variable-density flow and solute transport simulation

The groundwater system in southern Oahu, Hawaii consists of a thick, areally extensive freshwater lens overlying a zone of transition to a thick saltwater body. This system is analyzed in cross section with a variable-density groundwater flow and solute transport model on a regional scale. The simulation is difficult, because the coastal aquifer system has a saltwater transition zone that is broadly dispersed near the discharge area, but is very sharply defined inland. Steady-state simulation analysis of the transition zone in the layered basalt aquifer of southern Oahu indicates that a small transverse dispersivity is characteristic of horizontal regional flow. Further, in this system flow is generally parallel to isochlors and steady-state behavior is insensitive to the longitudinal dispersivity. Parameter analysis identifies that only six parameters control the complex hydraulics of the system: horizontal and vertical hydraulic conductivity of the basalt aquifer; hydraulic conductivity of the confining "caprock" layer; leakance below the caprock; specific yield; and aquifer matrix compressibility. The best-fitting models indicate the horizontal hydraulic conductivity is significantly greater than the vertical hydraulic conductivity. These models give values for specific yield and aquifer compressibility which imply a considerable degree of compressive storage in the water table aquifer.

Hawaii

Numerical simulation of advective-dispersive multisolute transport with sorption, ion exchange and equilibrium chemistry

A model was developed that can simulate the effect of certain chemical and sorption reactions simultaneously among solutes involved in advective-dispersive transport through porous media. The model is based on a methodology that utilizes physical-chemical relationships in the development of the basic solute mass-balance equations; however, the form of these equations allows their solution to be obtained by methods that do not depend on the chemical processes. The chemical environment is governed by the condition of local chemical equilibrium, and may be defined either by the linear sorption of a single species and two soluble complexation reactions which also involve that species, or binary ion exchange and one complexation reaction involving a common ion. Partial differential equations that describe solute mass balance entirely in the liquid phase are developed for each tenad (a chemical entity whose total mass is independent of the reaction process) in terms of their total dissolved concentration. These equations are solved numerically in two dimensions through the modification of an existing groundwater flow/transport computer code. (Author 's abstract)

Water-Resources Investigations Report

AQUIFEM-SALT; a finite-element model for aquifers containing a seawater interface

Described are modifications to AQUIFEM, a finite element areal ground-water flow model for aquifer evaluation. The modified model, AQUIFEM-SALT, simulates an aquifer containing a freshwater body that freely floats on seawater. Parts of the freshwater lens may be confined above and below by less permeable units. Theory, code modifications, and model verification are discussed. A modified input data list is included. This report is intended as a companion to the original AQUIFEM documentation. (USGS)

Water-Resources Investigations Report

A finite-element simulation model for saturated-unsaturated, fluid-density-dependent ground-water flow with energy transport or chemically- reactive single-species solute transport

SUTRA (Saturated-Unsaturated Transport) is a computer program which can be used to simulate the movement of fluid and the transport of either energy or dissolved substances in a subsurface environment. The model employs a two-dimensional hybrid finite-element and integrated-finite-difference method to approximate the governing equations that describe the two interdependent processes that are simulated by SUTRA: (1) fluid density-dependent saturated or unsaturated groundwater flow, and either (2a) transport of a solute in the groundwater, in which the solute may be subject to: equilibrium adsorption on the porous matrix, and both first-order and zero-order production or decay, or, (2b) transport of thermal energy in the groundwater and solid matrix of the aquifer. SUTRA provides, as the primary calculated results, fluid pressures and either solute concentrations or temperatures, as they vary with time, everywhere in the simulated subsurface system. SUTRA may also be used to simulate simpler subsets of the above process. SUTRA may be employed for areal and cross-sectional models of saturated groundwater flow systems, and for cross-sectional models of unsaturated zone flow. Solute transport simulation using SUTRA may be used to simulate natural or man-induced chemical transport, solute sorption, production and decay. SUTRA may be used for simulation of variable density leachate movement, and for cross-sectional simulation of salt-water intrusion in aquifers at near-well or regional scales, with either dispersed or relatively sharp transition zones between fresh water and salt water. SUTRA energy transport simulation may be employed to model thermal regimes in aquifers, subsurface heat conduction, aquifer thermal energy storage systems, geothermal reservoirs, thermal pollution of aquifers, and natural hydrogeologic convection systems. (USGS)

Water-Resources Investigations Report