Explore the vast range of titles available.

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'.

The paleomagnetic record is an archive of Earth’s geophysical history, informing reconstructions of ancient plate motions and probing the core via the geodynamo. We report a robust 3.25-billion-year-old (Ga) paleomagnetic pole from the East Pilbara Craton, Western Australia. Together with previous results from the East Pilbara between 3.34 and 3.18 Ga, this pole enables the oldest reconstruction of time-resolved lithospheric motions, documenting 160 My of both latitudinal drift and rotation at rates of at least 0.55°/My. Motions of this style, rate, and duration are difficult to reconcile with true polar wander or stagnant-lid geodynamics, arguing strongly for mobile-lid geodynamics by 3.25 Ga. Additionally, this pole includes the oldest documented geomagnetic reversal, reflecting a stably dipolar, core-generated Archean dynamo.

Subduction processes play a pivotal role in facilitating material exchange between the crust and mantle, contributing to the growth of continents. However, the onset and evolution of subduction remain hotly debated. Here, we developed a high‐dimensional machine learning (ML) model to use multiple compositional data (e.g., Nb/La, Nb, Ti, Nb/U, Pb/Nd, and Nb/Th) to distinguish arc‐type from non‐arc basalts worldwide, then applied this well‐trained ML model to identify and delineate the secular occurrence of Archean arc‐type basalts. Our findings indicate that basaltic arc magmatism can be traced back to only a few early Archean cratons. A noticeable increase in the prevalence of arc‐type basalts during the Mesoarchean and Neoarchean suggests the synchronous formation of Archean subduction across different cratons. This observation hints at transitioning from the predominantly stagnant geodynamic regime to the early, more mobile tectonic regime. Plain Language Summary Deep beneath the Earth's surface, subduction moves material from the crust to the mantle, generating magma and igneous activity. This process is crucial for building the continental crust. Despite its importance, how the subduction process evolved over billions of years remains debated. To better understand when and where subduction began, we trained a machine learning model to distinguish between arc (subduction‐related) and non‐arc (not subduction‐related) basalts to recognize the chemical “fingerprints” of subduction in these rocks. By applying the same model to ancient basalts, we could trace basaltic arc‐like magmatism back to several Archean eon cratons, known as cratons, early in Earth's history (during the Archean Eon). The number of arc‐type basalts increased during the Mesoarchean and late Archean, suggesting that subduction expanded to different regions for about 700 million years. The transition suggests a significant shift in how the Earth's crust moved and changed, from a primarily steady state to a more mobile tectonics as observed today. Our study sheds new light on the ancient history of our planet and helps us understand how the Earth's crust has evolved over billions of years. Key Points Machine learning can accurately distinguish between arc and non‐arc basalts worldwide Subduction‐related basaltic magmatism may have produced on several cratons during the early Archean Subduction‐related processes may have started in cratons worldwide during the late Archean

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 (O₂-whiffs) and glaciation on Archean Earth.

The Archaean (4.0–2.5 Ga) continental crust is mainly composed of granitoids, whose geochemical characteristics are a function of their formation mechanisms and components, as well as physical conditions of their source. Therefore, revealing changes in Archaean geodynamic processes requires understanding of geochemical changes in Archaean granitoids. This paper compares key geochemical signatures in granitoid occurrences from the Eoarchaean to Neoarchaean Eras and aims to highlight changes or variations in their geochemical signatures. The study is performed by exploring and comparing geochemical and geochronological datasets of Archaean granitoids compiled from literature. The results show that two end-members of sodic TTGs (tonalite–trondhjemite–granodiorite) occur throughout the Archaean: low- and high-HREE (heavy rare earth elements) types. A profound change in granitoid geochemistry occurred between 3.0 and 2.5 Ga when multi-source high-K calc-alkaline granitoid batholiths emerged, possibly indicating the onset of modern-type plate tectonics.

The relative depletion of high field strength elements (HFSE), such as Nb, Ta and Ti, on normalised trace-element plots is a geochemical proxy routinely used to fingerprint magmatic processes linked to Phanerozoic subduction. This proxy has increasingly been applied to ultramafic-mafic units in Archaean cratons, but as these assemblages have commonly been affected by high-grade metamorphism and hydrothermal alteration/metasomatism, the likelihood of element mobility is high relative to Phanerozoic examples. To assess the validity of HFSE anomalies as a reliable proxy for Archaean subduction, we here investigate their origin in ultramafic rocks from the Ben Strome Complex, which is a 7 km2 ultramafic-mafic complex in the Lewisian Gneiss Complex of NW Scotland. Recently interpreted as a deformed layered intrusion, the Ben Strome Complex has been subject to multiple phases of high-grade metamorphism, including separate granulite- and amphibolite-facies deformation events. Additional to bulk-rock geochemistry, we present detailed petrography, and major- and trace-element mineral chemistry for 35 ultramafic samples, of which 15 display negative HFSE anomalies. Our data indicate that the magnitude of HFSE anomalies in the Ben Strome Complex are correlated with light rare earth-element (LREE) enrichment likely generated during interaction with H2O and CO2-rich hydrothermal fluids associated with amphibolitisation, rather than primary magmatic (subduction-related) processes. Consequently, we consider bulk-rock HFSE anomalies alone to be an unreliable proxy for Archaean subduction in Archaean terranes that have experienced multiple phases of high-grade metamorphism, with a comprehensive assessment of element mobility and petrography a minimum requirement prior to assigning geodynamic interpretations to bulk-rock geochemical data.

Understanding the pressure conditions of the formation of tonalite‐trondhjemite‐granodiorite (TTG) series rocks is important, as these rocks are key constituents of the Archean continental crust. However, the detailed pressure history of Archean TTG gneisses remains ambiguous due to uncertainties surrounding their sources and the effects of fractional crystallization. To address this problem, we developed a machine learning‐based tool to predict the melting pressures of Archean TTG gneisses from their major element compositions. This method provides high accuracy, supports batch processing, and incorporates the effects of fractional crystallization, mineral segregation, and other factors, making it suitable for application to global data sets. Our results indicate that the Earth's oldest rocks record two key transitions from low to medium pressure conditions, suggesting that the formation of the earliest continental crust may not have required deep subduction‐related geodynamic processes. High‐pressure rocks and melting conditions first emerged around 3.8 Ga (up to 1.5 GPa), becoming more extensive globally after 3.5 Ga (up to 1.9–2.0 GPa). Fractional crystallization, source, and other petrogenetic factors significantly influence the identification of high‐pressure TTGs, suggesting that true global high‐pressure TTGs may make up only ∼8% of the total rock record, much lower than earlier estimates. Our findings suggest that after 3.5 Ga, increased crustal rigidity likely facilitated widespread high‐pressure melting and plate tectonic activity across numerous terranes. This research provides new insights into Archean geodynamic processes and demonstrates the power of machine learning in advancing our understanding of ancient Earth systems.

On the contemporary Earth, distinct plate tectonic settings are characterized by differences in heat flow that are recorded in metamorphic rocks as differences in apparent thermal gradients. In this study we compile thermal gradients [defined as temperature/pressure (T/P) at the metamorphic peak] and ages of metamorphism (defined as the timing of the metamorphic peak) for 456 localities from the Eoarchean to Cenozoic Eras to test the null hypothesis that thermal gradients of metamorphism through time did not vary outside of the range expected for each of these distinct plate tectonic settings. Based on thermal gradients, metamorphic rocks are classified into three natural groups: high dT/dP [>775°C/GPa, mean ∼1110°C/GPa (n = 199) rates], intermediate dT/dP [775-375°C/GPa, mean ∼575°C/GPa (n = 127)], and low dT/dP [<375°C/GPa, mean ∼255°C/GPa (n = 130)] metamorphism. Plots of T, P, and T/P against age demonstrate the widespread occurrence of two contrasting types of metamorphism-high dT/dP and intermediate dT/dP-in the rock record by the Neoarchean, the widespread occurrence of low dT/dP metamorphism in the rock record by the end of the Neoproterozoic, and a maximum in the thermal gradients for high dT/dP metamorphism during the period 2.3 to 0.85 Ga. These observations falsify the null hypothesis and support the alternative hypothesis that changes in thermal gradients evident in the metamorphic rock record were related to changes in geodynamic regime. Based on the observed secular changes, we postulate that the Earth has evolved through three geodynamic cycles since the Mesoarchean and has just entered a fourth. Cycle I began with the widespread appearance of paired metamorphism in the rock record, which was coeval with the amalgamation of widely dispersed blocks of protocontinental lithosphere into supercratons, and was terminated by the progressive fragmentation of the supercratons into protocontinents during the Siderian-Rhyacian (2.5 to 2.05 Ga). Cycle II commenced with the progressive reamalgamation of these protocontinents into the supercontinent Columbia and extended until the breakup of the supercontinent Rodinia in the Tonian (1.0 to 0.72 Ga). Thermal gradients of high dT/dP metamorphism rose around 2.3 Ga leading to a thermal maximum in the mid-Mesoproterozoic, reflecting insulation of the mantle beneath the quasi-integral continental lithosphere of Columbia, prior to the geographical reorganization of Columbia into Rodinia. This cycle coincides with the age span of most anorogenic magmatism on Earth and a scarcity of passive margins in the geological record. Intriguingly, the volume of preserved continental crust of Mesoproterozoic age is low relative to the Paleoproterozoic and Neoproterozoic Eras. These features are consistent with a relatively stable association of continental lithosphere between the assembly of Columbia and the breakup of Rodinia. The transition to Cycle III during the Tonian is marked by a steep decline in the thermal gradients of high dT/dP metamorphism to their lowest value and the appearance of low dT/dP metamorphism in the rock record. Again, thermal gradients for high dT/dP metamorphism show a rise to a peak at the end of the Variscides during the formation of Pangea, before another steep decline associated with the breakup of Pangea and the start of a fourth cycle at ca. 0.175 Ga. Although the mechanism by which subduction started and plate boundaries evolved remains uncertain, based on the widespread record of paired metamorphism in the Neoarchean we posit that plate tectonics was established globally during the late Mesoarchean. During the Neoproterozoic there was a change to deep subduction and colder thermal gradients, features characteristic of the modern plate tectonic regime.

A comprehensive study based on U–Pb and Hf isotope analyses of zircons from gneisses has been conducted along the western part (Babina area) of the E–W-trending Bundelkhand Tectonic Zone in the central part of the Archaean Bundelkhand Craton. 207Pb–206Pb zircon ages and Hf isotopic data indicate the existence of a felsic crust at ~ 3.59 Ga, followed by a second tectonothermal event at ~ 3.44 Ga, leading to calc-alkaline magmatism and subsequent crustal growth. The study hence suggests that crust formation in the Bundelkhand Craton occurred in a similar time-frame to that recorded from the Singhbhum and Bastar cratons of the North Indian Shield.

Although the thermal evolution of the mantle before c. 3.0 Ga remains unclear, since c. 3.0 Ga secular cooling has dominated over heat production--this is time's arrow. By contrast, the thermal history of the crust, which is preserved in the record of metamorphism, is more complex. Heat to drive metamorphism is generated by radioactive decay and viscous dissipation, and is augmented by the influx of heat from the mantle. Notwithstanding that reliable data are sparse before the Neoarchean, we use a dataset of temperature (T), pressure (P) and thermobaric ratio (T/P at the metamorphic 'peak'), and age of metamorphism (t, the timing of the metamorphic 'peak') for rocks from 564 localities ranging in age from the Cenozoic to the Eoarchean eras to interrogate the crustal record of metamorphism as a proxy for the heat budget of the crust through time. On the basis of T/P, metamorphic rocks are classified into three natural groups: high T/P type (T/P >775°C/GPa, mean T/P ∼1105°C/GPa), including common and ultrahigh-temperature granulites, intermediate T/P type (T/P between 775 and 375°C/GPa, mean T/P ∼575°C/GPa), including high-pressure granulites and medium- and high-temperature eclogites, and low T/P type (T/P <375°C/GPa, mean T/P ∼255°C/GPa), including blueschists, low-temperature eclogites and ultrahigh-pressure metamorphic rocks. A monotonic increase in the P of intermediate T/P metamorphism from the Neoarchean to the Neoproterozoic reflects strengthening of the lithosphere during secular cooling of the mantle - this is also time's arrow. However, temporal variation in the P of intermediate T/P metamorphism and in the moving means of T and T/P of high T/P metamorphism, combined with the clustered age distribution, demonstrate the cyclicity of collisional orogenesis and cyclic variations in the heat budget of the crust superimposed on secular cooling since c. 3.0 Ga -this is time's cycle. A first cycle began with the widespread appearance/survival of intermediate T/P and high T/P metamorphism in the Neoarchean rock record coeval with amalgamation of dispersed blocks of lithosphere to form protocontinents. This cycle was terminated by the fragmentation of the protocontinents into cratons in the early Paleoproterozoic, which signalled the start of a new cycle. The second cycle continued with the progressive amalgamation of the cratons into the supercontinent Columbia and extended until the breakup of the supercontinent Rodinia in the Neoproterozoic. This cycle represented a period of relative tectonic and environmental stability, and perhaps reduced subduction during at least part of the cycle. During most of the Proterozoic the moving means for both T and T/P of high T/P metamorphism exceeded the arithmetic means, reflecting insulation of the mantle beneath the quasi-integrated lithosphere of Columbia and, after a limited reorganisation, Rodinia. The third cycle began with the steep decline in thermobaric ratios of high T/P metamorphism to their lowest value, synchronous with the breakup of Rodinia and the formation of Pannotia, and the widespread appearance/preservation of low T/P metamorphism in the rock record. The thermobaric ratios for high T/P metamorphism rise to another peak associated with the Pan-African event, again reflecting insulation of the mantle. The subsequent steep decline in thermobaric ratios of high T/P metamorphism associated with the breakup of Pangea at c. 0.175 Ga may indicate the start of a fourth cycle. The limited occurrence of high and intermediate T/P metamorphism before the Neoarchean suggests either that suitable tectonic environments to generate these types of metamorphism were not widely available before then or that the rate of survival was low. We interpret the first cycle to record stabilisation of subduction and the emergence of a network of plate boundaries in a plate tectonics regime once the balance between heat production and heat loss changed in favour of secular cooling, possibly as early as c. 3.0 Ga in some areas. This is inferred to have been a globally linked system by the early Paleoproterozoic, but whether it remained continuous to the present is unclear. The second cycle was characterised by stability from the formation of Columbia to the breakup of Rodinia, generating higher than average T and T/P of high T/P metamorphism. The third cycle reflects colder collisional orogenesis and deep subduction of the continental crust, features that are characteristic of modern plate tectonics, which became possible once the average temperature of the asthenospheric mantle had declined to <100°C warmer than the present day after c. 1.0 Ga.