Surfactant-induced changes in gravity fingering of water through a light oil

Gravity-driven preferential flow (fingering) can greatly affect how one fluid displaces another in the subsurface. We have studied the internal properties of these preferential flow paths for water, with and without surfactants, infiltrating into oil saturated porous media using synchrotron X-rays, and miniature tensiometers to characterize fluid content and pressure relationships. We also used a light transmission technique to visualize overall flow pattern. Capillary pressure and water content decrease behind the front, similar to fingers in air-dry sand, with quantitative differences for five different surfactants with surface tensions ranging from 4–21 g/s². Using unstable flow theory, the finger widths, capillary pressure drops within the fingers, finger tip lengths, and finger splitting dynamics were scaled successfully with interfacial tension, fluid density, and the contact angle using the fingers in air–water systems as the reference.     

A multi-flowpath model for the interpretation of immiscible displacement experiments in heterogeneous soil columns

This work focuses on the phenomenon of the immiscible two-phase flow of water and oil in saturated heterogeneous soil columns. The goal is to develop a fast and reliable method for quantifying soil heterogeneities for incorporation into the relevant capillary pressure and relative permeability functions. Such data are commonly used as input data in simulators of contaminant transport in the subsurface. Rate-controlled drainage experiments are performed on undisturbed soil columns and the transient response of the axial distribution of water saturation is determined from electrical measurements. The transient responses of the axial distribution of water saturation and total pressure drop are fitted with the multi-flowpath model (MFPM) where the pore space is regarded as a system of parallel paths of different permeability. The MFPM enables us to quantify soil heterogeneity at two scales: the micro-scale parameters describe on average the effects of pore network heterogeneities on the two-phase flow pattern; the macro-scale parameters indicate the variability of permeability at the scale of interconnected pore networks. The capillary pressure curve is consistent with that measured with mercury intrusion porosimetry over the low pressure range. The oil relative permeability increases sharply at a very low oil saturation (< 10^− 3) and tends to a high end value. The water relative permeability decreases abruptly at a low oil saturation (~ 0.1), whereas the irreducible wetting phase saturation is quite high. The foregoing characteristics of the two-phase flow properties are associated with critical (preferential) flowpaths that comprise a very small percentage of the total pore volume, control the overall hydraulic conductivity, and are consistent with the very broad range of pore-length scales usually probed in soil porous matrix.     

Aquifer washing by micellar solutions: 1: Optimization of alcohol–surfactant–solvent solutions

Phase diagrams were used for the formulation of alcohol–surfactant–solvent and to identify the DNAPL (Dense Non Aqueous Phase Liquid) extraction zones. Four potential extraction zones of Mercier DNAPL, a mixture of heavy aliphatics, aromatics and chlorinated hydrocarbons, were identified but only one microemulsion zone showed satisfactory DNAPL recovery in sand columns. More than 90 sand column experiments were performed and demonstrate that: (1) neither surfactant in water, alcohol–surfactant solutions, nor pure solvent can effectively recover Mercier DNAPL and that only alcohol–surfactant–solvent solutions are efficient; (2) adding salts to alcohol–surfactant or to alcohol–surfactant–solvent solutions does not have a beneficial effect on DNAPL recovery; (3) washing solution formulations are site specific and must be modified if the surface properties of the solids (mineralogy) change locally, or if the interfacial behavior of liquids (type of oil) changes; (4) high solvent concentrations in washing solutions increase DNAPL extraction but also increase their cost and decrease their density dramatically; (5) maximum DNAPL recovery is observed with alcohol–surfactant–solvent formulations which correspond to the maximum solubilization in Zone C of the phase diagram; (6) replacing part of surfactant SAS by the alcohol n-butanol increases washing solution efficiency and decreases the density and the cost of solutions; (7) replacing part of n-butanol by the nonionic surfactant HOES decreases DNAPL recovery and increases the cost of solutions; (8) toluene is a better solvent than D-limonene because it increases DNAPL recovery and decreases the cost of solutions; (9) optimal alcohol–surfactant–solvent solutions contain a mixture of solvents in a mass ratio of toluene to D-limonene of one or two. Injection of 1.5 pore volumes of the optimal washing solution of n-butanol–SAS–toluene–D-limonene in water can recover up to 95% of Mercier DNAPL in sand columns. In the first pore volume of the washing solution recovered in the sand column effluent, the DNAPL is in a water-in-oil microemulsion lighter than the excess aqueous phase (Winsor Type II system), which indicates that part of the DNAPL was mobilized. In the next pore volumes, DNAPL is dissolved in a oil-in-water microemulsion phase and is mobilized in an excess oil phase lighter than the microemulsion (Winsor Type I system). The main drawback of this oil extraction process is the high concentration of ingredients necessary for DNAPL dissolution, which makes the process expensive. Because mobilization of oil seems to occur at the washing solution front, an injection strategy must be developed if there is no impermeable limit at the aquifer base. DNAPL recovery in the field could be less than observed in sand columns because of a smaller sweep efficiency related to field sand heterogeneities. The role of each component in the extraction processes in sand column as well as the Winsor system type have to be better defined for modeling purposes. Injection strategies must be developed to recover ingredients of the washing solution that can remain in the soil at the end of the washing process. ©1997 Elsevier Science B.V.     

One-dimensional simulation of solute transfer in saturated–unsaturated porous media using the discontinuous finite elements method

A one-dimensional transport model for simulating water flow and solute transport in homogeneous–heterogeneous, saturated–unsaturated porous media is presented. The model is composed of a combination of accurate numerical algorithms for solving the nonlinear Richard's and advection–dispersion equations (ADE). The mixed form of Richard's equation is solved using a standard finite element method (FEM) with primary variable switching. The transport equation is solved using operator splitting, with the discontinuous finite element method (DFE) for discretization of the advective term. A slope limiting procedure for DFE avoids numerical instabilities but creates very limited numerical dispersion for high Peclet numbers. An implicit finite differences scheme (FD) is used for the dispersive term.The unsaturated flow and transport model (Wamos-T) is applied to a variety of rigorous problems including transient flow, heterogeneous medium and abrupt variations of velocity in magnitude and direction due to time-varying boundary conditions. It produces accurate and mass-conservative solutions for a very large range of grid Peclet numbers. The Wamos-T model is a good and robust alternative for the simulation of mass transport in unsaturated domain.     

Multispectral image analysis method to determine dynamic fluid saturation distribution in two-dimensional three-fluid phase flow laboratory experiments

The need for measuring dynamic fluid saturation distribution in multi-dimensional three-fluid phase flow experiments is hampered by lack of appropriate techniques to monitor full field transient flow phenomena. There is no conventional technique able to measure dynamic three-fluid phase saturation at several array points of the flow field at the same time. A multispectral image analysis technique was developed to determine dynamic NAPL, water and air saturation distribution in two-dimensional three-fluid phase laboratory experiments. Using a digital near-infrared camera, images of sand samples with various degrees of NAPL, water and air saturation were taken, under constant lighting conditions and within three narrow spectral bands of the visible and near-infrared spectrum. It was shown that the optical density defined for the reflected luminous intensity was a linear function of the NAPL and the water saturation for each spectral band and for any two and three-fluid phase systems. This allowed the definition of dimensionless lump reflection coefficients for the NAPL and the water phase within each spectral band. Consequently, at any given time, two images taken within two different spectral bands provided two linear equations which could be solved for the water and the NAPL saturation. The method was applied to two-dimensional three-phase flow experiments, which were conducted to investigate the migration and the distribution of LNAPL in the vadose zone. The method was used to obtain continuous, quantitative and dynamic full field mapping of the NAPL saturation as well as the variation of the water and the air saturation during NAPL flow. The method provides a non-destructive and non-intrusive tool for studying multiphase flow for which rapid changes in fluid saturation in the entire flow domain is difficult to measure using conventional techniques.     

Geostatistical analysis of soil hydrologic properties in a field plot

A 91-m transect was set up in an irrigated field near Las Cruces, New Mexico to obtain soil water tension and water content data to investigate their spatial variability. A total of 455 sampling points were monitored along a grid consisting of 91 stations placed 1 m apart by 5 depths per station (at 0.3, 0.6, 0.9, 1.2 and 1.5 m below the surface). Post-irrigation tension and wetness measurements were recorded over 45 days at 11 time periods. Soil water tension was measured with tensiometers using a hand-held pressure transducer. A neutron probe was used to obtain volumetric water content. Using the observed wetness and tension data, unsaturated hydraulic conductivity values were derived (using a cubic spline function to estimate the gradient), and an exponential model was used to fit the calculated conductivity-tension curves to obtain hydraulic conductivity parameter values. The spatial and temporal variability of wetness, tension, saturated hydraulic conductivity and pore-size distribution parameters, and texture at the 0.3-m depth were examined using geostatistical techniques.The exponential model was found to inadequately describe the hydraulic conductivity/tension relationship for the full range of tension, particularly in the tension range near saturation. The derived values of the saturated hydraulic conductivity parameter were much greater than expected and do not correspond to reasonable saturated hydraulic conductivity values.All of the soil parameters studied exhibited large spatial variability horizontally and vertically in the field. Ranges of dependence determined from semivariogram analysis over the 44-day drainage period are 3–32 m for wetness, 6–34 m for soil water tension, 5–35 m for natural log of saturated hydraulic conductivity parameters, 5–11 m for pore-size distribution parameter, and 8–24 m for percent sand, silt and clay at the 0.3 m depth. An alternate hole-effect model is suggested to describe the texture semivariograms.It was determined that the variance of volumetric water content generally increased at each depth over the measured time periods, which is consistent with certain past field studies and a stochastic analysis of unsaturated flow in heterogeneous soils. Future research is recommended relating soil texture to soil hydrologic parameters with the goal of predicting soil behavior with less extensive sampling schemes.     

Detection of organic pollutants in sandy soils via TDR and eigendecomposition

In this study, organic pollutants in sandy soils are detected using a time domain relectrometry (TDR) probe system and eigendecomposition technique. As a demonstration, four types of organic fluids (motor oil, diesel oil, methanol and ethanol) were examined. Samples were prepared with different combinations of deionized water and organic fluid contents. For each experiment, reflected signals were captured by an oscilloscope, and their characteristics were identified using eigendecomposition technique. Autoregressive modeling and singular value decomposition methods were utilized for calculating the eigenvalues. The most significant eigenvalues were then identified based on their relative magnitude. Experimental results indicated that the presented system is sensitive to both water and pore fluid organic contents. For saturated conditions, signature curves were determined for identification of organic and/or water contents in soil pore fluids.     

Immiscible fluid flow in porous media: dielectric properties

The movement and retention of immiscible fluids in porous media was determined using a microwave dielectric measuring technique. The technique was used to monitor the composite dielectric constant of columns containing model soil materials to observe the effect of displacement of various saturating fluids by water and determine the retained immiscible components. Various combinations of displacing fluids, both miscible and immiscible, were studied, yielding time-dependent flow characteristics as well as the total retention factor. The flow circumstances were similar to what would be encountered in an accidental spill of chlorocarbons or crude oil.The retention of chlorocarbons and crude oil was determined for model soils consisting of sand or glass spheres. The storage capacity of the sand for immiscible fluids was dependent on the fluid initially saturating the porous medium. For chlorocarbons, it was increased in water-wet compared to dry sands while the inverse was true for crude oil. The retention of crude oil in the sand was dependent on the water flow velocity and was much larger than that for chlorocarbons. The crude oil could be displaced by solvents such as cyclohexane in the sand but the net volume occupied by the two immiscible fluids after water flushing did not change appreciably, i.e. the crude oil component was replaced by the solvent. Application of a miscible fluid-water mixture or surfactants improved the ability of water to displace the immiscible fluid from the sand column.     

Three-dimensional density-dependent flow and multicomponent reactive transport modeling of chlorinated solvent oxidation by potassium permanganate

A popular method for the treatment of aquifers contaminated with chlorinated solvents is chemical oxidation based on the injection of potassium permanganate (KMnO₄). Both the high density (1025 gL^− 1) and reactivity of the treatment solution influence the fate of permanganate (MnO₄) in the subsurface and affect the degree of contaminant treatment. The MIN3P multicomponent reactive transport code was enhanced to simulate permanganate-based remediation, to evaluate the pathways of MnO₄ utilization, and to assess the role of density contrasts for the delivery of the treatment solution. The modified code (MIN3P-D) provides a direct coupling between density-dependent fluid flow, solute transport, contaminant treatment, and geochemical reactions. The model is used to simulate a field trial of TCE oxidation in a sandy aquifer that is underlain by an aquitard. Three-dimensional simulations are conducted for a coupled reactive system comprised of ten aqueous components, two mineral phases, TCE (dissolved, adsorbed, and NAPL), reactive organic matter, and including ion exchange reactions. Model parameters are constrained by literature data and a detailed data set from the field site under investigation. The general spatial and transient evolution in observed concentrations of the oxidant, dissolved TCE, and reaction products are adequately reproduced by the simulations. The model elucidates the important role of density-induced flow and transport on the distribution of the treatment solution into NAPL containing regions located at the aquifer–aquitard interface. Model results further suggest that reactions that do not directly affect the stability of MnO₄ have a negligible effect on solution density and MnO₄ delivery.     

Sorption of polar and nonpolar organic contaminants by oil-contaminated soil

Sorption of nonpolar (phenanthrene and butylate) and polar (atrazine and diuron) organic chemicals to oil-contaminated soil was examined to investigate oil effects on sorption of organic chemicals and to derive oil–water distribution coefficients (K_oil). The resulting oil-contaminated soil–water distribution coefficients (K_d) for phenanthrene demonstrated sorption-enhancing effects at both lower and higher oil concentrations (C_oil) but sorption-reducing (competitive) effects at intermediate C_oil (approximately 1 g kg⁻¹). Rationalization of the different dominant effects was attempted in terms of the relative aliphatic carbon content which determines the accessibility of the aromatic cores to phenanthrene. Little or no competitive effect occurred for butylate because its sorption was dominated by partitioning. For atrazine and diuron, the changes in K_d at C_oil above approximately 1 g kg⁻¹ were negligible, indicating that the presently investigated oil has little or no effect on the two tested compounds even though the polarity of the oil is much less than soil organic matter (SOM). Therefore, specific interactions with the active groups (aromatic and polar domains) are dominantly responsible for the sorption of polar sorbates, and thus their sorption is controlled by available sorption sites. This study showed that the oil has the potential to be a dominant sorptive phase for nonpolar pollutants when compared to SOM, but hardly so for polar compounds. The results may aid in a better understanding of the role of the aliphatic and aromatic domains in sorption of nonpolar and polar organic pollutants.     

A field experiment with variable-suction multi-compartment samplers to measure the spatio-temporal distribution of solute leaching in an agricultural soil

Solutes spread out in time and space as they move downwards from the soil surface with infiltrating water. Solute monitoring in the field is often limited to observations of resident concentrations, while flux concentrations govern the movement of solutes in soils. A recently developed multi-compartment sampler is capable of measuring fluxes at a high spatial resolution with minimal disturbance of the local pressure head field. The objective of this paper is to use this sampler to quantify the spatial and temporal variation of solute leaching below the root zone in an agricultural field under natural rainfall in winter and spring. We placed two samplers at 31 and 25 cm depth in an agricultural field, leaving the soil above undisturbed. Each sampler contained 100 separate cells of 31 × 31 mm. Water fluxes were measured every 5 min for each cell. We monitored leaching of a chloride pulse under natural rainfall by frequently extracting the collected leachate while leaving the samplers buried in situ. This experiment was followed by a dye tracer experiment. This setting yielded information that widely surpassed the information that can be provided by separate anionic and dye tracer trials, and solute transport monitoring by coring or suction cups. The detailed information provided by the samplers showed that percolation at the sampling depth started much faster (approximately 3 h after the start of rainfall) in initially wet soil (pressure head above − 65 cm) than in drier soil (more than 14 h at pressure heads below − 80 cm). At any time, 25% of the drainage passed through 5–6% of the sampled area, reflecting the effect of heterogeneity on the flow paths. The amount of solute carried by individual cells varied over four orders of magnitude. The lateral concentration differences were limited though. This suggests a convective–dispersive regime despite the short vertical travel distance. On the other hand, the dilution index indicates a slight tendency towards stochastic–convective transport at this depth. There was no evidence in the observed drainage patterns and dye stained profiles of significant disturbance of the flow field by the samplers.     

Integrated biological responses of zebrafish (Danio rerio) to analyze water quality in regions under anthropogenic influence

This study analyzed water quality in regions around Patos lagoon (Southern Brazil) that are under anthropogenic pressure. Water samples were collected from five different sites, including one used as a source for human consumption (COR) and others known to be influenced by human activities (IP). Danio rerio (Teleostei, Cyprinidae) organisms were exposed for 24 h to these water samples, plus a control group. It was observed that: (1) reactive oxygen species levels were lower in COR and IP than in the control group; (2) glutamate-cysteine ligase (catalytic subunit) expression was higher in COR than in other sites; (3) exposure to all water samples affected long-term memory (LTM) when compared to control group. Thus, some water samples possess the ability to modulate the antioxidant system and to induce a decline in cognitive functions, as measured by LTM. The obtained results indicate that a combination of variables of different organization level (molecular, biochemical and behavioral) can be employed to analyze water quality in impacted regions.     

Thermal venting to recover less-volatile hydrocarbons from the unsaturated zone, 1. Theory

Thermal venting is a remediation technique suitable to the liquid unsaturated zone to enhance recovery of less-volatile residual hydrocarbon contaminants. Thermal venting is different to traditional soil venting because heated air instead of air at ambient conditions is applied to the contaminated zone. The vapor pressure of a less-volatile contaminants is typically increased by temperature causing the gas-phase concentrations to increase by three- to five-fold over a temperature increase of 20–30°C. The work described in this first paper provides the theoretical framework of analysis related to thermal venting. The analysis included nonisothermal gas flow, thermal energy transport and multicomponent mass transport in a multiphase porous medium. The transient gas flow analysis included the effect of temperature on fluid properties and gas compressibility. The heat energy transport analysis was performed under the thermodynamic equilibrium condition with phase-summed effective thermal properties. Multi-component mass transport was performed under local equilibrium for partitioning between phases. Model verification was performed to the extend possible using analytical and available experimental data for different physical processes. The second paper of this two-part series will demonstrate the applicability of thermal venting technique through numerical simulations of hypothetical laboratory and field-scale scenarios.     

Modelling tritium and phosphorus transport by preferential flow in structured soil

Subsurface solute transport through structured soil is studied by model interpretation of experimental breakthrough curves from tritium and phosphorus tracer tests in three intact soil monoliths. Similar geochemical conditions, with nearly neutral pH, were maintained in all the experiments. Observed transport differences for the same tracer are thus mainly due to differences in the physical transport process between the different monoliths. The modelling is based on a probabilistic Lagrangian approach that decouples physical and chemical mass transfer and transformation processes from pure and stochastic advection. Thereby, it enables explicit quantification of the physical transport process through preferential flow paths, honouring all independently available experimental information. Modelling of the tritium breakthrough curves yields a probability density function of non-reactive solute travel time that is coupled with a reaction model for linear, non-equilibrium sorption–desorption to describe the phosphorus transport. The tritium model results indicate that significant preferential flow occurs in all the experimental soil monoliths, ranging from 60–100% of the total water flow moving through only 25–40% of the total water content. In agreement with the fact that geochemical conditions were similar in all experiments, phosphorus model results yield consistent first-order kinetic parameter values for the sorption–desorption process in two of the three soil monoliths; phosphorus transport through the third monolith cannot be modelled because the apparent mean transport rate of phosphorus is anomalously rapid relative to the non-adsorptive tritium transport. The occurrence of preferential flow alters the whole shape of the phosphorus breakthrough curve, not least the peak mass flux and concentration values, and increases the transported phosphorus mass by 2–3 times relative to the estimated mass transport without preferential flow in the two modelled monoliths.     

基于剪切乳化机的采油污水再现配置与分散相破碎乳化行为

以剪切乳化机为主体,设计、加工了一套采油污水室内再现配置装置,探讨了含油量、剪切强度(与流量同比变化,下同)等不同配置条件下的再现采油污水粒度分布、乳化度、Zeta电位和粘度等性能参数以及乳化机前后的压力损失.结果表明,该采油污水再现配置装置,能满足不同性质采油污水的配置要求,剪切乳化机的压力损失随着剪切强度的增大而明显升高,而随含油量增大没有明显的变化趋势,呈现出基本稳定的状态.较大的剪切强度创造了湍流流体条件,强化了分散相颗粒之间的碰撞聚结行为,因而再现采油污水的乳化度随剪切强度增大而降低,这种降低的程度又受含油量的影响:低含油量时剪切强度的影响明显,高含油量时这种影响逐渐减弱并趋于平缓.鉴于高含油量时趋于形成大粒径油珠,再现采油污水的乳化度随含油量增大而呈现略微的下降趋势,这种下降趋势的程度同时受剪切强度的影响:低剪切强度时含油量的影响较为明显,而高剪切强度时的剪切条件削弱了含油量的影响.Zeta电位随着剪切强度增大而不断增大,使得分散相颗粒带电绝对值降低,乳化体系有脱稳的趋势,导致颗粒间的碰撞聚结机会变大,与乳化度降低的趋势相吻合.剪切强度和含油量对再现采油污水粘度的影响不大.     

Single- and dual-porosity modelling of flow in reclaimed mine soil cores with embedded lignitic fragments

Lignitic mine soils represent a typical two-scale dual-porosity medium consisting of a technogenic mixture of overburden sediments that include lignitic components as dust and as porous fragments embedded within a mostly coarse-textured matrix. Flow and transport processes in such soils are not sufficiently understood to predict the course of soil reclamation or of mine drainage. The objective of this contribution is to identify the most appropriate conceptual model for describing small-scale heterogeneity effects on flow on the basis of the physical structure of the system. Multistep flow experiments on soil cores are analyzed using either mobile–immobile or mobile–mobile type 1D dual-porosity models, and a 3D numerical model that considers a local-scale distribution of fragments. Simulations are compared with time series' of upward infiltration and matric potential heads measured at two depths using miniature tensiometers. The 3D and the 1D dual-permeability models yielded comparable results as long as pressure heads are in local equilibrium; however, could describe either the upward infiltration or the matric potential curves but not both at the same time. The mobile–immobile type dual-porosity model failed to describe the data. A simultaneous match with pressure heads and upward infiltration data could only be obtained with the 1D dual-permeability model (i.e., mobile–mobile) by assuming an additional restriction of the inter-domain water transfer. These results indicate that for unsaturated flow conditions at higher matric potential heads (i.e., here >− 40 hPa), water in a restricted part of the fragment domain must be more mobile as compared to water in the sandy matrix domain. Closer inspections of the pore system and first neutron radiographic imaging support the hypothesis that a more continuous pore region exists at these pressure heads in the vicinity of the lignitic fragments possibly formed by fragment contacts and a lignitic dust interface-region between the two domains. The results suggest that the small-scale structure is too complex as to be represented by weighted contributions of individual components alone.     

Macro-porosity and leaching of atrazine in tilled and orchard loamy soils

Atrazine is the most commonly detected herbicide in the groundwater. Leaching of atrazine largely depends on soil management practices. The aim of this study was to examine leaching of atrazine in tilled and orchard silty loam soils. The experimental objects included: conventionally tilled field (CT) with main tillage operations including pre-plow (10 cm) + harrowing, mouldboard ploughing (20 cm), and a 35 year-old apple orchard (OR) with a permanent sward. To determine leaching of atrazine soil columns of undisturbed structure were taken with steel cylinders of 21.5 cm diameter and 20 cm high from the depth of 0–20 cm. All columns were equilibrated at water content corresponding to field capacity (0.21 kg kg⁻¹). Atrazine suspended in distilled water was dripped uniformly onto the surface of each column. Then water was infiltrated and breakthrough times of leachates were recorded. Atrazine concentration in the leachates was determined by means of HPLC Waters. Macro-porosity and percolation rate were higher in OR than CT soil. Cumulative recovery % of the atrazine applied was 1.267% for OR and approximately one third more from the CT soil but the rate of leaching (per unit of time) was greater from the OR soil. The lower leaching under OR than CT can be due to a greater SOM and the presence of earthworm burrows with organic burrow linings that could adsorb atrazine and contribute to preferential flow allowing solutes to bypass parts whereas the greater rate of leaching due to a greater infiltration rate.The results indicate potential of management practices for minimizing atrazine leaching.     

Hydrological and geochemical factors affecting leachate composition in municipal solid waste incinerator bottom ash: Part I: The hydrology of Landfill Lostorf, Switzerland

The objective of the investigation of the municipal solid waste incinerator (MSWI) bottom ash landfill, Landfill Lostorf, was to determine the residence time of water in the landfill and the flow paths through the landfill. Over a period of 22 months, measurements of rainfall, landfill discharge and leachate electrical conductivity were recorded and tracer experiments made. Over the yearly period 1995, approximately 50% of the incident rainfall was measured in the discharge. An analysis of single rain events showed that in winter, 90–100% of rainfall was expressed in the landfill discharge, whereas in summer months, the value was between 9 and 40% depending on the intensity of the rain event. The response to rainfall was rapid. Within 30–100 h, approximately 50% of water discharged in response to a rain event had left the landfill. The discharge was less than 4 l/min for approximately 50% of the measurement periods. Qualitative tracer studies with fluorescein, pyranine and iodide clearly showed the existence of preferential flow paths. This was further substantiated by quantitative tracer studies of single rain events using / ratios and electrical conductivity measurements. The proportion of rainwater passing directly through the landfill was found to be between 20 and 80% in summer months and around 10% in winter months. The difference has been ascribed to the water content in the landfill. The average residence time of the water within the landfill has been estimated to be roughly 3 years and this water is the predominant component in the discharge over a yearly period.     

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.