Edward Michael Kwicklis

Temperature-Profile Methods for Estimating Percolation Rates in Arid Environments

Percolation rates are estimated using vertical temperature profiles from sequentially deeper vadose environments, progressing from sediments beneath stream channels, to expansive basin-fill materials, and finally to deep fractured bedrock underlying mountainous terrain. Beneath stream channels, vertical temperature profiles vary over time in response to downward heat transport, which is generally controlled by conductive heat transport during dry periods, or by advective transport during channel infiltration. During periods of stream-channel infiltration, two relatively simple approaches are possible: a heat-pulse technique, or a heat and liquid-water transport simulation code. Focused percolation rates beneath stream channels are examined for perennial, seasonal, and ephemeral channels in central New Mexico, with estimated percolation rates ranging from 100 to 2100 mm d−1. Deep within basin-fill and underlying mountainous terrain, vertical temperature gradients are dominated by the local geothermal gradient, which creates a profile with decreasing temperatures toward the surface. If simplifying assumptions are employed regarding stratigraphy and vapor fluxes, an analytical solution to the heat transport problem can be used to generate temperature profiles at specified percolation rates for comparison to the observed geothermal gradient. Comparisons to an observed temperature profile in the basin-fill sediments beneath Frenchman Flat, Nevada, yielded water fluxes near zero, with absolute values <10 mm yr−1. For the deep vadose environment beneath Yucca Mountain, Nevada, the complexities of stratigraphy and vapor movement are incorporated into a more elaborate heat and water transport model to compare simulated and observed temperature profiles for a pair of deep boreholes. Best matches resulted in a percolation rate near zero for one borehole and 11 mm yr−1 for the second borehole.

Hydrology of Yucca Mountain, Nevada

Yucca Mountain, located in southern Nevada in the Mojave Desert, is being considered as a geologic repository for high-level radioactive waste. Although the site is arid, previous studies indicate net infiltration rates of 5–10 mm yr−1 under current climate conditions. Unsaturated flow of water through the mountain generally is vertical and rapid through the fractures of the welded tuffs and slow through the matrix of the nonwelded tuffs. The vitric-zeolitic boundary of the nonwelded tuffs below the potential repository, where it exists, causes perching and substantial lateral flow that eventually flows through faults near the eastern edge of the potential repository and recharges the underlying groundwater system. Fast pathways are located where water flows relatively quickly through the unsaturated zone to the water table. For the bulk of the water a large part of the travel time from land surface to the potential repository horizon (∼300 m below land surface) is through the interlayered, low fracture density, nonwelded tuff where flow is predominately through the matrix. The unsaturated zone at Yucca Mountain is being modeled using a three-dimensional, dual-continuum numerical model to predict the results of measurements and observations in new boreholes and excavations. The interaction between experimentalists and modelers is providing confidence in the conceptual model and the numerical model and is providing researchers with the ability to plan further testing and to evaluate the usefulness or necessity of further data collection.

Preliminary development of the LBL/USGS three-dimensional site-scale model of Yucca Mountain, Nevada

A 3-D model of moisture flow within the unsaturated zone at Yucca Mountain is being developed at LBL in cooperation with USGS. This site-scale model covers an area of about 34 km{sup 2} and is bounded by major faults to the north, east, and west. The relatively coarse-grid model includes about 300 horizontal grid-blocks and 17 layers. Contour maps and isopach maps are presented defining different types of infiltration zones, and the spatial distribution of Tiva Canyon, Paintbrush, and Topopah Spring hydrogeological units. Matrix flow is approximated using the van Genuchten model, and the equivalent continuum approximation is used to account for fracture flow in the welded units. One-, two-, and three-dimensional simulations are conducted using the TOUGH2 computer program. Steady-state simulations are performed with various uniform and nonuniform infiltration rates; results are interpreted in terms of effect of fault characteristics on moisture flow distribution, and on the location and formation of preferential pathways.

The site-scale saturated zone flow model for Yucca Mountain: calibration of different conceptual models and their impact on flow paths

This paper presents several different conceptual models of the Large Hydraulic Gradient (LHG) region north of Yucca Mountain and describes the impact of those models on groundwater flow near the potential high-level repository site. The results are based on a numerical model of site-scale saturated zone beneath Yucca Mountain. This model is used for performance assessment predictions of radionuclide transport and to guide future data collection and modeling activities. The numerical model is calibrated by matching available water level measurements using parameter estimation techniques, along with more informal comparisons of the model to hydrologic and geochemical information. The model software (hydrologic simulation code FEHM and parameter estimation software PEST) and model setup allows for efficient calibration of multiple conceptual models. Until now, the Large Hydraulic Gradient has been simulated using a low-permeability, east-west oriented feature, even though direct evidence for this feature is lacking. In addition to this model, we investigate and calibrate three additional conceptual models of the Large Hydraulic Gradient, all of which are based on a presumed zone of hydrothermal chemical alteration north of Yucca Mountain. After examining the heads and permeabilities obtained from the calibrated models, we present particle pathways from the potential repository that record differences in the predicted groundwater flow regime. The results show that Large Hydraulic Gradient can be represented with the alternate conceptual models that include the hydrothermally altered zone. The predicted pathways are mildly sensitive to the choice of the conceptual model and more sensitive to the quality of calibration in the vicinity on the repository. These differences are most likely due to different degrees of fit of model to data, and do not represent important differences in hydrologic conditions for the different conceptual models.

Estimation of percolation flux from borehole temperature data at Yucca Mountain, Nevada

Temperature data from the unsaturated zone (UZ) at Yucca Mountain are analyzed to estimate percolation-flux rates and overall heat flux. A multilayer, one-dimensional analytical solution is presented for determining percolation flux from temperature data. Case studies have shown that the analytical solution agrees very well with results from the numerical code, TOUGH2. The results of the analysis yield percolation fluxes in the range from 0 to 20 mm/year for most of the deep boreholes. This range is in good agreement with the results of infiltration studies at Yucca Mountain. Percolation flux for the shallower boreholes, however, cannot be accurately determined from temperature data alone because large gas flow in the shallow system alters the temperature profiles. Percolation-flux estimates for boreholes located near or intersecting major faults are significantly higher than those for other boreholes. These estimates may be affected by gas flow in the faults.

The saturated zone at Yucca Mountain: An overview of the characterization and assessment of the saturated zone as a barrier to potential radionuclide migration

The US Department of Energy is pursuing Yucca Mountain, Nevada, for the development of a geologic repository for the disposal of spent nuclear fuel and high-level radioactive waste, if the repository is able to meet applicable radiation protection standards established by the US Nuclear Regulatory Commission and the US Environmental Protection Agency (EPA). Effective performance of such a repository would rely on a number of natural and engineered barriers to isolate radioactive waste from the accessible environment. Groundwater beneath Yucca Mountain is the primary medium through which most radionuclides might move away from the potential repository. The saturated zone (SZ) system is expected to act as a natural barrier to this possible movement of radionuclides both by delaying their transport and by reducing their concentration before they reach the accessible environment. Information obtained from Yucca Mountain Site Characterization Project activities is used to estimate groundwater flow rates through the site-scale SZ flow and transport model area and to constrain general conceptual models of groundwater flow in the site-scale area. The site-scale conceptual model is a synthesis of what is known about flow and transport processes at the scale required for total system performance assessment of the site. This knowledge builds on and is consistent with knowledge that has accumulated at the regional scale but is more detailed because more data are available at the site-scale level. The mathematical basis of the site-scale model and the associated numerical approaches are designed to assist in quantifying the uncertainty in the permeability of rocks in the geologic framework model and to represent accurately the flow and transport processes included in the site-scale conceptual model. Confidence in the results of the mathematical model was obtained by comparing calculated to observed hydraulic heads, estimated to measured permeabilities, and lateral flow rates calculated by the site-scale model to those calculated by the regional-scale flow model. In addition, it was confirmed that the flow paths leaving the region of the potential repository are consistent with those inferred from gradients of measured head and those independently inferred from water-chemistry data. The general approach of the site-scale SZ flow and transport model analysis is to calculate unit breakthrough curves for radionuclides at the interface between the SZ and the biosphere using the three-dimensional site-scale SZ flow and transport model. Uncertainties are explicitly incorporated into the site-scale SZ flow and transport abstractions through key parameters and conceptual models.

Upscaling retardation factor in hierarchical porous media with multimodal reactive mineral facies

Aquifer heterogeneity controls spatial and temporal variability of reactive transport parameters and has significant impacts on subsurface modeling of flow, transport, and remediation. Upscaling (or homogenization) is a process to replace a heterogeneous domain with a homogeneous one such that both reproduce the same response. To make reliable and accurate predictions of reactive transport for contaminant in chemically and physically heterogeneous porous media, subsurface reactive transport modeling needs upscaled parameters such as effective retardation factor to perform field-scale simulations. This paper develops a conceptual model of multimodal reactive mineral facies for upscaling reactive transport parameters of hierarchical heterogeneous porous media. Based on the conceptual model, covariance of hydraulic conductivity, sorption coefficient, flow velocity, retardation factor, and cross-covariance between flow velocity and retardation factor are derived from geostatistical characterizations of a three-dimensional unbounded aquifer system. Subsequently, using a Lagrangian approach the scale-dependent analytical expressions are derived to describe the scaling effect of effective retardation factors in temporal and spatial domains. When time and space scales become sufficiently large, the effective retardation factors approximate their composite arithmetic mean. Correlation between the hydraulic conductivity and the sorption coefficient can significantly affect the values of the effective retardation factor in temporal and spatial domains. When the temporal and spatial scales are relatively small, scaling effect of the effective retardation factors is relatively large. This study provides a practical methodology to develop effective transport parameters for field-scale modeling at which remediation and risk assessment is actually conducted. It does not only bridge the gap between bench-scale measurements to field-scale modeling, but also provide new insights into the influence of hierarchical mineral distribution on effective retardation factor.

Discrete Fracture Network Modeling and Simulation of Subsurface Transport for the Topopah Springs and Lava Flow Aquifers at Pahute Mesa, FY 15 Progress Report

This progress report for fiscal year 2015 (FY15) describes the development of discrete fracture network (DFN) models for Pahute Mesa. DFN models will be used to upscale parameters for simulations of subsurface flow and transport in fractured media in Pahute Mesa. The research focuses on modeling of groundwater flow and contaminant transport using DFNs generated according to fracture characteristics observed in the Topopah Spring Aquifer (TSA) and the Lava Flow Aquifer (LFA). This work will improve the representation of radionuclide transport processes in large-scale, regulatory-focused models with a view to reduce pessimistic bounding approximations and provide more realistic contaminant boundary calculations that can be used to describe the future extent of contaminated groundwater. Our goal is to refine a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, which tend to attenuate migration.