Lawrence Berkeley National Laboratory

Characterizing percolation patterns in unsaturated zones has posed a greater challenge to numerical modeling investigations than comparable saturated zone studies, because of the heterogeneous nature of unsaturated media as well as the great number of variables impacting unsaturated zone flow. This paper presents an integrated modeling methodology for quantitatively characterizing percolation patterns in the unsaturated zone of Yucca Mountain, Nevada, a proposed underground repository site for storing high-level radioactive waste. It takes into account the multiple coupled processes of air, water, heat flow and chemical isotopic transport in Yucca Mountain s highly heterogeneous, unsaturated fractured tuffs. The modeling approach integrates a wide variety of moisture, pneumatic, thermal, and isotopic geochemical field data into a comprehensive three-dimensional numerical model for modeling analyses. Modeling results are examined against different types of field-measured data and then used to evaluate different hydrogeological conceptual models and their results of flow patterns in the unsaturated zone. In particular, this integration model provides a much clearer understanding of percolation patterns and flow behavior through the unsaturated zone, both crucial issues in assessing repository performance. The integrated approach for quantifying Yucca Mountain s flow system is also demonstrated to provide a comprehensive modeling tool for characterizing flow and transport processes in complex subsurface systems.