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.     

Preliminary 3-D site-scale studies of radioactive colloid transport in the unsaturated zone at Yucca Mountain,Nevada

The U.S. Department of Energy (DOE) is actively investigating the technical feasibility of permanent disposal of high-level nuclear waste in a repository to be situated in the unsaturated zone (UZ) at Yucca Mountain (YM), Nevada. In this study we investigate, by means of numerical simulation, the transport of radioactive colloids under ambient conditions from the potential repository horizon to the water table. The site hydrology and the effects of the spatial distribution of hydraulic and transport properties in the Yucca Mountain subsurface are considered. The study of migration and retardation of colloids accounts for the complex processes in the unsaturated zone of Yucca Mountain, and includes advection, diffusion, hydrodynamic dispersion, kinetic colloid filtration, colloid straining, and radioactive decay. The results of the study indicate that the most important factors affecting colloid transport are the subsurface geology and site hydrology, i.e., the presence of faults (they dominate and control transport), fractures (the main migration pathways), and the relative distribution of zeolitic and vitric tuffs. The transport of colloids is strongly influenced by their size (as it affects diffusion into the matrix, straining at hydrogeologic unit interfaces, and transport velocity) and by the parameters of the kinetic-filtration model used for the simulations. Arrival times at the water table decrease with an increasing colloid size because of smaller diffusion, increased straining, and higher transport velocities. The importance of diffusion as a retardation mechanism increases with a decreasing colloid size, but appears to be minimal in large colloids.     

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.     

非饱和土壤轻油污染多相流研究进展

轻油泄漏到地下后 ,在非饱和土壤中水、气、油三相相互作用、相互影响 ,共同构成多相流体系统。在简要分析轻油污染问题的产生、危害以及轻油泄漏在地表以下运移方式的基础上 ,介绍并总结了近些年来国内外土壤轻油污染有关的多相流研究方面的进展 ,包括国外的模型研究成果、国内外试验研究方法和研究成果 ,并指出目前研究中存在的问题和需要进一步探讨的方向     

Transport and degradation of propyleneglycol and potassium acetate in the unsaturated zone

De-icing chemicals used during the winter season are potential pollutants for the groundwater underneath the new main airport of Norway. Several field experiments examining the transport and degradation of propyleneglycol (PG), potassium acetate (KAc) and non-reactive tracers were performed in a lysimeter trench under natural snowmelting conditions. Chemicals were applied underneath the snow cover and the transport in a heterogeneous coarse sandy soil was examined by extracting soil water from 30 or 40 suction cups placed at five depths between 0.4 and 2.4 m depth. Transport and degradation was analysed by spatial moment calculations. The de-icing chemicals showed the same basic displacement as chemically inactive tracers, an initial fast transport during the melting period followed by a period of stagnation throughout the summer season. PG seemed to be displaced to greater depths compared to non-reactive tracer after the first application. However, computer simulations of transport and degradation in a heterogeneous unsaturated soil showed that decreasing degradation constants with depth can generate a downward movement of the centre of mass without any flow occurring in the system. Potassium acetate showed some adsorption, with calculated retardation factors of approximately 1.3 and 1.2. The degradation rate constant for PG was calculated to be 0.015 day-1 in 1994 and increased to 0.047 day-1 in the second application made in 1995. The degradation rate constant for acetate was estimated to be 0.02 day-1. Increased manganese concentrations seem to be a good indicator of degradation of PG and Ac.     

Design and laboratory testing of a chamber device to measure total flux of volatile organic compounds from the unsaturated zone under natural conditions

To determine if an aquifer contaminated with volatile organic compounds (VOCs) has potential for natural remediation, all natural processes affecting the fate and transport of VOCs in the subsurface must be identified and quantified. This research addresses the quantification of air-phase volatile organic compounds (VOCs) leaving the unsaturated zone soil gas and entering the atmosphere-including the additional flux provided by advective soil-gas movement induced by barometric pumping. A simple and easy-to-use device for measuring VOC flux under natural conditions is presented. The vertical flux chamber (VFC) was designed using numerical simulations and evaluated in the laboratory. Mass-balance numerical simulations based on continuously stirred tank reactor equations (CSTR) provided information on flux measurement performance of several sampling configurations with the final chamber configuration measuring greater than 96% of model-simulated fluxes. A laboratory device was constructed to evaluate the flux chamber under both diffusion-only and advection-plus-diffusion transport conditions. The flux chamber measured an average of 82% of 15 diffusion-only fluxes and an average of 95% of 15 additional advection-plus-diffusion flux experiments. The vertical flux chamber has the capability of providing reliable measurement of VOC flux from the unsaturated zone under both diffusion and advection transport conditions.     

An integrated approach for assessing the migration behavior of chlorpyrifos and carbaryl in the unsaturated soil zone

Chlorpyrifos (O, O-diethyl O-3,5,6-trichloropyridin-2-yl phosphorothioate) and carbaryl (1-naphthyl methylcarbamate) are often applied concurrently as insecticides in food production. The aim of this study was to research their migration behavior in a real environment. We researched the leaching of both pesticides by setting up field lysimeters on a farm with the typical soil used in fruit production today. In order to analyze the variables involved in this process, we performed complementary adsorption studies, we performed complementary adsorption studies using batches and undisturbed soil laboratory columns for both compounds. The results for pesticide transport through the lysimeters showed that less than 1% of chlorpyrifos was recovered in the leachates, while almost 17% was recovered for carbaryl. Having completed the experiment in undisturbed laboratory columns, soil analysis showed that chlorpyrifos mainly remained in the first 5 cm, while carbaryl moved down to the lower sections. These results can be explained in view of the sorption coefficient values (K_D) obtained in horizons A and B for chlorpyrifos (393 and 184 L kg^?1) and carbaryl (3.1 and 4.2 L kg^?1), respectively. By integrating the results obtained in the different approaches, we were able to characterize the percolation modes of these pesticides in the soil matrix, thus contributing to the sustainable use of resources.     

Oxygenated gasoline release in the unsaturated zone, Part 2: Downgradient transport of ethanol and hydrocarbons

In the event of a gasoline spill containing oxygenated compounds such as ethanol and MTBE, it is important to consider the impacts these compounds might have on subsurface contamination. One of the main concerns commonly associated with ethanol is that it might decrease the biodegradation of aromatic hydrocarbon compounds, leading to an increase in the hydrocarbon dissolved plume lengths. The first part of this study (Part 1) showed that when gasoline containing ethanol infiltrates the unsaturated zone, ethanol is likely to partition to and be retained in the unsaturated zone pore water. In this study (Part 2), a controlled field test is combined with a two-dimensional laboratory test and three-dimensional numerical modelling to investigate how ethanol retention in the unsaturated zone affects the downgradient behaviour of ethanol and aromatic hydrocarbon compounds. Ethanol transport downgradient was extremely limited. The appearance of ethanol in downgradient wells was delayed and the concentrations were lower than would be expected based on equilibrium dissolution. Oscillations in the water table resulted in minor flushing of ethanol, but its effect could still be perceived as an increase in the groundwater concentrations downgradient from the source zone. Ethanol partitioning to the unsaturated zone pore water reduced its mass fraction within the NAPL thus reducing its anticipated impact on the fate of the hydrocarbon compounds. A conceptual numerical simulation indicated that the potential ethanol-induced increase in benzene plume length after 20 years could decrease from 136% to 40% when ethanol retention in the unsaturated zone is considered.     

Illicit crops and armed conflict as constraints on biodiversity conservation in the Andes region

Coca, once grown for local consumption in the Andes, is now produced for external markets, often in areas with armed conflict. Internationally financed eradication campaigns force traffickers and growers to constantly relocate, making drug-related activities a principal cause of forest loss. The impact on biodiversity is known only in general terms, and this article presents the first regional analysis to identify areas of special concern, using bird data as proxy. The aim of conserving all species may be significantly constrained in the Santa Marta and Perijá mountains, Darién, some parts of the Central Andes in Colombia, and between the middle Mara?ón and middle Huallaga valleys in Peru. Solutions to the problem must address the root causes: international drug markets, long-lasting armed conflict, and lack of alternative income for the rural poor.     

Strontium isotope constraints on the seasonal variation of the provenance of base cations in rain water at Kawakami, Central Japan

The ⁸⁷Sr/⁸⁶Sr ratios in monthly precipitation in the forested basin at Kawakami, central Japan, varied seasonally from 0.709 to 0.711 in spring to as low as 0.7062 ± 0.0004 in autumn over nine years from 1987 to 1995. The seasonal variation can be explained in terms of the mixing of three sources of Sr: sea salt (⁸⁷Sr/⁸⁶Sr 0.70917), soluble eolian minerals originating from deserts in continental Asia ( 0.711), and biogenic materials growing on soils derived from the volcanic rock substrates in the vicinity ( 0.706). It is estimated that the contribution of sea-salt Sr into the Kawakami rain Sr is relatively constant (10 ± 5%) and that more than half of the Sr is of biogenic origin throughout the year except during spring rains when 50% of the Sr is due to the dissolution of Ca-minerals from Asian dusts. The dominant contribution of eolian components to spring rains is consistent with the high concentration of ³H, which is typical of air masses from Asia.     

有机污染物在包气带中迁移转化试验研究

采集一定浓度的有机物废水 ,在试验室内进行了静态吸附、静态降解和动态土柱试验 ,对COD在包气带中迁移转化规律进行研究 ,提出了描述COD在包气带中迁移转化规律的数学模型。结果表明 :包气带对COD的吸附过程是线性的 ,可用亨利吸附模式s =Kdc +s0 表示 ,吸附系数Kd=0 .0 6 93;包气带对COD的降解曲线基本符合一级动力学方程c =c0 e-k1 t,降解系数k1=0 .0 4 99d-1;弥散试验测得弥散系数D =0 .0 0 2 4 2m2 /d。COD在包气带中的迁移转化过程是弥散、吸附、降解等多种作用共同作用的结果。     

Quantifying ozone uptake at the canopy level of spruce,pine and larch trees at the alpine timberline: an approach based on sap flow measurement

Micro-climatic and ambient ozone data were combined with measurements of sap flow through tree trunks in order to estimate whole-tree ozone uptake of adult Norway spruce (Picea abies), cembran pine (Pinus cembra), and European larch (Larix decidua) trees. Sap flow was monitored by means of the heat balance approach in two trees of each species during the growing season of 1998. In trees making up the stand canopy, the ozone uptake by evergreen foliages was significantly higher than by deciduous ones, when scaled to the ground area. However, if expressed per unit of whole-tree foliage area, ozone flux through the stomata into the needle mesophyll was 1.09, 1.18 and 1.40 nmol m(-2) s(-1) in Picea abies, Pinus cembra and Larix decidua, respectively. These fluxes are consistent with findings from measurements of needle gas exchange, published from the same species at the study site. It is concluded that the sap flow-based approach offers an inexpensive, spatially and temporally integrating way for estimating ozone uptake at the whole-tree and stand level, intrinsicly covering the effect of boundary layers on ozone flux.     

Colloid transport in unsaturated porous media: the role of water content and ionic strength on particle straining

Packed column and mathematical modeling studies were conducted to explore the influence of water saturation, pore-water ionic strength, and grain size on the transport of latex microspheres (1.1 microm) in porous media. Experiments were carried out under chemically unfavorable conditions for colloid attachment to both solid-water interfaces (SWI) and air-water interfaces (AWI) using negatively charged and hydrophilic colloids and modifying the solution chemistry with a bicarbonate buffer to pH 10. Interaction energy calculations and complementary batch experiments were conducted and demonstrated that partitioning of colloids to the SWI and AWI was insignificant across the range of the ionic strengths considered. The breakthrough curve and final deposition profile were measured in each experiment indicating colloid retention was highly dependent on the suspension ionic strength, water content, and sand grain size. In contrast to conventional filtration theory, most colloids were found deposited close to the column inlet, and hyper-exponential deposition profiles were observed. A mathematical model, accounting for time- and depth-dependent straining, produced a reasonably good fit for both the breakthrough curves and final deposition profiles. Experimental and modeling results suggest that straining--the retention of colloids in low velocity regions of porous media such as grain junctions--was the primary mechanism of colloid retention under both saturated and unsaturated conditions. The extent of stagnant regions of flow within the pore structure is enhanced with decreasing water content, leading to a greater amount of retention. Ionic strength also contributes to straining, because the number of colloids that are held in the secondary energy minimum increases with ionic strength. These weakly associated colloids are prone to be translated to stagnation regions formed at grain-grain junctions, the solid-water-air triple point, and dead-end pores and then becoming trapped.     

Seasonal and diurnal measurements of carbon monoxide and nonmethane hydrocarbons at Mt. Wilson,California: Indirect evidence of atomic Cl in the Los Angeles basin

We present a study of the seasonal and diurnal variability of carbon monoxide and selected volatile organic compounds in the Los Angeles area. Measurements were made during four different nine-day field campaigns in April/May, September, and November, 2007, and February, 2008, at the Mt. Wilson sampling site, which is located at an elevation of approximately 1700 m in the San Gabriel Mountains overlooking Pasadena and the Los Angeles basin. The results were used to characterize the Mt. Wilson site as a representative location for monitoring integrated Los Angeles basin emissions, and, by reference to carbon monoxide emissions, to estimate average annual emissions. The considerable seasonal variability of many hydrocarbons, in both their measured mixing ratios and their relationship to carbon monoxide, was indicative of variable source strengths. Most interestingly, perturbation of C₄ hydrocarbon ratios suggested an enhanced role for chlorine chemistry during the month of September, likely as the result of Los Angeles’ coastal location. Such coastal influence was confirmed by observations of enhanced mixing ratios of marine halocarbons, as well as air mass back trajectories.     

Movement and degradation of aldicarb residues in the saturated zone under citrus groves on the Florida ridge

A 1.7-ha section of citrus grove near Lake Hamilton was the site of a three-year field study designed to monitor the movement and degradation of the nematicide and insecticide aldicarb in the central ridge area of Florida. Soil cores were used to monitor the fate of aldicarb residues in the unsaturated zone and over 2,000 groundwater samples were collected from 174 monitoring wells to measure horizontal and vertical transport of aldicarb residues in the saturated zone. A simple saturated zone model was used to estimate the degradation rate of aldicarb residues and extrapolate findings to other ridge areas.The results of the study suggest that in the saturated zone aldicarb residues degrade at a rate corresponding to a half-life of approximately eight months. The predominantly horizontal movement of groundwater at this site limits aldicarb residues to the upper three to five meters of the saturated zone. Field data from this site together with unsaturated and saturated zone simulations suggest that in this area of Florida current restrictions on aldicarb used near potable wells are adequate to protect drinking water supplies.     

Microplastics in water,sediments, and fish at Alpine River,originating from the Hindu Kush Mountain,Pakistan: implications for conservation

Environmental Science and Pollution Research - Microplastics (MP) pollution is an emerging threat to life and the environment. These particles are not restricted to human-inhabited lands but also...     

Coupled effects of solution chemistry and hydrodynamics on the mobility and transport of quantum dot nanomaterials in the vadose zone

To investigate the coupled effects of solution chemistry and hydrodynamics on the mobility of quantum dot (QD) nanoparticles in the vadose zone, laboratory scale transport experiments involving single and/or sequential infiltrations of QDs in unsaturated and saturated porous media, and computations of total interaction and capillary potential energies were performed. As ionic strength increased, QD retention in the unsaturated porous media increased; however, this retention was significantly suppressed in the presence of a non-ionic surfactant in the infiltration suspensions as indicated by surfactant enhanced transport of QDs. In the vadose zone, the non-ionic surfactant limited the formation of QD aggregates, enhanced QD mobility and transport, and lowered the solution surface tension, which resulted in a decrease in capillary forces that not only led to a reduction in the removal of QDs, but also impacted the vadose zone flow processes. When chemical transport conditions were favorable (ionic strength of 5 × 10(-4)M and 5 × 10(-3)M, or ionic strengths of 5 × 10(-2)M and 0.5M with surfactant), the dominating phenomena controlling the mobility and transport of QDs in the vadose zone were meso-scale processes, where infiltration by preferential flow results in the rapid transport of QDs. When chemical transport conditions were unfavorable (ionic strength of 5 × 10(-2)M and 0.5M) the dominating phenomena controlling the mobility and transport of QDs in the vadose zone were pore-scale processes governed by gas-water interfaces (GWI) that impact the mobility of QDs. The addition of surfactant enhanced the transport of QDs both in favorable and unfavorable chemical transport conditions. The mobility and retention of QDs was controlled by interaction and capillary forces, with the latter being the most influential. GWI were found to be the dominant mechanism and site for QD removal compared with solid-water interfaces (SWI) and pore straining. Additionally, ripening phenomena were demonstrated to enhance QDs removal or retention in porous media and to be attenuated by the presence of surfactant.     

Effects of soil moisture and physical-chemical properties of organic pollutants on vapor-phase transport in the vadose zone

Vapor-phase transport of organic pollutants is one of the important pathways in the distribution and attenuation of volatile organic compounds in the vadose zone. In this study, the impact of vapor-phase partitioning and of the physical-chemical properties of organic pollutants on vapor-phase transport was assessed. An experimentally derived relationship to predict vapor sorption for a variety of soil types under varying soil moisture conditions was incorporated into the two-dimensional finite-element model, Vocwaste. The revised model was then used to simulate the transport of volatile organics. Vapor-phase partitioning in the model accounted for vapor uptake by sorption onto moist mineral surfaces as well as sorption at the liquid-solid interface and dissolution into soil water. Under dry conditions, vapor-phase sorption of volatile organic pollutants was shown to have a retarding effect on transport of organic vapors. However, for shallow, contaminated soils, volatilization was controlled by vapor diffusion, even under dry conditions where vapor-phase sorption was high. The influence of Henry's law constant and of the aqueous-phase (solid-liquid) partition coefficient for volatile organic pollutants was considered in the simulations. Volatilization of organic vapors was shown to be favored for contaminants with high Henry's law constants and low aqueous-phase partitioning coefficients. Because of the interdependence of these two physical-chemical properties, individual properties of the contaminant should not be considered in isolation in the evaluation of vapor transport.     

Matrix diffusion coefficients in volcanic rocks at the Nevada test site: influence of matrix porosity, matrix permeability, and fracture coating minerals

Diffusion cell experiments were conducted to measure nonsorbing solute matrix diffusion coefficients in forty-seven different volcanic rock matrix samples from eight different locations (with multiple depth intervals represented at several locations) at the Nevada Test Site. The solutes used in the experiments included bromide, iodide, pentafluorobenzoate (PFBA), and tritiated water ((3)HHO). The porosity and saturated permeability of most of the diffusion cell samples were measured to evaluate the correlation of these two variables with tracer matrix diffusion coefficients divided by the free-water diffusion coefficient (D(m)/D*). To investigate the influence of fracture coating minerals on matrix diffusion, ten of the diffusion cells represented paired samples from the same depth interval in which one sample contained a fracture surface with mineral coatings and the other sample consisted of only pure matrix. The log of (D(m)/D*) was found to be positively correlated with both the matrix porosity and the log of matrix permeability. A multiple linear regression analysis indicated that both parameters contributed significantly to the regression at the 95% confidence level. However, the log of the matrix diffusion coefficient was more highly-correlated with the log of matrix permeability than with matrix porosity, which suggests that matrix diffusion coefficients, like matrix permeabilities, have a greater dependence on the interconnectedness of matrix porosity than on the matrix porosity itself. The regression equation for the volcanic rocks was found to provide satisfactory predictions of log(D(m)/D*) for other types of rocks with similar ranges of matrix porosity and permeability as the volcanic rocks, but it did a poorer job predicting log(D(m)/D*) for rocks with lower porosities and/or permeabilities. The presence of mineral coatings on fracture walls did not appear to have a significant effect on matrix diffusion in the ten paired diffusion cell experiments.