Flow and transport in fractured tuff at Yucca Mountain: numerical experiments on fast preferential flow mechanisms

Recent discovery of bomb-related ³⁶Cl at depth in fractured tuff in the unsaturated zone at the Yucca Mountain candidate high-level waste (HLW) repository site has called into question the usual modeling assumptions based on the equivalent continuum model (ECM). A dual continuum model (DCM) for simulating transient flow and transport at Yucca Mountain is developed. In order to ensure properly converged flow solutions, which are used in the transport simulation, a new flow solution convergence criteria is derived. An extensive series of simulation studies is presented which indicates that rapid movement of solute through the fractures will not occur unless there are intense episodic infiltration events. Movement of solute in the environs of the repository is enhanced if the properties of the tuff layer at the repository horizon are modified from current best-estimate values. Due to a large advective–dispersive coupling between the matrix and fractures, the matrix acts as a major buffer which inhibits rapid transport along the fractures. Consequently, fast movement of solutes through the fractures to the repository depth can only be explained if the matrix–fracture coupling term is significantly reduced from a value that would be calculated on the basis of data currently available.     

Characterization of flow and transport processes within the unsaturated zone of Yucca Mountain, Nevada, under current and future climates

This paper presents a large-scale modeling study characterizing fluid flow and tracer transport in the unsaturated zone of Yucca Mountain, Nevada, a potential repository site for storing high-level radioactive waste. The study has been conducted using a three-dimensional numerical model, which incorporates a wide variety of field data and takes into account the coupled processes of flow and transport in the highly heterogeneous, unsaturated fractured porous rock. The modeling approach is based on a dual-continuum formulation of coupled multiphase fluid and tracer transport through fractured porous rock. Various scenarios of current and future climate conditions and their effects on the unsaturated zone are evaluated to aid in the assessment of the proposed repository's system performance using different conceptual models. These models are calibrated against field-measured data. Model-predicted flow and transport processes under current and future climates are discussed.     

A serendipitous, long-term infiltration experiment: Water and tritium circulation beneath the CAMBRIC trench at the Nevada Test Site

Underground nuclear weapons testing at the Nevada Test Site introduced numerous radionuclides that may be used subsequently to characterize subsurface hydrologic transport processes in arid climates. In 1965, a unique, 16-year pumping experiment designed to examine radionuclide migration away from the CAMBRIC nuclear test, conducted in the saturated zone beneath Frenchman Flat, Nevada, USA, gave rise to an unintended second experiment involving radionuclide infiltration through the vadose zone, as induced by seepage of pumping effluents beneath an unlined discharge trench. The combined experiments have been reanalyzed using a detailed, three-dimensional numerical model of transient, variably saturated flow and mass transport in a heterogeneous subsurface, tailored specifically for large-scale and efficient calculations. Simulations have been used to estimate tritium travel and residence times in various parts of the system for comparison with observations in wells. Model predictions of mass transport were able to clearly demonstrate radionuclide recycling behavior between the trench and pumping well previously suggested by isotopic age dating information; match travel time estimates for radionuclides moving between the trench, the water table, and monitoring and pumping wells; and provide more realistic ways in which to interpret the pumping well elution curves. Collectively, the results illustrate the utility of integrating detailed numerical modeling with diverse observational data in developing more accurate interpretations of contaminant migration processes.     

Assessing tungsten transport in the vadose zone: from dissolution studies to soil columns

This study investigates the dissolution, sorption, leachability, and plant uptake of tungsten and alloying metals from canister round munitions in the presence of model, well characterized soils. The source of tungsten was canister round munitions, composed mainly of tungsten (95%) with iron and nickel making up the remaining fraction. Three soils were chosen for the lysimeter studies while four model soils were selected for the adsorption studies. Lysimeter soils were representatives of the typical range of soils across the continental USA; muck-peat, clay-loamy and sandy-quartzose soil. Adsorption equilibrium data on the four model soils were modeled with Langmuir and linear isotherms and the model parameters were obtained. The adsorption affinity of soils for tungsten follows the order: Pahokee peat > kaolinite > montmorillonite > illite. A canister round munition dissolution study was also performed. After 24 d, the measured dissolved concentrations were: 61.97, 3.56, 15.83 mg L⁻¹ for tungsten, iron and nickel, respectively. Lysimeter transport studies show muck peat and sandy quartzose soils having higher tungsten concentration, up to 150 mg kg⁻¹ in the upper layers of the lysimeters and a sharp decline with depth suggesting strong retardation processes along the soil profile. The concentrations of tungsten, iron and nickel in soil lysimeter effluents were very low in terms of posing any environmental concern; although no regulatory limits have been established for tungsten in natural waters. The substantial uptake of tungsten and nickel by ryegrass after 120 d of exposure to soils containing canister round munition suggests the possibility of tungsten and nickel entering the food chain.     

Tritium and Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and ³⁶Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this “fast flow” in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891–989] and suggest that fast flow in fractures with minimal fracture–matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and ³⁶Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.     

Uranium transport around the reactor zone at Bangombé and Okélobondo (Oklo): examples of hydrogeological and geochemical model integration and data evaluation

The sites at Bangombé and Okélobondo (Oklo) in Gabon provide a unique opportunity to study the behaviour of products from natural nuclear reactions in the vicinity of reactor zones which were active around two billion years ago. The Commission of the European Communities initiated the Oklo Natural Analogue Programme. One of the principal aims was to study indications of present time migration of elements from the reactor zones under ambient conditions. The hydrogeological and hydrochemical data from the Oklo sites were modelled in order to better understand the geochemical behaviour of radionuclides in the natural system, by using independent models and by comparing the modelling outcome. Two modelling approaches were used: M3 code (hydrochemical mixing and mass balance model), developed by the Swedish Nuclear Fuel and Waste Management Company (SKB) and HYTEC (reactive transport model) developed by Ecole des Mines de Paris.Two different reactor zones were studied: Bangombé, a shallow site, the reactor being at 11 m depth, and OK84 at Okélobondo, situated at about 450 m depth, more comparable with a real repository location. This allowed the validation of modelling tools in two different sedimentary environments: one shallow, with a more homogeneous layering situated in an area of meteoric alteration, and the other offering the opportunity to study radionuclide migration from the reaction zone over a distance of 450 m through very heterogeneous sedimentary layers.The modeling results indicate that the chemical reactions retarding radionuclide transport are very different at the two sites. At Bangombé, the decomposition of organic material consumes oxygen and at Okélobondo the oxygen is consumed by inorganic reactions resulting, in both cases, in uranium retardation. Both modelling approaches (statistic with M3 code and deterministic with HYTEC code) could describe this situation.The goal of this exercise is to test codes which can help to describe and understand the processes taking place at the sites, validate the models with in situ data, and thus build confidence in the tools used for future site characterization. Ultimately, this allows identifying and selecting processes and parameters that can be used as input into repository performance assessment calculations and modelling exercises.     

Chemical fate,latitudinal distribution and long-range transport of cyclic volatile methylsiloxanes in the global environment: A modeling assessment

Cyclic volatile methylsiloxanes (cVMS) such as octamethycyclotetrasiloxane (D4), decamethycyclopentasiloxane (D5), and dodecamethylcyclohexasiloxane (D6) are widely used as intermediates in the synthesis of high-molecular weight silicone polymers or as ingredients in the formulation of personal care products. The global environmental fate, latitudinal distribution, and long range transport of those cVMS were analyzed by two multimedia chemical fate models using the best available physicochemical properties as inputs and known persistent organic pollutants (POPs) and highly persistent volatile organic chemicals (“fliers”) as reference. The global transport and accumulation characteristics of cVMS differ from those of typical POPs in three significant ways. First, a large fraction of the released cVMS tends to become airborne and is removed from the global environment by degradation in air, whereas known POPs have a tendency to be distributed and persistent in all media. Secondly, although cVMS can travel a substantial distance in the atmosphere, they have little potential for deposition to surface media in remote regions. This contrasts with a deposition potential of known POPs that exceeds that of cVMS by 4–5 orders of magnitude. Thirdly, cVMS have short global residence times with the majority of the global mass removed within 3 months of the end of release. Global residence times of POPs on the other hand are in years. The persistent fliers resemble the cVMS with respect to the first two attributes, but their global residence times are more like those of the POPs.