Statistical geochemistry reveals disruption in secular lithospheric evolution about 2.5 Gyr ago

Subaerial weathering drove stabilization of continents

Earth's silica-rich continental crust is unique among the terrestrial planets and is critical for planetary habitability. Cratons represent the most imperishable continental fragments and form about 50% of the continental crust of the Earth, yet the mechanisms responsible for craton stabilization remain enigmatic1. Large tracts of strongly differentiated crust formed between 3 and 2.5 billion years ago, during the late Mesoarchaean and Neoarchaean time periods2. This crust contains abundant granitoid rocks with elevated concentrations of U, Th and K; the formation of these igneous rocks represents the final stage of stabilization of the continental crust2,3. Here, we show that subaerial weathering, triggered by the emergence of continental landmasses above sea level, facilitated intracrustal melting and the generation of peraluminous granitoid magmas. This resulted in reorganization of the compositional architecture of continental crust in the Neoarchaean period. Subaerial weathering concentrated heat-producing elements into terrigenous sediments that were incorporated into the deep crust, where they drove crustal melting and the chemical stratification required to stabilize the cratonic lithosphere. The chain of causality between subaerial weathering and the final differentiation of Earth's crust implies that craton stabilization was an inevitable consequence of continental emergence. Generation of sedimentary rocks enriched in heat-producing elements, at a time in the history of the Earth when the rate of radiogenic heat production was on average twice the present-day rate, resolves a long-standing question of why many cratons were stabilized in the Neoarchaean period.

Enhanced U-Pb detrital zircon, Lu-Hf zircon, δ18O zircon, and Sm-Nd whole rock global databases

High-quality global isotopic databases provide Earth scientists with robust means for developing and testing a variety of geological hypotheses. Database design establishes the range of questions that can be addressed, and validation techniques can enhance data quality. Here, six validated global isotopic databases provide extensive records of analyses from U-Pb in detrital zircon, Lu-Hf in zircon, Sm-Nd from whole rocks, and δ18O in zircon. The U-Pb detrital zircon records are segregated into three independently sampled databases. Independent samples are critical for testing the replicability of results, a key requisite for gaining confidence in the validity of a hypothesis. An advantage of our updated databases is that a hypothesis developed from one of the global detrital zircon databases can be immediately tested with the other two independent detrital zircon databases to assess the replicability of results. The independent εHf(t) and εNd(t) values provide similar means of testing for replicable results. This contribution discusses database design, data limitations, and validation techniques used to ensure the data are optimal for subsequent geological investigations.

Lithospheric thickness records tectonic evolution by controlling metamorphic conditions

The lithosphere, as the outermost solid layer of our planet, preserves a progressively more fragmentary record of geological events and processes from Earth's history the further back in time one looks. Thus, the evolution of lithospheric thickness and its cascading impacts in Earth's tectonic system are presently unknown. Here, we track the lithospheric thickness history using machine learning based on global lithogeochemical data of basalt. Our results demonstrate that four marked lithospheric thinning events occurred during the Paleoarchean, early Paleoproterozoic, Neoproterozoic, and Phanerozoic with intermediate thickening scenarios. These events respectively correspond to supercontinent/supercraton breakup and assembly periods. Causality investigation further indicates that crustal metamorphic and deformation styles are the feedback of lithospheric thickness. Cross-correlation between lithospheric thickness and metamorphic thermal gradients records the transition from intraoceanic subduction systems to continental margin and intraoceanic in the Paleoarchean and Mesoarchean and a progressive emergence of large thick continents that allow supercontinent growth, which promoted assembly of the first supercontinent during the Neoarchean.

Changes in orogenic style and surface environment recorded in Paleoproterozoic foreland successions

The Earth's interior and surficial systems underwent dramatic changes during the Paleoproterozoic, but the interaction between them remains poorly understood. Rocks deposited in orogenic foreland basins retain a record of the near surface to deep crustal processes that operate during subduction to collision and provide information on the interaction between plate tectonics and surface responses through time. Here, we document the depositional-to-deformational life cycle of a Paleoproterozoic foreland succession from the North China Craton. The succession was deposited in a foreland basin following ca. 2.50-2.47 Ga Altaid-style arc-microcontinent collision, and then converted to a fold-and-thrust belt at ca. 2.0-1.8 Ga due to Himalayan-style continent-continent collision. These two periods correspond to the assembly of supercratons in the late Archean and of the Paleoproterozoic supercontinent Columbia, respectively, which suggests that similar basins may have been common at the periphery of other cratons. The multiple stages of orogenesis and accompanying tectonic denudation and silicate weathering, as recorded by orogenic foreland basins, likely contributed to substantial changes in the hydrosphere, atmosphere, and biosphere known to have occurred during the Paleoproterozoic.

Influence of the asthenosphere on earth dynamics and evolution

The existence of a thin, weak asthenospheric layer beneath Earth's lithospheric plates is consistent with existing geological and geophysical constraints, including Pleistocene glacio-isostatic adjustment, modeling of gravity anomalies, studies of seismic anisotropy, and post-seismic rebound. Mantle convection models suggest that a pronounced weak zone beneath the upper thermal boundary layer (lithosphere) may be essential to the plate tectonic style of convection found on Earth. The asthenosphere is likely related to partial melting and the presence of water in the sub-lithospheric mantle, further implying that the long-term evolution of the Earth may be controlled by thermal regulation and volatile recycling that maintain a geotherm that approaches the wet mantle solidus at asthenospheric depths.

Interplay between oceanic subduction and continental collision in building continental crust

Generation of continental crust in collision zones reflect the interplay between oceanic subduction and continental collision. The Gangdese continental crust in southern Tibet developed during subduction of the Neo-Tethyan oceanic slab in the Mesozoic prior to reworking during the India-Asia collision in the Cenozoic. Here we show that continental arc magmatism started with fractional crystallization to form cumulates and associated medium-K calc-alkaline suites. This was followed by a period commencing at ~70 Ma dominated by remelting of pre-existing lower crust, producing more potassic compositions. The increased importance of remelting coincides with an acceleration in the convergence rate between India and Asia leading to higher basaltic flow into the Asian lithosphere, followed by convergence deceleration due to slab breakoff, enabling high heat flow and melting of the base of the arc. This two-stage process of accumulation and remelting leads to the chemical maturation of juvenile continental crust in collision zones, strengthening crustal stratification.

Reconstructing Earth's atmospheric oxygenation history using machine learning

Reconstructing historical atmospheric oxygen (O₂) levels at finer temporal resolution is a top priority for exploring the evolution of life on Earth. This goal, however, is challenged by gaps in traditionally employed sediment-hosted geochemical proxy data. Here, we propose an independent strategy-machine learning with global mafic igneous geochemistry big data to explore atmospheric oxygenation over the last 4.0 billion years. We observe an overall two-step rise of atmospheric O₂ similar to the published curves derived from independent sediment-hosted paleo-oxybarometers but with a more detailed fabric of O₂ fluctuations superimposed. These additional, shorter-term fluctuations are also consistent with previous but less well-established suggestions of O₂ variability. We conclude from this agreement that Earth's oxygenated atmosphere may therefore be at least partly a natural consequence of mantle cooling and specifically that evolving mantle melts collectively have helped modulate the balance of early O₂ sources and sinks.

Evidence of secular variation in Archean crust formation in the Eastern Indian Shield

Understanding the dominant crustal accretion model in any Archean craton is the key to understanding the dominant geodynamic process responsible for early crust formation during the Hadean (> 4.0 Ga) and Archaean (4.0–2.5 Ga). The continental crust has been proposed to have formed through either horizontal/vertical accretion related to subduction or mantle plume tectonic processes. Here, the Moho depths and average crustal Vp/Vs ratios are modelled at 16 broadband stations in the Eastern Indian Shield (EIS) through HK stacking of radial P-receiver functions (PRFs). These modelled parameters are used to test both plume and subduction models, which might have played a key role in the crustal accretion of the EIS throughout the Archean. We observe a correlation between crustal age and composition within the ellipsoidal Paleoarchean cratonic domain in the Singhbhum-Odisha-Craton (SOC), which reveals an increase in age from the younger granitoid core of the SOC (with thinning of felsic crust) to the surrounding older greenstone belts (with thickening of felsic crust). A thinner mafic crust resulting from multiple magmatic events characterizes the neighbouring Meso-Proterozoic Chotanagpur Granitic Gneissic terrain (CGGT). The Common Conversion Point (CCP) image of radial PRFs reveals northward subduction of the Paleoarchean SOC below the Meso-Proterozoic CGGT.

Dichotomy in crustal melting on early Mars inferred from antipodal effect

The Martian crustal dichotomy (MCD) between the southern highlands and the northern lowlands is the planet’s most ancient crustal structure, but its origins and evolution remain enigmatic. Understanding of the MCD comes largely from present-day and shallow crustal constraints. Lacking ancient and deeper constraints, hypotheses for the origin of the MCD range from an early giant impact, partial melting from sustained mantle convection, or some combination. We investigate with seismological modeling the best-preserved case of the “antipodal effect”—energy from an impact that concentrates and induces uplift and fracturing promoting volcanism at its antipode—the Hellas crater and the Alba Patera volcano on Mars. The volcano is latitudinally offset ∼2° (∼119 km) from the expected antipode, and we explore whether the MCD can explain this deflection. Variations across the MCD in topography, thickness, and composition have only minor effects. Simulations capable of sufficiently decelerating southern surface waves require the presence of 2%–5% more partial melt in the southern highlands. As the age of impact ca. 4 billion years ago post-dates the formation of the MCD, our partial melting results thus imply that, with or without an early giant impact, the MCD was modified by mantle convection in order to supply enough heat for crustal melts for several hundreds of millions of years after Mars formation.

•

We model the effect of impact-induced seismic waves causing volcanism at its antipode

•

A lower southern hemisphere crustal velocity explains the observed Martian antipodes

•

The simulation reveals a hemispheric dichotomy in crustal melting on early Mars

The Deep-Time Digital Earth program: data-driven discovery in geosciences

Current barriers hindering data-driven discoveries in deep-time Earth (DE) include: substantial volumes of DE data are not digitized; many DE databases do not adhere to FAIR (findable, accessible, interoperable and reusable) principles; we lack a systematic knowledge graph for DE; existing DE databases are geographically heterogeneous; a significant fraction of DE data is not in open-access formats; tailored tools are needed. These challenges motivate the Deep-Time Digital Earth (DDE) program initiated by the International Union of Geological Sciences and developed in cooperation with national geological surveys, professional associations, academic institutions and scientists around the world. DDE’s mission is to build on previous research to develop a systematic DE knowledge graph, a FAIR data infrastructure that links existing databases and makes dark data visible, and tailored tools for DE data, which are universally accessible. DDE aims to harmonize DE data, share global geoscience knowledge and facilitate data-driven discovery in the understanding of Earth's evolution.

Episodic growth of felsic continents in the past 3.7 Ga

Continents form the most accessible parts of Earth, but their complex compositions make their origin difficult to investigate. A novel approach based on a comprehensive compilation of samarium-neodymium isotopic compositions of detrital sedimentary rocks is here used to unravel continental growth through time. This record reveals that continents were as felsic as today in the past 3.7 Ga (billion years) and that their growth was not continuous but episodic. Reworking of preexisting crust was a ubiquitous process during most of Earth history, but at least six periods of continental growth can be identified every 500 to 700 Ma (million years) in the past 3.7 Ga. This recurrence could be accounted for by changes in tectonic plate velocities favoring periods of rapid subduction and enhanced production of juvenile felsic crust.

A long-term record of early to mid-Paleozoic marine redox change

The extent to which Paleozoic oceans differed from Neoproterozoic oceans and the causal relationship between biological evolution and changing environmental conditions are heavily debated. Here, we report a nearly continuous record of seafloor redox change from the deep-water upper Cambrian to Middle Devonian Road River Group of Yukon, Canada. Bottom waters were largely anoxic in the Richardson trough during the entirety of Road River Group deposition, while independent evidence from iron speciation and Mo/U ratios show that the biogeochemical nature of anoxia changed through time. Both in Yukon and globally, Ordovician through Early Devonian anoxic waters were broadly ferruginous (nonsulfidic), with a transition toward more euxinic (sulfidic) conditions in the mid-Early Devonian (Pragian), coincident with the early diversification of vascular plants and disappearance of graptolites. This ~80-million-year interval of the Paleozoic characterized by widespread ferruginous bottom waters represents a persistence of Neoproterozoic-like marine redox conditions well into the Phanerozoic.

Thermal state and evolving geodynamic regimes of the Meso- to Neoarchean North China Craton

Constraining thickness and geothermal gradient of Archean continental crust are crucial to understanding geodynamic regimes of the early Earth. Archean crust-sourced tonalitic–trondhjemitic–granodioritic gneisses are ideal lithologies for reconstructing the thermal state of early continental crust. Integrating experimental results with petrochemical data from the Eastern Block of the North China Craton allows us to establish temporal–spatial variations in thickness, geothermal gradient and basal heat flow across the block, which we relate to cooling mantle potential temperature and resultant changing geodynamic regimes from vertical tectonics in the late Mesoarchean (~2.9 Ga) to plate tectonics with hot subduction in the early to late Neoarchean (~2.7–2.5 Ga). Here, we show the transition to a plate tectonic regime plays an important role in the rapid cooling of the mantle, and thickening and strengthening of the lithosphere, which in turn prompted stabilization of the cratonic lithosphere at the end of the Archean.

Constraining the thermal state of the lithosphere is crucial to understanding geodynamic regime in early Earth. Here the authors reconstruct ~2.9–2.5 Ga thermal structure of continental lithosphere of the North China Craton using TTG and propose a systematic Archean geodynamic evolution process.

Heterogeneous Hadean crust with ambient mantle affinity recorded in detrital zircons of the Green Sandstone Bed, South Africa

The nature of Earth's earliest crust and the processes by which it formed remain major issues in Precambrian geology. Due to the absence of a rock record older than ∼4.02 Ga, the only direct record of the Hadean is from rare detrital zircon and that largely from a single area: the Jack Hills and Mount Narryer region of Western Australia. Here, we report on the geochemistry of Hadean detrital zircons as old as 4.15 Ga from the newly discovered Green Sandstone Bed in the Barberton greenstone belt, South Africa. We demonstrate that the U-Nb-Sc-Yb systematics of the majority of these Hadean zircons show a mantle affinity as seen in zircon from modern plume-type mantle environments and do not resemble zircon from modern continental or oceanic arcs. The zircon trace element compositions furthermore suggest magma compositions ranging from higher temperature, primitive to lower temperature, and more evolved tonalite-trondhjemite-granodiorite (TTG)-like magmas that experienced some reworking of hydrated crust. We propose that the Hadean parental magmas of the Green Sandstone Bed zircons formed from remelting of mafic, mantle-derived crust that experienced some hydrous input during melting but not from the processes seen in modern arc magmatism.

Remnants of early Earth differentiation in the deepest mantle-derived lavas

The noble gas isotope systematics of ocean island basalts suggest the existence of primordial mantle signatures in the deep mantle. Yet, the isotopic compositions of lithophile elements (Sr, Nd, Hf) in these lavas require derivation from a mantle source that is geochemically depleted by melt extraction rather than primitive. Here, this apparent contradiction is resolved by employing a compilation of the Sr, Nd, and Hf isotope composition of kimberlites-volcanic rocks that originate at great depth beneath continents. This compilation includes kimberlites as old as 2.06 billion years and shows that kimberlites do not derive from a primitive mantle source but sample the same geochemically depleted component (where geochemical depletion refers to ancient melt extraction) common to most oceanic island basalts, previously called PREMA (prevalent mantle) or FOZO (focal zone). Extrapolation of the Nd and Hf isotopic compositions of the kimberlite source to the age of Earth formation yields a ¹⁴³Nd/¹⁴⁴Nd-¹⁷⁶Hf/¹⁷⁷Hf composition within error of chondrite meteorites, which include the likely parent bodies of Earth. This supports a hypothesis where the source of kimberlites and ocean island basalts contains a long-lived component that formed by melt extraction from a domain with chondritic ¹⁴³Nd/¹⁴⁴Nd and ¹⁷⁶Hf/¹⁷⁷Hf shortly after Earth accretion. The geographic distribution of kimberlites containing the PREMA component suggests that these remnants of early Earth differentiation are located in large seismically anomalous regions corresponding to thermochemical piles above the core-mantle boundary. PREMA could have been stored in these structures for most of Earth's history, partially shielded from convective homogenization.

Recycled selenium in hot spot-influenced lavas records ocean-atmosphere oxygenation

Oxygenation of Earth's oceans and atmosphere through time has consequences for subducted surface signatures that are now stored in the mantle. Here, we report significant mass-dependent selenium isotope variations in modern hot spot-influenced oceanic lavas. These variations are correlated with tracers of mantle source enrichment, which can only be explained by incorporation of abyssal pelagic sediments subducted from a redox-stratified mid-Proterozoic ocean. Selenium geochemical signatures of these sediments have mostly been preserved during long-term recycling and may therefore complement the global surface sediment record as ancient oxygen archives. Combined deep mantle and surface perspectives, together with emerging models for atmospheric oxygen based on selenium systematics, further imply a significantly oxygenated ocean-atmosphere system throughout the mid-Proterozoic.

Constraining crustal silica on ancient Earth

Accurately quantifying the composition of continental crust on Hadean and Archean Earth is critical to our understanding of the physiography, tectonics, and climate of our planet at the dawn of life. One longstanding paradigm involves the growth of a relatively mafic planetary crust over the first 1 to 2 billion years of Earth history, implying a lack of modern plate tectonics and a paucity of subaerial crust, and consequently lacking an efficient mechanism to regulate climate. Others have proposed a more uniformitarian view in which Archean and Hadean continents were only slightly more mafic than at present. Apart from complications in assessing early crustal composition introduced by crustal preservation and sampling biases, effects such as the secular cooling of Earth's mantle and the biologically driven oxidation of Earth's atmosphere have not been fully investigated. We find that the former complicates efforts to infer crustal silica from compatible or incompatible element abundances, while the latter undermines estimates of crustal silica content inferred from terrigenous sediments. Accounting for these complications, we find that the data are most parsimoniously explained by a model with nearly constant crustal silica since at least the early Archean.

The evolution of the continental crust and the onset of plate tectonics

The Earth is the only known planet where plate tectonics is active, and different studies have concluded that plate tectonics commenced at times from the early Hadean to 700 Ma. Many arguments rely on proxies established on recent examples, such as paired metamorphic belts and magma geochemistry, and it can be difficult to establish the significance of such proxies in a hotter, older Earth. There is the question of scale, and how the results of different case studies are put in a wider global context. We explore approaches that indicate when plate tectonics became the dominant global regime, in part by evaluating when the effects of plate tectonics were established globally, rather than the first sign of its existence regionally. The geological record reflects when the continental crust became rigid enough to facilitate plate tectonics, through the onset of dyke swarms and large sedimentary basins, from relatively high-pressure metamorphism and evidence for crustal thickening. Paired metamorphic belts are a feature of destructive plate margins over the last 700 Myr, but it is difficult to establish whether metamorphic events are associated spatially as well as temporally in older terrains. From 3.8-2.7 Ga, suites of high Th/Nb (subduction-related on the modern Earth) and low Th/Nb (non-subduction-related) magmas were generated at similar times in different locations, and there is a striking link between the geochemistry and the regional tectonic style. Archaean cratons stabilised at different times in different areas from 3.1-2.5 Ga, and the composition of juvenile continental crust changed from mafic to more intermediate compositions. Xenon isotope data indicate that there was little recycling of volatiles before 3 Ga. Evidence for the juxtaposition of continental fragments back to ~2.8 Ga, each with disparate histories highlights that fragments of crust were moving around laterally on the Earth. The reduction in crustal growth at ~ 3 Ga is attributed to an increase in the rates at which differentiated continental crust was destroyed, and that coupled with the other changes at the end of the Archaean are taken to reflect the onset of plate tectonics as the dominant global regime.

Seismological evidence for the earliest global subduction network at 2 Ga ago

The earliest evidence for subduction, which could have been localized, does not signify when plate tectonics became a global phenomenon. To test the antiquity of global subduction, we investigated Paleoproterozoic time, for which seismic evidence is available from multiple continents. We used a new high-density seismic array in North China to image the crustal structure that exhibits a dipping Moho bearing close resemblance to that of the modern Himalaya. The relict collisional zone is Paleoproterozoic in age and implies subduction operating at least as early as ~2 billion years (Ga) ago. Seismic evidence of subduction from six continents at this age is interpreted as the oldest evidence of global plate tectonics. The sutures identified can be linked in a plate network that resulted in the assembly of Nuna, likely Earth's first supercontinent. Global subduction by ~2 Ga ago can explain why secular planetary cooling was not appreciable until Proterozoic time.

Geochemical evidence for a widespread mantle re-enrichment 3.2 billion years ago: implications for global-scale plate tectonics

Progressive mantle melting during the Earth’s earliest evolution led to the formation of a depleted mantle and a continental crust enriched in highly incompatible elements. Re-enrichment of Earth’s mantle can occur when continental crustal materials begin to founder into the mantle by either subduction or, to a lesser degree, by delamination processes, profoundly affecting the mantle’s trace element and volatile compositions. Deciphering when mantle re-enrichment/refertilization became a global-scale process would reveal the onset of efficient mass transfer of crust to the mantle and potentially when plate tectonic processes became operative on a global-scale. Here we document the onset of mantle re-enrichment/refertilization by comparing the abundances of petrogenetically significant isotopic values and key ratios of highly incompatible elements compared to lithophile elements in Archean to Early-Proterozoic mantle-derived melts (i.e., basalts and komatiites). Basalts and komatiites both record a rapid-change in mantle chemistry around 3.2 billion years ago (Ga) signifying a fundamental change in Earth geodynamics. This rapid-change is recorded in Nd isotopes and in key trace element ratios that reflect a fundamental shift in the balance between fluid-mobile and incompatible elements (i.e., Ba/La, Ba/Nb, U/Nb, Pb/Nd and Pb/Ce) in basaltic and komatiitic rocks. These geochemical proxies display a significant increase in magnitude and variability after ~3.2 Ga. We hypothesize that rapid increases in mantle heterogeneity indicate the recycling of supracrustal materials back into Earth’s mantle via subduction. Our new observations thus point to a ≥ 3.2 Ga onset of global subduction processes via plate tectonics.

Argon constraints on the early growth of felsic continental crust

The continental crust is a major geochemical reservoir, the evolution of which has shaped the surface environment of Earth. In this study, we present a new model of coupled crust-mantle-atmosphere evolution to constrain the growth of continental crust with atmospheric 40Ar/36Ar. Our model is the first to combine argon degassing with the thermal evolution of Earth in a self-consistent manner and to incorporate the effect of crustal recycling and reworking using the distributions of crustal formation and surface ages. Our results suggest that the history of argon degassing favors rapid crustal growth during the early Earth. The mass of continental crust, highly enriched in potassium, is estimated to have already reached >80% of the present-day level during the early Archean. The presence of such potassium-rich, likely felsic, crust has important implications for tectonics, surface environment, and the regime of mantle convection in the early Earth.

Global geochemical fingerprinting of plume intensity suggests coupling with the supercontinent cycle

Plate tectonics and mantle plumes are two of the most fundamental solid-Earth processes that have operated through much of Earth history. For the past 300 million years, mantle plumes are known to derive mostly from two large low shear velocity provinces (LLSVPs) above the core-mantle boundary, referred to as the African and Pacific superplumes, but their possible connection with plate tectonics is debated. Here, we demonstrate that transition elements (Ni, Cr, and Fe/Mn) in basaltic rocks can be used to trace plume-related magmatism through Earth history. Our analysis indicates the presence of a direct relationship between the intensity of plume magmatism and the supercontinent cycle, suggesting a possible dynamic coupling between supercontinent and superplume events. In addition, our analysis shows a consistent sudden drop in MgO, Ni and Cr at ~3.2–3.0 billion years ago, possibly indicating an abrupt change in mantle temperature at the start of global plate tectonics.

The links between plate tectonics and deep mantle structure remain unclear. Here, the authors demonstrate that transition elements (Ni, Cr, and Fe/Mn) in basaltic rocks can be used as a tool to trace plume-related magmatism through Earth history, and their results indicate the presence of a direct relationship between the intensity of plume magmatism and the supercontinent cycle.

Global atmospheric oxygen variations recorded by Th/U systematics of igneous rocks

Scientists have been working on the dating of atmospheric oxygenation in Earth’s history for decades. However, most previous studies relied on evidence from sedimentary rocks. Here, we show that igneous rocks can also be linked with surface oxidation by a key geological process: plate subduction. We here make an attempt to apply the Th/U of worldwide arc igneous rocks as an indicator for the timing of atmospheric oxygenation over the Earth’s history. The results are coincident with the previously defined Great Oxidation Event and Neoproterozoic Oxygenation Event.

Atmospheric oxygen has evolved from negligible levels in the Archean to the current level of about 21% through 2 major step rises: The Great Oxidation Event (GOE) in the early Proterozoic and the Neoproterozoic Oxygenation Event (NOE) during the late Proterozoic. However, most previous methods for constraining the time of atmospheric oxygenation have relied on evidence from sedimentary rocks. Here, we investigate the temporal variations of the Th/U of arc igneous rocks since 3.0 billion y ago (Ga) and show that 2 major Th/U decreases are recorded at ca. 2.35 Ga and ca. 0.75 Ga, coincident with the beginning of the GOE and NOE. The decoupling of U from Th is predominantly caused by the significant rise of atmospheric oxygen. Under an increasingly oxidized atmosphere condition, more uranium in the surface environment became oxidized from the water-insoluble U⁴⁺ to the water-soluble U⁶⁺ valance and incorporated in the sea water and altered oceanic crust. Eventually, the subduction of this altered oceanic crust produced the low-Th/U signature of arc igneous rocks. Therefore, the sharp decrease of Th/U in global arc igneous rocks may provide strong evidence for the rise of atmospheric oxygen. We suggest that the secular Th/U evolution of arc igneous rocks could be an effective geochemical indicator recording the global-scale atmospheric oxygen variation.

Evolution of cellular metabolism and the rise of a globally productive biosphere

Metabolic processes in cells and chemical processes in the environment are fundamentally intertwined and have evolved in concert for most of Earth's existence. Here I argue that intrinsic properties of cellular metabolism imposed central constraints on the historical trajectories of biopsheric productivity and atmospheric oxygenation. Photosynthesis depends on iron, but iron is highly insoluble under the aerobic conditions produced by oxygenic photosynthesis. These counteracting constraints led to two major stages of Earth oxygenation. After a cyanobacteria-driven biospheric expansion near the Archean-Proterozoic boundary, productivity remained largely restricted to continental boundaries and shallow aquatic environments where weathering inputs made iron more accessible. The anoxic deep open ocean was rich in free iron during the Proterozoic, but this iron was largely inaccessible, partly because an otherwise nutrient-poor ocean was limiting to photosynthesis, but also because a photosynthetic expansion would have quenched its own iron supply. Near the Proterozoic-Phanerozoic boundary, bioenergetics innovations allowed eukaryotic photosynthesis to overcome these interconnected negative feedbacks and begin expanding into the deep open oceans and onto the continents, where nutrients are inherently harder to come by. Key insights into what drove the ecological rise of eukaryotic photosynthesis emerge from analyses of marine Synechococcus and Prochlorococcus, abundant marine picocyanobacteria whose ancestors colonized the oceans in the Neoproterozoic. The reconstructed evolution of this group reveals a sequence of innovations that ultimately produced a form of photosynthesis in Prochlorococcus that is more like that of green plant cells than other cyanobacteria. Innovations increased the energy flux of cells, thereby enhancing their ability to acquire sparse nutrients, and as by-product also increased the production of organic carbon waste. Some of these organic waste products had the ability to chelate iron and make it bioavailable, thereby indirectly pushing the oceans through a transition from an anoxic state rich in free iron to an oxygenated state with organic carbon-bound iron. Resulting conditions (and parallel processes on the continents) in turn led to a series of positive feedbacks that increased the availability of other nutrients, thereby promoting the rise of a globally productive biosphere. In addition to the occurrence of major biospheric expansions, the several hundred million-year periods around the Archean-Proterozoic and Proterozoic-Phanerozoic boundaries share a number of other parallels. Both epochs have also been linked to major carbon cycle perturbations and global glaciations, as well as changes in the nature of plate tectonics and increases in continental exposure and weathering. This suggests the dynamics of life and Earth are intimately intertwined across many levels and that general principles governed transitions in these coupled dynamics at both times in Earth history.

Continuous plate subduction marked by the rise of alkali magmatism 2.1 billion years ago

Over the Earth’s evolutionary history, the style of plate subduction has evolved through time due to the secular cooling of the mantle. While continuous subduction is a typical feature of modern plate tectonics, a stagnant-lid tectonic regime with localized episodic subduction likely characterized the early Earth. The timing of the transition between these two subduction styles bears important insights into Earth’s cooling history. Here we apply a statistical analysis to a large geochemical dataset of mafic rocks spanning the last 3.5 Ga, which shows an increasing magnitude of alkali basaltic magmatism beginning at ca. 2.1 Ga. We propose that the rapid rise of continental alkali basalts correlates with an abruptly decreasing degree of mantle melting resulting from the enhanced cooling of the mantle at ca. 2.1 Ga. This might be a consequence of the initiation of continuous subduction, which recycled increasing volumes of cold oceanic crust into the mantle.

The point in time in which the ancient stagnant-lid tectonic regime with episodic plate subduction changed to continuous plate subduction is still topic of scientific debate. Here, the authors present vast geochemical dataset covering the last 3.5 billion years of Earth’s history and date the onset of continuous subduction to about 2.1 billion years ago.

The causes of spatiotemporal variations in erupted fluxes and compositions along a volcanic arc

Decades of study on volcanic arcs have provided insight into the overarching processes that control magmatism, and how these processes manifest at individual volcanoes. However, the causes of ubiquitous and dramatic intra-arc variations in volcanic flux and composition remain largely unresolved. Investigating such arc-scale issues requires greater quantitative comparison of geophysical and geochemical data, linked through sets of common intensive variables. To work towards these goals, we use observed lava compositions to estimate the heat budget associated with Quaternary volcanism in the Cascades Arc and compare this to the heat required to produce the observed geophysical properties of the crust. Here we show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily related to variations in the flux of basalt into the crust, rather than variations in their crustal storage history. This approach shows promise for studying other large-scale frontier geologic problems in volcanic arcs.

The primary causes of dramatic variations in volcanic flux and composition along strike in subduction zones remain largely unknown. Here we use a promising new approach to show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily due to variations in the flux of basalt into the base of the crust, rather than crustal magma storage.

Neoproterozoic glacial origin of the Great Unconformity

The Great Unconformity, a profound gap in Earth's stratigraphic record often evident below the base of the Cambrian system, has remained among the most enigmatic field observations in Earth science for over a century. While long associated directly or indirectly with the occurrence of the earliest complex animal fossils, a conclusive explanation for the formation and global extent of the Great Unconformity has remained elusive. Here we show that the Great Unconformity is associated with a set of large global oxygen and hafnium isotope excursions in magmatic zircon that suggest a late Neoproterozoic crustal erosion and sediment subduction event of unprecedented scale. These excursions, the Great Unconformity, preservational irregularities in the terrestrial bolide impact record, and the first-order pattern of Phanerozoic sedimentation can together be explained by spatially heterogeneous Neoproterozoic glacial erosion totaling a global average of 3-5 vertical kilometers, along with the subsequent thermal and isostatic consequences of this erosion for global continental freeboard.

Titanium isotopes as a tracer for the plume or island arc affinity of felsic rocks

Indirect evidence for the presence of a felsic continental crust, such as the elevated ⁴⁹Ti/⁴⁷Ti ratios in Archean shales, has been used to argue for ongoing subduction at that time and therefore plate tectonics. However, rocks of intermediate to felsic compositions can be produced in both plume and island arc settings. The fact that Ti behaves differently during magma differentiation in these two geological settings might result in contrasting isotopic signatures. Here, we demonstrate that, at a given SiO₂ content, evolved plume rocks (tholeiitic) are more isotopically fractionated in Ti than differentiated island arc rocks (mainly calc-alkaline). We also show that the erosion of crustal rocks from whether plumes (mafic in average) or island arcs (intermediate in average) can all produce sediments having quite constant ⁴⁹Ti/⁴⁷Ti ratios being 0.1-0.3 per mille heavier than that of the mantle. This suggests that Ti isotopes are not a direct tracer for the SiO₂ contents of crustal rocks. Ti isotopes in crustal sediments are still a potential proxy to identify the geodynamical settings for the formation of the crust but only if combined with additional SiO₂ information.

When Proterozoic Crusts Became Thick: New Insights from Magma Petrology

The Earth’s continental crust represents the outermost envelope of the solid Earth, controlling exchanges within the geosphere and reflecting geodynamics processes. One of the fundamental issues of Earth Science aims to determine crustal thickness in past geodynamic environments in order to discuss the evolution of certain geodynamic processes through time. Despite presenting a continuing challenge, the evolution of crustal thickness during the last 3 billion years can be investigated using indirect clues yielded by the chemical signature of mafic magmas and associated ferromagnesian minerals (pyroxene, amphibole). Here, we present a new statistical assessment of a global database of magmatic and mineral chemical information. Analysis reveals the increasing occurrence of high-temperature pyroxenes and amphiboles growing in Ca-rich, Fe-poor magma since ~1 Ga, which contrasts with lower temperature conditions of minerals crystallization throughout the Meso- and Palaeoproterozoic times. This is interpreted to reflect temporal changes in the control of Earth’s crust on mantle-derived magma composition, related to changes in lithospheric thickness and mantle secular cooling. We propose that thick existing crust is associated with deeper, hotter magmatic reservoirs, potentially elucidating the mineral chemistry and the contrasting iron content between primary and derivative mafic magmas. Based on both the chemical and mineral information of mafic magma, an integrated approach provides qualitative estimates of past crustal thickness and associated magmatic systems. Our findings indicate that the Proterozoic was characterized by thicker crustal sections (>40–50 km) relative to the Phanerozoic and Archean (<35 km). This period of crustal thickening appears at the confluence of major changes on Earth, marked by the onset of mantle cooling and Plate Tectonics and the assembly of Columbia, the first supercontinent.

The inception of plate tectonics: a record of failure

The development of plate tectonics from a pre-plate tectonics regime requires both the initiation of subduction and the development of nascent subduction zones into long-lived contiguous features. Subduction itself has been shown to be sensitive to system parameters such as thermal state and the specific rheology. While generally it has been shown that cold-interior high-Rayleigh-number convection (such as on the Earth today) favours plates and subduction, due to the ability of the interior stresses to couple with the lid, a given system may or may not have plate tectonics depending on its initial conditions. This has led to the idea that there is a strong history dependence to tectonic evolution-and the details of tectonic transitions, including whether they even occur, may depend on the early history of a planet. However, intrinsic convective stresses are not the only dynamic drivers of early planetary evolution. Early planetary geological evolution is dominated by volcanic processes and impacting. These have rarely been considered in thermal evolution models. Recent models exploring the details of plate tectonic initiation have explored the effect of strong thermal plumes or large impacts on surface tectonism, and found that these 'primary drivers' can initiate subduction, and, in some cases, over-ride the initial state of the planet. The corollary of this, of course, is that, in the absence of such ongoing drivers, existing or incipient subduction systems under early Earth conditions might fail. The only detailed planetary record we have of this development comes from Earth, and is restricted by the limited geological record of its earliest history. Many recent estimates have suggested an origin of plate tectonics at approximately 3.0 Ga, inferring a monotonically increasing transition from pre-plates, through subduction initiation, to continuous subduction and a modern plate tectonic regime around that time. However, both numerical modelling and the geological record itself suggest a strong nonlinearity in the dynamics of the transition, and it has been noted that the early history of Archaean greenstone belts and trondhjemite-tonalite-granodiorite record many instances of failed subduction. Here, we explore the history of subduction failure on the early Earth, and couple these with insights from numerical models of the geodynamic regime at the time.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

Geological archive of the onset of plate tectonics

Plate tectonics, involving a globally linked system of lateral motion of rigid surface plates, is a characteristic feature of our planet, but estimates of how long it has been the modus operandi of lithospheric formation and interactions range from the Hadean to the Neoproterozoic. In this paper, we review sedimentary, igneous and metamorphic proxies along with palaeomagnetic data to infer both the development of rigid lithospheric plates and their independent relative motion, and conclude that significant changes in Earth behaviour occurred in the mid- to late Archaean, between 3.2 Ga and 2.5 Ga. These data include: sedimentary rock associations inferred to have accumulated in passive continental margin settings, marking the onset of sea-floor spreading; the oldest foreland basin deposits associated with lithospheric convergence; a change from thin, new continental crust of mafic composition to thicker crust of intermediate composition, increased crustal reworking and the emplacement of potassic and peraluminous granites, indicating stabilization of the lithosphere; replacement of dome and keel structures in granite-greenstone terranes, which relate to vertical tectonics, by linear thrust imbricated belts; the commencement of temporally paired systems of intermediate and high dT/dP gradients, with the former interpreted to represent subduction to collisional settings and the latter representing possible hinterland back-arc settings or ocean plateau environments. Palaeomagnetic data from the Kaapvaal and Pilbara cratons for the interval 2780-2710 Ma and from the Superior, Kaapvaal and Kola-Karelia cratons for 2700-2440 Ma suggest significant relative movements. We consider these changes in the behaviour and character of the lithosphere to be consistent with a gestational transition from a non-plate tectonic mode, arguably with localized subduction, to the onset of sustained plate tectonics.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

When crust comes of age: on the chemical evolution of Archaean, felsic continental crust by crustal drip tectonics

The secular evolution of the Earth's crust is marked by a profound change in average crustal chemistry between 3.2 and 2.5 Ga. A key marker for this change is the transition from Archaean sodic granitoid intrusions of the tonalite-trondhjemite-granodiorite (TTG) series to potassic (K) granitic suites, akin (but not identical) to I-type granites that today are associated with subduction zones. It remains poorly constrained as to how and why this change was initiated and if it holds clues about the geodynamic transition from a pre-plate tectonic mode, often referred to as stagnant lid, to mobile plate tectonics. Here, we combine a series of proposed mechanisms for Archaean crustal geodynamics in a single model to explain the observed change in granitoid chemistry. Numeric modelling indicates that upper mantle convection drives crustal flow and subsidence, leading to profound diversity in lithospheric thickness with thin versus thick proto-plates. When convecting asthenospheric mantle interacts with lower lithosphere, scattered crustal drips are created. Under increasing P-T conditions, partial melting of hydrated meta-basalt within these drips produces felsic melts that intrude the overlying crust to form TTG. Dome structures, in which these melts can be preserved, are a positive diapiric expression of these negative drips. Transitional TTG with elevated K mark a second evolutionary stage, and are blends of subsided and remelted older TTG forming K-rich melts and new TTG melts. Ascending TTG-derived melts from asymmetric drips interact with the asthenospheric mantle to form hot, high-Mg sanukitoid. These melts are small in volume, predominantly underplated, and their heat triggered melting of lower crustal successions to form higher-K granites. Importantly, this evolution operates as a disseminated process in space and time over hundreds of millions of years (greater than 200 Ma) in all cratons. This focused ageing of the crust implies that compiled geochemical data can only broadly reflect geodynamic changes on a global or even craton-wide scale. The observed change in crustal chemistry does mark the lead up to but not the initiation of modern-style subduction.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

The dependence of planetary tectonics on mantle thermal state: applications to early Earth evolution

For plate tectonics to operate on a planet, mantle convective forces must be capable of forming weak, localized shear zones in the lithosphere that act as plate boundaries. Otherwise, a planet's mantle will convect in a stagnant lid regime, where subduction and plate motions are absent. Thus, when and how plate tectonics initiated on the Earth is intrinsically tied to the ability of mantle convection to form plate boundaries; however, the physics behind this process are still uncertain. Most mantle convection models have employed a simple pseudoplastic model of the lithosphere, where the lithosphere 'fails' and develops a mobile lid when stresses in the lithosphere reach the prescribed yield stress. With pseudoplasticity high mantle temperatures and high rates of internal heating, conditions relevant for the early Earth, impede plate boundary formation by decreasing lithospheric stresses, and hence favour a stagnant lid for the early Earth. However, when a model for shear zone formation based on grain size reduction is used, early Earth thermal conditions do not favour a stagnant lid. While lithosphere stress drops with increasing mantle temperature or heat production rate, the deformational work, which drives grain size reduction, increases. Thus, the ability of convection to form weak plate boundaries is not impeded by early Earth thermal conditions. However, mantle thermal state does change the style of subduction and lithosphere mobility; high mantle temperatures lead to a more sluggish, drip-like style of subduction. This 'sluggish lid' convection may be able to explain many of the key observations of early Earth crust formation processes preserved in the geologic record. Moreover, this work highlights the importance of understanding the microphysics of plate boundary formation for assessing early Earth tectonics, as different plate boundary formation mechanisms are influenced by mantle thermal state in fundamentally different ways.This article is part of a discussion meeting issue 'Earth dynamics and the development of plate tectonics'.

A non-zircon Hf isotope record in Archean black shales from the Pilbara craton confirms changing crustal dynamics ca. 3 Ga ago

Plate tectonics and associated subduction are unique to the Earth. Studies of Archean rocks show significant changes in composition and structural style around 3.0 to 2.5 Ga that are related to changing tectonic regime, possibly associated with the onset of subduction. Whole rock Hf isotope systematics of black shales from the Australian Pilbara craton, selected to exclude detrital zircon components, are employed to evaluate the evolution of the Archean crust. This approach avoids limitations of Hf-in-zircon analyses, which only provide input from rocks of sufficient Zr-concentration, and therefore usually represent domains that already underwent a degree of differentiation. In this study, we demonstrate the applicability of this method through analysis of shales that range in age from 3.5 to 2.8 Ga, and serve as representatives of their crustal sources through time. Their Hf isotopic compositions show a trend from strongly positive εHf_initial values for the oldest samples, to strongly negative values for the younger samples, indicating a shift from juvenile to differentiated material. These results confirm a significant change in the character of the source region of the black shales by 3 Ga, consistent with models invoking a change in global dynamics from crustal growth towards crustal reworking around this time.

A record of deep-ocean dissolved O2 from the oxidation state of iron in submarine basalts

The oxygenation of the deep ocean in the geological past has been associated with a rise in the partial pressure of atmospheric molecular oxygen (O₂) to near-present levels and the emergence of modern marine biogeochemical cycles. It has also been linked to the origination and diversification of early animals. It is generally thought that the deep ocean was largely anoxic from about 2,500 to 800 million years ago, with estimates of the occurrence of deep-ocean oxygenation and the linked increase in the partial pressure of atmospheric oxygen to levels sufficient for this oxygenation ranging from about 800 to 400 million years ago. Deep-ocean dissolved oxygen concentrations over this interval are typically estimated using geochemical signatures preserved in ancient continental shelf or slope sediments, which only indirectly reflect the geochemical state of the deep ocean. Here we present a record that more directly reflects deep-ocean oxygen concentrations, based on the ratio of Fe³⁺ to total Fe in hydrothermally altered basalts formed in ocean basins. Our data allow for quantitative estimates of deep-ocean dissolved oxygen concentrations from 3.5 billion years ago to 14 million years ago and suggest that deep-ocean oxygenation occurred in the Phanerozoic (541 million years ago to the present) and potentially not until the late Palaeozoic (less than 420 million years ago).

Geochemical and mineralogical evidence that Rodinian assembly was unique

The mineralogy and geochemistry associated with Rodinian assembly (~1.3-0.9 Ga) are significantly different from those of other supercontinents. Compared to other supercontinents, relatively more Nb-bearing minerals, Y-bearing minerals, and zircons formed during Rodinian assembly, with corresponding enrichments of Nb, Y, and Zr concentrations in igneous rocks. By contrast, minerals bearing many other elements (e.g., Ni, Co, Au, Se, and platinum group elements) are significantly less abundant, without corresponding depletion of Ni and Co concentrations in igneous rocks. Here we suggest that the Nb, Y, and Zr enrichments in igneous rocks and relatively more occurrences of corresponding Nb-bearing minerals, Y-bearing minerals, and zircons result from significant non-arc magmatism during the mid-Proterozoic, while fewer occurrences of many other minerals suggest enhanced erosion of Rodinian volcanic arcs and orogens. The prolonged, extrovert assembly of Rodinia from thickened mid-Proterozoic continental crust via two-sided subduction can account for both the prevalence of non-arc magmatism and the enhanced erosion.

Recycled iron fuels new production in the eastern equatorial Pacific Ocean

Nitrate persists in eastern equatorial Pacific surface waters because phytoplankton growth fueled by nitrate (new production) is limited by iron. Nitrate isotope measurements provide a new constraint on the controls of surface nitrate concentration in this region and allow us to quantify the degree and temporal variability of nitrate consumption. Here we show that nitrate consumption in these waters cannot be fueled solely by the external supply of iron to these waters, which occurs by upwelling and dust deposition. Rather, a substantial fraction of nitrate consumption must be supported by the recycling of iron within surface waters. Given plausible iron recycling rates, seasonal variability in nitrate concentration on and off the equator can be explained by upwelling rate, with slower upwelling allowing for more cycles of iron regeneration and uptake. The efficiency of iron recycling in the equatorial Pacific implies the evolution of ecosystem-level mechanisms for retaining iron in surface ocean settings where it limits productivity.

The Neoarchaean surficial sulphur cycle: An alternative hypothesis based on analogies with 20th-century atmospheric lead

We revisit the S-isotope systematics of sedimentary pyrite in a shaly limestone from the ca. 2.52 Ga Gamohaan Formation, Upper Campbellrand Subgroup, Transvaal, South Africa. The analysed rock is interpreted to have been deposited in a water depth of ca. 50-100 m, in a restricted sub-basin on a drowning platform. A previous study discovered that the pyrites define a nonzero intercept δ³⁴ S_V-CDT -Δ³³ S data array. The present study carried out further quadruple S-isotope analyses of pyrite, confirming and expanding the linear δ³⁴ S_V-CDT -Δ³³ S array with an δ³⁴ S zero intercept at ∆³³ S ca. +5. This was previously interpreted to indicate mixing of unrelated S-sources in the sediment environment, involving a combination of recycled sulphur from sulphides that had originally formed by sulphate-reducing bacteria, along with elemental sulphur. Here, we advance an alternative explanation based on the recognition that the Archaean seawater sulphate concentration was likely very low, implying that the Archaean ocean could have been poorly mixed with respect to sulphur. Thus, modern oceanic sulphur systematics provide limited insight into the Archaean sulphur cycle. Instead, we propose that the 20th-century atmospheric lead event may be a useful analogue. Similar to industrial lead, the main oceanic input of Archaean sulphur was through atmospheric raindown, with individual giant point sources capable of temporally dominating atmospheric input. Local atmospheric S-isotope signals, of no global significance, could thus have been transmitted into the localised sediment record. Thus, the nonzero intercept δ³⁴ S_V-CDT -Δ³³ S data array may alternatively represent a very localised S-isotope signature in the Neoarchaean surface environment. Fallout from local volcanic eruptions could imprint recycled MIF-S signals into pyrite of restricted depositional environments, thereby avoiding attenuation of the signal in the subdued, averaged global open ocean sulphur pool. Thus, the superposition of extreme local S-isotope signals offers an alternative explanation for the large Neoarchaean MIF-S excursions and asymmetry of the Δ³³ S rock record.

Continental igneous rock composition: A major control of past global chemical weathering

The composition of igneous rocks in the continental crust has changed throughout Earth's history. However, the impact of these compositional variations on chemical weathering, and by extension on seawater and atmosphere evolution, is largely unknown. We use the strontium isotope ratio in seawater [(87Sr/86Sr)seawater] as a proxy for chemical weathering, and we test the sensitivity of (87Sr/86Sr)seawater variations to the strontium isotopic composition (87Sr/86Sr) in igneous rocks generated through time. We demonstrate that the 87Sr/86Sr ratio in igneous rocks is correlated to the epsilon hafnium (εHf) of their hosted zircon grains, and we use the detrital zircon record to reconstruct the evolution of the 87Sr/86Sr ratio in zircon-bearing igneous rocks. The reconstructed 87Sr/86Sr variations in igneous rocks are strongly correlated with the (87Sr/86Sr)seawater variations over the last 1000 million years, suggesting a direct control of the isotopic composition of silicic magmatism on (87Sr/86Sr)seawater variations. The correlation decreases during several time periods, likely reflecting changes in the chemical weathering rate associated with paleogeographic, climatic, or tectonic events. We argue that for most of the last 1000 million years, the (87Sr/86Sr)seawater variations are responding to changes in the isotopic composition of silicic magmatism rather than to changes in the global chemical weathering rate. We conclude that the (87Sr/86Sr)seawater variations are of limited utility to reconstruct changes in the global chemical weathering rate in deep times.

Diversity of burial rates in convergent settings decreased as Earth aged

The evolution and the growth of the continental crust is inextricably linked to the evolution of Earth's geodynamic processes. The detrital zircon record within the continental crust, as well as the isotopic composition of this crust, indicates that the amount of juvenile felsic material decreased with time and that in geologically recent times, the generation of new crust is balanced by recycling of the crust back into the mantle within subduction zones. However it cannot always have been so; yet the nature of the crust and the processes of crustal reworking in the Precambrian Earth are not well constrained. Here we use both detrital zircon ages and metamorphic pressure-temperature-time (P-T-t) information from metasedimentary units deposited in proposed convergent settings from Archaean, Proterozoic and Phanerozoic terrains to characterize the evolution of minimum estimates of burial rate (km.Ma(-1)) as a function of the age of the rocks. The demonstrated decrease in burial rate correlates positively with a progressive decrease in the production of juvenile felsic crust in the Archaean and Proterozoic. Burial rates are also more diverse in the Archaean than in modern times. We interpret these features to reflect a progressive decrease in the diversity of tectonic processes from Archaean to present, coupled with the emergence of the uniquely Phanerozoic modern-style collision.

Crustal thickness control on Sr/Y signatures of recent arc magmas: an Earth scale perspective

Arc magmas originate in subduction zones as partial melts of the mantle, induced by aqueous fluids/melts liberated by the subducted slab. Subsequently, they rise through and evolve within the overriding plate crust. Aside from broadly similar features that distinguish them from magmas of other geodynamic settings (e.g., mid-ocean ridges, intraplate), arc magmas display variably high Sr/Y values. Elucidating the debated origin of high Sr/Y signatures in arc magmas, whether due to mantle-source, slab melting or intracrustal processes, is instrumental for models of crustal growth and ore genesis. Here, using a statistical treatment of >23000 whole rock geochemical data, I show that average Sr/Y values and degree of maturation (MgO depletion at peak Sr/Y values) of 19 out of 22 Pliocene-Quaternary arcs correlate positively with arc thickness. This suggests that crustal thickness exerts a first order control on the Sr/Y variability of arc magmas through the stabilization or destabilization of mineral phases that fractionate Sr (plagioclase) and Y (amphibole ± garnet). In fact, the stability of these mineral phases is function of the pressure at which magma evolves, which depends on crustal thickness. The data presented show also that high Sr/Y Pliocene-Quaternary intermediate-felsic arc rocks have a distinct origin from their Archean counterparts.

The zircon archive of continent formation through time

The strong resilience of the mineral zircon and its ability to host a wealth of isotopic information make it the best deep-time archive of Earth's continental crust. Zircon is found in most felsic igneous rocks, can be precisely dated and can fingerprint magmatic sources; thus, it has been widely used to document the formation and evolution of continental crust, from pluton- to global-scale. Here, we present a review of major contributions that zircon studies have made in terms of understanding key questions involving the formation of the continents. These include the conditions of continent formation on early Earth, the onset of plate tectonics and subduction, the rate of crustal growth through time and the governing balance of continental addition v. continental loss, and the role of preservation bias in the zircon record.

Supplementary material:

A compilation used in this study of previously published detrital zircon U-Pb-Hf isotope data are available at http://www.geolsoc.org.uk/SUP18791