Available online at www.sciencedirect.com ScienceDirect

Geochimica et Cosmochimica Acta 209 (2017) 216–232 www.elsevier.com/locate/gca

Initiation of modern-style plate tectonics recorded in marine chemical

a,b, c a,b a,b Aaron M. Satkoski ⇑, Philip Fralick , Brian L. Beard , Clark M. Johnson

a University of Wisconsin-Madison, Department of Geoscience, 1215 West Dayton Street, Madison, WI 53706, United States b NASA Astrobiology Institute, United States c Department of , Lakehead University, Thunder Bay, ON PB7 5E1, Canada

Received 16 January 2017; accepted in revised form 14 April 2017; Available online 24 April 2017


The chemistry of the oceans in part reflects a balance between inputs from the and mantle. Traditionally, it has been thought that ocean chemistry was dominated by mantle sources, but recent work has suggested that continental during the Archean provided a much higher flux to the oceans than previously recognized. Here, we present new Rb-Sr and Sm-Nd isotope compositions on carbonate (dolomite and limestone) from the 2.94 Ga Red Lake and 2.80 Ga Steep Rock groups in the Superior Province, Canada to assess the potential impact continental weathering had on ocean chemistry during the Mesoarchean, a time when initiation of modern-style plate tectonics has been proposed to have occurred. The low Rb contents of all carbonate samples suggest that clastic contamination does not affect the Sr isotope compositions. Using O and Sr isotope modeling, we identified unaltered samples and estimate a 87Sr/86Sr ratio of 0.70173 for seawater at 2.94 Ga and 0.70182 at 2.80 Ga. Strontium isotope compositions from both Red Lake and Steep Rock indicate that seawater was significantly more radiogenic than contemporaneous mantle, and suggests that weathering of evolved continental crust was an important input to seawater. Continental weathering likely affected seawater chemistry through uplift of continental lithosphere during the initiation of modern-style plate tectonics at 3.2 Ga, a model that is contrary to those that suggest the Archean continents were small in extent and largely submerged. Initiation of modern-style plate tectonics and associated con- tinental weathering had an important effect on the biosphere, including increased nutrient delivery, as well as creation of eco- logical niches that allowed development of the first biologically produced shallow marine redox gradients. Ó 2017 Elsevier Ltd. All rights reserved.

Keywords: Plate tectonics; Archean; Carbonate; Strontium

1. INTRODUCTION evolved continental crust, of substantial thickness such that it was emergent, seems likely to have been tied to orogenesis The interplay between continental crustal and (e.g., Kemp and Hawkesworth, 2003). Multiple lines of evi- the beginning of plate tectonics has been a long-standing dence suggest that modern-style plate tectonics probably question in science (Condie and Pease, 2008, and ref- began to operate by 3.2–3.0 Ga (Shirey and Richardson, erences therein). Stabilization of significant quantities of 2011; Næraa et al., 2012; Dhuime et al., 2012, 2015; Satkoski et al., 2013), and if this coincided with stabiliza- tion of evolved, thick continental crust, particularly under conditions of high CO2 levels in the Archean atmosphere ⇑ Corresponding author at: University of Wisconsin-Madison, Department of Geoscience, 1215 West Dayton Street, Madison, WI (e.g., Kasting, 2010), there should have been a marked 53706, United States. increase in continental weathering fluxes to the oceans. Sup- E-mail address: [email protected] (A.M. Satkoski). port for this model is supplied by studies of clastic http://dx.doi.org/10.1016/j.gca.2017.04.024 0016-7037/Ó 2017 Elsevier Ltd. All rights reserved. A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 217 sequences suggesting that weathering conditions through- out the Archean were extreme (Hessler and Lowe, 2006 and references therein). Of the radiogenic isotope systems that respond to conti- nental crustal evolution, only the 87Rb-87Sr system is asso- ciated with an element that has a long residence time in the oceans, relative to ocean mixing timescales (Broecker and Peng, 1982; Hodell et al., 1990), and hence is amenable to tracing using marine chemical sediments. Historically, there is a strongly held notion that Archean ocean chemistry was almost entirely buffered by oceanic hydrothermal fluid cir- culation, with little influence from continental input (e.g., Shields, 2007; Pons et al., 2013). In part this reflects the influence of models that suggest continental crust was lar- gely submerged in the early to middle Archean (e.g., Shields and Veizer, 2002; Flament et al., 2013). In addition, however, Sr isotope studies of carbonates as a proxy for Archean ocean chemistry have been influential in inferring a large mantle contribution, but these have been hindered by the lack of pristine chemical sedimentary rocks in the geologic record throughout much of the Archean (e.g., Kamber, 2010). The wide scatter in 87Sr/86Sr ratios in Archean carbonates (Shields and Veizer, 2002) may reflect issues of siliciclastic contamination, diagenesis, and alter- ation by fluid-rock interaction. This scatter has led to the widespread assumption that the lowest 87Sr/86Sr ratios must be closest to seawater, and yet it can be shown that this assumption can lead to an under-estimate of seawater 87Sr/86Sr ratios if a local hydrothermal component was pre- sent during mineral formation or during post-depositional alteration. For example, using the relatively insoluble min- eral barite, which is likely resistant to post-depositional alteration, Satkoski et al. (2016) showed that the global ocean at 3.2 Ga had a significantly more radiogenic Sr iso- tope composition than contemporaneous mantle, suggest- ing that continental weathering affected seawater chemistry at that time. Although such an approach is promising, barite is only sporadically preserved in the Archean rock record, and hence provides limited snapshots of Archean ocean chemistry. In this contribution, we take a new look at the Sr isotope record for Archean carbonates through study of two excep- tionally well-preserved Mesoarchean sequences of 2.94 and 2.80 Ga age at Red Lake and Steep Rock in Superior Pro- vince, Canada. Combining Sr isotope data with REE and Nd isotope data, as well as stable isotope data, we can assess issues of siliciclastic contamination and post- depositional alteration, allowing us to isolate the seawater signals. The goal is to evaluate the potential impact of con- tinental weathering on ocean chemistry in a time period that immediately follows the proposed onset of modern- style plate tectonics at 3.2–3.0 Ga (e.g., Dhuime et al., 2015).

2. GEOLOGICAL BACKGROUND AND SAMPLES Fig. 1. Location maps for Red Lake and Steep Rock samples. (C) The numbered dots (1–3) show sampling sites for Red Lake 2.1. Red Lake carbonate – 2.94 Ga samples. (D) The numbered dots (4–6) show sampling sites for Steep Rock samples. The numbered dots (sampling locations) are The Red Lake carbonate platform is part of the Red linked to specific samples in Appendix 1. Lake in the Uchi Subprovince of the larger 218 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Superior Province, Canada (Tomlinson et al., 1998) Steep Rock carbonate has an age between 2801 and (Fig. 1). The Ball assemblage of the Red Lake greenstone 2780 Ma. We use 2800 Ma to represent the age of the car- belt is composed of mafic (pillowed) and felsic extrusive bonate samples analyzed as part of this study. The regional igneous rocks, iron formation, dolomite- beds, and metamorphic grade is lower facies (Wilks and stromatolitic carbonate (Pirie, 1981; Hofmann et al., Misbet, 1985), however, work by Veizer et al. (1982), 1985). The depositional environment for carbonate ranged Riding et al. (2014) and Fralick and Riding (2015) suggests from shallow marine, as highlighted by mounded stromato- that Steep Rock carbonate largely preserves its primary lites associated with flat-pebble breccias, to deeper marine chemistry. Detailed sample information can be found in for massive, slide brecciated and slumped carbonate layers Appendix 1. Based on texture, carbonate samples at Steep associated with iron formation and chert (McIntyre, 2014). Rock are divided into the following categories, (1) crystal The age of the is constrained by a rhyolitic fan, (2) void filling cement and (3) . Samples tuff(2940 + 2.4/ 1.7 Ma) below and a 2925 analyzed as part of this study were also studied by Riding À + 3.4/ 2.9 Ma rhyolitic flow above (Corfu and Wallace, et al. (2014) and Fralick and Riding (2015). À 1986). We use 2940 Ma to represent the age of the carbon- ate samples analyzed as part of this study as the sedimen- 3. METHODS tary rocks appear to lie conformably on the volcanic rocks. The regional metamorphic grade ranges from green- 3.1. Rb-Sr isotopes schist to lower amphibolite facies. Stromatolites in the Red Lake area include columnar, stratiform, and mounded mor- Approximately 10–40 mg of powder was spiked with an phologies (Hofmann et al., 1985; McIntyre, 2014), which enriched 87Rb-84Sr tracer solution and dissolved with 4 M suggests a great diversity in energy levels of the environ- HNO3 in clean Teflon on a hotplate for a minimum of ments and microbial at this time. Crystal fans occur 24 h. Strontium was separated from matrix elements Ò as layers interbedded with the stromatolites and as mounds, (including Rb) using Sr Spec resin in 4 M HNO3 (Beard which probably occupied areas further offshore (McIntyre, et al., 2013). Strontium was stripped from the resin using 2014). The crystals grew directly on the seafloor, reaching 2% HNO3. The resin was cleaned prior to use following 20 cm in height, and by analogy with Archean crystal fans using the method of Charlier et al. (2006). Rubidium was at other locations, were probably originally aragonite further separated from matrix elements using cation- (Sumner and Grotzinger, 2000). Detailed sample informa- exchange chromatography (BioRad AG-MP-50 100–200 tion can be found in Appendix 1. mesh) in 2.5 M HCl (Beard et al., 2013). Strontium samples were loaded onto Ta filaments with H3PO4 and isotopic 2.2. Steep Rock carbonate – 2.80 Ga analyses were made using a dynamic multi-collector routine on a VG Instruments Sector 54 thermal ionization mass The Steep Rock carbonate platform sampled here is part spectrometer (TIMS) at the University of Wisconsin- 88 11 of the larger 2.8 Ga Steep Rock Group of the central Wabi- Madison. Typical Sr ion intensities were 3 10À amps  goon Subprovince in the (Riding et al., and the Sr isotope ratios were exponentially normalized 2014; Fralick and Riding, 2015)(Fig 1). The Steep Rock using 86Sr/88Sr = 0.1194; this procedure removes any Group, at Steep Rock, contains a basal conglomerate mass-dependent Sr isotope variations that might exist in deposited in paleo-valleys incised into 3.0 Ga tonalite, the samples. Measurements of NIST SRM-987, E&A Sr  87 86 as well as a 500 m thick carbonate unit (Mosher Carbonate; and EN-1 Sr yielded an average Sr/ Sr ratio of focus of this study) that rests on the basal conglomerate or 0.710270 ± 0.000016 (n = 35, 2-SD), 0.708055 ± 0.000024 directly on tonalitic basement (Fralick et al., 2008). The (n = 8, 2-SD) and 0.709194 ± 0.000013 (n = 8, 2-SD) entire package represents a transgressive sequence that respectively. These uncertainties are comparable to those was deposited primarily as Ca-carbonate, with lesser estimated for initial 87Sr/86Sr ratios of ± 0.00002 based on amounts of dolomite, ankerite, and siderite (Riding et al., Sr and Rb (see below) isotope dilution analyses. Measured 2014). Laterally, the carbonate transitions into deep-water Sr isotope ratios are not normalized to any external stan- iron formation, and therefore the Steep Rock carbonate dard. Total procedural Sr blanks of 107, 15, and 16 pg were represents a platform that was constructed next to a deeper measured, which are negligible relative to samples. portion of the basin. Rubidium samples were loaded onto Ta filaments with The tonalitic basement underlying the Steep Rock car- H3PO4. Samples were analyzed in static-mode on Faraday bonate has an age of 3001.6 ± 1.7 Ma (Tomlinson et al., cups or using single Daly detector peak-hopping mode on 2003) and pyroclastics of the Dismal Ashrock overlying it a TIMS. Standard NIST SRM-984 (87Rb/85Rb = 0.386 are 2780.4 ± 1.4 Ma (Tomlinson et al., 2003), which repre- ± 0.002) was analyzed by both methods and produced an sents the minimum age of the Steep Rock Group, at Steep average 87Rb/85Rb = 0.3863 ± 0.0042 (n = 8, 2-SD) in static Rock Lake. The youngest detrital found in the basal mode and an average 87Rb/85Rb = 0.3804 ± 0.0020 (n = 9, conglomerate is 2779 ± 22 Ma, requiring the limestone to 2-SD) for single Daly detector peak-hopping. Relative to probably be younger than 2801 Ma. A correlative unit to NIST SRM-984, these measured ratios produce exponential the northeast (Lumby Lake belt) overlies volcanic rocks b factors of -0.0304 and 0.6247, respectively, which were with an age of 2828 ± 1 Ma. Using the age of the overlying applied to the measured 87Rb/85Rb ratios obtained on sam- volcanic rocks and the youngest detrital zircon in the basal ples. Total procedural Rb blanks of 58, 0.61 and 1.49 pg conglomerate, Fralick and Riding (2015) suggest that the were measured, which are negligible relative to samples. A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 219

3.2. Sm-Nd isotopes These uncertainties are comparable to those estimated for initial eNd(T) values at ±0.5 e-units based on Nd and Sm Approximately 500 mg of powdered carbonate was isotope dilution analyses. Measured Nd isotope ratios are spiked with an enriched 149Sm-150Nd tracer solution and not normalized to any external standard. A total procedu- dissolved with 4 M HNO3 in clean Teflon on a hotplate ral Nd blank of 20 pg was measured, which is negligible rel- for a minimum of 24 h. The bulk Rare Earth Elements ative to samples. (REEs) were separated from other matrix elements using Samarium was loaded onto single Re filaments with Si- cation-exchange chromatography (BioRad AG 50W X8 gel and H3PO4 and analyzed by TIMS in static mode as a 200–400 mesh resin) and 2.5 M HCl. The REEs are stripped metal (Sm+) using Faraday cups. Typical ion beams were 12 149 147 from the resin using 6 M HCl. The REEs were then sepa- 5 10À amps. The measured Sm/ Sm ratio was  147 152 rated from themselves using a 2-methylactic acid procedure exponentially corrected using Sm/ Sm = 0.5608. A (Johnson and Thompson, 1991). Neodymium was loaded total procedural Sm blank of 6 pg was measured, which is onto single Re filaments with Si-gel and H3PO4 and ana- negligible relative to samples. lyzed as NdO+ using a dynamic multi-collector routine on a VG Instruments Sector 54 TIMS at the University of 4. RESULTS Wisconsin-Madison. To enhance formation of NdO+ the pressure in the source was raised using an O2 gas bleed. 4.1. Oxygen isotopes Analyses were run with a 144Nd16O+ ion signal of between 11 0.74 and 1 10À amps. Measured data were exponen- Oxygen isotope compositions were determined on the  146 144 tially corrected using Nd/ Nd = 0.7219. Measurements same powders as used for Rb-Sr and Sm-Nd analyses, of UW AMES I, II, LaJolla and JNdi-1 standards yielded and are reported in McIntyre (2014) and Fralick and 143Nd/144Nd ratios of 0.512135 ± 0.000014 (n = 9, 2-SD), Riding (2015). The Red Lake samples have a range of 0.511968 ± 0.000012 (n = 4, 2-SD), 0.511849 ± 0.000022 d18O values from 14.7‰ to 20.3‰ (Fig. 2c). Dolomite sam- (n = 10, 2-SD), and 0.512106 ± 0.000021 (n = 6, 2-SD). ples have d18O values that range from 14.7‰ to 15.8‰. The The correlation between our measurements of LaJolla Ca-carbonate samples range from 18.0‰ to 20.3‰ and are and JNdi-1 is 1.000507, which lies within the range systematically higher than the dolomite samples. There is (1.000503 ± 15 1-SD) reported by Tanaka et al. (2000). no difference in d18O values between carbonate morpholo-

18 87 86 Fig. 2. MgO/CaO (A, D), Sr content in ppm (B, E) and d Ocarb (C, F) versus initial Sr/ Sr ratios for the Red Lake and Steep Rock samples, respectively. 220 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 gies. The Steep Rock samples have d18O values that range 2), which suggests the Sm-Nd relationship is not controlled from 19.7‰ to 22.4‰ (Fig. 2f), and there is no systematic by binary mixing (R2 = 0.0), and instead records small vari- difference in d18O values between carbonate morphologies. ations in initial Sm/Nd ratios and in situ growth.

4.2. Rb-Sr isotopes 4.4. REE + Y contents

Measured 87Rb/86Sr ratios for all samples are very low, REE and Y contents were measured on the same pow- less than 0.017, which results in very small corrections for ders as used for Rb-Sr and Sm-Nd analysis, and are in situ decay, <0.00071 for 87Sr/86Sr ratios. Although there reported in McIntyre (2014) and Fralick and Riding is a broad positive relation between measured 87Rb/86Sr (2015). Total REE contents for Red Lake samples range and 87Sr/86Sr ratios (Appendix 2), the slope of this correla- from 1.8–7.6 ppm and 1.6–19.7 ppm for Steep Rock sam- tion is far greater than that of the age of formation, indicat- ples. Elements with ‘‘SN” refer to those that have been nor- ing that 87Rb/86Sr-87Sr/86Sr relations do not record malized to PAAS (Post Archean Average Shale) values incorporation of a high Rb/Sr component (e.g., clays) at from Taylor and McLennan (1985). All samples have pos- * the time of formation. Data are discussed below in terms itive La anomalies, with Ce/Ce = CeSN/(0.5 LaSN + 0.5 87 86 of initial Sr/ Sr ratios, as calculated using the measured PrSN) less than unity; positive Eu anomalies, with Eu/ 87 86 87 86 * Rb/ Sr and Sr/ Sr ratios and age of formation. EuSN = EuSN/(2/3 SmSN + 1/3 TbSN) greater than unity; Initial 87Sr/86Sr ratios for Red Lake samples fall into and Y/Ho ratios greater than crustal rocks. two groups based on . Carbonate samples that have high MgO/CaO ratios (>0.1; dolomite), low Sr con- 5. DISCUSSION tents (<50 ppm), and low d18O values (<16‰) have high initial 87Sr/86Sr ratios that range from 0.70273 to 0.70337, 5.1. Assessing primary chemical signatures of carbonate whereas samples that have low MgO/CaO ratios (<0.1; cal- cite), high Sr contents (>200 ppm), and high d18O values The goal of the study is to constrain Archean seawater (>18‰) have significantly lower initial 87Sr/86Sr ratios that compositions as an indicator of continental weathering, range from 0.70173 to 0.70198 (Fig. 2a–c). but before this may be done, the effects of clastic contami- Similarly, initial 87Sr/86Sr ratios of Steep Rock samples nation and post-formation alteration must be determined. also vary based on mineralogy and composition (Fig. 2d– f). Samples that have MgO/CaO ratios <0.1 have initial 5.1.1. Clastic contamination 87Sr/86Sr ratios that vary from 0.70305 to 0.70182. Within Clastic contamination potentially affects Sr and Nd iso- the low MgO/CaO samples, initial 87Sr/86Sr ratios vary tope compositions of carbonates, as well as the REE + Y by carbonate morphology. On average, carbonates that contents. Here we used high-precision Rb contents, deter- have a crystal fan morphology have higher Sr contents, mined by isotope dilution, to quantify the extent of clastic consistent with the inferred aragonite precursor of the fans contamination, as Rb contents in carbonate are exclusively (Edwards et al., 2015 and references therein), and less derived from clastic erosion products. Carbonate that is radiogenic initial 87Sr/86Sr ratios than carbonate associated free of clastic material will typically have positive La and with microbial mats. Y anomalies in PAAS-normalized REE + Y patterns (Kamber and Webb, 2001), which may be quantified using * * 4.3. Sm-Nd isotopes Ce/Ce (Ce/Ce = CeSN/[0.5 LaSN + 0.5 PrSN]) and Y/Ho ratios, respectively. Relative enrichments in La and Y are Red Lake carbonates broadly define a due to the increased solution stability of these elements, 147Sm/144Nd-143Nd/144Nd errorchron with an apparent similar to the HREE (Nozaki et al., 1997; Alibo and age of 2596 ± 450 Ma and an eNd(2596) of -1.9 (initial Nozaki, 1999). The chondritic ratio for Y/Ho ranges from 143Nd/144Nd = 0.50918 ± 0.00037; MSWD = 15; Appendix 26 to 28 (Kamber and Webb, 2001), which is the range 2), which overlaps the accepted depositional age of these expected for bulk clastic contamination and crustal Ce/ samples (Section 2.1). No significant correlation exists on Ce* is expected to be unity. Siliciclastic contamination the same diagram if the samples are separated by MgO/ would therefore be expected to produce a negative correla- CaO ratios (calcite versus dolomite). The samples show tion between Y/Ho and Rb and a positive correlation no correlation (R2 = 0.1) when plotted on a 1/Nd versus between Ce/Ce* and Rb, which is not observed (Fig. 3). 143Nd/144Nd diagram (Appendix 2), which suggests the At the very low Rb contents measured in the samples, sili- Sm-Nd relations are not controlled by binary mixing, but ciclastic contamination is negligible, and this cannot record small differences in Sm/Nd ratios followed by explain the elevated initial 87Sr/86Sr ratios for some samples in situ growth since formation. (Fig. 3). Based on the trends in REEs and Sr isotopes, we Steep Rock carbonates similarly define a infer that siliciclastic contamination has not affected the 147Sm/144Nd-143Nd/144Nd errorchron with an apparent Sr and Nd isotope compositions. age of 2918 ± 410 Ma and an eNd(2918) of +2.6 (initial 143Nd/144Nd = 0.50898 ± 0.00029; MSWD = 2.9; Appen- 5.1.2. Fluid-rock alteration dix 2), which overlaps the accepted depositional age of these The high 87Sr/86Sr ratios for a number of carbonate samples (Section 2.2). Samples show no correlation when samples, if not due to clastic contamination (Fig. 3), may plotted on a 1/Nd versus 143Nd/144Nd diagram (Appendix reflect mixing with a high-87Sr/86Sr component after depo- A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 221

Fig. 3. LogRb concentration versus (A, D) Y/Ho, (B, E) Ce/Ce* and (C, F) initial 87Sr/86Sr for Red Lake and Steep Rock samples, respectively The arrow highlights the point along the mixing curve that represents 0.25% mixing of a PAAS composition and pristine carbonate. Full details of the modeling are in Appendix 3.

sition, as supported by the 87Rb/86Sr-87Sr/86Sr relations, with the interpretation that the crystal fans were originally which show data that plot far above an appropriate iso- aragonite, which has substantially higher Sr contents than chron relation (Appendix 2). Post-depositional mixing with calcite (Edwards et al., 2015 and references therein), and a high-87Sr/86Sr component would most likely occur hence more likely to retain its original Sr isotope composi- through fluid-rock interaction, which may be tested using tion during fluid-rock interaction. Variations in Sr/Ca O isotope compositions. Using the fluid-rock interaction ratios and d18O values can also identify primary carbonate model of McCulloch et al. (1981), samples that have high compositions (e.g., Edgar et al., 2015), and we interpret the initial 87Sr/86Sr ratios, which also tend to have low d18O Steep Rock samples with the highest Sr/Ca ratios and d18O values, can be explained by fluid-rock interaction over tem- values and lowest 87Sr/86Sr ratios to be closest to seawater peratures ranging from 100 to 300 °C(Fig. 4; Appendix 4). compositions (Appendix 5) and are shown as the least These results demonstrate that none of the dolomites from altered samples on Fig 4b. Red Lake (Fig. 4a) likely represent seawater Sr isotope Although O and Sr isotope compositions of the Red compositions; this is consistent with other studies that note Lake and Steep Rock carbonates have been affected by dolomite may not retain primary compositions (e.g., fluid-rock interaction, including dolomitization, it is unli- Kamber and Webb, 2001; Sena et al., 2014). Based on the kely these processes changed the REE + Y contents or modeling shown in Fig 4a, the crystal fan samples from Nd isotope compositions (e.g., Banner et al., 1988). The Red Lake are closest to seawater. For comparison of the high stability of REEs in carbonate reflects their substitu- crystal fans and microbialite from Steep Rock, the broadly tion for the Ca2+ ion in the calcite lattice (Zhong and lower d18O values and higher initial 87Sr/86Sr ratios of the Mucci, 1995) and, short of complete dissolution; the REEs microbialite samples also suggests that the crystal fans are will not be affected by later dolomitization (Webb and closer to seawater compositions (Fig. 4b). This is consistent Kamber, 2000). Recrystallization in the Red Lake and 222 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Fig. 4. Plots of d18O versus initial 87Sr/86Sr for (A) Red Lake and (B) Steep Rock samples. The dashed lines are modeled effects of water/rock interaction between an assumed pristine carbonate sample (d18O = 23‰, 87Sr/86Sr = 0.70173 for Red Lake and d18O = 23‰, 87Sr/86Sr = 0.70182 for Steep Rock). Full details of the modeling are in Appendix 4.

Steep Rock carbonates appears to be neomorphic HREE and has a positive La and Y anomaly (Ce/Ce* < 1, (McIntyre, 2014; Fralick and Riding, 2015), and thus the Y/Ho > 27). This LREE depletion has been observed in REE chemistry should be primary. Modeling fluid-rock the 3.45 Ga Strelley Pool carbonates, which suggests that interaction between seawater with a Nd content of the scavenging process has been operating since early in 6 2.8 10À ppm and carbonate with a Nd content 0.1 ppm the Archean (van Kranendonk et al., 2003). In modern sea- Â shows that the initial eNd(T) values only change by 0.01 water, REEs are also input into the oceans through oceanic 18 eNd units over the range of d Ocarb values observed as part hydrothermal fluids, but these are immediately scavenged of this study (Appendix 6). We therefore interpret all car- by Fe-Mn oxyhydroxides and therefore make a negligible bonates studied, including dolomite samples, to record pri- contribution to the overall seawater REE budget. In the mary REE compositions. Archean, however, when the heat flux from the interior of the Earth was greater and seawater contained low oxygen 5.2. The REE + Y characteristics of Archean seawater contents, REE budgets would have been dominated by hydrothermal REE, which would include a characteristic The primary REE input into the modern oceans is positive Eu anomaly and high Sm/Nd ratios (Danielson through continental weathering (Bau and Dulski, 1996). et al., 1992; Bau and Dulski, 1996). Alternatively, for car- Before entering the oceans, however, the LREE are scav- bonates formed in shallow-water, proximal settings where enged in estuaries due to adsorption onto the surfaces of continental input would be greater, we would expect high suspended or sinking particles (Webb and Kamber, 2000). Y/Ho ratios due to estuarine scavenging, decreased or no This scavenging results in a REE + Y pattern (PAAS nor- positive Eu anomalies, and low Sm/Nd ratios (Bau and malized) that is typically LREE depleted relative to the Dulski, 1996). A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 223

The two suites studied here appear to be distinct in terms the Sr isotope composition of the oceans is homogenous, of water mass mixing for the REE + Y. The range in REE even for restricted regions such as the Arctic Ocean + Y contents for the Red Lake carbonates suggests possible (Winter et al., 1997). It has been proposed, however, that mixtures between components (Fig. 5). The positive trend very low Sr contents measured in the 2.5 Ga Campbellrand * * observed between Sm/YbSN and Eu/EuSN (Eu/EuSN =- carbonates of South reflect a significantly lower sea- EuSN/[2/3 SmSN + 1/3 TbSN]), and negative trend between water Sr residence time in the (Kamber and Sm/YbSN and Y/Ho, suggest that seawater in the Red Lake Webb, 2001). In support of a short residence time for Sr, depositional environment represented a mixture of deep Kamber and Webb (2001) show a correlation (R2 = 0.72) (hydrothermal) and shallow water masses (Fig. 5c and d). between Sr and Nd isotopes in the 2.5 Ga Campbellrand The REE + Y patterns for Steep Rock carbonates, how- carbonates of South Africa. In modern seawater, the resi- ever, suggest a different interpretation (Fig. 6). When com- dence time of Nd is very short ( 500 yrs; Tachikawa  pared to a microbialite and carbonate from the et al., 2003), and assuming a similar residence time in the 3.45 Ga Strelley Pool carbonate occurrence, all Steep Rock late Archean, a reasonable assumption because Nd is not samples show similar features (Fig. 6a and b), which a redox-sensitive element, a correlation between Sr and include depleted LREE and positive La (Ce/Ce* < 1) and Nd isotopes could be taken as support for a short residence Y anomalies very similar to those observed for Red Lake time, and hence isotopic provinciality for Sr isotopes. samples. No clear correlations exist, however, between The eNd(T) values for the carbonates studied here vary * Sm/YbSN versus Eu/EuSN (Fig. 6c) and Y/Ho (Fig. 6d), significantly, particularly those of the Red Lake suite * suggesting they do not preserve mixing of deep (hydrother- (Appendix 7), and variations in Eu/EuSN, Y/Ho, and Sm/ mal) and shallow water masses. YbSN with Nd isotope compositions for Red Lake carbon- ates suggest mixing relations between a hydrothermal com- 5.3. Evolution of seawater Sr and Nd ponent and a shallow seawater component (Fig. 7). No such relations are seen in the data from Steep Rock In modern seawater, Sr has a residence time of 2.5 m.y., (Appendix 7), consistent with the lack of clear mixing rela- which is much longer than the ocean circulation time of tions among the REEs (Fig. 6), and the limited range in Nd 1500 yrs (Broecker and Peng, 1982), and hence at any time isotope compositions. It is difficult, however, to determine a

Fig. 5. (A) and (B) Normalized REE patterns for samples from Red Lake. For comparison a Holocene carbonate (Webb and Kamber, 2000) and carbonate (crystal fan and stromatolite) from the 3.45 Ga Strelley Pool (van Kranendonk et al., 2003). PAAS is post Archean average * shale from Taylor and McLennan (1985). (C) A plot of Sm/YbSN versus Eu/EuSN. Traditionally, the Eu anomaly is defined using Gd, however when samples have a Gd anomaly the Eu anomaly should be expressed as (Eu/[2/3Sm + 1/3 Tb] (Kamber and Webb, 2001). A * positive correlation between Sm/YbSN and Eu/EuSN could suggest mixing between shallow seawater and oceanic hydrothermal fluids. (D) A plot of Y/Ho versus Sm/YbSN. The observed negative trend can be explained by mixing of oceanic hydrothermal fluids and seawater. 224 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Fig. 6. Normalized REE patterns for samples from Steep Rock. For comparison a Holocene carbonate (Webb and Kamber, 2000) and carbonate (crystal fan and stromatolite) from the 3.45 Ga Strelley Pool (van Kranendonk et al., 2003). PAAS is post Archean average shale * from Taylor and McLennan (1985). (C) A plot of Sm/YbSN versus Eu/EuSN. Traditionally, the Eu anomaly is defined using Gd, however when samples have a Gd anomaly the Eu anomaly should be expressed as (Eu/[2/3Sm + 1/3 Tb] (Kamber and Webb, 2001). A positive * correlation between Sm/YbSN and Eu/EuSN could suggest mixing between shallow seawater and oceanic hydrothermal fluids, however, no definitive correlation exists. (D) A plot of Y/Ho versus Sm/YbSN. No definitive trend exists, suggesting that these samples do not preserve mixing between shallow seawater and deep oceanic hydrothermal fluids. clear relation between Sr and Nd isotope compositions for heterogeneity in seawater based on the Red Lake and Steep the Red Lake suite, despite the large range in initial Nd iso- Rock suites. We therefore take the inferred initial 87Sr/86Sr tope compositions because, as argued above, none of the ratios for Red Lake and Steep Rock carbonates, once the dolomites are considered to have preserved their primary effects of post-depositional changes are accounted for (iden- Sr isotope compositions. This leaves three samples of Ca- tified as the least altered samples of Fig 4 and listed in carbonate, and although there is a hint of a correlation of Table 1), to reflect ambient seawater in the Mesoarchean. 87 86 decreasing eNd(T) values and initial Sr/ Sr for these sam- ples (Fig. 7a), the range in initial 87Sr/86Sr is quite small. 5.4. Comparison to previously published Sr isotope data It is important to note that the Sr contents of the 2.5 Ga Campbellrand carbonates studied by Kamber and Webb The initial 87Sr/86Sr ratios (0.70173 and 0.70182) (2001) are exceptionally low, 10s of ppm Sr. The average inferred here for seawater at 2.94 and 2.80 Ga are within Sr content for Archean Ca-carbonate is 297 ppm, which is the range reported by Veizer et al. (1989) for 2.8 ± 0.2 Ga close to the average Sr content of the samples from this seawater (0.7025 ± 0.0015). Specifically, Veizer et al. study (328 ppm; Appendix 8). The high Sr contents from (1989) report 87Sr/86Sr ratios for Steep Rock with a range the majority of Archean carbonates imply that the from 0.70224 to 0.70178, with one sample at 0.70130. Com- 2.52 Ga carbonate studied by Kamber and Webb (2001) pared to the rest of the Archean carbonates (n = 39) pub- may not generally represent the chemistry of the Archean lished by Veizer et al. (1989), however, the least oceans. Whether the Sr residence time during the Archean radiogenic Sr sample (87Sr/86Sr = 0.70130) has a d18O value was truly similar to the REEs is still unknown, but the of 17.6‰, which is much lower than the average d18O value higher Sr contents of the samples measured in this study, measured in the current study. Therefore, the which is comparable to those of average Archean carbon- 87Sr/86Sr = 0.70130 carbonate likely represents a sample ate, suggests that there is no clear evidence for Sr isotope that underwent high-temperature water-rock interaction A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 225

87 86 * Fig. 7. Initial Nd isotope composition for Red Lake samples versus (A) initial Sr/ Sr (B) Sm/YbSN (C) Y/Ho and (D) Eu/EuSN. The arrows represent expected correlations to be observed if shallow seawater mixed with oceanic hydrothermal fluids. 226 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Table 1 6.1. Architecture of Archean crust Sr isotope compositions of Archean Reservoirs. Seawater Mantle 100 my 200 my 300 my Preserved Archean continental crust is primarily tona Age (Ga) 87Sr/86Sr 87Sr/86Srb Crustc Crustc Crustc lite-trondhjemite-granodiorite (TTG) in composition (e.g., 3.20 0.70139a 0.70075 0.70168 0.70261 0.70355 Moyen, 2011). Experimental results suggest that tonalite 2.94 0.70173 0.70095 0.70188 0.70281 0.70374 can be formed by directly melting a basaltic source 2.80 0.70182 0.70115 0.70208 0.70301 0.70394 (Foley, 2008), which is consistent with geologic studies using trace element and Hf isotopes (e.g., Satkoski et al., a Seawater isotope composition is from Satkoski et al. (2016). b Mantle isotope composition is from McCulloch (1994). 2013). Archean TTGs are unique in that they are typically c The Sr isotope composition of juvenile upper continental crust depleted in HREE, producing high La/YbCN ratios (‘‘CN” that was extracted from the mantle 100, 200 or 300 million is chondrite normalized), and show no EuCN anomaly (Eu/ prior to the age reported in the first column. The values are Eu* 1; Eu/Eu* = Eu /[pSm * Gd ]) when com- CN  CN CN CN CN calculated using mantle values from McCulloch (1994) and the Rb/ pared to Archean (e.g., Kleinhanns et al., 2003) Sr ratio of early Archean upper continental crust from Condie and and felsic igneous rocks (1993). (Moyen, 2011). These geochemical characteristics are ascribed to formation of TTG in the lower crust, at depths and hence does not reflect equilibrium with seawater. The > 30–40 km, below the plagioclase stability field and in  remaining samples from Veizer et al. (1989) have a range equilibrium with residual garnet (Kemp and in Sr and O isotope compositions that are similar to those Hawkesworth, 2003) and possibly amphibole (Rapp et al., reported here (Fig. 2). Of note, the large compilations of 2003), although amphibole does not by itself explain the Shields and Veizer (2002) and Prokoph et al. (2008) dis- very low HREE abundances. Extensive experimental stud- count the majority of Sr isotope compositions from Steep ies show that for mafic lithologies, dehydration melting in Rock of Veizer et al. (1989) in that they assume the lowest the lower crust breaks down plagioclase and produces gar- 87Sr/86Sr ratio most closely matches that of seawater, which net at depths of 35–45 km (e.g., Wolf and Wyllie, 1993; Sen is not always true (e.g., Satkoski et al., 2016). Importantly, and Dunn, 1994; Rapp and Watson, 1995), processes that * the seawater Sr isotope curve produced by Shields and can explain high La/YbCN ratios and Eu/EuCN 1. For-  Veizer (2002) is not defined by any actual data between mation of TTG at the base of thick basaltic crust is consis- 3.2 and 2.8 Ga, and essentially assumes mantle-like compo- tent with the extensive mantle melting that would be sitions are correct through this time interval. expected from the higher mantle potential temperatures in the Archean, relative to today (Davies, 1992). Using meta- 6. ARCHEAN CRUSTAL EVOLUTION – morphic petrology, a study of ten undisturbed Archean IMPLICATIONS FOR SEAWATER CHEMISTRY -greenstone belts indicates they formed in a crust that was, on average, 47 km thick, consistent with the The 87Sr/86Sr ratios inferred for seawater from this above arguments (Galer and Mezger, 1998). study, and that of Satkoski et al. (2016), indicate that sea- A compilation of REE data from Archean TTGs span- water was significantly more radiogenic than the mantle ning 3.9–2.5 Ga show no correlation between the La/YbCN 87 86 between 3.2 and 2.8 Ga (Table 1). The Sr/ Sr ratios of nor EuCN anomaly and time (Fig. 8), suggesting that crustal the continental crust during this time depends upon the thickness did not vary significantly over the whole of the average 87Rb/86Sr of the crust and average crustal age, Archean. Based on the stability fields of plagioclase and but assuming crustal ages between 100 and 300 my during garnet, we suggest an average continental crustal thickness this time, the Sr isotope compositions of seawater could between 35 and 45 km throughout the Archean as the best have ranged from 25% to 75% of crustal compositions rel- explanation for the REE contents, in line with previous  ative to the mantle (Table 1). These observations require a proposals noted above. Such a conclusion stands in con- re-assessment of an assumed mantle dominance on seawa- trast, however, to recent proposals (Dhuime et al., 2015) ter chemistry in the Archean (e.g., Shields and Veizer, that Archean continental crust was thin, between 15 and 2002; Kamber, 2010; Flament et al., 2013), which has found 25 km, and that temporal changes occurred in crustal thick- support from arguments that increased ridge formation and ness (Fig. 8). The proposal of Dhuime et al. (2015) is based hydrothermal circulation was high in the Archean due to an on scaling changes in Rb/Sr ratios in igneous rocks to an elevated mantle temperature (Kamber, 2015), features that assumed Rb/Sr-thickness relation in the modern crust. may provide efficient modes of heat loss on a silicate plan- Although the increases in Rb/Sr ratios with decreasing etary body (Stern, 2008). What is less understood is the pos- age in the model of Dhuime et al. (2015) is attractive in sibility that continental weathering had a significant impact the sense of providing evidence for expected increasing 87 86 on ocean chemistry during the Archean. The question we Sr/ Sr ratios with decreasing age, we feel that the scaling address in this section is that, given our current understand- used is incorrect and inconsistent with the evidence for ing of Archean crustal dynamics, how reasonable is it to thick Archean continental crust noted above. expect that Archean continental crust contributed substan- Thick Archean continental crust, by itself, does not tially to Archean seawater signals, and if this is reasonable, require high-elevation, emergent crust. If, for example, how does this conclusion require re-consideration of the Archean continental crust was relatively weak, high topog- of the Archean crust in terms of its thickness, rigid- raphy might not be supported (Rey and Coltice, 2008). ity, and composition. Weaker crust would tend to reduce continental hypsometry, A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 227

tion (e.g., Sclater et al., 1980), and higher rates of heat loss through Archean ocean environments is expected (Burke and Kidd, 1978). In fact, it has been estimated that Archean continental geothermal gradients were no more than 23 ° 1  1 C kmÀ , compared to a modern gradient of 17 °C kmÀ  (Burke and Kidd, 1978; Kramers et al., 2001, 2014), even when considering higher concentrations of radionuclides in the crust. Geologic evidence in support of a non- extreme continental geothermal gradient are (1) P-T condi- tions in metamorphic rocks (3.7 Ga-modern; Brown, 2007, 2008) show relatively constant continental geotherms over time, and (2) the preservation of coherent TTG crustal sec- tions that suggest the continental lithosphere was subject to less extreme geothermal gradients than where, for example, was generated, which likely only reflects localized high heat flow (Ernst, 2007).

6.2. Crust generation, weathering, and the start of modern- style plate tectonics

It has been suggested that Archean TTG can form in a subduction zone setting (e.g., Foley et al., 2002), or through melting of hydrous basaltic rocks deep in thickened por- tions of the crust (Smithies, 2000), possibly in oceanic pla- teaus (Zegers and van Keken, 2001). Of note, felsic continental crust can be generated from either of these models, suggesting that felsic continental crust could have existed in the . The existence of continental-like crust early in Earth history is supported by isotopic data from the Hadean Jack Hills zircon (e.g., Ushikubo et al., 2008). Most authors would agree that Archean tectonics Fig. 8. (A) The thickness of new crust versus age from Dhuime et al. (2015). The curve in (A) is calculated based on observation did not operate in exactly the same manner as modern tec- that the Rb/Sr ratio of igneous rocks in modern-day Central and tonics (see Hawkesworth et al., 2016 for a recent review). South America increases with crustal thickness. (B) Age versus Eu/ Questions have persisted on when the transition to a more * EuCN. On average, Archean TTG does not have a Eu anomaly modern-style of plate tectonics began to operate on Earth compared to post-Archean granitoids, which has been linked to (e.g., Condie and Kroner, 2008), but new research (utiliz- depth of melting (30–40 km). (C) Age versus La/Yb ratio. Many ing, in part, Hf isotopes in zircon) has suggested that Archean TTGs have low Yb contents and high La/Yb ratios, which modern-style plate tectonics began to operate by 3.2 Ga are attributed to melting deep in the crust (30–40 km) in equilib- (Næraa et al., 2012; Satkoski et al., 2013; van rium with garnet. Samples with a green square inside the white Kranendonk and Kirkland, 2016) and definitely by 3.0 Ga squares represent those TTG samples with normalized Yb values (Shirey and Richardson, 2011; Dhuime et al., 2012). The below 5. The consistently low Yb and high La/Yb ratios of TTGs suggests the presence of thick continental lithosphere throughout start of modern-style plate tectonics would have many the Archean. The blue line represents the highest La/Yb ratio implications for Earth, including the beginning of large- measured from all post-Archean granitoids. TTG data shown here scale orogenic activity (Shirey and Richardson, 2011; van are compiled from Moyen (2011). (For interpretation of the Kranendonk and Kirkland, 2016). This would have been references to color in this figure legend, the reader is referred to the a time of relatively hot mantle and punctuated continental web version of this article.) growth (van Kranendonk and Kirkland, 2016), both of which would increase the proportion of emergent continen- which is also dependent on the nature of the lithospheric tal crust in response to higher mantle potential temperature mantle (Griffin et al., 1998). A higher continental geotherm (Flament et al., 2008). has been invoked as a mechanism to reduce continental Given the evidence discussed above that Archean (4.0– topography, and some studies have argued that early to 2.5 Ga) continental crust was likely thick, and that conti- middle Archean continental crust was largely submerged nental crustal heat flow was moderate, which allows for a (Flament et al., 2013). There are, however, important argu- strong crust, and the fact that felsic continental crust could ments against higher continental geothermal gradients in exist throughout the entire Archean, there must have been a the Archean. Although it is true that overall heat flow from global-scale change in crustal dynamics that could cause a the Earth is expected to have been higher based on greater significant amount of continental crust to become emergent radiogenic heat production and a larger heat flux from the by 3.2 Ga as suggested by Satkoski et al. (2016), and con- heat of formation of the Earth, two-thirds of modern heat tinue to be emergent throughout the Archean (this study; flow is lost through oceanic crust and oceanic crust genera- Kamber and Webb, 2001). We suggest this change in crus- 228 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 tal dynamics was caused by initiation of modern-style plate seawater, which bear on continental weathering and its tectonics (Shirey and Richardson, 2011; Næraa et al., 2012; input to the oceans. Extremely high atmospheric CO2 con- Satkoski et al., 2013; van Kranendonk and Kirkland, 2016). centrations have been estimated for the early Archean to According to Dhuime et al. (2015) the beginning of large- compensate for lower solar luminosity (Kasting, 2010, scale orogenic activity due to the initiation of modern- 2014), which in turn would enhance high weathering inten- style plate tectonics would be marked by an increase of con- sities, as shown in high chemical alteration indices for early tinental detritus into the oceans. Archean clastic rocks (Hessler and Lowe, 2006 and refer- In Fig 9 we integrate temporal variations of crustal vol- ences therein), the presence of extensive first-cycle quart- ume, atmospheric CO2, crustal recycling, and Sr isotopes in zites (e.g., 3.2 Ga Moodies Group, Simpson et al., 2012),

Fig. 9. (A) The red curve plots age versus volumes of continental crust as determined by Dhuime et al. (2012). The black curves show the partial pressure of atmospheric CO2 that would be required to maintain a surface temperature of 15 and 0 °C given the differences in solar luminosity in the Archean (Kasting, 2010). The increased atmospheric CO2 relative to today would likely create a very aggressive global 87 86 weathering regime (Hessler and Lowe, 2006). (B) Age versus eHf (T) from zircon (Belousova et al., 2010). Age versus Sr/ Sr. The green (Red Lake) and yellow (Steep Rock) colored data are the best estimates of the seawater 87Sr/86Sr values at 2.94 and 2.80 Ga. The teal colored datum is from the Campbellrand carbonate (Kamber and Webb, 2001). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.) A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 229 and compositions of chemical sediments (BIF) that suggest photosynthesizers evolved early in Earth’s history, with evi- deep weathering of the Slave prior to 2.85 Ga dence for the oldest know marine oxygen gradient (and (Haugaard et al., 2016). Crustal growth curves suggest an hence ) at 3.2 Ga (Satkoski et al., 2015), as increase in continental volume from 40–60% relative to well as evidence for transient oxygen increases at 3.0 Ga  today during the (Fig. 9), and this is accom- (Planavsky et al., 2014) and 2.8 Ga (Riding et al., 2014). panied by an increase in crustal recycling, as recorded in Because the evolved continental crust is high in P, enhanced the isotopic compositions of detrital , including weathering of continental crust would produce a high flux marked decreases in eHf values (Fig. 9), as well as increases of P to the oceans (e.g., Ozaki and Tajika, 2013), and one in d18O values (Dhuime et al., 2012). We suggest that the can envision this as a trigger for bio-diversification by first sign of increased continental weathering fluxes to the 3.2 Ga based on the discussion above. While it may be  oceans coincides with this combination of high CO2- speculative that oxygenic photosynthesizers evolved in induced weathering, crustal volume, and uplifted continen- response to the initiation of modern-style plate tectonics tal crust due to the initiation of modern style plate tectonics and deep weathering of evolved, high-relief continental (Fig. 9). Although initiation of modern-style plate tectonics crust, it seems likely that large environmental changes does not necessarily require all of these components, our would have occurred in response to the large changes in interpretation is that the orogenic-related uplifts of crust crustal dynamics that are indicated by the revised Sr isotope that would occur when plate tectonics began is an attractive seawater curve for the Archean. mechanism for producing high weathering fluxes to the oceans, in addition to climate. Perhaps not coincidentally, ACKNOWLEDGMENTS this time period corresponds to the age of the first known redox gradient in the shallow marine environments, which We thank Bruno Dhuime, Kent Condie, and an anonymous has been interpreted to mark the emergence of oxygenic reviewer for constructive comments on the manuscript and Mark (Satkoski et al., 2015), and the role of Rehka¨mper for editorial handling and comments. This study was enhanced continental weathering during this time may have funded by the NASA Astrobiology Institute and NSF grant been in the delivery of nutrients to seawater, such as 1523697. PWF is supported by the Natural Sciences and Engineer- phosphorus. ing Research Council of Canada.

7. SUMMARY AND CONCLUSIONS APPENDIX A. SUPPLEMENTARY MATERIAL A re-evaluation of the Sr isotope curve for Archean sea- water places important constraints on crustal evolution, Supplementary data associated with this article can be changing the commonly held views on the Archean in terms found, in the online version, at http://dx.doi.org/10.1016/ of crustal compositions, thickness, and strength. Although j.gca.2017.04.024. some thermal models have characterized Archean crust as thin and weak, other geochemical data suggest that REFERENCES Mesoarchean continental thickness and strength was more similar to that of modern crust. This would indicate that Alibo D. S. and Nozaki Y. (1999) Rare earth elements in seawater: the continental lithosphere could support significant topog- particle association, shale-normalization, and Ce oxidation. raphy; topographic expression, therefore, would be con- Geochim. Cosmochim. Acta 63, 363–372. trolled by the amount of crust present and the global Banner J. L., Hanson G. N. and Meyers W. J. (1988) Rare earth uplift rate. During crustal convergence and orogenesis at element and Nd isotopic variations in regionally extensive 3.2–3.0 Ga, crust production and uplift to subaerial levels dolomites from the Burlington-Keokuk Formation (Mississip- would be expected to be accelerated. To keep surface tem- pian): implications for REE mobility during carbonate diage- nesis. J. . Petrol. 58, 415–432. perature above freezing during this time, atmospheric CO 2 Bau M. and Dulski P. (1996) Distribution of yttrium and rare-earth partial pressures needed to be between 0.01 and 0.1 bars elements in the Penge and Kuruman iron-formations, Trans- (Fig. 9), significantly higher than modern (0.0005 bar vaal Supergroup, South Africa. Precambr. Res. 79, 37–55. today; Kasting, 2010), which would create conditions for Beard B. L., Ludois J. M., Lapen T. J. and Johnson C. M. (2013) deep continental weathering, provided the continental crust Pre-4.0 billion weathering on Mars constrained by Rb-Sr was emergent. Considering only Archean continental on meteorite ALH84001. Earth Planet. Sci. dynamics and atmospheric composition it may be expected Lett. 361, 173–182. that a significant change in seawater chemistry should be Belousova E. A., Kostitsyn Y. A., Griffin W. L., Begg G. C., observed due to large amounts of continental runoffduring O’Reily S. Y. and Peason N. J. (2010) The growth of crustal uplift. continental crust: constraints from zircon Hf-isotope data. Lithos 119, 457–466. The new Sr isotope data presented here, combined with Broecker W. S. and Peng T. H. (1982) Tracers in the Sea. Eldigio data from Satkoski et al. (2016) and Kamber and Webb Press, Palisades, NY. (2001), suggest that seawater chemistry was affected by con- Brown M. (2007) Metamorphic conditions in orogenic belts: a tinental erosion starting at 3.2 Ga, which we interpret to record of secular change. Int. Geol. Rev. 49, 193–234. reflect the initiation of modern-style plate tectonics, which Brown M. (2008) Characteristic thermal regimes of plate tectonics continued for the rest of the Archean (Fig 9). Perhaps not and their metamorphic imprint throughout Earth history: when coincidentally, evidence continues to emerge that oxygenic did Earth first adopt a plate tectonics mode of behavior? 230 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Condie, K C. and Pease, V. (eds.). When Plate Tectonics Begin Galer S. J. G. and Mezger K. (1998) Metamorphism, denudation on Planet Earth? Geol. Soc. Am. Spec. Pap. 440, 97–128. and sea level in the Archean and cooling of the Earth. Burke K. and Kidd W. S. F. (1978) Were Archean continental Precambr. Res. 92, 389–412. geothermal gradients much steeper than those of today? Nature Griffin W. L., O’Reilly S. Y., Ryan C. G., Gaul O. and Ionov D. A. 272, 240–241. (1998) Secular variation in the composition of subcontinental Charlier B. L. A., Ginibre C., Morgan D., Nowell G. M., Pearson lithospheric mantle: geophysical and geodynamic implica- D. G., Davidson J. P. and Ottley C. J. (2006) Methods for the tionsBraun, J. et al. (ed.). Structure and Evolution of the microsampling and high-precision analysis of strontium and Australian . American Geophysical Union Geodynamics rubidium isotopes at single crystal scale for petrological and Series 26, 1–25. geochronological applications. Chem. Geol. 232, 114–133. Haugaard R., Ootes L., Creaser R. A. and Konhauser K. O. (2016) Condie K. C. (1993) Chemical composition and evolution of the The nature of Mesoarchaean seawater and continental weath- upper continental crust: contrasting results from surface sam- ering in 2.85 Ga , , NW ples and shales. Chem. Geol. 104, 1–37. Canada. Geochim. Cosmochim. Acta 194, 34–56. Condie C. C. and Kroner A. (2008) When did plate tectonics begin? Hawkesworth C. J., Cawood P. A. and Dhuime B. (2016) Tectonics Evidence from the geologic record. Geol. Soc. Am. Spec. Pap. and crustal evolution. GSA Today 26, 4–11. 440, 281–294. Hessler A. M. and Lowe D. R. (2006) Weathering and sediment Condie K. C. and Pease V. (2008) When did plate tectonics begin generation in the Archean: an integrated study of the evolution on planet Earth? Geol. Soc. Am. Spec. Pap. 440, 294. of siliciclastic sedimentary rocks of the 3.2 Ga Moodies Group, Corfu F. and Wallace H. (1986) U-Pb zircon ages from magmatism Barberton Greenstone Belt, South Africa. Precambr. Res. 151, in the Red Lake greenstone belt, northwestern Ontario. Can. J. 185–210. Earth Sci. 23, 967–977. Hodell D. A., Mead G. A. and Mueller P. A. (1990) Variation in Danielson A., Muller P. and Dulski P. (1992) The europium the strontium isotopic composition of seawater (8 Ma to anomalies in banded iron formations and the thermal history of present): implications for chemical weathering rates and the ocean crust. Chem. Geol. 97, 89–100. dissolved fluxes to the oceans. Chem. Geol. 80, 291–307. Davies G. F. (1992) On the emergence of plate tectonics. Geology Hofmann H. J., Thurston P. C. and Wallace H. (1985) Archean 20, 963–966. Stromatolites from the Uchi greenstone belt, northwestern Dhuime B., Hawkesworth C. J., Cawood P. A. and Storey C. D. Ontario. In Evolution of Archean Supracrustal Sequences, GAC (2012) A change in geodynamics of continental growth 3 billion Special Paper No. 28. GAC, St. Johns, Newfoundland, pp. 125– years ago. Science 335, 1334–1336. 132. Dhuime B., Wuestefeld A. and Hawkesworth C. J. (2015) Johnson C. M. and Thompson R. A. (1991) Isotopic composition Emergence of modern continental crust about 3 billion years of mafic volcanic rocks in the northern Rio Grande ago. Nat. Geosci. 8, 552–555. : evidence for contributions of ancient intraplate and Edgar K. M., Anognostou E., Pearson P. N. and Foster G. L. subduction magmatism to evolution of the lithosphere. J. (2015) Assessing the impact of diagenesis on d11B, d13C, d18O, Geophys. Res. 96, 13593–13608. Sr/Ca and B/Ca values in planktic foraminiferal calcite. Kamber B. S. (2010) Archean mafic-ultramafic volcanic landmasses Geochim. Cosmochim. Acta 166, 189–209. and their effect on ocean-atmosphere chemistry. Chem. Geol. Edwards C. T., Saltzman M. R., Leslie S. A., Bergstrom S. M., 274, 19–28. Sedlacek A. R. C., Howard A., Bauer J. A., Sweet W. C. and Kamber B. S. (2015) The evolving nature of terrestrial crust from Young S. A. (2015) Strontium isotope (87Sr/86Sr) stratigraphy the Hadean, through the Archaean, into the Proterozoic. of bulk carbonate: implications for preservation of Precambr. Res. 258, 48–82. primary seawater values. GSA Bull. 127, 1275–1289. Kamber B. S. and Webb G. E. (2001) The geochemistry of late Ernst W. G. (2007) Speculations on evolution of the terrestrial Archean microbial carbonate: implications for ocean chemistry lithosphere–asthenosphere system—plumes and plates. Gond- and continental erosion history. Geochim. Cosmochim. Acta 65, wana Res. 11, 38–49. 2509–2525. Flament N., Coltice N. and Rey P. F. (2008) A case for late- Kasting J. F. (2010) Faint young redux. Nature 464, 687–689. Archaean continental emergence from thermal evolution mod- Kasting J. F. (2014) Atmospheric composition of Hadean-early els and hypsometry. Earth Planet. Sci. Lett. 275, 326–336. Archean Earth: the importance of CO. GSA Spec. Pap. 504, 19– Flament N., Cotice N. and Rey P. F. (2013) The evolution of the 28. 87Sr/86Sr of marine carbonates does not constrain continental Kemp A. I. S. and Hawkesworth C. J. (2003) Granitic perspectives growth. Precambr. Res. 229, 177–188. on the generation and secular evolution of the continental crust. Fralick P., Hollings P. and King D. (2008) Stratigraphy, geochem- In Treatise on Geochemistry: The Crust (ed. R. L. Rudnick). istry, and depositional environments of Mesoarchean sedimen- Elsevier, Oxford, pp. 349–410. tary units in western Superior Province. Implications for Kleinhanns I. C., Kramers J. D. and Kamber B. S. (2003) generation of early crustCondie, K. C. and Pease, V. (eds.). Importance of water from Archaean granitoid petrology: a When Plate Tectonics Begin on Planet Earth? Geol. Soc. Am. comparative study of TTG and potassic granitoids from Spec. Pap. 440, 77–96. Barberton Mountain Land, South Africa. Contrib. Miner. Fralick P. and Riding R. (2015) Steep Rock Lake: sedimentology Petrol. 145, 377–389. and geochemistry of an Archean carbonate platform. Earth Sci. Kramers J. D., Hezen M. and Steidle L. (2014) Greenstone belts at Rev. 151, 132–175. the northernmost edge of the : timing of Foley S. (2008) A trace element perspective on Archean crust tectonic events and a possible crustal fluid source. Precambr. formation and on the presence or absence of Archean subduc- Res. 253, 96–113. tionCondie, K. C. and Pease, V. (eds.). When did Plate Tectonics Kramers J. D., Kreissig K. and Jones M. Q. W. (2001) Crustal heat Begin on Planet Earth. Geol. Soc. Am. Spec. Pap. 440, 31–50. production and style of metamorphism: a comparison between Foley S., Tiepolo M. and Vannucci R. (2002) Growth of early two Archean high grade provinces in the Limpopo Belt, continental crust controlled by melting of amphibolite in southern Africa. Precambr. Res. 112, 149–163. subduction zones. Nature 417, 837–840. A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232 231

McCulloch M. T. (1994) Primitive 87Sr/86Sr from an Archean Sclater J. G., Jaupart C. and Galson D. (1980) The heat flow barite and conjecture on the age and origin. Earth through the oceanic and continental crust and the heat loss of Planet. Sci. Lett. 126, 1–13. the earth. Rev. Geophys. 18, 269–312. McCulloch M. T., Gregory R. T., Wasserburg G. J. and Taylor H. Sen C. and Dunn T. (1994) Dehydration melting of a basaltic P. (1981) Sm-Nd, Rb-Sr, and 18O/16O isotopic systematics in an composition amphibolite at 1.5 and 2.0 GPa: implications for oceanic crustal section: evidence from the Samail . J. the origin of adakites. Contrib. Miner. Petrol. 117, 345–364. Geophys. Res. 86, 2721–2735. Sena C. M., John C. M., Jourdan A., Vandeginste V. and Manning McIntyre Ti (2014) Sedimentology and Geochemistry of Mesoarch- C. (2014) Dolomitization of lower peritidal carbon- ean Chemical Sediments of the Red Lake and Wallace Lake ates by modified seawater: constraints from clumped isotopic Greenstone Belts Honours Thesis. Lakehead University, Thun- paleothermometry, elemental chemistry, and strontium iso- der Bay, Ontario, Canada, 173 pp. topes. J. Sediment. Res. 84, 552–566. Moyen J. (2011) The composite Archaean grey : petrolog- Shields G. A. (2007) A normalized seawater strontium isotope ical significance, and evidence for a non-unique tectonic setting curve: possible implications for - for Archaean crustal growth. Lithos 123, 21–36. weathering rates and the further oxygenation of the Earth. Næraa T., Scherste´n Rosing M. T., Kemp A. I. S., Hoffmann J. E., eEarth 2, 35–42. Kokfelt T. F. and White- house M. J. (2012) Hafnium isotope Shields G. A. and Veizer J. (2002) marine carbonate evidence for a transition in the dynamics of continental growth isotope database: version 1.1. Geochem. Geophys. Geosyst. 3, 1– 3.2 Ga ago. Nature 485, 627–631. 12. Nozaki Y., Zhang J. and Amakawa H. (1997) The fractionation Shirey S. B. and Richardson S. H. (2011) Start of the between T and Ho in the marine environment. Earth Planet. at 3 Ga shown by diamonds from subcontinental mantle. Lett. 197, 329–340. Science 333, 434–436. Ozaki K. and Tajika E. (2013) Biogeochemical effects of atmo- Simpson E. L., Eriksson K. A. and Mueller W. U. (2012) 3.2 Ga spheric oxygen concentration, phosphorus weathering, and sea- eolian deposits from the Moodies Group, Barberton Green- level stand on oceanic redox chemistry: implications for stone Belt, South Africa: implications for the origin of first- greenhouse climates. Earth Planet. Sci. Lett. 373, 129–139. cycle quartz . Precambr. Res. 214–215, 185–191. Pirie J. (1981) Regional setting of deposits in the Red Lake Smithies R. H. (2000) The Archean tonalite-trondhjemite-granodi- area, northwestern Ontario. In Genesis of Archean Volcanic- orite (TTG) series is not an analogue of adakite. hosted Gold Deposits. Ontario Geological Survey, Miscellaneous Earth Planet. Sci. Lett. 182, 115–125. Paper No. 97. Ontario Geological Survey, Ontario, Canada, pp. Stern R. J. (2008) Modern-style plate tectonics began in Neopro- 71–93. terozoic time: an alternative interpretation of Earth’s tectonic Planavsky N. J., Asael D., Hofmann A., Reinhard C. T., Lalonde historyCondie, K. C. and Pease, V. (eds.). When Plate Tectonics S. V., Knudsen A., Wang X., Ossa Ossa F., Recoits E., Smith Begin on Planet Earth? Geol. Soc. Am. Spec. Pap. 440, 265–280. A. J. B., Beukes N. J., Bekker A., Johnson T. M., Konhauser Sumner D. Y. and Grotzinger J. P. (2000) Late Archean aragonite K. O., Lyons T. W. and Rouxel O. J. (2014) Evidence for precipitation; petrography, facies associations and environmen- oxygenic photosynthesis half a billion years before the Great tal significance. In Carbonate Sedimentation and Diagenesis in Oxidation Event. Nat. Geosci. 7, 283–286. the Evolving Precambrian World, vol. 33 (eds. J. P. Grotzinger Pons M. L., Fujii T., Rosing M., Quitte G., Telouk P. and and N. P. James). Society of Economic Paleontologists and Albarede F. (2013) A Zn isotope perspective on the rise of Mineralogists Special Publication, pp. 107–120. continents. 11, 201–214. Tachikawa K., Athias V. and Jeandel C. (2003) Neodymium Prokoph A., Shields G. A. and Veizer J. (2008) Compilation of budget in the modern ocean and paleo-oceanographic implica- time-series analysis of marine carbonate d18O, d13C, 87Sr/86Sr tions. J. Geophys. Res. 108, 3254. and d34S database through Earth history. Earth Sci. Rev. 87, Tanaka T., Togashi S., Kamioka H., Amakawa H., Kagami H., 113–133. Hamamoto T., Yuhara T., Orihashi Y., Yoneda S., Shimizu H., Rapp R. P. and Watson E. B. (1995) Dehydration melting of Kunimaru T., Takahashi K., Yanagi T., Nakano T., Fujimaki metabasalt at 8–32 kbar – implications for continental growth H., Shinjo R., Asahara Y., Tanimizu M. and Dragusanu C. and crust–mantle recycling. J. Petrol. 36, 891–931. (2000) Jndi-1: a neodymium isotopic reference in consistency Rapp R. P., Shimizu N. and Norman M. D. (2003) Growth of early with LaJolla neodymium. Chem. Geol. 168, 279–281. continental crust by partial melting of eclogite. Nature 425, Taylor S. R. and McLennan S. M. (1985) The Continental Crust: 605–609. Its Composition and Evolution. Blackwell, Cambridge, Mass. Rey P. F. and Coltice N. (2008) Neoarchean strengthening of the Tomlinson K. Y., Davis D. W., Stone D. and Hart T. R. (2003) U- lithosphere and the coupling of the Earth’s geochemical Pb age and Nd isotopic evidence for terrane development and reservoirs. Geology 36, 635–638. crustal recycling in the south-central Wabigoon Subprovince, Riding R., Fralick P. and Liang L. (2014) Identification of an Canada. Contrib. Miner. Petrol. 144, 684–702. Archean marine oxygen oasis. Precambr. Res. 251, 232–237. Tomlinson K. Y., Stevenson R. K., Hughes D. J., Hall R. P., Satkoski A. M., Beukes N. J., Li W., Beard B. L. and Johnson C. Thurston P. C. and Henry P. (1998) The Red Lake greenstone M. (2015) A redox-stratified ocean 3.2 billion years ago. Earth belt, Superior Province: evidence of plume-related magmatism Planet. Sci. Lett. 430, 43–53. at 3 Ga and evidence of an older enriched source. Precambr. Satkoski A. M., Bickford M. E., Samson S. D., Bauer R. L., Res. 89, 59–76. Mueller P. A. and Kamenov G. D. (2013) Geochemical and Hf- Ushikubo T., Noriko K. T., Cavosie A. J., Wilde S. A., Rudnick R. Nd isotopic constraints on the crustal evolution of Archean L. and Valley J. W. (2008) Lithium in Jack Hills zircons: rocks from the Minnesota River Valley, USA. Precambr. Res. evidence for extensive weathering of Earth’s earliest crust. Earth 224, 36–50. Planet. Sci. Lett. 272, 666–676. Satkoski A. M., Lowe D. R., Beard B. L., Coleman M. L. and Van Kranendonk M. J. and Kirkland C. L. (2016) Conditioned Johnson C. M. (2016) A high continental weathering flux into duality of the Earth system: geochemical tracing of the Paleoarchean seawater revealed by strontium isotope analysis cycle through Earth history. Earth Sci. Rev. of 3.26 Ga barite. Earth Planet. Sci. Lett. 454, 28–35. 160, 171–187. 232 A.M. Satkoski et al. / Geochimica et Cosmochimica Acta 209 (2017) 216–232

Van Kranendonk M. J., Webb G. E. and Kamber B. S. (2003) Winter B. L., Johnson C. M. and Clark D. L. (1997) Strontium, Geological and trace element evidence for a marine sedimentary neodymium and lead isotope variations of authigenic and environment of deposition and biogenicity of 3.45 Ga stroma- silicate sediment components from the Late Ceno- zoic Arctic tolitic carbonates in the , and support for a Ocean: implications for and the source of trace reducing Archaean ocean. Geobiology 1, 91–108. metals in seawater. Geochim. Cosmochim. Acta 19, 4181–4200. Veizer J., Compston W., Hoefs R. W. and Nielson H. (1982) Wolf M. B. and Wyllie P. J. (1993) Garnet growth during Mantle buffering of the early oceans. Naturwissenschaften 69, amphibolite anatexis – implications of a garnetiferous restite. 173–180. J. Geol. 101, 357–373. Veizer J., Hoefs J., Lowe D. R. and Thurston P. C. (1989) Zegers T. E. and van Keken P. E. (2001) Middle Archean continent Geochemistry of Precambrian carbonates: II. Archean green- formation by crustal delamination. Geology 29, 1083–1086. stone belts and Archean seawater. Geochim. Cosmochim. Acta Zhong S. and Mucci A. (1995) Partitioning or rare earth elements 53, 859–871. (REEs) between calcite and seawater solutions at 25 °C and Webb G. E. and Kamber B. S. (2000) Rare earth elements in 1 atm, and high dissolved REE concentrations. Geochim. Holocene reefal microbialites: a new shallow seawater proxy. Cosmochim. Acta 59, 443–453. Geochim. Cosmochim. Acta 64, 1557–1565. Wilks M. E. and Misbet E. G. (1985) Archean stromatolites from Associate editor: Mark Rehkamper the Steep Rock Group, northwest Ontario, Canada. Can. J. Earth Sci. 25, 370–391.