USGS ScienceSearch

Geology topics

Richard W. Healy

Publications and source records attributed to Richard W. Healy.

At least 19 recordsLinked to original sources

New capabilities in MT3D-USGS for simulating unsaturated-zone heat transport

Changes in climate and land use will alter groundwater heat transport dynamics in the future. These changes will in turn affect watershed processes (e.g., nutrient cycling) as well as watershed characteristics (e.g., distribution and persistence of cold-water habitat). Thus, groundwater flow and heat transport models at watershed scales that can characterize and quantify thermal impacts of surface temperature change on groundwater system temperatures are needed to forecast changes to groundwater-linked ecosystems in riparian zones, streams, and lakes. Including unsaturated zone processes has previously been shown to be important for properly determining the timing and magnitude of groundwater recharge (Hunt et al. 2008). Similarly, heat transport dynamics in the saturated-zone, as well as connected surface-water systems, can be appreciably influenced by unsaturated-zone processes; in this way the unsaturated zone forms an inextricable link between land surface where change occurs and the groundwater system that transmit that change. This paper presents new capabilities for the existing MT3D-USGS transport simulator by adding functionality for simulating heat transport through the unsaturated zone. New simulation capabilities are verified through comparison of simulation results with those of the variably-saturated heat transport simulator VS2DH under steady and transient conditions for both water and heat flow. The new capabilities are assessed using a number of conceptualizations and include evaluations of convective and conductive heat flow. These additional capabilities increase the utility for applied watershed-scale simulations, which in turn should facilitate more realistic characterizations of temperature change on thermally sensitive ecosystems, such as stream habitat.

Groundwater

VS2DRTI: Simulating heat and reactive solute transport in variably saturated porous media

Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid-rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult-to-use models. To address the need for a simple and easy-to-use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two-dimensional, constant-density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature-dependent cation exchange. VS2DRTI is freely available public domain software.

Groundwater

Chemical tracer methods

Tracers have a wide variety of uses in hydrologic studies: providing quantitative or qualitative estimates of recharge, identifying sources of recharge, providing information on velocities and travel times of water movement, assessing the importance of preferential flow paths, providing information on hydrodynamic dispersion, and providing data for calibration of water flow and solute-transport models (Walker, 1998; Cook and Herczeg, 2000; Scanlon et al ., 2002b). Tracers generally are ions, isotopes, or gases that move with water and that can be detected in the atmosphere, in surface waters, and in the subsurface. Heat also is transported by water; therefore, temperatures can be used to trace water movement. This chapter focuses on the use of chemical and isotopic tracers in the subsurface to estimate recharge. Tracer use in surface-water studies to determine groundwater discharge to streams is addressed in Chapter 4; the use of temperature as a tracer is described in Chapter 8. Following the nomenclature of Scanlon et al . (2002b), tracers are grouped into three categories: natural environmental tracers, historical tracers, and applied tracers. Natural environmental tracers are those that are transported to or created within the atmosphere under natural processes; these tracers are carried to the Earth’s surface as wet or dry atmospheric deposition. The most commonly used natural environmental tracer is chloride (Cl) (Allison and Hughes, 1978). Ocean water, through the process of evaporation, is the primary source of atmospheric Cl. Other tracers in this category include chlorine-36 ( 36 Cl) and tritium ( 3 H); these two isotopes are produced naturally in the Earth’s atmosphere; however, there are additional anthropogenic sources of them.

Book chapter

1DTempPro V2: new features for inferring groundwater/surface-water exchange

A new version of the computer program 1DTempPro extends the original code to include new capabilities for (1) automated parameter estimation, (2) layer heterogeneity, and (3) time-varying specific discharge. The code serves as an interface to the U.S. Geological Survey model VS2DH and supports analysis of vertical one-dimensional temperature profiles under saturated flow conditions to assess groundwater/surface-water exchange and estimate hydraulic conductivity for cases where hydraulic head is known.

Groundwater Monitoring & Remediation

Variably-saturated groundwater modeling for optimizing managed aquifer recharge using trench infiltration

Spreading-basin methods have resulted in more than 130 million cubic meters of recharge to the unconfined Navajo Sandstone of southern Utah in the past decade, but infiltration rates have slowed in recent years because of reduced hydraulic gradients and clogging. Trench infiltration is a promising alternative technique for increasing recharge and minimizing evaporation. This paper uses a variably saturated flow model to further investigate the relative importance of the following variables on rates of trench infiltration to unconfined aquifers: saturated hydraulic conductivity, trench spacing and dimensions, initial water-table depth, alternate wet/dry periods, and number of parallel trenches. Modeling results showed (1) increased infiltration with higher hydraulic conductivity, deeper initial water tables, and larger spacing between parallel trenches, (2) deeper or wider trenches do not substantially increase infiltration, (3) alternating wet/dry periods result in less overall infiltration than keeping the trenches continuously full, and (4) larger numbers of parallel trenches within a fixed area increases infiltration but with a diminishing effect as trench spacing becomes tighter. An empirical equation for estimating expected trench infiltration rates as a function of hydraulic conductivity and initial water-table depth was derived and can be used for evaluating feasibility of trench infiltration in other hydrogeologic settings

Utah

The water-energy nexus: An earth science perspective

Water availability and use are closely connected with energy development and use. Water cannot be delivered to homes, businesses, and industries without energy, and most forms of energy development require large amounts of water. The United States faces two significant and sometimes competing challenges: to provide sustainable supplies of freshwater for humans and ecosystems and to ensure adequate sources of energy for future generations. This report reviews the complex ways in which water and energy are interconnected and describes the earth science data collection and research that can help the Nation address these important challenges. The earth sciences have been a cornerstone in developing our current understanding of the water-energy nexus. A full understanding of the nexus, however, is limited by uncertainty in our knowledge of fundamental issues, such as the quantity of freshwater that is available, the amount of water that is used in energy development, the effects that emerging energy development technologies have on water quality and quantity, and the amount of energy required to treat and deliver freshwater. Enhanced data collection and research can improve our understanding of these important issues and thereby lay the groundwork for informed resource management. Relevant earth science issues analyzed and discussed herein include freshwater availability; water use; ecosystems health; assessment of saline water resources; assessment of fossil-fuel, uranium, and geothermal resources; subsurface injection of wastewater and carbon dioxide and related induced seismicity; climate change and its effect on water availability and energy production; byproducts and waste streams of energy development; emerging energy-development technologies; and energy for water treatment and delivery.

Circular

A framework for assessing water and proppant use and flowback water extraction associated with development of continuous petroleum resources

The U.S. Geological Survey is developing approaches for the quantitative assessment of water and proppant involved with possible future production of continuous petroleum deposits. The assessment approach is an extension of existing U.S. Geological Survey petroleum-assessment methods, and it aims to provide objective information that helps decision makers understand the tradeoffs inherent in resource-development decisions. This fact sheet provides an overview of U.S. Geological Survey assessments for quantities of water and proppant required for drilling and hydraulic fracturing and for flowback water extracted with petroleum; the report also presents the form of the intended assessment output information.

Fact Sheet

Estimating groundwater recharge

Understanding groundwater recharge is essential for successful management of water resources and modeling fluid and contaminant transport within the subsurface. This book provides a critical evaluation of the theory and assumptions that underlie methods for estimating rates of groundwater recharge. Detailed explanations of the methods are provided - allowing readers to apply many of the techniques themselves without needing to consult additional references. Numerous practical examples highlight benefits and limitations of each method. Approximately 900 references allow advanced practitioners to pursue additional information on any method. For the first time, theoretical and practical considerations for selecting and applying methods for estimating groundwater recharge are covered in a single volume with uniform presentation. Hydrogeologists, water-resource specialists, civil and agricultural engineers, earth and environmental scientists and agronomists will benefit from this informative and practical book. It can serve as the primary text for a graduate-level course on groundwater recharge or as an adjunct text for courses on groundwater hydrology or hydrogeology.

Book

Deep subsurface drip irrigation using coal-bed sodic water: part II. geochemistry

Waters with low salinity and high sodium adsorption ratios (SARs) present a challenge to irrigation because they degrade soil structure and infiltration capacity. In the Powder River Basin of Wyoming, such low salinity (electrical conductivity, EC 2.1 mS cm -1 ) and high-SAR (54) waters are co-produced with coal-bed methane and some are used for subsurface drip irrigation(SDI). The SDI system studied mixes sulfuric acid with irrigation water and applies water year-round via drip tubing buried 92 cm deep. After six years of irrigation, SAR values between 0 and 30 cm depth (0.5-1.2) are only slightly increased over non-irrigated soils (0.1-0.5). Only 8-15% of added Na has accumulated above the drip tubing. Sodicity has increased in soil surrounding the drip tubing, and geochemical simulations show that two pathways can generate sodic conditions. In soil between 45-cm depth and the drip tubing, Na from the irrigation water accumulates as evapotranspiration concentrates solutes. SAR values >12, measured by 1:1 water-soil extracts, are caused by concentration of solutes by factors up to 13. Low-EC (<0.7 mS cm -1 ) is caused by rain and snowmelt flushing the soil and displacing ions in soil solution. Soil below the drip tubing experiences lower solute concentration factors (1-1.65) due to excess irrigation water and also contains relatively abundant native gypsum (2.4 &plusmn; 1.7 wt.%). Geochemical simulations show gypsum dissolution decreases soil-water SAR to <7 and increases the EC to around 4.1 mS cm-1, thus limiting negative impacts from sodicity. With sustained irrigation, however, downward flow of excess irrigation water depletes gypsum, increasing soil-water SAR to >14 and decreasing EC in soil water to 3.2 mS cm-1. Increased sodicity in the subsurface, rather than the surface, indicates that deep SDI can be a viable means of irrigating with sodic waters.

Wyoming

Deep subsurface drip irrigation using coal-bed sodic water: part I. water and solute movement

Water co-produced with coal-bed methane (CBM) in the semi-arid Powder River Basin of Wyoming and Montana commonly has relatively low salinity and high sodium adsorption ratios that can degrade soil permeability where used for irrigation. Nevertheless, a desire to derive beneficial use from the water and a need to dispose of large volumes of it have motivated the design of a deep subsurface drip irrigation (SDI) system capable of utilizing that water. Drip tubing is buried 92 cm deep and irrigates at a relatively constant rate year-round, while evapotranspiration by the alfalfa and grass crops grown is seasonal. We use field data from two sites and computer simulations of unsaturated flow to understand water and solute movements in the SDI fields. Combined irrigation and precipitation exceed potential evapotranspiration by 300-480 mm annually. Initially, excess water contributes to increased storage in the unsaturated zone, and then drainage causes cyclical rises in the water table beneath the fields. Native chloride and nitrate below 200 cm depth are leached by the drainage. Some CBM water moves upward from the drip tubing, drawn by drier conditions above. Chloride from CBM water accumulates there as root uptake removes the water. Year over year accumulations indicated by computer simulations illustrate that infiltration of precipitation water from the surface only partially leaches such accumulations away. Field data show that 7% and 27% of added chloride has accumulated above the drip tubing in an alfalfa and grass field, respectively, following 6 years of irrigation. Maximum chloride concentrations in the alfalfa field are around 45 cm depth but reach the surface in parts of the grass field, illustrating differences driven by crop physiology. Deep SDI offers a means of utilizing marginal quality irrigation waters and managing the accumulation of their associated solutes in the crop rooting zone.

Wyoming

Modeling variably saturated subsurface solute transport with MODFLOW-UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated-zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume-averaged approach in which Lagrangian-based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF-MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably-Saturated Two-Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two- and three-dimensional simulations also were investigated to ensure unsaturated-saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large-scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF-MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three-dimensional variably saturated flow and transport simulations revealed UZF-MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide-spread use of both MODFLOW and MT3DMS, the added capability of unsaturated-zone transport in this familiar modeling framework stands to benefit a broad user-ship.

Ground Water

VS2DI: Model use, calibration, and validation

VS2DI is a software package for simulating water, solute, and heat transport through soils or other porous media under conditions of variable saturation. The package contains a graphical preprocessor for constructing simulations, a postprocessor for displaying simulation results, and numerical models that solve for flow and solute transport (VS2DT) and flow and heat transport (VS2DH). Flow is described by the Richards equation, and solute and heat transport are described by advection-dispersion equations; the finite-difference method is used to solve these equations. Problems can be simulated in one, two, or three (assuming radial symmetry) dimensions. This article provides an overview of calibration techniques that have been used with VS2DI; included is a detailed description of calibration procedures used in simulating the interaction between groundwater and a stream fed by drainage from agricultural fields in central Indiana. Brief descriptions of VS2DI and the various types of problems that have been addressed with the software package are also presented.

Transactions of the ASABE

Hydrologic and geochemical data collected near Skewed Reservoir, an impoundment for coal-bed natural gas produced water, Powder River Basin, Wyoming

The Powder River Structural Basin is one of the largest producers of coal-bed natural gas (CBNG) in the United States. An important environmental concern in the Basin is the fate of groundwater that is extracted during CBNG production. Most of this produced water is disposed of in unlined surface impoundments. A 6-year study of groundwater flow and subsurface water and soil chemistry was conducted at one such impoundment, Skewed Reservoir. Hydrologic and geochemical data collected as part of that study are contained herein. Data include chemistry of groundwater obtained from a network of 21 monitoring wells and three suction lysimeters and chemical and physical properties of soil cores including chemistry of water/soil extracts, particle-size analyses, mineralogy, cation-exchange capacity, soil-water content, and total carbon and nitrogen content of soils.

Wyoming

Influence of a thin veneer of low-hydraulic-conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells.

Colorado

Tracking solutes and water from subsurface drip irrigation application of coalbed methane-produced waters, Powder River Basin, Wyoming

One method to beneficially use water produced from coalbed methane (CBM) extraction is subsurface drip irrigation (SDI) of croplands. In SDI systems, treated CBM water (injectate) is supplied to the soil at depth, with the purpose of preventing the buildup of detrimental salts near the surface. The technology is expanding within the Powder River Basin, but little research has been published on its environmental impacts. This article reports on initial results from tracking water and solutes from the injected CBM-produced waters at an SDI system in Johnson County, Wyoming. In the first year of SDI operation, soil moisture significantly increased in the SDI areas, but well water levels increased only modestly, suggesting that most of the water added was stored in the vadose zone or lost to evapotranspiration. The injectate has lower concentrations of most inorganic constituents relative to ambient groundwater at the site but exhibits a high sodium adsorption ratio. Changes in groundwater chemistry during the same period of SDI operation were small; the increase in groundwater-specific conductance relative to pre-SDI conditions was observed in a single well. Conversely, groundwater samples collected beneath another SDI field showed decreased concentrations of several constituents since the SDI operation. Groundwater-specific conductance at the 12 other wells showed no significant changes. Major controls on and compositional variability of groundwater, surface water, and soil water chemistry are discussed in detail. Findings from this research provide an understanding of water and salt dynamics associated with SDI systems using CBM-produced water.

Wyoming

Numerical investigation of steady liquid water flow in a variably saturated fracture network

Numerical simulation was used to study steady liquid water movement in a 5-m by 5-m vertical section containing a hypothetical fracture network under conditions of variable saturation. The fracture network was assumed to be embedded within an impermeable rock matrix. Three variations of a network were considered. The “mixed” network consisted of two fracture sets, a subvertical set containing five 125 μm average aperture fractures and a subhorizontal set containing four 25 μm average aperture fractures. The other two networks had identical fracture orientation and contained either all 125 μm or all 25 μm average aperture fractures. The TOUGH simulator was used to calculate the total steady liquid water flux through the network, the flux through individual fracture segments, and the pressure head at each fracture segment. A unit hydraulic gradient was imposed on the network by applying fixed pressure head boundaries (ranging from −0.25 to 0.0 m of water) of equal value to the top and bottom. Saturation and permeability versus pressure head relations for the two sets of fractures were determined with the VSFRAC model, which assumed that aperture was variable within an individual fracture. Results showed that the spatial distributions of pressure head and flux within the network, as well as the location of the dominant pathways, depended strongly on the prescribed boundary pressure head. For the mixed network, both pressure head and flux tended to become more spatially uniform when the boundary pressure head approached the pressure head at which the permeability thickness products of the large- and small-aperture fractures are equal (the crossover pressure head). These results imply that for systems similar to the one considered here, interpretation of actual measurements of pressure head and flux may be quite complex, and that representation of variably saturated fracture networks as an equivalent continuum may be more valid for some ranges in pressure head than for others. Equivalent permeability as a function of pressure head was calculated for the fracture network, illustrating how information collected on individual fractures may be used to estimate the flow properties of rock at larger scales.

Water Resources Research

Water-budget methods

A water budget is an accounting of water movement into and out of, and storage change within, some control volume. Universal and adaptable are adjectives that reflect key features of water-budget methods for estimating recharge. The universal concept of mass conservation of water implies that water-budget methods are applicable over any space and time scales (Healy et al ., 2007). The water budget of a soil column in a laboratory can be studied at scales of millimeters and seconds. A water-budget equation is also an integral component of atmospheric general circulation models used to predict global climates over periods of decades or more. Water-budget equations can be easily customized by adding or removing terms to accurately portray the peculiarities of any hydrologic system. The equations are generally not bound by assumptions on mechanisms by which water moves into, through, and out of the control volume of interest. So water-budget methods can be used to estimate both diffuse and focused recharge, and recharge estimates are unaffected by phenomena such as preferential flow paths within the unsaturated zone. Water-budget methods represent the largest class of techniques for estimating recharge. Most hydrologic models are derived from a water-budget equation and can therefore be classified as water-budget models. It is not feasible to address all water-budget methods in a single chapter. This chapter is limited to discussion of the “residual” water-budget approach, whereby all variables in a water-budget equation, except for recharge, are independently measured or estimated and recharge is set equal to the residual. This chapter is closely linked with Chapter 3, on modeling methods, because the equations presented here form the basis of many models and because models are often used to estimate individual components in water-budget studies. Water budgets for streams and other surface-water bodies are addressed in Chapter 4. The use of soil-water budgets and lysimeters for determining potential recharge and evapotranspiration from changes in water storage is discussed in Chapter 5. Aquifer water-budget methods based on the measurement of groundwater levels are described in Chapter 6.

Book chapter