Next Article in Journal
Integrated Hydrologic and Hydrodynamic Models to Improve Flood Simulation Capability in the Data-Scarce Three Gorges Reservoir Region
Next Article in Special Issue
Phosphate Uptake by Precipitation in Model Animal Wastewaters: Adjusting Ionic Strength and Ionic Composition to Maximize Phosphorus Removal
Previous Article in Journal
Long-Term Consequences of Water Pumping on the Ecosystem Functioning of Lake Sekšu, Latvia
Previous Article in Special Issue
Impact of Microplastic Fibers from the Degradation of Nonwoven Synthetic Textiles to the Magdalena River Water Column and River Sediments by the City of Neiva, Huila (Colombia)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissolution-Desorption Dynamics of Strontium During Elution Following Evaporation: pH and Ionic Strength Effects

by
William C. Weaver
1,
Tohren C. G. Kibbey
2 and
Charalambos Papelis
1,*
1
Department of Civil Engineering, New Mexico State University, Las Cruces, NM 88003, USA
2
School of Civil Engineering and Environmental Science, University of Oklahoma, Norman, OK 73019, USA
*
Author to whom correspondence should be addressed.
Water 2020, 12(5), 1461; https://doi.org/10.3390/w12051461
Submission received: 5 March 2020 / Revised: 17 May 2020 / Accepted: 18 May 2020 / Published: 20 May 2020
(This article belongs to the Special Issue Water Quality Impacts of Contaminant Transport and Transformation)

Abstract

:
Radioactive strontium-90 (90Sr2+) is a fission byproduct of uranium and plutonium production, and therefore understanding its environmental fate is of particular importance for predicting the evolution of long-term risk from historical releases. The nonradioactive strontium cation, Sr2+, is a chemical analog for 90Sr2+ that is often used in studies designed to understand the environmental behaviors of 90Sr2+. The focus of this work was on understanding the dynamics of remobilization of strontium following evaporation to dryness in porous media. Evaporation is ubiquitous in the unsaturated zone, and has the potential to significantly impact the dynamics of transport by driving adsorption or precipitation on solid surfaces. For this work, a series of transport experiments were conducted examining the behavior of strontium over a range of pH values, ionic strengths, and concentrations. Saturated transport experiments were conducted, followed by experiments designed to examine the release and transport following evaporation to dryness. Results show increasing saturated retardation with increasing pH, decreasing ionic strength, and decreasing concentration, with the concentration exhibiting the strongest effect. Breakthrough curves at low concentrations were also found to be consistent with significant rate-limited desorption. Remobilization elution curves measured following evaporation to dryness exhibited the high initial effluent concentrations, exceeding the influent strontium concentration, most likely caused by the initial dissolution and accumulation of strontium by the advancing solution. Concentrations at later times were found to be largely consistent with the dynamics of saturated transport for the systems studied.

1. Introduction

The strontium cation, Sr2+, in its normal isotopic composition based on natural abundance, is a chemical analog to 90Sr2+, radioactive strontium-90, a fission byproduct of uranium and plutonium production [1,2,3,4,5,6]. Both uranium and plutonium are part of the nuclear fuel cycle and hence 90Sr is a byproduct of the nuclear fuel cycle, as well as a byproduct of nuclear weapon testing; strontium-90 is a pure beta emitter, as is its daughter product 90Y (half-life 64.2 h) [6]. Because of its half-life of 29.1 years, strontium-90 is one of the most important radioisotopes in low-level radioactive wastes, given that the half-lives of the other 16 major radioactive isotopes of strontium are shorter than 65 days [2,3]. Because of its chemical similarity to calcium, 90Sr can readily replace calcium in the human body; in addition, 90Sr has a tendency to deposit in bone and blood-forming tissue, causing chronic sicknesses, such as leukemia and anemia [1,2,3,4,5]. The United States Environmental Protection Agency (EPA) has published a maximum contaminant level (MCL) of 4 millirem per year for photon radioactivity and beta particles from man-made radionuclides in drinking water [3]. Therefore, the transport of the constituents of nuclear waste resulting from nuclear energy production and Cold War nuclear tests is a major environmental and human health concern at many United States Department of Energy (DOE) sites, such as Hanford, Oak Ridge National Laboratory, Idaho National Laboratory, Nevada National Security Site, and the Savannah River Site [4,5,7,8,9,10,11,12].
The transport of metals ions, often regarded as analogs of cationic radionuclides, has been studied extensively [13,14,15,16,17]. The effect of pH and other geochemical conditions on metal ion transport has been modeled using reactive transport and surface complexation models both under field [16] and laboratory conditions [14,15]. Specifically, the significant and varying effects of ionic strength on anion and cation sorption on oxide surfaces, as a function of the sorption mechanism, have long been known [18,19] and continue to be observed and reported in more recent studies [13]. However, the majority of these studies involve column experiments under saturated or unsaturated conditions and do not consider the effects of saturation and drying cycles on metal ion mobility.
The work conducted here aimed to explore the dynamics of strontium transport, with an emphasis on the mobilization from porous media after evaporation to dryness. Evaporation is ubiquitous in the unsaturated zone near the ground surface, so liquids entering the environment at the ground surface in dry environments may undergo cycles of evaporation followed by remobilization with rainfall infiltration events. Based on previous work with organic contaminants [20], it was hypothesized that evaporation/infiltration cycles would be likely to have a substantial impact on the transport of strontium.
Because of the importance of strontium as a component of nuclear waste, a substantial quantity of work has been conducted examining its transport in the environment in both saturated and unsaturated porous media [4,6,8,9,10,21,22,23,24]. Studies have included efforts aimed at understanding how strontium previously immobilized in the subsurface can be released by natural or engineered phenomena (e.g., [6,9,10,23,24]); however, little work has examined remobilization of strontium from porous media following evaporation to dryness. Note that as the water evaporates from a porous medium, the concentrations of nonvolatile compounds in porewater increase. This process drives adsorption to solid surfaces and can potentially also lead to surface precipitation once concentrations exceed solubility limits. The work described here was designed to explore how this phenomenon ultimately impacts the dynamics of release following resaturation of porous media. The work included a range of saturated and unsaturated transport experiments designed to explore this phenomenon, as well as an evaluation of the results of the work using a numerical model developed in previous work [20].

2. Materials and Methods

2.1. Materials

US Silica (Berkeley Springs, WV, USA) F-65 quartz sand was used for all experiments conducted for the work. The quartz sand is a foundry-grade and high-purity sand. Based on X-ray diffraction (XRD) analysis, only trace amounts of meionite, a tectosilicate, could be identified as a separate phase in addition to quartz. No other oxide or hydroxide phases were detectable. The Sauter mean diameter (diameter of a particle with the same volume to surface area ratio as the particle considered) was measured to be 0.215 mm, while 90% of the sand, by mass, had diameters between 0.106 and 0.295 mm. The specific surface area was determined by nitrogen adsorption and the Brunauer, Emmett, and Teller (BET) model using a Micromeritics (Norcross, GA, USA) ASAP 2050 Xtended Pressure Sorption Analyzer, and found to be 0.0519 m²·g−1. The pH of the point of zero net surface charge (pHpzc) of the sand was found to be 2 using an Anton Parr (Ashland, VA, USA) Electrokinetic Analyzer for Solid Surface Analysis: SurPASS 3.
Strontium (Sr) was purchased from Inorganic Ventures (Christiansburg, VA, USA) in the form of a SrCO3 certified reference material (10,000 µg·mL−1 in 1% v/v nitric acid, 99.9937% reported purity), and was diluted as necessary for all experiments. Ionic strength was adjusted in all solutions prepared for the work through addition of sodium nitrate (NaNO3).

2.2. Experimental Procedures

The experiments conducted for this work consisted of saturated breakthrough curves, as well as remobilization elution curves exploring how strontium is released after it has been retained by evaporative drying. The experiments conducted were based on similar systems and procedures described by Normile, Papelis, and Kibbey [20], adapted for the analysis of strontium.
All experiments were conducted in a cylindrical column measuring 1.91 cm for the inner diameter and 2.54 cm in length, having a total internal volume of 6.72 cm3. The column was made from transparent acrylic. Prior to filling the column, sand was rinsed three times with deionized water to remove fines, and then oven-dried at 110 °C for 24 h. Stainless steel no. 40 mesh screens were used at the top and bottom of the column to support water-wet capillary barrier membranes. Ahlstrom-Munksjö (Leominster, MA, USA) paper filters with a 1.5-µm pore size were used as water-wet capillary barriers to allow water flow while preventing air flow. Membranes were sealed to the column by placing O-rings within a slot machined into the top and bottom of the column.
Experiments were initiated by wet packing approximately 11 g of quartz sand into a column. The typical porosity of the sand column was approximately 0.37. Saturated transport experiments were conducted at a flow rate of 0.50 mL·min−1 (0.47 cm·min−1 porewater velocity) in a horizontal orientation. A pH- and ionic strength-adjusted background solution was pumped through the column for at least 5 pore volumes, removing residual quartz sand fines and stabilizing the column pH. The effluent pH was continuously measured using a Mettler Toledo (Columbus, OH, USA) Seven Excellence pH meter, and deviations from influent pH were small (less than a tenth of a pH unit). Background electrolyte ions were not measured. Especially at the higher ionic strengths, given the nature of the interactions of these major electrolyte ions with the quartz surface, it is not likely that any measurable fraction would be retained by the sand. It should also be noted that, as mentioned earlier, the sand was carefully rinsed three times before the experiments, so that any impurities that could interfere with the ions of interest were removed before the experiments were initiated. A bromide (conservative tracer) breakthrough curve was first obtained, and then a saturated breakthrough experiment was conducted using a solution with the specified strontium concentration, ionic strength, and pH. The effluent pH was monitored throughout each stage of all experiments.
Analysis of strontium concentrations in this and all subsequent experiment stages was conducted through the use of a PerkinElmer (Waltham, MA, USA) SCIEX ELAN DRC-e inductively coupled plasma mass spectrometer (ICP-MS). Effluent samples were collected at fixed intervals and diluted as necessary with a 1% nitric acid solution to a concentration appropriate for strontium determination. A calibration curve, with a correlation coefficient (R2) value of at least 0.999, was generated each time the ICP-MS was operated. Yttrium at 75 µg·L−1 was used as an internal standard. The ICP-MS sampling parameters included speed at 20 rpm, sample flush at 60 s, read delay at 15 s, and wash at 45 s.
Following completion of the saturated breakthrough curve, the column was partially drained to facilitate evaporation. In order to drain the column, the column was placed in a vertical orientation, and the top of the column was linked to a computer-controlled servo-pressure regulator, which was supplied air by an air compressor. The column was drained to approximately 40% saturation in 45 min by applying a ramped pressure input, using the servopressure regulator system described by Chen and Kibbey [25] and Chen et al. [26]. The drainage process was monitored by measuring and recording the column mass using a Mettler (Toledo, Columbus, OH, USA) Precision Balance PL602E with a computer interface, and the mass of strontium leaving the column was monitored by collecting samples at regular intervals for analysis by ICP-MS.
Following drainage, the column was opened by removing the top and allowed to evaporate to dryness over 24 h or more. The evaporation process was continuously monitored by recording the column mass as a function of time until no further change in mass could be observed.
Following evaporation, a remobilization elution experiment was conducted. The column was closed and the pH- and ionic strength-adjusted background electrolyte solution was introduced to the dry column. Once the column filled and solution began eluting from the column, samples were collected at fixed intervals and analyzed for the strontium concentration and pH. Flow was continued until the effluent strontium concentration reached the analytical detection limit.

2.3. Finite-Difference Elution Model

To better understand the results of the experiments, a custom-designed finite difference transport model was used. The model was developed in a previous work and is described in detail by Normile, Papelis, and Kibbey [20]. The purpose of the model is not necessarily to simulate observed elution behavior but rather to examine the extent to which the elution behavior can be predicted from rate-limited desorption, as opposed to dissolution dynamics. The model is based on a standard one-dimensional single-species transport model, Equation (1), with rate-limited adsorption/desorption, Equation (2), capable of simulating saturated transport but also modified to allow simulation of transport during the resaturation process, where a wetting front is advancing into the column:
C t = D 2 C x 2 v C x M V q t ,
q t = k r e v ( K d C q ) ,
where C is the is aqueous concentration (mg·L−1), v is the porewater velocity (m·s−1), D is the dispersion coefficient (m2·s−1) related to the dispersivity, α , by D = α · v , M is the mass of sand in the column, and V is the pore volume of the column. q is the mass adsorbed to the sand per kg of sand (mg·kg−1), krev is the desorption rate constant (s−1), and K d is the distribution coefficient (L·kg−1) [20].
Application of the model involves first using it to simulate a saturated breakthrough curve, for the purpose of determining a desorption rate constant (krev) needed to match the saturated breakthrough curve. The desorption rate constant, krev, is a first-order parameter that describes rate-limited desorption during transport. For a given Kd (adsorption distribution constant), a transport curve with a greater krev exhibits greater deviation from equilibrium transport behavior, i.e., more effects of rate limitation. The krev parameter is then used to simulate the elution behavior that would be expected from an initially dry column, assuming all transport behavior is the result of rate-limited desorption; deviation from the simulated behavior can be taken as an indicator that rate-limited dissolution is impacting the observed release. The elution model is formulated as an implicit finite-difference model and uses a moving boundary approach (i.e., expansion of the simulated region during the transport simulation) to simulate water movement into the column for resaturation and elution. All parameters, procedures, and solution methods for the model as used in this work were the same as those given by Normile, Papelis, and Kibbey [20].

3. Results and Discussion

Figure 1 shows saturated breakthrough curves for 10−4 M strontium at two different ionic strengths, 0.1 (Figure 1A) and 0.01 M (Figure 1B) NaNO3. In each subfigure, breakthrough curves at three pH values are shown. Note that at the higher ionic strength (0.1 M, Figure 1A), very little retardation is observed at any of the pH values and strontium behaves largely as a conservative tracer. This result is consistent with many previously reported studies of strontium-saturated transport, where authors have observed higher mobility in the presence of higher background electrolytes (e.g., [4,6,9,23]). The primary reason for the observed behavior is likely competition for sorption sites on the solid surface, with higher concentrations of electrolyte preventing the adsorption of strontium. At the lower ionic strength (0.01 M, Figure 1B), more substantial retardation is observed, and the behavior is pH dependent, with the greatest retardation occurring at pH 9, and the least at pH 5. Note that this behavior has also been observed by others (e.g., [4,6]), and is likely due to the increasingly negative surface charge of the solid as pH increases beyond the pHpzc (pH 2 for the quartz sand used here), thereby increasingly favoring the sorption of the strontium cation with increasing pH.
Figure 2 illustrates the effect of the strontium concentration and ionic strength on saturated breakthrough, showing breakthrough curves at pH 9 for strontium concentrations of 10−4 M (Figure 2A), 10−5 M (Figure 2B), and 10−6 M (Figure 2C). In each subfigure, breakthrough curves at three ionic strengths (0.1, 0.01, and 0.001 M NaNO3) are shown. In each case, the breakthrough of strontium is retarded to the greatest extent in the system with the lowest ionic strength, consistent with the results in Figure 1. However, the most significant difference between systems is observed between the different concentrations, with strontium at 10−6 M exhibiting retardation almost two orders of magnitude greater than that at 10−4 M (note that the horizontal axes of Figure 2A, B, and C cover different ranges). This result means that strontium exhibits strongly nonlinear adsorption behavior, with the lower concentrations having a substantially stronger adsorption affinity. Others have also observed nonlinear adsorption for strontium, with higher adsorption at lower concentrations (e.g., [22]).
The difference in adsorption affinities as a function of the concentration is even more apparent in Figure 3, which shows the breakthrough curves at pH 9, and ionic strength of 0.001 M for all three concentrations on the same horizontal axis. The figure highlights the order-of-magnitude differences in affinities between the three concentrations. This behavior can be interpreted in terms of the ratio between the concentration of strontium cations in solution and available sorption sites. Compared to the 10−6 M strontium concentration, when the strontium concentration is increased by 100-fold to 10−4 M, similar retardation for the two concentrations could only be observed if strontium sorption were not limited by sorption site availability. If, therefore, at the higher strontium concentrations, the finite number of sorption sites on the sand are exhausted, strontium retardation is expected to fall dramatically.
Retardation coefficients (Rf) were calculated for all saturated breakthrough curves by integrating to find the area above each C/C0 vs. pore volume breakthrough curve. The calculated values are shown in Table 1. Note that the trends in the preceding paragraphs are also apparent in the Rf values shown in Table 1.
Figure 4 shows the elution curves for the release of strontium following resaturation from an initially dry state. The curves in Figure 4 all correspond to pH 9 and show elution for ionic strengths of 0.1 (Figure 4A), 0.01 (Figure 4B), and 0.001 M (Figure 4C) NaNO3. Each curve shows elution from a column following a saturated breakthrough curve at the specified initial concentration, followed by drainage and then evaporation to dryness. Note that pore volumes on the horizontal axes in Figure 4 correspond to the number of pore volumes elapsed since water began to elute from the column following resaturation.
Several notable features are apparent in Figure 4. First, it is interesting to note that C/C0 at the start of most of the elution curves is greater than 1.0. The reason for this is that the process of resaturating the column causes strontium to be dissolved and carried along by the wetting front moving through the column, leading to a concentration magnification at the start of the elution curve as water passes through the column. This behavior was previously reported for a range of organic compounds [20]. It is also apparent from Figure 4 that conditions that favor greater retention in saturated transport also correspond to longer elution times. For example, at 0.1 M NaNO3 ionic strength (Figure 4A), the two higher concentrations elute almost immediately, consistent with their low saturated retardation factors (Table 1, Figure 1A), while the lowest concentration (10−6 M) system, which exhibits a higher retardation factor (Table 1), elutes over hundreds of pore volumes.
As with the saturated transport, the most significant factor influencing retention during elution is the concentration, with the lowest initial concentration (10−6 M) experiments exhibiting the greatest retention in both saturated transport and elution. Furthermore, it is interesting to note that in the systems studied, the number of pore volumes needed for concentrations during elution to approach zero is generally similar to that needed for C/C0 to approach 1.0 during saturated transport. For example, the 10−6 M elution curve in Figure 4C (I = 0.001 M) exhibits long tailing, but C/C0 is approaching zero by around 300 pore volumes, while for saturated transport with the same system (Figure 2C), C/C0 approaches 1.0 after approximately 350 pore volumes. In contrast, previous work with strongly adsorbing organic solutes found substantially longer tailing in elution curves compared with saturated breakthrough curves [20].
To further explore the mechanisms underlying the elution behavior observed in Figure 4, a finite difference elution model developed in a previous work [20] was used. The model is designed to highlight the extent to which elution results from rate-limited desorption. To apply the model, it must first be fit to saturated breakthrough curves to determine desorption rate constants. The fitting process involves determining an Rf value by integration of the saturated breakthrough curve (Table 1), and then using this value to simulate breakthrough while adjusting the desorption rate constant (krev) to match the observed behavior. Figure 5 shows the fits for saturated breakthrough curves at the three strontium concentrations (10−4, 10−5, and 10−6 M) at an ionic strength of 0.01 M NaNO3 and pH 9. Breakthrough curves in Figure 5 are shown on a semi-log scale for clarity so all three fits can be observed. Of particular note from the results in Figure 5 is the fact that lower concentrations correspond to significantly lower desorption rate constants, an indication that rate-limited desorption is increasingly important as concentration decreases for this system.
Figure 6 shows comparisons between model-predicted and measured remobilization elution curves for the same systems shown in Figure 5, with pore volumes shown on linear (Figure 6A) and semi-log (Figure 6B) scales. Note from Figure 6A that the trends predicted by the model generally follow those observed in the elution curves, with more rapid elution observed at the higher concentrations and more gradual elution for the lowest concentration. In particular, the long tailing observed at high pore volumes for the 10−6 M case is reasonably well-predicted by the model. This is a contrast with previous work with organic compounds [20], which found that elution curves for more strongly sorbing compounds were not particularly well-predicted by the model. In the organic systems studied in that earlier work, it was observed that the most strongly sorbing compounds did not exhibit particularly low krev values, yet still exhibited high initial eluting concentrations and long tailing similar to that observed for the 10−6 M concentration in Figure 6A. The difference between the model and experiment in that earlier work suggested that rate-limited dissolution was likely important for the high-sorbing organic systems. In this work, at least to a first approximation, rate-limited desorption does appear to capture many of the major features of elution. One notable exception to this is the fact that the model does not capture the high initial C/C0 values observed in the experiment, a discrepancy most apparent in Figure 6B, where the semi-log scale accentuates discrepancies at low pore volumes. In general, the model only predicts a high initial C/C0 for systems that exhibit low sorption affinity (Rf less than approximately 1.5), so would not be expected for any of the systems in Figure 6, all of which exhibit Rf values of 6 or greater. The high initial C/C0 observed in experiments means that at least some fraction of the strontium released during elution is more mobile than can be predicted from saturated transport behavior alone.

4. Summary and Conclusions

The results of this work show that the transport of strontium in a natural sand is strongly pH and ionic strength dependent. Higher ionic strengths lead to more rapid transport (lower adsorption), an observation that is likely the result of competition for adsorption sites on the sand surface between strontium and the background electrolyte cation. At the highest ionic strength studied (0.1 M NaNO3), retardation was negligible regardless of pH, and strontium transport behaved largely as a non-sorbing conservative tracer. At lower ionic strengths, increasing pH led to increasing retardation, likely a result of the increasing charge density on the solid surface at pH values farther above the pHpzc (approximately pH 2 for quartz sand). The most significant increases in retention were observed with a decreasing strontium concentration, with the lowest concentration studied (10−6 M) exhibiting two orders of magnitude the retardation of the highest concentration studied (10−4 M). This result is an indicator of a strongly nonlinear adsorption isotherm, most likely caused by sorption site limitation. Saturated transport was also found to exhibit significant rate limitation, with lower concentrations more strongly impacted by rate-limited desorption. The dynamics of remobilization of strontium after evaporation to dryness were found to be largely consistent with the dynamics of rate-limited saturated transport, with the exception of an initial high concentration release that could not be predicted by rate-limited desorption and was consistent with the presence of a highly mobile fraction of the retained strontium.

Author Contributions

Conceptualization, T.C.G.K. and C.P.; methodology, T.C.G.K. and C.P.; investigation, W.C.W., T.C.G.K., and C.P.; writing—original draft preparation, W.C.W.; writing—review and editing, T.C.G.K. and C.P.; supervision, T.C.G.K. and C.P.; project administration, T.C.G.K. and C.P.; funding acquisition, T.C.G.K. and C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation under grant nos. 1336083 and 1336579. Additional funding for this research was provided by the New Mexico Water Resources Research Institute.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huo, L.; Qian, T.; Hao, J.; Liu, H.; Zhao, D. Effect of water content on strontium retardation factor and distribution coefficient in Chinese loess. J. Radiol Prot. 2013, 33, 791–807. [Google Scholar] [CrossRef] [PubMed]
  2. Huo, L.; Qian, T.; Hao, J.; Zhao, D. Sorption and retardation of strontium in saturated Chinese loess: Experimental results and model analysis. J. Environ. Radioact. 2013, 116, 19–27. [Google Scholar] [CrossRef] [PubMed]
  3. Qian, T.; Li, S.; Ding, Q.; Wu, G.; Zhao, D. Two-dimensional numerical modeling of 90Sr transport in an unsaturated Chinese loess under artificial sprinkling. J. Environ. Radioact. 2009, 100, 422–428. [Google Scholar] [CrossRef] [PubMed]
  4. Um, W.; Papelis, C. Sorption mechanisms of Sr and Pb on zeolitized tuffs from the Nevada test site as a function of pH and ionic strength. Am. Mineral. 2003, 88, 2028–2039. [Google Scholar] [CrossRef]
  5. Wang, L.; Wu, J.Q.; Hull, L.C.; Schafer, A.L. Modeling Reactive Transport of Strontium-90 in a Heterogeneous, Variably Saturated Subsurface. Vadose Zone J. 2010, 9. [Google Scholar] [CrossRef]
  6. Wallace, S.H.; Shaw, S.; Morris, K.; Small, J.S.; Fuller, A.J.; Burke, I.T. Effect of groundwater pH and ionic strength on strontium sorption in aquifer sediments: Implications for 90Sr mobility at contaminated nuclear sites. Appl. Geochem. 2012, 27, 1482–1491. [Google Scholar] [CrossRef]
  7. Albarran, N.; Missana, T.; Garcia-Gutierrez, M.; Alonso, U.; Mingarro, M. Strontium migration in a crystalline medium: Effects of the presence of bentonite colloids. J. Contam. Hydrol. 2011, 122, 76–85. [Google Scholar] [CrossRef]
  8. Pace, M.N.; Mayes, M.A.; Jardine, P.M.; McKay, L.D.; Yin, X.L.; Mehlhorn, T.L.; Liu, Q.; Gurleyuk, H. Transport of Sr2+ and SrEDTA2- in partially-saturated and heterogeneous sediments. J. Contam. Hydrol. 2007, 91, 267–287. [Google Scholar] [CrossRef]
  9. Rod, K.A.; Um, W.; Flury, M. Transport of Strontium and Cesium in Simulated Hanford Tank Waste Leachate through Quartz Sand under Saturated and Unsaturated Flow. Environ. Sci. Technol. 2010, 44, 8089–8094. [Google Scholar] [CrossRef]
  10. Wang, G.H.; Um, W. Facilitated strontium transport by remobilization of strontium-containing secondary precipitates in Hanford Site subsurface. J. Hazard. Mater. 2013, 248, 364–370. [Google Scholar] [CrossRef]
  11. Li, Y.; Tian, S.; Qian, T. Transport and retention of strontium in surface-modified quartz sand with different wettability. J. Radioanal. Nucl. Chem. 2011, 289, 337–343. [Google Scholar] [CrossRef]
  12. Saunders, J.A.; Toran, L.E. Modeling of radionuclide and heavy metal sorption around low- and high-pH waste disposal sites at Oak Ridge, Tennessee. Appl. Geochem. 1995, 10, 673–684. [Google Scholar] [CrossRef]
  13. McNeece, C.J.; Hesse, M.A. Challenges in Coupling Acidity and Salinity Transport in Porous Media. Environ. Sci. Technol. 2017, 51, 11799–11808. [Google Scholar] [CrossRef] [PubMed]
  14. Muniruzzaman, M.; Rolle, M. Impact of multicomponent ionic transport on pH fronts propagation in saturated porous media. Water Resour. Res. 2015, 51, 6739–6755. [Google Scholar] [CrossRef] [Green Version]
  15. Muniruzzaman, M.; Rolle, M. Experimental investigation of the impact of compound-specific dispersion and electrostatic interactions on transient transport and solute breakthrough. Water Resour. Res. 2017, 53, 1189–1209. [Google Scholar] [CrossRef] [Green Version]
  16. Kent, D.B.; Wilkie, J.A.; Davis, J.A. Modeling the movement of a pH perturbation and its impact on adsorbed zinc and phosphate in a wastewater-contaminated aquifer. Water Resour. Res. 2007, 43. [Google Scholar] [CrossRef] [Green Version]
  17. Prigiobbe, V.; Bryant, S.L. pH-Dependent Transport of Metal Cations in Porous Media. Environ. Sci. Technol. 2014, 48, 3752–3759. [Google Scholar] [CrossRef]
  18. Hayes, K.F.; Leckie, J.O. Modeling ionic strength effects on cation adsorption at hydrous oxide/solution interfaces. J. Colloid Interface Sci. 1987, 115, 564. [Google Scholar] [CrossRef]
  19. Hayes, K.F.; Papelis, C.; Leckie, J.O. Modeling ionic strength effects on anion adsorption at hydrous oxide/solution interfaces. J. Colloid Interface Sci. 1988, 125, 717–726. [Google Scholar] [CrossRef]
  20. Normile, H.J.; Papelis, C.; Kibbey, T.C.G. Remobilization Dynamics of Caffeine, Ciprofloxacin, and Propranolol following Evaporation-Induced Immobilization in Porous Media. Environ. Sci. Technol. 2017, 51, 6082–6089. [Google Scholar] [CrossRef]
  21. Koss, V.; Kim, J.I. Modeling-of-strontium-sorption-and-speciation-in-a-natural-sediment-groundwater-system. J. Contam. Hydrol. 1990, 6, 267–280. [Google Scholar] [CrossRef]
  22. Szenknect, S.; Ardois, C.; Gaudet, J.P.; Barthes, V. Reactive transport of 85Sr in a chernobyl sand column: Static and dynamic experiments and modeling. J. Contam. Hydrol. 2005, 76, 139–165. [Google Scholar] [CrossRef] [PubMed]
  23. Hull, L.C.; Schafer, A.L. Accelerated transport of (90)Sr following a release of high ionic strength solution in vadose zone sediments. J. Contam. Hydrol. 2008, 97, 135–157. [Google Scholar] [CrossRef] [PubMed]
  24. Chang, H.S.; Um, W.; Rod, K.; Serne, R.J.; Thompson, A.; Perdrial, N.; Steefel, C.I.; Chorover, J. Strontium and Cesium Release Mechanisms during Unsaturated Flow through Waste-Weathered Hanford Sediments. Environ. Sci. Technol. 2011, 45, 8313–8320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chen, L.; Kibbey, T.C.G. Measurement of air-water interfacial area for multiple hysteretic drainage curves in an unsaturated fine sand. Langmuir 2006, 22, 6874–6880. [Google Scholar] [CrossRef] [PubMed]
  26. Chen, L.; Miller, G.A.; Kibbey, T.C.G. Rapid pseudo-static measurement of hysteretic capillary pressure-saturation relationships in unconsolidated porous media. Geotech Test. J. 2007, 30, 474–483. [Google Scholar]
Figure 1. Breakthrough curves for 10−4 M Sr in 0.1 M (A) and 0.01 M (B) NaNO3. Experiments conducted at pH 5, 7, and 9 are shown.
Figure 1. Breakthrough curves for 10−4 M Sr in 0.1 M (A) and 0.01 M (B) NaNO3. Experiments conducted at pH 5, 7, and 9 are shown.
Water 12 01461 g001
Figure 2. Breakthrough curves at pH 9 for Sr concentrations of 10−4 M (A), 10−5 M (B), and 10−6 M (C). Experiments conducted at ionic strengths of 0.1, 0.01, and 0.001 M NaNO3 are shown.
Figure 2. Breakthrough curves at pH 9 for Sr concentrations of 10−4 M (A), 10−5 M (B), and 10−6 M (C). Experiments conducted at ionic strengths of 0.1, 0.01, and 0.001 M NaNO3 are shown.
Water 12 01461 g002
Figure 3. Comparison between breakthrough curves at pH 9, ionic strength 0.001 M NaNO3 for Sr at concentrations of 10−4 M, 10−5 M and 10−6 M.
Figure 3. Comparison between breakthrough curves at pH 9, ionic strength 0.001 M NaNO3 for Sr at concentrations of 10−4 M, 10−5 M and 10−6 M.
Water 12 01461 g003
Figure 4. Remobilization elution curves showing the release of strontium at pH 9 following resaturation after the evaporation of porewater. Experiments correspond to ionic strengths of 0.1 (A), 0.01 (B), and 0.001 M (C) NaNO3. Curves are shown for elution after saturated breakthrough at initial strontium concentrations (C0) of 10−4 M, 10−5 M, and 10−6 M.
Figure 4. Remobilization elution curves showing the release of strontium at pH 9 following resaturation after the evaporation of porewater. Experiments correspond to ionic strengths of 0.1 (A), 0.01 (B), and 0.001 M (C) NaNO3. Curves are shown for elution after saturated breakthrough at initial strontium concentrations (C0) of 10−4 M, 10−5 M, and 10−6 M.
Water 12 01461 g004
Figure 5. Experimental and simulated saturated breakthrough curves at pH 9, ionic strength 0.01 M NaNO3 for strontium at concentrations of 10−4 M, 10−5 M, and 10−6 M. Simulations were fit to experimental curves by using measured Rf values (Table 1) and adjusting the desorption rate constant (krev) values.
Figure 5. Experimental and simulated saturated breakthrough curves at pH 9, ionic strength 0.01 M NaNO3 for strontium at concentrations of 10−4 M, 10−5 M, and 10−6 M. Simulations were fit to experimental curves by using measured Rf values (Table 1) and adjusting the desorption rate constant (krev) values.
Water 12 01461 g005
Figure 6. Experimental and simulated remobilization elution curves at pH 9, ionic strength 0.01 M NaNO3, showing the release of strontium at pH 9 following resaturation after the evaporation of porewater. Simulations correspond to the parameters determined from the saturated breakthrough curve fits in Figure 5.
Figure 6. Experimental and simulated remobilization elution curves at pH 9, ionic strength 0.01 M NaNO3, showing the release of strontium at pH 9 following resaturation after the evaporation of porewater. Simulations correspond to the parameters determined from the saturated breakthrough curve fits in Figure 5.
Water 12 01461 g006
Table 1. Measured retardation factor (Rf) from saturated breakthrough curves for all combinations of conditions studied. Values were determined by integrating to find the area above the dimensionless breakthrough curves.
Table 1. Measured retardation factor (Rf) from saturated breakthrough curves for all combinations of conditions studied. Values were determined by integrating to find the area above the dimensionless breakthrough curves.
pH = 5pH = 7pH = 9
Initial Sr2+, C0 (M)Initial Sr2+, C0 (M)Initial Sr2+, C0 (M)
I (M)10−410−510−610−410−510−610−410−510−6
0.11.51.640.31.41.977.51.42.474.3
0.012.68.293.04.310.7103.16.214.1126.7
0.0015.318.335.25.119.5121.87.128.5185.8

Share and Cite

MDPI and ACS Style

Weaver, W.C.; Kibbey, T.C.G.; Papelis, C. Dissolution-Desorption Dynamics of Strontium During Elution Following Evaporation: pH and Ionic Strength Effects. Water 2020, 12, 1461. https://doi.org/10.3390/w12051461

AMA Style

Weaver WC, Kibbey TCG, Papelis C. Dissolution-Desorption Dynamics of Strontium During Elution Following Evaporation: pH and Ionic Strength Effects. Water. 2020; 12(5):1461. https://doi.org/10.3390/w12051461

Chicago/Turabian Style

Weaver, William C., Tohren C. G. Kibbey, and Charalambos Papelis. 2020. "Dissolution-Desorption Dynamics of Strontium During Elution Following Evaporation: pH and Ionic Strength Effects" Water 12, no. 5: 1461. https://doi.org/10.3390/w12051461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop