Understanding the parameters that control colloid-mediated transport of radionuclides is important for the safe disposal of used nuclear fuel. We report an experimental and reactive transport ...modeling examination of americium transport in a groundwater–bentonite–fracture fill material system. A series of batch sorption and column transport experiments were conducted to determine the role of desorption kinetics from bentonite colloids in the transport of americium through fracture materials. We used fracture fill material from a shear zone in altered granodiorite collected from the Grimsel Test Site (GTS) in Switzerland and colloidal suspensions generated from FEBEX bentonite, a potential repository backfill material. The colloidal suspension (100 mg L−1) was prepared in synthetic groundwater that matched the natural water chemistry at GTS and was spiked with 5.5 × 10−10 M 241Am. Batch characterizations indicated that 97% of the americium in the stock suspension was adsorbed to the colloids. Breakthrough experiments conducted by injecting the americium colloidal suspension through three identical columns in series, each with mean residence times of 6 h, show that more than 95% of the bentonite colloids were transported through each of the columns, with modeled colloid filtration rates (kf) of 0.01–0.02 h−1. Am recoveries in each column were 55–60%, and Am desorption rate constants from the colloids, determined from 1-D transport modeling, were 0.96, 0.98, and 0.91 h−1 in the three columns, respectively. The consistency in Am recoveries and desorption rate constants in each column indicates that the Am was not associated with binding sites of widely-varying strengths on the colloids, as one binding site with fast kinetics represented the system accurately for all three sequential columns. Our data suggest that colloid-mediated transport of Am in a bentonite-fracture fill material system is unlikely to result in transport over long distance scales because of the ability of the fracture materials to rapidly strip Am from the bentonite colloids and the apparent lack of a strong binding site that would keep a fraction of the Am strongly-associated with the colloids.
Display omitted
•97% of 241Am (5.5 × 10−10 M) in a bentonite suspension was associated with bentonite colloids.•95–99% of bentonite colloids were transported through the fracture fill material columns.•Am desorption from bentonite was modeled/described with a single type of binding site.•Modeled Am desorption rate coefficients from bentonite colloids were 0.091–0.098 h−1.•Significant Am transport in an altered granite fracture with oxidizing pH 8 water is unlikely.
Understanding colloid transport in ground water is essential to assessing the migration of colloid-size contaminants, the facilitation of dissolved contaminant transport by colloids, in situ ...bioremediation, and the health risks of pathogen contamination in drinking water wells. Much has been learned through laboratory and field-scale colloid tracer tests, but progress has been hampered by a lack of consistent tracer testing methodology at different scales and fluid velocities. This paper presents laboratory and field tracer tests in fractured rock that use the same type of colloid tracer over an almost three orders-of-magnitude range in scale and fluid velocity. Fluorescently-dyed carboxylate-modified latex (CML) microspheres (0.19 to 0.98 micrometer diameter) were used as tracers in (1) a naturally fractured tuff sample, (2) a large block of naturally fractured granite, (3) a fractured granite field site, and (4) another fractured granite/schist field site. In all cases, the mean transport time of the microspheres was shorter than the solutes, regardless of detection limit. In all but the smallest scale test, only a fraction of the injected microsphere mass was recovered, with the smaller microspheres being recovered to a greater extent than the larger microspheres. Using existing theory, we hypothesise that the ohserved microsphere early arrival was due to volume exclusion and attenuation was due to aggregation and/or settling during transport. In most tests, microspheres were detected using flow cytometry, which proved to be an excellent method of analysis. CML microspheres appear to be useful tracers for fractured rock in forced gradient and short-term natural gradient tests, but longer residence times may result in small microsphere recoveries.
We use uranium (U) isotope ratios to detect and quantify the extent of natural U reduction in groundwater across a roll front redox gradient. Our study was conducted at the Smith Ranch-Highland in ...situ recovery (ISR) U mine in eastern Wyoming, USA, where economic U deposits occur in the Paleocene Fort Union formation. To evaluate the fate of aqueous U in and adjacent to the ore body, we investigated the chemical composition and isotope ratios of groundwater samples from the roll-front type ore body and surrounding monitoring wells of a previously mined area. The 238U/235U of groundwater varies by approximately 3‰ and is correlated with U concentrations. Fluid samples down-gradient of the ore zone are the most depleted in 238U and have the lowest U concentrations. Activity ratios of 234U/238U are ∼5.5 up-gradient of the ore zone, ∼1.0 in the ore zone, and between 2.3 and 3.7 in the down-gradient monitoring wells. High-precision measurements of 234U/238U and 238U/235U allow for development of a conceptual model that evaluates both the migration of U from the ore body and the extent of natural attenuation due to reduction. We find that the premining migration of U down-gradient of the delineated ore body is minimal along eight transects due to reduction in or adjacent to the ore body, whereas two other transects show little or no sign of reduction in the down-gradient region. These results suggest that characterization of U isotopic ratios at the mine planning stage, in conjunction with routine geochemical analyses, can be used to identify where more or less postmining remediation will be necessary.
Two cross-hole tracer tests involving the simultaneous injection of two nonsorbing solute tracers with different diffusion coefficients (bromide and pentafluorobenzoate) and one weakly sorbing solute ...tracer (lithium ion) were conducted in two different intervals at the C-wells complex near the site of a potential high-level nuclear waste repository at Yucca Mountain, NV. The tests were conducted to (1) test a conceptual radionuclide transport model for saturated, fractured tuffs near Yucca Mountain and (2) obtain transport parameter estimates for predictive modeling of radionuclide transport. The differences between the responses of the two nonsorbing tracers and the sorbing tracer (when normalized to injection masses) were consistent with a dual-porosity transport system in which matrix diffusion was occurring. The concentration attenuation of the sorbing tracer relative to the nonsorbing tracers suggested that diffusion occurred primarily into matrix pores, not simply into stagnant water within the fractures. The
K
d values deduced from the lithium responses were generally larger than
K
d values measured in laboratory batch sorption tests using crushed C-wells cores. This result supports the use of laboratory-derived
K
d values for predicting sorbing species transport at the site, as the laboratory
K
d values would result in underprediction of sorption and hence conservative transport predictions. The tracer tests also provided estimates of effective flow porosity and longitudinal dispersivity at the site. The tests clearly demonstrated the advantages of using multiple tracers of different physical and chemical characteristics to distinguish between alternative conceptual transport models and to obtain transport parameter estimates that are better constrained than can be obtained using only a single tracer or using multiple nonsorbing tracers without a sorbing tracer.
Characterizing the mobility of uranium and vanadium in groundwater with a hydraulic connection to surface water is important to inform the best management practices of former mill tailing sites. In ...this study, the recharge of river water to the unsaturated and saturated zones of a uranium-contaminated alluvial aquifer was simulated in a series of forced-gradient single- and multi-well injection-extraction tests. The injection fluid (river water) was traced with natural and artificial tracers that included halides, fluorobenzoates, lithium, and naphthalene sulfonate to characterize the potential mass transport mechanisms of uranium and vanadium. The extraction fluid (river water/groundwater mixture) was analyzed for the tracers, uranium, and vanadium. The results from the tracers indicated that matrix diffusion was likely negligible over the spatiotemporal scales of the tests as evident by nearly identical breakthrough curves of the halides and fluorobenzoates. In contrast, the breakthrough curves of lithium and naphthalene sulfonate indicated that sorption by cation exchange and sorption to organic matter, respectively, were potential mass transport mechanisms of uranium and vanadium. Uranium was mobilized in the saturated zone containing gypsum (gypsum-rich zone), the vadose zone (vadose-rich zone), and the saturated zone containing organic carbon (organic-rich zone) whereas vanadium was mobilized only in the saturated gypsum-rich zone. The mechanisms responsible for the mobilization of uranium and vanadium were likely dissolution of uranium- and vanadium-bearing minerals and/or desorption from the gypsum-rich zone, flushing of uranium from the vadose-rich zone, and desorption of uranium from the organic-rich zone due to the natural contrast in the geochemistry between the river water and groundwater. The experimental design of this study was unique in that it employed the use of multiple natural and artificial tracers coupled with a direct injection of native river water to groundwater. Here, these results demonstrated that natural recharge and flooding events at former mill tailing sites can mobilize uranium, and possibly vanadium, and contribute to persistent levels of groundwater contamination.
The objective of this work was to investigate flow and transport in a layered, variably saturated system consisting of both fractured rock and sedimentary material during focused infiltration from ...the surface. Two tracer tests were performed using the Vadose Zone Research Park (VZRP) at the Idaho National Laboratory (INL). The first test occurred under quasi-steady-state conditions and the second was initiated in a much drier system and thus provided information regarding flow and transport under transient conditions. A one-dimensional analytical model was used to fit breakthrough curves resulting from the two tracer tests. The results of this modeling provide insight into the nature of flow in the fractured basalt, surficial alluvium, and sedimentary interbeds that comprise the vadose zone of the eastern Snake River Plain. Flow through the fractured basalt is focused and preferential in nature, and multiple flow paths arise due to numerous fractures functioning as transmissive pathways in addition to flow splitting along geologic contacts. Flow velocities were significantly higher during the test with the wetter flow domain, presumably due to increases in hydraulic conductivity associated with higher water contents of the geologic materials. Perching was observed above the alluvium-basalt contact and above the lower boundary of a locally continuous sedimentary interbed. The perching behavior between the two contacts was fundamentally different; the perched layer above the alluvium-basalt contact was neither laterally extensive nor temporally persistent in the absence of infiltration from the surface. In contrast, the perched layer along the interbed was significantly thicker and gave rise to lateral flow over distances on the order of hundreds of meters. Vertical transport is shown to occur predominantly through the main bulk of the sedimentary material of the interbed; lateral flow appears to occur primarily in the fractured basalt directly above the interbed.
This paper presents a methodology for upscaling matrix‐material transport parameters in fractured‐flow dominated systems with multimodal reactive mineral facies. The upscaling method provides a ...theoretical and practical link between controlled experimental results at the laboratory/bench scale and multikilometer field scales at which contaminant remediation and risk assessment are actually conducted. As sorption reactions in matrix are in part determined by mineral properties, a new conceptual model is developed to reflect the hierarchical structure of reactive mineral facies at the microform, mesoform, and macroform scales. The conceptual model of hierarchical reactive matrix mineral facies is integrated with a dual‐porosity model for simulating diffusion of solutes out of fractures and sorption onto the matrix minerals. By assuming that sorption reactions primarily occur in the rock matrix, we develop a multimodal spatial random function for characterizing both the tortuosity (physical heterogeneity) and sorption coefficient (chemical heterogeneity) at different scales in the rock matrix. The effective tortuosity at the field scale is derived by volume averaging of mass transfer coefficient for a conservative species. Subsequently, using a sorbing species (e.g., uranium), we derive the equations for upscaling the sorption coefficients in a saturated, fractured‐rock system for field‐scale simulations. The effective field‐scale tortuosity and sorption coefficient are related to their mean, variance, integral scale, and domain size along a pathway through a three‐dimensional flow field. The upscaled values increase with the integral scale and are larger than their geometric mean. Simulations conducted with upscaled sorption coefficients and tortuousities are compared very well with high‐resolution Monte Carlo simulations capturing the parameter spatial variations. Results of this study can be extended to explore scale dependence of other important transport parameters for fractured‐rock solute transport.
A weakly sorbing cation, lithium, will be used as a reactive tracer in upcoming field tracer tests in the saturated alluvium south of Yucca Mountain, Nevada. One objective of the field tests is to ...determine how well field-scale reactive transport can be predicted using transport parameters derived from laboratory experiments. This paper describes several laboratory lithium batch sorption and column transport experiments that were conducted using ground water and alluvium obtained from the site of the planned field tests. In the batch experiments, isotherms were determined over 2.5 orders of magnitude of lithium concentrations, corresponding to the range expected in the field tests. In addition to measuring equilibrium lithium concentrations, concentrations of other cations, namely Na
+, K
+, and Ca
2+, were measured in the batch tests to determine Li
+-exchangeable equilibria. This information was used in conjunction with alluvium cation exchange capacity measurements to parameterize a three-component cation-exchange model (EQUIL) that describes lithium sorption in the alluvium system. This model was then applied to interpret the transport behavior of lithium ion in saturated alluvium column tests conducted at three different lithium bromide injection concentrations. The concentrations were selected such that lithium ion either dominated, accounted for a little over half, or accounted for only a small fraction of the total cation equivalents in the injection solution. Although tracer breakthrough curves differed significantly under each of these conditions, with highly asymmetric responses occurring at the highest injection concentrations, the three-component cation-exchange model reproduced the observed transport behavior of lithium and the other cations in each case with a similar set of model parameters. In contrast, a linear
K
d-type sorption model could only match the lithium responses at the lowest injection concentration. The three-component model will be used to interpret the field tests, with the expectation that it will help refine estimates of effective flow porosity, particularly if the lithium response curves are asymmetric.
An optical video microscopic system and image processing and data extraction and manipulation routines are developed for in situ detailed quantification of the deposition of colloids onto an ...arbitrary surface and determining their concentration distribution across the bulk suspension. The system produces a relatively large field of view (∼330×245μm) and utilizes dark-field light microscopy to visualize colloids as small as ∼0.3μm in diameter at the surface and in the bulk suspension with a sufficient resolution (∼0.5μm). On real-time basis, the routines automate various tasks from image capturing and processing to data manipulation, extraction, and display. The extracted data include: (i) surface concentration and flux of deposited, attached, and detached colloids, (ii) surface concentration and flux of arriving and departing colloids, (iii) distribution of colloids in the bulk suspension in the direction perpendicular to the deposition surface, and (iv) spatial and temporal distributions of deposited colloids. This article provides detailed description of the system and its image processing and data extraction and manipulation routines. Representative results of the deposition of 0.3-μm-diameter polystyrene colloids onto a glass surface, from a flowing suspension in a 0.02-cm-aperture parallel-plate channel, are presented and discussed.
The objective of this research was to investigate the effects of matrix diffusion on solute transport in fractured volcanic tuff. Two tuff cores were studied, one with a matrix porosity of 0.27 and ...the other with a porosity of 0.14. The matrix permeabilities of the cores were 4.7×10−15 and 7.8×10−19 m2, 5 and 9 orders of magnitude less than the respective fracture permeabilities. This suggested that the cores could be modeled as dual‐porosity systems with no flow in the matrix but significant solute storage capacity. Two types of tracer tests were conducted in each fractured core: (1) iodide was injected in separate experiments at different flow rates and (2) two tracers of different matrix diffusion coefficients (bromide and pentafluorobenzoate (PFBA)) were injected in another test. A difference in the maximum concentrations of the solutes and the extended tailing of the breakthrough curves were assumed to be indicative of diffusive mass transfer between the fracture and the porous matrix of the cores. Interpreting the results from both methods allowed the identification of matrix diffusion and dispersion effects within the fracture by simultaneously fitting the data sets (with known constraints) using a relatively simple conceptual model. Estimates of mass transfer coefficients for the fractured cores were also obtained.