USGS ScienceSearch

Geology topics

Richard W. Healy

Publications and source records attributed to Richard W. Healy.

30 records · Page 2Linked to original sources

Modeling methods

Simulation models are widely used in all types of hydrologic studies, and many of these models can be used to estimate recharge. Models can provide important insight into the functioning of hydrologic systems by identifying factors that influence recharge. The predictive capability of models can be used to evaluate how changes in climate, water use, land use, and other factors may affect recharge rates. Most hydrological simulation models, including watershed models and groundwater-flow models, are based on some form of water-budget equation, so the material in this chapter is closely linked to that in Chapter 2. Empirical models that are not based on a water-budget equation have also been used for estimating recharge; these models generally take the form of simple estimation equations that define annual recharge as a function of precipitation and possibly other climatic data or watershed characteristics. Model complexity varies greatly. Some models are simple accounting models; others attempt to accurately represent the physics of water movement through each compartment of the hydrologic system. Some models provide estimates of recharge explicitly; for example, a model based on the Richards equation can simulate water movement from the soil surface through the unsaturated zone to the water table. Recharge estimates can be obtained indirectly from other models. For example, recharge is a parameter in groundwater-flow models that solve for hydraulic head (i.e. groundwater level). Recharge estimates can be obtained through a model calibration process in which recharge and other model parameter values are adjusted so that simulated water levels agree with measured water levels. The simulation that provides the closest agreement is called the best fit, and the recharge value used in that simulation is the model-generated estimate of recharge.

Book chapter

Methods based on surface-water data

Streamflow data are commonly used to estimate recharge rates in humid and subhumid regions, in part because of the abundance of streamflow data and the availability of computer programs for analyzing those data. Most of the methods described in this chapter are easy to use, but application of any of the methods should be accompanied by a careful analysis of the underlying assumptions. The methods estimate exchange rates between groundwater and surface-water bodies. That exchange can represent focused recharge from a losing stream, or, as in the case of groundwater discharge to a stream, the exchange can reflect diffuse recharge that occurs over widespread areas. Some of these methods may be unfamiliar to groundwater hydrologists because they were not developed specifically for the study of groundwater recharge; instead, they were developed for purposes such as sizing of culverts and bridge openings, predicting low-flow rates in streams, or developing an understanding of stream-water quality and the ability of a stream to assimilate solutes and contaminants. The fact that base-flow or recharge estimates are generated as byproducts of these methods does not diminish the usefulness or applicability of the methods in recharge studies. Techniques presented herein include the stream water-budget method, seepage meters, Darcy methods, streamflow duration curves, traditional streamflow hydrograph analyses (including hydrograph separation and recession-curve displacement), and chemical and isotopic hydrograph separation techniques. Some of these methods are designed specifically for estimating focused recharge; others are for estimating diffuse recharge. Discussions are centered on groundwater movement to or from streams, but the principles discussed and the methods described are equally applicable for groundwater exchange with other surface-water bodies, such as lakes, reservoirs, and wetlands. Proper application of any method requires a good conceptual model of the hydrologic system and a solid understanding of underlying assumptions. Prior to presentation of individual methods, background discussions are given on the exchange of groundwater and surface water and on the relationship between base flow and recharge. These discussions illustrate assumptions inherent to the methods and provide some guidelines for assessing the validity of those assumptions.

Book chapter

Heat tracer methods

The flow of heat in the subsurface is closely linked to the movement of water (Ingebritsen et al ., 2006). As such, heat has been used as a tracer in groundwater studies for more than 100 years (Anderson, 2005). As with chemical and isotopic tracers (Chapter 7), spatial or temporal trends in surface and subsurface temperatures can be used to infer rates of water movement. Temperature can be measured accurately, economically, at high frequencies, and without the need to obtain water samples, facts that make heat an attractive tracer. Temperature measurements made over space and time can be used to infer rates of recharge from a stream or other surface water body (Lapham, 1989; Stonestrom and Constantz, 2003); measurements can also be used to estimate rates of steady drainage through depth intervals within thick unsaturated zones (Constantz et al ., 2003; Shan and Bodvarsson, 2004). Several thorough reviews of heat as a tracer in hydrologic studies have recently been published (Constantz et al ., 2003; Stonestrom and Constantz, 2003; Anderson, 2005; Blasch et al ., 2007; Constantz et al ., 2008). This chapter summarizes heat-tracer approaches that have been used to estimate recharge. Some clarification in terminology is presented here to avoid confusion in descriptions of the various approaches that follow. Diffuse recharge is that which occurs more or less uniformly across large areas in response to precipitation, infiltration, and drainage through the unsaturated zone. Estimates of diffuse recharge determined using measured temperatures in the unsaturated zone are referred to as potential recharge because it is possible that not all of the water moving through the unsaturated zone will recharge the aquifer; some may be lost to the atmosphere by evaporation or plant transpiration. Estimated fluxes across confining units in the saturated zone are referred to as interaquifer flow (Chapter 1). Focused recharge is that which occurs directly from a point or line source, such as a stream, on land surface. Focused recharge may vary widely in space and time. If the water table intersects a stream channel, estimates of stream loss are called actual recharge, or just recharge. If the water table lies below the stream channel, estimates are referred to as potential recharge. For simplicity, all vertical water fluxes are referred to as drainage throughout this chapter. Whether the estimated quantity represents actual or potential recharge or drainage depends on the circumstances of each individual study.

Book chapter

Infiltration from an impoundment for coal‐bed natural gas, Powder River Basin, Wyoming: Evolution of water and sediment chemistry

Development of coal‐bed natural gas (CBNG) in the Powder River Basin, Wyoming, has increased substantially in recent years. Among environmental concerns associated with this development is the fate of groundwater removed with the gas. A preferred water‐management option is storage in surface impoundments. As of January 2007, permits for more than 4000 impoundments had been issued within Wyoming. A study was conducted on changes in water and sediment chemistry as water from an impoundment infiltrated the subsurface. Sediment cores were collected prior to operation of the impoundment and after its closure and reclamation. Suction lysimeters were used to collect water samples from beneath the impoundment. Large amounts of chloride (12,300 kg) and nitrate (13,500 kg as N), most of which accumulated naturally in the sediments over thousands of years, were released into groundwater by infiltrating water. Nitrate was more readily flushed from the sediments than chloride. If sediments at other impoundment locations contain similar amounts of chloride and nitrate, impoundments already permitted could release over 48 × 10 6 kg of chloride and 52 × 10 6 kg of nitrate into groundwater in the basin. A solute plume with total dissolved solid (TDS) concentrations at times exceeding 100,000 mg/L was created in the subsurface. TDS concentrations in the plume were substantially greater than those in the CBNG water (about 2300 mg/L) and in the ambient shallow groundwater (about 8000 mg/L). Sulfate, sodium, and magnesium are the dominant ions in the plume. The elevated concentrations are attributed to cation‐exchange‐enhanced gypsum dissolution. As gypsum dissolves, calcium goes into solution and is exchanged for sodium and magnesium on clays. Removal of calcium from solution allows further gypsum dissolution.

Water Resources Research

Water Budgets: Foundations for Effective Water-Resources and Environmental Management

INTRODUCTION Water budgets provide a means for evaluating availability and sustainability of a water supply. A water budget simply states that the rate of change in water stored in an area, such as a watershed, is balanced by the rate at which water flows into and out of the area. An understanding of water budgets and underlying hydrologic processes provides a foundation for effective water-resource and environmental planning and management. Observed changes in water budgets of an area over time can be used to assess the effects of climate variability and human activities on water resources. Comparison of water budgets from different areas allows the effects of factors such as geology, soils, vegetation, and land use on the hydrologic cycle to be quantified. Human activities affect the natural hydrologic cycle in many ways. Modifications of the land to accommodate agriculture, such as installation of drainage and irrigation systems, alter infiltration, runoff, evaporation, and plant transpiration rates. Buildings, roads, and parking lots in urban areas tend to increase runoff and decrease infiltration. Dams reduce flooding in many areas. Water budgets provide a basis for assessing how a natural or human-induced change in one part of the hydrologic cycle may affect other aspects of the cycle. This report provides an overview and qualitative description of water budgets as foundations for effective water-resources and environmental management of freshwater hydrologic systems. Perhaps of most interest to the hydrologic community, the concepts presented are also relevant to the fields of agriculture, atmospheric studies, meteorology, climatology, ecology, limnology, mining, water supply, flood control, reservoir management, wetland studies, pollution control, and other areas of science, society, and industry. The first part of the report describes water storage and movement in the atmosphere, on land surface, and in the subsurface, as well as water exchange among these compartments. Our ability to measure these phenomena and inherent uncertainties in measurement techniques also are discussed. The latter part of the report presents a number of case studies that illustrate how water-budget studies are conducted, documents how human activities affect water budgets, and describes how water budgets are used to address water and environmental issues.

Circular

User's guide to the Variably Saturated Flow (VSF) process to MODFLOW

A new process for simulating three-dimensional (3-D) variably saturated flow (VSF) using Richards' equation has been added to the 3-D modular finite-difference ground-water model MODFLOW. Five new packages are presented here as part of the VSF Process--the Richards' Equation Flow (REF1) Package, the Seepage Face (SPF1) Package, the Surface Ponding (PND1) Package, the Surface Evaporation (SEV1) Package, and the Root Zone Evapotranspiration (RZE1) Package. Additionally, a new Adaptive Time-Stepping (ATS1) Package is presented for use by both the Ground-Water Flow (GWF) Process and VSF. The VSF Process allows simulation of flow in unsaturated media above the ground-water zone and facilitates modeling of ground-water/surface-water interactions. Model performance is evaluated by comparison to an analytical solution for one-dimensional (1-D) constant-head infiltration (Dirichlet boundary condition), field experimental data for a 1-D constant-head infiltration, laboratory experimental data for two-dimensional (2-D) constant-flux infiltration (Neumann boundary condition), laboratory experimental data for 2-D transient drainage through a seepage face, and numerical model results (VS2DT) of a 2-D flow-path simulation using realistic surface boundary conditions. A hypothetical 3-D example case also is presented to demonstrate the new capability using periodic boundary conditions (for example, daily precipitation) and varied surface topography over a larger spatial scale (0.133 square kilometer). The new model capabilities retain the modular structure of the MODFLOW code and preserve MODFLOW's existing capabilities as well as compatibility with commercial pre-/post-processors. The overall success of the VSF Process in simulating mixed boundary conditions and variable soil types demonstrates its utility for future hydrologic investigations. This report presents a new flow package implementing the governing equations for variably saturated ground-water flow, four new boundary condition packages unique to unsaturated flow, the Adaptive Time-Stepping Package for use with both the GWF Process and the new VSF Process, detailed descriptions of the input and output files for each package, and six simulation examples verifying model performance.

Techniques and Methods

Flow and storage in groundwater systems

The dynamic nature of groundwater is not readily apparent, except where discharge is focused at springs or where recharge enters sinkholes. Yet groundwater flow and storage are continually changing in response to human and climatic stresses. Wise development of groundwater resources requires a more complete understanding of these changes in flow and storage and of their effects on the terrestrial environment and on numerous surface-water features and their biota.

Science

Tritium and 14C concentrations in unsaturated-zone gases at test hole UZB-2, Amargosa Desert Research Site, 1994-98: A section in U.S. Geological Survey Toxic Substances Hydrology Program: Proceedings of the technical meeting, Charleston, South Carolina, March 8-12, 1999: Volume 3 (Part C) (WRI 99-4018C)>

Tritium concentrations have been determined yearly since April 1994 from water-vapor samples collected at test hole UZB-2. The hole was drilled about 100 m (meters) south of the southwest corner of a commercial burial site for low-level radioactive wastes in September 1993. UZB-2 is equipped with ten 2.5-cm (centimeters) diameter air ports permanently installed in the unsaturated zone between the depths of 5.5 and 108.8 m below land surface. Depth to ground water is about 110 m. Additional sampling ports were driven by hand to depths of 0.5, 1.0 and 1.5 m in May 1997. Initial samples of water vapor collected in April 1994 showed elevated tritium concentrations of more than 100 TU (tritium units) from all 10 air ports, with a maximum concentration of 762±10 TU from an air port at a depth of 24.1 m. Subsequent tritium concentrations increased in all air ports, although tritium concentrations at depths of less than 34.1 m have remained relatively constant since July 1995. The largest observed increase in tritium has been at a depth of 47.9 m. There, tritium concentration has increased from 198±5 TU in April 1994 to 2,570±30 TU in June 1998. Large increases also have been measured in samples collected from air ports at depths of 106.4 and 108.8 m, just above the water table. During September and October 1998, carbon dioxide samples were collected from all ten air ports in UZB-2 and at a depth of 1.5 m, and analyzed for radioactive carbon-14 ( 14 C). 14 C concentrations are highest in air ports at depths less than 6 m where they exceed 2,000 pmc (percent modern carbon). Concentrations decrease rapidly in air ports at depth and are about 20 pmc below 94.2 m. However, at 47.9 meters, the 14 C concentration is 205±1 pmc, which is 2 to 4 times higher than concentrations in air ports immediately above and below. This depth corresponds to the largest tritium increase in UZB-2. Concentrations of both tritium and 14 C are greater than what could be expected from atmospheric fallout. The distribution of tritium and 14 C likely represent a complex pattern of lateral and vertical transport through the unsaturated zone from buried wastes to UZB-2.

Nevada

Solution of the advection-dispersion equation in two dimensions by a finite-volume Eulerian-Lagrangian localized adjoint method

We extend the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) for solution of the advection-dispersion equation to two dimensions. The method can conserve mass globally and is not limited by restrictions on the size of the grid Peclet or Courant number. Therefore, it is well suited for solution of advection-dominated ground-water solute transport problems. In test problem comparisons with standard finite differences, FVELLAM is able to attain accurate solutions on much coarser space and time grids. On fine grids, the accuracy of the two methods is comparable. A critical aspect of FVELLAM (and all other ELLAMs) is evaluation of the mass storage integral from the preceding time level. In FVELLAM this may be accomplished with either a forward or backtracking approach. The forward tracking approach conserves mass globally and is the preferred approach. The backtracking approach is less computationally intensive, but not globally mass conservative. Boundary terms are systematically represented as integrals in space and time which are evaluated by a common integration scheme in conjunction with forward tracking through time. Unlike the one-dimensional case, local mass conservation cannot be guaranteed, so slight oscillations in concentration can develop, particularly in the vicinity of inflow or outflow boundaries.

Advances in Water Resources

Water and tritium movement through the unsaturated zone at a low-level radioactive-waste disposal site near Sheffield, Illinois, 1981-85

The movement of water and tritium through the unsaturated zone was studied at a low-level radioactive-waste disposal site near Sheffield, Bureau County, Illinois, from 1981 to 1985. Water and tritium movement occurred in an annual, seasonally timed cycle; recharge to the saturated zone generally occurred in the spring and early summer. Mean annual precipitation (1982-85) was 871 mm (millimeters); mean annual recharge to the disposal trenches (July 1982 through June 1984) was estimated to be 107 mm. Average annual tritium flux below the study trenches was estimated to be 3.4 mCi/yr (millicuries per year). Site geology, climate, and waste-disposal practices influenced the spatial and temporal variability of water and tritium movement. Of the components of the water budget, evapotranspiration contributed most to the temporal variability of water and tritium movement. Disposal trenches are constructed in complexly layered glacial and postglacial deposits that average 17 m (meters) in thickness and overlie a thick sequence of Pennsylvanian shale. The horizontal saturated hydraulic conductivity of the clayey-silt to sand-sized glacial and postglacial deposits ranges from 4.8x10 -1 to 3.4x10 4 mm/d (millimeters per day). A 120-m-long horizontal tunnel provided access for hydrologic measurements and collection of sediment and water samples from the unsaturated and saturated geologic deposits below four disposal trenches. Trench-cover and subtrench deposits were monitored with soil-moisture tensiometers, vacuum and gravity lysimeters, piezometers, and a nuclear soil-moisture gage. A cross-sectional, numerical ground-water-flow model was used to simulate water movement in the variably saturated geologic deposits in the tunnel area. Concurrent studies at the site provided water-budget data for estimating recharge to the disposal trenches. Vertical water movement directly above the trenches was impeded by a zone of compaction within the clayey-silt trench covers. Water entered the trenches primarily at the trench edges where the compacted zone was absent and the cover was relatively thin. Collapse holes in the trench covers that resulted from inadequate compaction of wastes within the trenches provided additional preferential pathways for surface-water drainage into the trenches; drainage into one collapse hole during a rainstorm was estimated to be 1,700 L (liters). Till deposits near trench bases induced lateral water and tritium movement. Limited temporal variation in water movement and small flow gradients (relative to the till deposits) were detected in the unsaturated subtrench sand deposit; maximum gradients during the spring recharge period averaged 1.62 mm/mm (millimeter per millimeter). Time-of-travel of water moving from the trench covers to below the trenches was estimated to be as rapid as 41 days (assuming individual water molecules move this distance in one recharge cycle). Tritium concentrations in water from the unsaturated zone ranged from 200 (background) to 10,000,000 pCi/L (picocuries per liter). Tritium concentrations generally were higher below trench bases (averaging 91,000 pCi/L) than below intertrench sediments (averaging 3,300 pCi/L), and in the subtrench Toulon Member of the Glasford Formation (sand) (averaging 110,000 pCi/L) than in the Hulick Till Member of the Glasford Formation (clayey silt) (averaging 59,000 pCi/L). Average subtrench tritium concentration increased from 28,000 to 100,000 pCi/L during the study period. Within the trench covers, there was a strong seasonal trend in tritium concentrations; the highest concentrations occurred in late summer when soil-moisture contents were at a minimum. Subtrench tritium movement occurred in association with the annual cycle of water movement, as well as independently of the cycle, in apparent response to continuous water movement through the subtrench sand deposits and to the deterioration of trench-waste containers. The increase in concen

Water Supply Paper

Variability of an unsaturated sand unit underlying a radioactive- waste trench

Properties of soils vary considerably within any field. This study was conducted to investigate the variability in properties of an unsaturated sand unit that lies at a depth of 13 m below land surface. Four-hundred-forty soil core samples, obtained from a 1.75 by 18-m horizontal plane within a sand unit underlying a waste trench, were used to describe the variability of moisture content, tritium concentration, and several physical properties. A simple model based on unit-gradient theory was used to calculate fluxes of water and tritium through the study plane. The effects of including spatial variability of properties in model calculations were examined. The sand was uniform in terms of median particle size, bulk density, and proposity, but several variables, including saturated hydraulic conductivity ( K s ) and tritium concentration, showed considerable variability. In general, correlation coefficients calculated between variables were low. Spatial correlation structures were identified for moisture content, tritium concentration, and a few particle-size variables, but were absent for K s and most other variables. Variances of K s and field-measured volumetric moisture content (θ) were lower than published values for soils. Model-calculated fluxes of water and tritium that accounted for spatial variability of properties were considerably higher than those that did not account for the variability. Model results were more sensitive to changes in θ and residual moisture content (each parameters of the relative hydraulic conductivity function) than to changes in K s .

Illinois