Thermochemical lithosphere differentiation and the origin of cratonic mantle.

Cratons record the early history of continental lithosphere formation, yet how they became the most enduring part of the lithosphere on Earth remains unknown1. Here we propose a mechanism for the formation of large volumes of melt-depleted cratonic lithospheric mantle (CLM) and its evolution to stable cratons. Numerical models show large decompression melting of a hot, early Earth mantle beneath a stretching lithosphere, where melt extraction leaves large volumes of depleted mantle at depth. The dehydrated, stiffer mantle resists further deformation, forcing strain migration and cooling, thereby assimilating depleted mantle into the lithosphere. The negative feedback between strain localization and stiffening sustains long-term diffused extension and emplacement of large amounts of depleted CLM. The formation of CLM at low pressure and its deeper re-equilibration reproduces the evolution of Archaean lithosphere constrained by depth-temperature conditions1,2, whereas large degrees of depletion3,4 and melt volumes5 in Archaean cratons are best matched by models with lower lithospheric strength. Under these conditions, which are otherwise viable for plate tectonics6,7, thermochemical differentiation effectively prevents yielding and formation of margins: rifting and lithosphere subduction are short lived and embedded in the cooling CLM as relict structures, reproducing the recycling and reworking environments that are found in Archaean cratons8,9. Although they undergo major melting and extensive recycling during an early stage lasting approximately 500 million years, the modelled lithospheres progressively differentiate and stabilize, and then recycling and reworking become episodic. Early major melting and recycling events explain the production and loss of primordial Hadean lithosphere and crust10, whereas later stabilization and episodic reworking provides a context for the creation of continental cratons in the Archaean era4,8.

Magmatic thickening of crust in non-plate tectonic settings initiated the subaerial rise of Earth's first continents 3.3 to 3.2 billion years ago.

When and how Earth's earliest continents-the cratons-first emerged above the oceans (i.e., emersion) remain uncertain. Here, we analyze a craton-wide record of Paleo-to-Mesoarchean granitoid magmatism and terrestrial to shallow-marine sedimentation preserved in the Singhbhum Craton (India) and combine the results with isostatic modeling to examine the timing and mechanism of one of the earliest episodes of large-scale continental emersion on Earth. Detrital zircon U-Pb(-Hf) data constrain the timing of terrestrial to shallow-marine sedimentation on the Singhbhum Craton, which resolves the timing of craton-wide emersion. Time-integrated petrogenetic modeling of the granitoids quantifies the progressive changes in the cratonic crustal thickness and composition and the pressure-temperature conditions of granitoid magmatism, which elucidates the underlying mechanism and tectonic setting of emersion. The results show that the entire Singhbhum Craton became subaerial ∼3.3 to 3.2 billion years ago (Ga) due to progressive crustal maturation and thickening driven by voluminous granitoid magmatism within a plateau-like setting. A similar sedimentary-magmatic evolution also accompanied the early (>3 Ga) emersion of other cratons (e.g., Kaapvaal Craton). Therefore, we propose that the emersion of Earth's earliest continents began during the late Paleoarchean to early Mesoarchean and was driven by the isostatic rise of their magmatically thickened (∼50 km thick), buoyant, silica-rich crust. The inferred plateau-like tectonic settings suggest that subduction collision-driven compressional orogenesis was not essential in driving continental emersion, at least before the Neoarchean. We further surmise that this early emersion of cratons could be responsible for the transient and localized episodes of atmospheric-oceanic oxygenation (O2-whiffs) and glaciation on Archean Earth.

Rates of generation and destruction of the continental crust: implications for continental growth.

Less than 25% of the volume of the juvenile continental crust preserved today is older than 3 Ga, there are no known rocks older than approximately 4 Ga, and yet a number of recent models of continental growth suggest that at least 60-80% of the present volume of the continental crust had been generated by 3 Ga. Such models require that large volumes of pre-3 Ga crust were destroyed and replaced by younger crust since the late Archaean. To address this issue, we evaluate the influence on the rock record of changing the rates of generation and destruction of the continental crust at different times in Earth's history. We adopted a box model approach in a numerical model constrained by the estimated volumes of continental crust at 3 Ga and the present day, and by the distribution of crust formation ages in the present-day crust. The data generated by the model suggest that new continental crust was generated continuously, but with a marked decrease in the net growth rate at approximately 3 Ga resulting in a temporary reduction in the volume of continental crust at that time. Destruction rates increased dramatically around 3 billion years ago, which may be linked to the widespread development of subduction zones. The volume of continental crust may have exceeded its present value by the mid/late Proterozoic. In this model, about 2.6-2.3 times of the present volume of continental crust has been generated since Earth's formation, and approximately 1.6-1.3 times of this volume has been destroyed and recycled back into the mantle.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'.

Zircons reveal the history of fluctuations in oxidation state of crustal magmatism and supercontinent cycle.

We apply a zircon redox index to a global compilation of detrital zircons to track the variation of oxidation state, expressed as ΔFMQ, through Earth's history. Those from I-type rocks, which comprise mantle and crustal igneous protoliths, including tonalite-trondhjemite-granodiorites (TTGs), generally have a high oxidation state (ΔFMQ > 0). In contrast, zircons from igneous rocks derived from supracrustal source rocks (S-type) are commonly reduced (ΔFMQ < 0). With the probability density function derived from the Gaussian-Kernel-Density-Estimation, we use the maximum likelihood estimation (MLE) to distinguish S-type from I-type zircons through Earth's history using zircon redox. Voluminous S-type magma production shows a ca. 600 Ma cyclicity that is closely related to the supercontinent cycle. We link a cyclic drop in redox values after 2.6 Ga to periodic S-type magma generation associated with burial and melting of metasedimentary rocks during supercontinent assembly and amalgamation. The ΔFMQ of the detrital zircons rise at ∼3.5 Ga followed by a consistent average ΔFMQ > 0 over the last 3 Ga. Given that the redox state of magmas is independent of crustal thickness and silica variation, and elevated values are likely more closely related to tectonic setting, we suggest that the consistent average ΔFMQ > 0 from ca. 3.5 Ga onwards relates to recycling of oceanic lithosphere back into the mantle in what eventually became established as subduction zones. The more reduced magmas associated with sedimentary sources, became established at 2.6 Ga, presumably in response to continental rocks rising above sea-level, and follow peaks of productivity associated with the supercontinent cycle.

A resource for automated search and collation of geochemical datasets from journal supplements.

This article presents a resource for automated search, extraction and collation of geochemical and geochronological data from the Figshare repository using web scraping code. To answer fundamental questions about the Earth's evolution, such as spatial and temporal evolution and interrelationships between the planet's solid and surficial reservoirs, researchers must utilize global geochemical datasets. Due to the volume of data being published, these datasets become quickly outdated. We present a resource that allows researchers to rapidly curate and update their own databases from existing published data. We use open-source Python code to web scrape the Figshare repository for journal supplementary files using the application programming interface, allowing for the collection and download of hundreds of supplementary files and metadata in minutes. Use of this web scraping tool is demonstrated here by collation of a zircon geochronology and chemistry database of >150,000 analyses. The database is consistent in reproducing trends in other published zircon compilations. Providing a resource for automated collection of Figshare data files will encourage data sharing and reuse.

Oxidation of Archean upper mantle caused by crustal recycling.

The redox evolution of Archean upper mantle impacted mantle melting and the nature of chemical equilibrium between mantle, ocean and atmosphere of the early Earth. Yet, the origin of these variations in redox remain controversial. Here we show that a global compilation of ∼3.8-2.5 Ga basalts can be subdivided into group B-1, showing modern mid-ocean ridge basalt-like features ((Nb/La)PM ≥ 0.75), and B-2, which are similar to contemporary island arc-related basalts ((Nb/La)PM < 0.75). Our V-Ti redox proxy indicates a more reducing upper mantle, and the results of both ambient and modified mantle obtained from B-1 and B-2 samples, respectively, exhibit a ∼1.0 log unit increase in their temporal evolution for most cratons. Increases in mantle oxygen fugacity are coincident with the changes in basalt Th/Nb ratios and Nd isotope ratios, indicating that crustal recycling played a crucial role, and this likely occurred either via plate subduction or lithospheric drips.

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.

Mariana-type ophiolites constrain the establishment of modern plate tectonic regime during Gondwana assembly.

Initiation of Mariana-type oceanic subduction zones requires rheologically strong oceanic lithosphere, which developed through secular cooling of Earth's mantle. Here, we report a 518 Ma Mariana-type subduction initiation ophiolite from northern Tibet, which, along with compilation of similar ophiolites through Earth history, argues for the establishment of the modern plate tectonic regime by the early Cambrian. The ophiolite was formed during the subduction initiation of the Proto-Tethys Ocean that coincided with slab roll-back along the southern and western Gondwana margins at ca. 530-520 Ma. This global tectonic re-organization and the establishment of modern plate tectonic regime was likely controlled by secular cooling of the Earth, and facilitated by enhanced lubrication of subduction zones by sediments derived from widespread surface erosion of the extensive mountain ranges formed during Gondwana assembly. This time also corresponds to extreme events recorded in climate and surface proxies that herald formation of the contemporary Earth.

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.

Crustal rejuvenation stabilised Earth's first cratons.

The formation of stable, evolved (silica-rich) crust was essential in constructing Earth's first cratons, the ancient nuclei of continents. Eoarchaean (4000-3600 million years ago, Ma) evolved crust occurs on most continents, yet evidence for older, Hadean evolved crust is mostly limited to rare Hadean zircons recycled into younger rocks. Resolving why the preserved volume of evolved crust increased in the Eoarchaean is key to understanding how the first cratons stabilised. Here we report new zircon uranium-lead and hafnium isotope data from the Yilgarn Craton, Australia, which provides an extensive record of Hadean-Eoarchaean evolved magmatism. These data reveal that the first stable, evolved rocks in the Yilgarn Craton formed during an influx of juvenile (recently extracted from the mantle) magmatic source material into the craton. The concurrent shift to juvenile sources and onset of crustal preservation links craton stabilisation to the accumulation of enduring rafts of buoyant, melt-depleted mantle.

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.

An Early Cretaceous subduction-modified mantle underneath the ultraslow spreading Gakkel Ridge, Arctic Ocean.

Earth's upper mantle, as sampled by mid-ocean ridge basalts (MORBs) at oceanic spreading centers, has developed chemical and isotopic heterogeneity over billions of years through focused melt extraction and re-enrichment by recycled crustal components. Chemical and isotopic heterogeneity of MORB is dwarfed by the large compositional spectrum of lavas at convergent margins, identifying subduction zones as the major site for crustal recycling into and modification of the mantle. The fate of subduction-modified mantle and if this heterogeneity transmits into MORB chemistry remains elusive. Here, we investigate the origin of upper mantle chemical heterogeneity underneath the Western Gakkel Ridge region in the Arctic Ocean through MORB geochemistry and tectonic plate reconstruction. We find that seafloor lavas from the Western Gakkel Ridge region mirror geochemical signatures of an Early Cretaceous, paleo-subduction zone, and conclude that the upper mantle can preserve a long-lived, stationary geochemical memory of past geodynamic processes.

Crustal growth and reworking: A case study from the Erguna Massif, eastern Central Asian Orogenic Belt.

Despite being the largest accretionary orogen on Earth, the record of crustal growth and reworking of individual microcontinental massifs within the Central Asian Orogenic Belt (CAOB) remain poorly constrained. Here, we focus on zircon records from granitoids in the Erguna Massif to discuss its crustal evolution through time. Proterozoic-Mesozoic granitoids are widespread in the Erguna Massif, and spatiotemporal variations in their zircon εHf(t) values and TDM2(Hf) ages reveal the crustal heterogeneity of the massif. Crustal growth curve demonstrates that the initial crust formed in the Mesoarchean, and shows a step-like pattern with three growth periods: 2.9-2.7, 2.1-1.9, and 1.7-0.5 Ga. This suggests that microcontinental massifs in the eastern CAOB have Precambrian basement, contradicting the hypothesis of significant crustal growth during the Phanerozoic. Phases of growth are constrained by multiple tectonic settings related to supercontinent development. Calculated reworked crustal proportions and the reworking curve indicate four reworking periods at 1.86-1.78 Ga, 860-720 Ma, 500-440 Ma, and 300-120 Ma, which limited the growth rate. These periods of reworking account for the crustal heterogeneity of the Erguna Massif.

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 εHfinitial 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.

Magmatic record of India-Asia collision.

New geochronological and geochemical data on magmatic activity from the India-Asia collision zone enables recognition of a distinct magmatic flare-up event that we ascribe to slab breakoff. This tie-point in the collisional record can be used to back-date to the time of initial impingement of the Indian continent with the Asian margin. Continental arc magmatism in southern Tibet during 80-40 Ma migrated from south to north and then back to south with significant mantle input at 70-43 Ma. A pronounced flare up in magmatic intensity (including ignimbrite and mafic rock) at ca. 52-51 Ma corresponds to a sudden decrease in the India-Asia convergence rate. Geological and geochemical data are consistent with mantle input controlled by slab rollback from ca. 70 Ma and slab breakoff at ca. 53 Ma. We propose that the slowdown of the Indian plate at ca. 51 Ma is largely the consequence of slab breakoff of the subducting Neo-Tethyan oceanic lithosphere, rather than the onset of the India-Asia collision as traditionally interpreted, implying that the initial India-Asia collision commenced earlier, likely at ca. 55 Ma.

A change in the geodynamics of continental growth 3 billion years ago.

Models for the growth of continental crust rely on knowing the balance between the generation of new crust and the reworking of old crust throughout Earth's history. The oxygen isotopic composition of zircons, for which uranium-lead and hafnium isotopic data provide age constraints, is a key archive of crustal reworking. We identified systematic variations in hafnium and oxygen isotopes in zircons of different ages that reveal the relative proportions of reworked crust and of new crust through time. Growth of continental crust appears to have been a continuous process, albeit at variable rates. A marked decrease in the rate of crustal growth at ~3 billion years ago may be linked to the onset of subduction-driven plate tectonics.