Description of hydrogeologic heterogeneity and evaluation of radionuclide transport at an underground nuclear test

Realistic models of contaminant transport in groundwater demand detailed characterization of the spatial distribution of subsurface hydraulic properties, while at the same time programmatic constraints may limit collection of pertinent hydraulic data. Fortunately, alternate forms of data can be used to improve characterization of spatial variability. We utilize a methodology that augments sparse hydraulic information (hard data) with more widely available hydrogeologic information to generate equiprobable maps of hydrogeologic properties that incorporate patterns of connected permeable zones. Geophysical and lithologic logs are used to identify hydrogeologic categories and to condition stochastic simulations using Sequential Indicator Simulation (SIS). The resulting maps are populated with hydraulic conductivity values using field data and Sequential Gaussian Simulation (SGS). Maps of subsurface hydrogeologic heterogeneity are generated for the purpose of examining groundwater flow and transport processes at the Faultless underground nuclear test, Central Nevada Test Area (CNTA), through large-scale, three-dimensional numerical modeling. The maps provide the basis for simulation of groundwater flow, while transport of radionuclides from the nuclear cavity is modeled using particle tracking methods. Sensitivity analyses focus on model parameters that are most likely to reduce the long travel times observed in the base case. The methods employed in this study have improved our understanding of the spatial distribution of preferential flowpaths at this site and provided the critical foundation on which to build models of groundwater flow and transport. The results emphasize that the impacts of uncertainty in hydraulic and chemical parameters are dependent on the radioactive decay of specific species, with rapid decay magnifying the effects of parameters that change travel time.     

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.     

Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation

This paper presents the results of a detailed field investigation that was performed for studying groundwater recharge processes and solute downward migration mechanisms prevailing in the unsaturated zone overlying a chalk aquifer in Belgium. Various laboratory measurements were performed on core samples collected during the drilling of boreholes in the experimental site. In the field, experiments consisted of well logging, infiltration tests in the unsaturated zone, pumping tests in the saturated zone and tracer tests in both the saturated and unsaturated zones. Results show that gravitational flows govern groundwater recharge and solute migration mechanisms in the unsaturated zone. In the variably saturated chalk, the migration and retardation of solutes is strongly influenced by recharge conditions. Under intense injection conditions, solutes migrate at high speed along the partially saturated fissures, downward to the saturated zone. At the same time, they are temporarily retarded in the almost immobile water located in the chalk matrix. Under normal recharge conditions, fissures are inactive and solutes migrate slowly through the chalk matrix. Results also show that concentration dynamics in the saturated zone are related to fluctuations of groundwater levels in the aquifer. A conceptual model is proposed to explain the hydrodispersive behaviour of the variably saturated chalk. Finally, the vulnerability of the chalk to contamination issues occurring at the land surface is discussed.     

Kinetic modeling of microbially-driven redox chemistry of radionuclides in subsurface environments: coupling transport, microbial metabolism and geochemistry

Microbial reactions play an important role in regulating pore water chemistry as well as secondary mineral distribution in many subsurface systems and, therefore, may directly impact radionuclide migration in those systems. This paper presents a general modeling approach to couple microbial metabolism, redox chemistry, and radionuclide transport in a subsurface environment. To account for the likely achievement of quasi-steady state biomass accumulations in subsurface environments, a modification to the traditional microbial growth kinetic equation is proposed. The conditions for using biogeochemical models with or without an explicit representation of biomass growth are clarified. Based on the general approach proposed in this paper, the couplings of uranium reactions with biogeochemical processes are incorporated into computer code BIORXNTRN Version 2.0. The code is then used to simulate a subsurface contaminant migration scenario, in which a water flow containing both uranium and a complexing organic ligand is recharged into an oxic carbonate aquifer. The model simulation shows that Mn and Fe oxyhydroxides may vary significantly along a flow path. The simulation also shows that uranium(VI) can be reduced and therefore immobilized in the anoxic zone created by microbial degradation.     

Numerical methods for describing chemical transport in the unsaturated zone of the subsurface

Finite-difference and finite-element methods of approximation have been extended to solve the one-dimensional nonlinear partial differential equations that describe the simultaneous transport of heat, moisture and chemical in the unsaturated zone. Especially for chemical transport, nodal spacing criteria are required to minimize numerical dispersion and oscillatory behavior in the solution vector for chemical concentration. Conservative criteria for nodal spacing for saturated flow can be used to set nodal spacing for unsaturated zone transport. When nodal spacing criteria are satisfied, for the same set of transport and boundary conditions, chemical concentration profiles calculated by the two numerical methods will be almost the same. A situation that is simulated very well with one-dimensional models, is the application of chemicals to land surfaces. To compare and contrast the characteristics of solutions given by the two numerical methods, moisture content, temperature and chemical concentration profiles for a 75-day period after application in the unsaturated zone are calculated for two representative types of organic chemicals. In the first, the chemical is very slowly degraded in the subsurface environment but strongly sorbed to soil surfaces. In the second, the chemical is rapidly degraded but weakly sorbed to soil surfaces. Because of differences in sorption coefficients and mechanisms of degradation, for the same set of hydrodynamic properties of the subsurface, the weakly sorbed chemical is more widely distributed throughout the unsaturated zone, whereas the strongly sorbed chemical stays very close to where it is put initially with little penetration into the subsurface. Satisfying nodal spacing criteria minimizes the impact of the method of approximation on the calculated solutions of the transport equations. For better model predictive performance, however, there are needs for more fundamental information on processes governing transport in the subsurface.     

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.     

Predicting natural attenuation of xylene in groundwater using a numerical model.

The aquifer beneath an abandoned refinery in the Lower Rhine area, Germany, was contaminated with a number of different mineral oil products. Groundwater sampling in the area around the former xylene plant revealed that a xylene plume had developed in the underlying groundwater, and moreover, that there is strong evidence for in situ microbial xylene degradation with oxygen, nitrate, sulfate and ferric iron as electron acceptors. In order to prevent further xylene spreading, three pumping wells extracting contaminated water were installed downgradient of the spill zone. The numerical reactive transport code Transport Biochemisty Chemistry (TBC) was applied to this situation to quantify the relation of microbial degradation to xylene removal by the pumping wells. It could be shown that the unamended in situ degradation was an appreciable xylene removal process that contributed to about one-third to the total xylene removal (degradation plus extraction). A further objective of the model application was to predict xylene spreading under regional flow conditions, i.e. without operation of the three pumping wells, to consider the possible effects of natural xylene attenuation. To accomplish this, the model calibrated for the situation with operating wells was transferred to the hydraulic situation of regional flow while retaining the parameters of the biochemical model. It turned out that the xylene plume that is expected to develop downgradient of the source area will be limited to an extension of not more than 1000 m. An interesting feature of the simulations results was that xylene degradation under iron-reducing conditions, which was of minor importance for the situation with operating pumping wells, becomes the dominant degradation mechanism under regional flow conditions. Moreover, iron reduction will be the key process in controlling plume evolution. The model application illustrates that multi-species reactive transport models are needed to adequately transfer reactive processes from one hydraulic situation to another, while single species models are not suited for this predictive task.     

The concentration, apparent molecular weight and chemical reactivity of silica from groundwater in southern Nevada

Sorption of radionuclides, metals and organic compounds to colloidal particles has been suggested to increase the mobility of these pollutants in groundwater. Because silicates and alumino-silicates can be important components of groundwater colloids, we have conducted a study to characterize the nature of silica in various springs and wells in Southern Nevada and to determine the extent that silica may be associated with colloidal particles that can participate in pollutant transport. The total silica content was measured using inductively coupled plasma emission spectroscopy (ICP). In addition, reactive silica was measured using the silica molybdate colorimetric technique. The apparent molecular weight of the silica was investigated using split-flow-lateral-transport-thin-cell-fractionation (SPLITT) which can readily distinguish between colloidal and low molecular weight associations. This study indicates that silica does not tend to form stable inorganic colloids in Southern Nevada groundwaters but exists as low molecular weight species. However, water from one of the test facilities on the Nevada Test Site (NTS) did contain stable siliceous colloids that could have important implications for the modeling the transport of radionuclides at this site.     

Influence of hydrological and geochemical processes on the transport of chelated metals and chromate in fractured shale bedrock

Field-scale processes governing the transport of chelated radionuclides in groundwater remain conceptually unclear for highly structured, heterogeneous environments. The objectives of this research were to provide an improved understanding and predictive capability of the hydrological and geochemical mechanisms that control the transport behavior of chelated radionuclides and metals in anoxic subsurface environments that are complicated by fracture flow and matrix diffusion. Our approach involved a long-term, steady-state natural gradient field experiment where nonreactive Br- and reactive 57Co(II)EDTA2- 109CdEDTA2-, and 51Cr(VI) were injected into a fracture zone of a contaminated fractured shale bedrock. The spatial and temporal distribution of the tracer and solutes was monitored for 500 days using an array of groundwater sampling wells instrumented within the fast-flowing fracture regime and a slower flowing matrix regime. The tracers were preferentially transported along strike-parallel fractures coupled with the slow diffusion of significant tracer mass into the bedrock matrix. The chelated radionuclides and metals were significantly retarded by the solid phase with the mechanisms of retardation largely due to redox reactions and sorption coupled with mineral-induced chelate-radionuclide dissociation. The formation of significant Fe(III)EDTA byproduct that accompanied the dissociation of the radionuclide-chelate complexes was believed to be the result of surface interactions with biotite which was the only Fe(III)-bearing mineral phase present in these Fe-reducing environments. These results counter current conceptual models that suggest chelated contaminants move conservatively through Fe-reducing environments since they are devoid of Fe-oxyhydroxides that are known to aggressively compete for chelates in oxic regimes. Modeling results further demonstrated that chelate-radionuclide dissociation reactions were most prevalent along fractures where accelerated weathering processes are expected to expose more primary minerals than the surrounding rock matrix. The findings of this study suggest that physical retardation mechanisms (i.e. diffusion) are dominant within the matrix regime, whereas geochemical retardation mechanisms are dominant within the fracture regime.     

Laboratory migration experiments with radionuclides and natural colloids in a granite fracture

Natural colloids in groundwater could facilitate radionuclide transport, provided the colloids are mobile, are present in sufficient concentrations and can adsorb radionuclides. This paper describes the results of a laboratory migration study carried out with combinations of radionuclides and natural colloids within a fracture in a large granite block to experimentally determine the impact of colloids on radionuclide transport. The 85Sr used in this study is an example of a moderately sorbing radionuclide, while the 241Am is typical of a strongly sorbed radionuclide with very low solubility. The natural colloids used in this study were isolated from granite groundwater from Atomic Energy of Canada (AECL) Underground Research Laboratory (URL), and consisted of mostly 1-10 nm organic colloids, along with lesser amounts of 10-450 nm colloids (organics and aluminosilicates). The measured coefficients for radionuclide sorption onto these colloids were between 3 x 10(2) and 1 x 10(3) ml/g for 85Sr, and between 7 x 10(4) and 7 x 10(5) mg/l for 241Am. The 85Sr sorption on the natural colloids appeared to be reversible. Migration experiments in the granite block were carried out by establishing a flow field between two boreholes (out of a total of nine) intersecting a main horizontal fracture. These experiments showed that dissolved 85Sr behaved as a moderately sorbing tracer, while dissolved 241Am was completely adsorbed by the fracture surfaces and showed no evidence of transport. However, when natural colloids were injected together with dissolved 241Am, a small amount of 241Am transport was observed, demonstrating the ability of natural colloids to facilitate the transport of radionuclides with low solubility. Natural colloids had only a minor effect on the transport of 85Sr. In a separate experiment to test the effect of higher colloid concentrations on 85Sr migration, synthetic colloids were produced from Avonlea bentonite. The introduction of a relatively high concentration of bentonite colloids actually reduced 85Sr transport because, compared to natural colloids, the bentonite colloids were less mobile and they sorbed 85Sr more strongly.     

A triple-continuum approach for modeling flow and transport processes in fractured rock

This paper presents a triple-continuum conceptual model for simulating flow and transport processes in fractured rock. Field data collected from the unsaturated zone of Yucca Mountain, a repository site of high-level nuclear waste, show a large number of small-scale fractures. The effect of these small fractures has not been considered in previous modeling investigations within the context of a continuum approach. A new triple-continuum model (consisting of matrix, small-fracture, and large-fracture continua) has been developed to investigate the effect of these small fractures. This paper derives the model formulation and discusses the basic triple-continuum behavior of flow and transport processes under different conditions, using both analytical solutions and numerical approaches. The simulation results from the site-scale model of the unsaturated zone of Yucca Mountain indicate that these small fractures may have an important effect on radionuclide transport within the mountain.     

Flow and transport in unsaturated fractured rock: effects of multiscale heterogeneity of hydrogeologic properties

The heterogeneity of hydrogeologic properties at different scales may have different effects on flow and transport processes in a subsurface system. A model for the unsaturated zone of Yucca Mountain, Nevada, is developed to represent complex heterogeneity at two different scales: (1) layer scale corresponding to geologic layering and (2) local scale. The layer-scale hydrogeologic properties are obtained using inverse modeling, based on the available measurements collected from the Yucca Mountain site. Calibration results show a significant lateral and vertical variability in matrix and fracture properties. Hydrogeologic property distributions in a two-dimensional, vertical cross-section of the site are generated by combining the average layer-scale matrix and fracture properties with local-scale perturbations generated using a stochastic simulation method. The unsaturated water flow and conservative (nonsorbing) tracer transport through the cross-section are simulated for different sets of matrix and fracture property fields. Comparison of simulation results indicates that the local-scale heterogeneity of matrix and fracture properties has a considerable effect on unsaturated flow processes, leading to fast flow paths in fractures and the matrix. These paths shorten the travel time of a conservative tracer from the source (repository) horizon in the unsaturated zone to the water table for small fractions of total released tracer mass. As a result, the local-scale heterogeneity also has a noticeable effect on global tracer transport processes, characterized by an average breakthrough curve at the water table, especially at the early arrival time of tracer mass. However, the effect is not significant at the later time after 20% tracer mass reaches the water table. The simulation results also verify that matrix diffusion plays an important role in overall solute transport processes in the unsaturated zone at Yucca Mountain.     

Identification and Characterization of Potential Discharge Areas for Radionuclide Transport by Groundwater from a Nuclear Waste Repository in Sweden

This paper describes solute transport modeling carried out as a part of an assessment of the long-term radiological safety of a planned deep rock repository for spent nuclear fuel in Forsmark, Sweden. Specifically, it presents transport modeling performed to locate and describe discharge areas for groundwater potentially carrying radionuclides from the repository to the surface where man and the environment could be affected by the contamination. The modeling results show that topography to large extent determines the discharge locations. Present and future lake and wetland objects are central for the radionuclide transport and dose calculations in the safety assessment. Results of detailed transport modeling focusing on the regolith and the upper part of the rock indicate that the identification of discharge areas and objects considered in the safety assessment is robust in the sense that it does not change when a more detailed model representation is used.     

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.     

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.     

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.     

Variably saturated flow and multicomponent biogeochemical reactive transport modeling of a uranium bioremediation field experiment

Three-dimensional, coupled variably saturated flow and biogeochemical reactive transport modeling of a 2008 in situ uranium bioremediation field experiment is used to better understand the interplay of transport and biogeochemical reactions controlling uranium behavior under pulsed acetate amendment, seasonal water table variation, spatially variable physical (hydraulic conductivity, porosity) and geochemical (reactive surface area) material properties. While the simulation of the 2008 Big Rusty acetate biostimulation field experiment in Rifle, Colorado was generally consistent with behaviors identified in previous field experiments at the Rifle IFRC site, the additional process and property detail provided several new insights. A principal conclusion from this work is that uranium bioreduction is most effective when acetate, in excess of the sulfate-reducing bacteria demand, is available to the metal-reducing bacteria. The inclusion of an initially small population of slow growing sulfate-reducing bacteria identified in proteomic analyses led to an additional source of Fe(II) from the dissolution of Fe(III) minerals promoted by biogenic sulfide. The falling water table during the experiment significantly reduced the saturated thickness of the aquifer and resulted in reactants and products, as well as unmitigated uranium, in the newly unsaturated vadose zone. High permeability sandy gravel structures resulted in locally high flow rates in the vicinity of injection wells that increased acetate dilution. In downgradient locations, these structures created preferential flow paths for acetate delivery that enhanced local zones of TEAP reactivity and subsidiary reactions. Conversely, smaller transport rates associated with the lower permeability lithofacies (e.g., fine) and vadose zone were shown to limit acetate access and reaction. Once accessed by acetate, however, these same zones limited subsequent acetate dilution and provided longer residence times that resulted in higher concentrations of TEAP reaction products when terminal electron donors and acceptors were not limiting. Finally, facies-based porosity and reactive surface area variations were shown to affect aqueous uranium concentration distributions with localized effects of the fine lithofacies having the largest impact on U(VI) surface complexation. The ability to model the comprehensive biogeochemical reaction network, and spatially and temporally variable processes, properties, and conditions controlling uranium behavior during engineered bioremediation in the naturally complex Rifle IFRC subsurface system required a subsurface simulator that could use the large memory and computational performance of a massively parallel computer. In this case, the eSTOMP simulator, operating on 128 processor cores for 12h, was used to simulate the 110-day field experiment and 50 days of post-biostimulation behavior.     

Transport of Sr2+ and SrEDTA2- in partially-saturated and heterogeneous sediments

Strontium-90 has migrated deep into the unsaturated subsurface beneath leaking storage tanks in the Waste Management Areas (WMA) at the U.S. Department of Energy's (DOE) Hanford Reservation. Faster than expected transport of contaminants in the vadose zone is typically attributed to either physical hydrologic processes such as development of preferential flow pathways, or to geochemical processes such as the formation of stable, anionic complexes with organic chelates, e.g., ethylenediaminetetraacetic acid (EDTA). The goal of this paper is to determine whether hydrological processes in the Hanford sediments can influence the geochemistry of the system and hence control transport of Sr(2+) and SrEDTA(2-). The study used batch isotherms, saturated packed column experiments, and an unsaturated transport experiment in an undisturbed core. Isotherms and repacked column experiments suggested that the SrEDTA(2-) complex was unstable in the presence of Hanford sediments, resulting in dissociation and transport of Sr(2+) as a divalent cation. A decrease in sorption with increasing solid:solution ratio for Sr(2+) and SrEDTA(2-) suggested mineral dissolution resulted in competition for sorption sites and the formation of stable aqueous complexes. This was confirmed by detection of MgEDTA(2-), MnEDTA(2-), PbEDTA(2-), and unidentified Sr and Ca complexes. Displacement of Sr(2+) through a partially-saturated undisturbed core resulted in less retardation and more irreversible sorption than was observed in the saturated repacked columns, and model results suggested a significant reservoir (49%) of immobile water was present during transport through the heterogeneous layered sediments. The undisturbed core was subsequently disassembled along distinct bedding planes and subjected to sequential extractions. Strontium was unequally distributed between carbonates (49%), ion exchange sites (37%), and the oxide (14%) fraction. An inverse relationship between mass wetness and Sr suggested that sandy sediments of low water content constituted the immobile flow regime. Our results suggested that the sequestration of Sr(2+) in partially-saturated, heterogeneous sediments was most likely due to the formation of immobile water in drier regions having low hydraulic conductivities.     

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.     

Modeling colloid transport for performance assessment

The natural system is expected to contribute to isolation at the proposed high-level nuclear waste (HLW) geologic repository at Yucca Mountain, NV (YM). In developing performance assessment (PA) computer models to simulate long-term behavior at YM, colloidal transport of radionuclides has been proposed as a critical factor because of the possible reduced interaction with the geologic media. Site-specific information on the chemistry and natural colloid concentration of saturated zone groundwaters in the vicinity of YM is combined with a surface complexation sorption model to evaluate the impact of natural colloids on calculated retardation factors (RF) for several radioelements of concern in PA. Inclusion of colloids into the conceptual model can reduce the calculated effective retardation significantly. Strongly sorbed radionuclides such as americium and thorium are most affected by pseudocolloid formation and transport, with a potential reduction in RF of several orders of magnitude. Radioelements that are less strongly sorbed under YM conditions, such as uranium and neptunium, are not affected significantly by colloid transport, and transport of plutonium in the valence state is only moderately enhanced. Model results showed no increase in the peak mean annual total effective dose equivalent (TEDE) within a compliance period of 10,000 years, although this is strongly dependent on container life in the base case scenario. At longer times, simulated container failures increase and the TEDE from the colloidal models increased by a factor of 60 from the base case. By using mechanistic models and sensitivity analyses to determine what parameters and transport processes affect the TEDE, colloidal transport in future versions of the TPA code can be represented more accurately.