skip to main content


Title: Emerging investigator series: radium accumulation in carbonate river sediments at oil and gas produced water discharges: implications for beneficial use as disposal management
In the western U.S., produced water from oil and gas wells discharged to surface water augments downstream supplies used for irrigation and livestock watering. Here we investigate six permitted discharges on three neighboring tributary systems in Wyoming. During 2013–16, we evaluated radium activities of the permitted discharges and the potential for radium accumulation in associated stream sediments. Radium activities of the sediments at the points of discharge ranged from approximately 200–3600 Bq kg −1 with elevated activities above the background of 74 Bq kg −1 over 30 km downstream of one permitted discharge. Sediment as deep as 30 cm near the point of discharge had radium activities elevated above background. X-ray diffraction and targeted sequential extraction of radium in sediments indicate that radium is likely coprecipitated with carbonate and, to a lesser extent, sulfate minerals. PHREEQC modeling predicts radium coprecipitation with aragonite and barite, but over-estimates the latter compared to observations of downstream sediment, where carbonate predominates. Mass-balance calculations indicate over 3 billion Bq of radium activity ( 226 Ra + 228 Ra) is discharged each year from five of the discharges, combined, with only 5 percent of the annual load retained in stream sediments within 100 m of the effluent discharges; the remaining 95 percent of the radium is transported farther downstream as sediment-associated and aqueous species.  more » « less
Award ID(s):
1640634
NSF-PAR ID:
10093913
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Environmental Science: Processes & Impacts
Volume:
21
Issue:
2
ISSN:
2050-7887
Page Range / eLocation ID:
324 to 338
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Oil and gas (O&G) extraction generates large volumes of produced water (PW) in regions that are often water-stressed. In Wyoming, generators are permitted under the National Pollutant Discharge Elimination System (NPDES) program to discharge O&G PW for beneficial use. In one Wyoming study region, downstream of the NPDES facilities exist naturally occurring wetlands referred to herein as produced water retention ponds (PWRPs). Previously, it was found that dissolved radium (Ra) and organic contaminants are removed within 30 km of the discharges and higher-resolution sampling was required to understand contaminant attenuation mechanisms. In this study, we sampled three NPDES discharge facilities, five PWRPs, and a reference background wetland not impacted by O&G PW disposal. Water samples, grab sediments, sediment cores and vegetation were collected. No inorganic PW constituents were abated through the PWRP series but Ra was shown to accumulate within PWRP grab sediments, upwards of 2721 Bq kg −1 , compared to downstream sites. Ra mineral association with depth in the sediment profile is likely controlled by the S cycle under varying microbial communities and redox conditions. Under anoxic conditions, common in wetlands, Ra was available as an exchangeable ion, similar to Ca, Ba and Sr, and S was mostly water-soluble. 226 Ra concentration ratios in vegetation samples, normalizing vegetation Ra to sediment Ra, indicated that ratios were highest in sediments containing less exchangeable 226 Ra. Sequential leaching data paired with redox potentials suggest that oxic conditions are necessary to contain Ra in recalcitrant sediment minerals and prevent mobility and bioavailability. 
    more » « less
  2. Green Lake is the deepest natural inland lake in Wisconsin, with a maximum depth of about 72 meters. In the early 1900s, the lake was believed to have very good water quality (low nutrient concentrations and good water clarity) with low dissolved oxygen (DO) concentrations occurring in only the deepest part of the lake. Because of increased phosphorus (P) inputs from anthropogenic activities in its watershed, total phosphorus (TP) concentrations in the lake have increased; these changes have led to increased algal production and low DO concentrations not only in the deepest areas but also in the middle of the water column (metalimnion). The U.S. Geological Survey has routinely monitored the lake since 2004 and its tributaries since 1988. Results from this monitoring led the Wisconsin Department of Natural Resources (WDNR) to list the lake as impaired because of low DO concentrations in the metalimnion, and they identified elevated TP concentrations as the cause of impairment. As part of this study by the U.S. Geological Survey, in cooperation with the Green Lake Sanitary District, the lake and its tributaries were comprehensively sampled in 2017–18 to augment ongoing monitoring that would further describe the low DO concentrations in the lake (especially in the metalimnion). Empirical and process-driven water-quality models were then used to determine the causes of the low DO concentrations and the magnitudes of P-load reductions needed to improve the water quality of the lake enough to meet multiple water-quality goals, including the WDNR’s criteria for TP and DO. Data from previous studies showed that DO concentrations in the metalimnion decreased slightly as summer progressed in the early 1900s but, since the late 1970s, have typically dropped below 5 milligrams per liter (mg/L), which is the WDNR criterion for impairment. During 2014–18 (the baseline period for this study), the near-surface geometric mean TP concentration during June–September in the east side of the lake was 0.020 mg/L and in the west side was 0.016 mg/L (both were above the 0.015-mg/L WDNR criterion for the lake), and the metalimnetic DO minimum concentrations (MOMs) measured in August ranged from 1.0 to 4.7 mg/L. The degradation in water quality was assumed to have been caused by excessive P inputs to the lake; therefore, the TP inputs to the lake were estimated. The mean annual external P load during 2014–18 was estimated to be 8,980 kilograms per year (kg/yr), of which monitored and unmonitored tributary inputs contributed 84 percent, atmospheric inputs contributed 8 percent, waterfowl contributed 7 percent, and septic systems contributed 1 percent. During fall turnover, internal sediment recycling contributed an additional 7,040 kilograms that increased TP concentrations in shallow areas of the lake by about 0.020 mg/L. The elevated TP concentrations then persisted until the following spring. On an annual basis, however, there was a net deposition of P to the bottom sediments. Empirical models were used to describe how the near-surface water quality of Green Lake would be expected to respond to changes in external P loading. Predictions from the models showed a relatively linear response between P loading and TP and chlorophyll-a (Chl-a) concentrations in the lake, with the changes in TP and Chl-a concentrations being less on a percentage basis (50–60 percent for TP and 30–70 percent for Chl-a) than the changes in P loading. Mean summer water clarity, quantified by Secchi disk depths, had a greater response to decreases in P loading than to increases in P loading. Based on these relations, external P loading to the lake would need to be decreased from 8,980 kg/yr to about 5,460 kg/yr for the geometric mean June–September TP concentration in the east side of the lake, with higher TP concentrations than in the west side, to reach the WDNR criterion of 0.015 mg/L. This reduction of 3,520 kg/yr is equivalent to a 46-percent reduction in the potentially controllable external P sources (all external sources except for precipitation, atmospheric deposition, and waterfowl) from those measured during water years 2014–18. The total external P loading would need to decrease to 7,680 kg/yr (a 17-percent reduction in potentially controllable external P sources) for near-surface June–September TP concentrations in the west side of the lake to reach 0.015 mg/L. Total external P loading would need to decrease to 3,870–5,320 kg/yr for the lake to be classified as oligotrophic, with a near-surface June–September TP concentration of 0.012 mg/L. Results from the hydrodynamic water-quality model GLM–AED (General Lake Model coupled to the Aquatic Ecodynamics modeling library) indicated that MOMs are driven by external P loading and internal sediment recycling that lead to high TP concentrations during spring and early summer, which in turn lead to high phytoplankton production, high metabolism and respiration, and ultimately DO consumption in the upper, warmer areas of the metalimnion. GLM–AED results indicated that settling of organic material during summer might be slowed by the colder, denser, and more viscous water in the metalimnion and thus increase DO consumption. Based on empirical evidence from a comparison of MOMs with various meteorological, hydrologic, water quality, and in-lake physical factors, MOMs were lower during summers, when metalimnetic water temperatures were warmer, near-surface Chl-a and TP concentrations were higher, and Secchi depths were lower. GLM–AED results indicated that the external P load would need to be reduced to about 4,060 kg/yr, a 57-percent reduction from that measured in 2014–18, to eliminate the occurrence of MOMs less than 5 mg/L during more than 75 percent of the years (the target provided by the WDNR). Large reductions in external P loading are expected to have an immediate effect on the near-surface TP concentrations and metalimnetic DO concentrations in Green Lake; however, it may take several years for the full effects of the external-load reduction to be observed because internal sediment recycling is an important source of P for the following spring. 
    more » « less
  3. Unconventional oil and gas residual solid wastes are generally disposed in municipal waste landfills (RCRA Subtitle D), but they contain valuable raw materials such as proppant sands. A novel process for recovering raw materials from hydraulic fracturing residual waste is presented. Specifically, a novel hydroacoustic cavitation system, combined with physical separation devices, can create a distinct stream of highly concentrated sand, and another distinct stream of clay from the residual solid waste by the dispersive energy of cavitation conjoined with ultrasonics, ozone and hydrogen peroxide. This combination cleaned the sand grains, by removing previously aggregated clays and residues from the sand surfaces. When these unit operations were followed by a hydrocyclone and spiral, the solids could be separated by particle size, yielding primarily cleaned sand in one flow stream; clays and fine particles in another; and silts in yet a third stream. Consequently, the separation of particle sizes also affected radium distribution – the sand grains had low radium activities, as lows as 0.207 Bq g −1 (5.6 pCi g −1 ). In contrast, the clays had elevated radium activities, as high as 1.85–3.7 Bq g −1 (50–100 pCi g −1 ) – and much of this radium was affiliated with organics and salts that could be separated from the clays. We propose that the reclaimed sand could be reused as hydraulic fracturing proppant. The separation of sand from silt and clay could reduce the volume and radium masses of wastes that are disposed in landfills. This could represent a significant savings to facilities handling oil and gas waste, as much as $100 000–300 000 per year. Disposing the radium-enriched salts and organics downhole will mitigate radium release to the surface. Additionally, the reclaimed sand could have market value, and this could represent as much as a third of the cost savings. Tests that employed the toxicity characteristic leaching protocol (TCLP) on these separated solids streams determined that this novel treatment diminished the risk of radium mobility for the reclaimed sand, clays or disposed material, rendering them better suited for landfilling. 
    more » « less
  4. Hotspot tracks (chains of seamounts, ridges, and other volcanic structures) provide important records of plate motions, as well as mantle geodynamics, magma flux, and mantle source compositions. The Tristan-Gough-Walvis Ridge (TGW) hotspot track, extending from the active volcanic islands of Tristan da Cunha and Gough through a province of guyots and then along Walvis Ridge to the Etendeka flood basalt province, forms one of the most prominent and complex global hotspot tracks. The TGW hotspot track displays a tight linear age progression in which ages increase from the islands to the flood basalts (covering ~135 My). Unlike Pacific tracks, which are often simple, nearly linear chains of seamounts, the TGW track is alternately a steep-sided narrow ridge, an oceanic plateau, subparallel linear ridges and chains of seamounts (most are flat-topped guyots). The track displays isotopic zonation over the last ~70 My. The zonation appears near the middle of the track just before it splits into two to three chains of ridge- and guyot-type seamounts. Walvis Ridge, forming the older part of the track, is also overprinted with age-progressive late-stage volcanism, which was emplaced ~30–40 My after the initial eruptions and has a distinct isotopic composition. The plan for Expedition 391 was to drill at six sites, three along Walvis Ridge and three in the seamounts of the Guyot Province, to collect igneous rocks to better understand the formation of volcanic edifices, the temporal and geochemical evolution of the hotspot, and the variation in paleolatitudes at which the volcanic edifices formed. After a delay of 18 days to address a shipboard Coronavirus (COVID-19) outbreak, Expedition 391 proceeded to drill at four of the proposed sites: three sites on Walvis Ridge around Valdivia Bank, an ocean plateau within the ridge, and one site on the lower flank of a guyot in the Center track of the Guyot Province, a ridge located between the Tristan subtrack (which extends from the end of Walvis Ridge to the islands of Tristan da Cunha) and the Gough subtrack (which extends from Walvis Ridge to Gough Island). The first hole was drilled at Site U1575, located on a low portion of the northeastern Walvis Ridge just north of Valdivia Bank. At this location, 209.9 m of sediments and 122.4 m of igneous basement were cored. The sediments ranged in age from Late Pleistocene (~0.43–1.24 Ma) to Late Cretaceous (Campanian; 72–78 Ma). The igneous basement comprised 10 submarine lava units consisting of pillow, lobate, sheet, and massive lava flows, the thickest of which was ~21 m. Most lavas are tholeiitic, but some alkalic basalts were recovered. A portion of the igneous succession consists of low-Ti basalts, which are unusual because they appear in the Etendeka flood basalts but have not been previously found on Walvis Ridge. Two holes were drilled at Site U1576 on the west flank of Valdivia Bank. The first of these holes was terminated because a bit jammed shortly after entering the igneous basement. Hole U1576A recovered a remarkable ~380 m thick sedimentary section consisting mostly of chalk covering a nearly complete sequence from Late Pleistocene (~0.43–1.24 Ma) to Late Cretaceous (Campanian; ~79–81.38 Ma). These sediments display short and long cyclic color changes that imply astronomically forced and longer term paleoenvironmental changes. The igneous basement recovered in Hole U1576B yielded 11 submarine lava units (total thickness = ~65 m). The flows range from pillows to massive flows with compositions varying from tholeiitic basalt to basaltic andesite, only the second occurrence of the latter composition recovered from the TGW track thus far. These units are separated by seven sedimentary chalk units that range 0.1–11.6 m in thickness, implying a long-term interplay of sedimentation and lava eruptions. These intercalated sediments revealed Upper Cretaceous (Campanian) ages of ~77–79 Ma for the upper two interbeds and ~79–81.38 Ma for the lower beds. Coring at Site U1577, on the extreme eastern flank of Valdivia Bank, penetrated a 154.8 m thick sedimentary section ranging from the Paleocene (Thanetian; ~58.8 Ma) to Upper Cretaceous (Campanian; ~81.43–83.20 Ma). Igneous basement coring progressed only 39.1 m below the sediment/basalt contact, recovering three massive submarine tholeiitic basalt lava flows that are 4.1, 15.5, and >19.1 m thick, respectively. Paleomagnetic data from Sites U1577 and U1576 indicate that the former volcanic basement formed just before the end of the Cretaceous Normal Superchron and the latter during Chron 33r, shortly afterward. Biostratigraphic and paleomagnetic data suggest that Valdivia Bank becomes younger from east to west. Site U1578, located on a Center track guyot, provided a long and varied igneous section. After coring through 184.3 m of pelagic carbonate sediments mainly consisting of Eocene and Paleocene chalk (~55.64–63.5 Ma), Hole U1578A cored 302.1 m of igneous basement. Basement lavas are largely pillows but are interspersed with sheet and massive flows. Lava compositions are mostly alkalic basalts with some hawaiite. Several intervals contain abundant olivine (some fresh), and some of the pillow stacks consist of basalt with remarkably high Ti content. The igneous sequence is interrupted by 10 sedimentary interbeds consisting of chalk and volcaniclastics and ranging 0.46–10.19 m in thickness. Investigations of toothpick samples from the intercalated sediments were examined, each revealing the same age range of ~63.5–64.81 Ma (lower Paleocene; Danian). Paleomagnetic data display a change in basement magnetic polarity ~100 m above the base of the hole. Combining magnetic stratigraphy with biostratigraphic data, the igneous section is inferred to span >1 My. Nearly 7 months after Expedition 391, JOIDES Resolution transited from Cape Town to the north Atlantic. During this transit (Expedition 397T), 7.9 days of ship time were used to drill two holes (U1584A and U1585A) at sites on the Gough and Tristan tracks that had been omitted because of COVID-19–related time loss on the earlier cruise. For both, coring was begun only a short distance above the igneous basement to save time. The 75.2 m thick section drilled in Hole U1584A contains two sedimentary units: clay-rich carbonate sediments overlie a pumice-dominated volcaniclastic deposit containing basalt fragments. Because the goal was to core basalt and the base of the volcaniclastic deposit was not imaged in the seismic profile, the hole was terminated early to save operation time for the next site. In Hole U1585A, coring penetrated a 273.5 m thick sediment section overlying an 81.2 m thick pile of massive basalt flows. The sediment section is divided into four units: The uppermost unit consists of nannofossil chalk; The two intermediate units contain alternating chalk and volcaniclastic sediments containing several breccia units; and The lowermost unit consists of volcanic breccia containing juvenile blocks, bombs, and accretionary lapilli. This thick sedimentary section documents a transition from shallow-water volcanism to open-ocean sedimentation as the seamount subsided. The thick underlying basalt section is made up of four sparsely to highly phyric massive flows, the thickest of which is >43 m thick. Samples of these units are mostly basalt with a few trachybasalts and one trachyandesite. Although the igneous penetration was less than planned, coring during Expeditions 391 and 397T obtained samples that clearly will lead to an improved understanding of the evolution of the TGW hotspot and its track. Reasonable recovery of fresh basalt in some holes provides ample samples for geochemical, geochronologic, and paleomagnetic studies. Good recovery of Late Cretaceous and early Cenozoic chalk successions provides samples for paleoenvironmental study. 
    more » « less
  5. Abstract

    Ice streams that flow into Ross Ice Shelf are underlain by water-saturated sediments, a dynamic hydrological system, and subglacial lakes that intermittently discharge water downstream across grounding zones of West Antarctic Ice Sheet (WAIS). A 2.06 m composite sediment profile was recently recovered from Mercer Subglacial Lake, a 15 m deep water cavity beneath a 1087 m thick portion of the Mercer Ice Stream. We examined microbial abundances, used 16S rRNA gene amplicon sequencing to assess community structures, and characterized extracellular polymeric substances (EPS) associated with distinct lithologic units in the sediments. Bacterial and archaeal communities in the surficial sediments are more abundant and diverse, with significantly different compositions from those found deeper in the sediment column. The most abundant taxa are related to chemolithoautotrophs capable of oxidizing reduced nitrogen, sulfur, and iron compounds with oxygen, nitrate, or iron. Concentrations of dissolved methane and total organic carbon together with water content in the sediments are the strongest predictors of taxon and community composition. δ¹³C values for EPS (−25 to −30‰) are consistent with the primary source of carbon for biosynthesis originating from legacy marine organic matter. Comparison of communities to those in lake sediments under an adjacent ice stream (Whillans Subglacial Lake) and near its grounding zone provide seminal evidence for a subglacial metacommunity that is biogeochemically and evolutionarily linked through ice sheet dynamics and the transport of microbes, water, and sediments beneath WAIS.

     
    more » « less