USGS Science⌕ Search

Geology topics

R. W. Healy

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

At least 37 records · Page 2Linked to original sources

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↗

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 millimeters; mean annual recharge to the disposal trenches (July 1982 through June 1984) was estimated to be 107 millimeters. Average annual tritium flux below the study trenches was estimated to be 3.4 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 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 millimeters per day. A 120-meter-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 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 millimeters 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 concentrations of tritium with increasing distance from the trench bases in the sand unit indicates that water movement through the unit may be more pronounced than indicated by the pressure-head data. Localized, preferential flow paths may have gone undetected by the monitoring instruments used in the study.

Illinois↗

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↗

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↗

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↗

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↗

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↗

Documentation of computer program VS2D to solve the equations of fluid flow in variably saturated porous media

This report documents FORTRAN computer code for solving problems involving variably saturated single-phase flow in porous media. The flow equation is written with total hydraulic potential as the dependent variable, which allows straightforward treatment of both saturated and unsaturated conditions. The spatial derivatives in the flow equation are approximated by central differences, and time derivatives are approximated either by a fully implicit backward or by a centered-difference scheme. Nonlinear conductance and storage terms may be linearized using either an explicit method or an implicit Newton-Raphson method. Relative hydraulic conductivity is evaluated at cell boundaries by using either full upstream weighting, the arithmetic mean, or the geometric mean of values from adjacent cells. Nonlinear boundary conditions treated by the code include infiltration, evaporation, and seepage faces. Extraction by plant roots that is caused by atmospheric demand is included as a nonlinear sink term. These nonlinear boundary and sink terms are linearized implicitly. The code has been verified for several one-dimensional linear problems for which analytical solutions exist and against two nonlinear problems that have been simulated with other numerical models. A complete listing of data-entry requirements and data entry and results for three example problems are provided. (USGS)

Water-Resources Investigations Report↗

Simulation of trickle irrigation, an extension to the US Geological Survey's computer program VS2D

A method is presented for simulating water movement through unsaturated porous media in response to a constant rate of application from a surface source. Because the rate at which water can be absorbed by soil is limited, the water will pond; therefore the actual surface area over which the water is applied may change with time and in general will not be known beforehand. An iterative method is used to determine the size of this ponded area at any time. This method will be most useful for simulating trickling irrigation, but also may be of value for simulating movement of water is soils as the result of an accidental spill. The method is an extension to the finite difference computer program VS2D developed by the U.S. Geological Survey, which simulates water movement through variably saturated porous media. The simulated region can be a vertical, 2-dimensional cross section for treatment of a surface line source or an axially symmetric, 3-dimensional cylinder for a point source. Five test problems, obtained from the literature , are used to demonstrate the ability of the method to accurately match analytical and experimental results. (Author 's abstract)

Water-Resources Investigations Report↗

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↗

Hydrogeology, ground-water flow, and tritium movement at low-level radioactive-waste disposal site near Sheffield, Illinois

Groundwater flow and tritium movement are described at and near a low-level radioactive waste disposal site near Sheffield, Illinois. Flow in the shallow aquifer is confined to three basins that ultimately drain into a stripmine lake. Most of the flow from the site is through a buried, pebbly sandfilled channel. Remaining flow is toward alluvium of an existing stream. Conceptual flow models for the two largest basins are used to improve definition of flow velocity and direction. Flow velocities range from about 25 to 2,500 ft/yr. Tritium was found in all three basins. The most extensive migration of tritium is coincident with buried channel. Tritium concentrations ranged from detection level to more than 300 nanocuries/L.

Illinios↗

Hydrogeology of a low-level radioactive-waste disposal site near Sheffield, Illinois

The Sheffield low-level radioactive-waste facility is located on 20 acres of rolling terrain 3 miles southwest of Sheffield, Illinois. The shallow hydrogeologic system is composed of glacial sediments. Pennsylvania shale and mudstone bedrock isolate the regional aquifers below from the hydrogeologic system in the overlying glacial deposits. Pebbly sand underlies 67 percent of the site. Two ground-water flow paths were identified. The primary path conveys ground water from the site to the east through the pebbly-sand unit; a secondary path conveys ground water to the south and east through less permeable material. The pebbly-sand unit provides an underdrain that eliminates the risk of water rising into the trenches. Digital computer model results indicate that the pebbly-sand unit controls ground-water movement. Tritium found migrating in ground water in the southeast corner of the site travels approximately 25 feet per year. A group of water samples from wells which contained the highest tritium concentrations had specific conductivities, alkalinities, hardness, and chloride, sulfate, calcium, and magnesium contents higher than normal for local shallow ground water. (USGS)

Water-Resources Investigations Report↗

River mileages and drainage areas for Illinois streams—Volume 1, Illinois except Illinois River basin

River mileages are presented for points of interest on Illinois streams draining 10 square miles or more. Points of interest include bridges, dams, gaging stations, county lines, hydrologic unit boundaries, and major tributaries. Drainage areas are presented for selected sites, including total drainage area for any stream draining at least 100 square miles. The report is contained in two volumes. This volume (Volume I) includes Illinois except Illinois River basin.

Illinois↗

River mileages and drainage areas for Illinois streams—Volume 2, Illinois River basin

River mileages are presented for points of interest on Illinois streams draining 10 square miles or more. Points of interest include bridges, dams, gaging stations, county lines, hydrologic unit boundaries, and major tributaries. Drainage areas are presented for selected sites, including total drainage area for any stream draining at least 100 square miles. The report is contained in two volumes. This volume (Volume II) includes Illinois River basin.

Illinois↗

Chemical analyses of surface water in Illinois, 1975-77--Volume I, Des Plaines River basin and Lake Michigan

Samples of surface water were collected and analyzed by the Illinois Environmental Protection Agency. The results from water years 1975 to 1977 are presented in three volumes. The history of sampling and analytical methods used during that period are summarized. Stream discharge data from records of the U.S. Geological Survey are included for all sites where samples were collected at gaging stations or near enough that reliable discharge estimates could be made. Volume I includes the Des Plaines River basin and Lake Michigan.

Illinois↗

Chemical analyses of surface water in Illinois, 1975-77--Volume II, Illinois River basin and Mississippi River tributaries north of Illinois River basin

Samples of surface water were collected and analyzed by the Illinois Environmental Protection Agency. The results from water years 1975 to 1977 are presented in three volumes. The history of sampling and analytical methods used during that period are summarized. Stream discharge data from records of the U.S. Geological Survey are included for all sites where samples were collected at gaging stations or near enough that reliable discharge estimates could be made. Volume II includes the Illinois River basin and Mississippi River tributaries north of Illinois River basin.

Illinois↗

Chemical analyses of surface water in Illinois, 1975-77--Volume III, Ohio River tributaries and Mississippi River tributaries south of Illinois River basin

Samples of surface water were collected and analyzed by the Illinois Environmental Protection Agency. The results from water years 1975 to 1977 are presented in three volumes. The history of sampling and analytical methods used during that period are summarized. Stream discharge data from records of the U.S. Geological Survey are included for all sites where samples were collected at gaging stations or near enough that reliable discharge estimates could be made. Volume III includes Ohio River tributaries and Mississippi River tributaries south of Illinois River basin.

Illinois↗