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

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

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↗
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↗
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

Solution of the advection-dispersion equation by a finite-volume eulerian-lagrangian local adjoint method

A finite-volume Eulerian-Lagrangian local adjoint method for solution of the advection-dispersion equation is developed and discussed. The method is mass conservative and can solve advection-dominated ground-water solute-transport problems accurately and efficiently. An integrated finite-difference approach is used in the method. A key component of the method is that the integral representing the mass-storage term is evaluated numerically at the current time level. Integration points, and the mass associated with these points, are then forward tracked up to the next time level. The number of integration points required to reach a specified level of accuracy is problem dependent and increases as the sharpness of the simulated solute front increases. Integration points are generally equally spaced within each grid cell. For problems involving variable coefficients it has been found to be advantageous to include additional integration points at strategic locations in each well. These locations are determined by backtracking. Forward tracking of boundary fluxes by the method alleviates problems that are encountered in the backtracking approaches of most characteristic methods. A test problem is used to illustrate that the new method offers substantial advantages over other numerical methods for a wide range of problems.

Conference Paper↗
USGS1989

Seepage through a hazardous-waste trench cover

Water movement through a waste-trench cover under natural conditions at a low-level radioactive waste disposal site in northwestern Illinois was studied from July 1982 to June 1984, using tensiometers, a moisture probe, and meteorological instruments. Four methods were used to estimate seepage: the Darcy, zero-flux plane, surface-based water-budget, and groundwater-based water-budget methods. Annual seepage estimates ranged from 48 to 216 mm (5-23% of total precipitation), with most seepage occurring in spring. The Darcy method, although limited in accuracy by uncertainty in hydraulic conductivity, was capable of discretizing seepage in space and time and indicated that seepage varied by almost an order of magnitude across the width of the trench. Lowest seepage rates occurred near the center of the cover, where seepage was gradual. Highest rates occurred along the edge of the cover, where seepage was highly episodic, with 84% of the total there being traced to wetting fronts from 28 individual storms. Limitations of the zero-flux-plane method were severe enough for the method to be judged inappropriate for use in this study.Water movement through a waste-trench cover under natural conditions at a low-level radioactive waste disposal site in northwestern Illinois was studied from July 1982 to June 1984, using tensiometers, a moisture probe, and meteorological instruments. Four methods were used to estimate seepage: the Darcy, zero-flux plane, surface-based water-budget, and groundwater-based water-budget methods. Annual seepage estimates ranged from 48 to 216mm (5-23% of total precipitation), with most seepage occurring in spring. The Darcy method, although limited in accuracy by uncertainty in hydraulic conductivity, was capable of discretizing seepage in space and time and indicated that seepage varied by almost an order of magnitude across the width of the trench. Lowest seepage rates occurred near the center of the cover, where seepage was gradual. Highest rates occurred along the edge of the cover, where seepage was highly episodic, with 84% of the total there being traced to wetting fronts from 28 individual storms. Limitations of the zero-flux-plane method were severe enough for the method to be judged inappropriate for use in this study.

Illinois↗
USGS1989

A numerical solution for the diffusion equation in hydrogeologic systems

The documentation of a computer code for the numerical solution of the linear diffusion equation in one or two dimensions in Cartesian or cylindrical coordinates is presented. Applications of the program include molecular diffusion, heat conduction, and fluid flow in confined systems. The flow media may be anisotropic and heterogeneous. The model is formulated by replacing the continuous linear diffusion equation by discrete finite-difference approximations at each node in a block-centered grid. The resulting matrix equation is solved by the method of preconditioned conjugate gradients. The conjugate gradient method does not require the estimation of iteration parameters and is guaranteed convergent in the absence of rounding error. The matrixes are preconditioned to decrease the steps to convergence. The model allows the specification of any number of boundary conditions for any number of stress periods, and the output of a summary table for selected nodes showing flux and the concentration of the flux quantity for each time step. The model is written in a modular format for ease of modification. The model was verified by comparison of numerical and analytical solutions for cases of molecular diffusion, two-dimensional heat transfer, and axisymmetric radial saturated fluid flow. Application of the model to a hypothetical two-dimensional field situation of gas diffusion in the unsaturated zone is demonstrated. The input and output files are included as a check on program installation. The definition of variables, input requirements, flow chart, and program listing are included in the attachments. (USGS)

Water-Resources Investigations Report↗
USGS1989

Evapotranspiration and microclimate at a low-level radioactive-waste disposal site in northwestern Illinois

From July 1982 through June 1984, a study was made of the evapotranspiration and microclimate at a low-level radioactive-waste disposal site near Sheffield, Bureau County, Illinois. Vegetation at the site consists of mixed pasture grasses, primarily awnless brome (Bromus inermis) and red clover (Trifoleum pratense). Three methods were used to estimate evapotranspiration: (1) an energy budget with the Bowen ratio, (2) an aerodynamic profile, and (3) a soil-based water budget. For the aerodynamic-profile method, sensible-heat flux was estimated by a profile equation and evapotranspiration was then calculated as the residual in the energy-balance equation. Estimates by the energy-budget and aerodynamic-profile methods were computed from hourly data and then summed by days and months. Yearly estimates (for March through November) by these methods were in close agreement: 648 and 626 millimeters, respectively. Daily estimates reach a maximum of about 6 millimeters. The water-budget method produced only monthly estimates based on weekly or biweekly soil-moisture content measurements. The yearly evapotranspiration estimated by this method (which actually included only the months of April through October) was 655 millimeters. The March-through-November average for the three methods of 657 millimeters was equivalent to 70 percent of total precipitation. Continuous measurements were made of incoming and reflected shortwave radiation, incoming and emitted longwave radiation, net radiation, soil-heat flux, soil temperature, horizontal windspeed, and wet- and dry-bulb air temperature. Windspeed and air temperature were measured at heights of 0.5 and 2.0 meters (and also at 1.0 meter after September 1983). Soilmoisture content of the soil zone was measured with a gamma-attenuation gage. Annual precipitation (938 millimeters) and average temperature (10.8 degrees Celsius) at the Sheffield site were virtually identical to long-term averages from nearby National Weather Service stations. Solar radiation averaged 65 percent of that normally expected under clear skies. Net radiation averaged 70.1 watts per square meter and was highest in July and negative during some winter months. Wind direction varied but was predominately south-southeasterly. Wind speed at the 2-meter height averaged 3.5 meters per second and was slightly higher in winter months than the rest of the year. The amount of water stored within the soil zone was greatest in early spring and least in late summer. Seasonal and diurnal trends of evapotranspiration rates mirrored those of net radiation; July was usually the month with the highest evapotranspiration rate. The ratio of sensible- to latentheat fluxes (commonly called the Bowen ratio) for the 2-year study period was 0.38, as averaged from the three methods. Monthly Bowen ratios fluctuated somewhat but averaged about 0.35 for late spring through summer. In fall, the ratio declined to zero or to slightly negative values. When the ratio was negative, the latent-heat flux was slightly greater than the net radiation because of additional energy supplied by' the cooling soil and air. Evapotranspiration calculated by the three methods averaged 75 percent of potential evapotranspiration, as estimated by the Penman equation. There was no apparent seasonal trend in the relation between actual and potential evapotranspiration rates.

Water Supply Paper↗
USGS1988

Tritium migration from a low-level radioactive-waste disposal site near Chicago, Illinois

This paper describes the results of a study to determine the geologic and hydrologic factors that control migration of tritium from a closed, low-level radioactive-waste disposal site. The disposal site, which operated from 1943 to mid1949, contains waste generated by research activities at the world's first nuclear reactors. Tritium has migrated horizontally at least 1,300 feet northward in glacial drift and more than 650 feet in the underlying dolomite. Thin, gently sloping sand layers in an otherwise clayey glacial drift are major conduits for ground-water flow and tritium migration in a perched zone beneath the disposal site. Tritium concentrations in the drift beneath the disposal site exceed 100,000 nanocuries per liter. Regional horizontal joints in the dolomite are enlarged by solution and are the major conduits for ground-water flow and tritium migration in the dolomite. A weathered zone at the top of the dolomite also is a pathway for tritium migration. The maximum measured tritium concentration in the dolomite is 29.4 nanocuries per liter. Fluctuations of tritium concentration in the dolomite are the result of dilution by seasonal recharge from the drift.

Water Supply Paper↗
USGS1987

Evapotranspiration and microclimate at a low-level radioactive-waste disposal site in northwestern Illinois

From July 1982 through June 1984, a study was made of the microclimate and evapotranspiration at a low-level radioactive-waste disposal site near Sheffield, Bureau County, Illinois. Vegetation at the site consists of mixed pasture grasses, primarily brome (Bromus inermis) and red clover (Trifoleum pratense). Three methods were used to estimate evapotranspiration: (1) an energy-budget with the Bowen ratio, (2) an aerodynamic-profile, and (3) a soil-based water-budget. For the aerodynamic-profile method, sensible-heat flux was estimated by a profile equation and evapotranspiration was then calculated as the residual in the energy-balance equation. Estimates by the energy-budget and aerodynamic-profile methods were computed from hourly data, then summed by days and months. Yearly estimates for March through November, by these methods, were quite close--648 and 626 millimeters, respectively. Daily estimates range up to a maximum of about 6 millimeters. The water-budget method produced only monthly estimates based on weekly or biweekly soil-moisture content measurements. The yearly evapotranspiration estimated by this method (which actually included only the months of April through October) was 655 millimeters. The March-through-November average for the three methods of 657 millimeters was equivalent to 70 percent of precipitation. Continuous measurements were made of incoming and reflected shortwave radiation, incoming and emitted longwave radiation, net radiation, soil-heat flux, soil temperature, horizontal windspeed, and wet- and dry-bulb air temperature. Windspeed and air temperature were measured at heights of 0.5 and 2.0 meters (and also at 1.0 meter after September 1983). Soil-moisture content of the soil zone was measured with a gamma-attenuation gage. Annual precipitation (938 millimeters) and average temperature (10.8 degrees Celsius) were virtually identical to long-term averages from nearby National Weather Service stations. Solar radiation averaged 65 percent of that normally expected under clear skies. Net radiation averaged 70.1 watts per square meter and was highest in July and negative during some winter months. Wind direction varied but was predominately out of the south-southeast. Wind speed at the 2-meter height averaged 3.5 meters per second and was slightly higher in winter months than the rest of the year. The amount of water stored within the soil zone was greatest in early spring and least in late summer. Seasonal and diurnal trends in evapotranspiration rates mirrored those in net radiation; July was usually the month with the highest rate. The ratio of sensible- to latent-heat fluxes (commonly called the Bowen ratio) for the 2-year period was 0.38, as averaged from the three methods. Monthly Bowen ratios fluctuated somewhat but averaged about 0.35 for late spring through summer. In fall, the ratio declined to zero or to slightly negative values. When the ratio was negative, the latent-heat flux was slightly greater than the net radiation because of additional energy supplied by the cooling soil and air. Evapotranspiration calculated by the three methods averaged 75 percent of potential evapotranspiration, as estimated by the Penman equation. There was no apparent seasonal trend in the relation between actual and potential evapotranspiration rates.

Illinois↗
USGS1986

Concepts and data-collection techniques used in a study of the unsaturated zone at a low-level radioactive-waste disposal site near Sheffield, Illinois

A study of water and radionuclide movement through the unsaturated zone is being conducted at the low level radioactive waste disposal site near Sheffield, Illinois. Included in the study are detailed investigations of evapotranspiration, movement of water through waste trench covers, and movement of water and radionuclides (dissolved and gaseous) from the trenches. An energy balance/Bowen ratio approach is used to determine evapotranspiration. Precipitation, net radiation, soil-heat flux, air temperature and water vapor content gradients, wind speed, and wind direction are measured. Soil water tension is measured with tensiometers which are connected to pressure transducers. Meteorological sensors and tensiometers which are connected to pressure transducers. Meteorological sensors and tensiometers are monitored with automatic data loggers. Soil moisture contents are measured through small-diameter access tubes with neutron and gamma-ray attenuation gages. Data beneath the trenches are obtained through a 130-meter-long tunnel which extends under four of the trenches. Water samples are obtained with suction lysimeters, and samples of the geologic material are obtained with core tubes. These samples are analyzed for radiometric and inorganic chemistry. Gas samples are obtained from gas piezometers and analyzed for partial pressures of major constituents, Radon-222, tritiated water vapor, and carbon-14 dioxide. (USGS)

Water-Resources Investigations Report↗