Mantle plumes and episodic crustal growth

Destabilization of deep oxidized mantle drove the Great Oxidation Event

The rise of Earth's atmospheric O₂ levels at ~2.4 Ga was driven by a shift between increasing sources and declining sinks of oxygen. Here, we compile recent evidence that the mantle shows a significant increase in oxidation state leading to the Great Oxidation Event (GOE), linked to sluggish upward mixing of a deep primordial oxidized layer. We simulate this scenario by implementing a new rheological model for this oxidized, bridgmanite-enriched viscous material and demonstrate slow mantle mixing in simulations of early Earth's mantle. The eventual homogenization of this layer may take ~2 Ga, in line with the timing of the observed mantle redox shift, and would result in the increase in upper mantle oxidation of >1 log(fO₂) unit. Such a shift would alter the redox state of volcanic degassing products to more oxidized species, removing a major sink of atmospheric O₂ and allowing oxygen levels to rise at ~2.4 Ga.

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.

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.

Building cratonic keels in Precambrian plate tectonics

The ancient cores of continents (cratons) are underlain by mantle keels-volumes of melt-depleted, mechanically resistant, buoyant and diamondiferous mantle up to 350 kilometres thick, which have remained isolated from the hotter and denser convecting mantle for more than two billion years. Mantle keels formed only in the Early Earth (approximately 1.5 to 3.5 billion years ago in the Precambrian eon); they have no modern analogues^1-4. Many keels show layering in terms of degree of melt depletion^5-7. The origin of such layered lithosphere remains unknown and may be indicative of a global tectonics mode (plate rather than plume tectonics) operating in the Early Earth. Here we investigate the possible origin of mantle keels using models of oceanic subduction followed by arc-continent collision at increased mantle temperatures (150-250 degrees Celsius higher than the present-day values). We demonstrate that after Archaean plate tectonics began, the hot, ductile, positively buoyant, melt-depleted sublithospheric mantle layer located under subducting oceanic plates was unable to subduct together with the slab. The moving slab left behind craton-scale emplacements of viscous protokeel beneath adjacent continental domains. Estimates of the thickness of this sublithospheric depleted mantle show that this mechanism was efficient at the time of the major statistical maxima of cratonic lithosphere ages. Subsequent conductive cooling of these protokeels would produce mantle keels with their low modern temperatures, which are suitable for diamond formation. Precambrian subduction of oceanic plates with highly depleted mantle is thus a prerequisite for the formation of thick layered lithosphere under the continents, which permitted their longevity and survival in subsequent plate tectonic processes.

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.

Episodic zircon age spectra mimic fluctuations in subduction

Decades of geochronological work have shown the temporal distribution of zircon ages to be episodic on billion-year timescales and seemingly coincident with the lifecycle of supercontinents, but the physical processes behind this episodicity remain contentious. The dominant, end-member models of fluctuating magmatic productivity versus selective preservation of zircon during times of continental assembly have important and very different implications for long-term, global-scale phenomena, including the history of crustal growth, the initiation and evolution of plate tectonics, and the tempo of mantle outgassing over billions of years. Consideration of this episodicity has largely focused on the Precambrian, but here we analyze a large collection of Phanerozoic zircon ages in the context of global, full-plate tectonic models that extend back to the mid-Paleozoic. We scrutinize two long-lived and relatively simple active margins, and show that along both, a relationship between the regional subduction flux and zircon age distribution is evident. In both cases, zircon age peaks correspond to intervals of high subduction flux with a ~10–30 Ma time lag (zircons trailing subduction), illuminating a possibly intrinsic delay in the subduction-related magmatic system. We also show that subduction fluxes provide a stronger correlation to zircon age distributions than subduction lengths do, implying that convergence rates play a significant role in regulating the volume of melting in subduction-related magmatic systems, and thus crustal growth.

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

Earth's punctuated tectonic evolution: cause and effect

Peaks in the Precambrian preserved crustal record are associated with major volcanic, tectonic and climatic events. These include addition of juvenile continental crust, voluminous high-temperature volcanism, massive mantle depletion, widespread orogeny and mineralization, large apparent polar wander velocity spikes, and subsequent palaeomagnetic intensity increases. These events impinge on the glaciation record, atmospheric and ocean chemistry, and on the rise of oxygen. Here we summarize and assess a number of geodynamic models that have been proposed to explain the observed episodicity in the Precambrian record. We find that episodic behaviour from nonlinear slab-driven models best explains the observed record. Examples of such slab-driven systems include mantle avalanches or episodic subduction. In these cases, rapid descent of slabs into the mantle drives fast plate motions and convergence at the surface. This is accompanied by large-scale upwellings of deep hot mantle, which contribute to voluminous volcanism. Further modelling will determine the relative importance of each mechanism, and reinforce the fundamental contribution of the mantle to the evolution of Earth's surface systems.

Temporal relations between mineral deposits and global tectonic cycles

Mineral deposits are heterogeneously distributed in both space and time, with variations reflecting tectonic setting, evolving environmental conditions, as in the atmosphere and hydrosphere, and secular changes in the Earth's thermal history. The distribution of deposit types whose settings are tied to plate margin processes (e.g. orogenic gold, volcanic-hosted massive sulphide, Mississippi valley type Pb–Zn deposits) correlates well with the supercontinent cycle, whereas deposits related to intra-cratonic settings and mantle-driven igneous events, such as Ni–Cu–PGE deposits, lack a clear association. The episodic distribution of deposits tied to the supercontinent cycle is accentuated by selective preservation and biasing of rock units and events during supercontinent assembly, a process that encases the deposit within the assembled supercontinent and isolates it from subsequent removal and recycling at plate margins.

Arc-continent collisions, sediment recycling and the maintenance of the continental crust

Subduction zones are both the source of most new continental crust and the locations where crustal material is returned to the upper mantle. Globally the total amount of continental crust and sediment subducted below forearcs currently lies close to 3.0 Armstrong Units (1 AU=1 km3a−1), of which 1.65 AU comprises subducted sediments and 1.33 AU tectonically eroded forearc crust, compared with an average ofc. 0.4 AU lost during subduction of passive margins during Cenozoic continental collision. Margins may retreat in a wholesale, steady-state mode, or in a slower way involving the trenchward erosion of the forearc coupled with landward underplating, such as seen in the central and northern Andean margins. Tephra records of magmatism evolution from Central America indicate pulses of recycling through the roots of the arc. While this arc is in a state of long-term mass loss this is achieved in a discontinuous fashion via periods of slow tectonic erosion and even sediment accretion interrupted by catastrophic erosion events, probably caused by seamount subduction. Crustal losses into subduction zones must be balanced by arc magmatism and we estimate global average melt production rates to be 96 and 64 km3Ma−1km−1in oceanic and continental arc, respectively. Critical to maintaining the volume of the continental crust is the accretion of oceanic arcs to continental passive margins. Mass balancing across the Taiwan collision zones suggests that almost 90% of the colliding Luzon Arc crust is accreted to the margin of Asia in that region. Rates of exhumation and sediment recycling indicate that the complete accretion process spans only 6–8 Ma. Subduction of sediment in both erosive and inefficient accretionary margins provides a mechanism for returning continental crust to the upper mantle. Sea level governs rates of continental erosion and thus sediment delivery to trenches, which in turn controls crustal thicknesses over 107–109years. Tectonically thickened crust is reduced to normal values (35–38 km) over time scales of 100–200 Ma.

Accretionary orogens through Earth history

Accretionary orogens form at intraoceanic and continental margin convergent plate boundaries. They include the supra-subduction zone forearc, magmatic arc and back-arc components. Accretionary orogens can be grouped into retreating and advancing types, based on their kinematic framework and resulting geological character. Retreating orogens (e.g. modern western Pacific) are undergoing long-term extension in response to the site of subduction of the lower plate retreating with respect to the overriding plate and are characterized by back-arc basins. Advancing orogens (e.g. Andes) develop in an environment in which the overriding plate is advancing towards the downgoing plate, resulting in the development of foreland fold and thrust belts and crustal thickening. Cratonization of accretionary orogens occurs during continuing plate convergence and requires transient coupling across the plate boundary with strain concentrated in zones of mechanical and thermal weakening such as the magmatic arc and back-arc region. Potential driving mechanisms for coupling include accretion of buoyant lithosphere (terrane accretion), flat-slab subduction, and rapid absolute upper plate motion overriding the downgoing plate. Accretionary orogens have been active throughout Earth history, extending back until at least 3.2 Ga, and potentially earlier, and provide an important constraint on the initiation of horizontal motion of lithospheric plates on Earth. They have been responsible for major growth of the continental lithosphere through the addition of juvenile magmatic products but are also major sites of consumption and reworking of continental crust through time, through sediment subduction and subduction erosion. It is probable that the rates of crustal growth and destruction are roughly equal, implying that net growth since the Archaean is effectively zero.

Helium isotopic evidence for episodic mantle melting and crustal growth

The timing of formation of the Earth's continental crust is the subject of a long-standing debate, with models ranging from early formation with little subsequent growth, to pulsed growth, to steadily increasing growth. But most models do agree that the continental crust was extracted from the mantle by partial melting. If so, such crustal extraction should have left a chemical fingerprint in the isotopic composition of the mantle. The subduction of oceanic crust and subsequent convective mixing, however, seems to have largely erased this record in most mantle isotopic systems (for example, strontium, neodymium and lead). In contrast, helium is not recycled into the mantle because it is volatile and degasses from erupted oceanic basalts. Therefore helium isotopes may potentially preserve a clearer record of mantle depletion than recycled isotopes. Here I show that the spectrum of 4He/3He ratios in ocean island basalts appears to preserve the mantle's depletion history, correlating closely with the ages of proposed continental growth pulses. The correlation independently predicts both the dominant 4He/3He peak found in modern mid-ocean-ridge basalts, as well as estimates of the initial 4He/3He ratio of the Earth. The correspondence between the ages of mantle depletion events and pulses of crustal production implies that the formation of the continental crust was indeed episodic and punctuated by large, potentially global, melting events. The proposed helium isotopic evolution model does not require a primitive, undegassed mantle reservoir, and therefore is consistent with whole mantle convection.

Evolution of the continental crust

The continental crust covers nearly a third of the Earth's surface. It is buoyant--being less dense than the crust under the surrounding oceans--and is compositionally evolved, dominating the Earth's budget for those elements that preferentially partition into silicate liquid during mantle melting. Models for the differentiation of the continental crust can provide insights into how and when it was formed, and can be used to show that the composition of the basaltic protolith to the continental crust is similar to that of the average lower crust. From the late Archaean to late Proterozoic eras (some 3-1 billion years ago), much of the continental crust appears to have been generated in pulses of relatively rapid growth. Reconciling the sedimentary and igneous records for crustal evolution indicates that it may take up to one billion years for new crust to dominate the sedimentary record. Combining models for the differentiation of the crust and the residence time of elements in the upper crust indicates that the average rate of crust formation is some 2-3 times higher than most previous estimates.

Episodic growth of the Gondwana supercontinent from hafnium and oxygen isotopes in zircon

It is thought that continental crust existed as early as 150 million years after planetary accretion, but assessing the rates and processes of subsequent crustal growth requires linking the apparently contradictory information from the igneous and sedimentary rock records. For example, the striking global peaks in juvenile igneous activity 2.7, 1.9 and 1.2 Gyr ago imply rapid crustal generation in response to the emplacement of mantle 'super-plumes', rather than by the continuous process of subduction. Yet uncertainties persist over whether these age peaks are artefacts of selective preservation, and over how to reconcile episodic crust formation with the smooth crustal evolution curves inferred from neodymium isotope variations of sedimentary rocks. Detrital zircons encapsulate a more representative record of igneous events than the exposed geology and their hafnium isotope ratios reflect the time since the source of the parental magmas separated from the mantle. These 'model' ages are only meaningful if the host magma lacked a mixed or sedimentary source component, but the latter can be diagnosed by oxygen isotopes, which are strongly fractionated by rock-hydrosphere interactions. Here we report the first study that integrates hafnium and oxygen isotopes, all measured in situ on the same, precisely dated detrital zircon grains. The data reveal that crust generation in part of Gondwana was limited to major pulses at 1.9 and 3.3 Gyr ago, and that the zircons crystallized during repeated reworking of crust formed at these times. The implication is that the mechanisms of crust formation differed from those of crustal differentiation in ancient orogenic belts.

Controls on the heterogeneous distribution of mineral deposits through time

Mineral deposits exhibit heterogeneous distributions, with each major deposit type showing distinctive, commonly unique, temporal patterns. These reflect a complex interplay between formational and preservational forces that, in turn, largely reflect changes in tectonic processes and environmental conditions in an evolving Earth. The major drivers were the supercontinent cycle and evolution from plume-dominated to modern-style plate tectonics in a cooling Earth. Consequent decrease in the growth rate of continental crust, and change from thick, buoyant sub-continental lithospheric mantle (SCLM) in the Precambrian to thinner, negatively buoyant SCLM in the Phanerozoic, led to progressive decoupling of formational and preservational processes through time. This affected the temporal patterns of deposit types including orogenic gold, porphyry and epithermal deposits, volcanic hosted massive sulphide (VHMS), palaeoplacer Au, iron oxide, copper gold (IOCG), platinum group elements (PGE), diamond and probably massive sulphide SEDEX deposits. Sedimentary mineral deposits mined for redox-sensitive metals show highly anomalous temporal patterns in which specific deposit types are restricted to particular times in Earth history. In particular, palaeoplacer uranium, banded iron formation (BIF) and BIF-associated manganese carbonates that formed in the early Precambrian do not reappear in younger basins. The most obvious driver is progressive oxidation of the atmosphere, with consequent long-term changes in the hydrosphere and biosphere, the latter influencing the temporal distribution and peak development of deposits such as Mississippi Valley types (MVT), hosted in biogenic sedimentary rocks.

Zoned mantle convection

We review the present state of our understanding of mantle convection with respect to geochemical and geophysical evidence and we suggest a model for mantle convection and its evolution over the Earth's history that can reconcile this evidence. Whole-mantle convection, even with material segregated within the D" region just above the core-mantle boundary, is incompatible with the budget of argon and helium and with the inventory of heat sources required by the thermal evolution of the Earth. We show that the deep-mantle composition in lithophilic incompatible elements is inconsistent with the storage of old plates of ordinary oceanic lithosphere, i.e. with the concept of a plate graveyard. Isotopic inventories indicate that the deep-mantle composition is not correctly accounted for by continental debris, primitive material or subducted slabs containing normal oceanic crust. Seismological observations have begun to hint at compositional heterogeneity in the bottom 1000 km or so of the mantle, but there is no compelling evidence in support of an interface between deep and shallow mantle at mid-depth. We suggest that in a system of thermochemical convection, lithospheric plates subduct to a depth that depends - in a complicated fashion - on their composition and thermal structure. The thermal structure of the sinking plates is primarily determined by the direction and rate of convergence, the age of the lithosphere at the trench, the sinking rate and the variation of these parameters over time (i.e. plate-tectonic history) and is not the same for all subduction systems. The sinking rate in the mantle is determined by a combination of thermal (negative) and compositional buoyancy and as regards the latter we consider in particular the effect of the loading of plates with basaltic plateaux produced by plume heads. Barren oceanic plates are relatively buoyant and may be recycled preferentially in the shallow mantle. Oceanic plateau-laden plates have a more pronounced negative buoyancy and can more easily founder to the very base of the mantle. Plateau segregation remains statistical and no sharp compositional interface is expected from the multiple fate of the plates. We show that the variable depth subduction of heavily laden plates can prevent full vertical mixing and preserve a vertical concentration gradient in the mantle. In addition, it can account for the preservation of scattered remnants of primitive material in the deep mantle and therefore for the Ar and (3)He observations in ocean-island basalts.

Late cretaceous polar wander of the pacific plate: evidence of a rapid true polar wander event

We reexamined the Late Cretaceous-early Tertiary apparent polar wander path for the Pacific plate using 27 paleomagnetic poles from seamounts dated by (40)Ar/(39)Ar geochronology. The path shows little motion from 120 to 90 million years ago (Ma), northward motion from 79 to 39 Ma, and two groups of poles separated by 16 to 21 degrees with indistinguishable mean ages of 84 +/- 2 Ma. The latter phenomenon may represent a rapid polar wander episode (3 to 10 degrees per million years) whose timing is not adequately resolved with existing data. Similar features in other polar wander paths imply that the event was a rapid shift of the spin axis relative to the mantle (true polar wander), which may have been related to global changes in plate motion, large igneous province eruptions, and a shift in magnetic field polarity state.

Lead and Helium Isotope Evidence from Oceanic Basalts for a Common Deep Source of Mantle Plumes

Linear arrays in lead isotope space for mid-ocean ridge basalts (MORBs) converge on a single end-member component that has intermediate lead, strontium, and neodymium isotope ratios compared with the total database for oceanic island basalts (OIBs) and MORBs. The MORB data are consistent with the presence of a common mantle source region for OIBs that is sampled by mantle plumes. 3He/4He ratios for MORBs show both positive and negative correlation with the 206Pb/204Pb ratios, depending on the MORB suite. These data suggest that the common mantle source is located in the transition zone region. This region contains recycled, oceanic crustal protoliths that incorporated some continental lead before their subduction during the past 300 to 2000 million years.

Mantle Convection with Plates and Mobile, Faulted Plate Margins

A finite-element formulation of faults has been incorporated into time-dependent models of mantle convection with realistic rheology, continents, and phase changes. Realistic tectonic plates naturally form with self-consistent coupling between plate and mantle dynamics. After the initiation of subduction, trenches rapidly roll back with subducted slabs temporarily laid out along the base of the transition zone. After the slabs have penetrated into the lower mantle, the velocity of trench migration decreases markedly. The inhibition of slab penetration into the lower mantle by the 670-kilometer phase change is greatly reduced in these models as compared to models without tectonic plates.