Search USGS⌕ Search

Geology topics

C.I. Voss

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

At least 19 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↗

Airborne electromagnetic imaging of discontinuous permafrost

The evolution of permafrost in cold regions is inextricably connected to hydrogeologic processes, climate, and ecosystems. Permafrost thawing has been linked to changes in wetland and lake areas, alteration of the groundwater contribution to streamflow, carbon release, and increased fire frequency. But detailed knowledge about the dynamic state of permafrost in relation to surface and groundwater systems remains an enigma. Here, we present the results of a pioneering ∼1,800 line-kilometer airborne electromagnetic survey that shows sediments deposited over the past ∼4 million years and the configuration of permafrost to depths of ∼100 meters in the Yukon Flats area near Fort Yukon, Alaska. The Yukon Flats is near the boundary between continuous permafrost to the north and discontinuous permafrost to the south, making it an important location for examining permafrost dynamics. Our results not only provide a detailed snapshot of the present-day configuration of permafrost, but they also expose previously unseen details about potential surface – groundwater connections and the thermal legacy of surface water features that has been recorded in the permafrost over the past ∼1,000 years. This work will be a critical baseline for future permafrost studies aimed at exploring the connections between hydrogeologic, climatic, and ecological processes, and has significant implications for the stewardship of Arctic environments.

Alaska↗

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↗

Controls on groundwater flow in the Bengal Basin of India and Bangladesh: Regional modeling analysis

Groundwater for domestic and irrigation purposes is produced primarily from shallow parts of the Bengal Basin aquifer system (India and Bangladesh), which contains high concentrations of dissolved arsenic (exceeding worldwide drinking water standards), though deeper groundwater is generally low in arsenic. An essential first step for determining sustainable management of the deep groundwater resource is identification of hydrogeologic controls on flow and quantification of basin-scale groundwater flow patterns. Results from groundwater modeling, in which the Bengal Basin aquifer system is represented as a single aquifer with higher horizontal than vertical hydraulic conductivity, indicate that this anisotropy is the primary hydrogeologic control on the natural flowpath lengths. Despite extremely low hydraulic gradients due to minimal topographic relief, anisotropy implies large-scale (tens to hundreds of kilometers) flow at depth. Other hydrogeologic factors, including lateral and vertical changes in hydraulic conductivity, have minor effects on overall flow patterns. However, because natural hydraulic gradients are low, the impact of pumping on groundwater flow is overwhelming; modeling indicates that pumping has substantially changed the shallow groundwater budget and flowpaths from predevelopment conditions. ?? Springer-Verlag 2009.

Hydrogeology Journal↗

Comparison of alternative representations of hydraulic-conductivity anisotropy in folded fractured-sedimentary rock: Modeling groundwater flow in the Shenandoah Valley (USA)

A numerical representation that explicitly represents the generalized three-dimensional anisotropy of folded fractured-sedimentary rocks in a groundwater model best reproduces the salient features of the flow system in the Shenandoah Valley, USA. This conclusion results from a comparison of four alternative representations of anisotropy in which the hydraulic-conductivity tensor represents the bedrock structure as (model A) anisotropic with variable strikes and dips, (model B) horizontally anisotropic with a uniform strike, (model C) horizontally anisotropic with variable strikes, and (model D) isotropic. Simulations using the US Geological Survey groundwater flow and transport model SUTRA are based on a representation of hydraulic conductivity that conforms to bedding planes in a three-dimensional structural model of the valley that duplicates the pattern of folded sedimentary rocks. In the most general representation, (model A), the directions of maximum and medium hydraulic conductivity conform to the strike and dip of bedding, respectively, while the minimum hydraulic-conductivity direction is perpendicular to bedding. Model A produced a physically realistic flow system that reflects the underlying bedrock structure, with a flow field that is significantly different from those produced by the other three models.

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↗

Testing density-dependent groundwater models: Two-dimensional steady state unstable convection in infinite, finite and inclined porous layers

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA––A model for saturated–unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4 π 2 and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature.

Advances in Water Resources↗

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↗