Geology ReportsSearch

Geology topics

Joseph P. Rousseau

Publications and source records attributed to Joseph P. Rousseau.

10 recordsLinked to original sources

A comparison of U.S. Geological Survey three-dimensional model estimates of groundwater source areas and velocities to independently derived estimates, Idaho National Laboratory and vicinity, Idaho

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, evaluated a three-dimensional model of groundwater flow in the fractured basalts and interbedded sediments of the eastern Snake River Plain aquifer at and near the Idaho National Laboratory to determine if model-derived estimates of groundwater movement are consistent with (1) results from previous studies on water chemistry type, (2) the geochemical mixing at an example well, and (3) independently derived estimates of the average linear groundwater velocity. Simulated steady-state flow fields were analyzed using backward particle-tracking simulations that were based on a modified version of the particle tracking program MODPATH. Model results were compared to the 5-microgram-per-liter lithium contour interpreted to represent the transition from a water type that is primarily composed of tributary valley underflow and streamflow-infiltration recharge to a water type primarily composed of regional aquifer water. This comparison indicates several shortcomings in the way the model represents flow in the aquifer. The eastward movement of tributary valley underflow and streamflow-infiltration recharge is overestimated in the north-central part of the model area and underestimated in the central part of the model area. Model inconsistencies can be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones. Sources of water at well NPR-W01 were identified using backward particle tracking, and they were compared to the relative percentages of source water chemistry determined using geochemical mass balance and mixing models. The particle tracking results compare reasonably well with the chemistry results for groundwater derived from surface-water sources (-28 percent error), but overpredict the proportion of groundwater derived from regional aquifer water (108 percent error) and underpredict the proportion of groundwater derived from tributary valley underflow from the Little Lost River valley (-74 percent error). These large discrepancies may be attributed to large contrasts in hydraulic conductivity between hydrogeologic zones and (or) a short-circuiting of underflow from the Little Lost River valley to an area of high hydraulic conductivity. Independently derived estimates of the average groundwater velocity at 12 well locations within the upper 100 feet of the aquifer were compared to model-derived estimates. Agreement between velocity estimates was good at wells with travel paths located in areas of sediment-rich rock (root-mean-square error [RMSE] = 5.2 feet per day [ft/d]) and poor in areas of sediment-poor rock (RMSE = 26.2 ft/d); simulated velocities in sediment-poor rock were 2.5 to 4.5 times larger than independently derived estimates at wells USGS 1 (less than 14 ft/d) and USGS 100 (less than 21 ft/d). The models overprediction of groundwater velocities in sediment-poor rock may be attributed to large contrasts in hydraulic conductivity and a very large, model-wide estimate of vertical anisotropy (14,800).

Idaho

Steady-state and transient models of groundwater flow and advective transport, Eastern Snake River Plain aquifer, Idaho National Laboratory and vicinity, Idaho

Three-dimensional steady-state and transient models of groundwater flow and advective transport in the eastern Snake River Plain aquifer were developed by the U.S. Geological Survey in cooperation with the U.S. Department of Energy. The steady-state and transient flow models cover an area of 1,940 square miles that includes most of the 890 square miles of the Idaho National Laboratory (INL). A 50-year history of waste disposal at the INL has resulted in measurable concentrations of waste contaminants in the eastern Snake River Plain aquifer. Model results can be used in numerical simulations to evaluate the movement of contaminants in the aquifer. Saturated flow in the eastern Snake River Plain aquifer was simulated using the MODFLOW-2000 groundwater flow model. Steady-state flow was simulated to represent conditions in 1980 with average streamflow infiltration from 1966-80 for the Big Lost River, the major variable inflow to the system. The transient flow model simulates groundwater flow between 1980 and 1995, a period that included a 5-year wet cycle (1982-86) followed by an 8-year dry cycle (1987-94). Specified flows into or out of the active model grid define the conditions on all boundaries except the southwest (outflow) boundary, which is simulated with head-dependent flow. In the transient flow model, streamflow infiltration was the major stress, and was variable in time and location. The models were calibrated by adjusting aquifer hydraulic properties to match simulated and observed heads or head differences using the parameter-estimation program incorporated in MODFLOW-2000. Various summary, regression, and inferential statistics, in addition to comparisons of model properties and simulated head to measured properties and head, were used to evaluate the model calibration. Model parameters estimated for the steady-state calibration included hydraulic conductivity for seven of nine hydrogeologic zones and a global value of vertical anisotropy. Parameters estimated for the transient calibration included specific yield for five of the seven hydrogeologic zones. The zones represent five rock units and parts of four rock units with abundant interbedded sediment. All estimates of hydraulic conductivity were nearly within 2 orders of magnitude of the maximum expected value in a range that exceeds 6 orders of magnitude. The estimate of vertical anisotropy was larger than the maximum expected value. All estimates of specific yield and their confidence intervals were within the ranges of values expected for aquifers, the range of values for porosity of basalt, and other estimates of specific yield for basalt. The steady-state model reasonably simulated the observed water-table altitude, orientation, and gradients. Simulation of transient flow conditions accurately reproduced observed changes in the flow system resulting from episodic infiltration from the Big Lost River and facilitated understanding and visualization of the relative importance of historical differences in infiltration in time and space. As described in a conceptual model, the numerical model simulations demonstrate flow that is (1) dominantly horizontal through interflow zones in basalt and vertical anisotropy resulting from contrasts in hydraulic conductivity of various types of basalt and the interbedded sediments, (2) temporally variable due to streamflow infiltration from the Big Lost River, and (3) moving downward downgradient of the INL. The numerical models were reparameterized, recalibrated, and analyzed to evaluate alternative conceptualizations or implementations of the conceptual model. The analysis of the reparameterized models revealed that little improvement in the model could come from alternative descriptions of sediment content, simulated aquifer thickness, streamflow infiltration, and vertical head distribution on the downgradient boundary. Of the alternative estimates of flow to or from the aquifer, only a 20 percent decrease in

Scientific Investigations Report

Field Methods and Quality-Assurance Plan for Quality-of-Water Activities, U.S. Geological Survey, Idaho National Laboratory, Idaho

Water-quality activities conducted by the staff of the U.S. Geological Survey (USGS) Idaho National Laboratory (INL) Project Office coincide with the USGS mission of appraising the quantity and quality of the Nation's water resources. The activities are conducted in cooperation with the U.S. Department of Energy's (DOE) Idaho Operations Office. Results of the water-quality investigations are presented in various USGS publications or in refereed scientific journals. The results of the studies are highly regarded, and they are used with confidence by researchers, regulatory and managerial agencies, and interested civic groups. In its broadest sense, quality assurance refers to doing the job right the first time. It includes the functions of planning for products, review and acceptance of the products, and an audit designed to evaluate the system that produces the products. Quality control and quality assurance differ in that quality control ensures that things are done correctly given the 'state-of-the-art' technology, and quality assurance ensures that quality control is maintained within specified limits.

Open-File Report

Hydraulic Characteristics of Bedrock Constrictions and Evaluation of One- and Two-Dimensional Models of Flood Flow on the Big Lost River at the Idaho National Engineering and Environmental Laboratory, Idaho

A 1.9-mile reach of the Big Lost River, between the Idaho National Engineering and Environmental Laboratory (INEEL) diversion dam and the Pioneer diversion structures, was investigated to evaluate the effects of streambed erosion and bedrock constrictions on model predictions of water-surface elevations. Two one-dimensional (1-D) models, a fixed-bed surface-water flow model (HEC-RAS) and a movable-bed surface-water flow and sediment-transport model (HEC-6), were used to evaluate these effects. The results of these models were compared to the results of a two-dimensional (2-D) fixed-bed model [Transient Inundation 2-Dimensional (TRIM2D)] that had previously been used to predict water-surface elevations for peak flows with sufficient stage and stream power to erode floodplain terrain features (Holocene inset terraces referred to as BLR#6 and BLR#8) dated at 300 to 500 years old, and an unmodified Pleistocene surface (referred to as the saddle area) dated at 10,000 years old; and to extend the period of record at the Big Lost River streamflow-gaging station near Arco for flood-frequency analyses. The extended record was used to estimate the magnitude of the 100-year flood and the magnitude of floods with return periods as long as 10,000 years. In most cases, the fixed-bed TRIM2D model simulated higher water-surface elevations, shallower flow depths, higher flow velocities, and higher stream powers than the fixed-bed HEC-RAS and movable-bed HEC-6 models for the same peak flows. The HEC-RAS model required flow increases of 83 percent [100 to 183 cubic meters per second (m3/s)], and 45 percent (100 to 145 m3/s) to match TRIM2D simulations of water-surface elevations at two paleoindicator sites that were used to determine peak flows (100 m3/s) with an estimated return period of 300 to 500 years; and an increase of 13 percent (150 to 169 m3/s) to match TRIM2D water-surface elevations at the saddle area that was used to establish the peak flow (150 m3/s) of a paleoflood with a return period of 10,000 years. A field survey of the saddle area, however, indicated that the elevation of the lowest point on the saddle area was 1.2 feet higher than indicated on the 2-ft contour map that was used in the TRIM2D model. Because of this elevation discrepancy, HEC-RAS model simulations indicated that a peak flow of at least 210 m3/s would be needed to initiate flow across the 10,000-year old Pleistocene surface. HEC-6 modeling results indicated that to compensate for the effects of streambed scour, additional flow increases would be needed to match HEC-RAS and TRIM2D water-surface elevations along the upper and middle reaches of the river, and to compensate for sediment deposition, a slight decrease in flows would be needed to match HEC-RAS water-surface elevations along the lower reach of the river. Differences in simulated water-surface elevations between the TRIM2D and the HEC-RAS and HEC-6 models are attributed primarily to differences in topographic relief and to differences in the channel and floodplain geometries used in these models. Topographic differences were sufficiently large that it was not possible to isolate the effects of these differences on simulated water-surface elevations from those attributable to the effects of supercritical flow, streambed scour, and sediment deposition.

Scientific Investigations Report

A conceptual model of ground-water flow in the eastern Snake River Plain aquifer at the Idaho National Laboratory and vicinity with implications for contaminant transport

Ground-water flow in the west-central part of the eastern Snake River Plain aquifer is described in a conceptual model that will be used in numerical simulations to evaluate contaminant transport at the Idaho National Laboratory (INL) and vicinity. The model encompasses an area of 1,940 square miles (mi2) and includes most of the 890 mi2 of the INL. A 50-year history of waste disposal associated with research activities at the INL has resulted in measurable concentrations of waste contaminants in the aquifer. A thorough understanding of the fate and movement of these contaminants in the subsurface is needed by the U.S. Department of Energy to minimize the effect that contaminated ground water may have on the region and to plan effectively for remediation. Three hydrogeologic units were used to represent the complex stratigraphy of the aquifer in the model area. Collectively, these hydrogeologic units include at least 65 basalt-flow groups, 5 andesite-flow groups, and 61 sedimentary interbeds. Three rhyolite domes in the model area extend deep enough to penetrate the aquifer. The rhyolite domes are represented in the conceptual model as low permeability, vertical pluglike masses, and are not included as part of the three primary hydrogeologic units. Broad differences in lithology and large variations in hydraulic properties allowed the heterogeneous, anisotropic basalt-flow groups, andesite-flow groups, and sedimentary interbeds to be grouped into three hydrogeologic units that are conceptually homogeneous and anisotropic. Younger rocks, primarily thin, densely fractured basalt, compose hydrogeologic unit 1; younger rocks, primarily of massive, less densely fractured basalt, compose hydrogeologic unit 2; and intermediate-age rocks, primarily of slightly-to-moderately altered, fractured basalt, compose hydrogeologic unit 3. Differences in hydraulic properties among adjacent hydrogeologic units result in much of the large-scale heterogeneity and anisotropy of the aquifer in the model area, and differences in horizontal and vertical hydraulic conductivity in individual hydrogeologic units result in much of the small-scale heterogeneity and anisotropy of the aquifer in the model area. The inferred three-dimensional geometry of the aquifer in the model area is very irregular. Its thickness generally increases from north to south and from west to east and is greatest south of the INL. The interpreted distribution of older rocks that underlie the aquifer indicates large changes in saturated thickness across the model area. The boundaries of the model include physical and artificial boundaries, and ground-water flows across the boundaries may be temporally constant or variable and spatially uniform or nonuniform. Physical boundaries include the water-table boundary, base of the aquifer, and northwest mountain-front boundary. Artificial boundaries include the northeast boundary, southeast-flowline boundary, and southwest boundary. Water flows into the model area as (1) underflow (1,225 cubic feet per second (ft3/s)) from the regional aquifer (northeast boundary-constant and nonuniform), (2) underflow (695 ft3/s) from the tributary valleys and mountain fronts (northwest boundary-constant and nonuniform), (3) precipitation recharge (70 ft3/s) (constant and uniform), streamflow-infiltration recharge (95 ft3/s) (variable and nonuniform), wastewater return flows (6 ft3/s) (variable and nonuniform), and irrigation-infiltration recharge (24 ft3/s) (variable and nonuniform) across the water table (water-table boundary-variable and nonuniform), and (4) upward flow across the base of the aquifer (44 ft3/s) (uniform and constant). The southeast-flowline boundary is represented as a no-flow boundary. Water flows out of the model area as underflow (2,037 ft3/s) to the regional aquifer (southwest boundary-variable and nonuniform) and as ground-water withdrawals (45 ft3/s) (water table boundary-variable and nonuniform). Ground-water flow i

Scientific Investigations Report

Historical development of the U.S. Geological Survey hydrologic monitoring and investigative programs at the Idaho National Engineering and Environmental Laboratory, Idaho, 1949 to 2001

This report is a summary of the historical development, from 1949 to 2001, of the U.S. Geological Survey's (USGS) hydrologic monitoring and investigative programs at the Idaho National Engineering and Environmental Laboratory. The report covers the USGS's water-level monitoring program, water-quality sampling program, geophysical program, geologic framework program, drilling program, modeling program, surface-water program, and unsaturated-zone program. The report provides physical information about the wells and information about the frequencies of sampling and measurement. Summaries of USGS published reports with U.S. Department of Energy (DOE) report numbers also are provided in an appendix. This report was prepared by the USGS in cooperation with the DOE.

Open-File Report

Review of the transport of selected radionuclides in the interim risk assessment for the Radioactive Waste Management Complex, Waste Area Group 7 Operable Unit 7-13/14, Idaho National Engineering and Environmental Laboratory, Idaho

The U.S. Department of Energy (DOE) requested that the U.S. Geological Survey conduct an independent technical review of the Interim Risk Assessment (IRA) and Contaminant Screening for the Waste Area Group 7 (WAG-7) Remedial Investigation, the draft Addendum to the Work Plan for Operable Unit 7-13/14 WAG-7 comprehensive Remedial Investigation and Feasibility Study (RI/FS), and supporting documents that were prepared by Lockheed Martin Idaho Technologies, Inc. The purpose of the technical review was to assess the data and geotechnical approaches that were used to estimate future risks associated with the release of the actinides americium, uranium, neptunium, and plutonium to the Snake River Plain aquifer from wastes buried in pits and trenches at the Subsurface Disposal Area (SDA). The SDA is located at the Radioactive Waste Management Complex in southeastern Idaho within the boundaries of the Idaho National Engineering and Environmental Laboratory. Radionuclides have been buried in pits and trenches at the SDA since 1957 and 1952, respectively. Burial of transuranic wastes was discontinued in 1982. The five specific tasks associated with this review were defined in a ?Proposed Scope of Work? prepared by the DOE, and a follow-up workshop held in June 1998. The specific tasks were (1) to review the radionuclide sampling data to determine how reliable and significant are the reported radionuclide detections and how reliable is the ongoing sampling program, (2) to assess the physical and chemical processes that logically can be invoked to explain true detections, (3) to determine if distribution coefficients that were used in the IRA are reliable and if they have been applied properly, (4) to determine if transport model predictions are technically sound, and (5) to identify issues needing resolution to determine technical adequacy of the risk assessment analysis, and what additional work is required to resolve those issues.

Idaho

Use of submersible pressure transducers in water-resources investigations

Submersible pressure transducers, developed in the early 1960s, have made the collection of water-level and pressure data much more convenient than former methods. Submersible pressure transducers, when combined with electronic data recorders have made it possible to collect continuous or nearly continuous water-level or pressure data from wells, piezometers, soil-moisture tensiometers, and surface water gages. These more frequent measurements have led to an improved understanding of the hydraulic processes in streams, soils, and aquifers. This manual describes the operational theory behind submersible pressure transducers and provides information about their use in hydrologic investigations conducted by the U.S. Geological Survey.

Techniques of Water-Resources Investigations

Estimating the magnitude of the 100-year peak flow in the Big Lost River at the Idaho National Engineering and Environmental Laboratory, Idaho

Accurate estimates of peak flows in the Big Lost River at the Idaho National Engineering and Environmental Laboratory (INEEL) are needed to assist planners and managers with evaluating possible effects of flooding on facilities at the INEEL. A large difference of 4,350 cubic feet per second (ft 3 /s) between two previous estimates of the magnitude of the 100-year peak flow in the Big Lost River near the western boundary of the INEEL prompted the present study. Regression models that compared annual peak flows and attenuation of annual peak flows between successive gaging stations for the same flow event were used to estimate the magnitude of the 100-year peak flow in the Big Lost River. The 100-year peak flow of 4,790 ft 3 /s at the Howell Ranch gaging station was used as the starting point for this analysis. This estimate was determined by using a three-parameter log-Pearson Type III distribution as outlined in “Guidelines for Determining Flood Flow Frequency” (Bulletin 17B by the Interagency Advisory Committee on Water Data). The regression models indicated that, in the reach of the Big Lost River between Howell Ranch and Mackay Reservoir, downstream peak flows are lower than upstream peak flows. Peak-flow attenuation values for this reach of the river decreased nonlinearly as the magnitude of the peak flow increased. Extrapolation of the trend resulted in an attenuation estimate of 13 percent for this reach relative to the 100-year peak flow at the Howell Ranch gaging station. In the lower reach of the Big Lost River between Mackay Reservoir and Arco, downstream peak flows are also lower than upstream peak flows. However, in contrast to the upper reach, peak-flow attenuation values decreased linearly as the magnitude of the peak flow increased. Extrapolation of the data indicated that peak-flow attenuations in this reach of the river approach zero for flows approaching the 100-year peak-flow estimate immediately upstream and downstream from Mackay Reservoir. A regression model of annual maximum daily mean flows between Arco and the INEEL diversion dam indicated that the attenuation values in this reach of the river are nearly the same for all flows of record. Extrapolation of the linear regression of these values resulted in an attenuation estimate of 10 percent. Seepage measurements made during 1951–53 also resulted in a loss estimate of approximately 10 percent. This attenuation value, combined with the values from analyses of the upstream reaches, resulted in an estimate of the 100-year peak flow for the Big Lost River immediately upstream from the INEEL diversion dam of 3,750 ft 3 /s; upper and lower 95-percent confidence limits were 6,250 ft 3 /s and 1,300 ft 3 /s, respectively. Localized rainfall, even of high intensity, is not likely to produce large peak flows at the INEEL because of high loss rates (infiltration, bank storage, and channel storage) along much of the stream channel. The relatively short flow durations resulting from rainstorms historically have not provided sufficient volumes of water to satisfy local storage demands (bank and channel storage). Only after these storage demands are met do the loss rates decrease enough for significant peak flows to reach the INEEL site. An uncertain component of the present analysis is the effect of seismic activity on the 100-year peak-flow estimate. Analysis of the effect of the magnitude 7.3 Borah Peak earthquake in 1983 on normal flow conditions in the Big Lost River suggests that the joint occurrence of a large earthquake and a 100-year peak flow could significantly increase the magnitude of the peak flow at the INEEL.

Idaho

Kilometer-scale rapid transport of naphthalene sulfonate tracer in the unsaturated zone at the Idaho National Engineering and Environmental Laboratory

To investigate possible long-range flow paths through the interbedded basalts and sediments of a 200-m-thick unsaturated zone, we applied a chemical tracer to seasonally filled infiltration ponds on the Snake River Plain in Idaho. This site is near the Subsurface Disposal Area for radioactive and other hazardous waste at the Idaho National Engineering and Environmental Laboratory. Within 4 mo, we detected tracer in one of 13 sampled aquifer wells, and in eight of 11 sampled perched-water wells as far as 1.3 km away. These detections show that (i) low-permeability layers in the unsaturated zone divert some flow horizontally, but do not prevent rapid transport to the aquifer; (ii) horizontal convective transport rates within the unsaturated zone may exceed 14 m d −1 , perhaps through essentially saturated basalt fractures, tension cracks, lava tubes, or rubble zones; and (iii) some perched water beneath the Subsurface Disposal Area derives from episodic surface water more than 1 km away. Such rapid and far-reaching flow may be common throughout the Snake River Plain, and possibly occurs in other locations that have a geologically complex unsaturated zone and comparable sources of infiltrating water.

Idaho