Steady-State and Transient Groundwater Flow and Advective Transport, Eastern Snake River Plain Aquifer, Idaho National Laboratory and Vicinity, Idaho
Abstract
Three-dimensional steady-state and transient models of groundwater flow and advective transport through the fractured basalts and interbedded sediments of the Eastern Snake River Plain (ESRP) aquifer were developed by the U.S. Geological Survey in cooperation with the U.S. Department of Energy. The model domain covers an area of 1,940 square miles that includes most 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 aquifer. Numerical models simulated 1980 steady-state conditions and transient flow for 1980-95. In the transient model, streamflow infiltration was the major stress. The models were calibrated using the parameter-estimation program incorporated in MODFLOW-2000. The steady-state model reasonably simulated the observed water-table altitude and gradients. Simulation of transient conditions reproduced changes in the flow system resulting from episodic infiltration from the Big Lost River. Analysis of simulations shows that flow is (1) dominantly horizontal through interflow zones in basalt, vertical anisotropy resulting from contrasts in hydraulic conductivity of different 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. Particle-tracking simulations were used to evaluate how simulated groundwater flow paths and travel times differ between the steady-state and transient flow models, and how well model-derived groundwater flow directions and velocities compare to independently-derived estimates. Particle tracking also was used to simulate the growth of tritium plumes originating at two INL facilities over a 16 year period under steady-state and transient flow conditions (1953-68). The shape, dimensions, and areal extent of these plumes were compared to a map of the plumes for 1968 from tritium releases beginning in 1952. Collectively, the particle-tracking simulations indicate that groundwater flow paths and velocities, based on uncalibrated estimates of porosity, are influenced by the dynamic character of the water table and the large contrasts in the hydraulic properties of the media, primarily hydraulic conductivity. Simulation results also indicate that temporal changes in the local hydraulic gradient can account for some of the observed dispersion of contaminants in the aquifer near the major sources of contamination and perhaps the majority of the observed dispersion several miles downgradient of these facilities. The distance downgradient of the facilities where simulated particle plumes were able to reasonably reproduce the 1968 tritium plume extended only to the boundary separating sediment-rich from sediment-poor aquifer layers about 4 mi downgradient of the contaminant source. Particle plumes simulated beyond this boundary were narrow and long, and did not reasonably reproduce the shape, dimensions, or position of the leading edge of the tritium plume; however, few data were available to characterize its true areal extent and shape.
- Publication:
-
AGU Fall Meeting Abstracts
- Pub Date:
- December 2009
- Bibcode:
- 2009AGUFM.H21H..05F
- Keywords:
-
- 1829 HYDROLOGY / Groundwater hydrology