Search

This paper presents a study of solute transport through ground water in the saturated zone and the resulting breakthrough curves (BTCs), using a field-scale numerical model that incorporates the processes of advection, dispersion, matrix diffusion in fractured volcanic formations, sorption, and colloid-facilitated transport. Such BTCs at compliance boundaries are often used as performance measures for a site. The example considered here is that of the saturated zone study prepared for the Yucca Mountain license application. The saturated zone at this site occurs partly in volcanic, fractured rock formations and partly in alluvial formations. This paper presents a description of the site and the ground water flow model, the development of the conceptual model of transport, model uncertainties, model validation, and the influence of uncertainty in input parameters on the downstream BTCs at the Yucca Mountain site., (Copyright (c) 2010 Elsevier B.V. All rights reserved.)

Two papers recently published in the Journal of Contaminant Hydrology by Marshall et al. [Marshall, B.D., Neymark, L.A., Peterman, Z.E., 2003. Estimation of past seepage volumes from calcite distribution in the Topopah Spring Tuff, Yucca Mountain, Nevada. J. Contam. Hydrol. 62-63, 237-247] and Xu et al. [Xu, T., Sonnenthal, E., Bodvarsson, G., 2003. A reaction-transport model for calcite precipitation and evaluation of infiltration fluxes in unsaturated fractured rock. J. Contam. Hydrol. 64, 113-127] attempt to assess past volumes of seepage and infiltration fluxes through the vadose zone of Yucca Mountain, Nevada, on the basis of the modeling of the spatial distribution of secondary calcite. In this commentary, we argue that the employed methodology is not viable. In addition, the thermal boundary conditions used in simulations do not correspond to the temperatures of the mineral forming fluids established on the basis of the fluid inclusion studies.

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain., (Copyright 2004 Elsevier B.V.)

A numerical model was used to simulate the flow and transport processes at the Busted Butte Field Test Site for the purpose of quantifying the effects of hydrogeologic conditions beneath the potential Yucca Mountain repository horizon. In situ experiments were conducted on a 10 x 10 x 7 m block comprising a layered Topopah Springs/Calico Hills formation with two imbedded faults. Tracer solution was continuously injected in eight parallel boreholes arranged on two horizontal planes. Twelve collection boreholes were emplaced perpendicular to the injection holes and were both horizontal and inclined. Solution samples were collected regularly using a sampling assembly consisting of an inverted membrane and sorbing-paper sampling pads. Comparisons between measurements and predictions show that, except for the occasional drops of concentrations observed in the field, the current model is able to capture the general characteristics of the system with varying levels of agreement using laboratory-measured mean hydraulic properties. Simulation results and field observations revealed a capillary-driven flow in the system. Good quantitative agreement is generally observed for near-field boreholes, however, this agreement deteriorates and the simulated solute concentration is underestimated at boreholes farther away from the injection points. Increasing the spatial resolution of the simulation improves the model predictions only to a limited extent. Scaling issues may need to be considered to describe flow and transport events, as the travel distance becomes large.

Several numerical codes have been used to simulate radionuclide transport in fractured rock systems. The validation of such numerical codes can be accomplished by comparison of numerical simulations against appropriate analytical solutions. In this paper, we present analytical solutions for the reactive transport of N-member radionuclide chains (i.e., multiple species of radionuclides and their daughter species) through a discrete fracture in a porous rock matrix applying a system decomposition approach. We consider the transport of N-member radionuclide chains in a single-fracture-matrix system as a starting point to simulate more realistic and complex systems. The processes considered are advection along the fracture, lateral diffusion in the matrix, radioactive decay of multiple radionuclides, and adsorption in both the fracture and matrix. Different retardation factors can be specified for the fracture and matrix. However, all species are assumed to share the same retardation factors for the fracture and matrix, respectively. Although a daughter species may penetrate farther along the fracture than its parent species when a constant-concentration boundary condition is applied, our results indicate that all species retain the same transport speed in the fracture if a pulse of the first species is released into the fracture. This solution scheme provides a way to validate numerical computer codes of radionuclide transport in fractured rock, such as those being used to assess the performance of a potential nuclear-waste repository at Yucca Mountain.

This paper presents a numerical study on the response of the unsaturated zone (UZ) system of Yucca Mountain to heat generated from decaying radioactive wastes emplaced at the proposed repository. The modeling study is based on the current thermal-hydrological (TH) mountain-scale model, which uses a locally refined 2D north-south cross-section and dual-permeability numerical approach. The model provides a prediction of the mountain-scale TH response under the thermal-load scenario of 1.45 kW/m, while accounting for future climatic changes and the effects of drift ventilation. The TH simulation results show that ventilation of the repository drifts has a large impact on thermal-hydrologic regimes and moisture-flow conditions at the repository. In both cases, with and without ventilation, the TH model predicts dry or reduced liquid saturation near the drifts for over 1000 years, during which liquid flux through the drifts is reduced to either zero or less than the ambient flux. Without ventilation, the model predicts higher temperatures at the repository, but no major moisture redistribution in the UZ except in the areas very near the heated drifts.

This paper presents the application of parallel computing techniques to large-scale modeling of fluid flow in the unsaturated zone (UZ) at Yucca Mountain, Nevada. In this study, parallel computing techniques, as implemented into the TOUGH2 code, are applied in large-scale numerical simulations on a distributed-memory parallel computer. The modeling study has been conducted using an over-1-million-cell three-dimensional numerical model, which incorporates a wide variety of field data for the highly heterogeneous fractured formation at Yucca Mountain. The objective of this study is to analyze the impact of various surface infiltration scenarios (under current and possible future climates) on flow through the UZ system, using various hydrogeological conceptual models with refined grids. The results indicate that the 1-million-cell models produce better resolution results and reveal some flow patterns that cannot be obtained using coarse-grid modeling models.

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.

Because the continuum approach is relatively simple and straightforward to implement, it has been commonly used in modeling flow and transport in unsaturated fractured rock. However, the usefulness of this approach can be questioned in terms of its adequacy for representing fingering flow and transport in unsaturated fractured rock. The continuum approach thus needs to be evaluated carefully by comparing simulation results with field observations directly related to unsaturated flow and transport processes. This paper reports on such an evaluation, based on a combination of model calibration and prediction, using data from an infiltration test carried out in a densely fractured rock within the unsaturated zone of Yucca Mountain, Nevada. Comparisons between experimental and modeling results show that the continuum approach may be able to capture important features of flow and transport processes observed from the test. The modeling results also show that matrix diffusion may have a significant effect on the overall transport behavior in unsaturated fractured rocks, which can be used to estimate effective fracture-matrix interface areas based on tracer transport data. While more theoretical, numerical, and experimental studies are needed to provide a conclusive evaluation, this study suggests that the continuum approach is useful for modeling flow and transport in unsaturated, densely fractured rock.

A weakly sorbing cation, lithium, will be used as a reactive tracer in upcoming field tracer tests in the saturated alluvium south of Yucca Mountain, Nevada. One objective of the field tests is to determine how well field-scale reactive transport can be predicted using transport parameters derived from laboratory experiments. This paper describes several laboratory lithium batch sorption and column transport experiments that were conducted using ground water and alluvium obtained from the site of the planned field tests. In the batch experiments, isotherms were determined over 2.5 orders of magnitude of lithium concentrations, corresponding to the range expected in the field tests. In addition to measuring equilibrium lithium concentrations, concentrations of other cations, namely Na(+), K(+), and Ca(2+), were measured in the batch tests to determine Li(+)-exchangeable equilibria. This information was used in conjunction with alluvium cation exchange capacity measurements to parameterize a three-component cation-exchange model (EQUIL) that describes lithium sorption in the alluvium system. This model was then applied to interpret the transport behavior of lithium ion in saturated alluvium column tests conducted at three different lithium bromide injection concentrations. The concentrations were selected such that lithium ion either dominated, accounted for a little over half, or accounted for only a small fraction of the total cation equivalents in the injection solution. Although tracer breakthrough curves differed significantly under each of these conditions, with highly asymmetric responses occurring at the highest injection concentrations, the three-component cation-exchange model reproduced the observed transport behavior of lithium and the other cations in each case with a similar set of model parameters. In contrast, a linear K(d)-type sorption model could only match the lithium responses at the lowest injection concentration. The three-component model will be used to interpret the field tests, with the expectation that it will help refine estimates of effective flow porosity, particularly if the lithium response curves are asymmetric.

Performance assessment and design evaluation of the proposed repository at Yucca Mountain are facilitated by a thermohydrologic modeling tool that simultaneously accounts for processes occurring at a scale of a few tens of centimeters around individual waste packages and emplacement drifts, and accounts for processes at the multi-kilometer scale of the mountain. The most straightforward approach is to account for the 3-D drift- and mountain-scale dimensionality all within a single monolithic thermohydrologic model. This approach is too computationally expensive to be a viable simulation tool capable of addressing all waste-package locations in the repository. The Multiscale Thermohydrologic Model (MSTHM) is a computationally efficient alternative to addressing these modeling issues. In this paper, we describe the principal calculation stages to predict temperature, relative humidity, and liquid-saturation, as well as other thermohydrologic variables, in the drifts and in the host rock. Using a three-drift repository example (which is a scaled-down version of the proposed repository), we demonstrate the validity of the MSTHM approach against a nested monolithic thermohydrologic model.

This paper describes the development and use of a particle-tracking model to perform radionuclide-transport simulations in the unsaturated zone at Yucca Mountain, Nevada. The goal of the effort was to produce a computational model that can be coupled to the project's calibrated 3D site-scale flow model so that the results of that effort could be incorporated directly into the Total System Performance Assessment (TSPA) analyses. The transport model simulates multiple species (typically 20 or more) with complex time-varying and spatially varying releases from the potential repository. Water-table rise, climate-change scenarios, and decay chains are additional features of the model. A cell-based particle-tracking method was employed that includes a dual-permeability formulation, advection, longitudinal dispersion, matrix diffusion, and colloid-facilitated transport. This paper examines the transport behavior of several key radionuclides through the unsaturated zone using the calibrated 3D unsaturated flow fields. Computational results illustrate the relative importance of fracture flow, matrix diffusion, and lateral diversion on the distribution of travel times from the simulated repository to the water table for various climatic conditions. Results also indicate rapid transport through fractures for a portion of the released mass. Further refinement of the model will address several issues, including conservatism in the transport model, the assignment of parameters in the flow and transport models, and the underlying assumptions used to support the conceptual models of flow and transport in the unsaturated zone at Yucca Mountain.

Understanding thermally driven coupled hydrological, mechanical, and chemical processes in unsaturated fractured tuff is essential for evaluating the performance of the potential radioactive waste repository at Yucca Mountain, Nevada. The Drift Scale Test (DST), intended for acquiring such an understanding of these processes, has generated a huge volume of temperature and moisture redistribution data. Sophisticated thermal-hydrological (TH) conceptual models have yielded a good fit between simulation results and those measured data. However, some uncertainties in understanding the TH processes associated with the DST still exist. This paper evaluates these uncertainties and provides quantitative estimates of the range of these uncertainties. Of particular interest for the DST are the uncertainties resulting from the unmonitored loss of vapor through an open bulkhead of the test. There was concern that the outcome from the test might have been significantly altered by these losses. Using alternative conceptual models, we illustrate that predicted mean temperatures from the DST are within 1 degrees C of the measured mean temperatures through the first 2 years of heating. The simulated spatial and temporal evolution of drying and condensation fronts is found to be qualitatively consistent with measured saturation data. Energy and mass balance computation shows that no more than 13% of the input energy is lost because of vapor leaving the test domain through the bulkhead. The change in average saturation in fractures is also relatively small. For a hypothetical situation in which no vapor is allowed to exit through the bulkhead, the simulated average fracture saturation is not qualitatively different enough to be discerned by measured moisture redistribution data. This leads us to conclude that the DST, despite the uncertainties associated with open field testing, has provided an excellent understanding of the TH processes.