Original Research ARTICLE
Hydrothermal and Cold Spring Water and Primary Productivity Effects on Magnesium Isotopes: Lake Myvatn, Iceland
- 1LOGIC (London Geochemistry and Isotope Centre), Institute of Earth and Planetary Sciences, University College London and Birkbeck, University of London, London, United Kingdom
- 2Department of Earth Sciences, Durham University, Durham, United Kingdom
- 3Earth and Life Institute, Université catholique de Louvain, Louvain-la-Neuve, Belgium
- 4Institute of Earth Sciences, University of Iceland, Reykjavik, Iceland
- 5Department of Earth Sciences, Oxford University, Oxford, United Kingdom
- 6Myvatn Research Station, Skútustaðahreppur, Iceland
- 7Institute of Life- and Environmental Sciences, University of Iceland, Reykjavik, Iceland
Lake Myvatn, Iceland, is one of the most biologically productive lakes in the northern hemisphere, despite seasonal ice cover. Hydrothermal and groundwater springs make up the dominant source to this lake, and we investigate their Mg isotope ratio to assess the effect of mid-ocean ridge hydrothermal springs, which are the primary modern sink of seawater magnesium. We also examine a time series in the only outflow from this lake, the Laxa River, to assess the effects of seasonal primary productivity on Mg isotopes. In the hydrothermal waters, there is a clear distinction between cold waters (largely unfractionated from primary basalt) and relatively hot waters, which exhibit over 1‰ fractionation, with consequences for the oceanic mass balance if the hydrothermal removal of Mg is not fully quantitative. The outflow Mg isotopes are similar to basalts (δ26Mg = −0.2 to −0.3) during winter but reach a peak of ∼0‰ in August. This fractionation corresponds to calcite precipitation during summer in Lake Myvatn, preferentially taking up light Mg isotopes and driving the residual waters isotopically heavy as observed, meaning that overall the lake is a CO2 sink.
Chemical weathering of continental silicate rocks is the primary long-term drawdown process of atmospheric CO2 (Walker et al., 1981; Berner et al., 1983; Chamberlin, 1899). Dissolution of Ca–Mg silicates dissolves CO2, transports it to the oceans as bicarbonate via rivers, where it, and Ca and Mg, is sequestered into marine carbonates (Berner, 2003). Continental weathering also provides many critical nutrients to the coastal oceans (e.g., Fe, P), which fertilize organic carbon growth (Gíslason et al., 2006). Finally, continental weathering and erosion processes deliver particulate clays and reactive iron to the oceans, which enhance organic carbon burial (Kennedy and Wagner, 2011; Lalonde et al., 2012; Kennedy et al., 2014; Hawley et al., 2017). In short, weathering controls the climate both on short timescales (via organic carbon growth and burial) and long timescales (via carbonate precipitation). Chemical weathering also affects the carbonate saturation state of the oceans on millennial timescales, affecting its pH and ability to store CO2 (Archer et al., 2000).
However, there are a number of still open questions on these processes, including the ultimate controlling mechanism of silicate weathering [where both climate (via temperature and runoff) and supply of fresh primary silicates are likely important (Raymo and Ruddiman, 1992; West et al., 2005)], and the rate at which any weathering-related feedback can mitigate climate change. One reason why these questions have not been resolved is that no truly unambiguous tracers of weathering, both in the present and in the geological record, have been determined. For example, traditional radiogenic isotopes (e.g., 87Sr/86Sr or 188Os/187Os) have been used to quantify palaeo-weathering (Peucker-Ehrenbrink et al., 1995; McArthur et al., 2001; Bickle et al., 2005). However, they suffer from a strong lithological control (i.e., the isotope ratio of rivers is controlled by the isotope ratio of rocks being weathered) and often cannot distinguish between weathering of silicates and carbonates (Palmer and Edmond, 1992; Oliver et al., 2003). For this reason, other isotopic tracers have been sought, often from the “non-traditional” stable metal isotopes. In this study, we not only focus on magnesium isotopes but also compare them to silicon and lithium isotopes previously measured in the same samples.
Magnesium is of interest because it is an element directly involved in the long-term carbon cycle. It is weathered from rocks on the continents and then sequestered in the oceans into carbonates (primarily dolomites and high-Mg calcites), as well as substituting for Ca in hydrothermal interactions with basalt at mid-ocean ridges. The reactions that control this process largely involve greenschist-facies metamorphism of oceanic basalts, where Mg is basically lost into high-temperature secondary minerals, such as olivine reacting to chlorite or plagioclase reacting to albite + epidote (Humphris and Thompson, 1978; Berner et al., 1983; Albarede and Michard, 1986; Elderfield and Schultz, 1996; Gaillardet et al., 1999; Holland, 2005; Tipper et al., 2006b; Higgins and Schrag, 2015). Hence, weathering is the primary source of Mg to the oceans, while the two main sinks are carbonates and hydrothermal fluids, with additional minor sinks during low-temperature clay formation (Tipper et al., 2006b).
Magnesium has three stable isotopes (24Mg, 25Mg, and 26Mg), with isotope ratios reported as δ26Mg, which is the deviation of the 26/24 ratio from the standard DSM-3 (Galy et al., 2003). Simple dissolution imparts no Mg isotope fractionation (Wimpenny et al., 2010, 2014). During continental weathering, Mg isotopes in surface waters are strongly controlled by the dissolution mixing of isotopically light carbonates and heavy silicates (Tipper et al., 2006a, 2008; Brenot et al., 2008; Wombacher et al., 2011; Pogge von Strandmann et al., 2019b,c). Furthermore, there is additional isotopic fractionation during the formation of secondary minerals from silicate weathering, where these minerals tend to prefer the heavy Mg isotopes (Tipper et al., 2006a, b, 2008, 2010, 2012a; Pogge von Strandmann et al., 2008; Chapela Lara et al., 2017; Hindshaw et al., 2020). The same direction of fractionation occurs during Mg sorption onto mineral surfaces (Huang et al., 2012; Opfergelt et al., 2012, 2014; Pogge von Strandmann et al., 2012). Finally, the uptake of Mg by plants causes variable isotope fractionation (Black et al., 2006; Ra and Kitagawa, 2007; Bolou-Bi et al., 2012).
While a considerable body of literatures exists on Mg isotope behavior in rivers and occasionally groundwaters, much less is known about Mg isotope fractionation in continental hydrothermal fluids (Tipper et al., 2006b; Huang et al., 2018; Oelkers et al., 2019; Pogge von Strandmann et al., 2019c). This is potentially also important for the removal of Mg from seawater. While it has been generally considered that the hydrothermal removal of Mg is quantitative (and hence causes no isotope fractionation) (Holland, 2005; Higgins and Schrag, 2015), totally quantitative reactions are rarely observed in geochemistry and hence raises the interesting possibility of considerable isotope fractionation occurring. Equally, little is known about Mg isotope fractionation during uptake by phytoplankton, important for the budgets of all surface waters.
In this study, we present Mg isotope data from high- and low-temperature groundwaters that dominantly source Lake Myvatn, Iceland. The high-temperature, hydrothermal, waters have undergone water–rock interactions at temperatures up to ∼300°C (cooled during transport to 17–44°C), while the cold groundwaters are at temperatures of 3–7°C. These samples allow the comparison of the effect of temperature on groundwaters from the same area, and particularly the impact on Mg isotope behavior. We also report Mg isotope data from a time series from the single outflow from Lake Myvatn, the Laxa River. The lake has significant seasonal changes in primary productivity that strongly affect its water chemistry and provides an opportunity to study the effects of phytoplankton (including a large component of diatoms) on isotope fractionation.
Lake Myvatn (65°35′ N, 17°00′ W) is a shallow eutrophic lake in the northeast of Iceland located close to the Arctic Circle (Figure 1). The lake’s altitude is 278 m above sea level, and it has a maximum depth of ∼4.2 m, with an area of ∼37 km2. During the summer (ice-free) season (mid-May to late-September), the entire water column is well-mixed. The primary water sources to the lake are groundwaters, with the southern basin (Sydrifloi) fed by cold springs, and the northern basin (Ytrifloi) fed by hot springs (Pogge von Strandmann et al., 2016). The heat in the hot springs stems from the Namafjall and Krafla geothermal fields (Kristmannsdottir and Armannsson, 2004). Only one river (the Laxa River) drains the lake, and the water in the lake has a residence time of about 27 days (Olafsson, 1979a).
Figure 1. Location map of samples within Lake Myvatn, Iceland. The sample numbers correspond to the final numbers of each sample name in Table 1. The monitoring for 2000–2001 was made at the outlet Geirastadaskudur (Gst.).
The underlying geology of the lake is basaltic and ∼2,300 years old, formed by a major basaltic volcanic eruption. Since then, the local geology has been modified by more recent eruptions, including the Krafla Volcano’s 1725–1729 and 1975–1984 eruptions (Thordarson and Hoskuldsson, 2002; Hauptfleisch and Einarsson, 2012).
Despite winter ice cover, the lake is one of the most biologically productive in the northern hemisphere. This is due to high nutrient delivery from the groundwater sources, causing very high seasonal productivity, mainly associated with diatoms, green algae, and cyanobacteria (Olafsson, 1979b; Gislason et al., 2004; Thorbergsdottir and Gislason, 2004; Opfergelt et al., 2011). Biological activity in the lake also controls the concentrations of many dissolved metals. When the lake is covered by ice, sediment pore waters become enriched in nutrients relative to the lake waters by several orders of magnitude. These nutrients are then subsequently released by bioturbation, as well as sediment resuspension and diffusion, in the ice-free summer, leading to high nutrient fluxes and hence primary productivity (Gislason et al., 2004).
Samples were taken from both the “cold” and “hot” springs sourcing the lake in August 2009, and a monthly time series spanning March 2000 to March 2001 was taken from the Laxa River draining the lake following sampling protocols outlined in Opfergelt et al. (2011).
Materials and Methods
Major and trace element concentrations were measured by Opfergelt et al. (2011) and were determined by inductively coupled plasma mass spectrometry (ICP-MS) and ion chromatography. At each site, pH and temperature were measured in situ (Opfergelt et al., 2011). For Mg isotope analysis, sufficient water was dried down to attain ∼1 μg of Mg. This was passed through a two-step cation exchange chromatography procedure, using AG50W X-12 resin, and using 2 M HNO3 as an eluent. Isotopic analyses were conducted on a Thermo Neptune multicollector ICP-MS at the Bristol Isotope Group, normalizing to the standard DSM-3.
This procedure has been described in detail before (Pogge von Strandmann et al., 2011), and resulting standards have been shown in multiple compilations and interlaboratory comparisons (Foster et al., 2010; Pogge von Strandmann et al., 2011; Teng et al., 2015; Shalev et al., 2018). During the analyses of this study’s waters, seawater was used as a standard, with a δ26Mg = −0.82 ± 0.05‰ (n = 3), which agrees with the long-term running average of over 10 years (Foster et al., 2010) and where δ26Mg = (26Mg/24Mgsample/26Mg/24MgDSM–3 - 1) × 1,000.
The solution elemental concentrations were used to calculate mineral saturation states using PHREEQC and a combination of its standard thermodynamic database, plus THERMODEM and MINTEQ databases (Parkhurst and Appelo, 1999). These are reported as the logarithmic saturation index, SI, where positive numbers indicate supersaturation, and negative numbers indicate undersaturation.
Hydrothermal and Cold Groundwaters
All results are listed in Table 1, and the sample locations are shown in Figure 1. The groundwaters that are the primary source to Lake Myvatn are undersaturated with respect to the primary basaltic minerals of glass, olivine, pyroxene and plagioclase, as well as minerals such as calcite. Secondary minerals such as kaolinite are supersaturated, while smectites are supersaturated at relatively low pH (8–9) but undersaturated at higher pH. Secondary minerals that are known to remove Mg in seawater hydrothermal settings, such as chlorite or Mg-rich smectite (Haymon and Kastner, 1986), are also supersaturated.
Magnesium concentrations vary between 74 and 160 μmol/l, on average 50% or more lower than river waters from Iceland (Gíslason et al., 1996). The “cold” groundwaters (<7.5°C) are on average slightly less concentrated than the “hot” groundwaters (106 vs. 126 μmol/L, respectively).
Mg isotope ratios (δ26Mg) vary between −1.36 and −0.16‰, with the cold waters being on average significantly isotopically heavier (average, −0.40‰) compared to the hot groundwaters (−1.01‰). Overall, the Mg isotope ratios negatively covary with elemental concentrations such as Na, and overlap with, but extend to lower values than Icelandic river waters (Pogge von Strandmann et al., 2008). There is no evidence of mixing between different Mg sources (e.g., seawater via precipitation) in any of the samples (Figure 2).
Figure 2. Mixing diagram [1/Mg (in μmol/l) vs. δ26Mg] for the groundwaters, outflow time series, and seawater (representing Icelandic precipitation). No mixing trends are visible.
Laxa River Outflow
The Mg concentration of the outflow time series is overall slightly higher than most of the groundwater springs, with values ranging from 130 to 170 μmol/L. There is no obvious seasonal trend, although the highest concentrations are in the summer. Mg isotope ratios, in contrast, are strongly seasonal, with δ26Mg values of −0.31 to −0.24‰ in the winter (ice-covered) months, increasing to a maximum of −0.03‰ in August, within the range of Icelandic river waters (Pogge von Strandmann et al., 2008). This yearly peak is similar to that observed for Al concentrations and inversely to Si concentrations, but not other nutrients like NO3, which exhibits two minima (Opfergelt et al., 2011).
Sample MY12 was collected from the middle of the Lake at the end of August 2009, and its δ26Mg (0.00‰) is within analytical uncertainty of the outflow sample collected in August 2000 (−0.03‰), suggesting that the broad time series patterns observed in 2000–2001 are repeatable yearly cycles.
Hydrothermal and Cold Groundwaters
The primary modern sink of Mg from seawater is exchange of Mg for Ca during high-temperature basalt–water interactions (Holland, 2005). At the same time, the low-temperature alteration of the oceanic crust is likely to be at the very least a minor sink of Mg from seawater (Tipper et al., 2006b; Huang et al., 2018). The Myvatn groundwaters provide the opportunity to examine Mg isotope fractionation at various temperatures of basalt–water interaction.
Mg concentrations in the low-temperature springs (<7.5‘C) are quite variable, ranging from 70 to 150 μmol/L. However, the higher temperature springs show a negative covariation between [Mg] and temperature, such that the warmest samples have low [Mg] (Figure 3A). This is a broadly similar behavior to Mg concentrations during seawater–basalt interactions at mid-ocean ridges (Albarede and Michard, 1986; Mottl and Wheat, 1994; Elderfield and Schultz, 1996). There, [Mg] drops from concentrations of > 90% of that of seawater to < 10% of seawater once temperatures are > 40°C. By 90°C, Mg concentrations are < 2% of that of seawater (Mottl and Wheat, 1994).
Figure 3. (A) The effect of temperature on Mg concentrations during hydrothermal removal. At low temperatures, there is no effect, but as temperatures exceed 30°C, Mg concentrations both in Myvatn and seawater samples (Mottl and Wheat, 1994) rapidly decrease (note different scale for seawater samples). (B) This effect causes significant Mg isotope fractionation in the Myvatn hydrothermal springs. The trend line is logarithmic to correspond to stable isotope fractionation theory.
While overall concentrations are significantly lower in the Myvatn springs, because the original solution is not seawater, the overall relationship is similar and suggests that these results may be able to inform on Mg isotope fractionation occurring during Mg removal from seawater. It is therefore significant that the Mg isotope ratios of these springs show a clear relationship with temperature (Figure 3B). The coldest groundwater (∼3°C) has a δ26Mg value that is indistinguishable from primary basalt [−0.29 vs. −0.23‰, respectively (Pogge von Strandmann et al., 2008, 2011, 2012; Opfergelt et al., 2014)]. At slightly higher temperatures (>3, but < 7.5°C), there is more fractionation (δ26Mg = −0.68 to −0.16‰), while at higher temperatures, Mg isotope ratios in the solution continue to decrease, reaching a value of -1.36‰ at 44°C. There are two general processes that could be causing greater Mg isotope fractionation at higher temperatures: (1) standard temperature-dependent stable isotope fractionation, which scales according to 1,000 × ln(α) as a function of 106/T2 (where T is temperature in K, and αrock/soln is the isotopic fractionation factor between solution and solid, which, for geological processes, approximates the Δ7Li between two phases) (Faure, 1986); (2) removal of Mg into high-temperature secondary minerals, such as chlorite or smectite, which is the process by which most of the present-day Mg is removed from seawater (Humphris and Thompson, 1978).
Figure 4 shows that there is a covariation between the Mg isotope fractionation and temperature (shown as 1/T2), such that at higher temperature, the solution δ26Mg is more fractionated from the basaltic starting composition. Such a relationship is, in effect, opposite from what would be expected from stable isotope fractionation theory, where more fractionation occurs at lower temperature (Urey, 1947). It therefore seems likely that, while there is a temperature effect occurring, it is not dominantly due to simple stable isotope fractionation effects, suggesting that high-temperature secondary mineral formation may be playing a role.
In the case of basalt–seawater hydrothermal reactions, Mg is removed by high-temperature secondary mineral formation [such as chlorite (Berner et al., 1983; Albarede and Michard, 1986; Elderfield and Schultz, 1996; Gaillardet et al., 1999; Holland, 2005; Tipper et al., 2006b; Higgins and Schrag, 2015), which is supersaturated in these groundwaters]. The same effect can be observed by the negative covariation between Δ26Mg and the Mg/Na ratio (Figure 5). As temperature increases, more Mg is removed, removing Mg preferentially over Na [which is the most soluble major cation in basalt weathering (Gíslason et al., 1996)], and causing greater Mg isotope fractionation. The best-fit Rayleigh fractionation relationship requires an αrock–solution (isotopic fractionation factor) value of 1.00045 (while for equilibrium fractionation, this would be 1.00095, although given the potential of Mg depletion during hydrothermal reactions, a Rayleigh relationship would be expected to be more accurate). The Rayleigh fractionation is given by:
Figure 5. Mg isotope fractionation in the hydrothermal springs as a function of Mg/Na ratios. The black dotted line represents a best-fit line to the data, while the red line represents the Rayleigh fractionation trend fitted to the best-fit line. The blue line is the best-fit equilibrium fractionation line. The italicized numbers represent the required α values.
where δi is the isotope ratio (in delta notation) of the initial material, δf is that of the product, and f is the fraction of material remaining. Model estimates for the high-temperature hydrothermal removal of Mg have α values of 1.0000 (i.e., no isotopic fractionation) and values of 1.0007 for low-temperature clays (Higgins and Schrag, 2015). Experimental interactions between saponite (Mg-rich smectite) and seawater at 250 and 290°C report fractionation factors of 1.000345–1.000422 (Voigt et al., 2020), similar to our fractionation values.
The effect of hydrothermal fractionation on Mg isotopes can also be observed with a comparison to lithium isotopes. For Li, the low-temperature groundwaters are highly isotopically fractionated from the original basalt (Figure 6), while the high-temperature waters are relatively unfractionated. In effect, this represents standard stable isotope fractionation behavior, plus the effect that, in hot waters, the secondary minerals that fractionate Li isotopes are more soluble and hence precipitate less (Pogge von Strandmann et al., 2017). The positive covariation between δ7Li and δ26Mg in these waters (Figure 6), where clay formation should cause a negative covariation because clays tends to preferentially take up light Li but heavy Mg isotopes, shows that temperature-dependent stable isotope fractionation is not the controlling process for Mg isotopes.
Figure 6. Mg isotope ratios as a function of Li isotope ratios (Pogge von Strandmann et al., 2016). Fractionation-controlling mechanisms are clearly different in these two systems (see text for details).
Such a relationship could have interesting effects on the oceanic Mg isotope mass balance if the removal of Mg through hydrothermal interactions is not fully quantitative. Although it is debated whether the modern oceans are in concentration or isotopic equilibrium for Mg (Holland, 2005; Tipper et al., 2006b; Pogge von Strandmann et al., 2014; Higgins and Schrag, 2015; Gothmann et al., 2017; Voigt et al., 2020), if we assume for simplicity’s sake that they are, the mass balance can be inverted to determine the effect of isotope fractionation caused by a not fully quantitative Mg removal associated with the Rayleigh fractionation observed in our Myvatn samples. We use the original basic steady-state oceanic Mg mass balance of Tipper et al. (2006b):
where f is the fraction the source or sink flux for rivers, dolomite, hydrothermal, and the hypothetical “hydrothermal return,” which designates the efficiency of the hydrothermal sink, make up of the total source or sink. δ is the δ26Mg of the source or of seawater (SW), and Δ is the isotopic fractionation imposed by the sinks. This is then solved for three separate individual parameters: (i) the dolomite Δ26Mg [set to −1.7‰ when not varied (Pogge von Strandmann et al., 2014)], (ii) the riverine δ26Mg [set to −1.09‰ when not varied (Tipper et al., 2006b)], and (iii) the amount of the total sink comprised by dolomite [set to 15% when not varied (Tipper et al., 2006b)]. For example, as summarized in Figure 7, to maintain a steady-state seawater δ26Mg of −0.82‰ (Foster et al., 2010), if the hydrothermal sink efficiency for Mg is 90% (with a fractionation of −1.27‰, as determined from the Myvatn Rayleigh relationship), then dolomite must cause a fractionation of −0.96‰. In contrast, if that sink efficiency is 99.9%, then the dolomite fractionation must be −1.78‰.
Figure 7. Models on a steady-state modern seawater Mg budget, assuming that the hydrothermal sink is not fully efficient, and causes isotope fractionation as observed in our samples. See text for details.
Equally, if the dolomite fractionation were constant, then the amount of Mg that dolomite removes from seawater would have to decrease with increasing hydrothermal removal efficiency. Finally, if both the dolomite sink and Mg isotope fractionation were constant, then the global riverine δ26Mg would have to change to maintain a constant seawater δ26Mg, as shown in Figure 7. In other words, different parameters of the global ocean Mg mass balance can be significantly affected by the efficiency of Mg removal by the hydrothermal sink if, as shown here, there is an isotopic fractionation associated with that sink.
Laxa River Outflow
Previous work on these outflow samples show that two phytoplankton blooms occurred during the year, the first during May–June and the second during August–September, representing the spring and autumn diatom blooms (Opfergelt et al., 2011). During these periods, NO3 concentrations are at a minimum, as nutrients are being scavenged for plankton growth. In general, the groundwater is thought to act as a constant source of dissolved nutrients, split between cold springs (∼80% of input) and hot springs (∼20%). Thus, the phytoplankton blooms are controlled by the buildup of nutrients during the ice-covered season, and the spring bloom then occurs immediately after ice breakup. A second bloom then occurs in the autumn, corresponding to the second decline in dissolved nitrate (Table 1).
Silicon isotopes are strongly affected by these blooms: diatoms preferentially take up the light Si isotope, driving the lake water isotopically heavy. This results in Si isotopes showing a double peak during these periods. In contrast, Si concentrations only show a single minimum, with values declining throughout spring and increasing again after the July minimum (Figures 8, 9) (Opfergelt et al., 2011).
Figure 8. Calcite (A), forsterite, and smectite (B) mineral saturation indices for the Laxa River outflow time series.
Figure 9. Parameters for the Laxa River outflow time series: (A) pH and temperature; (B) Si and Li isotopes (Opfergelt et al., 2011; Pogge von Strandmann et al., 2016); (C) Mg isotopes and Mg concentrations. The error bars represent the 2sd analytical uncertainty. The dark gray bars represent the spring and autumn diatom blooms, while the light gray bars represent the ice-free months.
In contrast, lithium is an element largely unaffected by any biology, including phytoplankton. Li concentrations are unaffected by any of the phenomenon that affect other elements in these samples (Pogge von Strandmann et al., 2016). Equally, Li isotope ratios also exhibit broadly constant values throughout the year, showing that δ7Li is unaffected by plankton growth (Figure 9). Li isotopes are normally strongly affected by secondary mineral formation (Hindshaw et al., 2019; Pogge von Strandmann et al., 2019a). PHREEQC calculations show that the saturation state of several minerals changes dramatically during the year. For example, olivine is strongly undersaturated during the winter but approaches saturation in summer (Figure 8). Equally, secondary minerals such as smectite are supersaturated during the winter but become undersaturated in the summer. It is, however, unclear how much this change represents actual changes in mineral saturation and how much is element utilization by plankton (misinterpreted by PHREEQC as changes in saturation). The invariant Li isotopes suggest that silicate secondary mineral formation does not change much during the year, implying that plankton, rather than mineral precipitation, are primarily changing elemental concentrations during the course of the year.
Mg isotopes do not exhibit a double peak like Si isotopes or NO3 concentrations (Figure 9). Instead, they more closely follow trends like pH, which also only show a single maximum (Opfergelt et al., 2011). This implies that Mg isotopes are not being fractionated by the phytoplankton (of which diatoms make up approximately two-thirds), which is also suggested by the relatively unchanging Mg concentrations, which are clearly not limiting plankton growth. This presents a contrast to coccolith growth experiments, where the phytoplankton preferentially took up light Mg isotopes (Ra and Kitagawa, 2007). On the other hand, it has been suggested that plants have little overall effect on riverine Mg isotope ratios (Tipper et al., 2012b), although this effect does appear to depend on location and stream size (Pogge von Strandmann et al., 2019b). This may present a difference between laboratory growth experiments and natural samples, and of course, the diatoms that make up a large proportion of the Myvatn phytoplankton may also not fractionate Mg isotopes.
Equally, it does not seem that Mg isotopes are being controlled by the precipitation of silicate secondary minerals. If the PHREEQC saturation calculations are correct, maximum precipitation of minerals like Mg smectite should occur in the winter months (Figure 8). However, smectites, and most other secondary silicates, preferentially take up heavy Mg isotopes, driving residual waters isotopically light (Tipper et al., 2006a, 2008; Opfergelt et al., 2012; Pogge von Strandmann et al., 2012; Chapela Lara et al., 2017; Hindshaw et al., 2020). Thus, if precipitation of these minerals were controlling the Mg isotopes, it would be expected that, during winter, Mg isotope ratios would be lower than that of the initial basalt (∼-0.25‰). Instead, during the winter, the outflow Mg isotope ratios are within analytical uncertainty of basalt, implying only dissolution with little secondary mineral formation.
However, a secondary mineral that is suggested by PHREEQC to be undersaturated during the winter and supersaturated during the summer is calcite. From June to September, calcite precipitates in Lake Myvatn, not only suggested by our saturation calculations but also confirmed by precipitation boxes placed on the bottom of the lake (Thorbergsdottir and Gislason, 2004). Calcite preferentially takes up light Mg isotopes, driving residual waters isotopically heavy (Immenhauser et al., 2010; Wombacher et al., 2011; Mavromatis et al., 2012; Saulnier et al., 2012; Saenger and Wang, 2014). The calcite saturation state corresponds to Mg isotopes in Lake Myvatn, where values are only distinguishable from primary basalt from June to September, when they become isotopically heavier. A calcitic control on Mg isotope ratios is also suggested by a direct positive covariation between the calcite SI and δ26Mg (r2 = 0.79; Figure 10).
Figure 10. Mg isotope ratios in the Laxa River outflow time series as a function of the saturation index of calcite.
Calcite precipitation has previously been observed driving riverine Mg isotope ratios heavy in Iceland during highly supersaturated outflow from the 2010 Eyjafjallajökull eruption (Pogge von Strandmann et al., 2019c). In this case, calcite SI reached values of > 1.5, travertine precipitation was observed, and δ26Mg changed by over 2.5‰. Here at Myvatn, calcite never becomes quite as supersaturated (maximum SI = 0.8), and δ26Mg changes considerably less. This may also be exacerbated by the much larger water reservoir in Lake Myvatn and the Laxa River, compared to the smaller streams samples from Eyjafjallajökull.
Given that the groundwater influx rate into Lake Myvatn is fairly constant at 33 m3/s, it follows that the outflow rate through the Laxa River is broadly the same rate (Olafsson, 1979a; Opfergelt et al., 2011). Assuming that Mg isotopes are solely being fractionated by calcite formation, with a known fractionation factor of α = 0.9965 (Mavromatis et al., 2013; Saenger and Wang, 2014; Pogge von Strandmann et al., 2019c), up to 8% of dissolved Mg is lost into calcite (Pogge von Strandmann et al., 2019c). Following the calculations detailed in Pogge von Strandmann et al. (2019c), using the Laxa outflow rate and the volume of Lake Myvatn (∼9.25 × 107 m3), the calcite precipitation rate reaches a maximum of ∼53 t/month (∼6 mmol/month/m3) in August. Over the course of a year, this results in the precipitation of ∼240 t/year calcite. However, during the winter months, calcite is undersaturated, and much of this calcite may redissolve, as also predicted elsewhere (Thorbergsdottir and Gislason, 2004). Overall, this calcite precipitation is an order of magnitude less than the Si sequestration by diatom growth of ∼7,500 t/year (Opfergelt et al., 2011), showing that the phytoplankton (including diatoms) in Myvatn can have a significantly larger effect on the concentration of many elements compared to calcite precipitation. In comparison, the stream draining the Eyjafjallajökull eruption exhibited a Mg isotope-based calcite precipitation rate of 3,200 t/year (Pogge von Strandmann et al., 2019c).
The precipitation of that amount of calcite would sequester ∼105 t/year of atmospheric CO2, in lake-stored carbonate. Eventual dissolution of that calcite would release alkalinity, which the Laxa River would transport to the oceans, where it would form marine carbonate, again sequestering CO2 on the timescale of the completion of the carbonate weathering reaction, which is on the order of a few 1000 years (Archer et al., 2000). Thus, even if the carbonate that forms within Lake Myvatn is redissolved, it still acts as a medium-term CO2 sink, meaning that, overall, the lake is a location of carbon sequestration even when not considering any potential burial of organic carbon.
This study analyzed Mg isotope ratios in hydrothermal and groundwaters of different temperatures that form the primary input into Lake Myvatn, Iceland. There is a clear trend between their δ26Mg and temperature, such that δ26Mg values decrease with increasing temperature. Hence, as Mg is removed, Mg isotopes become increasingly fractionated, with a fractionation factor of α ∼1.00045. Hydrothermal solutions are the main sink of Mg from modern seawater. While they are generally considered to be a quantitative sink (i.e., causing no isotope fractionation), if they are not fully efficient, our data suggest that they will cause significant fractionation, with consequences for the ocean mass balance.
This study also examined Mg isotopes from a time series of the Laxa River, the only outflow from Lake Myvatn. This lake is one of the most productive in the northern hemisphere and has a strong seasonal growth rate dependence that affects nutrients like NO3 and Si. Mg isotopes do not exhibit the double peak typical of the two seasonal growth cycles but do exhibit a single peak, reaching maximum values in August. This behavior appears to correspond to the saturation behavior of calcite, which preferentially takes up light Mg isotopes driving residual water isotopically heavy. Hence, calcite is precipitated in the lake during summer and then likely dissolved again during winter when calcite is undersaturated, following which it will contribute alkalinity to the oceans. Myvatn is thus a CO2 sink even without considering the formation and burial of its organic carbon.
Data Availability Statement
All datasets generated for this study are included in the article/supplementary material.
PP analyzed the samples, interpreted the data, and wrote the manuscript. KB, SO, AE, EE, and SG provided samples and local knowledge. MM assisted with analyses. All co-authors assisted with interpretation.
PP and the analyses were funded by NERC grant NE/I020571/2, and PP and MM are also supported by ERC Consolidator grant 682760 CONTROLPASTCO2.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the two reviewers LB and MD for their comments.
Albarede, F., and Michard, A. (1986). Transfer of continental Mg, S, O and U to the Mantle through Hydrothermal alteration of the Oceanic-Crust. Chem. Geol. 57, 1–15. doi: 10.1016/0009-2541(86)90090-2
Berner, R. A., Lasaga, A. C., and Garrels, R. M. (1983). The Carbonate-Silicate Geochemical Cycle and Its Effect on Atmospheric Carbon-Dioxide over the Past 100 Million Years. Am. J. Sci. 283, 641–683. doi: 10.2475/ajs.283.7.641
Bickle, M. J., Chapman, H. J., Bunbury, J., Harris, N. B. W., Fairchild, I. J., Ahmad, T., et al. (2005). Relative contributions of silicate and carbonate rocks to riverine Sr fluxes in the headwaters of the Ganges. Geochim. Cosmochim. Acta 69, 2221–2240. doi: 10.1016/j.gca.2004.11.019
Black, J. R., Yin, Q. Z., and Casey, W. H. (2006). An experimental study of magnesium-isotope fractionation in chlorophyll-a photosynthesis. Geochim. Cosmochim. Acta 70, 4072–4079. doi: 10.1016/j.gca.2006.06.010
Bolou-Bi, E. B., Vigier, N., Poszwa, A., Boudot, J.-P., and Dambrine, E. (2012). Effects of biogeochemical processes on magnesium isotope variations in a forested catchment in the Vosges Mountains (France). Geochim. Cosmochim. Acta 87, 341–355. doi: 10.1016/j.gca.2012.04.005
Brenot, A., Cloquet, C., Vigier, N., Carignan, J., and France-Lanord, C. (2008). Magnesium isotope systematics of the lithologically varied Moselle river basin. France. Geochim. Cosmochim. Acta 72, 5070–5089. doi: 10.1016/j.gca.2008.07.027
Chapela Lara, M., Buss, H. L., Pogge von Strandmann, P. A. E., Schuessler, J. A., and Moore, O. W. (2017). The influence of critical zone processes on the Mg isotope budget in a tropical, highly weathered andesitic catchment. Geochim. Cosmochim. Acta 202, 77–100. doi: 10.1016/j.gca.2016.12.032
Gaillardet, J., Dupre, B., Louvat, P., and Allegre, C. J. (1999). Global silicate weathering and CO2 consumption rates deduced from the chemistry of large rivers. Chem. Geol. 159, 3–30. doi: 10.1016/s0009-2541(99)00031-5
Galy, A., Yoffe, O., Janney, P. E., Williams, R. W., Cloquet, C., Aldard, O., et al. (2003). Magnesium isotope heterogeneity of the isotopic standard SRM980 and new reference materials for magnesium-isotope-ratio measurements. J. Anal. Atomic Spectr. 18, 1352–1356.
Gíslason, S. R., Arnorsson, S., and Armannsson, H. (1996). Chemical weathering of basalt in southwest Iceland: effects of runoff, age of rocks and vegetative/glacial cover. Am. J. Sci. 296, 837–907. doi: 10.2475/ajs.296.8.837
Gislason, S. R., Eiriksdottir, E. S., and Olafsson, J. (2004). Chemical composition of interstitial water and diffusive fluxes within the diatomaceous sediment in Lake Myvatn. Iceland. Aqua. Ecol. 38, 163–175. doi: 10.1023/b:aeco.0000032063.15765.e4
Hauptfleisch, U., and Einarsson, A. (2012). Age of the Younger Laxá Lava and Lake Myvatn, northern Iceland, determined by AMS radiocarbon dating. Radiocarbon 54, 155–164. doi: 10.2458/azu_js_rc.v54i2.15970
Hawley, S. M., Pogge von Strandmann, P. A. E., Burton, K. W., Williams, H. M., and Gislason, S. R. (2017). Continental weathering and terrestrial (oxyhydr)oxide export: comparing glacial and non-glacial catchments in Iceland. Chem. Geol. 462, 55–66. doi: 10.1016/j.chemgeo.2017.04.026
Haymon, R. M., and Kastner, M. (1986). The formation of high temperature clay minerals from basalt alteration during hydrothermal discharge on the East Pacific Rise axis at 21°N. Geochim. Cosmochim. Acta 50, 1933–1939. doi: 10.1016/0016-7037(86)90249-8
Higgins, J. A., and Schrag, D. P. (2015). The Mg isotopic composition of Cenozoic seawater – evidence for a link between Mg-clays, seawater Mg/Ca, and climate. Earth Planet. Sci. Lett. 416, 73–81. doi: 10.1016/j.epsl.2015.01.003
Hindshaw, R. S., Tosca, R., Tosca, N. J., and Tipper, E. T. (2020). Experimental constraints on Mg isotope fractionation during clay formation: implications for the global biogeochemical cycle of Mg. Earth Planet. Sci. Lett. 531:115980. doi: 10.1016/j.epsl.2019.115980
Huang, K. J., Teng, F. Z., Wei, G. J., Ma, J. L., and Bao, Z. Y. (2012). Adsorption- and desorption-controlled magnesium isotope fractionation during extreme weathering of basalt in Hainan Island. China. Earth Planet. Sci. Lett. 35, 73–83. doi: 10.1016/j.epsl.2012.10.007
Huang, K.-J., Tosca, R., Gout, T. L., Farnan, L., Tosca, N. J., and Tipper, E. T. (2018). Magnesium isotopic composition of altered oceanic crust and the global Mg cycle. Geochim. Cosmochim. Acta 238, 357–373. doi: 10.1016/j.gca.2018.07.011
Immenhauser, A., Buhl, D., Richter, D., Niedermay, A., Riechelmann, D., Dietzel, M., et al. (2010). Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments. Geochim. Cosmochim. Acta 74, 4346–4364. doi: 10.1016/j.gca.2010.05.006
Kennedy, M. J., Löhr, S. C., Fraser, S. A., and Baruch, E. T. (2014). Direct evidence for organic carbon preservation as clay-organic nanocomposites in a Devonian black shale; from deposition to diagenesis. Earth Planet. Sci. Lett. 388, 59–70. doi: 10.1016/j.epsl.2013.11.044
Kennedy, M. J., and Wagner, T. (2011). Clay mineral continental amplifier for marine carbon sequestration in a greenhouse ocean. Proc. Natl. Acad. Sci. U.S.A. 108, 9776–9781. doi: 10.1073/pnas.1018670108
Mavromatis, V., Gautier, Q., Bosc, O., and Schott, J. (2013). Kinetics of Mg partition and Mg stable isotope fractionation during its incorporation in calcite. Geochim. Cosmochim. Acta 114, 188–203. doi: 10.1016/j.gca.2013.03.024
Mavromatis, V., Pearce, C. R., Shirokova, L. S., Bundeleva, I. A., Pokrovsky, O. S., Benezeth, P., et al. (2012). Magnesium isotope fractionation during hydrous magnesium carbonate precipitation with and without cyanobacteria. Geochim. Cosmochim. Acta 76, 161–174. doi: 10.1016/j.gca.2011.10.019
McArthur, J. M., Howarth, R. J., and Bailey, T. R. (2001). Strontium isotope stratigraphy: LOWESS version 3: best fit to the marine Sr-isotope curve for 0-509 Ma and accompanying look-up table for deriving numerical age. J. Geol. 109, 155–170. doi: 10.1086/319243
Oelkers, E. H., Butcher, R., Strandmann, P. E. A. P., Schuessler, J. A., Blankenburg, F. V., Snæbjörnsdóttir, S. O., et al. (2019). Using stable Mg isotope signatures to assess the fate of magnesium during the in situ mineralisation of CO2 and H2S at the CarbFix site in SW-Iceland. Geochim. Cosmochim. Acta 245, 542–555. doi: 10.1016/j.gca.2018.11.011
Oliver, L., Harris, N., Bickle, M. J., Chapman, H. J., Dise, N., and Horstwood, M. S. A. (2003). Silicate weathering rates decoupled from the 87Sr/86Sr ratio of the dissolved load during Himalayan erosion. Chem. Geol. 201, 119–139. doi: 10.1016/s0009-2541(03)00236-5
Opfergelt, S., Burton, K. W., George, R. B., West, A. J., Guicharnaud, R., Sigfusson, B., et al. (2014). Magnesium retention on the soil exchange complex controlling Mg isotope variations in soils, soil solutions and vegetation in volcanic soils, Iceland. Geochim. Cosmochim. Acta 125, 110–130. doi: 10.1016/j.gca.2013.09.036
Opfergelt, S., Eiriksdottir, E., Burton, K. W., Einarsson, A., Siebert, C., Gislason, S., et al. (2011). Quantifying the impact of freshwater diatom productivity on silicon isotopes and silicon fluxes: Lake Myvatn, Iceland. Earth Planet. Sci. Lett. 305, 73–82. doi: 10.1016/j.epsl.2011.02.043
Opfergelt, S., Georg, R. B., Delvaux, B., Cabidoche, Y.-M., Burton, K. W., Halliday, A. N., et al. (2012). Mechanisms of magnesium isotope fractionation in volcanic soil weathering sequences. Guadeloupe. Earth Planet. Sci. Lett. 341, 176–185. doi: 10.1016/j.epsl.2012.06.010
Parkhurst, D. L., and Appelo, C. A. J. (1999). User’s guide to PHREEQC (version 2) - a Computer Program For Speciation, Batch-Reaction, One-Dimensional Transport, And Inverse Geochemical Calculations. Reston: USGS.
Pogge von Strandmann, P. A. E., Burton, K. W., James, R. H., van Calsteren, P., and Gislason, S. R. (2008). The influence of weathering processes on riverine magnesium isotopes in a basaltic terrain. Earth Planet. Sci. Lett. 276, 187–197. doi: 10.1016/j.epsl.2008.09.020
Pogge von Strandmann, P. A. E., Burton, K. W., Opfergelt, S., Eiríksdóttir, E. S., Murphy, M. J., and Gislason, S. R. (2016). The effect of hydrothermal spring weathering processes and primary productivity on lithium isotopes: Lake Myvatn. Iceland. Chem. Geol. 445, 4–13. doi: 10.1016/j.chemgeo.2016.02.026
Pogge von Strandmann, P. A. E., Elliott, T., Marschall, H., Coath, C., Lai, Y. J., Jeffcoate, A. B., et al. (2011). Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths. Geochim. Cosmochim. Acta 75, 5247–5268. doi: 10.1016/j.gca.2011.06.026
Pogge von Strandmann, P. A. E., Forshaw, J., and Schmidt, D. N. (2014). Modern and Cenozoic records of seawater magnesium from foraminiferal Mg isotopes. Biogeosciences 11, 5155–5168. doi: 10.5194/bg-11-5155-2014
Pogge von Strandmann, P. A. E., Fraseret, W. T., Hammond, S. J., Tarbuck, G., Wood, I. G., Oelkers, E. H., et al. (2019a). Experimental determination of Li isotope behaviour during basalt weathering. Chem. Geol. 517, 34–43. doi: 10.1016/j.chemgeo.2019.04.020
Pogge von Strandmann, P. A. E., Hendry, K. R., Hatton, J. E., and Robinson, L. F. (2019b). The Response of Magnesium, Silicon, and Calcium Isotopes to rapidly uplifting and weathering terrains: South Island, New Zealand. Front. Earth Sci. 7:240. doi: 10.3389/feart.2019.00240
Pogge von Strandmann, P. A. E., Olsson, J., Luu, T.-H., Gislason, S. R., and Burton, K. W. (2019c). Using Mg isotopes to estimate natural calcite compositions and precipitation rates during the 2010 Eyjafjallajökull eruption. Front. Earth Sci. 7:6. doi: 10.3389/feart.2019.00006
Pogge von Strandmann, P. A. E., Opfergelt, S., Lai, Y. J., Sigfusson, B., Gislason, S., Burton, K. W., et al. (2012). Lithium, magnesium and silicon isotope behaviour accompanying weathering in a basaltic soil and pore water profile in Iceland. Earth Planet. Sci. Lett. 33, 11–23. doi: 10.1016/j.epsl.2012.05.035
Pogge von Strandmann, P. A. E., Vaks, A., Bar-Matthews, M., Ayalon, A., Jacob, E., and Henderson, G. E. (2017). Lithium isotopes in speleothems: temperature-controlled variation in silicate weathering during glacial cycles. Earth Planet. Sci. Lett. 469, 64–74. doi: 10.1016/j.epsl.2017.04.014
Saenger, C., and Wang, Z. (2014). Magnesium isotope fractionation in biogenic and abiogenic carbonates: implications for paleoenvironmental proxies. Quat. Sci. Rev. 90, 1–21. doi: 10.1016/j.quascirev.2014.01.014
Saulnier, S., Rollion-Bard, C., Vigier, N., and Chaussidon, M. (2012). Mg isotope fractionation during calcite precipitation: an experimental study. Geochim. Cosmochim. Acta 91, 75–91. doi: 10.1080/10256016.2018.1430692
Shalev, N., Farkas, J., Fietzke, J., Novák, M., Schuessler, J. A., von Strandmann, P. A. E. P., et al. (2018). Mg Isotope inter-laboratory comparison of reference materials from earth-surface low-temperature environments. Geostand. Geoanal. Res. 42, 205–221. doi: 10.1111/ggr.12208
Teng, F. Z., Yin, Q.-Z., Ullmann, C. V., Chakrabarti, R., von Strandmann, P. A. E. P., Yang, W., et al. (2015). Interlaboratory comparison of magnesium isotopic compositions of 12 felsic to ultramafic igneous rock standards analyzed by MC-ICPMS. Geochem. Geophys. Geosyst. 16, 3197–3209. doi: 10.1002/2015gc005939
Tipper, E. T., Calmels, D., Gaillardet, J., Louvat, P., Capmas, F., and Dubacq, B. (2012a). Positive correlation between Li and Mg isotope ratios in the river waters of the Mackenzie Basin challenges the interpretation of apparent isotopic fractionation during weathering. Earth Planet. Sci. Lett. 333, 35–45. doi: 10.1016/j.epsl.2012.04.023
Tipper, E. T., Lemarchand, E., Hindshaw, R. S., Reynolds, B. C., and Bourdon, B. (2012b). Seasonal sensitivity of weathering processes: hints from magnesium isotopes in a glacial stream. Chem. Geol. 31, 80–92. doi: 10.1016/j.chemgeo.2012.04.002
Tipper, E. T., Gaillardet, J., Louvat, P., Capmas, F., and White, A. F. (2010). Mg isotope constraints on soil pore-fluid chemistry: evidence from Santa Cruz. California. Geochim. Cosmochim. Acta 74, 3883–3896. doi: 10.1016/j.gca.2010.04.021
Tipper, E. T., Galy, A., and Bickle, M. (2008). Calcium and magnesium isotope systematics in rivers draining the Himalaya-Tibetan-Plateau region: Lithological or fractionation control? Geochim. Cosmochim. Acta 72, 1057–1075. doi: 10.1016/j.gca.2007.11.029
Tipper, E. T., Galy, A., and Bickle, M. J. (2006a). Riverine evidence for a fractionated reservoir of Ca and Mg on the continents: implications for the oceanic Ca cycle. Earth Planet. Sci. Lett. 247, 267–279. doi: 10.1016/j.epsl.2006.04.033
Tipper, E. T., Galy, A., Gailardet, J., Elderfield, H., and Carder, E. A. (2006b). The magnesium isotope budget of the modern ocean: constraints from riverine magnesium isotope ratios. Earth Planet. Sci. Lett. 250, 241–253. doi: 10.1016/j.epsl.2006.07.037
Voigt, M., Pearce, C. R., Fries, D. M., Baldermann, A., and Oelkers, E. H. (2020). Magnesium isotope fractionation during hydrothermal seawater-basalt interaction. Geochim. Cosmochim. Acta 272, 21–35. doi: 10.1016/j.gca.2019.12.026
Wimpenny, J., Colla, C. A., Yin, Q. Z., Rustad, J. R., and Casey, W. H. (2014). Investigating the behaviour of Mg isotopes during the formation of clay minerals. Geochim. Cosmochim. Acta 128, 178–194. doi: 10.1016/j.gca.2013.12.012
Wimpenny, J., Gíslason, S. R., James, R. H., Gannoun, A., Pogge von Strandmann, P. A. E., and Burton, K. W. (2010). The behaviour of Li and Mg isotopes during primary phase dissolution and secondary mineral formation in basalt. Geochim. Cosmochim. Acta 74, 5259–5279. doi: 10.1016/j.gca.2010.06.028
Keywords: weathering (IGC: D3/D5/D6), isotope geochemistry, carbonate chemistry, hydrothermal spring, groundwater (G.W.), phytoplankton
Citation: Pogge von Strandmann PAE, Burton KW, Opfergelt S, Eiríksdóttir ES, Murphy MJ, Einarsson A and Gislason SR (2020) Hydrothermal and Cold Spring Water and Primary Productivity Effects on Magnesium Isotopes: Lake Myvatn, Iceland. Front. Earth Sci. 8:109. doi: 10.3389/feart.2020.00109
Received: 02 January 2020; Accepted: 24 March 2020;
Published: 30 April 2020.
Edited by:Paul Frederick Dennis, University of East Anglia, United Kingdom
Reviewed by:Martin Dietzel, Graz University of Technology, Austria
Luke Bridgestock, University of Oxford, United Kingdom
Copyright © 2020 Pogge von Strandmann, Burton, Opfergelt, Eiríksdóttir, Murphy, Einarsson and Gislason. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Philip A. E. Pogge von Strandmann, email@example.com