Search USGS⌕ Search

USGS · 70073530

Evolution of the conceptual model of unsaturated zone hydrology at Yucca Mountain, Nevada

Abstract

Yucca Mountain is an arid site proposed for consideration as the United States’ first underground high-level radioactive waste repository. Low rainfall (approximately 170 mm/yr) and a thick unsaturated zone (500–1000 m) are important physical attributes of the site because the quantity of water likely to reach the waste and the paths and rates of movement of the water to the saturated zone under future climates would be major factors in controlling the concentrations and times of arrival of radionuclides at the surrounding accessible environment. The framework for understanding the hydrologic processes that occur at this site and that control how quickly water will penetrate through the unsaturated zone to the water table has evolved during the past 15 yr. Early conceptual models assumed that very small volumes of water infiltrated into the bedrock (0.5–4.5 mm/yr, or 2–3 percent of rainfall), that much of the infiltrated water flowed laterally within the upper nonwelded units because of capillary barrier effects, and that the remaining water flowed down faults with a small amount flowing through the matrix of the lower welded, fractured rocks. It was believed that the matrix had to be saturated for fractures to flow. However, accumulating evidence indicated that infiltration rates were higher than initially estimated, such as infiltration modeling based on neutron borehole data, bomb-pulse isotopes deep in the mountain, perched water analyses and thermal analyses. Mechanisms supporting lateral diversion did not apply at these higher fluxes, and the flux calculated in the lower welded unit exceeded the conductivity of the matrix, implying vertical flow of water in the high permeability fractures of the potential repository host rock, and disequilibrium between matrix and fracture water potentials. The development of numerical modeling methods and parameter values evolved concurrently with the conceptual model in order to account for the observed field data, particularly fracture flow deep in the unsaturated zone. This paper presents the history of the evolution of conceptual models of hydrology and numerical models of unsaturated zone flow at Yucca Mountain, Nevada (Flint, A.L., Flint, L.E., Kwicklis, E.M., Bodvarsson, G.S., Fabryka-Martin, J.M., 2001. Hydrology of Yucca Mountain. Reviews of Geophysics in press). This retrospective is the basis for recommendations for optimizing the efficiency with which a viable and robust conceptual model can be developed for a complex site.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 35.4976° to 37.501° latitude; -117.2379° to -115.4938° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Alan L. Flint, Lorraine E. Flint, Gudmundur S. Bodvarsson, Edward M. Kwicklis, June Fabryka-Martin. 2001. Evolution of the conceptual model of unsaturated zone hydrology at Yucca Mountain, Nevada. https://doi.org/10.1016/s0022-1694(01)00358-4

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Inverse modeling for post-wildfire hydrologic process understanding over snow-dominated watersheds

Wildfire-induced changes to vegetation and soil properties can alter snowpack dynamics and hydrologic processes in ways that complicate predictions of post-wildfire water availability. To better understand these effects, we coupled the Parameter ESTimation in C++ (PEST++) software with the Weather Research & Forecasting Hydrologic model (WRF-Hydro) across three burned, snow-dominated watersheds in the western United States. We incorporated three burn-severity-based land cover categories into WRF-Hydro and used an iterative Ensemble Smoother (iES) within PEST++ to estimate model parameters controlling snow accumulation and melt, infiltration, and streamflow generation. This framework supported three complementary objectives: (1) assessing parameter sensitivities under pre- and post-wildfire conditions, (2) comparing pre- and post-wildfire parameter ensembles to identify wildfire-driven shifts in hydrologic processes, and (3) using counterfactual simulations to isolate the hydrologic response to wildfire from interannual climate variability. Our findings indicate that model parameter sensitivity varied by domain and shifted following wildfire. Maximum snow albedo and snow albedo decay were most influential in the higher-elevation East Troublesome domain, while the rain-snow partitioning temperature threshold dominated in the lower-elevation Chippy Creek and Las Conchas domains. Following wildfire, streamflow sensitivity shifted away from snowpack parameters toward soil hydraulic parameters, consistent with wildfire-driven reductions in evapotranspiration and increases in soil moisture. Incorporating wildfire effects led to complex changes in simulated snowpack accumulation and streamflow timing, generally leading to earlier snowmelt freshet initiation and cessation and reduced late season flows. Simulations were more sensitive to interannual climate variations than wildfire land cover changes. These results highlight both the potential and challenges of ensemble-based inverse modeling for understanding post-wildfire hydrologic change in snow-dominated watersheds.

Colorado, Montana, New Mexico↗

Reservoir releases and land cover interact to drive event-scale nitrate export in a large agricultural basin

Understanding the drivers of nitrate export in rivers is critical for developing effective nutrient management strategies. However, few studies have explored event-scale drivers of export in large river basins with human modifications like reservoirs. Here, we analyzed nitrate concentration-discharge (C-Q) relationships from 215 events at the outlet of the Kansas River Basin, USA (155,690 km 2 ) from 2014 to 2022 to (1) characterize event-scale nitrate export behaviors in a large agricultural basin, and (2) determine how different event characteristics are linked to these behaviors. We found that C-Q behaviors varied greatly, with 60% of events producing nitrate enrichment (n = 130) and 40% of events producing nitrate dilution (n = 85). These behaviors were correlated with complex spatial interlinkages between climate and land cover: across most of the basin, nitrate enrichment was correlated with drier antecedent conditions, but in wetter areas with higher proportions of urban and forested land cover, enrichment was more strongly correlated with precipitation magnitude/intensity. This difference in hydroclimatic controls on nitrate export might be related to differential distributions of nitrate sources within these land covers. Upstream of major reservoirs, however, neither variable was strongly correlated with C-Q behavior, suggesting that nitrate attenuation within reservoirs decouples event-scale concentration signals in upstream waters from those downstream. Reservoir outflows had variable impacts on C-Q behavior, reflecting reservoir-specific variations in nitrate attenuation efficiency. Together, these results identify specific complex interactions between hydroclimate, land cover, reservoir positioning, and individual reservoir properties that control event-scale nitrate export from large basins.

Colorado, Kansas, Nebraska↗

Deep learning error post-processing improves stochastic watershed modeling

Hydrologic extremes, including floods and droughts, pose substantial societal risks that are expected to intensify with climate change. Deterministic watershed models (DWMs) remain a mainstay for modeling these extremes, but lack explicit representation of uncertainty, limiting their utility for risk-informed planning. Stochastic watershed models (SWMs) address this limitation by generating ensembles of streamflow via models of observed DWM residuals. However, most SWMs struggle with the complex dependence between DWM residuals and the underlying hydrologic state, which can complicate stochastic simulations under nonstationary climates. Deep learning (DL) models, whether used as standalone models or post-processors for process-based DWMs, offer a pathway to address this challenge by reducing conditional dependence. In this study, we evaluate SWMs applied to seven models: three process-based models (PRMS, Hymod, and HBV), their hybrid process-DL counterparts, and a pure DL DWM, focusing on daily simulations and extremes under both historical conditions and synthetic climate change scenarios. Results for a case study watershed in Massachusetts show that SWMs applied to hybrid or pure DL DWMs consistently outperform those applied to process-based DWMs. However, an SWM applied to the pure DL model exhibits weaknesses at low flows for this study basin, underscoring the value of hybrid approaches. Extending the analysis across 73 additional basins demonstrates that these improvements are robust and generalizable statewide. This work highlights the potential of a DL-enhanced stochastic watershed modeling framework to advance hydrologic risk prediction under changing climate conditions, offering a scalable methodology for integrating uncertainty into watershed modeling for long-term planning.

Journal of Hydrology↗