USGS Science⌕ Search

Geology topics

Gordon W. Rattray

Publications and source records attributed to Gordon W. Rattray.

24 records · Page 2Linked to original sources

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

Open-File Report↗

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↗

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↗

Radiochemical and chemical constituents in water from selected wells and springs from the southern boundary of the Idaho National Laboratory to the Hagerman Area, Idaho, 2003

The U.S. Geological Survey and the Idaho Department of Water Resources, in cooperation with the U.S. Department of Energy, sampled water from 14 sites as part of an ongoing study to monitor the water quality of the eastern Snake River Plain aquifer between the southern boundary of the Idaho National Laboratory (INL) and the Burley-Twin Falls-Hagerman area. The State of Idaho, Department of Environmental Quality, Division of INL Oversight and Radiation Control cosampled with the U.S. Geological Survey and the Idaho Department of Water Resources and their analytical results are included in this report. The samples were collected from four domestic wells, two dairy wells, two springs, four irrigation wells, one observation well, and one stock well and analyzed for selected radiochemical and chemical constituents. Two quality-assurance samples, sequential replicates, also were collected and analyzed. None of the concentrations of radiochemical or organic-chemical constituents exceeded the maximum contaminant levels for drinking water established by the U.S. Environmental Protection Agency. However, the concentration of one inorganic-chemical constituent, nitrate (as nitrogen), in water from site MV-43 was 20 milligrams per liter which exceeded the maximum contaminant level for that constituent. Of the radiochemical and chemical concentrations analyzed for in the replicate-sample pairs, 267 of the 270 pairs (with 95 percent confidence) were statistically equivalent.

Open-File Report↗

Radiochemical and Chemical Constituents in Water from Selected Wells and Springs from the Southern Boundary of the Idaho National Engineering and Environmental Laboratory to the Hagerman Area, Idaho, 2002

The U.S. Geological Survey, Idaho Department of Water Resources, and the State of Idaho INEEL Oversight Program, in cooperation with the U.S. Department of Energy, sampled water from 17 sites as part of the sixth round of a long-term project to monitor water quality of the eastern Snake River Plain aquifer from the southern boundary of the Idaho National Engineering and Environmental Laboratory to the Hagerman area. The samples were collected from eight irrigation wells, three domestic wells, one stock well, one dairy well, one commercial well, one observation well, and two springs and analyzed for selected radiochemical and chemical constituents. One quality-assurance sample, a sequential replicate, also was collected and analyzed. Many of the radionuclide and inorganic-constituent concentrations were greater than the reporting levels and most of the organic-constituent concentrations were less than the reporting levels. However, none of the reported radiochemical- or chemical-constituent concentrations exceeded the maximum contaminant levels for drinking water established by the U.S. Environmental Protection Agency. Statistical evaluation of the replicate sample pair indicated that, with 95 percent confidence, 132 of the 135 constituent concentrations of the replicate pair were equivalent.

Idaho↗