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.     

Comparison of the solution behaviour of a pyrite-calcite mixture in batch and unsaturated sand column

A successful application of reaction transport algorithms to calculate the chemical evolution of natural systems requires accurate methods to compute the rates of mineral/fluid surface reactions. Regarding the transport of radio-nuclides in mining dumps the dissolution of minerals is of special importance. Using a kinetic rate law of the mineral dissolution verified for unsaturated conditions will allow a realistic modelling of the mineral weathering in the environment. Dissolution rates of minerals in an aqueous solution are determined by several characteristics. These are surface reaction rates, morphology of the mineral's surface and, in case it is the unsaturated zone, the degree of the water saturation. For this process, the quantity of the particle surfaces which are in contact with percolating water is most decisive. In order to study the differences of mineral dissolution under saturated and unsaturated conditions batch and column experiments were carried out with a pyrite-calcite mixture. The experimental results were verified by calculations. Comparing the dissolution in batch with those in the column experiment, which was performed with a water flow velocity of 0.64 cm/day and was analyzed in the region of a water saturation of 0.11, one can conclude that only a small portion of about 5% of the grain surface is chemically reactive in this unsaturated flow.     

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.