Humic colloid-borne migration of uranium in sand columns

Column experiments were carried out to investigate the influence of humic colloids on subsurface uranium migration. The columns were packed with well-characterized aeolian quartz sand and equilibrated with groundwater rich in humic colloids (dissolved organic carbon (DOC): 30 mg dm(-3)). U migration was studied under an Ar/1% CO2 gas atmosphere as a function of the migration time, which was controlled by the flow velocity or the column length. In addition, the contact time of U with groundwater prior to introduction into a column was varied. U(VI) was found to be the dominant oxidation state in the spiked groundwater. The breakthrough curves indicate that U was transported as a humic colloid-borne species with a velocity up to 5% faster than the mean groundwater flow. The fraction of humic colloid-borne species increases with increasing prior contact time and also with decreasing migration time. The migration behavior was attributed to a kinetically controlled association/dissociation of U onto and from humic colloids and also a subsequent sorption of U onto the sediment surface. The column experiments provide an insight into humic colloid-mediated U migration in subsurface aquifers.     

Migration of uranium(IV)/(VI) in the presence of humic acids in quartz sand: a laboratory column study

The migration behavior of U(IV) and U(VI) in the presence of humic acid was studied in a quartz sand system. Laboratory column experiments were performed using humic acid, U(VI) in humic acid absence, U(IV) and U(VI) in humic acid presence, and for comparison a conservative tracer. In experiments using humic acid, both redox species of U migrate nearly as fast as the conservative tracer. Humic acid accelerates the U(VI) breakthrough compared to the humic acid-free system. There are strong indications for a similar effect on the U(IV) transport. At the same time, a part of U(IV) and U(VI) associated with the humic acid is immobilized in the quartz sand due to humic colloid filtration thus producing a delaying effect. Tailing at a low concentration level was observed upon tracer elution. The experimental breakthrough curves were described by reactive transport modeling using equations for equilibrium and kinetic reactions. The present study demonstrates that humic acids can play an important role in the migration of actinides. As natural organic matter is ubiquitous in aquifer systems, the humic colloid-borne transport of actinides is of high relevance in performance assessment.     

Development of a scalable model for predicting arsenic transport coupled with oxidation and adsorption reactions

Understanding the fundamentals of arsenic adsorption and oxidation reactions is critical for predicting its transport dynamics in groundwater systems. We completed batch experiments to study the interactions of arsenic with a common MnO2(s) mineral, pyrolusite. The reaction kinetics and adsorption isotherm developed from the batch experiments were integrated into a scalable reactive transport model to facilitate column-scale transport predictions. We then completed a set of column experiments to test the predictive capability of the reactive transport model. Our batch results indicated that the commonly used pseudo-first order kinetics for As(III) oxidation reaction neglects the scaling effects with respect to the MnO2(s) concentration. A second order kinetic equation that explicitly includes MnO2(s) concentration dependence is a more appropriate kinetic model to describe arsenic oxidation by MnO2(s) minerals. The arsenic adsorption reaction follows the Langmuir isotherm with the adsorption capacity of 0.053micromol of As(V)/g of MnO2(s) at the tested conditions. The knowledge gained from the batch experiments was used to develop a conceptual model for describing arsenic reactive transport at a column scale. The proposed conceptual model was integrated within a reactive transport code that accurately predicted the breakthrough profiles observed in multiple column experiments. The kinetic and adsorption process details obtained from the batch experiments were valuable data for scaling to predict the column-scale reactive transport of arsenic in MnO2(s)-containing sand columns.     

Actinide geochemistry: From the molecular level to the real system

Geochemical processes leading to either mobilization or retention of radionuclides in an aquifer system are significantly influenced by their interaction with rock, sediment and colloid surfaces. Therefore, a sound safety assessment of nuclear waste disposal requires the elucidation and quantification of those processes. State-of-the-art analytical techniques as e.g. laser- and X-ray spectroscopy are increasingly applied to study solid–liquid interface reactions to obtain molecular level speciation insight.We have studied the sorption of trivalent lanthanides and actinides onto aluminium oxides, hydroxides and purified clay minerals by the time-resolved laser fluorescence spectroscopy and X-ray-absorption spectroscopy. Chemical constitution and structure of surface bound actinides are proposed based on spectroscopic information. Open questions still remain with regard to the exact nature of mineral surface ligands and the mineral/water interface. Similarities of spectroscopic data obtained for M(III) sorbed onto γ-alumina, and clay minerals suggest the formation of very comparable inner-sphere surface complexes such as S–O–An(III)(OH)_x^(2 − x)(H₂O)_5 − x at pH > 5. Those speciation data are found consistent with those predicted by surface complexation modelling. The applicability of data obtained for pure mineral phases to actinide sorption onto heterogeneously composed natural clay rock is examined by experiments and by geochemical modelling. Good agreement of experiment and model calculations is found for U(VI) and trivalent actinide/lanthanide sorption to natural clay rock. The agreement of spectroscopy, geochemical modelling and batch experiments with natural rock samples and purified minerals increases the reliability in model predictions.The assessment of colloid borne actinide migration observed in various laboratory and field studies calls for detailed information on actinide–colloid interaction. Kinetic stabilization of colloid bound actinides can be due to inclusion into inorganic colloid matrix or by macromolecular rearrangement in case of organic, humic/fulvic like colloids. Only a combination of spectroscopy, microscopy and classical batch sorption experiments can help to elucidate the actinide–colloid interaction mechanisms and thus contribute to the assessment of colloids for radionuclide migration.     

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.     

The lanthanum precipitation method. Part 2: quantification of the conditional interaction constant between technetium(IV) and humic substances

The interaction between colloidal Tc(IV) species and colloidal Gorleben humic substances (HS) was quantified after application of the La-precipitation method on supernatant solutions obtained under various experimental conditions but at constant ionic strength of the Gorleben groundwater (0.04M). The determined interaction constant LogK(HS) (2.3+/-0.3) remained unchanged over a large range of Tc(IV) and HS concentrations and was independent of the pH of the original supernatant solution (pH range 6-10), Tc(IV)-HS loading (10(-3)-10(-6)molTcg(-1) HS) and the nature of the reducing surface (Magnetite, Pyrite and Gorleben sand) used for the pertechnetate reduction. The LogK(HS) value determined by the La-precipitation method is lower than the LogK value obtained from a previous study where the interaction between colloidal Tc(IV) species and Gorleben humic substances was quantified using a modified Schubert approach (2.6+/-0.3). The La-precipitation method allows to accurately determine the amount of Tc(IV) associated with HS but leads to a (small) overestimation of the free inorganic Tc(IV) species.     

Transport of reactive colloids and contaminants in groundwater: effect of nonlinear kinetic interactions

Transport of reactive colloids in groundwater may enhance the transport of contaminants in groundwater. Often, the interpretation of results of transport experiments is not a simple task as both reactions of colloids with the solid matrix and reactions of contaminants with the solid matrix and mobile and immobile colloids may be time dependent and nonlinear. Further colloid transport properties may differ from solute transport properties. In this paper, a one-dimensional model for coupled and contaminant in a porous medium (COLTRAP) is presented together with simulation results. Calculated breakthrough curves (BTC's) during contamination and decontamination show systematically the effect of nonlinear and kinetic interactions on contaminant transport in the presence of reactive colloids, and the effect of colloid transport properties that differ from solute transport properties. It is shown that in case of linear kinetic reactions, the rate of exchange of mobile and immobile colloids have a large impact on the shape of BTC's even if the solid matrix is saturated with respect to colloids. BTC's during the contamination and decontamination phase have identical shapes in this case. Moreover, the slow reactions of contaminants and colloids may lead to unretarded breakthrough of contaminants. Independent of reaction rates, nonlinear reactions lead to BTC's that are steeper during contamination than in the linear case. A characteristic aspect of nonlinear sorption is that shapes of BTC's differ during the contamination and decontamination phase. It has been observed that shapes of some of the simulated adsorption and desorption curves are similar as shapes found in experiments reported in literature. This stresses the importance of incorporating both kinetics and nonlinearity in models for coupled colloid and contaminant transport and the capability of COLTRAP to interpret experimental results. Finally, to figure out whether nonlinear processes play a role, it is very important to consider both contamination and decontamination in transport experiments.     

Influence of trivalent electrolytes on the humic colloid-borne transport of contaminant metals: competition and flocculation effects

With the objective to assess the relevance of competitive effects in respect of the humic colloid-borne migration of actinides in case of release, the influence of Al(III) on humate complexation of La(III) as an analogue of trivalent actinides was investigated for various humic materials by using 140La as a radioactive tracer, allowing measurements in very dilute systems to simulate realistic settings. Generally, competition by aluminium is not detectable unless the metal-loading capacity of the humic colloids is nearly exhausted. For average contents of organic carbon, a threshold value of 10(-6) M Al(III) can be specified. The metal exchange turned out to be kinetically hindered. Effects on co-adsorption of La(III) and humic acid were found to be less important. Immobilization by the concomitantly induced flocculation process outweighs the role of displacement effects. Comparative studies on complexation and flocculation of humic acid with Al(III), Ga(III), In(III), Sc(III), Y(III), and La(III) were undertaken in order to evaluate the influence of specific properties apart from ion charge and to characterize the mechanism of flocculation. In spite of considerable variations in the binding affinities among these metals, it can be inferred that the possibility of significant competitive effects in natural aquatic systems is confined to Al(III). Complex stabilities and flocculation efficiencies proved to be interrelated. Precipitation is thus attributed to homocoagulation of humic colloids induced by charge compensation, which is further supported by flocculation experiments with Al(III) depending on pH, ionic strength, and humic acid concentration.     

Migration of colloids in discretely fractured porous media: effect of colloidal matrix diffusion

Numerical simulations of colloid transport in discretely fractured porous media were performed to investigate the importance of matrix diffusion of colloids as well as the filtration and remobilization of colloidal particles in both the fractures and porous matrix. To achieve this objective a finite element numerical code entitled COLDIFF was developed. The processes that COLDIFF takes into account include advective-dispersive transport of colloids, filtration and remobilization of colloidal particles in both fractures and porous matrix, and diffusive interactions of colloids between the fractures and porous matrix. Three sets of simulations were conducted to examine the importance of parameters and processes controlling colloid migration. First, a sensitivity analysis was performed using a porous block containing a single fracture to determine the relative importance of various phenomenological coefficients on colloid transport. The primary result of the analysis showed that the porosity of the matrix and the process of colloid filtration in fractures play important roles in controlling colloid migration. Second, simulations were performed to replicate and examine the results of a laboratory column study using a fractured shale saprolite. Results of this analysis showed that the filtration of colloidal particles in the porous matrix can greatly affect the tailing of colloid concentrations after the colloid source was removed. Finally, field-scale simulations were performed to examine the effect of matrix porosity, fracture filtration and fracture remobilization on long-term colloid concentration and migration distance. The field scale simulations indicated that matrix diffusion and fracture filtration can significantly reduce colloid migration distance. Results of all three analyses indicated that in environments where porosity is relatively high and colloidal particles are small enough to diffuse out of fractures, the characteristics of the porous matrix that affect colloid transport become more important than those of the fracture network. Because the properties of the fracture network tend to have greater uncertainty due to difficulties in their measurement relative to those of the porous matrix, prediction uncertainties associated with colloid transport in discretely fractured porous media may be reduced.     

Sensitivity analysis of radionuclide migration in compacted bentonite: a mechanistic model approach

Mechanistic model calculations for the migration of Cs, Ra, Am and Pb in compacted bentonite have been carried out to evaluate sensitivities with respect to different parameter variations. A surface chemical speciation/electric double layer model is used to calculate: (i) porewater composition and radionuclide speciation in solution and at the bentonite surface, yielding the distribution of mobile and sorbed species and (ii) interaction of diffusing species with negatively charged pore walls to obtain diffusion parameters. The basic scenario considers the interaction of compacted bentonite with a fresh-type groundwater; variations include the presence of bentonite impurities and saline groundwater. It is shown that these scenarios result in significant variations of porewater composition that affect migration via three mechanisms that can partly compensate each other: (1) effects on sorption through radionuclide complexation in solution, and competition of major cations for surface sites; (2) changes in radionuclide solution speciation leading to different diffusing species under different conditions; (3) effects on diffusion through changes in the electric double layer properties of the clay pores as a function of ionic strength.     

Chromium (III) recovery from waste acid solution by ion exchange processing using Amberlite IR-120 resin: batch and continuous ion exchange modelling

The use of ion exchange technology was studied to remove chromium (III) from acidic waste solution by Amberlite IR-120 resin. Batch and column experimental tests were conducted to provide data for theoretical models and verify the system performance of the adsorption process. Results of batch equilibrium tests indicated that Langmuir isotherm describes well the adsorption process, whereas experimental data also provide evidence that, under the present experimental conditions, chromium (III) adsorption by Amberlite IR-120 resin is film-diffusion controlled; on the other hand, the theoretical model used in the present investigation was found to predict reasonably well the ion exchange breakthrough performance.     

The lanthanum precipitation method. Part 1: a new method for technetium(IV) speciation in humic rich natural groundwater

A new and quick method for direct speciation of Tc(IV) in humic rich solutions, based on the induced aggregation of humic substances in the presence of the trivalent cation La(3+), is presented. This method (the "La-precipitation method") allows flocculating all the humic substances and also the Tc(IV) associated with humic substances. The method is tested on solutions containing Tc(IV) and Gorleben humic substances. The influence of different parameters (humic substance concentration, Tc concentration, reaction time and pH) is investigated on the observed free Tc(IV) concentration after precipitation of all humic substances. None of these parameters had a (significant) influence on the observed Tc(IV) concentration in solution after addition of La(3+) to Tc(IV)-HS containing solutions. It is therefore proposed that the method can be used to separate the Tc(IV) bound to humic substances from the free inorganic Tc species in solution.     

Soil column experiments used as a means to assess transport, sorption, and biodegradation of pesticides in groundwater

Soil column experiments are used to investigate the fate of three pesticides of high, intermediate, and low solubility in groundwater: N- phosphonomethyl glycine (glyphosate); O,O-diethyl-S-[(ethylthio)methyl]phosphorodithioate (phorate); (2,4-dichlorophenoxy)acetic acid (2,4-D). Feed solutions are prepared by adding each pesticide (100 mg/L glyphosate, 50 micro g/L phorate, 50 mg/L 2,4-D) along with conservative tracer, KBr, in synthetic groundwater. The concentration of the pesticides in effluents is detected by ion chromatography (glyphosate, 2,4-D) and GC-FID (phorate). The Br(-) breakthrough curves are employed to estimate the dispersion coefficient and mean pore velocity in each column. Solute transport and reactive models accounting for equilibrium/non-equilibrium sorption and biodegradation are coupled with inverse modeling numerical codes to estimate the kinetic parameters for all pesticides.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

Selectivity in the complexation of actinides by humic substances

The interactions of a range of actinide elements (Th, U, Np, Pu, Am) with humic substances from the Needle's Eye natural analogue site were studied by gel permeation chromatography. Bulk humic substances were isolated by ammonia extraction, followed by dialysis against distilled water and freeze-drying. The gel permeation results suggest that Needle's Eye humic substances can be fractionated into three incompletely resolved fractions with average molecular weights determined by analytical ultracentrifugation around 49 000 for Fraction 1, around 14 700 for Fraction 2 and around 8000 for Fraction 3. Although there are significant differences between the organic matter elution patterns in individual gel permeation experiments, presumably due to differences in column packing, these are much smaller than the differences between metal ions. The uranium that is naturally present in these humic substances is largely bound in the late-eluting fraction. Spikes of the early actinides, including Np and Pu in controlled valency states, have been added to the humic substances, and gel permeation of the spiked humic substances shows that the three humic fractions vary greatly in their effectiveness and selectivity as ligands for early actinides.     

Soil column experiments used as a means to assess transport,sorption, and biodegradation of pesticides in groundwater

Soil column experiments are used to investigate the fate of three pesticides of high, intermediate, and low solubility in groundwater: N- phosphonomethyl glycine (glyphosate); O,O-diethyl-S-[(ethylthio)methyl]phosphorodithioate (phorate); (2,4-dichlorophenoxy)acetic acid (2,4-D). Feed solutions are prepared by adding each pesticide (100 mg/L glyphosate, 50 μ g/L phorate, 50 mg/L 2,4-D) along with conservative tracer, KBr, in synthetic groundwater. The concentration of the pesticides in effluents is detected by ion chromatography (glyphosate, 2,4-D) and GC-FID (phorate). The Br^? breakthrough curves are employed to estimate the dispersion coefficient and mean pore velocity in each column. Solute transport and reactive models accounting for equilibrium/non-equilibrium sorption and biodegradation are coupled with inverse modeling numerical codes to estimate the kinetic parameters for all pesticides.     

Groundwater flow, multicomponent transport and biogeochemistry: development and application of a coupled process model

A research tool for modeling the reactive flow and transport of groundwater contaminants in multiple dimensions is presented. Arbitrarily complex coupled kinetic–equilibrium heterogeneous reaction networks, automatic code generation, transfer-function based solutions, parameter estimation, high-resolution methods for advection, and robust solvers for the mixed kinetic–equilibrium chemistry are some of the features of reactive flow and transport (RAFT) that make it a versatile research tool in the modeling of a wide variety of laboratory and field experiments. The treatment of reactions is quite general so that RAFT can be used to model biological, adsorption/desorption, complexation, and mineral dissolution/precipitation reactions among others. The integrated framework involving automated code generation and parameter estimation allows for the development, characterization, and evaluation of mechanistic process models. The model is described and used to solve a problem in competitive adsorption that illustrates some of these features. The model is also used to study the development of an in situ Fe(II)-zone by encouraging the growth of an iron-reducing bacterium with lactate as the electron donor. Such redox barriers are effective in sequestering groundwater contaminants such as chromate and TCE.     

Database uncertainty as a limiting factor in reactive transport prognosis

The effect of uncertainties in thermodynamic databases on prediction performances of reactive transport modeling of uranium (VI) is investigated with a Monte Carlo approach using the transport code TReaC. TReaC couples the transport model to the speciation code PHREEQC by a particle tracking method. A speciation example is given to illustrate the effect of uncertainty in thermodynamic data on the predicted solution composition. The transport calculations consequently show the prediction uncertainty resulting from uncertainty in thermodynamic data. A conceptually simple scenario of elution of uranium from a sand column is used as an illustrating example. Two different cases are investigated: a carbonate-enriched drinking water and an acid mine water associated with uranium mine remediation problems. Due to the uncertainty in the relative amount of positively charged and neutral solution species, the uncertainty in the thermodynamic data also infers uncertainty in the retardation behavior. The carbonated water system shows the largest uncertainties in speciation calculation. Therefore, the model predictions of total uranium solubility have a broad range. The effect of data uncertainty in transport prediction is further illustrated by a prediction of the time when eluted uranium from the column exceeds a threshold value. All of these Monte Carlo transport calculations consume large amounts of computing time.     

Biogeochemical influence on transport of chromium in manganese sediments: experimental and modeling approaches

Hexavalent chromium (Cr(VI)) was reduced to immobile and nontoxic Cr(III) by a dissimilatory metal reducing bacteria, Shewanella alga Simidu (BrY-MT) ATCC 55627. A series of kinetic batch and dynamic column experiments were conducted to provide an understanding of Cr(VI) reduction by the facultative anaerobe BrY-MT. Reduction of Cr(VI) was rapid (within 1 h) in columns packed with quartz sand and bacteria, whereas Cr(VI) reduction by BrY-MT was delayed (57 h) in the presence of beta-MnO2-coated sand. A mathematical model was developed and evaluated against data obtained from column experiments. The model takes into account (1) advective-dispersive transport of Cr(III), Cr(VI), lactate, and protein (mobile and immobile bacteria); (2) first-order kinetic adsorption of Cr(III) and lactate; (3) conversion of solid phase beta-MnO2 to solid phase MnOOH due to oxidation of Cr(III); (4) dual-Monod kinetics, where Cr(VI) is the electron acceptor and lactate is the electron donor. The breakthrough data for Cr(III), Cr(VI), lactate, and protein (mobile and immobile bacteria) were fitted simultaneously. The breakthrough data are well described by the mathematical model that considers the above processes. This result demonstrates the ability of the coupled hydrobiogeochemical model to simulate chromium transport in complex reactive systems.     

Contaminant transport in groundwater in the presence of colloids and bacteria: Model development and verification

Colloids and bacteria (microorganisms) naturally exist in groundwater aquifers and can significantly impact contaminant migration rates. A conceptual model is first developed to account for the different physiochemical and biological processes, reaction kinetics, and different transport mechanisms of the combined system (contaminant–colloids–bacteria). All three constituents are assumed to be reactive with the reactions taking place between each constituent and the porous medium and also among the different constituents. A general linear kinetic reaction model is assumed for all reactive processes considered. The mathematical model is represented by fourteen coupled partial differential equations describing mass balance and reaction processes. Two of these equations describe colloid movement and reactions with the porous medium, four equations describe bacterial movement and reactions with colloids and the porous medium, and the remaining eight equations describe contaminant movement and its reactions with bacteria, colloids, and the porous medium. The mass balance equations are numerically solved for two-dimensional groundwater systems using a third-order, total variance-diminishing scheme (TVD) for the advection terms. Due to the complex coupling of the equations, they are solved iteratively each time step until a convergence criterion is met. The model is tested against experimental data and the results are favorable.