Skip to main content

Surface complexation model for strontium sorption to amorphous silica and goethite

Abstract

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 Sr2+ and SrOH+ complexes on the β-plane and a monodentate Sr2+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 Sr2+ 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 Sr2+ 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.

Strontium surface complexation equilibrium constants determined in this study combined with other alkaline earth surface complexation constants are used to recalibrate a predictive model based on Born solvation and crystal-chemistry theory. The model is accurate to about 0.7 log K units. More studies are needed to determine the dependence of alkaline earth sorption on ionic strength and dissolved carbonate and sulfate concentrations for the development of a robust surface complexation database to estimate alkaline earth sorption in the environment.

Background

Ion sorption to mineral and amorphous solids has long been recognized as a process that controls the composition of trace elements in water. This process is particularly important for the transport of contaminants in the Earth's surface environment where sorption may retard transport by removing the contaminant from a mobile aqueous phase to a more stationary solid phase. Efforts to describe sorption in complex geological settings has evolved from a purely empirical approach in which distribution coefficients (Kd) are a measure of the total amount of specific ion between the solid and aqueous phases for a complex solution and solid matrix specific to a contaminated site. Although this approach provides a direct measure of the ability of the solid matrix to sequester the contaminant from a specific solution, its empirical nature does not allow it to be applied outside of the specific parameters of the contaminated site. Another approach measures thermodynamic surface complexation constants which describe sorption as a series of specific reactions between dissolved ions and surface sites. In principle, thermodynamic data from several single mineral and element experiments can be combined to build a model that represents the complex systems found in nature, especially when coupled with aqueous speciation, mineral solubility, and kinetic databases. However an internally consistent surface complexation database for a wide range of ions and solids found in natural waters that capture the surface charge is still lacking [1]. Databases tend to adopt surface complexation models that account for surface charge if only one solid is available for uptake [2] and non-electrostatic models that ignore surface charge if multiple solids are available for uptake [3–5].

It is also important to be able to estimate adsorption constants for reactions between aqueous components and substrates for which data are lacking both so that surface complexation models can be applied to complex geochemical systems and to build a robust database. Towards this effort, Sverjensky and colleagues have applied the Born solvation and crystal-chemistry theory together with a site-occupancy standard state to develop a predictive triple layer surface complexation model for surface protonation, alkali, alkaline earth, heavy metal and anion sorption for aluminum, iron, manganese, silica, and titanium oxides/hydroxides [6–16].

We illustrate the need for a predictive surface complexation model by considering the role that strontium sorption may play for the safe disposal of radioactive waste. 90Sr is one of several fission products that are concentrated in nuclear weapon and energy reprocessing waste and may interact with several different oxides depending on the waste form and disposal environment. At the Hanford (Washington, USA) site, caustic liquid waste with high 90Sr concentrations was disposed in tanks buried below the subsurface. Some of these tanks have leaked into the subsurface, where the migration of strontium is dependent on both its interaction with natural minerals and the reaction products formed from the interaction of the waste liquid with the subsurface fluvial-glacial sedimentary deposits [17–20]. Future disposal of 90Sr may include solid waste forms of cement, glass, or ceramics [21–24]. The long-term disposal of 90Sr depends not only on the stability of the waste form, but also on the sorption of strontium leached from the waste form to possible secondary phases, such as calcite, amorphous silica, iron hydroxides and rutile produced by the degradation of the waste form and corrosion of steel canisters containing the waste.

In this paper new strontium sorption data to amorphous silica and goethite collected over a range of total strontium concentrations, pH, and dissolved carbonate concentrations are described using a surface complexation model that builds on and further calibrates Sverjensky's [16] predictive model for alkaline earth sorption. Surface complexation reactions are constrained with structural information inferred from spectroscopic analysis of strontium at the mineral-solution interface (this study, [25–28]).

Experimental methods

Starting materials

For experiments conducted in the absence of CO2, all reagents were prepared with freshly distilled and deionized water collected under a nitrogen atmosphere using a portable microwave still. The water was then transferred to a nitrogen atmosphere glove box and used to make FeCl3, SrCl2, KOH, and NaOH stock solutions from solids that were purged for 20 to 30 minutes under a nitrogen stream and weighed in the nitrogen atmosphere glove box.

For experiments conducted in the presence of CO2, stock solutions were made by dissolving reagent grade NaCl, SrCl2, FeCl3 and Fe(NO3)3 solids in distilled and deionized water. Commercial high purity NaOH, HNO3, and HCl stock solutions were also used to adjust the pH of the sorption experiments. Solutions were stored in sealed containers and were not continually exposed to the atmosphere. We define dissolved carbonate as the sum of dissolved aqueous carbon species in this paper. Sources for dissolved carbonate in the sorption experiments include diffusion of atmospheric CO2 when the suspensions were prepared, dissolved carbonate present in the NaOH stock solution used to adjust solution pH, and possibly carbonate sorbed to goethite when synthesized at atmospheric pCO2.

The amorphous silica used in the sorption experiments was synthetic silica gel (Mallinckrodt Silica), 100–200 mesh, (lot # 6512), with an average pore diameter of 150 Å. The gel was repeatedly cleaned ultrasonically with distilled and deionized water until the suspension yielded a clear supernatant after 10 minutes of settling. Cleaned gel was dried at 40°C for 24 hours and stored in a plastic container at room temperature. Surface area was 277 m2 g-1 determined by BET nitrogen gas adsorption. No effort was made to exclude atmospheric CO2 in the cleaning and drying procedure for amorphous silica.

Four lots of goethite were synthesized following protocols outlined in Schwertmann and Cornell [29]. For CO2-free sorption experiments, goethite (Lot 1) was synthesized using KOH and FeCl3•6H20 in a nitrogen atmosphere from reagents dissolved in CO2-free distilled and deionized water. After initial formation of iron hydroxide, the suspension was purged with nitrogen for 60 hours at 70°C to transform the hydroxide to goethite. It was then rinsed repeatedly to remove chloride using dialysis tubing and CO2-free distilled and deionized water. Goethite was dried under a nitrogen stream at 40°C and stored in a nitrogen atmosphere glove box. Mineralogy was confirmed by XRD. For sorption experiments prepared in atmospheric CO2, goethite (Lots 2 and 3) was synthesized following the same protocol except that no effort was made to exclude CO2. For Lot 4, goethite was prepared from Fe(NO3)3 instead of FeCl3•6H2O and no effort was made to exclude CO2. Surface areas determined by BET nitrogen gas adsorption were: Lot 1 = 37.9 m2 g-1; Lot 2 = 38.2 m2 g-1; Lot 3 = 37.9 m2 g-1; and Lot 4 = 27.7 m2 g-1. An average surface area of 37.8 m2 g-1 was used to model the strontium sorption to goethite in NaCl solutions.

Sorption experiments

Strontium sorption was measured in amorphous silica and goethite suspensions prepared in the presence and absence of atmospheric CO2 at 25°C from pH 6 to 10. Tables 1, 2, 3, 4, 5, 6 list total surface area, sorption atmosphere, ionic strength, and initial and final solution composition for each experiment. Amorphous silica or goethite was mixed with a freshly prepared SrCl2/NaCl or Sr(NO3)2/NaNO3 solution of the desired concentration in polycarbonate test tubes. After the pH was adjusted, the tubes were sealed, shaken vigorously by hand, and then reacted for 2 or 14 days in a constant temperature orbital-shaker water bath at 200 rpm. The 14-day experiments were conducted in goethite suspensions to see if additional reaction was needed to precipitate strontium carbonate at higher pH. At the end of the experiment, the final pH of each solution was measured, a sample (2.5 ml) was taken, filtered (4.1 nm pore size), acidified with high purity HCl or HNO3 to prevent SrCO3 precipitation, and analyzed for strontium by inductively coupled plasma atomic emission spectrometry (ICP-AES) (detection limit = 10-7 molal; precision ± 2%). For many of the experiments prepared in the presence of atmospheric CO2, total dissolved carbonate was measured from a filtered sample using a carbon analyzer with an IR detector (detection limit = 5 × 10-5 molal). With this technique, dissolved carbon is purged with 11 N phosphoric acid and nitrogen gas. The resulting concentrations are under saturated with respect to atmospheric CO2 at higher pH. Additional control experiments with no solid present were done to check for strontium sorption to vessel walls and for precipitation of SrCO3(s). Strontium sorption uncertainty is calculated from the analytical uncertainty of the initial and final solution concentrations and the uncertainty associated with a small amount of strontium inherent in the substrate measured from substrate control experiments. For the experiments prepared in the absence of CO2, preparation, sampling, and reagent storage were done in a nitrogen atmosphere glove box.

Table 1 Solution analyses for strontium sorption to amorphous silica in nitrogen atmosphere, I = 0.1 M NaCl and T = 25°C.
Table 2 Solution analyses for strontium sorption to amorphous silica in air, I = 0.1 M NaCl and T = 25°C.
Table 3 Solution analyses for strontium sorption to amorphous silica in air, I = 0.005 M NaCl and T = 25°C.
Table 4 Solution analyses for strontium sorption to goethite in nitrogen atmosphere, I = 0.1 M NaCl and T = 25°C.
Table 5 Solution analyses for strontium sorption to goethite in air, I = 0.1 M NaCl and T = 25°C.
Table 6 Solution analyses for strontium sorption to goethite in air, I = 0.1 M NaNO3 and T = 25°C.

EXAFS sorption samples

Strontium sorption to amorphous silica experiments from pH 8 to 10 with dissolved CO2 and initial strontium concentrations of 10-3 M were analyzed with EXAFS. After reaction and centrifugation, supernatant liquids were removed and sorption samples were loaded as wet pastes into teflon sample holders with Kapton windows just prior to XAS analysis. For sorption samples collected at cryogenic temperatures, wet samples were quenched by immersion in liquid nitrogen and then placed in a helium cryostat in the beamline hutch.

For strontium sorption samples, EXAFS spectra were collected on wiggler beamline IV-3 at the Stanford Synchrotron Radiation Laboratory (SSRL). The incident beam was detuned to 50–70% of the maximum incoming intensity to reject higher-order harmonic reflections. The mid-point of the absorption edge of SrCO3(s) reference compound (set to 16105 eV) was used for energy calibration. Spectra were collected in fluorescence mode using a 13-element germanium array detector. For each sorption sample, 20–40 scans were collected to achieve an adequate signal.

Spectra in the EXAFS region were analyzed with the program EXAFSPAK [30]. Reference phase shift and amplitude functions used in non-linear least-squares fitting of experimental spectra were calculated using the ab initio program FEFF6 [31–33]. In non-linear least-squares fits of the sorption sample spectra, bond distance (R), backscatterer number (N), and the disorder or Debye-Waller term (σ2) were treated as adjustable parameters. The difference between theoretical and experimental threshold energies (ΔE0) was treated as a single adjustable parameter for all Sr-backscatterer shells [34]. Least-squares fits were performed on both filtered spectra of individual peaks in the radial structure functions (RSF) and on normalized χ(k) spectra with no significant differences in fit results. Detailed analyses of crystalline and hydrated strontium reference compounds and strontium in aqueous solution at ambient and cryogenic temperatures are given in our previous study [25]. These references allowed us to constrain adjustable EXAFS fitting parameters and to estimate errors in fit parameters based on empirical analysis (rather than using only the statistical errors derived from the least-squares fit). Our previous study [25], showed that anharmonic vibrational disorder of oxygen-ligated strontium compounds can be neglected because the third cummulate term (C3) of the EXAFS phase shift function is generally not significant above the error in fitted EXAFS distances (i.e., R ± 0.02 Å with and without C3). Also, we showed that rapid quenching of sorption samples for data collection at low temperature does not appear to introduce any new features into EXAFS spectra when compared with room temperature spectra. Based on our previous strontium EXAFS analyses, S02 was fixed at 0.92 and estimated empirical errors in fit parameters for first-shell Sr-O analysis are: R ± 0.02 Å; N ± 1 for N in the range of 6–12; σ2 ± 25% [25, 26].

Geochemical calculations

GEOSURF [35] and FITEQL [36] were used to fit specific surface complexation reactions to the experimental data. GEOSURF is tied to a thermodynamic database that automatically accounts for aqueous speciation and ionic strength during the sorption simulation. We used GEOSURF to fit strontium sorption to amorphous silica and goethite in carbonate-free suspensions and to amorphous silica in the presence of dissolved carbonate because SrCO3(aq) is negligible and dissolved carbonate is not known to sorb to amorphous silica to the best of our knowledge. We used FITEQL to fit strontium sorption data to goethite in the presence of dissolved carbonate to account for the significant carbonate sorption to the goethite surface because FITEQL allows the input of measured carbonate concentrations at each titration point for the speciation calculation. The pH dependence of measured dissolved carbonate concentrations in goethite suspensions was used to estimate dissolved carbonate concentrations when modeling strontium sorption for those experiments in which dissolved carbonate was not measured. The extended-Debye Hückel equation was used to correct for ionic strength effects. Aqueous equilibrium constants used in the calculations are listed in Table 7[37].

Table 7 Thermodynamic equilibrium constants used to account for aqueous speciation for strontium sorption model [37].

Results

EXAFS analysis

Absorption spectra were collected for samples of strontium sorbed to amorphous silica from solutions of 10-3 M SrCl2 and 0.1 M NaCl with dissolved CO2 from pH 8 to 10. Normalized χ(k) EXAFS spectra and Fourier transforms of the spectra are shown in Figure 1. Numerical fit results are given in Table 8. For all sorption samples collected at low temperature, there is only a single shell of oxygen backscatterers with a Sr-O distance of 2.60 ± 0.02 Å and a coordination number of 10 ± 1. Compared to a spectrum of strontium sorbed to silica gel collected at room temperature reported in our previous study [26], the fitted Sr-O distance at room temperature is slightly shorter (2.57 ± 0.02 Å) than that derived for low-temperature spectra. This small distance contraction was also noted for room- and low-temperature spectra of strontium sorbed to kaolinite [26] and probably results from a small anharmonic effect [25]. For strontium sorption on both silica gel and kaolinite, the fitted Sr-O distance is slightly shorter than the Sr-O distance (2.65 ± 0.02 Å) determined for aqueous Sr2+ in a 10-3 M SrCl2 solution [25].

Figure 1
figure 1

Normalized EXAFS spectra and corresponding radial structural functions for strontium sorbed to amorphous silica and for Sr2+(aq), strontianite (SrCO3(s)), and a Ca/Sr zeolite models. Solid lines are the data and dashed lines are fits to the data. For all sorption samples collected at low temperature, there is only a single shell of oxygen backscatterers indicated by one peak in the radial structural functions compared to multiple peaks for the model compounds reflecting Sr-Sr, Sr-C and Sr-Si/Al backscatters.

Table 8 EXAFS fit results for strontium sorbed to amorphous silica and goethite (R, N, σ2 and ΔE0 adjusted in fits).

Comparison of EXAFS sorption spectra indicates no change in strontium coordination with increasing pH and sorption. There is no evidence for silica backscatterers from the substrate, nor is there evidence for carbon or strontium backscatterers indicative of strontianite precipitation or other multi-nuclear sorption complexes. Evidence for scattering from atoms beyond the first coordination shell would be seen in multiple sine-wave oscillations in normalized spectra (i.e., "beat" patterns or shoulders on primary sine waves). Figure 1 compares the strontium sorption spectra to reference spectra for crystalline strontium carbonate (SrCO3(s)) and strontium in the calcium zeolite mineral heulandite (≈ 4500 ppm strontium substitution in the calcium site). In the calcium site in heulandite, strontium is eight-coordinated by oxygen, with three ligands of framework oxygen atoms from the mineral surface and five ligands of water extending into the zeolite channel. As shown here and in our previous study of strontium in zeolites [25], scattering from aluminum or silicon atoms in the zeolite framework is apparent in spectra up to a distance of 4.15 Å from central strontium because of direct bonding to the zeolite framework. Thus, strontium in zeolites is a good analog structure for inner-sphere complexation of strontium on silica gel, if it occurs. Although aluminum and silicon are relatively light backscatterers, they are easily identified as backscatterers in zeolite when strontium is partially dehydrated and bonds to the framework structure. Likewise, precipitation of strontianite is readily observed by backscattering from carbon and strontium atoms at distances of ≈ 4 Å or less, but is not seen in the sorption sample spectra, even for samples in which reacting solutions were supersaturated with respect to strontianite (Figure 2a). In our previous study of strontium reference compounds [25], we found that backscatterer atoms beyond the first oxygen coordination shell were apparent in normalized spectra when fitted values of σ2 (the Debye-Waller disorder parameter) were below ≈ 0.025 Å 2 (for N < 12 and R > 3 Å). We did not collect EXAFS spectra on samples of strontium sorbed to amorphous silica in the absence of CO2 because the bulk sorption behavior was the same with CO2 present. Nor did we collect EXAFS spectra on strontium sorption in 5 × 10-3 M NaCl solutions.

Figure 2
figure 2

ΔGr for final solution compositions with respect to strontianite (SrCO3) in which both dissolved strontium and carbonate were measured from amorphous silica (A) and goethite (B) suspensions. Solutions are supersaturated when ΔGr is > 0 and undersaturated when ΔGr < 0.

Figure 3 and Table 8 reproduce EXAFS spectra on strontium sorption to goethite with dissolved carbonate [26] to show changes in bonding at the surface as a function of solution pH. These results show that strontium forms a surface precipitate at pH 8.5, but not at higher pH where the solutions were more supersaturated with respect to strontianite for total strontium concentrations of 10-3 M (Figure 2b). For solutions with pH above 8.5 only Sr-O backscatters were detected indicating that strontium retains all or part of its hydration sphere when sorbed to the surface. This behavior was attributed to a maximum sorption of carbonate on goethite near pH 8.5 that nucleated a SrCO3-type surface precipitate, and decreasing carbonate sorption at higher pH that resulted in formation of hydrated Sr surface complexes [26].

Figure 3
figure 3

Normalized EXAFS spectra and corresponding radial structural functions for strontium sorbed to goethite and for Sr2+(aq) and SrCO3(s) models (reproduced from Sahai et al., 2000). Solid lines are the data and dashed lines are fits to the data. At pH 8.5 multiple peaks in the radial structural functions reflect Sr-O, Sr-C, and Sr-Sr backscatters that match SrCO3 (s). At higher pH, there is only a single shell of oxygen backscatterers indicated by one peak in the radial structural functions.

Macroscopic sorption experiments

Amorphous silica

The pH dependence of strontium sorption to amorphous silica in ~7 × 10-7 to 10-3 M SrCl2 and 0.1 M NaCl solutions with and without dissolved carbonate is shown in Figure 4 and Tables 1 and 2. Strontium sorption was near zero below pH 7.0 and increased with pH to about 80 % of the initial strontium concentration at pH 9.5 for all experiments. There was no measurable effect of dissolved CO2 on strontium sorption as shown by the similar pH dependence of sorption in suspensions with and without dissolved carbonate. Carbonate sorption to amorphous silica was not measured. We infer that strontium carbonate precipitation was negligible because it was not detected by EXAFS, even though final solution compositions for experiments with total strontium concentrations of 10-3 and 10-4 M are supersaturated or approach saturation with respect to strontianite (SrCO3) with increasing pH (Figure 2).

Figure 4
figure 4

Comparison of experimental and model results for strontium sorption to amorphous silica in 7 × 10-7 to 10-3 M SrCl2 and 0.1 M NaCl solutions with and without dissolved carbonate (A-E), and in 10-4 M SrCl2 and 0.005 M NaCl solutions with dissolved carbonate (F). Symbols are the experimental data and solid lines are fits to the data.

Strontium sorption to amorphous silica in 10-4 M SrCl2 and 5 × 10-3 M NaCl solutions prepared in the presence of dissolved carbonate is shown in Figure 4F and Table 3. These experiments were conducted to investigate the effect ionic strength on sorption. The strontium sorption increases at lower ionic strength as shown by a shift in the midpoint of the sorption edge (50% sorbed Sr) from pH 9.0 in 0.1 M NaCl solutions (Fig. 4B) to pH 7.2 in 5 × 10-3 M NaCl solutions (Fig. 4F). A similar shift in the strontium sorption edge has been observed for strontium sorbed to amorphous silica in 0.1 M and 0.01 M NaNO3, NaCl, and NaClO4 background electrolytes [28].

Goethite

The pH dependence of strontium sorption to goethite in 10-5 to 10-3 M SrCl2 and 0.1 M NaCl solutions without dissolved carbonate are shown in Figure 5(A–C) and Table 4. Similar to amorphous silica, strontium sorption to goethite has a broad pH sorption edge from pH 7 to 10. Enhanced strontium sorption to goethite between pH 7 and 8 with total strontium concentrations of 10-3 M was observed compared to lower total strontium concentrations of 10-4 and 10-5 M. The enhanced sorption at Sr = 10-3 M was reproduced in duplicate experiments and does not appear to be an experimental artifact.

Figure 5
figure 5

Comparison of experimental and model results for strontium sorption to goethite in 10-5 to 10-3 M SrCl2 and 0.1 M NaCl solutions without dissolved carbonate (A-C) and in 10-6 to 10-4 M SrCl2 and 0.1 M NaCl solutions with dissolved carbonate (D-F). Symbols are the experimental data and solid lines are fits to the data.

The pH dependence of strontium sorption to goethite with dissolved carbonate in 10-6 to 10-4 M Sr in 0.1 M NaCl (aged for 2 and 14 days) and 0.1 M NaNO3 (aged for 2 days) solutions is shown in Figures 5(D–F) and 6 and Tables 5 and 6. Similar to strontium sorption in carbonate-free systems, strontium sorption exhibits a broad sorption edge from pH 7 to 10. In the presence of dissolved carbonate, strontium carbonate precipitation has been observed at high total strontium concentrations [26] and may form strontium carbonate surface complexes because dissolved carbonate is known to sorb to goethite [38–42]. In our studies, strontium carbonate precipitation at the goethite surface appears to be negligible for Sr of 10-6 to 10-4 M. We assume that surface precipitation of SrCO3 (s) is minimal in these experiments because it was not detected in an EXAFS spectra from a strontium sorption sample prepared in the presence of atmospheric CO2 with pH = 9.9 and total Sr ~ 10-4 M where final solutions are highly supersaturated with respect to strontianite (Figures 2 and 3). It is possible that strontium carbonate precipitates at the goethite surface at lower pH even though the extent of supersaturation is less, because surface precipitates have been observed at pH 8.5 but not at higher pH in solutions with total Sr = 10-3 M [26].

Figure 6
figure 6

Total carbonate predicted at the goethite surface versus pH (A) in 0.1 M NaCl suspensions using measured carbonate concentrations (B) and strontium carbonate complexation model (Table 9).

Discussion

The strontium surface complexation model for amorphous silica and goethite presented here builds on a triple layer alkaline earth surface complexation model developed for oxides and hydroxides [16]. The model defines a limited number of surface complexation reactions that describe alkaline earth sorption over a wide range of ionic strength, pH, surface coverage, and oxide type by collectively fitting experimental data which investigated a limited set of these variables [43–54]. We model strontium sorption to amorphous silica and goethite using the same set of surface species to contribute to an internally consistent model for alkaline earth sorption that can be used for a wide range of solution and substrate compositions. Alkaline earth sorption is modeled primarily as a series of surfaces complexes placed on the β- or diffuse planes to account for the broad sorption edge, the dependence of sorption on ionic strength and available spectroscopic data. Although these complexes are not strictly identified as outer-sphere complexes by Sverjensky [16], their placement in the β- or diffuse planes is consistent with EXAFS analysis, which show that sorbed strontium remains hydrated at amorphous silica and goethite surfaces (this study, [26, 28]). The stoichiometry of surface complexes is further constrained by requiring the formation of some tetradentate complexes because strontium was found to bond to four surface oxygens on rutile with X-ray standing wave spectroscopy [27]. Although the exact stoichiometry could not be identified with spectroscopy, Sverjensky [16] found that the surface species (>SOH)2(>SiO-)2_MOH+ or (>SOH)2(>SiO-)2_M2+, where M stands for any alkaline earth, captured the broad pH dependence for much of the alkaline earth sorption.

Strontium surface complexation reactions were fit to the experimental data using equilibrium constants from Sverjensky [15] to account for surface protonation and sorption of the background electrolyte. All equilibrium constants were adjusted in accordance with the site-occupancy standard state [14]. Table 9 reports the total number of sites Ns, solid concentration Cs, surface area SA, and equilibrium constants normalized to 1 M standard state (K0) used in the fitting programs and the site-occupancy standard state (Kθ) used in the Born solvation and crystal-chemistry analysis. The site-occupancy standard state allows sorption experiments to be compared with one another independent of surface area, site density, and solid concentration. All complexes are mononuclear unless stated otherwise. The fitted constants were found to be accurate within an uncertainty of ± 0.3 log K [16].

Table 9 Equilibrium constants for alkaline earth sorption, where Log K0 is the equilibrium constant referenced to a 1 M standard state, Log Kθ is the equilibrium constant referenced to a site-occupancy standard state, Ns is the site density (# nm-2), As is the surface area (m2g-1), and Cs is the solid concentration (g L-1). The site-occupancy standard state assumes reference Ns = 10 sites nm-2 and As = 10 m2g-1.

Strontium sorption to amorphous silica

The fitted and experimental results are shown in Figure 4 for the sorption of strontium to amorphous silica in solutions with total Sr ranging from ~7 × 10-7 to 10-3 M, 0.1 M NaCl, and with and without dissolved carbonate, and for one experiment with total Sr = 10-4 M, 0.005 M NaCl, and dissolved carbonate. Strontium sorption to amorphous silica in the presence and absence of dissolved carbonate can be described with two tetradentate strontium complexes on the β-plane and one monodentate strontium complex on the diffuse plane (Table 9):

4 >SOH + Sr2+ = (>SOH)2(>SO-)2_ Sr2+ + 2H+ (1)

4 >SOH + Sr2+ +H2O = (>SOH)2(>SO-)2_ SrOH+ + 3H+ (2)

>SOH + Sr2+ = >SOH...Sr2+ (3)

In 0.1 M NaCl solutions, tetradentate Sr2+ and SrOH+ complexes capture the gradual sorption edge with increasing pH. The tetradentate Sr2+ complex, (>SOH)2(>SO-)2_Sr2+, dominates at near neutral pH where sorption is minimal. As pH and percent sorption increase, the tetradentate hydrolyzed SrOH+ complex, (>SOH)2(>SO-)2_ SrOH+, accounts for most of the strontium sorbed to the surface. In solutions with lower ionic strength (0.005 M NaCl), the sorption edge shifts to lower pH and a monodenate Sr2+ complex on the diffuse plane, >SOH...Sr2+, is needed to fit the data in addition to the tetradentate strontium complexes. Inclusion of other outersphere complexes (such as >SOH_ Sr2+, tetradentate SrCl+, and tetradentate SrClOH) failed to capture the enhanced strontium uptake at lower ionic strength. The dominance of the β-plane Sr2+ complexes at lower ionic strength over the β-plane SrOH+ complexes reflects the interplay between the charge in the β-plane and the charge of the surface complexes. As the ionic strength increases the β-plane SrOH+ complex becomes more dominant. The position of the surface complexes on the β- and diffuse planes suggests that strontium retains some or all of its waters of hydration at the amorphous silica surface and is an outer-sphere complex. The designation of sorbed strontium as outer-sphere is supported by the shift in the sorption edge from pH = 7.2 at I = 0.005 M NaCl to pH = 9 at I = 0.1 M NaCl (Figure 4B,F) and by EXAFS data showing only Sr-O bonding and coordination similar to aqueous Sr from pH 8.5 to 9.9 (Figure 1).

The strontium surface complexation model presented here is consistent with Sverjensky's [16] model fit for other alkaline earth sorption data for amorphous silica. Sverjensky [16] described Ca [50] and Mg [47] sorption data over a range of ionic strengths (0.001 to 0.1 N for Ca and 0.005 to 0.05 N for Mg) for a single surface coverage for each cation. The Ca model used the same type of complexes used to describe strontium sorption [>SOH...Ca2+, (>SOH)2(>SO-)2_ Ca2+, and (>SOH)2(>SO-)2_ CaOH+], and the Mg model used two of the three reactions [>SOH... Mg2+ and (>SOH)2(>SO-)2_ Mg2+].

Strontium sorption to goethite

Carbonate-free system

The fitted and experimental results are shown in Figure 5(A–C) and Table 8 for the sorption of strontium to goethite in carbonate-free solutions with total Sr ranging from 10-5 to 10-3 M and 0.1 M NaCl. In the absence of dissolved carbonate, strontium sorption to goethite can be fit with monodentate and tetradentate SrOH+ complexes and a tetradentate binuclear Sr2+ complex all on the β-plane:

>SOH + Sr2+ + H2O = >SO-_ SrOH+ + 2H+ (4)

4 >SOH + Sr2+ +H2O = (>SOH)2(>SO-)2_ SrOH+ + 3H+ (5)

4 >SOH + 2Sr2+ = (>SOH)2(>SO-)2_ Sr22+ + 2H+ (6)

Combination of reactions 4 and 5 fit most of the experimental sorption data for total strontium concentrations of 10-5 to 10-3 M, however they do not capture the enhanced strontium sorption between pH 7 and 8.5 for total strontium concentrations of 10-3 M (nor did the addition of tetradentate Sr+ or tetradentate SrOHCl complexes). We fit the data with a binuclear tetradentate strontium complex because the abundance of the binuclear strontium complex falls off with decreasing total strontium concentrations due to its second-order dependence on Sr2+(Equation 6). Unfortunately, there is no spectroscopic data in this pH range to confirm the presence of a polynuclear strontium complex. If a binuclear strontium complex forms at high surface coverage on goethite, it appears to be fairly unique for alkaline earth cation sorption. Polynuclear complexes were not needed to fit Ca and Mg data [45] with similar surface coverage [16]. Classification of strontium as an outer-sphere complex and its placement on the β-plane for all surface complexes is also consistent with EXAFS spectra which show only Sr-O bonding and coordination similar to aqueous strontium indicating that sorbed strontium retains waters of hydration at the mineral surface in carbonate free solutions [26].

The strontium model is consistent with surface complexation reactions fitted to other alkaline earth sorption data to goethite, but could be better refined with data collected over a range of ionic strength. In addition to the mono- and tetradentate MOH+ complexes, Sverjensky [16] found that tetradentate MgOHCl or BaOHNO3 surface complexes were needed to fit sorption data at high ionic strength of 0.5 N [45, 49] and that a tetradentate M2+ complex was needed to fit sorption of Ca at low surface coverage and low ionic strength [45, 52].

Carbonate system

The fitted and experimental results are shown in Figure 5(D–F) for the sorption of strontium to goethite with total Sr ranging from 10-6 to 10-4 M and 0.1 M NaCl in suspensions with dissolved carbonate. The strontium surface complexation model developed here is based on sorption data from experiments with total Sr from 10-6 to 10-4 M to avoid possible precipitation of strontium carbonate from experiments with the total Sr ~ 10-3 M (data not shown). The base model consists of the tetradentate and monodenate SrOH+ complexes (Equations 4 and 5) and carbonate surface complexation reactions and constants from Villalobos and Leckie [41] adjusted in accordance with the site-occupancy standard state [14] (Table 9) to account for carbonate sorption to goethite:

>SOH + CO32- + H+ = >SO-0.2_COO-0.8 + H2O (7)

>SOH + CO32- + H+ + Na+ = >SOCOONa + H2O (8)

>SOH + CO32- + 2H+ = >SOCOOH + H2O (9)

The binuclear tetradentate strontium surface complex (Equation 6) was not included in the calculations because its overall contribution is minimal at lower surface coverages (total Sr from 10-6 to 10-4 M) and because it led to convergence problems within FITEQL. Fits to the strontium sorption data require the addition of two monodentate strontium carbonate complexes. In the reactions below we maintain the carbonate stoichiometry and charge distribution between the 0- and β-planes as modeled by Villalobos and Leckie [41] for the strontium carbonate complexes.

>SOH + CO32- + Sr2+ = >SO0.2-_COOSrOH0.2+ + H2O (10)

>SOH + CO32- + H+ + Sr2+ = >SO0.2-_COOSr1.2+ + H2O (11)

Strontium carbonate surface complexes dominate the goethite surface from pH 6 to 10, with the non-carbonate SrOH+ complexes becoming important only at higher pH where less carbonate sorbs to the surface.

The strontium sorption model provides some insight into the pH dependence of SrCO3 precipitation at the goethite surface observed using EXAFS [26]. Strontium carbonate precipitate has been identified in the presence of goethite at pH 8.5, but not at pH greater than 8.7 in samples with the same surface loading (total Sr = 10-3 M) and for one sample at pH 9.9 with total strontium = 10-4 M. (Figure 3). Sahai et al [26] concluded that the absence of a strontium carbonate precipitate at pH above 8.7 was due to insufficient carbonate on the surface to nucleate the precipitate, based on a decrease in carbonate sorption on the goethite surface from pH 6 to 10 (using constants from VanGeen et al. [39]). Figure 6A is an example of the amount of carbonate predicted to sorb to the goethite surface using our strontium surface complexation model and measured dissolved carbonate. Strontium carbonate complexes comprise at most 3% of the total amount of carbonate on the surface. The amount of carbonate on the surface has a complex dependence on solution pH. Unlike the linear increase of dissolved carbonate concentrations with pH shown in Figure 6B, the amount of carbonate on the surface increases from pH 6 to 8, reaches a maximum from pH 8 to 9, and sharply decreases with increasing pH where the aqueous carbonate complexes dominate. From pH 6 to 9, there is 10 to 4 times more carbonate on the surface than dissolved in solution, illustrating the high affinity of the goethite surface for carbonate in near-neutral pH suspensions, where surface precipitation of strontium carbonate has been observed [26]. The stoichiometry of the carbonate complexes in our model further suggests that the >SO0.2-_COOSr1.2+ may be a precursor to surface precipitation, where as the hydrolyzed carbonate complex is not. Thus strontium carbonate precipitation at the goethite surface is inhibited despite having solutions that are supersaturated at higher pH.

Extension of strontium surface complexation model to strontium sorption to goethite in NaNO3electrolyte with dissolved carbonate

Figure 7 compares measured and predicted strontium sorption to goethite in suspensions containing dissolved carbonate, 0.1 M NaNO3 and total Sr ranging from 10-6 to 10-4 M. In this figure we show the total amount of strontium sorbed assuming an uncertainty of ± 0.3 log K for the strontium sorption constants (Equations 4, 5, 10, 11). Dissolved carbonate concentrations were estimated from concentrations measured as a function of pH in the NaCl experiments because they were not measured in the NaNO3 experiments. The surface complexation model is identical to the strontium carbonate surface complexation model described above with the exception that NO3- replaces Cl- for the complexation of the background electrolyte with the surface (Table 9). All constants have been adjusted in accordance with the site-occupancy standard state [14], where the total number of sites equals 16.4 nm-2, the solid concentration equals 10 g L-1 and the BET surface area of 27.7 m2g-1.

Figure 7
figure 7

Comparison of experimental and model results for strontium sorption to goethite in 10-6 to 10-4 M Sr(NO3)2 and 0.1 M NaNO3 solutions with dissolved carbonate (A-C). Symbols are the experimental data, solid lines are predictions made with the strontium surface complexation model (Table 9), dashed lines represent ± 0.3 uncertainty in log K for all strontium surface complexation reactions.

Application of the strontium carbonate surface complexation model developed in 0.1 M NaCl goethite suspensions to experiments conducted in 0.1 M NaNO3 goethite suspensions is a test of the site-occupancy standard state which allows sorption experiments with varying solid concentrations and specific surface areas to be compared. The solid concentration of the NaNO3 experiments was 25% of the solid concentration of the NaCl experiments and the goethite was synthesized from nitrate starting materials resulting in goethite with a surface area that was 75% of the goethite surface area synthesized from chloride starting materials. The strontium carbonate model slightly underpredicts the strontium sorption in NaNO3 goethite suspensions assuming an uncertainty of ± 0.3 log K suggested for surface complexation models [16].

Predictive model

Sverjensky and colleagues have applied Born solvation and crystal-chemistry theory using a site-occupancy standard state to develop a predictive model for surface protonation, alkali, heavy metal, alkaline earth, and anion sorption [6–16]. In this section we augment and re-calibrate Sverjenksy's [16] model for the prediction of alkaline earth speciation sorbed on oxide surfaces by including the fitted strontium surface complexation reactions for goethite and amorphous silica from this study.

Application of Born solvation and crystal-chemistry theory to metal sorption assumes that the standard Gibbs free energy of sorption (ΔGθr,m) depends on contributions from ion solvation (ΔGθs,m), electrostatic interactions between the sorbing ion and the surface sites (ΔGθai,m), and contributions specific to the sorbing ion (ΔGθii,m):

ΔGθr,m = ΔGθs,m + ΔGθai,m + ΔGθii,m (12)

such that a given surface complexation equilibrium constant (Log Kθr,m) can be expressed as:

Log Kθr,m = -ΔΩr,m/RT*(1/εs) - Bm(s/rm) + log K"ii,m (13)

The first term on the right hand side of equation 13 accounts for ion solvation, where ΔΩr,m is the Born solvation coefficient for the rth reaction containing the metal m and εs is the dielectric constant for the solid. The second term accounts for the repulsive interaction between the sorbing ion and near surface species, where s is the Pauling's bond strength of the metal-oxygen bonds in the bulk mineral, rm is the distance the sorbing ion is repulsed by the underlying metal in the solid due to short-range electrostatic interactions, and Bm is a constant characteristic of the surface reaction. The final term represents interactions intrinsic to the sorbing ion as well as solvation contributions from the interfacial dielectric constant and the electrostatic attractive interactions. Equation 13 can be reduced to:

Log Kθr,m = -ΔΩr,m/RT*(1/εs)+log K''ii,m (14)

for a given surface reaction if the repulsive interactions in the electrostatic term are minimal. Linear regressions of log Kθr,m versus 1/εs yield a slope equal to -ΔΩr,m/RT and a y-intercept equal to log K"ii,m (or Bm(s/rm) + log K"ii,m if repulsion interactions are important). This regression serves as a fundamental calibration for the predictive model, because it allows the equilibrium constant for a given surface reaction to be estimated for solids of varying dielectric constants.

In the absence of enough data to calibrate Equation 14, surface equilibrium constants can be calculated with ΔΩr,m and log K"ii,m derived from two additional regressions (Equations 15 and 17). Unknown values of ΔΩr,m can be estimated from the absolute solvation coefficient of a given surface complex, Ωabs r, and the effective electrostatic radius of the sorbing ion, Re,m. The Ωabs r is estimated from a regression of known ΔΩr,m versus Re,m:

ΔΩr,m= η/Re,m - Ωabs r (15)

where η = 166.027 kcal Å mol-1 [8]. Re,m is a function of the hydrated radii, rm,hydr, [55] and an empirical constant specific to the sorbed species, γm:

Re,m = rm,hydr + γm (16)

In the regression analysis γm is a variable used to produce a theoretical slope equal to 1 for ΔΩr,m versus η/Re,m. Similarly, unknown log K"ii,m for a given surface complex can be estimated from a linear regression of known values of log K"ii,m (the y-intercept in Equation 14 if repulsive interactions between the near surface species and sorbing ion are not important) versus the ion radius, rx,m, because log K"ii,m is assumed to be intrinsic to the sorbing ion and independent of differing solid properties. The resulting regression is:

Log K''ii,m = slope * rx,m + y-intercept (17)

Model Calibration

In this section we analyze the new strontium surface reactions for amorphous silica and goethite together with other alkaline earth surface reactions to re-calibrate Sverjensky's [16] predictive model for alkaline earth sorption. Strontium surface reactions on goethite and amorphous silica from this study are re-written as:

2>SOH + 2>SO- + Sr2+ + OH- = (>SOH)2(>SO-)2_SrOH+ (tet_SrOH+) (18)

2>SOH + 2>SO- + Sr2+ = (>SOH)2(>SO-)2_Sr2+ (tet_Sr2+) (19)

>SO- + Sr2+ + OH- = >SO-_SrOH+ (20)

>SOH + Sr2+ = >SOH...Sr2+ (21)

and the equilibrium constants used in the fits are converted from the 1 M standard state (Log K0) to the site-occupancy standard state (Log Kθ) using the following equations to be consistent with stoichiometry presented by Sverjensky [16]:

Log Kθtet_SrOH+ = Log K0tet_SrOH+ + 2pHzpc + ΔpKon + log Cs3(NsAs)4/(N†A†)2 + 14 (18a)

Log Kθtet_Sr2+ = Log K0tet_Sr2+ + 2pHzpc + ΔpKon + log Cs3(NsAs)4/(N†A†)2 (19a)

Log Kθ>SO-_SrOH+ = Log K0>SO-_SrOH+ + pHzpc + ΔpKon + log (NsAs)/(N†A†) + 14 (20a)

Log Kθ>SOH...Sr2+ = Log K0>SOH...Sr2+ + log (NsAs)/(N†A†) (21a)

We cannot evaluate the binuclear and strontium carbonate surface reactions with this model, because data are insufficient for calibration (Equations 6–11). Table 10 lists the surface equilibrium constants (Equation 18a-21a) from fits to Sr sorption to goethite and amorphous silica from this study and from fits for other alkaline earth cations [16]. The regression of log Kθr,m versus 1/εs for the formation of Sr, Ca, Mg, and Ba surface reactions are linear with correlation coefficients of 0.938 ≤ R ≤ 0.997 with at least three data points (Figure 8, Table 10). We include regressions for Ca, Mg and Ba in Figure 8 because the slopes and y-intercepts used in the calibration of the model are slightly different than those reported in Sverjensky [16] and to provide the reader with snapshot of the data used to calibrate the model. Table 11 lists ΔΩr,m and log K"ii,m derived for Sr and for Ca, Mg, and Ba from the regressions in Figure 8 (Equation 14), Ωabs r and fitting parameter γi from regressions of ΔΩr,m and η/Re,m in Figure 9 (Equation 15 and 16). Our calculations differ from Sverjensky [16] only in our use of alkaline earth hydration radii from Robinson and Stokes [55] for rm,hydr. In some cases Sverjensky [16] used rm,hydr as a fitting parameter to achieve reasonable linear regressions. We chose not do this because it introduces a second fitting parameter in addition to γi, which is adjusted to yield a theoretical slope equal to one in Equation 15. Additionally, we observe identical values of ΔΩr,m for the formation of (>SOH)2(>SO-)2_CaOH+ and (>SOH)2(>SO-)2_SrOH+ which is consistent with the two species having the same effective hydration radius. There is a greater difference in ΔΩr,m for the formation of (>SOH)2(>SO-)2_M2+ and >SO-_MOH+ for Ca and Sr; however it is only slightly greater than an uncertainty of ± 1.5 ΔΩr,m [16].

Table 10 Equilibrium constants (log Kθ) for alkaline earth surface complexation reactions.
Table 11 Regression slopes and intercepts from Figure 8 used to calibrate predictive alkaline earth surface complexation model.
Figure 8
figure 8

Log Kθfor (>SOH)2(>SO-)2_MOH+, (>SOH)2(>SO-)2_M2+, >SO-_MOH+, and >SOH...M2+ versus the dielectric constants for rutile (R), goethite (G), γ-alumina (A), quartz (Q), and amorphous silica (AS). Linear regressions are shown as lines of the same color as the data. Values of log Kθr,m for reactions involving Sr and goethite or Sr and amorphous silica are from this study and are shown as open symbols, all other log Kθr,m are from Sverjensky 2006 and are shown as solid symbols. The ''tet'' prefix refers to tetradentate surface sites.

Figure 9
figure 9

Born solvation coefficient (ΔΩr,M) versus the effective electrostatic radii (η/Re,m) for Sr, Ba, Ca and Mg for A: (>SOH)2(>SO-)2_MOH+, B: (>SOH)2(>SO-)2_M2+, and C: >SO-_MOH+. Values for ΔΩr,M are derived from regression slopes in Figure 8.

Figure 10 is a plot of the log K"ii for (>SOH)2(>SO-)2_MOH+, (>SOH)2(>SO-)2_M2+ and >SO-_MOH+ versus ion radii (Equation 17). The linear trends (R > 0.97) support the notion that repulsive interactions between alkaline earths and metals in the underlying substrate are minimal and the y-intercept in Figure 9 represents log K"ii,m for all three surface species. Sverjensky [16] estimated the repulsive interactions for M2+ surface reactions, but assumed they were insignificant for surface reactions involving MOH+. Results from EXAFS analyses indicate that strontium sorbs mostly as a hydrated surface complex, which is consistent with a large separation between the sorbed cation and metal cations in the substrate and thus minimal repulsion.

Figure 10
figure 10

Regression of log K"ii,m versus rx,m for the formation of (>SOH)2(>SO-)2_MOH+, (>SOH)2(>SO-)2_M2+, and >SO-_MOH+. Log K"ii,m are derived from regression y-intercepts in Figure 9. The "tet" prefix refers to the tetranuclear surface sites.

Our strontium surface reaction model fits nicely within the constraints of Sverjensky's [16] larger alkaline earth sorption model based on the regression analysis (Figures 8, 9, 10). The good correlation between Log Kθr,Sr versus 1/εs for the formation of >SO-_SrOH+ and (>SOH)2(>SO-)2_SrOH+ surface complexes suggests that equilibrium constants for these reactions can be predicted for many solids with varying dielectric constants. It is also possible to estimate equilibrium constants for (>SOH)2(>SO-)2_Sr2+ for other solids even though the regression is based only on two solids, because γ-alumina (εs = 10.4) and amorphous silica (εs = 4.6) span a fairly wide range in dielectric constant. It is not possible to estimate the formation of >SOH...Sr for a wide range of solids, because data are available only for amorphous silica (I = 0.005 N NaCl). Estimates of ΔΩr,m for the formation of the diffuse layer species (>SOH...M2+) from ΔΩr,m and a theoretical slope equal to one did not reproduce the limited data for this reaction. This species appears to be a very important complex for the uptake of alkaline earths on amorphous silica (and perhaps quartz) at low ionic strength, but it doesn't appear to be important for other oxides and hydroxides based on the available data.

Figure 11 compares the difference between equilibrium constants for the formation of (>SOH)2(>SO-)2_MOH+, (>SOH)2(>SO-)2_M2+ and >SO-_MOH+ fitted to sorption data and predicted directly by the regression of Log Kθr,m vs 1/εs (Equation 14) and by substitution of estimated values of ΔΩr,m and log K"ii,m into Equation 14 from regressions in Equations 15 and 17. We compare predicted constants from regression of Equation 14 if the surface reaction is calibrated with three or more different solids. On average the regression over predicts Δlog Kθr,m by 0.1 ± 0.6 (1σ). Prediction of all fitted equilibrium constants using regressed values for ΔΩr,m and log K"ii,m (Equation 15 and 17) substituted into Equation 14 yield an average log Kθ= 0.2 ± 0.7 (1σ). The overall uncertainty of the predicted model appears to be about twice that of log Kθvalues fitted to sorption data (log Kθ± 0.3).

Figure 11
figure 11

Comparison of log Kθr,m values fitted to sorption data and predicted from Born solvation and crystal chemistry theory using (A) regressions of log Kθr,m = -ΔΩr,m/RT*(1/εs) + log K''ii,m (Equation 14) and (B) substitution of estimated ΔΩr,m and log K''ii,m from Equations 15 and 17 to solve for log Kθr,m = ΔΩr,M/RT*(1/εs) + log K''ii,m (Equation 14). Colored boxes correspond to standard deviation of 0.6 and 0.7 log K units in A and B respectively.

Conclusion

Strontium sorption to amorphous silica and goethite can be modeled as a series of outer-sphere SrOH+ and Sr2+ complexes at tetradentate and monodentate surface sites. Reaction stoichiometry for strontium sorption is consistent with that used to model sorption of other alkaline earth metals [16] and allows strontium sorption to be evaluated over a wide range of solids in waters of varying composition. Surface equilibrium constants fit the sorption data to ± 0.3 log K units over a wide range of strontium surface coverage (total Sr = ~10-6 to 10-3M) in the presence and absence of dissolved carbonate. There are two key differences between strontium sorption to amorphous silica and goethite. Amorphous silica requires the formation of Sr2+ at the diffuse plane to account for enhanced sorption at low ionic strength, where as goethite does not. Dissolved carbonate does not appear to sorb to amorphous silica or impact the uptake of strontium to its surface, where as significant amounts of carbonate sorb to goethite and suggest the formation of strontium carbonate surface complexes to account for much of strontium uptake to goethite.

The overall alkaline earth model together with its predictive capability suggests that an additive approach can be used to describe sorption reactions in complex geochemical environments. The regression analysis done here suggests that alkaline earth sorption is largely a function of the solvation of the sorbing cation with minimal contributions between the sorbing cation and metals in the substrate as would be expected for outer-sphere sorption. Although the model can predict equilibrium constants for three non-carbonate surface reactions to within ± 0.7 log K, calibration is still fairly limited. There is a need for experimental data over wider range in ionic strength to determine the importance of β-plane MOHL and diffuse-plane M2+ surface species. It is also important to determine the role that carbonate and sulfate play on metal sorption to iron hydroxides and other oxides because both anions are abundant in the Earth's surface environment and will play a large role in the mobility of contaminants.

References

  1. Ochs M, Davis JA, Olin M, Payne TE, Tweed CJ, Askarieh MM, Altmann S: Use of thermodynamic sorption models to derive radionuclide Kd values for performance assessment: Selected results and recommendations of the NEA sorption project. Radiochimica Acta. 2006, 94: 799-785. 10.1524/ract.2006.94.9-11.779.

    Article  Google Scholar 

  2. Dzombak DA, Morel FMM: Surface Complexation Modeling: Hydrous Ferric Oxide. 1990, John Wiley and Sons

    Google Scholar 

  3. Davis JA, Coston JA, Kent DB, Fuller CC: Application of the surface complexation concept to complex mineral assemblages. Environmental Science and Technology. 1998, 32: 2820-2828. 10.1021/es980312q.

    Article  Google Scholar 

  4. Zavarin M, Bruton C: A non-electrostatic surface complexation approach to modeling radionuclide migration at the Nevada Test Site I. Iron hydroxides and calcite. Lawrence Livermore National Laboratory UCRL-TR-208673. 2004

    Google Scholar 

  5. Zavarin M, Bruton C: A non-electrostatic surface complexation approach to modeling radionuclide migration at the Nevada Test Site II. Aluminosilicates. Lawrence Livermore National Laboratory. UCRL-TR-208672. 2004

    Google Scholar 

  6. Sverjensky DA: Physical surface-complexation models for sorption at the mineral-water interface. Nature. 1993, 364: 776-780. 10.1038/364776a0.

    Article  Google Scholar 

  7. Sverjensky DA: Zero-point-of-charge prediction from crystal chemistry and solvation theory. Geochimica et Cosmochimica Acta. 1994, 58: 3123-3129. 10.1016/0016-7037(94)90184-8.

    Article  Google Scholar 

  8. Sverjensky DA, Sahai N: Theoretical prediction of single-site surfaceprotonation equilibrium constants for oxides and silicates in water. Geochimica et Cosmochimica Acta. 1996, 60: 3773-3797. 10.1016/0016-7037(96)00207-4.

    Article  Google Scholar 

  9. Sahai N, Sverjensky DA: Evaluation of internally consistent parameters for the triple-layer model by the systematic analysis of oxide surface titration data. Geochimica et Cosmochimica Acta. 1997, 61: 2801-2826. 10.1016/S0016-7037(97)00128-2.

    Article  Google Scholar 

  10. Sahai N, Sverjensky DA: Solvation and electrostatic model for specific electrolyte adsorption. Geochimica et Cosmochimica Acta. 1997, 61: 2827-2848. 10.1016/S0016-7037(97)00127-0.

    Article  Google Scholar 

  11. Criscenti LJ, Sverjensky DA: The role of electrolyte anions (ClO4-, NO3-, and Cl-) in divalent metal (M2+)adsorption on oxide and hydroxide surfaces in salt solutions. American Journal of Science. 1999, 299: 828-899. 10.2475/ajs.299.10.828.

    Article  Google Scholar 

  12. Criscenti LJ, Sverjensky DA: A single-site model for divalent transition and heavy metal adsorption over a range of metal concentrations. Journal of Colloid and Interface Science. 2002, 253: 329-352. 10.1006/jcis.2002.8529.

    Article  Google Scholar 

  13. Sverjensky DA: Interpretation and prediction of triple-layer model capacitances and the structure of the oxide-electrolyte-water interface. Geochimica et Cosmochimica Acta. 2001, 65: 3643-3655. 10.1016/S0016-7037(01)00709-8.

    Article  Google Scholar 

  14. Sverjensky DA: Standard states for the activities of mineral surface sites and species. Geochimica et Cosmochimica Acta. 2003, 67: 17-28. 10.1016/S0016-7037(02)01074-8.

    Article  Google Scholar 

  15. Sverjensky DA: Prediction of surface charge on oxides in salt solutions: Revisions for 1:1 (M+ L-) electrolytes. Geochimica et Cosmochimica Acta. 2005, 69: 225-257. 10.1016/j.gca.2004.05.040.

    Article  Google Scholar 

  16. Sverjensky DA: Prediction of the speciation of alkaline earths adsorbed on mineral surfaces in salt solutions. Geochimica et Cosmochimica Acta. 2006, 70: 2427-2453. 10.1016/j.gca.2006.01.006.

    Article  Google Scholar 

  17. Serne RJ, Bjornstad BN, Horton DG, Lanigan DC, Lindenmeier CW, Lindberg MJ, Clayton RE, LeGore VL, Geiszler KN, Baum SR, Valenta MM, Kutnyakov IV, Vickerman TS, Orr RD, Brown CF: Characterization of vadose zone sediments below the T Tank farm: Boreholes C4104 C4105, 299-W10-196 and RCRA borehole 299-W11-39. Pacific Northwest National Laboratory, PNNL-14849. 2004

    Google Scholar 

  18. Chorover J, Choi S, Amistadi MK, Karthikeyan KG, Crosson G, Mueller KT: Linking cesium and strontium uptake to kaolinite weathering in simulated tank waste leachate. Environmental Science and Technology. 2003, 37: 2200-2208. 10.1021/es025980x.

    Article  Google Scholar 

  19. Choi S, Amistadi MK, Chorover J: Clay mineral weathering and contaminant dynamics in a caustic aqueous system 1. Wet chemistry and aging effects. Geochimica et Cosmochimica Acta. 2005, 69: 4425-4436. 10.1016/j.gca.2005.04.003.

    Article  Google Scholar 

  20. Choi S, O'Day PA, Rivera NA, Mueller KT, Vairavamurthy MA, Seraphin S, Chorover J: Strontium speciation during reaction of kaolinite with simulated tank-waste leachate: Bulk and microfocused EXAFS analysis. Environmental Science and Technology. 2006, 40: 2608-2614. 10.1021/es051869q.

    Article  Google Scholar 

  21. Bao Y, Grutzeck MW, Jantzen CM: Preparation and properties of hydroceramic waste forms made with simulated Hanford low-activity waste. J Am Ceram Soc. 2005, 88: 3287-3302. 10.1111/j.1551-2916.2005.00775.x.

    Article  Google Scholar 

  22. Zhang Z, Blackford MG, Lumpkin GR, Smith KL, Vance ER: Aqueous dissolution of perovskite (CaTiO3): Effects of surface damage and [Ca2+] in the leachant. J Mater Res. 2005, 20: 2462-2473. 10.1557/jmr.2005.0294.

    Article  Google Scholar 

  23. Grambow B: Nuclear Waste Glasses – How durable?. Elements. 2006, 2: 357-364. 10.2113/gselements.2.6.357.

    Article  Google Scholar 

  24. Nyman M, Hobbs DT: A family of peroxo-titantate materials tailored for optimal strontium and actinide sorption. Chem Mater. 2006, 18: 6425-6435. 10.1021/cm061797h.

    Article  Google Scholar 

  25. O'Day PA, Newville M, Neuhoff PS, Sahai N, Carroll SA: X-ray absorption spectroscopy of strontium (II) coordination I. Static and thermal disorder in crystalline, hydrated, and precipitated solids and in aqueous solution. Journal of Colloid and Interface Science. 2000, 222: 184-197. 10.1006/jcis.1999.6621.

    Article  Google Scholar 

  26. Sahai N, Carroll SA, Roberts S, O'Day PA: X-ray adsorption spectroscopy of strontium(II) coordination II. Sorption and precipitation at kaolinite, amorphous silica, and goethite surfaces. Journal of Colloid and Interface Science. 2000, 222: 198-212. 10.1006/jcis.1999.6562.

    Article  Google Scholar 

  27. Fenter P, Cheng L, Rihs S, Machesky M, Bedzyk MJ, Sturchio NC: Electrical double-layer structure at the rutile-water interface as observed in situ with small-period x-ray standing waves. Journal of Colloid and Interface Science. 2000, 225: 154-165. 10.1006/jcis.2000.6756.

    Article  Google Scholar 

  28. Chen CC, Coleman ML, Katz LE: Bridging the gap between macroscopic and spectroscopic studies of metal ion sorption at the oxide/water interface: Sr(II), Co(II), and Pb(II) sorption to quartz. Environmental Science and Technology. 2006, 40: 142-148. 10.1021/es050356g.

    Article  Google Scholar 

  29. Schwertmann U, Cornell RM: Iron Oxides in the Laboratory: Preparation and Characterization. 1991, VCH, New York

    Google Scholar 

  30. George GN, Pickering IJ: EXAFSPAK: A Suite of Computer Programs for Analysis of X-Ray Absorption Spectra. 1993, Stanford Synchrotron Radiation Laboratory

    Google Scholar 

  31. Rehr JJ, Albers RC, Zabinsky SI: High-order multiple scattering calculations of X-ray absorption fine structure. Phys Rev Lett. 1992, 69: 3397-10.1103/PhysRevLett.69.3397.

    Article  Google Scholar 

  32. Rehr JJ: Recent developments in multiple scattering calculations of XAFS and XANES. Jpn J Appl Phys. 1993, 32: 8-12. 10.1143/JJAP.32.L8.

    Article  Google Scholar 

  33. Zabinsky SL, Rehr JJ, Ankudinov A, Albers RC, Eller MJ: Multiple-scattering calculations of x-ray-absorption spectra. Phys Rev B. 1995, 52: 2995-3009. 10.1103/PhysRevB.52.2995.

    Article  Google Scholar 

  34. O'Day PA, Rehr JJ, Zabinsky SI, Brown GE: Extended x-ray-absorption fine-structure (EXAFS) analysis of disorder and multiple-scattering in complex crystalline solids. J Am Chem Soc. 1994, 116: 2938-2949. 10.1021/ja00086a026.

    Article  Google Scholar 

  35. Sahai N, Sverjensky DA: GEOSURF: A Computer Program for Modeling Adsorption on Mineral Surfaces from Aqueous Solution. Computers and Geosciences. 1998, 24: 853-873. 10.1016/S0098-3004(97)00142-8.

    Article  Google Scholar 

  36. Herbelin AL, Westall JC: FITEQL: A Computer Program for Determination of Chemical Equilibrium Constants from Experimental Data. Oregon State University, Report 99-01. 1999

    Google Scholar 

  37. Johnson JW, Oelkers EH, Helgeson HC: SUPCRT92: A software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000°C. Computers and Geosciences. 1992, 18: 899-947. 10.1016/0098-3004(92)90029-Q.

    Article  Google Scholar 

  38. Zeltner WA, Anderson MA: Surface charge development at the goethite/aqueous solution interface: Effects of CO2adsorption. Langmuir. 1988, 4: 469-474. 10.1021/la00080a039.

    Article  Google Scholar 

  39. Van Geen A, Robertson AP, Leckie JO: Complexation of carbonate species at the goethite surface: Implications for adsorption of metal ions in natural waters. Geochimica et Cosmochimica Acta. 1994, 58: 2073-2086. 10.1016/0016-7037(94)90286-0.

    Article  Google Scholar 

  40. Villalobos M, Leckie JO: Carbonate adsorption on goethite under closed and open CO2 conditions. Geochimica et Cosmochimica Acta. 2000, 64: 3787-3802. 10.1016/S0016-7037(00)00465-8.

    Article  Google Scholar 

  41. Villalobos M, Leckie JO: Surface complexation modeling and FTIR study of carbonate adsorption to goethite. Journal of Colloid and Interface Science. 2001, 235: 15-32. 10.1006/jcis.2000.7341.

    Article  Google Scholar 

  42. Villalobos M, Trotz MA, Leckie JO: Variability in goethite surface site density: evidence from proton and carbonate sorption. Journal of Colloid and Interface Science. 2003, 268: 273-287. 10.1016/j.jcis.2003.07.044.

    Article  Google Scholar 

  43. Vuceta J: Adsorption of Pb(II) and Cu(II) on α-quartz from aqueous solutions: Influence of pH, ionic strength, and complexing ligands. Ph.D., California Institute of Technology. 1976

    Google Scholar 

  44. Huang C, Stumm W: Specific adsorption of cations on hydrous γ-Al2O3. Journal of Colloid and Interface Science. 1973, 43: 409-420. 10.1016/0021-9797(73)90387-1.

    Article  Google Scholar 

  45. Balistrieri LS, Murray JW: The surface chemistry of goethite (αFeOOH) in major ion seawater. American Journal of Science. 1981, 281: 788-806.

    Article  Google Scholar 

  46. Riese AC: Adsorption of radium and thorium onto quartz and kaolinite: A comparison of solution/surface equilibrium models. Ph. D., Colorado School of Mines. 1983

    Google Scholar 

  47. Kent DB, Kastner M: Mg2+ removal in the system Mg2+-amorphous SiO2-H2O by adsorption and Mg-hydroxysilicate precipitation. Geochimica et Cosmochimica Acta. 1985, 49: 1123-1136. 10.1016/0016-7037(85)90003-1.

    Article  Google Scholar 

  48. Jang HM, Fuerstenau DW: The specific adsorption of alkaline-earth cations at the rutile/water interface. Colloids and Surfaces. 1986, 21: 235-257. 10.1016/0166-6622(86)80094-4.

    Article  Google Scholar 

  49. Hayes KF, Leckie JO: Modeling ionic strength effects on cation adsorption at hydrous oxide/solution interfaces. Journal of Colloid and Interface Science. 1987, 115: 564-572. 10.1016/0021-9797(87)90078-6.

    Article  Google Scholar 

  50. Kosmulski M: Co-adsorption on mono- and multivalent ions on silica and alumina. Ber Bunsenges Phys Chem. 1994, 98: 1062-1067.

    Article  Google Scholar 

  51. Ali MA: The influence of simple organic acids on sorption of inorganic ions at the goethite/water interface. Ph. D., University of Pittsburgh. 1994

    Google Scholar 

  52. Ali MA, Dzombak DA: Effects of simple organic acids on sorption of Cu2+ and Ca2+ on goethite. Geochimica et Cosmochimica Acta. 1996, 60: 291-304. 10.1016/0016-7037(95)00385-1.

    Article  Google Scholar 

  53. Machesky M, Wesolowski DJ, Palmer DA, Ichiro-Hayashi K: Potentiometric titrations of rutile suspensions to 250°C. Journal of Colloid and Interface Science. 1988, 200: 298-309. 10.1006/jcis.1997.5401.

    Article  Google Scholar 

  54. Ridley MK, Machesky ML, Wesolowski DJ, Palmer DA: Calcium adsorption at the rutile-water interface: A potentiometric study in NaCl media to 250°C. Geochimica Cosmochimica Acta. 1999, 63: 3087-3096. 10.1016/S0016-7037(99)00236-7.

    Article  Google Scholar 

  55. Robinson RA, Stokes RH: Electrolyte Solutions. 1959, Butterworths, London, 126-

    Google Scholar 

Download references

Acknowledgements

We thank Catherine Hayden, Phil Toretto and Nita Sahai for their contributions early in this study. We also value the comments of two anonymous reviewers. This work was supported by the Department of Energy (DOE) Environmental Management Science Program (Grant #55249) and by subcontract (B335246) from Lawrence Livermore National Laboratory to Arizona State University. Work was done (partially) at SSRL, which is operated by the DOE, Office of Basic Energy Sciences. SAC and LJC acknowledge support from the U.S. Department of Energy, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. Sandia National Laboratories is a multi-program laboratory operated by the Sandia Corporation, a Lockheed Martin Company for the United States Department of Energy's National Nuclear Security Administration under Contract DE-AC04-94AL85000.

This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Susan A Carroll.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Carroll, S.A., Roberts, S.K., Criscenti, L.J. et al. Surface complexation model for strontium sorption to amorphous silica and goethite. Geochem Trans 9, 2 (2008). https://doi.org/10.1186/1467-4866-9-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1467-4866-9-2

Keywords