Sorption of Eu(III)/Cm(III) on Ca-montmorillonite and Na-illite. Part 2: Surface complexation modelling

Sorption edges and isotherms for Eu(III) uptake on Ca-montmorillonite and Na-illite in 0.066 mol/L Ca(ClO₄)₂ and 0.1 mol/L NaClO₄ background electrolytes, respectively, were modelled using a quasi-mechanistic sorption model (the two site protolysis non electrostatic surface complexation and cation exchange (2SPNE SC/CE) model). For both clay minerals the Eu sorption edges could be quantitatively modelled in the pH range ∼3 to ∼10 using cation exchange reactions for Eu³⁺/Na⁺ and Eu³⁺/Ca²⁺ and three surface complexation reactions on the strong sorption sites forming ≡S^SOEu²⁺, ≡S^SOEuOH⁺ and ≡S^SOEu(OH)₂° inner sphere complexes which appear successively with increasing pH. Time resolved laser fluorescence spectroscopy (TRLFS) measurements of Cm(III) loaded Ca-montmorillonite and Na-illite were available from Part 1 of this work. De-convolution of the normalised fluorescence spectra measured at different pH values indicated three distinct Cm surface complexes, Cm complexes 1, 2 and 3 for both clay minerals, in agreement with model predictions, but with different distribution functions for the individual species. Under the assumption that Eu and Cm exhibit essentially the same hydrolysis and sorption behaviour, the Eu surface complexation constants were used to predict surface species distribution functions for Cm under the same experimental conditions used in the TRLFS measurements. Comparison of modelled and experimentally deduced species distributions indicated that for both clay minerals peak heights and widths of the three peaks did not correspond particularly well. It is shown that the calculated species distribution functions are sensitive to the values of the hydrolysis constants used in the calculations, whereas modelling the sorption edge measurements by applying the 2SPNE SC/CE approach is much less sensitive. By modifying the values of the hydrolysis constants within their uncertainty range and re-modelling the sorption edges, considerably better correspondence between the modelled and TRLFS species distribution functions was found. In particular, peak positions, heights and widths for the model predicted peaks for the ≡S^SOCm²⁺ and ≡S^SOCmOH⁺ species distribution, and those for Cm complexes 1 and 2 derived from TRLFS, were found to be very close for both clay minerals. However, discrepancies were still apparent between the profile for the calculated ≡S^SOEu(OH)₂° surface species and the Cm complex 3 species, especially in the case of Na-illite.     

Modelling the sorption of Mn(II), Co(II), Ni(II), Zn(II), Cd(II), Eu(III), Am(III), Sn(IV), Th(IV), Np(V) and U(VI) on montmorillonite: Linear free energy relationships and estimates of surface binding constants for some selected heavy metals and actinides

In solution thermodynamics, and more recently in surface chemistry, it is well established that relationships can be found between the free energies of formation of aqueous or surface metal complexes and thermodynamic properties of the metal ions or ligands. Such systematic dependencies are commonly termed linear free energy relationships (LFER). A 2 site protolysis non-electrostatic surface complexation and cation exchange (2SPNE SC/CE) model has been used to model “in house” and literature sorption edge data for eleven elements: Mn(II), Co(II), Ni(II), Zn(II), Cd(II), Eu(III), Am(III), Sn(IV), Th(IV), Np(V) and U(VI) to provide surface complexation constants for the strong sites on montmorillonite. Modelling a further 4 sets of sorption isotherms for Ni(II), Zn(II), Eu(III) and U(VI) provided complexation constants for the weak sites. The protolysis constants and site capacities derived for the 2SPNE SC/CE model in previous work were fixed in all of the calculations. Cation exchange was modelled simultaneously to provide selectivity coefficients. Good correlations between the logarithms of strong ^SK_x−1 and weak ^W1K_x−1 site binding constants on montmorillonite and the logarithm of the aqueous hydrolysis constants ^OHK_x were found which could be described by the following equations: Strong (≡S^SOH) sites:

SlogKX−1=8.1±0.3+(0.90±0.02)log^OHKX     

Sorption of Eu(III)/Cm(III) on Ca-montmorillonite and Na-illite. Part 1: Batch sorption and time-resolved laser fluorescence spectroscopy experiments

Sorption of Cm(III) and Eu(III) at trace concentrations onto Ca-montmorillonite (SWy-1) and Na-illite (Illite du Puy) has been studied under anaerobic conditions by batch sorption experiments and time-resolved laser fluorescence spectroscopy (TRLFS). Comparison of the results from spectroscopic and batch sorption experiments with Cm and Eu indicates the existence of outer-sphere complexes at pH <4 in the experiments with Na-illite (0.25 g/L solid; 2.5 × 10⁻⁷ mol/L Cm; 0.1 mol/L NaClO₄). In the case of Ca-montmorillonite, (0.25 g/L solid, 2.5 × 10⁻⁷ mol/L Cm or 10⁻⁶ mol/L Eu, 0.066 mol/L Ca(ClO₄)₂), Cm/Eu outer-sphere complexes do not form at significant levels due to the Ca²⁺ competition for the clay mineral cation-exchange sites. TRLFS spectra indicate the formation of inner-sphere surface complexes at pH >5 for both clay minerals. Five H₂O/OH⁻ molecules remain in the first metal ion coordination sphere of the sorbed Eu/Cm. Measured fluorescence lifetimes of sorbed Eu/Cm and peak deconvolution of Cm-spectra are consistent with the formation of surface complexes of the form ≡S-O-Eu/Cm(OH)_x^(2−x)(H₂O)_5−x. At pH ≥ 12 Cm becomes incorporated into a surface precipitate at the Ca-montmorillonite surface presumably composed of Ca(OH)₂ or calcium silicate hydrate. A dramatic shift of the fluorescence emission band by more than 20 nm and a clear increase in the fluorescence lifetime suggests the almost complete displacement of coordinated H₂O and OH⁻. The pH dependent Eu sorption data obtained in batch experiments are consistent with spectroscopic data on Eu and Cm within experimental uncertainties thus demonstrating the validity of Eu as a homologue for trivalent actinides. Parameterization of a two-site protolysis nonelectrostatic surface complexation and cation exchange model using the batch sorption data and spectroscopic results is discussed in Part 2 of this work.     

Sorption modelling on illite Part I: Titration measurements and the sorption of Ni, Co, Eu and Sn

In this study the physico-chemical, titration and sorption characteristics of Na-illite du Puy (Na-illite) have been measured and modelled. Samples of illite, collected in the region of le Puy-en-Velay, France, were purified and conditioned to the Na-form and physico-chemically characterised. Potentiometric titrations on suspensions of the Na-illite were carried out using a batch backtitration technique in 0.01, 0.1 and 0.5 M NaClO₄ background electrolytes from pH∼3 to ∼11.5 in an inert atmosphere glove box. The supernatant solutions from each titration experiment in each series were analysed for K, Mg, Ca, Sr, Si, Al, Fe and Mn. The titration data were modelled in terms of the protolysis of two amphoteric edge sites (S^W1OH and S^W2OH) without an electrostatic term. Sorption edges (solid/liquid distribution ratios versus pH at trace sorbate concentrations and constant ionic strength) were determined for the transitions metals Ni(II) and Co(II), the lanthanide Eu(III), and the heavy metal Sn(IV) on Na-illite with NaClO₄ as the background electrolyte under anoxic conditions (CO₂ ? 2 ppm, O₂ ? 2 ppm). The study thus encompasses a broad range of metals with different thermodynamic characteristics and with valence states ranging from II to IV. The results from the modelling of the titration data, in combination with a non electrostatic surface complexation and cation exchange sorption model were applied to quantitatively describe the uptake characteristics of the metals listed above on Na-illite. Since sorption edges were measured at trace concentrations, metal uptake was modelled as occurring on strong type sites (S^SOH) only. This sorption model, the two site protolysis non electrostatic surface complexation and cation exchange model (2SPNE SC/CE model) had been previously developed and used to describe metal uptake on montmorillonite.     

Sorption modelling on illite. Part II: Actinide sorption and linear free energy relationships

Sorption edge data for Ni(II), Co(II), Eu(III) and Sn(IV) [Bradbury M. H. and Baeyens B. (2009) Sorption modelling on illite. Part I: titration measurements and sorption of Ni(II), Co(II), Eu(III) and Sn(IV), Part I] on purified Na-Illite du Puy are available from some previous work, and some new measurements for Am(III), Th(IV), Pa(V) and U(VI) are presented here. All of these sorption edge measurements have been modelled with a 2 site protolysis non-electrostatic surface complexation and cation exchange (2SPNE SC/CE) sorption model for which the site types, site capacities and protolysis constants were fixed [Bradbury M. H. and Baeyens B. (2009), Part I]. In addition, two further data sets for the sorption of Am(III) and Np(V) on Illite du Puy, obtained from the literature, were also modelled in this work. Thus, surface complexation constants for the strong sites in the 2SPNE SC/CE sorption model for nine metals with valence states from II to VI have been obtained. A linear relationship between the logarithm of strong site metal binding constants, ^SK_x−1, and the logarithm of the corresponding aqueous hydrolysis stability constant, ^OHK_x, extending over nearly 35 orders of magnitude is established here for illite for these nine metals. Such correlations are often termed linear free energy relationships (LFER), and although they are quite common in aqueous phase chemistry, they are much less so in surface chemistry, especially over this large range. The LFER for illite could be described by the equation: where, “x” is an integer. A similar relationship has been previously obtained for montmorillonite, thus LFERs relating to the sorption on two of the most important clay minerals present in natural systems have been established. Such an LFER approach is an extremely useful tool for estimating surface complexation constants for metals in a chemically consistent manner. It provides a means of obtaining sorption values for radionuclides for which there are no measured values and thus allows gaps in missing sorption data to be filled. An ultimate goal of this approach is to develop a thermodynamic sorption database. This could then be used in radioactive waste management performance assessment studies to calculate sorption in natural systems, and thereby replace the current usage of single solid liquid distribution coefficients (K_d values) to describe radionuclide uptake. Finally, with the data now available, the 2SPNE SC/CE sorption model can be ported into reactive transport models allowing radionuclide migration to be calculated under spatially and temporally changing conditions.     

Modeling of Cs+ diffusion and retention in the DI-A2 experiment (Mont Terri). Uncertainties in sorption and diffusion parameters

In the DI-A2 experiment several non-reactive and reactive tracers were injected as a pulse in a packed-off borehole in the Opalinus Clay. Unlike the previous DI-A1 test, the design of the Teflon filter in the injection borehole forced the water to flow through the filter and the open space between the filter and the borehole wall (the filter itself did not act as a diffusion barrier between the circulating solution and the rock). The decrease in tracer concentration in the liquid phase was monitored during a period of a year. Afterwards, the borehole section was overcored and the tracer profiles in the rock were analyzed. A main interest of this experiment was to understand the chemical behavior of sorbing tracers: Cs⁺ (stable), ⁸⁵Sr²⁺, ⁶⁰Co²⁺ and Eu³⁺ (stable). The complete dataset (except for Eu³⁺ because of strong sorption to experimental equipment) was analyzed in a previous study with a 2D diffusion–reaction model and the derived diffusion and sorption parameters were compared with laboratory data. As in DI-A1, a difference by a factor of about 2 for sorption (magnitude of the Freundlich isotherm) was obtained between in situ and laboratory batch sorption experiments.Recent experimental and modeling studies have shown equivalent Cs⁺ sorption on intact and disaggregated Opalinus Clay samples. In view of these developments, new modeling of Cs⁺ diffusion and retention in the DI-A2 experiment has been performed using CrunchFlow. The calculations include transport by diffusion and a multisite cation exchange model to account for the retention of Cs⁺. The new results show that upscaling of Cs⁺ sorption from laboratory to field is no longer required. However, a difference in sorption by a factor of about 2 is still explained by the use of different versions of the same cation exchange model (a small difference in the selectivity coefficient for one type of site). This uncertainty in sorption leads to an uncertainty in the effective diffusion coefficient (D_e) for Cs⁺, also by a factor of 2 (2–4 × 10⁻¹⁰ m²/s). Clearly, the values of D_e obtained are correlated with the strength of sorption in the model, with stronger sorption leading to larger D_e values. Discrimination between the two versions of the exchange model is not possible when using only the results of the in situ test. Additionally, during early times (t < 10 days) the drop in Cs⁺ concentration in the circulation system is slower than expected. Due to the experimental setup, this slow decrease in concentration cannot be caused by the filter in the contact between borehole and rock. Poor mixing in the circulation system could explain this effect.     

Geochemical and hydrogeochemical evolution of groundwater in Ferdows area, northeast of Iran

The chemical analysis of 19 water wells in Ferdows area, Northeastern Iran, has been evaluated to determine the hydrogeochemical processes and ion concentration background in the region. In the study area, the order of cation and anion abundance is Na⁺ > Ca²⁺ > Mg²⁺ > K⁺ and Cl^? > SO ₄ ^?2  > HCO₃ ^? > NO₃ ^?, respectively, and the dominating hydrochemical types are Na–Cl. Most metal concentrations in water depend on the mineral solubility, and pH, Eh, and salinity of the solution. Their ΣREE concentrations showed excellent correlations with parameters such as TDS and pH. North American Shale Composite (NASC)-normalized REE patterns are enriched in the HREEs relative to the LREEs for all groundwaters. They all have positive Eu anomalies (Eu/Eu* = 0.752–3.934) and slightly negative Ce anomalies (Ce/Ce* = 0.019–1.057). Reduction–oxidation, complexation, desorption, and re-adsorption alter groundwater REE concentrations and fractionation patterns. The positive Eu anomalies in groundwaters are probably due to preferential mobilization of Eu²⁺ relative to the trivalent REEs in the reducing condition.     

Surface complexation model for strontium sorption to amorphous silica and goethite

Strontium sorption to amorphous silica and goethite was measured as a function of pH and dissolved strontium and carbonate concentrations at 25°C. Strontium sorption gradually increases from 0 to 100% from pH 6 to 10 for both phases and requires multiple outer-sphere surface complexes to fit the data. All data are modeled using the triple layer model and the site-occupancy standard state; unless stated otherwise all strontium complexes are mononuclear. Strontium sorption to amorphous silica in the presence and absence of dissolved carbonate can be fit with tetradentate Sr²⁺ and SrOH⁺ complexes on the β-plane and a monodentate Sr²⁺complex on the diffuse plane to account for strontium sorption at low ionic strength. Strontium sorption to goethite in the absence of dissolved carbonate can be fit with monodentate and tetradentate SrOH⁺ complexes and a tetradentate binuclear Sr²⁺ species on the β-plane. The binuclear complex is needed to account for enhanced sorption at hgh strontium surface loadings. In the presence of dissolved carbonate additional monodentate Sr²⁺ and SrOH⁺ carbonate surface complexes on the β-plane are needed to fit strontium sorption to goethite. Modeling strontium sorption as outer-sphere complexes is consistent with quantitative analysis of extended X-ray absorption fine structure (EXAFS) on selected sorption samples that show a single first shell of oxygen atoms around strontium indicating hydrated surface complexes at the amorphous silica and goethite surfaces.     

Modelling Zn(II) sorption onto clayey sediments using a multi-site ion-exchange model

In environmental studies, it is necessary to be able to predict the behaviour of contaminants in more or less complex physico-chemical contexts. The improvement of this prediction partly depends on establishing thermodynamic models that can describe the behaviour of these contaminants and, in particular, the sorption reactions on mineral surfaces. In this way, based on the mass action law, it is possible to use surface complexation models and ion exchange models. Therefore, the aim of this study is (i) to develop an ion-exchange model able to describe the sorption of transition metal onto pure clay minerals and (ii) to test the ability of this approach to predict the sorption of these elements onto natural materials containing clay minerals (i.e. soils/sediments) under various chemical conditions. This study is focused on the behaviour of Zn(II) in the presence of clayey sediments. Considering that clay minerals are cation exchangers containing multiple sorption sites, it is possible to interpret the sorption of Zn(II), as well as competitor cations, by ion-exchange equilibria with the clay minerals. This approach is applied with success to interpret the experimental data obtained previously in the Zn(II)–H⁺–Na⁺–montmorillonite system. The authors’ research team has already studied the behaviour of Na⁺, K⁺, Ca²⁺ and Mg²⁺ versus pH in terms of ion exchange onto pure montmorillonite, leading to the development of a thermodynamic database including the exchange site concentrations associated with montmorillonite and the selectivity coefficients of Na⁺, K⁺, Ca²⁺, Mg²⁺, and Zn²⁺ versus H⁺.     

Investigation of the different binding edge sites for Zn on montmorillonite using P-EXAFS - The strong/weak site concept in the 2SPNE SC/CE sorption model

The 2 site protolysis non electrostatic surface complexation and cation exchange (2SPNE SC/CE) sorption model has been used over the past decade or so to quantitatively describe the uptake of metals with oxidation states from II to VI on 2:1 clay minerals; montmorillonite and illite. One of the main features in this model is that there are two broad categories of amphoteric edge sorption sites; the so called strong (S^SOH) and weak (S^W1OH) sites. Because of their different sorption characteristics, it was expected that the coordination environments of the surface complexes on the two site types would be different. Zn isotherm data on two montmorillonites, Milos and STx-1, were measured and modelled using the 2SPNE SC/CE sorption model. The results were used to define the most favourable experimental conditions under which Zn sorption was either dominated by the strong (S^SOH, ∼2 mmol kg⁻¹) or by the weak sites (S^W1OH, ∼40 mmol kg⁻¹). Highly oriented self-supporting films were prepared for polarised extended X-ray absorption fine structure (P-EXAFS) investigations.Montmorillonites often contain Zn incorporated in the clay matrix. The Zn bound in this form was quantified and the results from the analysis of the P-EXAFS spectra were taken into account in the interpretation of the spectra measured at low Zn loadings (∼2 mmol kg⁻¹) and medium Zn loadings (∼30 mmol kg⁻¹). The Zn spectra on the “strong sites” exhibited a pronounced angular dependency and formed surface complexes in the continuity of the Al-octahedral sheets at the montmorillonite edges. In contrast, the Zn “weak site” spectra showed only a weak angular dependency. The spectroscopic evidence indicates the existence of two distinct groups of edge surface binding sites which is consistent with a multi-site sorption model and in particular with the strong/weak site concept intrinsic to the 2SPNE S/CE sorption model.     

Estimates of the second dissociation constant of H2S from the surface sulfidation of crystalline sulfur

The adsorption of hydrogen sulfide (Γ_H2S) and protons (Γ_H⁺) on the surface of crystalline sulfur was investigated experimentally in H₂S-bearing solutions at temperatures of 25, 50, and 70°C, NaCl concentrations of 0.1 and 0.5 mol/dm⁻³ and log C_H⁺ values in the range −2.3 to −5. At all temperatures, the dominant process on the surface of the sulfur was deprotonation, and the average values of Γ_H2S were very close to the highest values determined for Γ_H⁺. This finding, combined with the lack of detectable proton adsorption in H₂S-free solutions, suggests that proton adsorption/desorption on the surface of sulfur occurs through formation of ≡ S − H₂S complexes in the presence of H₂S.We propose that this complexation represents sulfidation of the sulfur surface, a process analogous to hydroxylation of oxide surfaces, and that the sulfidation can be described by the reaction: ≡ S + H₂S = ≡SSH₂⁰ β° The deprotonation of the ≡ SH° complex occurs via the reaction: ≡ SSH₂⁰ = ≡SSH⁻ + H⁺ β⁻ Values of 2.9, 2.8, and 2.9 (± 0.23) were obtained for −log β⁻ at 25, 50, and 70°C, respectively. These data were employed to estimate the second dissociation constant for hydrogen sulfide in aqueous solutions using the extrapolation method proposed by Schoonen and Barnes (1988) and yielded corresponding values for the constant of 17.4 ± 0.3, 15.7, and 14.5, respectively. The value for 25°C is in very good agreement with the experimentally determined values of Giggenbach (1971) at 17 ± 0.1; Meyer et al. (1983) at 17 ± 1; Licht and Manassen (1987) at 17.6 ± 0.3; and Licht et al. (1990) at 17.1 ± 0.3.     

Sorption of Cs to micaceous subsurface sediments from the Hanford site, USA

The sorption of Cs⁺ was investigated over a large concentration range (10⁻⁹−10⁻² mol/L) on subsurface sediments from a United States nuclear materials site (Hanford) where high-level nuclear wastes (HLW) have been accidentally released to the vadose zone. The sediment sorbs large amounts of radiocesium, but expedited migration has been observed when HLW (a NaNO₃ brine) is the carrier. Cs⁺ sorption was measured on homoionic sediments (Na⁺, K⁺, Ca²⁺) with electrolyte concentrations ranging from 0.01 to 1.0 mol/L. In Na⁺ electrolyte, concentrations were extended to near saturation with NaNO_3(s) (7.0 mol/L). The sediment contained nonexpansible (biotite, muscovite) and expansible (vermiculite, smectite) phyllosilicates. The sorption data were interpreted according to the frayed edge-planar site conceptual model. A four-parameter, two-site (high- and low-affinity) numeric ion exchange model was effective in describing the sorption data. The high-affinity sites were ascribed to wedge zones on the micas where particle edges have partially expanded due to the removal of interlayer cations during weathering, and the low-affinity ones to planar sites on the expansible clays. The electrolyte cations competed with Cs⁺ for both high- and low-affinity sites according to the trend K⁺ >> Na⁺ ≥ Ca²⁺. At high salt concentration, Cs⁺ adsorption occurred only on high-affinity sites. Na⁺ was an effective competitor for the high-affinity sites at high salt concentrations. In select experiments, silver-thiourea (AgTU) was used as a blocking agent to further isolate and characterize the high-affinity sites, but the method was found to be problematic. Mica particles were handpicked from the sediment, contacted with Cs⁺_(aq), and analyzed by electron microprobe to identify phases and features important to Cs⁺ sorption. The microprobe study implied that biotite was the primary contributor of high-affinity sites because of its weathered periphery. The poly-phase sediment exhibited close similarity in ion selectivity to illite, which has been well studied, although its proportion of high-affinity sites relative to the cation exchange capacity (CEC) was lower than that of illite. Important insights are provided on how Na⁺ in HLW and indigenous K⁺ displaced from the sediments may act to expedite the migration of strongly sorbing Cs⁺ in subsurface environments.     

Geochemical and thermodynamic aspects of sorption of strontium on kaolinite dominated clay samples at Kalpakkam

The mobility of strontium in subsurface is largely influenced by sorption on to clay minerals. In the present study, kaolinite clay samples collected from the Kalpakkam nuclear plant site were employed to understand the sorption characteristics of strontium by batch method. The effect of several parameters such as time, strontium ion concentration, pH, temperature and ionic strength was investigated. The kinetic studies suggested pseudo-second-order mechanism. The experimental sorption data was fitted to Langmuir adsorption model for obtaining the sorption capacity of the sorbent. The maximum sorption capacity was 5.77 mg/g at 298 K and was found to increase with an increase in temperature. It was observed that the distribution coefficient (K _d) of strontium on clay increased as the pH of the solution increased. The distribution coefficient was found to decrease with an increase in concentration of Na⁺ and Ca²⁺ ions. This variation of K _d suggests that cation exchange is the predominant sorption process. It was also observed that sorption process is endothermic. The thermodynamic parameters such as ∆G ⁰, ∆H ⁰ and ∆S ⁰ were calculated. The negative values obtained for ∆G ⁰ indicated that the sorption of strontium on clay was spontaneous at all studied concentrations. ∆G ⁰ becomes more negative with an increase in temperature, suggests that the sorption process is more favorable at higher temperatures.     

Experimental study of europium (III) coprecipitation with calcite

Partitioning of Eu(III) in calcite, CaCO₃, was evaluated with the aim of collecting data on partition coefficients and to enhance understanding of the incorporation mechanisms. This information will aid in the interpretation of geological processes from rare Earth element (REE) data and in the use of Eu(III) as a chemical analogue for the trivalent actinides, particularly Am(III) and Cm(III). Coprecipitation experiments were carried out by the constant addition method at 25°C and P_CO2 = 1 atm. Eu(III) was strongly partitioned from the solution into calcite. For dilute solid solutions (X_Eu < 0.001), Eu partition coefficients were estimated to be 770 ± 290 and found to be independent of calcite precipitation rate in the range of 0.02 to 2.7 nmol mg⁻¹ min⁻¹. This could be explained by the approximately equal values of the Eu partition and adsorption coefficients. Several solid solution models were tested. A vacancy model for Eu₂(CO₃)₃-CaCO₃ is consistent with the experimental results and constraints on geometry for Eu fit in the calcite lattice. For low Eu content, vacancy density is independent of Eu concentration in the solid so logarithm of the ion activity product, log (Eu)²(CO₃²⁻)³, depends linearly on log X_Eu². The fit of the data to such a model is good evidence that Eu(III) is taken up as a true solid solution, not simply by physical trapping. A model using EuOHCO₃-CaCO₃ is also consistent with the uptake stoichiometry, but EuOH²⁺ substitution for Ca²⁺ would be expected to distort the calcite structure more than is compatible with such a high K_D. Several other models, including EuNa(CO₃)₂-CaCO₃, were abandoned because their stoichiometric relationships did not fit the experimental data.     

Uptake of uranium and trace elements in pyrite (FeS2) suspensions

Pyrite dissolution and interaction with Fe^(II), Co^(II), Eu^(III) and U^(VI) have been studied under anoxic conditions by solution chemistry and spectroscopic techniques. Aqueous data show a maximal cation uptake above pH 5.5. Iron (II) uptake can explain the non-stoichiometric [S]_aq/[Fe]_aq ratios often observed during dissolution experiments. Protonation data corrected for pyrite dissolution resulted in a proton site density of 9 ± 3 sites nm⁻². Concentration isotherms for Eu^(III) and U^(VI) sorption on pyrite indicate two different behaviours which can be related to the contrasted redox properties of these elements. For Eu^(III), sorption can be explained by the existence of a unique site with a saturation concentration of 1.25 × 10⁻⁶ mol g⁻¹. In the U^(VI) case, sorption seems to occur on two different sites with a total saturation concentration of 4.5 × 10⁻⁸ mol g⁻¹. At lower concentration, uranium reduction occurs, limiting the concentration of dissolved uranium to the solubility of UO₂(s).Scanning electron microscopy and micro-Raman spectrometry of U^(VI)-sorbed pyrite indicate a heterogeneous distribution of U at the pyrite surface and a close association with oxidized S. X-ray photoelectron spectroscopy confirms the partial reduction of U and the formation of a hyperstoichiometric UO_2+x(s). Our results are consistent with a chemistry of the pyrite surface governed not by Fe^(II)-bound hydroxyl groups, but by S groups which can either sorb cations and protons, or sorb and reduce redox-sensitive elements such as U^(VI).     

Surface complexation of arsenic(V) to iron(III) (hydr)oxides: structural mechanism from ab initio molecular geometries and EXAFS spectroscopy

Arsenic(V), as the arsenate (AsO₄)³⁻ ion and its conjugate acids, is strongly sorbed to iron(III) oxides (α-Fe₂O₃), oxide hydroxides (α-,γ-FeOOH) and poorly crystalline ferrihydrite (hydrous ferric oxide). The mechanism by which arsenate complexes with iron oxide hydroxide surfaces is not fully understood. There is clear evidence for inner sphere complexation but the nature of the surface complexes is controversial. Possible surface complexes between AsO₄ tetrahedra and surface FeO₆ polyhedra include bidentate corner-sharing (²C), bidentate edge-sharing (²E) and monodentate corner-sharing (¹V). We predicted the relative energies and geometries of AsO₄-FeOOH surface complexes using density functional theory calculations on analogue Fe₂(OH)₂(H₂O)_nAsO₂(OH)₂³⁺ and Fe₂(OH)₂(H₂O)_nAsO₄⁺ clusters. The bidentate corner-sharing complex is predicted to be substantially (55 kJ/mole) more favored energetically over the hypothetical edge-sharing bidentate complex. The monodentate corner-sharing (¹V) complex is very unstable. We measured EXAFS spectra of 0.3 wt. % (AsO₄)³⁻ sorbed to hematite (α-Fe₂O₃), goethite(α-FeOOH), lepidocrocite(γ-FeOOH) and ferrihydrite and fit the EXAFS directly with multiple scattering. The phase-shift-corrected Fourier transforms of the EXAFS spectra show peaks near 2.85 and 3.26 Å that have been attributed by previous investigators to result from ²E and ²C complexes. However, we show that the peak near 2.85 Å appears to result from As-O-O-As multiple scattering and not from As-Fe backscatter. The observed 3.26 Å As-Fe distance agrees with that predicted for the bidentate corner-sharing surface (²C) complex. We find no evidence for monodentate (¹V) complexes; this agrees with the predicted high energies of such complexes.     

Fe sorption onto nontronite (NAu-2)

The sorption of ferrous iron to a clay mineral, nontronite (NAu-2, a ferruginous smectite), was investigated under strictly anoxic conditions as a function of pH (3-10), Fe²⁺ concentration (0.01-50 mM), equilibration time (1-35 days), and ionic strength (0.01-0.5 M NaClO₄). The surface properties of NAu-2 were independently characterized to determine its fixed charge and amphoteric site density in order to interpret the Fe²⁺ sorption data. Fe²⁺ sorption to NAu-2 was strongly dependent on pH and ionic strength, reflecting the coupled effects of Fe²⁺ sorption through ion exchange and surface complexation reactions. Fe²⁺ sorption to NAu-2 increased with increasing pH from pH 2.5 to 4.5, remained constant from pH 4.5 to 7.0, increased again with further increase of pH from pH 7.0 to 8.5, and reached a maximum above pH 8.5. The Fe²⁺ sorption below pH 7.0 increased with decreasing ionic strength. The differences of Fe²⁺ sorption at different ionic strengths, however, diminished with increasing equilibration time. The Fe²⁺ sorption from pH 4.5 to 7.0 increased with increasing equilibration time up to 35 days and showed stronger kinetic behavior in higher ionic strength solutions. The kinetic uptake of Fe²⁺ onto NAu-2 is consistent with a surface precipitation mechanism although our measurements were not able to identify secondary precipitates. An equilibrium model that integrates ion exchange, surface complexation and aqueous speciation reactions reasonably well describes the Fe²⁺ sorption data as a function of pH, ionic strength, and Fe²⁺ concentration measured at 24 h of equilibration. Model calculations show that the species Fe(OH)⁺ was required to describe Fe²⁺ sorption above pH 8.0 satisfactorily. Overall, this study demonstrated that Fe²⁺ sorption to NAu-2 is affected by complex equilibrium and kinetic processes, likely caused by surface precipitation reactions.     

Laser Ablation ICPMS study of trace element partitioning between plagioclase and basaltic melts: an experimental approach

Plagioclase-melt partition coefficients (D) for 34 trace elements at natural concentration levels were determined experimentally in a natural MORB composition at atmospheric pressure using thin Pt-wire loops. Experiments were carried out at three temperatures (1,220, 1,200, and 1,180°C), and at three different oxygen fugacities (fO₂ = IW, QFM, air) in order to assess the effect of fO₂ on the partitioning of elements with multiple valence (Fe, Eu, Cr). Run products were analyzed by laser-ablation ICP-MS. Most trace element Ds increase slightly as temperature decreases, except for D _Zr, D _Fe, D _Eu and D _Cr that vary systematically with fO₂. Applying the Lattice Strain Model to our data suggests the presence of Fe²⁺ entirely in the octahedral site at highly to moderate reducing conditions, while Fe³⁺ was assigned wholly to the tetrahedral site of the plagioclase structure. Furthermore, we provide a new quantitative framework for understanding the partitioning behaviour of Eu, which occurs as both 2+ and 3+ cations, depending on fO₂and confirm the greater compatibility of Eu²⁺, which has an ionic radius similar to Sr, relative to Eu³⁺ in plagioclase and the higher Eu²⁺/ Eu³⁺ under reducing conditions. For petrogenetic basaltic processes, a combined fractionation of Eu²⁺–Sr and Fe–Mg by plagioclase has considerable potential as an oxybarometer for natural magmatic rocks.     

Diffusion and sorption of Cs, I and HTO in samples of the argillaceous Wakkanai Formation from the Horonobe URL, Japan: Clay-based modeling approach

Diffusion and sorption behaviors of cationic Cs⁺, anionic I⁻ and neutral HTO in samples of the Wakkanai Formation from the Horonobe underground research laboratory (URL), Japan, were investigated as a function of ionic strength (I) of synthetic groundwater by through-diffusion and batch sorption experiments and mechanistic modeling. The effective diffusivities (D_e) measured by through-diffusion experiments showed cation excess and anion exclusion effects, which were strongly dependent on I; D_e for Cs⁺ decreased as I increased, D_e for I⁻ showed the opposite dependency and D_e for HTO showed no dependence. The sorption of Cs⁺ measured by through-diffusion and batch sorption experiments were described by Freundlich isotherms with consistent parameters and decreased with I as a result of competitive ion exchange.Diffusion and sorption behaviors were interpreted by assuming the clay components of illite and smectite control diffusion and sorption mechanisms. The component additive (CA) sorption model, which includes illite and smectite contents and their ion exchange constants, provided a reasonable account of the Cs⁺ sorption trends measured as functions of I and Cs concentration. The diffusion model was developed by coupling the electrical double layer (EDL) model, describing the change of ionic concentrations (cation excess and anion deficit) and viscoelectric effects caused by electrostatic interaction at negatively charged clay surfaces, and a simplified pore model assuming one type of pore shape and includes their size distribution. When averaging the electrostatic effects by using the pore surface area distribution, the model could predict the cation excess and anion exclusion effects, and its dependence on I reasonably well. This result implies the nanoscale pores dominating the pore surface area can strongly impact on ionic diffusion in argillaceous rocks. The clay-based modeling approach described here provides a useful tool to predict ionic diffusion and sorption in argillaceous rocks.