Rate of K–Mg–Fe sequestration by glauconite in marine settings. Reverse weathering reactions to produce authigenic clay minerals in the ocean can significantly impact the ratio of element (re)cycling vs element sequestration in marine sediments. Such element sequestration reactions are commonly slower (105 to 106 yr) compared to clay mineral formation by continental weathering (103 to 105 yr) or hydrothermal-driven fluid-rock interactions at mid-oceanic ridges (102 to 104 yr). Nevertheless, previous studies have identified ‘hot spot’ areas that favor clay mineral precipitation in marine sediments, such as mangrove forests, deltas and estuaries (Michalopoulos and Aller, 1995; Cuadros et al., 2016), low- to high-temperature hydrothermal sites (Huang et al., 2018; Shalev et al., 2019; Voigt et al., 2020; Laureijs et al., 2021) and shallow-water shelfal settings characterized by mega-condensed sedimentation (Baldermann et al., 2017; Bansal et al., 2020b). In such surroundings, clay retrograde reactions can behave as important controls on element sequestration and marine elemental output fluxes at or close to the sediment-seawater interface.
Herein, we focus on constraining the sequestration rates for K and Mg, as these elements are taken up directly from seawater or seawater-derived pore fluids by glauconite, and on Fe, because this element is an essential micro-nutrient in the ocean and plays a key role in the sedimentary Fe redox cycle, in particular when the rate of pyrite formation is low (Raiswell, 2006; Baldermann et al., 2013; López-Quirós et al., 2019). We do not consider Na and Ca, as they are barely incorporated in glauconite (Table S1), and also exclude Al and Si from further consideration, as these elements are virtually ‘absent’ in modern ocean water. This is attributed to the low solubility of most Al- and Si-bearing solids and the strong concentration limitations arising from the activity of silicifying organisms, such as radiolarians, diatoms and siliceous sponges, in particular in the surface oceans (Hoffmann et al., 2008; Kranzler et al., 2019).
Assuming (i) a sedimentation rate of 130 cm·Myr− 1 for the M. dixoni Zone (Wilmsen, 2007) and (ii) a density of 2.7 g·cm−³ for the glauconitized strata (Logvinenko, 1982), and considering (iii) the K–Mg–Fe contents of the glauconites (Table S1) and (iv) the abundance of glauconite across the studied sedimentary sequence (Fig. 3a,b), element-specific sequestration rates associated with shallow-water glauconite formation can be calculated (Fig. 3e and Table S3). The sequestration rates range from 0.7 to 42 mmol K·cm−²·kyr− 1, 0.4 to 23 mmol Mg·cm−²·kyr− 1 and 1.1 to 65 mmol Fe·cm−²·kyr− 1, reflecting the extremely low to very high abundances (1 vs 70 wt.%) of glauconite in the profile.
The ‘low rate’ of element sequestration may reflect reverse weathering processes taking place in depositional systems that are characterized by high sedimentation rates and a short residence time at the sediment-seawater interface. Authigenic clay precipitation is inefficient under such conditions (Amorosi, 2012; Chattoraj et al., 2016) so that element sequestration is reduced. The ‘high rate’ of element sequestration may represent transgressive systems that are characterized by low to close to zero overall sedimentation and a prolonged residence time at the sediment-seawater interface, which favor glauconite formation (Amorosi, 2012; Banerjee et al., 2020).
Considering an average glauconite content of 5–10 wt.% for the ~ 2.8 m thick glauconitized interval at Langenstein, which is representative of Phanerozoic glauconite deposits (7 ± 4 wt.%; Banerjee et al., 2016), the sequestration rate for glauconite formation in the palaeo-shelf environment of the Upper Cretaceous averages 7 ± 5 mmol K·cm−²·kyr− 1, 3 ± 1 mmol Mg·cm−²·kyr− 1 and 14 ± 13 mmol Fe·cm−²·kyr− 1. For comparison, the sequestration rates for glauconite-smectite and glauconite forming in the modern deep-sea sediments of the Ivory Coast (Ghana Marginal Ridge) were determined as 0.017 mmol K·cm−²·kyr− 1, 0.024 mmol Mg·cm−²·kyr− 1 and 0.080 mmol Fe·cm−²·kyr− 1 (partly recalculated from Baldermann et al., 2015), which is, on average, ~ 100-400-times lower than shelfal sequestration rates. However, the glauconite content (2.5 wt.%, on average), the K2O concentration (2.9 wt.%, on average, reflecting the ‘nascent’ stage) and the rate of glauconite formation (~ 5-times lower) are substantially lower in deep-water settings than in the shelf regions. The lower rate of glauconite formation in deep-water settings is due to the low temperature (~ 5 vs ~ 25°C), the reduced supply and sedimentary reflux of Al3+ ions, silicic acid and organic matter, the latter controlling local redox restrictions (semi-confined micromilieu vs redoxcline) that predetermine the availability and the speciation of Fe (Fe2+ vs Fe3+), which is the rate-determining factor for glauconite formation (Baldermann et al., 2013; López-Quirós et al., 2019).
Effect of glauconite formation on global marine K–Mg–Fe palaeo-fluxes. The source-sink relations of chemical elements in the modern ocean are relatively well constrained (Raiswell, 2006; Henderson et al., 2007; Arvidson et al., 2013; Dale et al., 2015; Jeandel and Oelkers, 2015; Cuadros et al., 2016; Huang et al., 2018; Shalev et al., 2019; Somes et al., 2021). Model estimates indicate that the marine K, Mg and Fe cycle is mainly controlled by flux variations of riverine (and groundwater), atmospheric, glacial and hydrothermal inputs vs silicate and carbonate mineral deposition on the ocean floor. However, the role of authigenic clay formation on the fluxes in the marine sedimentary K–Mg–Fe budgets of the modern and ancient oceans remains enigmatic.
Although the K–Mg–Fe sequestration rates reported for glauconite formation at Langenstein (shallow-water; this study) and the Ivory Coast (deep-water; Baldermann et al., 2015) may not be directly transferrable to all other marine settings that accumulated glauconite through time and space, (i) the mode of glauconite formation (Fe-smectite-to-glauconite reaction), (ii) the micro-environment (faecal pellets and foraminifera), (iii) the timing (105-106 yr), (iv) the composition (Fe-rich, highly evolved vs nascent), (v) the abundance (5–10 wt.% vs 2–3 wt.%) and (vi) the depositional environment (warm shallow shelf vs cool deep-sea; low sedimentation rate) at the two localities are representative of the range expected for many modern and past glauconite-forming environments (Odin and Matter, 1981; Baldermann et al., 2015; Banerjee et al., 2016; Baldermann et al., 2017; Bansal et al., 2017; López-Quirós et al., 2019; Banerjee et al., 2020; Bansal et al., 2020a,b; Choudhury et al., 2021). Furthermore, the Langenstein glauconites share many similarities with other Mesozoic to Cenozoic glauconite deposits, such as a similar abundance (7 ± 4 wt.% in the marine rock record from the Triassic to the Holocene time; Fig. 4a; Banerjee et al., 2016), a comparable chemical composition (7 ± 1 wt.% K2O, 4 ± 1 wt.% MgO and 23 ± 2 wt.% TFe; Fig. 4b; Banerjee et al., 2016), and a low overall sedimentation rate (1 to 5 m·Myr− 1; Fig. 4c; Chattoraj et al., 2016). If we assume that the sedimentation rate (1.3 m·Myr− 1) and the sediment density (2.7 g·cm−³) on the palaeo-shelf at Langenstein are representative of the global shelf through geological time, we can compute K–Mg–Fe palaeo-sequestration rates for shallow-water glauconite formation for Mesozoic and Cenozoic times (Fig. 4d-f; Table S4).
It is evident that glauconite formation significantly contributed to K–Mg–Fe sequestration in shallow marine sediments throughout Earth`s history, averaging 4 ± 3 mmol K·cm−²·kyr− 1, 4 ± 2 mmol Mg·cm−²·kyr− 1 and 10 ± 6 mmol Fe·cm−²·kyr− 1, respectively. We note that K–Mg–Fe sequestration by glauconite formation also happened in the ‘older’ sediments of the Archaean, Proterozoic and early Cambrian, but this elemental uptake can be barely quantified given that sedimentary archives of this time are scarce and that most of the ‘old’ glauconites are at least partly altered to illite or chlorite minerals (Bansal et al., 2020b; Rafiei et al., 2020). K–Mg–Fe sequestration by green clay authigenesis was of minor importance in the Late Ordovician, Early Silurian and Late Devonian, which corresponded to major glacial events, where glauconite formation is inefficient (Banerjee et al., 2016). On the contrary, high glauconite abundances and related high K–Mg–Fe sequestration rates are evident, for example, in the Neogene (glauconite sands from the Chatham Rise, Southwest Pacific), at the Paleogene-Eocene Transition (glauconite sands from the continental margins of the northern hemisphere; Upper and Lower Greensands of England) and in the Cretaceous (New Jersey, Maryland and Delaware Greensands; greensand giants from the Duwi group, Egypt; Bakchar glauconite deposit, Western Siberia) (cf. Figure 4d-f). These periods record glauconite deposits of huge economic or geological value (Banerjee et al., 2016; Banerjee et al., 2020).
Using average elemental sequestration rates per geological period (Fig. 5a) and corresponding occurrences of glauconite on the shelf (Fig. 4a), as well as calculated low and high estimates of the shallow ocean areas over time (Fig. 5b; Cao et al., 2017; Scotese and Wright, 2018), we can compute K–Mg–Fe palaeo-fluxes (Tmol·yr− 1) associated with green clay authigenesis that progressed on the world`s shelf area (defined here as 0-200 m water depth) over time (Fig. 5c-e; Table S5). Based on comparison with global K–Mg–Fe fluxes of the modern and past ocean, we propose that the obtained ‘high’ K–Mg–Fe palaeo-fluxes are overestimated (i.e., the shelf areas reported by Cao et al. (2017) are overestimated) and that the ‘low’ K–Mg–Fe palaeo-fluxes (i.e., the shelf areas reported by Scotese and Wright (2018)) better portray the average elemental burial related to green clay authigenesis per geological period, which averages 0.09 ± 0.12 Tmol K·yr− 1, 0.08 ± 0.11 Tmol Mg·yr− 1 and 0.20 ± 0.28 Tmol Fe·yr− 1 during the Triassic to Holocene. The ratio of the K–Mg–Fe palaeo-fluxes associated with glauconite formation in the past vs modern ocean (Fig. 5d) indicates further that the elemental fluxes were much higher from the Jurassic to the Oligocene compared to the modern ocean, averaging a factor of 2.4 ± 3.4, which we attribute to the warm and sea level highstand ‘greenhouse’ conditions prior to the Eocene-Oligocene transition. The lower K–Mg–Fe palaeo-fluxes ever since the Oligocene are caused by the decrease of the shallow-water shelf areas and seawater temperature with the onset of the first Southern Hemisphere glaciation (~ 34 Myr ago) and then Northern Hemisphere glaciation (~ 5 Myr ago), where glauconite formation is reduced (Banerjee et al., 2016). Although these estimates have a relatively high uncertainty and need to be better constrained in future work, it is evident that green clay authigenesis greatly affected the global marine K, Mg and in particular Fe cycle throughout Earth`s history.
K–Mg–Fe sequestration by green clay authigenesis in modern marine settings. To the best of our knowledge, elemental output fluxes attributed to widespread glauconite formation taking place at the shallow shelf and in the deep-sea (defined here as > 2000 m water depth) of the modern oceans have not been determined yet and are not fully accounted for in earth system models (Arvidson et al., 2013; Baldermann et al., 2015; Cuadros et al., 2016; Haley et al., 2017; Rahman et al., 2017; Isson and Planavsky, 2018; Abbott et al., 2019). We recognize that some of the glauconite deposits of the Quaternary and Holocene are of para-autochthonous or detrital origin, representing reworked glauconites of the Neogene or even older ages (Odin and Matter, 1981; Baldermann et al., 2013; Banerjee et al., 2016). As a first-order approximation to calculate the present-day global K–Mg–Fe output fluxes attributed to green clay authigenesis, we use (i) published glauconite contents in the shallow and deep-marine sediments of the Holocene (5.6 wt.% vs 2.5 wt.%; Banerjee et al., 2016; Logvinenko, 1982), (ii) the total areas of the modern shelf and deep-sea regions (27.12·1012 m² vs 302.5·1012 m²; Dale et al., 2015) and (iii) K–Mg–Fe sequestration rates by Holocene shallow-water and deep-water glauconite formation (3.0–4.0 mmol K·cm−²·kyr− 1, 2.4–4.1 mmol Mg·cm−²·kyr− 1 and 7.4–8.8 mmol Fe·cm−²·kyr− 1, Fig. 5d-f, vs 0.017 mmol K·cm−²·kyr− 1, 0.024 mmol Mg·cm−²·kyr− 1 and 0.080 mmol Fe·cm−²·kyr− 1; Baldermann et al., 2015). Upscaling predicts global K–Mg–Fe output fluxes associated with glauconite formation of ~ 0.05–0.06 Tmol K·yr− 1, ~ 0.04–0.06 Tmol Mg·yr− 1 and ~ 0.11–0.14 Tmol Fe·yr− 1 at the shallow shelf, and of ~ 0.001 Tmol K·yr− 1, ~ 0.002 Tmol Mg·yr− 1 and ~ 0.006 Tmol Fe·yr− 1 for the deep-sea glauconites (Fig. 6). The calculated K–Mg–Fe fluxes have a relatively high uncertainty, but the estimates are well within global marine K–Mg–Fe fluxes published in the literature (Dale et al., 2015; Sun et al., 2016; Huang et al., 2018; Shalev et al., 2019). However, we assume the herein reported K–Mg–Fe fluxes are conservative estimates, if considering that the abundance of glauconite in shallow marine sediments was much higher in the majority of the Phanerozoic due to sea level lowstands and higher seawater temperatures under warm ‘greenhouse’ conditions vs modern ‘icehouse’ conditions that significantly reduce the burial fluxes linked to green clay authigenesis (Banerjee et al., 2020).
K sequestration by deep-water glauconite formation is barely significant and accounts for ~ 0.1 % removal of the total dissolved riverine K influx (~ 1.51 Tmol·yr− 1; Garrels and Mackenzie, 1971) and of the K supplied to the ocean by hydrothermal alteration of the modern oceanic crust (~ 1.51 Tmol·yr− 1; Berner and Berner, 2012). Contrary, shallow-water glauconite formation can play an important role in the global K cycle (Fig. 6), as it fixes ~ 3–4 % of the total oceanic K inventory that is sourced from riverine and hydrothermal fluxes, or ~ 3 % of the K that is removed from the ocean through low-temperature alteration of ocean-floor basalts (~ 1.99 Tmol·yr− 1; Sayles, 1979; Sun et al., 2016). Hence, K sequestration by glauconite formation at the shelf is at the same order of magnitude as K sequestration by Fe-illite formation taking place in the mangrove forests worldwide (~ 0.02–0.08 Tmol·yr− 1; Cuadros et al., 2016). An exception to this amount of K sequestration occurs when K is available locally from e.g., a K-feldspar substrate altering to glauconite, as reported in Proterozoic oceans (Banerjee et al., 2015).
The same conclusions can be drawn for the marine Mg cycle (Fig. 6): Glauconite formation at the shelf consumes ~ 1 % of the terrestrial Mg flux (~ 5.51 Tmol·yr− 1) that is brought to the ocean via continental weathering of Mg-bearing carbonates and silicates (Fig. 6) and constitutes ~ 1–2 % to the marine Mg sink that is associated with low- and high-temperature alteration of the oceanic crust (~ 2.71 Tmol·yr− 1; Huang et al., 2018). Further, Mg sequestration by glauconite is equivalent to ~ 1–4 % of the estimated Mg sink by enigmatic (yet hidden) dolomite deposits (~ 1.48–2.88 Tmol·yr− 1; Shalev et al., 2019) or ~ 19–29 % of the Mg consumption by authigenic clays forming in the Amazon deltaic sediments (~ 0.21 Tmol·yr− 1; Sun et al., 2016). Contrary, Mg sequestration associated with deep-water glauconite formation accounts for only ~ 0.3 % of the low-temperature alteration flux (~ 0.66 Tmol·yr− 1; Higgins and Schrag, 2012; Sun et al., 2016; Shalev et al., 2019).
In addition, glauconite acts as an important sink for Fe (Fig. 6), with shallow-water glauconite formation accounting for up to ~ 13–29 % removal of the dissolved and particulate riverine flux of highly reactive Fe to the ocean (~ 0.48–0.86 Tmol·yr− 1 and ~ 0.63 Tmol·yr− 1; Raiswell, 2006). Although Fe uptake by deep-water glauconite formation is less significant (Fig. 6), it is still equivalent to ~ 2 % removal by the total hydrothermal alteration flux (~ 0.25 Tmol·yr− 1) or ~ 3 % and ~ 30 % removal by the glacial flux (~ 0.20 Tmol·yr− 1) and the atmospheric dust flux (~ 0.02 Tmol·yr− 1; Raiswell, 2006). Even though oxidation and scavenging processes are the first-order controls for the benthic Fe fluxes in the ocean (Dale et al., 2015), we conclude that glauconite formation taking place in shallow and deep marine settings is an important but currently overlooked mechanism and elemental sink. These glauconite formation settings can, at least locally, affect the pore water inventory of Fe2+, K+ and Mg2+ ions that might otherwise be available for biogeochemical use for the ocean and the marine sediments.
We conclude that fast retrograde clay mineral reactions, which occur to wide extent on the ocean floor, are of great significance in the marine K, Mg and in particular Fe cycle and have to be considered in earth system models of the present and past marine element cycles. The elemental burial fluxes attributed to green clay authigenesis were significantly higher under the sea level highstand and ‘greenhouse’ conditions in the majority of the Phanerozoic compared to the modern sea level lowstand and icehouse’ conditions, which suppress glauconite formation. It is now up to future studies (i) to estimate how K–Mg–Fe sequestration through glauconite formation impacts the isotopic composition of the ocean, the pore water reservoir and the modern marine (deep-sea) sediments, and (ii) to assess the impact of climate change through time on the elemental burial fluxes attributed with reverse silicate weathering and green clay authigenesis.