Search USGS⌕ Search

Geology topics

R. W. Healy

Publications and source records attributed to R. W. Healy.

At least 19 recordsLinked to original sources
USGS2011

Groundwater chemistry near an impoundment for produced water, Powder River Basin, Wyoming, USA

The Powder River 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 the large amounts of groundwater extracted during CBNG production. Most of this produced water is disposed of in unlined surface impoundments. A 6-year study of groundwater flow and water chemistry at one impoundment, Skewed Reservoir, has produced the most detailed data set for any impoundment in the Basin. Data were collected from a network of 21 observation wells and three suction lysimeters. A groundwater mound formed atop bedrock within initially unsaturated, unconsolidated deposits underlying the reservoir. Heterogeneity in physical and chemical properties of sediments resulted in complex groundwater flow paths and highly variable groundwater chemistry. Sulfate, bicarbonate, sodium, and magnesium were the dominant ions in all areas, but substantial variability existed in relative concentrations; pH varied from less than 3 to more than 9, and total dissolved solids concentrations ranged from less than 5000 to greater than 100,000 mg/L. Selenium was a useful tracer of reservoir water; selenium concentrations exceeded 300 μg/L in samples obtained from 18 of the 24 sampling points. Groundwater travel time from the reservoir to a nearby alluvial aquifer (a linear distance of 177 m) was calculated at 474 days on the basis of selenium concentrations. The produced water is not the primary source of solutes in the groundwater. Naturally occurring salts and minerals within the unsaturated zone, dissolved and mobilized by infiltrating impoundment water, account for most of the solute mass in groundwater. Gypsum dissolution, cation-exchange, and pyrite oxidation appear to be important reactions. The complex geochemistry and groundwater flow paths at the study site underscore the difficulty in assessing effects of surface impoundments on water resources within the Powder River Basin.

Wyoming↗
USGS2010

A finite-volume Eulerian-Lagrangian Localized Adjoint Method for solution of the advection-dispersion equation

A new mass-conservative method for solution of the one-dimensional advection-dispersion equation is derived and discussed. Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods, in terms of accuracy and efficiency, for solute transport problems that are dominated by advection. For dispersion-dominated problems, the performance of the method is similar to that of standard methods. Like previous ELLAM formulations, FVELLAM systematically conserves mass globally with all types of boundary conditions. FVELLAM differs from other ELLAM approaches in that integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking, as used by most characteristic methods, of characteristic lines intersecting inflow boundaries. FVELLAM extends previous ELLAM results by obtaining mass conservation locally on Lagrangian space-time elements. Details of the integration, tracking, and boundary algorithms are presented. Test results are given for problems in Cartesian and radial coordinates.

Water Resources Research↗
USGS2008

Water movement within the unsaturated zone in four agricultural areas of the United States

Millions of tons of agricultural fertilizer and pesticides are applied annually in the USA. Due to the potential for these chemicals to migrate to groundwater, a study was conducted in 2004 using field data to calculate water budgets, rates of groundwater recharge and times of water travel through the unsaturated zone and to identify factors that influence these phenomena. Precipitation was the only water input at sites in Indiana and Maryland; irrigation accounted for about 80% of total water input at sites in California and Washington. Recharge at the Indiana site (47.5 cm) and at the Maryland site (31.5 cm) were equivalent to 51 and 32%, respectively, of annual precipitation and occurred between growing seasons. Recharge at the California site (42.3 cm) and Washington site (11.9 cm) occurred in response to irrigation events and was about 29 and 13% of total water input, respectively. Average residence time of water in the unsaturated zone, calculated using a piston-flow approach, ranged from less than 1 yr at the Indiana site to more than 8 yr at the Washington site. Results of bromide tracer tests indicate that at three of the four sites, a fraction of the water applied at land surface may have traveled to the water table in less than 1 yr. The timing and intensity of precipitation and irrigation were the dominant factors controlling recharge, suggesting that the time of the year at which chemicals are applied may be important for chemical transport through the unsaturated zone. Copyright ?? 2008 by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America. All rights reserved.

Journal of Environmental Quality↗
USGS2008

Variations in pesticide leaching related to land use, pesticide properties, and unsaturated zone thickness

Pesticide leaching through variably thick soils beneath agricultural fields in Morgan Creek, Maryland was simulated for water years 1995 to 2004 using LEACHM (Leaching Estimation and Chemistry Model). Fifteen individual models were constructed to simulate five depths and three crop rotations with associated pesticide applications. Unsaturated zone thickness averaged 4.7 m but reached a maximum of 18.7 m. Average annual recharge to ground water decreased from 15.9 to 11.1 cm as the unsaturated zone increased in thickness from 1 to 10 m. These point estimates of recharge are at the lower end of previously published values, which used methods that integrate over larger areas capturing focused recharge in the numerous detention ponds in the watershed. The total amount of applied and leached masses for five parent pesticide compounds and seven metabolites were estimated for the 32-km 2 Morgan Creek watershed by associating each hectare to the closest one-dimensional model analog of model depth and crop rotation scenario as determined from land-use surveys. LEACHM parameters were set such that branched, serial, first-order decay of pesticides and metabolites was realistically simulated. Leaching is predicted to be greatest for shallow soils and for persistent compounds with low sorptivity. Based on simulation results, percent parent compounds leached within the watershed can be described by a regression model of the form e −depth ( a ln t ½ −b ln K OC ) where t 1/2 is the degradation half-life in aerobic soils, K OC is the organic carbon normalized sorption coefficient, and a and b are fitted coefficients ( R 2 = 0.86, p value = 7 × 10 −9 ).

Journal of Environmental Quality↗
USGS2008

Simulating water, solute, and heat transport in the subsurface with the VS2DI software package

The software package VS2DI was developed by the U.S. Geological Survey for simulating water, solute, and heat transport in variably saturated porous media. The package consists of a graphical preprocessor to facilitate construction of a simulation, a postprocessor for visualizing simulation results, and two numerical models that solve for flow and solute transport (VS2DT) and flow and heat transport (VS2DH). The finite-difference method is used to solve the Richards equation for flow and the advection-dispersion equation for solute or heat transport. This study presents a brief description of the VS2DI package, an overview of the various types of problems that have been addressed with the package, and an analysis of the advantages and limitations of the package. A review of other models and modeling approaches for studying water, solute, and heat transport also is provided. ?? Soil Science Society of America. All rights reserved.

Vadose Zone Journal↗
USGS2007

Water balance at a low-level radioactive-waste disposal site

The water balance at a low-level radioactive-waste disposal site in northwestern Illinois was studied from July 1982 through June 1984. Continuous data collection allowed estimates to be made for each component of the water-balance equation independent of other components. The average annual precipitation was 948 millimeters. Average annual evapotranspiration was estimated at 637 millimeters, runoff was 160 millimeters, change in water storage in a waste-trench cover was 24 millimeters, and deep percolation was 208 millimeters. The magnitude of the difference between precipitation and all other components (81 millimeters per year) indicates that, in a similar environment, the water-budget method would be useful in estimating evapotranspiration, but questionable for estimation of other components. Precipitation depth and temporal distribution had a very strong effect on all other components of the water-balance equation. Due to the variability of precipitation from year to year, it appears that two years of data are inadequate for characterization of the long-term average water balance at the site.

Water Resources Bulletin↗
USGS2007

Effects of topography and soil properties on recharge at two sites in an agricultural field

Field experiments were conducted from 1992 to 1995 to estimate ground water recharge rates at two sites located within a 2.7-hectare agricultural field. The field lies in a sand plain setting in central Minnesota and is cropped continuously in field corn. The sites are located at a topographically high (upland) site and a topographically low (lowland) site in an effort to quantify the effects of depression focusing of recharge. Three site-specific methods were used to estimate recharge rates: well hydrograph analysis, chlorofluorocarbon age dating, and an unsaturated zone water balance. All three recharge methods indicated that recharge rates at the lowland site (annual average of all methods of 29 cm) exceeded those at the upland site (annual average of 18 cm). On an annual basis, estimates by the individual methods ranged from 12 to 44 percent of precipitation at the upland site and from 21 to 83 percent at the lowland site. The difference in recharge rates between the sites is primarily attributed to depression focusing of surface water runon at the lowland site. However, two other factors were also important: the presence of thin lamellae at the upland site, and coarser textured soils below a depth of 1.5 m at the lowland site.

Minnesota↗
USGS2007

Factors influencing ground-water recharge in the eastern United States

Ground-water recharge estimates for selected locations in the eastern half of the United States were obtained by Darcian and chloride-tracer methods and compared using statistical analyses. Recharge estimates derived from unsaturated-zone (RUZC) and saturated-zone (RSZC) chloride mass balance methods are less variable (interquartile ranges or IQRs are 9.5 and 16.1 cm/yr, respectively) and more strongly correlated with climatic, hydrologic, land use, and sediment variables than Darcian estimates (IQR = 22.8 cm/yr). The unit-gradient Darcian estimates are a nonlinear function of moisture content and also reflect the uncertainty of pedotransfer functions used to estimate hydraulic parameters. Significance level is <0.001 for nearly all explanatory variables having correlations with RUZC of <-0.3 or >0.3. Estimates of RSZC were evaluated using analysis of variance, multiple comparison tests, and an exploratory nonlinear regression (NLR) model. Recharge generally is greater in coastal plain surficial aquifers, fractured crystalline rocks, and carbonate rocks, or in areas with high sand content. Westernmost portions of the study area have low recharge, receive somewhat less precipitation, and contain fine-grained sediment. The NLR model simulates water input to the land surface followed by transport to ground water, depending on factors that either promote or inhibit water infiltration. The model explains a moderate amount of variation in the data set (coefficient of determination = 0.61). Model sensitivity analysis indicates that mean annual runoff, air temperature, and precipitation, and an index of ground-water exfiltration potential most influence estimates of recharge at sampled sites in the region. Soil characteristics and land use have less influence on the recharge estimates, but nonetheless are significant in the NLR model. ?? 2006 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2007

Comparison of local- to regional-scale estimates of ground-water recharge in Minnesota, USA

Regional ground-water recharge estimates for Minnesota were compared to estimates made on the basis of four local- and basin-scale methods. Three local-scale methods (unsaturated-zone water balance, water-table fluctuations (WTF) using three approaches, and age dating of ground water) yielded point estimates of recharge that represent spatial scales from about 1 to about 1000 m 2 . A fourth method (RORA, a basin-scale analysis of streamflow records using a recession-curve-displacement technique) yielded recharge estimates at a scale of 10&ndash;1000s of km 2 . The RORA basin-scale recharge estimates were regionalized to estimate recharge for the entire State of Minnesota on the basis of a regional regression recharge (RRR) model that also incorporated soil and climate data. Recharge rates estimated by the RRR model compared favorably to the local and basin-scale recharge estimates. RRR estimates at study locations were about 41% less on average than the unsaturated-zone water-balance estimates, ranged from 44% greater to 12% less than estimates that were based on the three WTF approaches, were about 4% less than the age dating of ground-water estimates, and were about 5% greater than the RORA estimates. Of the methods used in this study, the WTF method is the simplest and easiest to apply. Recharge estimates made on the basis of the UZWB method were inconsistent with the results from the other methods. Recharge estimates using the RRR model could be a good source of input for regional ground-water flow models; RRR model results currently are being applied for this purpose in USGS studies elsewhere.

Minnesota↗
USGS2002

Using groundwater levels to estimate recharge

Accurate estimation of groundwater recharge is extremely important for proper management of groundwater systems. Many different approaches exist for estimating recharge. This paper presents a review of methods that are based on groundwater-level data. The water-table fluctuation method may be the most widely used technique for estimating recharge; it requires knowledge of specific yield and changes in water levels over time. Advantages of this approach include its simplicity and an insensitivity to the mechanism by which water moves through the unsaturated zone. Uncertainty in estimates generated by this method relate to the limited accuracy with which specific yield can be determined and to the extent to which assumptions inherent in the method are valid. Other methods that use water levels (mostly based on the Darcy equation) are also described. The theory underlying the methods is explained. Examples from the literature are used to illustrate applications of the different methods.

Hydrogeology Journal↗
USGS2000

Chamber measurement of surface-atmosphere trace gas exchange: Numerical evaluation of dependence on soil interfacial layer, and source/sink products

We employed a three-dimensional finite difference gas diffusion model to simulate the performance of chambers used to measure surface-atmosphere tace gas exchange. We found that systematic errors often result from conventional chamber design and deployment protocols, as well as key assumptions behind the estimation of trace gas exchange rates from observed concentration data. Specifically, our simulationshowed that (1) when a chamber significantly alters atmospheric mixing processes operating near the soil surface, it also nearly instantaneously enhances or suppresses the postdeployment gas exchange rate, (2) any change resulting in greater soil gas diffusivity, or greater partitioning of the diffusing gas to solid or liquid soil fractions, increases the potential for chamber-induced measurement error, and (3) all such errors are independent of the magnitude, kinetics, and/or distribution of trace gas sources, but greater for trace gas sinks with the same initial absolute flux. Finally, and most importantly, we found that our results apply to steady state as well as non-steady-state chambers, because the slow rate of gas diffusion in soil inhibits recovery of the former from their initial non-steady-state condition. Over a range of representative conditions, the error in steady state chamber estimates of the trace gas flux varied from -30 to +32%, while estimates computed by linear regression from non-steadystate chamber concentrations were 2 to 31% too small. Although such errors are relatively small in comparison to the temporal and spatial variability characteristic of trace gas exchange, they bias the summary statistics for each experiment as well as larger scale trace gas flux estimates based on them.

Journal of Geophysical Research D: Atmospheres↗
USGS2000

Treatment of internal sources in the finite-volume ELLAM

The finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) is a mass-conservative approach for solving the advection-dispersion equation. The method has been shown to be accurate and efficient for solving advection-dominated problems of solute transport in ground water in 1, 2, and 3 dimensions. Previous implementations of FVELLAM have had difficulty in representing internal sources because the standard assumption of lowest order Raviart-Thomas velocity field does not hold for source cells. Therefore, tracking of particles within source cells is problematic. A new approach has been developed to account for internal sources in FVELLAM. It is assumed that the source is uniformly distributed across a grid cell and that instantaneous mixing takes place within the cell, such that concentration is uniform across the cell at any time. Sub-time steps are used in the time-integration scheme to track mass outflow from the edges of the source cell. This avoids the need for tracking within the source cell. We describe the new method and compare results for a test problem with a wide range of cell Peclet numbers.

Conference Paper↗
USGS1997

Hydrogeologic and water-quality data used to evaluate the effects of focused recharge on ground-water quality near Princeton, Minnesota, 1991-95

This study was part of the Management Systems Evaluation Area (MSEA) Program, a multi-scale, inter-agency initiative to evaluate the effects of agricultural systems on water quality in the midwest corn belt. The research was part of the U.S. Geological Survey Toxics Substances Hydrology Program. The research area was located in the Anoka Sand Plain about 5 kilometers southwest of Princeton, Minnesota. This report presents temperature, precipitation, soil-moisture, water-quality, and other data collected during 1991-95 to evaluate the effects of focused recharge on ground-water quality at the Princeton MSEA.

Minnesota↗
USGS1996

Factors affecting tritium and 14carbon distributions in the unsaturated zone near the low-level radioactive-waste burial site south of Beatty, Nevada

Interpretations of the distributions of tritiated water vapor (HTO V ) and 14 carbon dioxide gas ( 14 CO 2 ) concentrations in the unsaturated zone adjacent to the low-level radioactive-waste burial site south of Beatty, Nevada, suggest that observed concentrations of 14 CO 2 could be explained by either diffusive or advective transport of the radioactive gas from the site. The distribution of HTO V cannot be explained by vapor transport, either by diffusive or advective mechanisms. Thus, liquid transport appears to have played a role in moving HTO V to well UZB-2. Although the process by which this occurred cannot be determined from available data (and indeed may never be known), it is likely that liquid wastes disposed directly into the trenches during the period from 1962 to 1975 contributed to the offsite contamination. Liquid transport may have been enhanced by precipitation and runoff into open trenches that resulted in the occasional accumulation of ponded water in the trenches and flow along preferential pathways in the underlying unsaturated zone

Nevada↗
USGS1992

Effects of low-level radioactive-waste disposal on water chemistry in the unsaturated zone at a site near Sheffield, Illinois, 1982-84

A 1982-84 field study defined the chemistry of water collected from the unsaturated zone at a low-level radioactive-waste disposal site near Sheffield, Bureau County, Illinois. Chemical data were evaluated to determine the principal naturally occurring geochemical reactions in the unsaturated zone and to evaluate waste-induced effects on pore-water chemistry. Samples of precipitation, unsaturated-zone pore water, and saturated-zone water were analyzed for specific conductance, pH, alkalinity, major cations and anions, dissolved organic carbon, gross alpha and beta radiation, and tritium. Little change in concentration of most major constituents in the unsaturated-zone water was observed with respect to depth or distance from disposal trenches. Tritium and dissolved organic carbon concentrations were, however, dependent on proximity to trenches. The primary reactions, both on- site and off-site, were carbonate and clay dissolution, cation exchange, and the oxidation of pyrite. The major difference between on-site and off-site inorganic water chemistry resulted from the removal of the Roxana Silt and the Radnor Till Member of the Glasford Formation from on-site. Off-site, the Roxana Silt contributed substantial quantities of sodium to solution from montmorillonite dissolution and associated cation-exchange reactions. The Radnor Till Member provided exchange surfaces for magnesium. Precipitation at the site had an ionic composition of calcium zinc sulfate and an average pH of 4.6. Within 0.3 meter of the land surface, infiltrating rain water or snowmelt changed to an ionic canposition of calcium sulfate off-site and calcium bicarbonate on-site and had an average pH of 7.9; below that depth, pH averaged 7.5 and the ionic composition generally was calcium magnesium bicarbonate. Alkalinity and specific conductance differed primarily according to composition of geologic materials. Tritium concentrations ranged from 0.2 (detection limit) to 1,380 nanocuries per liter. The methods of constructing, installing, and sampling with lysimeters were evaluated to ensure data reliability. These evaluations indicate that, with respect to most constituents, the samples retrieved from the lysimeters accurately represented pore-water chemistry.

Open-File Report↗
USGS1992

Effects of low-level radioactive-waste disposal on water chemistry in the unsaturated zone at a site near Sheffield, Illinois, 1982-84

A 1982-84 field study defined the chemistry of water collected from the unsaturated zone at a low-level radioactive-waste disposal site near Sheffield, Bureau County, Ill. Chemical data were evaluated to determine the principal, naturally occurring geochemical reactions in the unsaturated zone and to evaluate waste-induced effects on pore-water chemistry. Samples of precipitation, unsaturated-zone pore water, and saturated-zone water were analyzed for specific conductance, pH, alkalinity, major cations and anions, dissolved organic carbon, gross alpha and beta radiation, and tritium. Little change in concentration of most major constituents in the unsaturated-zone water was observed with respect to depth or distance from disposal trenches. Tritium and dissolved organic carbon concentrations were, however, dependent on proximity to trenches. The primary reactions, both on-site and off-site, were carbonate and clay dissolution, cation exchange, and the oxidation of pyrite. The major difference between on-site and off-site inorganic water chemistry resulted from the removal of the Roxana Silt and the Radnor Till Member of the Glasford Formation from on-site. Off-site, the Roxana Silt contributed substantial quantities of sodium to solution from montmorillonite dissolution and associated cation-exchange reactions. The Radnor Till Member provided exchange surfaces for magnesium. Precipitation at the site had an ionic composition of calcium zinc sulfate and an average pH of 4.6. Within 0.3 meter of the land surface, infiltrating rainwater or snowmelt changed to an ionic composition of calcium sulfate off-site and calcium bicarbonate on-site and had an average pH of 7.9; below that depth, pH averaged 7.5 and the ionic composition generally was calcium magnesium bicarbonate. Alkalinity and specific conductance differed primarily according to composition of geologic materials. Tritium concentrations ranged from 0.2 (detection limit) to 1,380 nanocuries per liter. The methods of constructing, installing, and sampling with lysimeters were evaluated to ensure data reliability. These evaluations indicate that, with respect to most constituents, the samples retrieved from the lysimeters accurately represented pore-water chemistry.

Illinois↗