Rise of major subaerial landmasses about 3.0 to 2.7 billion years ago
Affiliations | Corresponding Author | Cite as | Funding information- Share this article
-
Article views:626Cumulative count of HTML views and PDF downloads.
- Download Citation
- Rights & Permissions
top
Abstract
Figures
Figure 1 K/La of fresh and hydrothermally altered basalts, exemplified by DSDP and ODP sites 417A, 417D, 418A, 801C, 843B, 1140D and 1149B (see Supplementary Information). Mean normal mid-ocean ridge basalts (N-MORB) and enriched mid-ocean ridge basalts (E-MORB) are from Sun and McDonough (1989). | Figure 2 (a) Secular mean K/La evolution of global continental mafic volcanic and plutonic rocks. Error bars show 1 s.e.m. uncertainties. Modern geological reservoirs are shown as grey arrows with numbers: (1) upper continental crust, (2) continental arc basalts, (3) lower continental crust, (4) primitive mantle, (5) E-MORB, (6) ocean island basalts (OIB), (7) N-MORB (Sun and McDonough, 1989; Kelemen et al., 2014; Rudnick and Gao, 2014). Secular changes in (b) normalised seawater 87Sr/86Sr curve (Shields et al., 2007), (c) the juvenile continental crustal thickness (Dhuime et al., 2015), (d) the upper continental crustal mass (Tang et al., 2016), and (e) atmospheric O2 level (Lyons et al., 2014) are also shown for comparison. Global detrital zircon U-Pb age distribution (Voice et al., 2011) is plotted at the top of the figure, with the periods of supercontinents assembly (Hawkesworth et al., 2017). | Figure 3 Schematic diagrams illustrating submarine versus subaerial states of landmasses in the Archean era. (a) Largely submarine state of proto-landmasses at pre-3.0 Ga. (b) Rise of subaerial landmasses during ∼3.0–2.7 Ga, likely driven by initiation of plate tectonics and assembly of the first supercontinent. These tectono-magmatic events can create high elevations as happens at present day (c) Tibet and (d) Andes plateaus (see Supplementary Information). |
Figure 1 | Figure 2 | Figure 3 |
top
Introduction
Determining when major subaerial landmasses first emerged and how they rose is the key to understand secular evolution of Earth’s ecosystem, but remains challenging. Static models considering multiple parameters (e.g., crustal thickness, density and mantle potential temperature) suggest the rise of major subaerial landmasses from the late Archean to the Neoproterozoic (Flament et al., 2008
Flament, N., Coltice, N., Rey, P.F. (2008) A case for late-Archaean continental emergence from thermal evolution models and hypsometry. Earth and Planetary Science Letters 275, 326–336.
; Lee et al., 2017Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
). Exposed landmasses have been traced back to as early as ∼3.5 Ga (e.g., detrital sediments in the East Pilbara Craton; Buick et al., 1995Buick, R., Thornett, J.R., McNaughton, N.J., Smith, J.B., Barley, M.E., Savage, M. (1995) Record of emergent continental crust ∼3.5 billion years ago in the Pilbara craton of Australia. Nature 375, 574–577.
; Campbell and Davies, 2017Campbell, I.H., Davies, D.R. (2017) Raising the continental crust. Earth and Planetary Science Letters 460, 112–122.
), as indicated by geological records of uplift-subsidence cycles (e.g., in the East Pilbara Craton during ∼3.5–3.2 Ga; Campbell and Davies, 2017Campbell, I.H., Davies, D.R. (2017) Raising the continental crust. Earth and Planetary Science Letters 460, 112–122.
). Early subaerial landmasses may have experienced subsidence and eventually submerged below sea level again, due to crustal erosion and thermal relaxation of continental lithospheric mantle root (Lee et al., 2017Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
). Accordingly, it remains unclear to what extent sporadic geological records can reflect the relative proportion of subaerial landmasses with time.An increase in subaerial large igneous provinces (LIPs) and a stepwise change of trip-oxygen-isotope (Δ17O) composition of shales indicate major landmasses could have rapidly emerged above sea level at the Archean-Proterozoic boundary (Kump and Barley, 2007
Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
; Bindeman et al., 2018Bindeman, I.N., Zakharov, D.O., Palandri, J., Greber, N.D., Dauphas, N., Retallack, G.J., Hofmann, A., Lackey, J.S., Bekker, A. (2018) Rapid emergence of subaerial landmasses and onset of a modern hydrologic cycle 2.5 billion years ago. Nature 557, 545–548.
). Most intriguingly, it was coincident with Earth’s first blast of molecular oxygen (i.e. the Great Oxidation Event (GOE)) and has been proposed as a key driving force for the GOE, probably by increasing oxidised subaerial volcanic gases (Kump and Barley, 2007Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
; Gaillard et al., 2011Gaillard, F., Scaillet, B., Arndt, N.T. (2011) Atmospheric oxygenation caused by a change in volcanic degassing pressure. Nature 478, 229–232.
) and decreasing O2 sinks (e.g., H2, CH4) of weathering and alteration (Smit and Mezger, 2017Smit, M.A., Mezger, K. (2017) Earth’s early O2 cycle suppressed by primitive continents. Nature Geoscience 10, 788–792.
). Both the number of subaerial LIPs and the Δ17O values of shales, however, appear too muted to reveal the emergence of subaerial landmasses. Here we show the mean K/La of global continental mafic volcanic rocks as a novel proxy to unravel the emergence history of subaerial landmasses.top
Rationale & Results
As highly incompatible elements, K and La do not fractionate from each other significantly during mafic magma processes. In contrast, K and La behaviour are very different during weathering and alteration at the surface. K is water-soluble and can thereby be easily removed during chemical weathering, while La is largely insoluble and thus immobile in this process. Weathered subaerial mafic rocks thus on average have depleted K/La (Ma et al., 2007
Ma, J.-L., Wei, G.-J., Xu, Y.-G., Long, W.-G., Sun, W.-D. (2007) Mobilization and re-distribution of major and trace elements during extreme weathering of basalt in Hainan Island, South China. Geochimica et Cosmochimica Acta 71, 3223–3237.
). With input from riverine runoff and high temperature (>100 or 150 °C) submarine hydrothermal processes (Hofmann and Harris, 2008Hofmann, A., Harris, C. (2008) Silica alteration zones in the Barberton greenstone belt: A window into subseafloor processes 3.5–3.3 Ga ago. Chemical Geology 257, 221–239.
; Staudigel, 2014Staudigel, H. (2014) 4.16 – Chemical Fluxes from Hydrothermal Alteration of the Oceanic Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 583–606.
), K is highly enriched relative to La in seawater. For submarine magma eruptions, the K/La of altered mafic volcanic rocks can be strongly elevated during low temperature hydrothermal alteration (<100 or 150 °C) with seawater (e.g., hydrothermally altered basalts; Fig. 1) (Hofmann and Harris, 2008Hofmann, A., Harris, C. (2008) Silica alteration zones in the Barberton greenstone belt: A window into subseafloor processes 3.5–3.3 Ga ago. Chemical Geology 257, 221–239.
; Staudigel, 2014Staudigel, H. (2014) 4.16 – Chemical Fluxes from Hydrothermal Alteration of the Oceanic Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 583–606.
; Supplementary Information). Therefore, mafic volcanic rocks in submarine eruptions are expected to possess systematically higher K/La compared to subaerial eruptions. Accordingly, the mean K/La of global continental mafic volcanic rocks could be an ideal proxy to track the relative proportion of submarine to subaerial eruptions over time.The weighted bootstrap re-sampling method with Monte Carlo simulations was applied to reveal secular changes in the K/La of less differentiated (SiO2: 43–51 wt. %) global continental mafic volcanic and plutonic rocks (Fig. S-1) (Keller and Schoene, 2012
Keller, C.B., Schoene, B. (2012) Statistical geochemistry reveals disruption in secular lithospheric evolution about 2.5 Ga ago. Nature 485, 490–493.
; Supplementary Information). Mean K/La values of continental igneous rocks were calculated in 100 Ma bins with uncertainty of 1 standard error of the mean (1 s.e.m.) (Fig. 2a, Table S-1). Mafic plutons intruded at depths and thus have been least affected by superficial seawater-rock interactions (Fig. S-2). Accordingly, these rocks can serve as a magmatic background, which yields a nearly constant mean K/La since 3.6 Ga (582 ± 92, 1 standard deviation; 1 s.d.), whereas the mean K/La of mafic volcanic rocks exhibits a substantial change through time. What is striking in mafic volcanic rocks is much higher mean K/La at pre-3.0 Ga—yielding an average value of 2125 ± 309 (1 s.d.)—that progressively decline to near the magmatic background (584 ± 83, 1 s.e.m.) at 2.7 Ga. Three transient positive excursions of K/La have been identified for post-2.7 Ga mafic volcanic rocks, i.e. at 0.4, 0.8 and 2.6 Ga, coincident with global tectono-magmatism lulls (Voice et al., 2011Voice, P.J., Kowalewski, M., Eriksson, K.A. (2011) Quantifying the timing and rate of crustal evolution: global compilation of radiometrically dated detrital zircon grains. Journal of Geology 119, 109–126.
; Hawkesworth et al., 2017Hawkesworth, C.J., Cawood, P.A., Dhuime, B., Kemp, T.I.S. (2017) Earth’s continental lithosphere through time. Annual Review of Earth and Planetary Sciences 45, 169–198.
). Without considering these excursions, the mean K/La of post-2.7 Ga mafic volcanic rocks yield an average of 536 ± 74 (1 s.d.) statistically lower than the magmatic background (t test, p < 0.01). It indicates distinguishable K depletion relative to La in post-2.7 Ga volcanic rocks relative to their intrusive equilibrants (Fig. 2a).top
Discussion
With the development of plate tectonics (Tang et al., 2016
Tang, M., Chen, K., Rudnick, R.L. (2016) Archean upper crust transition from mafic to felsic marks the onset of plate tectonics. Science 351, 372–375.
), differentiation in subduction zones and recycling of crustal materials created substantial mantle heterogeneity. On the one hand, influx of more K relative to La into the mantle wedge via aqueous fluids derived from subducted sediments and oceanic slabs (Kelemen et al., 2014Kelemen, P.B., Hanghøj, K., Greene, A.R. (2014) 4.21 – One View of the Geochemistry of Subduction-Related Magmatic Arcs, with an Emphasis on Primitive Andesite and Lower Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 749–806.
), tends to produce arc basalts with elevated K/La (e.g., average of continental arc basalts; Fig. 2a). On the other hand, recycling of dehydrated crustal materials can account for some intra-plate basalts with depleted K/La (e.g., average K/La values of Deccan Traps and Siberia LIPs are ∼201 and ∼161, respectively; Farmer, 2014Farmer, G.L. (2014) 4.3 – Continental Basaltic Rocks. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 75–110.
). However, development of mantle reservoirs with differentiated K/La (e.g., OIB, N-MORB, E-MORB; Sun and McDonough, 1989Sun, S.-S., McDonough, W.F. (1989) Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes. Geological Society, London, Special Publications 42, 313–345.
) appears not to have significantly changed the mean K/La of global mafic magmas with time, as illustrated by nearly constant mean K/La of mafic plutons since 3.6 Ga (Fig. 2a).Crustal contamination can elevate K/La of continental mafic magmas to various extents by interaction with felsic continental crust with a higher K/La (Fig. 2a). On the contrary, K2O/Na2O of felsic continents increases since ∼3.0 Ga reaching a modern value at ∼2.5 Ga (Fig. S-3). Contamination by felsic continental crust thus cannot explain the nearly contemporaneous decline in the mean K/La of mafic volcanic rocks. Further, mean K/La values of global mafic volcanic rocks at pre-2.7 Ga are too high (∼1039–2580; Table S-1) to be explained by crustal contamination, as the mean K/La of the modern upper continental crust is only ∼750 (Rudnick and Gao, 2014
Rudnick, R.L., Gao, S. (2014) 4.1 – Composition of the Continental Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 1–51.
; Fig. 2a). Archean mafic volcanic and plutonic samples could both have been subjected to metamorphism. However, the nearly invariable mean K/La of mafic plutonic rocks over time indicates that metamorphism may have not significantly changed mean K/La of global mafic rocks.The mean K/La of mafic volcanic rocks prior to ∼3.0 Ga is too high to be subjected to any magmatic and metamorphic origin. Instead, such high values, to our knowledge, can only be explained by influence of widespread low temperature seawater-rock interactions on these mafic volcanic rocks (Figs. 1, S-2). Therefore, the remarkably elevated mean K/La of pre-3.0 Ga mafic volcanic rocks reveals magma eruptions predominantly at a submarine environment (Fig. 3a), supporting the consensus that most early Archean landmasses were still below sea level (Bindeman et al., 2018
Bindeman, I.N., Zakharov, D.O., Palandri, J., Greber, N.D., Dauphas, N., Retallack, G.J., Hofmann, A., Lackey, J.S., Bekker, A. (2018) Rapid emergence of subaerial landmasses and onset of a modern hydrologic cycle 2.5 billion years ago. Nature 557, 545–548.
). The mantle-like seawater 87Sr/86Sr compositions in pre-3.0 Ga (Shields, 2007Shields, G.A. (2007) A normalised seawater strontium isotope curve: possible implications for Neoproterozoic-Cambrian weathering rates and the further oxygenation of the Earth. eEarth 2, 35–42.
; Fig. 2b) may also reflect limited erosion of proto-landmasses and a largely submarine state for early Archean landmasses.Potassium mobility could be quite different during anoxic versus oxic alteration of basalts (Jagoutz, 2012
Jagoutz, O. (2012) Were ancient granitoid compositions influenced by contemporaneous atmospheric and hydrosphere oxidation states? Terra Nova 25, 95–101.
). Saponite forms in anoxic altered basalts while celadonite forms in oxic altered basalts. The former mineral contains a much lower K2O content than the latter one. Therefore, transition from widespread anoxic to oxic alteration across the GOE would predict a systematic increase in the K/La of global mafic volcanic rocks, which is inconsistent with the reverse trend observed in Figure 2a. Instead, the progressive decrease in mean K/La of mafic volcanic rocks during ∼3.0–2.7 Ga may reflect a gradual decline in the proportion of magmas that erupted at submarine environments and experienced hydrothermal alteration, in other words, the emergence of subaerial landmasses (Fig. 3b). After this transition, the relatively steady mean K/La of mafic volcanic rocks suggests that, on the whole, the area of subaerial landmasses has not significantly changed since ∼2.7 Ga. Statistical analysis on the mean K/La of global mafic volcanic and plutonic rocks thus provides a compelling record for the rise of major subaerial landmasses during ∼3.0–2.7 Ga, likely to a present day level since ∼2.7 Ga (Figs. 2a, 3b).The gradual rise of major subaerial landmasses during ∼3.0–2.7 Ga can be supported by multiple geological observations. Emergence of large subaerial landmasses initially occurred at ∼3.0 Ga at the Kaapvaal craton (e.g., Pongola-Wiwatersrand Basin; with an area of hundreds of square kilometres of exposed surface; Kump and Barley, 2007
Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
; Korenaga et al., 2017Korenaga, J., Planavsky, N.J., Evans, D.A.D. (2017) Global water cycle and the coevolution of the Earth’s interior and surface environment. Philosophical Transactions of the Royal Society A-Mathematical Physical and Engineering Sciences 375, 20150393.
) that then extended to various cratons at ∼2.7 Ga (e.g., Kalgoorlie-Kambalda, Abitibi, Slave, Dharwar, Zimbabwe, Wyoming, Kaapvaal and Pilbara; Kump and Barley, 2007Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
; Campbell and Davies, 2017Campbell, I.H., Davies, D.R. (2017) Raising the continental crust. Earth and Planetary Science Letters 460, 112–122.
). Increasing seawater 87Sr/86Sr isotopic ratios since ∼3.0 Ga may reflect more extensive chemical weathering of exposed landmasses at that time (Fig. 2b). The earliest mature detrital zircon age distributions at ∼3.0–2.8 Ga also indicates widespread rise of subaerial continents at that time (Reimink et al., 2021Reimink, J.R., Davies, J.H.F.L., Ielpi, A. (2021) Global zircon analysis records a gradual rise of continental crust throughout the Neoarchean. Earth and Planetary Science Letters 554, 116654.
).The area of global subaerial landmasses in Earth’s history is regulated by water volume in the hydrosphere, topography differentiation driven by tectono-magmatic events (likely what happens at present day Tibet and Andes plateaus; Fig. 3c,d) as well as processes that eliminate topography difference, e.g., weathering erosion and subsidence due to thermal relaxation of continental lithospheric mantle root (Lee et al., 2017
Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
). More than 50 % water in the hydrosphere has been suggested to be returned back to the mantle through subduction in the past (Korenaga et al., 2017Korenaga, J., Planavsky, N.J., Evans, D.A.D. (2017) Global water cycle and the coevolution of the Earth’s interior and surface environment. Philosophical Transactions of the Royal Society A-Mathematical Physical and Engineering Sciences 375, 20150393.
). Here we consider that the decline in the total water volume of the hydrosphere was minor during the rise of major subaerial landmasses, as Archean subduction zones, if they existed, were hot and deep subduction of water should have been prevented by prevailing slab melting (Rapp et al., 2003Rapp, R.P., Shimizu, N., Norman, M.D. (2003) Growth of early continental crust by partial melting of eclogite. Nature 425, 605–609.
).Crustal thickness plays an important role in controlling continental elevations (Lee et al., 2017
Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
). Rapid growth of crustal thickness as well as felsic components since ∼3.0 Ga (Tang et al., 2016Tang, M., Chen, K., Rudnick, R.L. (2016) Archean upper crust transition from mafic to felsic marks the onset of plate tectonics. Science 351, 372–375.
; Dhuime et al., 2015Dhuime, B., Wuestefeld, A., Hawkesworth, C.J. (2015) Emergence of modern continental crust about 3 billion years ago. Nature Geoscience 8, 552–555.
) should have increased continental elevations. However, static models with consideration of the crustal thickness, density and mantle potential temperature predict the rise of major subaerial landmasses until the Neoproterozoic (Lee et al., 2017Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
). Instead, the rise of major subaerial landmasses during ∼3.0–2.7 Ga can be explained by initiation of plate tectonics (∼3.0 Ga; Tang et al., 2016Tang, M., Chen, K., Rudnick, R.L. (2016) Archean upper crust transition from mafic to felsic marks the onset of plate tectonics. Science 351, 372–375.
) and assembly of the Earth’s first supercontinent (Kenorland), which was marked by the first burst of felsic magmatism at ∼2.7–2.5 Ga (Voice et al., 2011Voice, P.J., Kowalewski, M., Eriksson, K.A. (2011) Quantifying the timing and rate of crustal evolution: global compilation of radiometrically dated detrital zircon grains. Journal of Geology 119, 109–126.
; Hawkesworth et al., 2017Hawkesworth, C.J., Cawood, P.A., Dhuime, B., Kemp, T.I.S. (2017) Earth’s continental lithosphere through time. Annual Review of Earth and Planetary Sciences 45, 169–198.
; Fig. 2). Continuous subduction with the episodic assembly of supercontinents could have provided persistent mountain building processes globally, which must have largely overcome weathering erosion and thermal subsidence to keep the area of global subaerial landmasses nearly constant since ∼2.7 Ga. The three K/La positive excursions of mafic volcanic rocks may have recorded several periods with a shrunken area of global subaerial landmasses, which can be explained by weaker orogeny during tectono-magmatic lulls (Fig. 2a).Submarine landmasses universal in the early Archean, indicated by the extremely high mean K/La of mafic volcanic rocks (Figs. 2a, 3a), with widespread submarine hydrothermal processes and serpentinisation of ultramafic-mafic rocks could have buffered a reduced atmosphere (Smit et al., 2017
Smit, M.A., Mezger, K. (2017) Earth’s early O2 cycle suppressed by primitive continents. Nature Geoscience 10, 788–792.
) and provided proper environments (e.g., N2, H+, temperature) for early life (Sleep et al., 2011Sleep, N.H., Bird, D.K., Pope, E.C. (2011) Serpentinite and the dawn of life. Philosophical Transactions of the Royal Society B-Biological Sciences 366, 2857–2869.
). For instance, H2-rich fluids were likely to support the prosperity of methanogens and promote complex organic compounds (biotic or prebiotic life) in pristine seas (Hofmann and Harris, 2008Hofmann, A., Harris, C. (2008) Silica alteration zones in the Barberton greenstone belt: A window into subseafloor processes 3.5–3.3 Ga ago. Chemical Geology 257, 221–239.
). Intriguingly, the emergence of subaerial landmasses has been proposed as a key driving force for the GOE (Kump and Barley, 2007Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
; Gaillard et al., 2011Gaillard, F., Scaillet, B., Arndt, N.T. (2011) Atmospheric oxygenation caused by a change in volcanic degassing pressure. Nature 478, 229–232.
). A “time lag” is identified here between the rise of major subaerial landmasses and the initiation of the GOE (Lyons et al., 2014Lyons, T.W., Reinhard, C.T., Planavsky, N.J. (2014) The rise of oxygen in Earth’s early ocean and atmosphere. Nature 506, 307–315.
; Fig. 2e). This time lag may not be explained by the delayed response of the surficial systems to evolution of the solid Earth, if a steady atmospheric oxygen level could be established at a time scale <50 Ma (Laakso and Schrag, 2014Laakso, T.A., Schrag, D.P. (2014) Regulation of atmospheric oxygen during the Proterozoic. Earth and Planetary Science Letters 388, 81–91.
). Accordingly, the emergence of subaerial landmasses during ∼3.0–2.7 Ga most likely acted as a precondition rather than a trigger for the GOE. This study suggests alternative mechanisms, e.g., the evolution of oxygen photosynthesis, sedimentary burial of organic carbon (Lyons et al., 2014Lyons, T.W., Reinhard, C.T., Planavsky, N.J. (2014) The rise of oxygen in Earth’s early ocean and atmosphere. Nature 506, 307–315.
) and/or deep oxygen cycles driven by volatiles (He et al., 2019He, Y.-S., Meng, X.-N., Ke, S., Wu, H.-J., Zhu, C.-W., Teng, F.-Z., Hoefs, J., Huang, J., Yang, W., Xu, L.-J., Hou, Z.-Q., Ren, Z.-Y., Li, S.-G. (2019) A nephelinitic component with unusual δ56Fe in Cenozoic basalts from eastern China and its implications for deep oxygen cycle. Earth and Planetary Science Letters 512, 175–183.
), have eventually triggered the GOE.top
Author Contributions
Y-SH designed the project. C-TL compiled the datasets and conducted the MATLAB coding. Both authors interpreted the results and prepared the manuscript.
top
Acknowledgements
The authors appreciate C.B. Keller for his generous and helpful suggestions about theory of the weighted bootstrap resampling method. K. Qu and Y.-Z. Ge provided help with coding. This work was supported by the National Science Foundation of China (Nos. 41688103 and 41730214), the Fundamental Research Funds for the Central Universities (3-7-5-2019-07), the 111 Project of the Ministry of Science and Technology, China (Grant No. BP0719021) and State Key Lab of Geological Processes and Mineral Resources. This is CUGB petrogeochemical contribution No. PGC-201565.
Editor: Cin-Ty Lee
top
References
Bindeman, I.N., Zakharov, D.O., Palandri, J., Greber, N.D., Dauphas, N., Retallack, G.J., Hofmann, A., Lackey, J.S., Bekker, A. (2018) Rapid emergence of subaerial landmasses and onset of a modern hydrologic cycle 2.5 billion years ago. Nature 557, 545–548.
Show in context
An increase in subaerial large igneous provinces (LIPs) and a stepwise change of trip-oxygen-isotope (Δ17O) composition of shales indicate major landmasses could have rapidly emerged above sea level at the Archean-Proterozoic boundary (Kump and Barley, 2007; Bindeman et al., 2018).
View in article
Therefore, the remarkably elevated mean K/La of pre-3.0 Ga mafic volcanic rocks reveals magma eruptions predominantly at a submarine environment (Fig. 3a), supporting the consensus that most early Archean landmasses were still below sea level (Bindeman et al., 2018).
View in article
Buick, R., Thornett, J.R., McNaughton, N.J., Smith, J.B., Barley, M.E., Savage, M. (1995) Record of emergent continental crust ∼3.5 billion years ago in the Pilbara craton of Australia. Nature 375, 574–577.
Show in context
Exposed landmasses have been traced back to as early as ∼3.5 Ga (e.g., detrital sediments in the East Pilbara Craton; Buick et al., 1995; Campbell and Davies, 2017), as indicated by geological records of uplift-subsidence cycles (e.g., in the East Pilbara Craton during ∼3.5–3.2 Ga; Campbell and Davies, 2017).
View in article
Campbell, I.H., Davies, D.R. (2017) Raising the continental crust. Earth and Planetary Science Letters 460, 112–122.
Show in context
Exposed landmasses have been traced back to as early as ∼3.5 Ga (e.g., detrital sediments in the East Pilbara Craton; Buick et al., 1995; Campbell and Davies, 2017), as indicated by geological records of uplift-subsidence cycles (e.g., in the East Pilbara Craton during ∼3.5–3.2 Ga; Campbell and Davies, 2017).
View in article
Exposed landmasses have been traced back to as early as ∼3.5 Ga (e.g., detrital sediments in the East Pilbara Craton; Buick et al., 1995; Campbell and Davies, 2017), as indicated by geological records of uplift-subsidence cycles (e.g., in the East Pilbara Craton during ∼3.5–3.2 Ga; Campbell and Davies, 2017).
View in article
Emergence of large subaerial landmasses initially occurred at ∼3.0 Ga at the Kaapvaal craton (e.g., Pongola-Wiwatersrand Basin; with an area of hundreds of square kilometres of exposed surface; Kump and Barley, 2007; Korenaga et al., 2017) that then extended to various cratons at ∼2.7 Ga (e.g., Kalgoorlie-Kambalda, Abitibi, Slave, Dharwar, Zimbabwe, Wyoming, Kaapvaal and Pilbara; Kump and Barley, 2007; Campbell and Davies, 2017).
View in article
Dhuime, B., Wuestefeld, A., Hawkesworth, C.J. (2015) Emergence of modern continental crust about 3 billion years ago. Nature Geoscience 8, 552–555.
Show in context
Secular changes in (b) normalised seawater 87Sr/86Sr curve (Shields et al., 2007), (c) the juvenile continental crustal thickness (Dhuime et al., 2015), (d) the upper continental crustal mass (Tang et al., 2016), and (e) atmospheric O2 level (Lyons et al., 2014) are also shown for comparison.
View in article
Rapid growth of crustal thickness as well as felsic components since ∼3.0 Ga (Tang et al., 2016; Dhuime et al., 2015) should have increased continental elevations.
View in article
Farmer, G.L. (2014) 4.3 – Continental Basaltic Rocks. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 75–110.
Show in context
On the other hand, recycling of dehydrated crustal materials can account for some intra-plate basalts with depleted K/La (e.g., average K/La values of Deccan Traps and Siberia LIPs are ∼201 and ∼161, respectively; Farmer, 2014).
View in article
Flament, N., Coltice, N., Rey, P.F. (2008) A case for late-Archaean continental emergence from thermal evolution models and hypsometry. Earth and Planetary Science Letters 275, 326–336.
Show in context
Static models considering multiple parameters (e.g., crustal thickness, density and mantle potential temperature) suggest the rise of major subaerial landmasses from the late Archean to the Neoproterozoic (Flament et al., 2008; Lee et al., 2017).
View in article
Gaillard, F., Scaillet, B., Arndt, N.T. (2011) Atmospheric oxygenation caused by a change in volcanic degassing pressure. Nature 478, 229–232.
Show in context
Most intriguingly, it was coincident with Earth’s first blast of molecular oxygen (i.e. the Great Oxidation Event (GOE)) and has been proposed as a key driving force for the GOE, probably by increasing oxidised subaerial volcanic gases (Kump and Barley, 2007; Gaillard et al., 2011) and decreasing O2 sinks (e.g., H2, CH4) of weathering and alteration (Smit and Mezger, 2017).
View in article
Intriguingly, the emergence of subaerial landmasses has been proposed as a key driving force for the GOE (Kump and Barley, 2007; Gaillard et al., 2011).
View in article
Hawkesworth, C.J., Cawood, P.A., Dhuime, B., Kemp, T.I.S. (2017) Earth’s continental lithosphere through time. Annual Review of Earth and Planetary Sciences 45, 169–198.
Show in context
Three transient positive excursions of K/La have been identified for post-2.7 Ga mafic volcanic rocks, i.e. at 0.4, 0.8 and 2.6 Ga, coincident with global tectono-magmatism lulls (Voice et al., 2011; Hawkesworth et al., 2017).
View in article
Global detrital zircon U-Pb age distribution (Voice et al., 2011) is plotted at the top of the figure, with the periods of supercontinents assembly (Hawkesworth et al., 2017).
View in article
Instead, the rise of major subaerial landmasses during ∼3.0–2.7 Ga can be explained by initiation of plate tectonics (∼3.0 Ga; Tang et al., 2016) and assembly of the Earth’s first supercontinent (Kenorland), which was marked by the first burst of felsic magmatism at ∼2.7–2.5 Ga (Voice et al., 2011; Hawkesworth et al., 2017; Fig. 2).
View in article
He, Y.-S., Meng, X.-N., Ke, S., Wu, H.-J., Zhu, C.-W., Teng, F.-Z., Hoefs, J., Huang, J., Yang, W., Xu, L.-J., Hou, Z.-Q., Ren, Z.-Y., Li, S.-G. (2019) A nephelinitic component with unusual δ56Fe in Cenozoic basalts from eastern China and its implications for deep oxygen cycle. Earth and Planetary Science Letters 512, 175–183.
Show in context
This study suggests alternative mechanisms, e.g., the evolution of oxygen photosynthesis, sedimentary burial of organic carbon (Lyons et al., 2014) and/or deep oxygen cycles driven by volatiles (He et al., 2019), have eventually triggered the GOE.
View in article
Hofmann, A., Harris, C. (2008) Silica alteration zones in the Barberton greenstone belt: A window into subseafloor processes 3.5–3.3 Ga ago. Chemical Geology 257, 221–239.
Show in context
With input from riverine runoff and high temperature (>100 or 150 °C) submarine hydrothermal processes (Hofmann and Harris, 2008; Staudigel, 2014), K is highly enriched relative to La in seawater.
View in article
For submarine magma eruptions, the K/La of altered mafic volcanic rocks can be strongly elevated during low temperature hydrothermal alteration (<100 or 150 °C) with seawater (e.g., hydrothermally altered basalts; Fig. 1) (Hofmann and Harris, 2008; Staudigel, 2014; Supplementary Information).
View in article
For instance, H2-rich fluids were likely to support the prosperity of methanogens and promote complex organic compounds (biotic or prebiotic life) in pristine seas (Hofmann and Harris, 2008).
View in article
Jagoutz, O. (2012) Were ancient granitoid compositions influenced by contemporaneous atmospheric and hydrosphere oxidation states? Terra Nova 25, 95–101.
Show in context
Potassium mobility could be quite different during anoxic versus oxic alteration of basalts (Jagoutz, 2012).
View in article
Kelemen, P.B., Hanghøj, K., Greene, A.R. (2014) 4.21 – One View of the Geochemistry of Subduction-Related Magmatic Arcs, with an Emphasis on Primitive Andesite and Lower Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 749–806.
Show in context
On the one hand, influx of more K relative to La into the mantle wedge via aqueous fluids derived from subducted sediments and oceanic slabs (Kelemen et al., 2014), tends to produce arc basalts with elevated K/La (e.g., average of continental arc basalts; Fig. 2a).
View in article
Modern geological reservoirs are shown as grey arrows with numbers: (1) upper continental crust, (2) continental arc basalts, (3) lower continental crust, (4) primitive mantle, (5) E-MORB, (6) ocean island basalts (OIB), (7) N-MORB (Sun and McDonough, 1989; Kelemen et al., 2014; Rudnick and Gao, 2014).
View in article
Keller, C.B., Schoene, B. (2012) Statistical geochemistry reveals disruption in secular lithospheric evolution about 2.5 Ga ago. Nature 485, 490–493.
Show in context
The weighted bootstrap re-sampling method with Monte Carlo simulations was applied to reveal secular changes in the K/La of less differentiated (SiO2: 43–51 wt. %) global continental mafic volcanic and plutonic rocks (Fig. S-1) (Keller and Schoene, 2012; Supplementary Information).
View in article
Korenaga, J., Planavsky, N.J., Evans, D.A.D. (2017) Global water cycle and the coevolution of the Earth’s interior and surface environment. Philosophical Transactions of the Royal Society A-Mathematical Physical and Engineering Sciences 375, 20150393.
Show in context
Emergence of large subaerial landmasses initially occurred at ∼3.0 Ga at the Kaapvaal craton (e.g., Pongola-Wiwatersrand Basin; with an area of hundreds of square kilometres of exposed surface; Kump and Barley, 2007; Korenaga et al., 2017) that then extended to various cratons at ∼2.7 Ga (e.g., Kalgoorlie-Kambalda, Abitibi, Slave, Dharwar, Zimbabwe, Wyoming, Kaapvaal and Pilbara; Kump and Barley, 2007; Campbell and Davies, 2017).
View in article
More than 50 % water in the hydrosphere has been suggested to be returned back to the mantle through subduction in the past (Korenaga et al., 2017).
View in article
Kump, L.R., Barley, M.E. (2007) Increased subaerial volcanism and the rise of atmospheric oxygen 2.5 billion years ago. Nature 448, 1033–1036.
Show in context
An increase in subaerial large igneous provinces (LIPs) and a stepwise change of trip-oxygen-isotope (Δ17O) composition of shales indicate major landmasses could have rapidly emerged above sea level at the Archean-Proterozoic boundary (Kump and Barley, 2007; Bindeman et al., 2018).
View in article
Most intriguingly, it was coincident with Earth’s first blast of molecular oxygen (i.e. the Great Oxidation Event (GOE)) and has been proposed as a key driving force for the GOE, probably by increasing oxidised subaerial volcanic gases (Kump and Barley, 2007; Gaillard et al., 2011) and decreasing O2 sinks (e.g., H2, CH4) of weathering and alteration (Smit and Mezger, 2017).
View in article
Emergence of large subaerial landmasses initially occurred at ∼3.0 Ga at the Kaapvaal craton (e.g., Pongola-Wiwatersrand Basin; with an area of hundreds of square kilometres of exposed surface; Kump and Barley, 2007; Korenaga et al., 2017) that then extended to various cratons at ∼2.7 Ga (e.g., Kalgoorlie-Kambalda, Abitibi, Slave, Dharwar, Zimbabwe, Wyoming, Kaapvaal and Pilbara; Kump and Barley, 2007; Campbell and Davies, 2017).
View in article
Intriguingly, the emergence of subaerial landmasses has been proposed as a key driving force for the GOE (Kump and Barley, 2007; Gaillard et al., 2011).
View in article
Laakso, T.A., Schrag, D.P. (2014) Regulation of atmospheric oxygen during the Proterozoic. Earth and Planetary Science Letters 388, 81–91.
Show in context
This time lag may not be explained by the delayed response of the surficial systems to evolution of the solid Earth, if a steady atmospheric oxygen level could be established at a time scale <50 Ma (Laakso and Schrag, 2014).
View in article
Lee, C.-T.A., Caves, J., Jiang, H.-H., Cao, W.-R., Lenardic, A., McKenzie, N.R., Shorttle, O., Yin, Q.-Z., Dyer, B. (2017) Deep mantle roots and continental emergence: implications for whole-earth elemental cycling, long-term climate, and the Cambrian explosion. International Geology Review 4, 431–448.
Show in context
Static models considering multiple parameters (e.g., crustal thickness, density and mantle potential temperature) suggest the rise of major subaerial landmasses from the late Archean to the Neoproterozoic (Flament et al., 2008; Lee et al., 2017).
View in article
Early subaerial landmasses may have experienced subsidence and eventually submerged below sea level again, due to crustal erosion and thermal relaxation of continental lithospheric mantle root (Lee et al., 2017).
View in article
The area of global subaerial landmasses in Earth’s history is regulated by water volume in the hydrosphere, topography differentiation driven by tectono-magmatic events (likely what happens at present day Tibet and Andes plateaus; Fig. 3c,d) as well as processes that eliminate topography difference, e.g., weathering erosion and subsidence due to thermal relaxation of continental lithospheric mantle root (Lee et al., 2017).
View in article
Crustal thickness plays an important role in controlling continental elevations (Lee et al., 2017).
View in article
However, static models with consideration of the crustal thickness, density and mantle potential temperature predict the rise of major subaerial landmasses until the Neoproterozoic (Lee et al., 2017).
View in article
Lyons, T.W., Reinhard, C.T., Planavsky, N.J. (2014) The rise of oxygen in Earth’s early ocean and atmosphere. Nature 506, 307–315.
Show in context
A “time lag” is identified here between the rise of major subaerial landmasses and the initiation of the GOE (Lyons et al., 2014; Fig. 2e).
View in article
This study suggests alternative mechanisms, e.g., the evolution of oxygen photosynthesis, sedimentary burial of organic carbon (Lyons et al., 2014) and/or deep oxygen cycles driven by volatiles (He et al., 2019), have eventually triggered the GOE.
View in article
Secular changes in (b) normalised seawater 87Sr/86Sr curve (Shields et al., 2007), (c) the juvenile continental crustal thickness (Dhuime et al., 2015), (d) the upper continental crustal mass (Tang et al., 2016), and (e) atmospheric O2 level (Lyons et al., 2014) are also shown for comparison.
View in article
Ma, J.-L., Wei, G.-J., Xu, Y.-G., Long, W.-G., Sun, W.-D. (2007) Mobilization and re-distribution of major and trace elements during extreme weathering of basalt in Hainan Island, South China. Geochimica et Cosmochimica Acta 71, 3223–3237.
Show in context
Weathered subaerial mafic rocks thus on average have depleted K/La (Ma et al., 2007).
View in article
Rapp, R.P., Shimizu, N., Norman, M.D. (2003) Growth of early continental crust by partial melting of eclogite. Nature 425, 605–609.
Show in context
Here we consider that the decline in the total water volume of the hydrosphere was minor during the rise of major subaerial landmasses, as Archean subduction zones, if they existed, were hot and deep subduction of water should have been prevented by prevailing slab melting (Rapp et al., 2003).
View in article
Reimink, J.R., Davies, J.H.F.L., Ielpi, A. (2021) Global zircon analysis records a gradual rise of continental crust throughout the Neoarchean. Earth and Planetary Science Letters 554, 116654.
Show in context
The earliest mature detrital zircon age distributions at ∼3.0–2.8 Ga also indicates widespread rise of subaerial continents at that time (Reimink et al., 2021).
View in article
Rudnick, R.L., Gao, S. (2014) 4.1 – Composition of the Continental Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 1–51.
Show in context
Modern geological reservoirs are shown as grey arrows with numbers: (1) upper continental crust, (2) continental arc basalts, (3) lower continental crust, (4) primitive mantle, (5) E-MORB, (6) ocean island basalts (OIB), (7) N-MORB (Sun and McDonough, 1989; Kelemen et al., 2014; Rudnick and Gao, 2014).
View in article
Further, mean K/La values of global mafic volcanic rocks at pre-2.7 Ga are too high (∼1039–2580; Table S-1) to be explained by crustal contamination, as the mean K/La of the modern upper continental crust is only ∼750 (Rudnick and Gao, 2014; Fig. 2a).
View in article
Shields, G.A. (2007) A normalised seawater strontium isotope curve: possible implications for Neoproterozoic-Cambrian weathering rates and the further oxygenation of the Earth. eEarth 2, 35–42.
Show in context
Secular changes in (b) normalised seawater 87Sr/86Sr curve (Shields et al., 2007), (c) the juvenile continental crustal thickness (Dhuime et al., 2015), (d) the upper continental crustal mass (Tang et al., 2016), and (e) atmospheric O2 level (Lyons et al., 2014) are also shown for comparison.
View in article
The mantle-like seawater 87Sr/86Sr compositions in pre-3.0 Ga (Shields, 2007; Fig. 2b) may also reflect limited erosion of proto-landmasses and a largely submarine state for early Archean landmasses.
View in article
Sleep, N.H., Bird, D.K., Pope, E.C. (2011) Serpentinite and the dawn of life. Philosophical Transactions of the Royal Society B-Biological Sciences 366, 2857–2869.
Show in context
Submarine landmasses universal in the early Archean, indicated by the extremely high mean K/La of mafic volcanic rocks (Figs. 2a, 3a), with widespread submarine hydrothermal processes and serpentinisation of ultramafic-mafic rocks could have buffered a reduced atmosphere (Smit et al., 2017) and provided proper environments (e.g., N2, H+, temperature) for early life (Sleep et al., 2011).
View in article
Smit, M.A., Mezger, K. (2017) Earth’s early O2 cycle suppressed by primitive continents. Nature Geoscience 10, 788–792.
Show in context
Most intriguingly, it was coincident with Earth’s first blast of molecular oxygen (i.e. the Great Oxidation Event (GOE)) and has been proposed as a key driving force for the GOE, probably by increasing oxidised subaerial volcanic gases (Kump and Barley, 2007; Gaillard et al., 2011) and decreasing O2 sinks (e.g., H2, CH4) of weathering and alteration (Smit and Mezger, 2017).
View in article
Submarine landmasses universal in the early Archean, indicated by the extremely high mean K/La of mafic volcanic rocks (Figs. 2a, 3a), with widespread submarine hydrothermal processes and serpentinisation of ultramafic-mafic rocks could have buffered a reduced atmosphere (Smit et al., 2017) and provided proper environments (e.g., N2, H+, temperature) for early life (Sleep et al., 2011).
View in article
Staudigel, H. (2014) 4.16 – Chemical Fluxes from Hydrothermal Alteration of the Oceanic Crust. In: Holland, H.D., Turekian, K.K. (Eds.) Treatise on Geochemistry. Second Edition, Elsevier, Oxford, 583–606.
Show in context
With input from riverine runoff and high temperature (>100 or 150 °C) submarine hydrothermal processes (Hofmann and Harris, 2008; Staudigel, 2014), K is highly enriched relative to La in seawater.
View in article
For submarine magma eruptions, the K/La of altered mafic volcanic rocks can be strongly elevated during low temperature hydrothermal alteration (<100 or 150 °C) with seawater (e.g., hydrothermally altered basalts; Fig. 1) (Hofmann and Harris, 2008; Staudigel, 2014; Supplementary Information).
View in article
Sun, S.-S., McDonough, W.F. (1989) Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes. Geological Society, London, Special Publications 42, 313–345.
Show in context
Mean normal mid-ocean ridge basalts (N-MORB) and enriched mid-ocean ridge basalts (E-MORB) are from Sun and McDonough (1989).
View in article
Modern geological reservoirs are shown as grey arrows with numbers: (1) upper continental crust, (2) continental arc basalts, (3) lower continental crust, (4) primitive mantle, (5) E-MORB, (6) ocean island basalts (OIB), (7) N-MORB (Sun and McDonough, 1989; Kelemen et al., 2014; Rudnick and Gao, 2014).
View in article
However, development of mantle reservoirs with differentiated K/La (e.g., OIB, N-MORB, E-MORB; Sun and McDonough, 1989) appears not to have significantly changed the mean K/La of global mafic magmas with time, as illustrated by nearly constant mean K/La of mafic plutons since 3.6 Ga (Fig. 2a).
View in article
Tang, M., Chen, K., Rudnick, R.L. (2016) Archean upper crust transition from mafic to felsic marks the onset of plate tectonics. Science 351, 372–375.
Show in context
Secular changes in (b) normalised seawater 87Sr/86Sr curve (Shields et al., 2007), (c) the juvenile continental crustal thickness (Dhuime et al., 2015), (d) the upper continental crustal mass (Tang et al., 2016), and (e) atmospheric O2 level (Lyons et al., 2014) are also shown for comparison.
View in article
With the development of plate tectonics (Tang et al., 2016), differentiation in subduction zones and recycling of crustal materials created substantial mantle heterogeneity.
View in article
Rapid growth of crustal thickness as well as felsic components since ∼3.0 Ga (Tang et al., 2016; Dhuime et al., 2015) should have increased continental elevations.
View in article
Instead, the rise of major subaerial landmasses during ∼3.0–2.7 Ga can be explained by initiation of plate tectonics (∼3.0 Ga; Tang et al., 2016) and assembly of the Earth’s first supercontinent (Kenorland), which was marked by the first burst of felsic magmatism at ∼2.7–2.5 Ga (Voice et al., 2011; Hawkesworth et al., 2017; Fig. 2).
View in article
Voice, P.J., Kowalewski, M., Eriksson, K.A. (2011) Quantifying the timing and rate of crustal evolution: global compilation of radiometrically dated detrital zircon grains. Journal of Geology 119, 109–126.
Show in context
Three transient positive excursions of K/La have been identified for post-2.7 Ga mafic volcanic rocks, i.e. at 0.4, 0.8 and 2.6 Ga, coincident with global tectono-magmatism lulls (Voice et al., 2011; Hawkesworth et al., 2017).
View in article
Global detrital zircon U-Pb age distribution (Voice et al., 2011) is plotted at the top of the figure, with the periods of supercontinents assembly (Hawkesworth et al., 2017).
View in article
Instead, the rise of major subaerial landmasses during ∼3.0–2.7 Ga can be explained by initiation of plate tectonics (∼3.0 Ga; Tang et al., 2016) and assembly of the Earth’s first supercontinent (Kenorland), which was marked by the first burst of felsic magmatism at ∼2.7–2.5 Ga (Voice et al., 2011; Hawkesworth et al., 2017; Fig. 2).
View in article
top
Supplementary Information
The Supplementary Information includes:
- Materials and Methods
- Tables S-1 to S-5
- Figures S-1 to S-3
- Supplementary Information References
Download Tables S-2 to S-5 (Excel).
Download the Supplementary Information (PDF).