Macroscale Solid-Phase Uranium Associations
A total acid digestion extraction method was used to measure total sediment-associated uranium. The digestion results of the solids indicate that the uranium concentrations range from 0.02 to 22 mg/kg in sediment. A summary of uranium concentrations (mg/kg) with depths (ft) for borings 852, 859, 860, and 858 is provided in Fig. 2. Boring 852 is the background with a uranium concentration of 2 mg/kg above the water table and 1 mg/kg below the water table.
Borings 859 and 860 are in the source zone with above background uranium concentrations of 5.5 mg/kg at 8.5 ft depth in fine sand and 2.7 mg/kg in coarse sand at 15.5 ft depth, respectively, in the saturated zone. The tailings (in addition to 3.9 ft of sub-pile material) from the unsaturated zone of the former tailings area were removed as a part of the DOE reclamation program (U.S. Department of Energy 2016); hence, the unsaturated zones of these two borings are not included in this paper.
Boring 858 is in the plume zone with a uranium concentration of 0.72 mg/kg (below background) in a coarse sand layer at 13 ft depth below the water table. Although the concentration is below the background, it has been selected for this study to understand the uranium transport mechanism between the upgradient source zone to the downgradient plume zone in the saturated layers. This study focused only on the saturated layers. Hence, the unsaturated layer of boring 858 was not studied despite uranium being above the background concentration (7.8 mg/kg). The high uranium concentration (5 mg/kg) at the bottom of boring 858 is due to the higher clay content in that portion of the bedrock.
Thus, the macroscale uranium concentration data from total digestions indicate significant variation in the solid-phase uranium concentration between the former tailings area-source zone and plume zone (as shown with depths, stratigraphy, and water table location in Fig. 2). Based on the total uranium digestions; it is the finer sediments in the saturated zone that contain more uranium than the coarser fractions, similar to other sediments exposed to uranium mill tailings (Davis et al. 2004). The sediments to be analyzed (red rectangle in Fig. 2) in the microscopic methods were selected from the depths based on the solid-phase uranium concentrations from the total acid digestion method.
Microscale Uranium Associations: Uranium Hosts
The microscale analysis of uranium in this section includes findings from fission-track radiography, SEM and EDS. The combined microscopic analyses show solid-phase uranium associations from the upgradient source zone to the downgradient plume zone. Three types of uranium hosts in the contaminated sediments are identified in this study. In the source zone, two distinctly different hosts for uranium are identified in aquifer sediments: Al hydroxide and ferric hydroxide coatings. In the plume zone, uranium is hosted by apatite. The amorphous coatings and apatite minerals do not contain sufficient uranium to be detected using EDS in the electron microprobe. However, fission-track analysis collected in the same thin sections exhibited precise spatial distribution of uranium that worked as a guide in SEM/EDS analysis. It is to be noted that the backscattered electron (BSE) image, fission-track radiography image, transmitted light image, and elemental maps are generated from the same ROI.
Al hydroxide coating
Al-rich coating: The aquifer sediments from boring 859 (Fig. 1 and Fig. 2) located in the source zone contain Al-rich precipitate (Fig. 3 and Fig. 4), as described below. Figures 3 and 4 depict two different sample points from the same thin section (Online Resource Fig. S1) prepared from boring 859. The Al-rich precipitates are found to be a few microns in width and up to a few tens of microns long (Fig. 3). SEM and EDS allowed us to show that these occur both as individual precipitates (Fig. 3) and at the edges of coatings of feldspar (Fig. 3) and quartz (Online Resource Fig. S4). The dried gel appearance with a distinct shrinkage crack (Fig. 3 and Fig. 4) and a consistent color (Online Resource Fig. S3) indicate that the precipitate is Al-rich rather than clay that has a molted appearance (Johnson et al. 2021). The stronger signal of the Al elemental map than the Si map (Fig. 3 and Fig. 4), and the higher Al peak than Si, with trace amounts of Ca, S, and P in the EDS results (Fig. 3, Fig. 4, and Online Resource Fig. S11) at the region of the precipitate support that Al is the primary element in the precipitate and not a component of clay.
Identification of the phase of Al-rich coating: The exact phases of the Al-rich precipitate could not be identified using SEM-EDS alone. However, based on a study by White et al. (1984) in Riverton, Wyoming, and EDS data with high Al peaks than that of Si (Fig. 3 and Fig. 4), it can be postulated that the Al-rich precipitate is Al hydroxide. According to White et al. (1984), the low-pH process waters caused the dissolution of calcite and the production of CO2 during contact with the underlying calcite-rich aquifer. White et al. (1984) hypothesized that the buffering of the mill tailing effluents in the groundwater caused significant Al and Fe hydroxide precipitation, indicated by the dramatic decrease of dissolved Al and Fe concentrations at near-neutral pH conditions.
Al hydroxide as a uranium host: Fission-track technology pinpointed trace uranium associations with corresponding hosts in the thin sections by overlapping the fission-track maps with the BSE images and EDS elemental maps. The fission tracks of uranium occur along the Al hydroxide precipitate that revealed the co-association of uranium with Al hydroxide precipitates and coatings around feldspar and quartz (Fig. 3 and Fig. 4). Fe-rich grains (Fig. 3), and iron oxide (Fig. 4) are also observed in proximity to Al hydroxide precipitate. Corresponding image (b) and Fe elemental map in Fig. 3 and Fig. 4 show no uranium fission tracks with the Fe-rich grains. Ferric hydroxide phases are known to have a high affinity towards uranium and are assumed to be one of the primary sorbing phases of uranium (Duff and Amrhein 1996; Walter et al. 2003; Davis et al. 2004). This suggests that Al hydroxide coating is a preferential reservoir for trace-concentration uranium.
Mechanism with Al hydroxide coating: Uranium is known to coprecipitate with a variety of hydroxides and carbonates, including aluminum hydroxides (Luo et al. 2009). Boring 859 contains a low carbonate concentration (alkalinity total as CaCO3 is 314 mg/L) and a higher pH (>5) that are favorable to form uranium coprecipitate with Al hydroxide (Luo et al. 2009). According to Luo et al. (2020), this precipitate is stable against dissolution. Moreover, the fact that uranium is associated with Al hydroxide precipitates, but not Fe-rich phases indicates that the residual uranium in boring 859 sediments is also unlikely to release by desorption from Fe oxides under ambient conditions. However, the long-term stability of coprecipitated cases depends on weathering of soil minerals (Noubactep et al. 2004). Thus, the coprecipitated uranium might release into groundwater if the Al hydroxide-coated feldspar or quartz weathers. Nevertheless, the significantly higher sedimentary uranium concentration (5.5 mg/kg) (Fig. 2) than the dissolved uranium concentration (~0.09 mg/L) (U.S. Department of Energy 2016) indicates a relatively slower uranium release.
Ferric hydroxide coating
Fe-rich coating: The aquifer sediments from boring 860 (Fig. 1 and Fig. 2), which is also located in the source zone, contain Fe-rich coating (Fig. 5 and Fig. 6) in contrast to boring 859 which contain Al-rich coating. Figures 5 and 6 depict two different sample points (Online Resource Fig. S14) from the same thin section prepared from boring 860. The Fe-rich coatings are a few microns thick (Fig. 4 and Fig. 5) and appear around the edges of feldspar (Fig. 5 and 6). At least two morphologies of the Fe-rich coating are observed in the two sample points- purely amorphous (Fig. 5) and particulates (Fig. 6). A stronger Fe signal and weaker Al and Si signals observed in the coatings (Fig. 5 and Fig. 6, Online Resource Fig. S18, and S19) depict the coating as Fe-rich. However, the morphology of feldspar (Fig. 5 and Fig. 6) indicates that it is chemically weathered around the edges. A significant amount of Al and Si with high Fe peaks (Fig. 5 and Online Resource Fig. S21) on the coating around weathered feldspar indicates the possibility of Fe-rich coating being intermixed with weathered kaolinite from feldspar.
Identification of phases of the Fe-rich coating: Goldberg. (1989) showed that Fe-oxides carry a sufficient positive charge at low pH and can precipitate on clay surfaces. These coatings are stable at high pHs and form oxide phases separate from clays. Charlet et al. (2002) found that the surface coatings on clay particles frequently encountered in soil particles are very likely a product of heterogeneous oxidative precipitation of iron by a variety of natural electron acceptors. Moreover, the ferric hydroxide precipitates are prominent under former tailings areas (White et al. 1984). Although the exact phases of the Fe-rich coating could not be identified with the methods used in this study, however, based on the above discussion, it can be postulated that the Fe-rich coating with clay are ferric hydroxides.
Ferric hydroxide as a uranium host: The fission-track analysis (Fig. 5 and Fig. 6) and BSE images (Fig. 5 and 6) show a relatively high density of uranium fission tracks occurring with ferric hydroxide coating around feldspar. Although ferric hydroxide coating is intermixed with kaolinite, uranium is found to have more affinity towards ferric hydroxide coating than kaolinite. It is more evident when the ferric hydroxide coating occurred as particles (Fig. 6). Corresponding images (a), (b), and (c) in Fig. 6 show that uranium fission tracks are more concentrated where ferric hydroxide particles are apparent (Fig. 6). Kaolinite is known to exhibit much greater uranium sorption due to more available aluminol sites (Bachmaf and Merkel 2011). However, the comparative visualization in the density of fission tracks (Fig. 6) shows that the trace uranium has more affinity towards ferric hydroxide coating than clay. The corresponding image (a) and Fe elemental map in Fig. 5 and Fig. 6 show that no uranium fission tracks are associated with Fe when it is incorporated into an aluminosilicate (Fig. 5) or a Fe-rich crystalline grain (Fig. 6). Although, the grains are in proximity to the ferric hydroxide coating. The morphology and uneven distribution of the elements on the grain in Fig. 5 indicate a crack with a weathering surface that might have been exposed when the thin section was cut. The above discussion and the direct comparative visualization of uranium co-occurrence in the admixture of clays (Fig. 5 and Fig. 6) suggest that in a low temperature environment, the ferric hydroxide coating plays a significant role as a uranium reservoir.
Mechanism with ferric hydroxide coating: Deposition of Fe-rich particles on kaolinite (Fig. 5 and Fig. 6) would increase the surface area independently of the underlying material (Arias, 1995). Hence a significant amount of uranium might get sorbed onto ferric hydroxide coating on feldspar. Typically, when uranium is adsorbed onto aged ferric hydroxide coating, the mechanism is ion exchange which is almost entirely reversible (Noubactep et al. 2004). Thus, the remobilization of sorbed uranium into groundwater is possible via concentration-driven desorption. The aqueous uranium concentrations remain ~0.77 mg/L which is significantly above the mean concentration level (MCL) of 0.044 mg/l (U.S. Department of Energy 2016). Meanwhile, the solid-phase uranium is relatively close to the background concentration of 2 mg/kg (Fig. 2) in boring 860, indicating a significant release of uranium from the solid-phase to the aqueous phase. Concentration-driven desorption of uranium from amorphous ferric hydroxide coating is likely responsible for the continuous release of uranium in the 860 boring.
Apatite
The aquifer sediments from boring 858 located in the plume zone (Fig. 1 and Fig. 2) contain a calcium phosphate mineral incorporated into quartz (Fig. 7). The EDS spectrum 1 in fig. 6 with a high Ca and P peak with traces of F, but no U peak indicates that the mineral is apatite but not autunite. The morphology of apatite (Fig. 7) shows that it is present as a crystalline grain at the edges of quartz rather than an amorphous coating. Although the uranium peak is not identified in EDS spectra generated from apatite, a relatively high density of uranium fission tracks is found to be associated with apatite (Fig. 7). Apatite is one of the most promising phases to immobilize uranium for the long term by accommodating a significant amount of U(VI) in the apatite structure (Rakovan et al. 2002, Dangelmayr et al. 2022). Mackinawite (Fig. 7) is also observed in quartz fractures, suggesting an anoxic environment below the water table. However, uranium fission tracks are found only with apatite. Since uranium did not occur with reduced iron minerals and did not precipitate as autunite (Fig. 6), apatite is the main sink for uranium in the saturated layer of the plume zone. Since the sample has a low uranium concentration in sediments (Fig. 2), it indicates the insignificant presence of apatite-like uranium sorbing phases, thus limiting the retardation of uranium in the saturated zone. However, the uranium concentration (~0.643 mg/L) in groundwater is above MCL. Such elevated uranium concentration can be explained by uranium transport from the source zone and plume stability in groundwater.
Effect on the transport behavior of uranium
The groundwater plume at the Riverton, Wyoming processing site persists to this day despite the completion of surface remediation in 1989 with the removal of the tailings. The persistence of the high groundwater uranium concentrations is hypothesized to result from the secondary sources in the saturated zone of the source zone. In this study, unique hosts of solid-phase uranium were identified in the saturated layer of source zone-1. Al-hydroxide coatings on feldspar and quartz, and 2. ferric hydroxide coatings intermixed with clay. The two separate sinks for uranium in the two locations under the former tailings area in the source zone would exhibit different mobilization mechanisms and require different management strategies. Thus, identifying major sinks for uranium in natural sediment is essential for risk assessment. Uranium coprecipitated into Al hydroxides is relatively stable against dissolution and insensitive to concentration-driven desorption. Meanwhile, the highly sorptive ferric hydroxide coatings increase the overall sorption capacity of the host rock and sorb a significant amount of uranium. Eventually, the concentration-driven desorption process might release considerable amounts of uranium from the solid phase. However, the uranium release from Al hydroxide coating is hypothesized to be slower than from the ferric hydroxide coating. These secondary sources of the former tailings area were not considered in the original transport model, which misled the uranium transport prediction (Dam et al. 2015). The dual hosts in the uranium-contaminated sediments increase complexity in the modeling of an already complex natural environment. Moreover, the association of uranium with tailings-derived Al and Fe precipitates is not well documented in modeling databases or literature. To our knowledge, Allard et al. (1999) is the only study that have analyzed the uranium mechanism in a quaternary system (U-Si-Al-Fe). Allard et al. (1999) showed that the sorption structure of uranium in Fe-rich gels was not the same as in the Si/Al-rich gels. Thus, there is a need to identify the elements associated with uranium at the field scale to avoid an incorrect correlation between uranium and other elements, which might lead to a misleading interpretation of uranium mobility.
The release of uranium from the source zone transport downward in the plume zone, which results in groundwater concentrations exceeding regulatory limits. Although uranium is found to be associated with stable apatite, the uranium-associated minerals that might significantly retard uranium transport were also not observed in this layer. Thus, minimum retardation of uranium is observed in the saturated layer of the plume zone, which aids in plume persistence.