Rock matrix diffusivity determinations by in-situ electrical conductivity measurements

A fast method to determine rock matrix diffusion properties directly in the bedrock would be valuable in the investigation of a possible site for disposal of nuclear waste. An "effective diffusivity borehole log" would provide important information on the variability of this entity over the area studied. As opposed to traditional matrix diffusion laboratory experiments, electrical conductivity measurements are fast, inexpensive and also easy to carry out in-situ. In this study, electrical resistivity data from borehole logging, as well as from measurements on the actual core, is evaluated with the purpose of extracting matrix diffusivity data. The influence of migration of ions in the electrical double layer, which can be of great importance in low ionic strength pore water, is also considered in evaluating the in-situ data to accurately determine the effective pore diffusivity. The in-situ data compare fairly well to those measured in the rock core.     

An axisymmetric diffusion experiment for the determination of diffusion and sorption coefficients of rock samples

Diffusion anisotropy is a critical property in predicting migration of substances in sedimentary formations with very low permeability. The diffusion anisotropy of sedimentary rocks has been evaluated mainly from laboratory diffusion experiments, in which the directional diffusivities are separately estimated by through-diffusion experiments using different rock samples, or concurrently by in-diffusion experiments in which only the tracer profile in a rock block is measured. To estimate the diffusion anisotropy from a single rock sample, this study proposes an axisymmetric diffusion test, in which tracer diffuses between a cylindrical rock sample and a surrounding solution reservoir. The tracer diffusion between the sample and reservoir can be monitored from the reservoir tracer concentrations, and the tracer profile could also be obtained after dismantling the sample. Semi-analytical solutions are derived for tracer concentrations in both the reservoir and sample, accounting for an anisotropic diffusion tensor of rank two as well as the dilution effects from sampling and replacement of reservoir solution. The transient and steady-state analyses were examined experimentally and numerically for different experimental configurations, but without the need for tracer profiling. These experimental configurations are tested for in- and out-diffusion experiments using Koetoi and Wakkanai mudstones and Shirahama sandstone, and are scrutinized by a numerical approach to identify favorable conditions for parameter estimation. The analysis reveals the difficulty in estimating diffusion anisotropy; test configurations are proposed for enhanced identifiability of diffusion anisotropy. Moreover, it is demonstrated that the axisymmetric diffusion test is efficient in obtaining the sorption parameter from both steady-state and transient data, and in determining the effective diffusion coefficient if isotropic diffusion is assumed. Moreover, measuring reservoir concentrations in an axisymmetric diffusion experiment coupled with tracer profiling may be a promising approach to estimate of diffusion anisotropy of sedimentary rocks.     

Prediction of some in situ tracer tests with sorbing tracers using independent data

Some recent converging tracer tests with sorbing tracers at the Asp? Hard Rock Laboratory in Sweden, the TRUE tests, have been predicted using only laboratory data and hydraulic data from borehole measurements. No model parameters were adjusted to obtain a better fit with the experiments. The independent data were fracture frequency and transmissivity data obtained in the field and laboratory data on sorption and matrix diffusion. Transmissivity measurements in five boreholes in the rock volume containing the region surrounding the injection and collection points show that there is a high frequency of water conducting fractures. Of 162 packed off sections with 0.5 m packer distances, 112 were found to have a transmissivity above the detection limit. The specific flow-wetted surface (FWS) of the rock mass could be estimated from these data. The transmissivities were found to be reasonably well described by a lognormal distribution. Laboratory data on diffusion and sorption properties together with the hydraulic data were used to "predict" the residence time distribution (RTD) of the sorbing tracers. The results were compared with the experimental breakthrough curves. In these experiments, the water residence time is very small compared to the residence time of the sorbing tracers due to their diffusion and sorption within the rock matrix. We thus could neglect the influence of the water residence time in our predictions. Therefore, no information on water residence times or on "dispersion" was needed. The dispersion of the sorbing tracers is caused by the different sorbing tracer residence times in different pathways. The sorbing tracer residence time is determined by the ratio of flowrate to the flow-wetted surface in the different pathways and not by the water residence time. Assuming a three-dimensional flow pattern and using the observed fracture frequency and flowrate distribution, breakthrough curves for three strongly sorbing tracers were predicted. Only the laboratory data, the transmissivity measurements and the pumping flowrate were used in the predictions. No information on the water residence time as obtained by the nonsorbing tracers was used. The predictions were surprisingly accurate.     

Characterisation of HTO diffusion properties by an in situ tracer experiment in Opalinus clay at Mont Terri

A long-term single borehole diffusion experiment using tritiated water as tracer was carried out in Opalinus clay, an argillaceous rock formation that is accessible at the Mont Terri Underground Research Laboratory, situated in the Swiss Jura. The tracer was diluted in reconstituted formation water and introduced into a packed-off section of a borehole located in saturated rock. Pressure in this interval was maintained equal to the pore pressure of the surrounding rock in order to prevent any hydraulic gradient around the borehole and to avoid advective transport processes. The evolution of the tracer concentration in the injection system was monitored over time. After 1 year of diffusion, the claystone surrounding the interval was retrieved by overcoring the whole borehole and packer system, and by an adjacent oblique borehole. Compressed air was used as drilling fluid to reduce rock disturbances. The recovered overcore was sampled along profiles perpendicular to the borehole wall with a view to determining the tracer-concentration profiles in the rock. To avoid further evaporation of tritiated water, subsamples were immediately transferred into polyethylene bottles and disaggregated by adding a known amount of tracer-free water. Fifteen profiles were determined and showed a decreasing tracer concentration with distance into the rock. The pore-water contents were constant along those profiles, confirming that only very little water was lost during overcoring operations. The evolution of tritium-tracer concentration in the injection system over time and in situ profiles were interpreted with a 3-D numerical simulation of the experiment. That allowed for the identification of the transport parameters (orthotropic diffusion tensor and porosity) by minimising the relative quadratic error between the experimental and simulated data. The fitting is good and the results are consistent with data obtained on drill-core samples. The result of tritiated water is discussed regarding (1) the potential effect of mechanical and/or chemical disturbances around the injection borehole and (2) the specific behaviour of tritiated water.     

Tracer diffusion coefficients in sedimentary rocks: correlation to porosity and hydraulic conductivity

Matrix diffusion is an important transport process in geologic materials of low hydraulic conductivity. For predicting the fate and transport of contaminants, a detailed understanding of the diffusion processes in natural porous media is essential. In this study, diffusive tracer transport (iodide) was investigated in a variety of geologically different limestone and sandstone rocks. Porosity, structural and mineralogical composition, hydraulic conductivity, and other rock properties were determined. The effective diffusion coefficients were measured using the time-lag method. The results of the diffusion experiments indicate that there is a close relationship between total porosity and the effective diffusion coefficient of a rock (analogous to Archie's Law). Consequently, the tortousity factor can be expressed as a function of total porosity. The relationship fits best for thicker samples (> 1.0 cm) with high porosities (> 20%), because of the reduced influence of heterogeneity in larger samples. In general, these correlations appear to be a simple way to determine tortuosity and the effective diffusion coefficient from easy to determine rock porosity values.     

Interpretation of injection-withdrawal tracer experiments conducted between two wells in a large single fracture

Tracer experiments conducted using a flow field established by injecting water into one borehole and withdrawing water from another are often used to establish connections and investigate dispersion in fractured rock. As a result of uncertainty in the uniqueness of existing models used for interpretation, this method has not been widely used to investigate more general transport processes including matrix diffusion or advective solute exchange between mobile and immobile zones of fluid. To explore the utility of the injection-withdrawal method as a general investigative tool and with the intent to resolve the transport processes in a discrete fracture, two tracer experiments were conducted using the injection-withdrawal configuration. The experiments were conducted in a fracture which has a large aperture (>500 microm) and horizontally pervades a dolostone formation. One experiment was conducted in the direction of the hydraulic gradient and the other in the direction opposite to the natural gradient. Two tracers having significantly different values of the free-water diffusion coefficient were used. To interpret the experiments, a hybrid numerical-analytical model was developed which accounts for the arcuate shape of the flow field, advection-dispersion in the fracture, diffusion into the matrix adjacent to the fracture, and the presence of natural flow in the fracture. The model was verified by comparison to a fully analytical solution and to a well-known finite-element model. Interpretation of the tracer experiments showed that when only one tracer, advection-dispersion, and matrix diffusion are considered, non-unique results were obtained. However, by using multiple tracers and by accounting for the presence of natural flow in the fracture, unique interpretations were obtained in which a single value of matrix porosity was estimated from the results of both experiments. The estimate of porosity agrees well with independent measurements of porosity obtained from core samples. This suggests that: (i) the injection-withdrawal method is a viable tool for the investigation of general transport processes provided all relevant experimental conditions are considered and multiple conservative tracers are used; and (ii) for the conditions of the experiments conducted in this study, the dominant mechanism for exchange of solute between the fracture and surrounding medium is matrix diffusion.     

Diffusion of tritiated water and 22Na+ through non-degraded hardened cement pastes

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and 22Na(+), considered to be conservative tracers, have been carried out in triplicates in a glove box under a controlled nitrogen atmosphere. Each experiment consisted of a through-diffusion test followed by an out-diffusion test. The experimental data were inversely modelled applying an automated Marquardt-Levenberg procedure. The analysis of the through-diffusion data allowed the extraction of values for the effective diffusion coefficients, D(e), and the rock capacity factor, alpha. Good agreement between measured and calculated tracer breakthrough curves was achieved using both a simple diffusion model without sorption and a diffusion/linear sorption model. The best-fit K(d)-values were found to be consistent with R(d)-values measured in previous batch-sorption experiments. The best-fit values from the through-diffusion tests were then used to predict the results of subsequent out-diffusion experiments. Good agreement between experimental data and predictions was achieved only for the case of linear sorption. Isotopic exchange can only partially account for both the amount of tracer taken up in the batch-sorption tests and the measured retardation in the diffusion experiments and, hence, additional mechanisms have to be invoked to explain the data.     

Modelling the migration of contaminants through variably saturated dual-porosity, dual-permeability chalk

In the Hesbaye region in Belgium, tracer tests performed in variably saturated fissured chalk rocks presented very contrasting results in terms of transit times, according to artificially controlled water recharge conditions prevailing during the experiments. Under intense recharge conditions, tracers migrated across the partially or fully saturated fissure network, at high velocity in accordance with the high hydraulic conductivity and low effective porosity (fracture porosity). At the same time, a portion of the tracer was temporarily retarded in the almost immobile water located in the matrix. Under natural infiltration conditions, the fissure network remained inactive. Tracers migrated downward through the matrix, at low velocity in relation with the low hydraulic conductivity and the large porosity of the matrix. Based on these observations, Brouyère et al. (2004a) [Brouyère, S., Dassargues, A., Hallet, V., 2004a. Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation, J. Contam. Hydrol., 72 (1-4), 135-164, doi: 10.1016/j.conhyd.2003.10.009] proposed a conceptual model in order to explain the migration of solutes in variably saturated, dual-porosity, dual-permeability chalk. Here, mathematical and numerical modelling of tracer and contaminant migration in variably saturated fissured chalk is presented, considering the aforementioned conceptual model. A new mathematical formulation is proposed to represent the unsaturated properties of the fissured chalk in a more dynamic and appropriate way. At the same time, the rock water content is partitioned between mobile and immobile water phases, as a function of the water saturation of the chalk rock. The groundwater flow and contaminant transport in the variably saturated chalk is solved using the control volume finite element method. Modelling the field tracer experiments performed in the variably saturated chalk shows the adequacy and usefulness of the new conceptual, mathematical and numerical model.     

Modeling tracer diffusion in fractured and unfractured, unsaturated, porous media

A proposed tracer diffusion test for the Exploratory Shaft Facility at Yucca Mountain, NV, is modeled. For the proposed test, a solution containing conservative tracers will be introduced into a borehole in the geologic medium of interest. The tracers will diffuse and advect from the saturated source region into the unsaturated matrix in the surrounding tuff. After some time, the borehole is to be overcored, and tracer concentrations in the fluid will be measured in the core as a function of distance from emplacement. The data will be used to evaluate diffusive behavior and to derive effective diffusion coefficients for the tracers in the specific tuff. Numerical simulations are used to study the effects of effective diffusion coefficient, porosity, saturation, and fracturing on tracer transport. Results are reported for numerical simulations of tests in the Topopah Spring Member and the Tuff of Calico Hills, which have significantly different porosities and saturations. The simulations make the following predictions: The spread of tracer during the test will be sensitive to the effective diffusion coefficient of the tracer. Tracer will diffuse farther in the Topopah Spring Member than in the Tuff of Calico Hills because of the former's lower porosity and saturation. Tracer transport by advection into the Topopah Spring Member will be greater than that into the Tuff of Calico Hills because of capillary effects. While advection will be a significant mechanism for tracer penetration into the Topopah Spring tuff, it will be less significant for tracer penetration into the Calico Hills tuff. The proximity of a single vertical fracture to the source region determines its effects on tracer transport, especially if the fracture diverts fluid flowing from the source region into the matrix.     

Modelling of nonreactive tracer dipole tests in a shear zone at the Grimsel test site

The investigation of the migration of a high pH plume in a fractured shear zone is foreseen by a long-term experiment at the Grimsel rock laboratory. In order to characterise the initial conditions for the long-term experiment and to evaluate an optimal hydraulic in situ set-up, several dipole experiments with nonreacting tracers have been performed. The dipole experiments differ in geometry, pumping rates and orientation to the background water flow. Several single and double-porosity models have been applied to fit the results of these dipole tracer tests in order to extract values for some transport parameters and discriminate for certain transport processes. A two-dimensional porous medium approach was successfully used to fit tracer breakthrough curves measured for a dipole experiment. A model based on a one-dimensional dual porous medium approach was also successful, although the applied hydraulic dipole, with similar injection and extraction rates, suggests the existence of an extended two-dimensional flow field. For the two-dimensional porous medium approach, tracer breakthrough could only be fitted with a complex flow field geometry within the heterogeneous fractured shear zone. The heterogeneity was generated by heterogeneous porosity and hydraulic permeability distributions. Predictions for further dipole geometries and a sorbing tracer have been calculated by means of both models using the flow and transport parameters deduced from fits for a single dipole experiment. This allows for comparison with the measured breakthrough of sorbing tracers. The foreseen experiment with sorbing (radionuclide) tracers will help decide on the appropriate approach that should be used to describe such dipole experiments in this shear zone. Additionally, the migration and spreading of a solution with high pH has been calculated taking into account mineral dissolution and precipitation in a two-dimensional porous medium approach in order to estimate the amount and character of the mineral reactions induced by the interaction between the high pH solution and the rock.     

Contribution of the diffusive exchange method to the characterization of pore-water in consolidated argillaceous rocks

Consolidated argillaceous rocks, potential host rocks for nuclear waste disposal, display particular properties (low water content, very small pore size) that require specific methods, especially for the characterization of their pore-water. The diffusive exchange method was applied to the claystone of Tournemire by taking into account the effect of swelling or micro-cracks induced by drilling. The good agreement of diffusion parameters with those obtained from the through-diffusion method showed that the extent of this disturbance was low. This method led to stable isotope content values in pore-water in good agreement with those measured in the fracture fluids located in the vicinity, validating the method. On the contrary, the vacuum distillation technique gave a systematic depletion in heavy isotopes of water, probably induced by an incomplete extraction of pore-water. For the chloride tracer, the determination of the accessible porosity on each sample allowed the estimate of chloride concentration of pore-water to be improved, by taking into account all the petrophysical heterogeneities.     

Use of tracer tests to evaluate the impact of enhanced-solubilization flushing on in-situ biodegradation

Tracer tests were conducted to evaluate the effect of a complexing sugar flush (CSF) on in-situ biodegradation potential at a site contaminated by jet fuel, solvents, and other organic compounds. Technical-grade hydroxypropyl-beta-cyclodextrin was used during the CSF study, which was conducted in a hydraulically isolated cell emplaced in a surficial aquifer. In-situ biodegradation potential was assessed with the use of tracer tests, which were conducted prior to and immediately following the CSF study. Ethanol, hexanol, and benzoate were used as the biodegradable tracers, while bromide was used as a nonreactive tracer. The results indicate that the biodegradation of benzoate was similar for both tracer tests. Conversely, the biodegradation of ethanol (23% increase) and hexanol (41% increase) was greater for the post-CSF tracer test. In addition, analysis of core samples collected from within the test cell indicates that the population density of aerobic jet-fuel degraders increased in the vicinity of the injection wells during the CSF. These results indicate that the cyclodextrin flush did not deleteriously affect the indigenous microbial community. This study illustrates that tracer tests can be used to evaluate the impact of remediation activities on in-situ biodegradation potential.     

Diffusion of tritiated water and Na through non-degraded hardened cement pastes

Diffusion experiments through hardened cement pastes (HCP) using tritiated water (HTO) and ²²Na⁺, considered to be conservative tracers, have been carried out in triplicates in a glove box under a controlled nitrogen atmosphere. Each experiment consisted of a through-diffusion test followed by an out-diffusion test.The experimental data were inversely modelled applying an automated Marquardt–Levenberg procedure. The analysis of the through-diffusion data allowed the extraction of values for the effective diffusion coefficients, D_e, and the rock capacity factor, α. Good agreement between measured and calculated tracer breakthrough curves was achieved using both a simple diffusion model without sorption and a diffusion/linear sorption model. The best-fit K_d-values were found to be consistent with R_d-values measured in previous batch-sorption experiments.The best-fit values from the through-diffusion tests were then used to predict the results of subsequent out-diffusion experiments. Good agreement between experimental data and predictions was achieved only for the case of linear sorption.Isotopic exchange can only partially account for both the amount of tracer taken up in the batch-sorption tests and the measured retardation in the diffusion experiments and, hence, additional mechanisms have to be invoked to explain the data.     

Analysis of a gas-phase partitioning tracer test conducted in an unsaturated fractured-clay formation

The gas-phase partitioning tracer method was used to estimate non-aqueous phase liquid (NAPL), water, and air saturations in the vadose zone at a chlorinated-solvent contaminated field site in Tucson, AZ. The tracer test was conducted in a fractured-clay system that is the confining layer for the underlying regional aquifer. Three suites of three tracers were injected into wells located 14, 24, and 24 m from a single, central extraction well. The tracers comprised noble gases (traditionally thought to be nonsorbing), alkanes (primarily water partitioning), perfluorides (primarily NAPL partitioning), and halons (both NAPL and water partitioning). Observations of vacuum response were consistent with flow in a fractured system. The halon tracers exhibited the greatest amount of retardation, and helium and the perfluoride tracers the least. The alkane tracers were unexpectedly more retarded than the perfluoride tracers, indicating low NAPL saturations and high water saturations. An NAPL saturation of 0.01, water saturation of 0.215, and gas saturation of 0.775 was estimated based on analysis of the suite of tracers comprising helium, perfluoromethylcyclohexane and dibromodifluoromethane, which was considered to be the most robust set. The estimated saturations compare reasonably well to independently determined values.     

Tracer test for the measurement of gas diffusion and non-aqueous phase liquid (NAPL) saturation in soil

During soil bioremediation, the diffusion of oxygen into the soil is an important prerequisite for aerobic biodegradation, and the decrease of petroleum products is the ultimate goal. Both processes need to be monitored. The aim of this work was to develop a gas tracer test that yields information on both, gas diffusion and residual saturation with non-aqueous phase liquids (NAPLs) in unsaturated soil heaps. One conservative tracer (methane) and 4 partitioning gas tracers (diethylether, methyl tert-butyl ether, chloroform and n-heptane) were injected as vapors into laboratory columns filled with unsaturated sand with increasing NAPL saturation. Breakthrough curves of gaseous compounds were measured at two points and compared to analytical solutions of an analytical diffusive-reactive transport equation. By fitting of methane data, robust results for effective diffusivity (tortuosity) were obtained. NAPL saturation was most accurately measured by the moderately water soluble tracers (ethers and chloroform). The hydrophobic tracer n-heptane did not partition into water-immersed NAPL. An easy and accurate way to assess air-NAPL partitioning constants from gas chromatography retention times is furthermore reported. It is concluded that gas tracer tests have the potential for measuring two important properties in soil bioremediation systems easily and quickly.     

Partitioning Gas Tracer Tests for Measurement of Water in Municipal Solid Waste

Abstract

A key component in the operation of almost all bioreactor landfills is the addition of water to maintain optimal moisture conditions. To determine how much water is needed and where to add it, in situ methods are required to measure water within solid waste. Existing technologies often result in measurements of unknown accuracy, because of the variability of solid waste materials and time-dependent changes in packing density, both of which influence most measurement methods. To overcome these problems, a new technology recently developed by hydrologists for measuring water in the vadose zone—the partitioning gas tracer test—was tested. In this technology, the transport behavior of two gas tracers within solid waste is used to measure the fraction of the void space filled with water. One tracer is conservative and does not react with solids or liquids, while a second tracer partitions into the water and is separated from the conservative tracer during transport. This technology was tested in four different solid waste packings and was capable of determining the volumetric water content to within 48% of actual values, with most measurement errors less than 15%. This technology and the factors that affect its applicability to landfills are discussed in this paper.     

Monitoring nutrients fate after digestate spreading in a short rotation buffer area

One of the main sources of reactive nitrogen pollution is animal manure. The disposal of digestate (material remaining after the anaerobic digestion of a biodegradable feedstock) in agricultural soils could solve both the problems of soil fertilization and waste removal, but the fate of digestate in the environment must be assessed carefully before its massive utilization. To investigate whether digestate could be safely employed as a soil fertilizer, an agricultural field located in Monastier di Treviso (Northern Italy) and characterized by the presence of low hydraulic conductivity clay soils, was selected to be amended with bovine digestate. The experimental site was intensively monitored by a three-dimensional array of probes recording soil water content, temperature, and electrical conductivity, to solve the water and bulk mass fluxes in the unsaturated zone. High-resolution soil coring allowed the characterization of soil water composition over two hydrological years. Chloride, found in high concentrations in the digestate, was used as environmental tracer to track the fate of the percolating water. The study concluded that digestate could be confidently employed in short rotation buffer areas at an average rate of 195 ± 26 kg-N/ha/year in low hydraulic conductivity soils not affected by diffuse fracturing during dry periods.     

Evaluation of groundwater flow patterns around a dual-screened groundwater circulation well

Dual-screened groundwater circulation wells (GCWs) can be used to remove contaminant mass and to mix reagents in situ. GCWs are so named because they force water in a circular pattern between injection and extraction screens. The radial extent, flux and direction of the effective flow of this circulation cell are difficult to measure or predict. The objective of this study is to develop a robust protocol for assessing GCW performance. To accomplish this, groundwater flow patterns surrounding a GCW are assessed using a suite of tools and data, including: hydraulic head, in situ flow velocity, measured hydraulic conductivity data from core samples, chemical tracer tests, contaminant distribution data, and numerical flow and transport models. The hydraulic head data show patterns that are consistent with pumping on a dual-screened well, however, many of the observed changes are smaller than expected. In situ thermal perturbation flow sensors successfully measured horizontal flow, but vertical flow could not be determined with sufficient accuracy to be useful in mapping flow patterns. Two types of chemical tracer tests were utilized at the site and showed that much of the flow occurs within a few meters of the GCW. Flow patterns were also assessed based on changes in contaminant (trichloroethylene, TCE) concentrations over time. The TCE data clearly showed treated water moving away from the GCW at shallow and intermediate depths, but the circulation of that water back to the well, except very close to the well, was less clear. Detailed vertical and horizontal hydraulic conductivities were measured on 0.3 m-long sections from a continuous core from the GCW installation borehole. The measured vertical and horizontal hydraulic conductivity data were used to construct numerical flow and transport models, the results of which were compared to the head, velocity and concentration data. Taken together, the field data and modeling present a fairly consistent picture of flow and transport around the GCW. However, the time and expense associated with conducting all of those tests would be prohibitive for most sites. As a consequence, a sequential protocol for GCW characterization is presented here in which the number of tools used can be adjusted to meet the needs of individual sites. While not perfect, we believe that this approach represents the most efficient means for evaluating GCW performance.     

Field investigation into unsaturated flow and transport in a fault: model analyses

Results of a fault test performed in the unsaturated zone of Yucca Mountain, Nevada, were analyzed using a three-dimensional numerical model. The fault was explicitly represented as a discrete feature and the surrounding rock was treated as a dual-continuum (fracture-matrix) system. Model calibration against seepage and water-travel-velocity data suggests that lithophysal cavities connected to fractures can considerably enhance the effective fracture porosity and therefore retard water flow in fractures. Comparisons between simulation results and tracer concentration data also indicate that matrix diffusion is an important mechanism for solute transport in unsaturated fractured rock. We found that an increased fault-matrix and fracture-matrix interface areas were needed to match the observed tracer data, which is consistent with previous studies. The study results suggest that the current site-scale model for the unsaturated zone of Yucca Mountain may underestimate radionuclide transport time within the unsaturated zone, because an increased fracture-matrix interface area and the increased effective fracture porosity arising from lithophysal cavities are not considered in the current site-scale model.