skip to main content


This content will become publicly available on May 9, 2025

Title: Tracing the Environmental Effects of Mineral Fertilizer Application with Trace Elements and Strontium Isotope Variations
Fertilizer utilization is critical for global food security. This study examines the occurrence of trace elements (TEs) and Sr isotope (87Sr/86Sr) variations in phosphate rocks and mineral fertilizers from a sample collection representative of major phosphate producing countries. We show high concentrations of several TEs in phosphate rocks (n=76) and their selective enrichment in phosphate fertilizers (n=40) of specific origin. Consistent with the concentrations in parent phosphate rocks, phosphate fertilizers from the U.S. and Middle East have substantially higher concentrations of U, Cd, Cr, V, and Mo than in fertilizers from China and India. Yet, fertilizers from China and India generally have higher concentrations of As. The 87Sr/86Sr in phosphate fertilizers directly mimic the composition of their source phosphate rocks, with distinctive higher ratios in fertilizers from China and India (0.70955–0.71939) relative to phosphate fertilizers from U.S. and Middle East (0.70748–0.70888). Potash fertilizers have lower Sr and TEs and higher 87Sr/86Sr (0.72017–0.79016), causing higher 87Sr/86Sr in mixed NPK-fertilizers. Selective extraction (Mehlich III) of soils from an experimental agricultural site shows relative enrichment of potentially plant-available P, Sr, and TEs in topsoil, which is associated with Sr isotope variation towards the 87Sr/86Sr of the local utilized phosphate fertilizer.  more » « less
Award ID(s):
2305946 2305947
NSF-PAR ID:
10507951
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
American Chemical Society Publications
Date Published:
Journal Name:
Environmental Science & Technology Letters
ISSN:
2328-8930
Subject(s) / Keyword(s):
["Phosphate, Potash, Soil and Water Pollution, Agriculture"]
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Phosphate rock bears both geologically and environmentally significant information. Rare earth elements and yttrium (i.e., REY) characteristics have been commonly used for reconstructing the redox conditions of depositional environments of and the effects of post-depositional diagenetic alteration on phosphate rock. In addition, phosphate rock is typically enriched in a range of trace elements such as uranium (U) and cadium (Cd) that can be dispersed as contaminants into the environment with phosphate mining and phosphate fertilizer application. Here we report the lead (Pb) isotope compositions combined with Pb and REY concentrations of both global sedimentary and igneous phosphate rocks, aiming to evalute the geological origin of phosphate rocks over time and the potential of using them for environmental tracing. Phopshate rocks samples analyzed in this study were sourced from major economic phosphate deposits in the world, including China, Southern Tethys (e.g., Morocco, Tunisia, Israel), the U.S., India, South Africa and Russia. Our results show a wide range of 208Pb/204Pb (35.70 to 60.58), 207Pb/204Pb (15.20 to 18.25), and 206Pb/204Pb (16.369 to 71.806) ratios in phosphate rocks, with sedimentary phosphate rocks being significantly more radiogenic than igneous rocks. The majority of the sedimentray phosphate rocks show a notable isotopic overprinting by non-radiogenic terrestrial Pb, except for those from Israel and Morocco that have the most radiogenic Pb isotope compositions. Correspondingly, phosphate rocks with more radiogenic Pb isotope ratios show relatively pristine seawater REY features, likely suggesting their preservation of the original oxic seawater conditions and/or minimal diagenetic alteration. In contrast, phosphate rocks with less radiogenic Pb isotope compositions show REY indications for more anoxic seawater redox conditions and/or greater diagenetic alteration. We further evaluate the potential utility of Pb isotopes for tracing the associated contamination with phosphate rock mining and fertilizer application in the environment. In most cases, the radiogenic Pb isotope composition of phosphate rocks and corresponding P-fertilizers is distinctive from both natural crustal Pb and major anthropogenic Pb sources (e.g., Pb ore deposits and pesticides), which provides a great advantage for applying Pb isotopes as environmental tracers for metal(loid) contamination from phosphate sources. The combination of Pb isotope ratios and REY proxies could further constrain the Pb source discrimination. Overall, this study provides new Pb isotopic and REY geochemical data on global phosphate rocks and fertilizers, which lays the groundwork for future regional and local studies on both their geological and environmental implications. 
    more » « less
  2. null (Ed.)
    Here we report the strontium isotope ratios (87Sr/86Sr) from pore fluids collected during International Ocean Discovery Program (IODP) Expedition 355. Ratios from Sites U1456 (N = 21) and U1457 (N = 20) are reported. Pore fluid 87Sr/86Sr is a useful tool to establish fluid-rock reactions, sources of Sr, and fluid mixing. The measured 87Sr/86Sr of the pore fluid has significant variations at both sites, and three distinct zones are identified. At Site U1456, 87Sr/86Sr starts at values similar to that of modern seawater (~0.7092) from near the seafloor down to ~100 meters below seafloor (mbsf). Over this interval, Sr concentration increases, whereas Ca decreases (Zone 1). Below 100 mbsf, 87Sr/86Sr values increase to a max of ~0.7100 at ~224 mbsf, and Sr concentrations decrease (Zone 2). Isotopic values then gradually decrease to ~0.7085, with an increase in Sr concentrations (Zone 3). Site U1457 displays the same general trends in pore fluid 87Sr/86Sr composition; however, there are distinct differences. First, Zone 1 occurs over a shorter interval (~0–54 mbsf) due to a lower sedimentation rate, and it also has lower Sr concentrations compared to Site U1456. Additionally, 87Sr/86Sr reaches a higher peak value in Zone 2 than at Site U1456. Finally, in Zone 3 the maximum Sr concentration reached is significantly lower than that at Site U1456. 
    more » « less
  3. null (Ed.)
    Oceanic island basalts are targeted for geochemical study because they provide a direct window into mantle composition and a wealth of information on the dynamics and timescales associated with Earth mixing. Previous studies mainly focused on the shield volcanic stage of oceanic islands and the more fusible, enriched mantle components that are easily distinguished in those basalts. Mantle depleted compositions are typically more difficult to resolve unless large amounts of this material participated in mantle melting (e.g., mid-ocean ridges), or unique processes allow for their compositions to be erupted undiluted, such as very small degrees of melting of a source with minimal fusible enriched components (e.g., rejuvenated basalts) or as xenoliths (e.g., abyssal peridotites). Mantle depleted components, defined here as material with low time-integrated Rb/Sr (low 87Sr/86Sr) and high time-integrated Sm/Nd and Lu/Hf ratios (high 143Nd/144Nd and 176Hf/177Hf) relative to primitive mantle, derive from a potentially very large volume reservoir (up to 80% of the mantle), and therefore need adequate characterization in order estimate the composition of the Earth and mantle-derived melts. This review focuses on mantle depleted compositions in oceanic island basalts using the Hawaiian-Emperor chain as a case study. The Hawaiian-Emperor chain is the ∼6000 km long geological record of the deeply sourced Hawaiian mantle plume, active for>81 Myr. Hawaiian volcanism evolves through four volcanic stages as a volcano traverses the Hawaiian plume: alkalic preshield, tholeiitic shield (80–90% volcano volume), alkalic postshield (∼1%), and silica undersaturated rejuvenated (< 0.1%). We report Pb-Sr-Nd-Hf isotope compositions and trace element concentrations of three rejuvenated Northwest Hawaiian Ridge basalts and compare them to an exhaustive compiled dataset of basalts from the Hawaiian Islands to the Emperor Seamounts. The Northwest Hawaiian Ridge (NWHR) includes 51 volcanoes spanning ∼42 m.y. between the bend in the Hawaiian-Emperor chain and the Hawaiian Islands where there is no high-precision isotopic data published on the rejuvenated-stage over ∼47% of the chain. NWHR and Hawaiian Island rejuvenated basalts are geochemically similar, indicating a consistent source for rejuvenated volcanism over ∼12.5 million years. In contrast, shield-stage basalts from the oldest Emperor Seamounts are more depleted in isotopic composition (i.e., higher 176Hf/177Hf, and 143Nd/144Nd with lower 87Sr/86Sr and 208Pb*/206Pb*) and trace element concentrations (i.e., much lower concentrations of highly incompatible elements) than all other depleted Hawaiian basalts younger than the bend, including NWHR rejuvenated basalts. The strongly depleted source for the oldest Emperor Seamounts (> 70 Ma) was likely related to interaction with the Kula-Pacific-Izanagi mid-ocean ridge spreading system active near the Hawaiian plume in the Late Cretaceous. In contrast, the incompatible trace element ratios of NWHR rejuvenated basalts require a distinct source in the Hawaiian mantle plume that was imprinted by ancient (> 1 Ga) partial melting, likely ancient recycled oceanic lithosphere. This review of the geochemistry of Hawaiian depleted components documents the need for the sampling of multiple distinctive depleted compositions, each preferentially melted during specific periods of Hawaiian plume activity. This suggests that the composition of depleted components can evolve during the lifetime of the mantle plume, as observed for enriched components in the Hawaiian mantle plume. Changes in the composition of depleted components are dominantly controlled by the upper mantle tectonic configurations at the time of eruption (i.e., proximity to a mid-ocean ridge), as this effect overwhelms the signal imparted by potentially sampling different lower mantle components through time. 
    more » « less
  4. The Sierra Nevada Batholith (SNB) records copious Mesozoic magmatism and is an important touchstone for understanding crustal growth at continental convergent margins. Recent research in the SNB has focused on defining magmatic cyclicity and arc “flare ups” based on the ages, magma production rates, and radiogenic isotope heterogeneities of the plutonic and volcanic rocks found throughout the batholith. Two main intervals at ca. 170–148 Ma and ca. 125–85 Ma delivered >95% of the magmas in the exposed plutonic bulk in the SNB and suggest elevated emplacement rates and hotter-than-usual magmas, though the Cretaceous is by far the most productive era and the most promising for understanding the factors modulating magmatic flux. The mid-Cretaceous of the Sierra (ca. 105–98 Ma) saw the appearance of conspicuous, high-silica (>65 wt.% SiO2; average ~71%) granitic plutons of similar chemical nature that span a large geographic area, breaking the well-established west-to-east “younging” trend found in the more common rocks of intermediate compositions. This study focuses on thirteen of these high-silica granites: the Bullfrog, Independence, McGann, Rawson Creek, and Spook Plutons of the eastern Sierra; and the Shaver Intrusive Suite, Grant Grove, Case Mountain, Coyote Pass, Dennison Peak, and Frys Point Plutons of the western/central Sierra. Whole rock geochemistry, zircon trace elements, and radiogenic isotope ratios (Sr and Nd) in these high-silica granites show some transitional patterns with other contemporaneous and geographically related plutons of intermediate compositions, suggesting fractionation trajectories; however, some distinct dissimilarities are observed, including: 1) elevated, but highly varied initial 87Sr/86Sr ratios, 2) elevated fluorine in granites, and 3) hotter apparent zircon saturation conditions. These geochemical data, hotter conditions, and higher flux suggest that mantle conditions favored more crustal melting and crustal source input than at any other time in the Cretaceous. We conclude that the granitic outburst of the mid-Cretaceous was a flare up like no other. 
    more » « less
  5. Jacobson, A. (Ed.)
    Ocean anoxic events (OAE) are characterized by increased organic content of marine sediment on a global scale with accompanying positive excursions in sedimentary organic and inorganic d 13C values. To sustain the increased C exports and burial required to explain the C isotope excursion, increased supplies of nutrients to the oceans are often invoked during ocean anoxic events. The potential source of nutrients in these events is investigated in this study for Oceanic Anoxic Event 2, which spans the Cenomanian-Turonian boundary. Massive eruptions of one or more Large Igneous Provinces (LIPs) are the proposed trigger for OAE 2. The global warming associated with volcanogenic loading of carbon dioxide to the atmosphere has been associated with increased continental weathering rates during OAE 2, and by extension, enhanced nutrient supplies to the oceans. Seawater interactions with hot basalts at LIP eruption sites can further deliver ferrous iron and other reduced metals to seawater that can stimulate increased productivity in surface waters and increased oxygen demand in deep waters. The relative importance of continental and submarine weathering drivers of expanding ocean anoxia during OAE 2 are difficult to disentangle. In this paper, a box model of the marine Sr cycle is used to constrain the timing and relative magnitudes of changes in the continental weathering and hydrothermal Sr fluxes to the oceans during OAE 2 using a new high-resolution record of seawater 87Sr/86Sr ratios preserved in a marl-limestone succession from the Iona-1 core collected from the Eagle Ford Formation in Texas. The results show that seawater 87Sr/86Sr ratios change synchronously with Os isotope evidence for the onset of massive LIP volcanism 60 kyr before the positive C isotope excursion that traditionally marks the onset of OAE 2. The higher temporal resolution of the seawater Sr isotope record presented in this study warrants a detailed quantitative analysis of the changes in continental weathering and hydrothermal Sr inputs to the oceans during OAE 2. Using an ocean Sr box model, it is found that increasing the continental weathering Sr flux by  1.8-times captures the change in seawater 87Sr/86Sr recorded in the Iona-1 core. The increase in the continental weathering flux is smaller than the threefold increase estimated by studies of seawater Ca isotope changes during OAE 2, suggesting that hydrothermal forcing may have played a larger role in the development of ocean anoxic events than previously considered. 
    more » « less