1
Title: Two-billion-year-old evaporites capture Earth’s great oxidation
1
Authors: C.L. Blättler1*, M.W. Claire2,3,4, A.R. Prave2, K. Kirsimäe5, J.A. Higgins1, P.V.
2
Medvedev6, A.E. Romashkin6, D.V. Rychanchik6, A.L. Zerkle2,3, K. Paiste7, T. Kreitsmann5, I.L.
3
Millar8, J.A. Hayles9, H. Bao10, A.V. Turchyn11, M.R. Warke2, A. Lepland12,7,5,13
4
Affiliations:
5
1Department of Geosciences, Princeton University, Princeton, NJ 08544, USA.
6
2School of Earth and Environmental Sciences, University of St Andrews, St Andrews KY16
7
9AL, Scotland/UK.
8
3Centre for Exoplanet Science, University of St Andrews, St Andrews KY16 9AL, Scotland/UK.
9
4Blue Marble Space Institute of Science, 1001 4th Ave., Suite 3201 Seattle, WA 98154, USA.
10
5University of Tartu, Department of Geology, 50411 Tartu, Estonia.
11
6Institute of Geology, Karelian Science Centre, Pushkinskaya 11, 185610 Petrozavodsk, Russia.
12
7Centre for Arctic Gas Hydrate, Environment and Climate, Department of Geosciences, UiT The
13
Arctic University of Norway, 9037 Tromsø, Norway.
14
8NERC Isotope Geosciences Laboratory, British Geological Survey, Keyworth, Nottingham
15
NG12 5GG, UK.
16
9Department of Earth Science, Rice University, 6100 Main St., Houston, TX 77005, USA.
17
10Department of Geology and Geophysics, E235 Howe-Russell Geoscience Complex, Louisiana
18
State University, Baton Rouge, Louisiana 70803, USA.
19
11Department of Earth Sciences, University of Cambridge, Downing Street, Cambridge CB2
20
3EQ, UK.
21
12Geological Survey of Norway, 7491 Trondheim, Norway.
22
13Tallinn University of Technology, Institute of Geology, 19086 Tallinn, Estonia.
23
*Correspondence to: blattler@princeton.edu.
24
This is the author's version of the work. It is posted here by permission of the AAAS for personal use, not for redistribution. The definitive version was published in Science , (), doi:
10.1126/science.aar2687"
2 Abstract:
25
Major changes in atmospheric and ocean chemistry occurred in the Paleoproterozoic Era
26
(2.5–1.6 billion years ago). Increasing oxidation dramatically changed Earth’s surface, but few
27
quantitative constraints exist on this important transition. This study describes the
28
sedimentology, mineralogy, and geochemistry of a two-billion-year-old and ~800 meter-thick
29
evaporite succession from the Onega Basin in Russian Karelia. The deposit consists of a basal
30
unit dominated by halite (~100 meters) followed by anhydrite-magnesite (~500 meters) and
31
dolomite-magnesite (~200 meters) dominated units. The evaporite minerals provide a robust
32
constraint that marine sulfate concentrations were at least 10 mmol/kg, representing an oxidant
33
reservoir equivalent to over 20% of the modern ocean-atmosphere oxidizing capacity. These
34
results show that substantial amounts of surface oxidant accumulated during this critical
35
transition in Earth’s oxygenation.
36
37
One Sentence Summary:
38
A well preserved, extensive evaporite deposit reveals that a large oxidant reservoir in the form of
39
marine sulfate accumulated two billion years ago.
40
41
Main Text:
42
The geological record preserves evidence of Earth’s dynamic surface oxygenation
43
(reviews by 1, 2), but quantifying this history remains a challenge. The presence/absence of
44
redbeds, banded iron formations, and detrital grains of pyrite and uraninite (1, 3) qualitatively
45
indicate increasing oxidation during the Paleoproterozoic, and the disappearance of large-
46
3
magnitude mass-independent fractionation (MIF) of sulfur isotopes at 2.4–2.3 Ga (4, 5) shows
47
that the atmosphere exceeded a redox threshold of ~1 ppm pO2 (6, 7). However, this limit
48
reflects only a tiny fraction of the potential surface oxidant budget and does not capture
49
subsequent redox changes in the Earth system. Today, marine sulfate ([SO42-
(aq)] = 28 mmol/kg)
50
constitutes one of the largest surface oxidant reservoirs, equivalent to almost twice the modern
51
atmospheric O2 inventory. Therefore, quantitative bounds on marine sulfate concentrations are
52
essential for constraining the net electron balance and accumulation of oxidants on Earth’s
53
surface.
54
Sedimentary evaporite minerals are one of the best archives of ancient seawater
55
chemistry (e.g. 8, 9) and specific isotopic signals (e.g. 10, 11). Unfortunately, most Precambrian
56
evaporite deposits consist of pseudomorphic replacements (12), and until recently, the oldest
57
known preserved evaporitic halite and bedded sulfates dated from ~830 Ma (13, 14) and ~1.2 Ga
58
(15), respectively. This study presents analyses from a remarkably preserved ~2.0 Ga marine
59
evaporite succession bearing carbonates, sulfates, halites, and bittern salts. This succession was
60
discovered during the 2007–2009 drilling of the Onega Parametric Hole (OPH), which
61
intersected 2.9 km of Paleoproterozoic sedimentary and volcanic rocks and 600 m of Archean
62
gneiss in the Onega Basin, Karelia, Russia (16, 17). By extending the record of extensive marine
63
evaporites by almost a billion years, core samples from the OPH provide a unique window into
64
surface conditions and redox balance in the aftermath of the initial rise of oxygen on Earth.
65
The interval of the OPH studied here lies between 2940 and 2115 m depth and
66
corresponds to the ~2.0 Ga Tulomozero Formation (age discussed in supplementary materials).
67
In other cores and outcrop exposures, this formation contains abundant pseudomorphic
68
replacements of evaporite minerals (18, 19), but original evaporite minerals are only preserved in
69
4
the OPH where they define three major units. Unit 1 (2940–2833 m; 40% average core recovery
70
in cored intervals) comprises dark red-pink, recrystallized halite with intraclasts of anhydrite,
71
magnesite, and mudstone (Fig. 1A,B) and ~10% various magnesium- and potassium-sulfate salts
72
(Fig. 2). Unit 2 (2833–2330 m; 56% average core recovery in cored intervals) consists of
73
decimeter- to meter-scale interlayered anhydrite, magnesite, and laminated dark-gray to red
74
mudstone (Fig. 1C) with minor glauberite, gypsum, and halite in its lower part. Unit 3 (2330–
75
2115 m; 44% average core recovery in cored intervals) is typified by pink-tan, commonly
76
microbially laminated dolostone (Fig. 1D), laminated red-brown-gray mudstone, and variable
77
amounts of magnesite; quartz and dolomite pseudomorphs of calcium-sulfate minerals occur
78
throughout the lower half of this unit, forming laths, nodules, discs, swallow-tail crystals, and
79
chicken-wire fabric (Fig. S2).
80
Considering the Tulomozero Formation in the OPH and its development across the
81
18,000 km2 of the Onega Basin (18, 19), the interpreted depositional setting is a restricted marine
82
embayment with sabkha/coastal plain and shallow-marine environments (Fig. S3). The OPH
83
succession captures a decreasing degree of evaporation from a state of halite and magnesium-
84
/potassium-sulfate saturation (Unit 1), through calcium-sulfate deposition (Unit 2), and then to
85
more open marine dolostone precipitation (Unit 3). The extent, thickness (>800 m in the OPH
86
core), and mineral sequence of the evaporite succession are comparable to Phanerozoic evaporite
87
basins.
88
Isotope geochemistry of the OPH evaporites presents an opportunity to investigate
89
Earth’s ocean-atmosphere system at ~2.0 Ga (Fig. 2). Quadruple sulfur isotope analyses of
90
samples from Units 1 and 2 reveal ∆33S and ∆36S values indistinguishable from zero, confirming
91
an oxic atmosphere (6, 7) and production of the Tulomozero Formation sulfate well after that
92
5
atmospheric transition occurred. Triple oxygen isotope measurements of sulfates yield resolvably
93
negative ∆17O values. Oxygen MIF derived from atmospheric O3/O2 photochemistry (20) cannot
94
be ruled out, but the small magnitude ∆17O signals preclude a unique interpretation (see
95
supplementary materials) and a quantification of pO2 is not possible.
96
The mass-dependent behavior of sulfur and calcium isotopes provides compositional
97
constraints on ancient seawater. Sulfate δ34S values lie between 5 and 7‰, except for the
98
uppermost sample in Unit 2 (further discussion in supplementary materials). Given the small 34S
99
enrichment during sulfate evaporite formation (10), the seawater sulfate δ34S composition is
100
estimated to have been 4–6‰ during deposition of Units 1 and 2. The homogeneous sulfur
101
isotopic composition across ~400 m of OPH stratigraphy and the composition and sheer volume
102
of evaporite minerals suggest that the OPH evaporites must have derived from seawater and
103
preserve robust, primary isotopic signals. Additionally, the presence of halite and highly soluble
104
bittern salts in Unit 1 argues against interaction with large volumes of fluid and supports the
105
interpretation of primary isotopic ratios for the major mineral-forming elements. Calcium isotope
106
ratios show a clear stratigraphic relationship following mineralogical trends, with the highest
107
δ44/40Ca values in Unit 1, decreasing values in Unit 2, and even lower values in Unit 3. Three
108
samples in Unit 3 with the lowest δ44/40Ca values (–1.6 to –1.4‰ relative to modern seawater)
109
also exhibit relatively high strontium content, with one sample containing minor relict aragonite;
110
these observations indicate that the bulk sediment likely contained primary aragonite which has
111
now largely been converted to other carbonate minerals (see supplementary materials).
112
Excluding those samples associated with aragonite, where mineralogy rather than evaporitic
113
processes is the first-order control on calcium isotope ratios, the increase in δ44/40Ca values in the
114
more highly evaporated facies is consistent with evaporite precipitation driving isotopic
115
6
distillation of calcium, the magnitude of which is sensitive to the initial composition of seawater
116
(21).
117
An estimate for seawater sulfate concentrations at ~2.0 Ga can be derived from the
118
observed sequence, mineralogy, and calcium isotope ratios of the OPH evaporites (Fig. 3).
119
Constraints were assessed by comparison to batch evaporation simulations with varying initial
120
ion concentrations. The relative concentrations of calcium and sulfate are the principal controls
121
governing the precipitation sequence that is expressed in the OPH, as in modern evaporites, of
122
carbonates (Unit 3), followed by calcium sulfates (Unit 2), halite, and eventually magnesium
123
sulfates (Unit 1). As such, the OPH preserves a reversed evaporite sequence where the degree of
124
evaporation decreases stratigraphically upwards, progressing from the most evolved brine at the
125
base of the Tulomozero Formation towards more open marine conditions. During a forward
126
evaporite sequence, calcium precipitates as sulfate minerals and minor carbonate with an isotopic
127
fractionation, such that calcium in the remaining brine becomes enriched in the heavier isotopes
128
through Rayleigh distillation. If sufficient sulfate is present to remove a large fraction of the
129
original calcium content of seawater, later calcium-bearing phases can record large δ44/40Ca
130
enrichments (21). The ~1‰ δ44/40Ca range captured in the OPH between shallow-marine
131
carbonates in Unit 3 and halite-hosted anhydrite in Unit 1 therefore places a lower limit on
132
sulfate concentrations. Together with mineralogical constraints, and assuming modern
133
concentrations of other major ions and a conservative interpretation of δ44/40Ca values (but see
134
supplementary materials for further discussion), the minimum sulfate concentration consistent
135
with these observations is ~10 mmol/kg (Fig. 3).
136
The OPH core provides quantitative evidence that marine sulfate concentrations at ~2.0
137
Ga were at least a third those of modern seawater. This constraint validates assertions that a large
138
7
Paleoproterozoic sulfate reservoir existed (18, 22) and increases fourfold the previous lower
139
bound of 2.5 mmol/kg sulfate, derived from observing that gypsum evaporites precipitated
140
before halite (22, 23). Although the ancient ocean volume is unknown, 10 mmol/kg sulfate in a
141
modern-sized ocean represents an oxidant reservoir equivalent to 23% of the present ocean-
142
atmosphere oxidizing capacity (or 62% of the present atmospheric O2 inventory). The growth of
143
such a reservoir from < 200 µmol/kg sulfate in the Neoarchean (24) would account for a redox
144
imbalance of at least 8–24×1010 mol/yr in equivalent moles of O2 produced or organic carbon
145
buried over 100–300 Myr. Compared to estimates of modern organic carbon burial (5×1012
146
mol/yr (25)), the accumulation of such a sizable fluid oxidant reservoir within the given time
147
constraints can be explained by either a large and rapid decline in reductant sinks (i.e. sulfide) or
148
a prolonged 2–5% imbalance over 108-year timescales. In either case, the observations suggest
149
that a sustained increase in net O2 production occurred in the Paleoproterozoic.
150
The geologically rapid growth of a massive sulfate reservoir, with or without a
151
commensurate increase in atmospheric O2, also has implications for feedbacks between the
152
global biogeochemical cycles of O2 and CO2 and Earth’s climate. In particular, the oxidation of
153
large amounts of reduced sulfur requires additional sources of carbon to offset the inferred
154
organic carbon burial (the initial source of oxygen) and prevent catastrophic cooling (26).
155
Models for Earth’s oxidation must balance these considerations as well as the new evidence for a
156
substantial oxidant reservoir in the form of marine sulfate. Finally, although sulfate likely
157
represented the largest oxidant reservoir during deposition of the Tulomozero Formation, its
158
concentration subsequently decreased such that evaporites at ~1.9 Ga no longer precipitated
159
gypsum before halite (12, 27) and a fundamental change in the sedimentary sulfur isotopic
160
composition occurred (28). The apparent transient accumulation of surface oxidants is not yet
161
8
well understood (2, 29) but implies a protracted re-organization of the global redox budget on the
162
timescales of sedimentary recycling of pyrite and organic carbon (i.e. hundreds of millions of
163
years). Regardless of the mechanisms involved, the observations presented here from the OPH
164
core document a remarkably large oxidant pool at ~2.0 Ga – a pivotal new constraint on the
165
history of Earth’s oxidation.
166
9 References and Notes:
167
1. D. E. Canfield, The early history of atmospheric oxygen: Homage to Robert M. Garrels. Annu. Rev.
168
Earth Planet. Sci. 33, 1 (2005).
169
2. T. W. Lyons, C. T. Reinhard, N. J. Planavsky, The rise of oxygen in Earth's early ocean and 170
atmosphere. Nature 506, 307 (2014).
171
3. H. Holland, in Early Life on Earth: Nobel Symposium, 84. (Columbia University Press, 1994).
172
4. J. Farquhar, H. Bao, M. Thiemens, Atmospheric influence of Earth's earliest sulfur cycle. Science 173
289, 756 (2000).
174
5. G. Luo et al., Rapid oxygenation of Earth’s atmosphere 2.33 billion years ago. Science Advances 2, 175
e1600134 (2016).
176
6. A. A. Pavlov, J. F. Kasting, Mass-independent fractionation of sulfur isotopes in Archean sediments:
177
Strong evidence for an anoxic Archean atmosphere. Astrobiology 2, 27 (2002).
178
7. K. Zahnle, M. Claire, D. Catling, The loss of mass‐independent fractionation in sulfur due to a 179
Palaeoproterozoic collapse of atmospheric methane. Geobiology 4, 271 (2006).
180
8. T. K. Lowenstein, M. N. Timofeeff, S. T. Brennan, L. A. Hardie, R. V. Demicco, Oscillations in 181
Phanerozoic seawater chemistry: Evidence from fluid inclusions. Science 294, 1086 (2001).
182
9. J. Horita, H. Zimmermann, H. D. Holland, Chemical evolution of seawater during the Phanerozoic:
183
Implications from the record of marine evaporites. Geochim. Cosmochim. Acta 66, 3733 (2002).
184
10. G. E. Claypool, W. T. Holser, I. R. Kaplan, H. Sakai, I. Zak, The age curves of sulfur and oxygen 185
isotopes in marine sulfate and their mutual interpretation. Chem. Geol. 28, 199 (1980).
186
11. H. Bao, J. R. Lyons, C. Zhou, Triple oxygen isotope evidence for elevated CO2 levels after a 187
Neoproterozoic glaciation. Nature 453, 504 (2008).
188
12. M. C. Pope, J. P. Grotzinger, Paleoproterozoic Stark Formation, Athapuscow Basin, northwest 189
Canada: Record of cratonic-scale salinity crisis. J. Sediment. Res. 73, 280 (2003).
190
13. N. Spear et al., Analyses of fluid inclusions in Neoproterozoic marine halite provide oldest 191
measurement of seawater chemistry. Geology 42, 103 (2014).
192
14. J. F. Lindsay, Upper Proterozoic evaporites in the Amadeus basin, central Australia, and their role in 193
basin tectonics. Geol. Soc. Am. Bull. 99, 852 (1987).
194
15. L. C. Kah, T. W. Lyons, J. T. Chesley, Geochemistry of a 1.2 Ga carbonate-evaporite succession, 195
northern Baffin and Bylot Islands: Implications for Mesoproterozoic marine evolution. Precambrian 196
Research 111, 203 (2001).
197
16. A. F. Morozov et al., Rock salt mass in the Paleoproterozoic sequence of the Onega Trough in 198
Karelia (from the Onega Parametric Well data). Doklady Earth Sciences 435, 1483 (2010).
199
10
17. V. A. Krupenik, K. Y. Sveshnikova, “Correlation of the Onega Parametric Hole with the reference 200
sections of the Onega Structure” in The Onega Palaeoproterozoic structure (Geology, tectonics, deep 201
structure and minerogeny). L. V. Glushanin, N. V. Sharov, V. V. Shchiptsov, Eds. (Institute of 202
Geology, Karelian Research Centre RAS, Petrozavodsk, 2011), pp. 190–195 (in Russian).
203
18. V. A. Melezhik, A. E. Fallick, D. V. Rychanchik, A. B. Kuznetsov, Palaeoproterozoic evaporites in 204
Fennoscandia: Implications for seawater sulphate, the rise of atmospheric oxygen and local 205
amplification of the δ13C excursion. Terra Nova 17, 141 (2005).
206
19. A. T. Brasier et al., Coastal sabkha dolomites and calcitised sulphates preserving the Lomagundi- 207
Jatuli carbon isotope signal. Precambrian Research 189, 193 (2011).
208
20. X. Cao, H. Bao, Dynamic model constraints on oxygen-17 depletion in atmospheric O2 after a 209
snowball Earth. PNAS 110, 14546 (2013).
210
21. C. L. Blättler, J. A. Higgins, Calcium isotopes in evaporites record variations in Phanerozoic seawater 211
SO4 and Ca. Geology 42, 711 (2014).
212
22. M. Reuschel et al., Isotopic evidence for a sizeable seawater sulfate reservoir at 2.1 Ga. Precambrian 213
Research 192, 78 (2012).
214
23. S. Schröder, A. Bekker, N. J. Beukes, H. Strauss, H. S. van Niekerk, Rise in seawater sulphate 215
concentration associated with the Paleoproterozoic positive carbon isotope excursion: Evidence from 216
sulphate evaporites in the ~2.2–2.1 Gyr shallow-marine Lucknow Formation, South Africa. Terra 217
Nova 20, 108 (2008).
218
24. K. S. Habicht, M. Gade, B. Thamdrup, P. Berg, D. E. Canfield, Calibration of sulfate levels in the 219
Archean ocean. Science 298, 2372 (2002).
220
25. R. A. Berner, D. E. Canfield, A new model for atmospheric oxygen over Phanerozoic time. Am. J.
221
Sci. 289, 333 (1989).
222
26. A. Bachan, L. R. Kump, The rise of oxygen and siderite oxidation during the Lomagundi Event.
223
PNAS 112, 6562 (2015).
224
27. J. P. Grotzinger, J. F. Kasting, New constraints on Precambrian ocean composition. J. Geol. 101, 235 225
(1993).
226
28. N. J. Planavsky, A. Bekker, A. Hofmann, J. D. Owens, T. W. Lyons, Sulfur record of rising and 227
falling marine oxygen and sulfate levels during the Lomagundi event. PNAS 109, 18300 (2012).
228
29. A. Bekker, H. D. Holland, Oxygen overshoot and recovery during the early Paleoproterozoic. Earth 229
Planet. Sci. Lett. 317, 295 (2012).
230
30. V. A. Krupenik, A. M. Akhmedov, K. Y. Sveshnikova, “Isotopic composition of carbon, oxygen and 231
sulphur in the Ludicovian and Jatulian rocks” in The Onega Palaeoproterozoic structure (Geology, 232
tectonics, deep structure and minerogeny). L. V. Glushanin, N. V. Sharov, V. V. Shchiptsov, Eds.
233
(Institute of Geology, Karelian Research Centre RAS, Petrozavodsk, 2011), pp. 250–255 (in 234
Russian).
235
11
31. V. A. Melezhik et al., Eds., The Palaeoproterozoic of Fennoscandia as Context for the 236
Fennoscandian Arctic Russia – Drilling Early Earth Project, vol. 1 of Reading the Archive of Earth's 237
Oxygenation (Springer-Verlag, 2013).
238
32. O. A. Eusipco, I. V. Neronova, “Geophysical studies of the borehole” in The Onega 239
Palaeoproterozoic structure (Geology, tectonics, deep structure and minerogeny). L. V. Glushanin, 240
N. V. Sharov, V. V. Shchiptsov, Eds. (Institute of Geology, Karelian Research Centre RAS, 241
Petrozavodsk, 2011), pp. 237-247 (in Russian).
242
33. O. A. Eusipco, I. V. Neronova, N. V. Sharov, Geophysical Study of the Onega Parametric Borehole, 243
Lennex Corporation Report. (Nobel Press, 2014), pp. 62 (in Russian).
244
34. P. V. Medvedev, V. V. Makarikhin, A. I. Golubev, D. V. Rychanchik, N. N. Trofimov, “The Jatuli”
245
in The Onega Palaeoproterozoic structure (Geology, tectonics, deep structure and minerogeny). L. V.
246
Glushanin, N. V. Sharov, V. V. Shchiptsov, Eds. (Institute of Geology, Karelian Research Centre 247
RAS, Petrozavodsk, 2011), pp. 52–67 (in Russian).
248
35. Y. A. Karpychev, Variations in the sedimentation in Kara Bogaz Gol Bay related to sea level 249
fluctuations during the Novocaspian time. Oceanology 47, 857 (2007).
250
36. A. N. Kosarev, A. G. Kostianoy, I. S. Zonn, Kara-Bogaz-Gol Bay: Physical and chemical evolution.
251
Aquatic Geochemistry 15, 223 (2009).
252
37. L. R. Kump et al., Isotopic evidence for massive oxidation of organic matter following the Great 253
Oxidation Event. Science 334, 1694 (2011).
254
38. A. V. Stepanova, A. V. Samsonov, A. N. Larionov, The final episode of middle Proterozoic 255
magmatism in the Onega structure: Data on trans-Onega dolerites. Proceedings of the Karelian 256
Research Centre of the Russian Academy of Sciences: Precambrian Geology 1, 3 (1998).
257
39. G. V. Ovchinnikova et al., Pb-Pb age of Jatulian carbonate rocks: the Tulomozero Formation of 258
southeast Karelia. Stratigraphy and Geological Correlation 15, 359 (2007).
259
40. I. S. Puchtel et al., Petrology of mafic lavas within the Onega plateau, central Karelia: Evidence for 260
2.0 Ga plume-related continental crustal growth in the Baltic Shield. Contrib. Mineral. Petrol. 130, 261
134 (1998).
262
41. A. P. Martin et al., Multiple Palaeoproterozoic carbon burial episodes and excursions. Earth Planet.
263
Sci. Lett. 424, 226 (2015).
264
42. J. C. Taylor, Computer programs for standardless quantitative analysis of minerals using the full 265
powder diffraction profile. Powder Diffraction 6, 2 (1991).
266
43. H. Bao, Purifying barite for oxygen isotope measurement by dissolution and reprecipitation in a 267
chelating solution. Analytical Chemistry 78, 304 (2006).
268
44. H. Bao, M. H. Thiemens, Generation of O2 from BaSO4 using a CO2−laser fluorination system for 269
simultaneous analysis of δ18O and δ17O. Analytical Chemistry 72, 4029 (2000).
270
45. H. Bao, X. Cao, J. A. Hayles, Triple oxygen isotopes: fundamental relationships and applications.
271
Annu. Rev. Earth Planet. Sci. Lett. 44, 463 (2016).
272
12
46. P. W. Crockford et al., Triple oxygen and multiple sulfur isotope constraints on the evolution of the 273
post-Marinoan sulfur cycle. Earth Planet. Sci. Lett. 435, 74 (2016).
274
47. D. E. Canfield, R. Raiswell, J. T. Westrich, C. M. Reaves, R. A. Berner, The use of chromium 275
reduction in the analysis of reduced inorganic sulfur in sediments and shales. Chem. Geol. 54, 149 276
(1986).
277
48. H. G. Thode, J. Monster, H. B. Dunford, Sulphur isotope geochemistry. Geochim. Cosmochim. Acta 278
25, 159 (1961).
279
49. G. Antler, A. V. Turchyn, V. Rennie, B. Herut, O. Sivan, Coupled sulfur and oxygen isotope insight 280
into bacterial sulfate reduction in the natural environment. Geochim. Cosmochim. Acta 118, 98 281
(2013).
282
50. C. L. Blättler, N. R. Miller, J. A. Higgins, Mg and Ca isotope signatures of authigenic dolomite in 283
siliceous deep-sea sediments. Earth Planet. Sci. Lett. 419, 32 (2015).
284
51. E. D. Young, A. Galy, H. Nagahara, Kinetic and equilibrium mass-dependent isotope fractionation 285
laws in nature and their geochemical and cosmochemical significance. Geochim. Cosmochim. Acta 286
66, 1095 (2002).
287
52. G. Caro, D. A. Papanastassiou, G. J. Wasserburg, 40K–40Ca isotopic constraints on the oceanic 288
calcium cycle. Earth Planet. Sci. Lett. 296, 124 (2010).
289
53. A. Heuser, A. Eisenhauer, The calcium isotope composition (δ44/40Ca) of NIST SRM 915b and NIST 290
SRM 1486. Geostandards and Geoanalytical Research 32, 311 (2008).
291
54. A. D. Jacobson, M. G. Andrews, G. O. Lehn, C. Holmden, Silicate versus carbonate weathering in 292
Iceland: New insights from Ca isotopes. Earth Planet. Sci. Lett. 416, 132 (2015).
293
55. D. L. Parkhurst, C. A. J. Appelo, “User's guide to PHREEQC (Version 2): A computer program for 294
speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations.” USGS 295
Water-Resources Investigations Report 99-4259 (1999).
296
56. S. Das Sharma, D. J. Patil, K. Gopalan, Temperature dependence of oxygen isotope fractionation of 297
CO2 from magnesite-phosphoric acid reaction. Geochim. Cosmochim. Acta 66, 589 (2002).
298
57. E. D. Young, A. Galy, The isotope geochemistry and cosmochemistry of magnesium. Rev. Mineral.
299
Geochem. 55, 197 (2004).
300
58. J. A. Karhu, H. D. Holland, Carbon isotopes and the rise of atmospheric oxygen. Geology 24, 867 301
(1996).
302
59. A. B. Kuznetsov et al., Sr isotopic composition of Paleoproterozoic 13C-rich carbonate rocks: The 303
Tulomozero Formation, SE Fennoscandian Shield. Precambrian Research 182, 300 (2010).
304
60. V. A. Melezhik, A. E. Fallick, P. V. Medvedev, V. V. Makarikhin, Palaeoproterozoic magnesite:
305
Lithological and isotopic evidence for playa/sabkha environments. Sedimentology 48, 379 (2001).
306
61. J. K. Warren, Evaporites through time: Tectonic, climatic and eustatic controls in marine and 307
nonmarine deposits. Earth-Science Reviews 98, 217 (2010).
308
13
62. C. E. Harvie, J. H. Weare, The prediction of mineral solubilities in natural waters: The Na–K–Mg–
309
Ca–Cl–SO4–H2O system from zero to high concentration at 25°C. Geochim. Cosmochim. Acta 44, 310
981 (1980).
311
63. R. J. Spencer, L. A. Hardie, Control of seawater composition by mixing of river waters and mid- 312
ocean ridge hydrothermal brines. Fluid-mineral interactions: A tribute to HP Eugster: Geochemical 313
Society Special Publication 19, 409 (1990).
314
64. T. K. Lowenstein, M. N. Timofeeff, V. M. Kovalevych, J. Horita, The major-ion composition of 315
Permian seawater. Geochim. Cosmochim. Acta 69, 1701 (2005).
316
65. T. K. Lowenstein, Origin of depositional cycles in a Permian "saline giant": The Salado (McNutt 317
zone) evaporites of New Mexico and Texas. Geol. Soc. Am. Bull. 100, 592 (1988).
318
66. R. Langbein, “The Zechstein sulphates: the state of the art” in The Zechstein Facies in Europe, T. M.
319
Peryt, Ed. (Springer-Verlag, 1987), pp. 143-188.
320
67. V. A. Melezhik, A. E. Fallick, P. V. Medvedev, V. V. Makarikhin, Extreme 13Ccarb enrichment in ca.
321
2.0 Ga magnesite–stromatolite–dolomite–'red beds' association in a global context: A case for the 322
world-wide signal enhanced by a local environment. Earth-Science Reviews 48, 71 (1999).
323
68. C. E. Harvie, N. Møller, J. H. Weare, The prediction of mineral solubilities in natural waters: The Na- 324
K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O system to high ionic strengths at 25°C. Geochim.
325
Cosmochim. Acta 48, 723 (1984).
326
69. C. M. Bethke, “The Geochemist's Workbench® Release 7.0 (four volumes)” (Hydrogeology 327
Program, University of Illinois, Urbana, Illinois, 2007).
328
70. J. Huang, Y. Xiao, Mg-Sr isotopes of low-δ26Mg basalts tracing recycled carbonate species:
329
Implication for the initial melting depth of the carbonated mantle in Eastern China. International 330
Geology Review 58, 1350 (2016).
331
71. J. Schott, V. Mavromatis, T. Fujii, C. R. Pearce, E. H. Oelkers, The control of carbonate mineral Mg 332
isotope composition by aqueous speciation: Theoretical and experimental modeling. Chem. Geol.
333
445, 120 (2016).
334
72. T. M. Hensley, “Calcium isotopic variation in marine evaporites and carbonates: Applications to Late 335
Miocene Mediterranean brine chemistry and Late Cenozoic calcium cycling in the oceans”, thesis, 336
University of California, San Diego (2006).
337
73. N. Gussone et al., Calcium isotope fractionation in calcite and aragonite. Geochim. Cosmochim. Acta 338
69, 4485 (2005).
339
74. W. T. Holser, “Mineralogy of evaporites” in Marine minerals, R. G. Burns, Ed. (Mineralogical 340
Society of America, 1979), vol. 6, pp. 211–294.
341
75. C. L. Blättler, G. M. Henderson, H. C. Jenkyns, Explaining the Phanerozoic Ca isotope history of 342
seawater. Geology 40, 843 (2012).
343
76. J. Horita, T. J. Friedman, B. Lazar, H. D. Holland, The composition of Permian seawater. Geochim.
344
Cosmochim. Acta 55, 417 (1991).
345
14
77. R. J. Spencer, Sulfate minerals in evaporite deposits. Rev. Mineral. Geochem. 40, 173 (2000).
346
78. I. Halevy, A. Bachan, The geologic history of seawater pH. Science 355, 1069 (2017).
347
79. D. A. Fike, A. S. Bradley, C. V. Rose, Rethinking the ancient sulfur cycle. Annu. Rev. Earth Planet.
348
Sci. Lett. 43, 593 (2015).
349
350
Acknowledgements:
351
Data presented in this study are available in the Supplementary Materials. Funding sources:
352
Simons Foundation (SCOL 339006 to C.L.B.), European Research Council (ERC Horizon 2020
353
grant 678812 to M.C.), Research Council of Norway (RCN Centres of Excellence funding
354
scheme project 223259 to K.P. and A.L.), Estonian Science Agency (PUT696 to K.K., A.L.,
355
K.P., T.K.). Thanks to T.H. Bui for assistance with sulfur isotope analyses and B.A. Wing for
356
helpful discussions. Core material from the Onega Parametric Hole is maintained by the Institute
357
of Geology, Karelian Science Centre, Petrozavodsk.
358
359
Fig. 1. Representative Tulomozero Formation evaporite rocks in the OPH. A. Halite cored
360
intervals (box length is 1 m; 2900 m depth); note rounded gray intraclasts consisting of
361
mudstone, anhydrite and magnesite. B. Halite with felted anhydrite grains and anhydrite-
362
magnesite intraclasts (2854 m depth). C. Magnesite-anhydrite (white) and halite overlain sharply
363
by laminated red-gray mudstone with desiccation cracks; magnesite-anhydrite bed infilling
364
compacted desiccation crack at top of image (2528 m depth). D. Laminated fine-grained
365
dolostone (2304 m depth). White bars are 1 cm in length.
366
367
15
Fig. 2. Interpretive stratigraphy of the Tulomozero Formation in the OPH and associated
368
mineralogical and geochemical data. Calcium isotope data for samples influenced by former
369
aragonite shown by open symbols with crosses. Additional δ34S from (30); cc–calcite, arag–
370
aragonite, dol–dolomite. Methods described in supplementary materials.
371
372
Fig. 3. Constraints on seawater chemistry during deposition of the Tulomozero Formation.
373
Circles show batch evaporation simulations with variable calcium and sulfate concentrations; all
374
other ions as in modern seawater. Small, empty circles indicate failure to precipitate gypsum
375
prior to halite. Color of filled circles indicates fraction of initial calcium (fCa) removed at halite
376
saturation. Conversion of fCa to δ44/40Ca range is based on a Rayleigh distillation model with α =
377
0.99905 (see supplementary materials). M and S identify composition of modern and estimated
378
Silurian seawater (8), respectively. Lines indicate constraints from OPH observations (arrows
379
give directionality), and the blue region shows the range of seawater compositions consistent
380
with these constraints.
381
Tulomozero FormationArchean gneissZaonega and younger Fms
dark gray = cored intervals of Tulomozero Formation OPH
2 km
3 km Unit 1 halite-rich interval:
restricted marine embayment Unit 2 anhydrite, magnesite, mudstone with evaporite pseudomorphs:
sabkha-marine embayment Unit 3
dolostone, magnesite, mudstone with evaporite fabrics:
sabkha-coastal plain to shallow- marine
2500
2600 2200
2300
2700
2800
2900
0 50 100 50 100
depth (m)
evaporite mineralogy (%)
insoluble fraction (%)
analytical error:
1.5 1.0 0.5 0 44/40Ca (‰ SW)
0 5 10
34S (‰ VCDT)
0.2 0 0.2
33S (‰)
2200
2300
2400
2500
2600
2700
2800
2900 17O (‰)
depth (m)
carbonates cc+arag+dol magnesite Ca sulfates halite Mg-K sulfates
this study Krupenik et al., 2011
(2)
0.2 0.1 0
0 5 10 15 20 25 30 0
5 10 15 20 25 30 35
Gypsum
bef ore halit
e
Mg sulf
ates occur Ca
isot ope
const raint
SO4(mmol/kg)
Ca (mmol/kg)
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
0.2 0.4 0.6 0.8 1.0 1.5 2.0
0.0 44/40Ca range (‰) fCa
M S
1
Supplementary Materials for
Two-billion-year-old Evaporites Capture Earth’s Great Oxidation
C.L. Blättler, M.W. Claire, A.R. Prave, K. Kirsimäe, J.A. Higgins, P.V. Medvedev, A.E.
Romashkin, D.V. Rychanchik, A.L. Zerkle, K. Paiste, T. Kreitsmann, I.L. Millar, J.A.
Hayles, H. Bao, A.V. Turchyn, M.R. Warke, A. Lepland correspondence to: blattler@princeton.edu
This PDF file includes:
Materials and Methods Supplementary Text Figs. S1 to S14 Table S1
Other Supplementary Materials for this manuscript include the following:
Data tables S1 to S2
2 Materials and Methods
Geological context
The Paleoproterozoic sedimentary-volcanic succession in the Onega Basin (Fig. S1) is moderately deformed into open folds cut by high-angle faults and metamorphosed to lower greenschist facies during the 1.85 Ga Svecofennian orogeny (31). In the Onega Parametric Hole (OPH), the Tulomozero Formation sits above Archean gneiss and below the conformably overlying Zaonega Formation, composed of siliceous dolomarl-
rhythmite and sulfidic fine-grained siliciclastic and carbonaceous rocks. The lower contact with Archean rocks was not recovered, but downhole geophysical data were used to identify contacts in intervals without core recovery and indicate that the halite-
dominated unit sits directly on gneiss (32). Seismic reflection data (33) show that the Tulomozero Formation halite and associated salts form wedge-shaped bodies that thicken west-northwest (from ca. 150 to 400 m over 2 km) and that reflectors have uniformly low-angle dips, indicating that the evaporite units have not undergone significant mass flowage.
Although robust 3-D control for the units dominated by halite and anhydrite- magnesite are lacking (they are only known from the OPH core), the lithofacies
associations displayed in the Tulomozero Formation across the Onega Basin indicate that the major evaporite facies was fringed by sabkha-coastal plains and stromatolite reefs (observed in outcrop, e.g. (34); Fig. S2). Additional sedimentological features observed in the Tulomozero Formation in the OPH that support a shallow, tidally influenced marine sabkha environment include common desiccation cracks and quartz- and
dolomite-pseudomorphed sulfate minerals in the form of laths, swallow-tail-terminating crystals, cauliflower nodules and lenses, and chicken-wire fabric, as well as tidal bundles (cm-scale rhythmic thickening and thinning of layers) and combined-flow ripples with reactivation surfaces and mud drapes. The stromatolite reefs may have acted as barriers restricting exchange of water between the embayment and ocean. Garabogazköl Bay, Turkmenistan offers a partial modern analogue to envisage such a paleogeography (Fig.
S3). There, a slender bedrock ridge cut by a narrow inlet separates the Bay from the Caspian Sea; the inlet allows an influx from the Sea with salinities of 1–3% into the Bay where evaporation then raises concentrations to 30–35% in the salt-precipitating areas distal to the inlet (35, 36).
Age constraints
The age of the Tulomozero Formation is not precisely constrained. Carbonate carbon isotope ratios in the OPH from this study (see below) and previous work (30), as well as other cores in the Onega Basin (18, 19, 37), show that the Tulomozero Formation records heavy δ13C values, in line with the Paleoproterozoic Lomagundi-Jatuli positive carbonate δ13C excursion. The formation is older than 1956 ± 5 Ma, the age of igneous dykes (38) that intrude the conformably overlying Zaonega Formation, but a maximum age for deposition remains debated: a Pb-Pb age of 2090 ± 70 Ma has been obtained from Tulomozero Formation carbonates (39), but zircons from a putative lava underlying the Formation yielded U-Pb ages of 1976 ± 9 Ma (40) and 1975.3 ± 2.8 Ma (41). Thus, until resolved, the most judicious depositional age to consider for the Tulomozero Formation is ca. 2.0 Ga.
3 X-ray diffractometry
Fifty-nine samples were taken for X-ray diffractometry (XRD); these were pulverized by hand with an agate mortar and pestle, and unoriented preparations were made on low-background sample holders of Si single-crystal wafers with high-index surface orientation. Preparations were measured on a Bruker D8 Advance using a Johannson-type Vario 1 focusing primary monochromator providing pure CuKα1 radiation to minimize peak overlaps. Scans were recorded with a LynxEye detector in the range 2θ = 2–70°, with a step size of 0.012° and counting time of 1 second per step.
Mineralogical compositions (Data table S1) were modeled and interpreted using the Rietveld algorithm-based codes Topaz and Siroquant 3.0 (42); the relative error of quantification is better than 10% for major phases (>5 wt%) and better than 20% for minor phases (<5 wt%).
Triple oxygen isotope and sulfur isotope analysis
Triple oxygen isotope analyses of sulfate were performed at Louisiana State University. Raw samples were treated with 3 M HCl for the extraction of acid-soluble sulfate (gypsum, anhydrite, carbonate-associated sulfate). To the resulting solution, ~4 mL of BaCl2-saturated 3 M HCl was added, resulting in the quantitative precipitation of BaSO4 which was then washed in distilled and deionized water until the washing solution was pH-neutral. The resulting BaSO4 powders were purified by repeated diethylene- triaminepentaacetic acid (DTPA) dissolution and reprecipitation (43). This purification procedure is critical for the removal of nitrate. Purified BaSO4 samples were converted to O2 by heating with a CO2-laser under BrF5. The O2 from laser fluorination was
analyzed for its triple oxygen isotope composition on a Thermo Finnigan MAT 253 (44).
All reported Δ17O values are the average of triplicate measurements and were computed using the logarithmic delta-prime notation with λ = 0.52 to facilitate comparison with published data: Δ17O = δʹ17O – 0.52 · δʹ18O, where δʹxO = ln(δxO/1000 + 1) · 1000. The precision on Δ17O values is better than ± 0.05‰.
The triple oxygen isotope composition is additionally reported using λ = 0.5305 (Data table S2). This alternate reference frame (45) is considered more appropriate for small magnitude Δ17O anomalies. The correction is performed using the more accurate δʹ18O from TCEA results (high temperature conversion elemental analyzer): Δ17O0.5305 = (0.52 · δʹ18OTCEA + Δ17O0.52) – 0.5305 · δʹ18OTCEA. Following this correction, the OPH samples average Δ17O0.5305 = –0.19 ± 0.03‰ (1 s.d.). Such small magnitude Δ17O values can have several origins, including but not restricted to stratospheric photochemical effects involving O2 and O3, which are in turn functions of the photosynthetic oxygen flux, oxygen consumption via respiration and weathering, the total atmospheric oxygen reservoir size (pO2), and its relative size to pCO2 (20). Additional processes such as low- temperature equilibrium fractionation and microbial sulfate cycling are also possible sources of the Δ17O signal.
Multiple sulfur isotope analyses were conducted on these same purified sulfate samples at McGill University, following an established methodology (46). Sulfate samples were purified with a chromium reduction procedure and converted to H2S gas through reaction with a Thode reduction solution (47, 48). The sulfide gas was captured as ZnS, converted to AgS2, and finally converted to SF6 gas through reaction with excess
4 F2. The resulting SF6 was purified cryogenically and chromatographically, and its
isotopic composition was quantified via gas source mass spectrometry on a Thermo Finnigan MAT 253. Long-term reproducibility (2σ) for δ34S values is better than
± 0.2‰, and for Δ33S and Δ36S values is ± 0.02‰ and ± 0.4‰, respectively. Triple- isotope and Δ-Δ plots (Fig. S4) demonstrate the relationships among the multiple sulfur isotopes.
Additional sulfate δ34S measurements were performed at the University of Cambridge in the Godwin Laboratory, following previously described methods (49).
Samples were dissolved in deionized water over 4–5 hours and sulfate was precipitated from the supernatant as barite by adding a saturated barium chloride solution. The barite precipitate was rinsed with deionized water, then dried and combusted at 1030 °C in a Flash Element Analyzer (EA) with the resulting sulfur dioxide (SO2) analyzed by continuous flow gas source isotope ratio mass spectrometry (Thermo Delta V Plus) and δ34S values corrected to NBS 127 = 20.3‰. The error for these δ34Ssulfate analyses is 0.3‰ based on replicate analyses of standards.
Calcium isotope analysis
Calcium isotope ratios were measured at Princeton University on fractions of the powdered samples used for XRD. The bulk powders were dissolved in 0.1 M buffered acetic acid–ammonium hydroxide solution (pH of ~4.5) for four hours. The soluble phase was then diluted and introduced into an automated ion chromatograph (Thermo- Dionex ICS-5000+) with a fraction collection system to isolate calcium for isotopic analysis (50). The isotope ratios therefore reflect the sum of calcium from all soluble calcium-bearing phases in a given sample, predominantly anhydrite and dolomite with minor gypsum, calcite, and aragonite (Data table S1).
Mass spectrometry was performed on a Thermo Neptune Plus multi-collector inductively coupled plasma mass spectrometer (ICP-MS) with an ESI Apex-IR sample introduction system (50). Corrections for isobaric strontium interferences were based on the double-charged 87Sr peak; polyatomic interferences (predominantly argon hydride, carbon dioxide, and nitrous oxide) were avoided with sufficient mass resolution. Sample- standard bracketing was used to calculate δ44/42Ca values relative to modern seawater. A triple-isotope plot shows the expected mass dependence among 44Ca, 43Ca, and 42Ca (Fig.
S5). Reported δ44/40Ca values are calculated following kinetic fractionation laws (51) and assuming no radiogenic 40Ca excess, which is consistent with previous work showing a lack of resolvable 40Ca excess in modern seawater or ancient marine precipitates (52);
alternatively, these can be considered initial δ44/40Ca values. External precision on δ44/40Ca values is ± 0.14‰ (2σ), derived from the long-term reproducibility of the
carbonate standard SRM-915b (n=155). The δ44/40Ca value obtained for SRM-915b is –1.15‰, which is indistinguishable from published values (53, 54).
Evaporation simulations
Evaporation simulations were conducted with the PHREEQC program (55) by the U.S. Geological Survey (https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/).
Simulations were initialized with varying concentrations of carbonate alkalinity, calcium, chloride, potassium, magnesium, sodium, sulfate, and trace strontium and set in
equilibrium with pCO2 = 10–1.5 atm (100× pre-industrial, see below for further
5 explanation) at a temperature of 25 °C with no charge balance imposed (i.e. net positive or negative charges were permitted and remained constant throughout a single
simulation). The minerals in Table S1 were permitted to precipitate if the respective saturations were attained. Note that calcite was the only carbonate mineral included;
dolomite and magnesite precipitation were suppressed because of known kinetic barriers (see below for discussion of magnesite). The pitzer.dat database provided with
PHREEQC was used for thermodynamic speciation and equilibrium calculations during the removal of over 99% of the initial water, and no back-reactions were allowed between fluid and solid phases.
Simulated evaporation of modern seawater yields a mineral assemblage of calcite, celestine, gypsum, glauberite, halite, polyhalite, bloedite, epsomite, kainite, kieserite, carnallite, and bischofite (listed in order of appearance with progressive evaporation).
Simulated evaporation of estimated Silurian seawater (8) yields a distinctly different assemblage, notably lacking magnesium sulfate minerals, of calcite, gypsum, halite, anhydrite, sylvite, carnallite, bischofite, and MgCl2 tetrahydrate. The constraints in Fig. 3 (main text) were established using simulations to determine the degree of evaporation (= initial water content / final water content) for the initial appearance of gypsum and halite, the presence/absence of select evaporite minerals, and the fraction of calcium removed at the onset of halite precipitation (see below for discussion of interpreting calcium isotopes and sensitivity tests of simulation behavior). The results using PHREEQC differ notably from simple calculations relying on stoichiometric
relationships among the major ions because of variable ion activities in concentrated brines.
Supplementary Text
Supporting geochemical data: Methods
Additional geochemical data, including carbonate carbon isotope ratios (δ13C), magnesium isotope ratios (δ26Mg), strontium isotope ratios (87Sr/86Sr), and trace element concentrations, provide additional context for the OPH (Data table S2). Carbon isotope ratios in calcium carbonate, dolomite, and magnesite were measured on micro-drilled powdered samples on a Thermo Delta V Advantage continuous flow isotope ratio mass spectrometer at Tartu University. Magnesite and dolomite-magnesite samples were reacted for 96 h at 70 °C to allow for the complete reaction of magnesite (method modified from (56)). Dolomite samples with <1 wt% magnesite content, as determined by XRD analysis, were reacted for 2.5 h at 70 °C. Carbon isotope ratios are reported relative to Vienna Pee Dee Belemnite (VPDB) with reproducibility (2σ) better than
± 0.1‰.
Magnesium isotopes and trace elements were measured at Princeton University on the same dissolved fractions used for calcium isotope analysis (see above). For
magnesium isotopes, automated ion chromatography (Thermo-Dionex ICS-5000+) was used to isolate magnesium, and isotopic ratios were then obtained by multi-collector ICP- MS (Thermo Neptune Plus) with wet aerosol introduction through a cyclonic spray chamber (50). Sample-standard bracketing in low resolution was used to calculate δ26Mg values relative to the standard DSM3. A triple-isotope plot shows the expected mass dependence among 24Mg, 25Mg, and 26Mg (Fig. S6). External precision on δ26Mg values is ± 0.09‰ (2σ), derived from the long-term reproducibility of the standard Cambridge-1
6 (n=76). The δ26Mg value obtained for Cambridge-1 is –2.56‰, which is
indistinguishable from published values (57).
Trace elements were analyzed on a quadrupole ICP-MS at Princeton University (Thermo iCAP). Samples were diluted so that major cations were no greater than 10 ppm, and scandium was added as an internal standard. External standard calibration curves bracketed the concentrations of each element of interest in the sample solutions, yielding precision of approximately ± 5%. Reported trace element concentrations reflect those of the acetic-acid-soluble component only (sum of carbonate, sulfate, and chloride evaporite minerals), not whole rock concentrations.
Strontium isotopes were analyzed at the NERC Isotope Geosciences Laboratory (British Geological Survey). The powdered samples were first leached in dilute (10%) acetic acid at 50 °C for 30 minutes to remove secondary carbonate phases. Residual material was washed repeatedly in water, then dissolved in distilled 6 M HCl at 50 °C for one hour. The dissolved material was removed after centrifuging, dried down, and converted to nitrate form using 1 mL 16 M HNO3. Strontium was separated using Sr- SPEC ion exchange resin, loaded onto outgassed single Re-filaments, and analyzed in a Thermo Scientific Triton thermal ionization mass spectrometer operated in multi-
dynamic mode with 100 measurement cycles for each sample. Eleven measurements of the strontium standard NBS-987 across the interval of analysis yield 87Sr/86Sr = 0.710256 (± 0.000007). Five seawater standards analyzed in the same period yield 87Sr/86Sr = 0.709182 (± 0.000004). The sample 87Sr/86Sri ratios have been corrected for the age of the Tulomozero Formation (2.0 Ga), and precision (2σ) for all samples is better than
± 0.000011.
Supporting geochemical data: Results
The additional geochemical datasets (Fig. S7 and S8) are generally consistent with the observations and interpretations presented in the main text. The magnesium isotope data are discussed below. Carbonate δ13C values from dolomite, magnesite, and rarely calcium carbonate range from +8.1 to +13.9‰, except for one sample in Unit 2 showing extreme 13C enrichment at +18.8‰. These values are characteristic of carbonates precipitated during the Lomagundi-Jatuli positive carbonate δ13C excursion (58).
Elevated strontium concentrations are found in three carbonate samples that also have very low δ44/40Ca values, one of which contains 4% aragonite (detected by XRD).
The combined geochemical evidence suggests that aragonite was likely the precursor carbonate mineral for these three samples, and that they were partially/completely
converted to calcite, dolomite, and magnesite. The sample with detectable aragonite also contains the highest proportion of calcite (55%) and highest strontium content (690 ppm) of any OPH sample, indicating relatively closed system (sediment-buffered) conditions during early marine neomorphism and recrystallization.
Relatively high strontium is also found in association with the anhydrite-magnesite of Unit 2, and to a lesser degree in Unit 1, where it scales with the abundance of calcium sulfate minerals (Fig. S9). This relationship suggests that where strontium is not
associated with aragonite, it may either be substituted into calcium sulfate minerals or exist as strontium sulfate (celestine), a trace phase predicted to co-precipitate with gypsum and anhydrite in evaporation simulations. Manganese is enriched in the