Evolution and diversity of subduction zones controlled by slab width

Onset of double subduction controls plate motion reorganisation

Face-to-face double subduction systems, in which two oceanic plates subduct toward each other, are essential elements of plate tectonics. Two subduction zones in such systems are typically uneven in age and their spatially and temporally variable dynamics remain enigmatic. Here, with 2D numerical modelling, we demonstrate that the onset of the younger subduction zone strongly changes the dynamics of the older subduction zone. The waxing younger subduction may gradually absorb plate convergence from the older one, resulting in older subduction waning featured by the dramatic decrease in subduction rate and trench retreat. The dynamical transformation of subduction predominance alters the intraplate stress and mantle flow, regulating the relative motion among the three different plates. The process of waxing and waning of subduction zones controls plate motion reorganisation, providing a reference to interpret the past, present, and future evolution of several key double subduction regions found on the modern Earth.

Vertical tearing of subducting plates controlled by geometry and rheology of oceanic plates

Lateral non-uniform subduction is impacted by continuous plate segmentation owing to vertical tearing of the subducting plate. However, the dynamics and physical controls of vertical tearing remain controversial. Here, we employed 3D numerical models to investigate the effects of trench geometry (offset by a transform boundary) and plate rheology (plate age and the magnitude of brittle/plastic strain weakening) on the evolution of shear stress-controlled vertical tearing within a homogenous subducting oceanic plate. Numerical results suggest that the trench offset geometry could result in self-sustained vertical tearing as a narrow shear zone within the intact subducting oceanic plate, and that this process of tearing could operate throughout the entire subduction process. Further, the critical trench offset length for the maturation of vertical tearing is impacted by plate rheology. Comparison between numerical modelling results and natural observations suggests that vertical tearing attributed to trench offset geometry is broadly developed in modern subduction and collision systems worldwide.

Kimberlite eruptions driven by slab flux and subduction angle

Kimberlites are sourced from thermochemical upwellings which can transport diamonds to the surface of the crust. The majority of kimberlites preserved at the Earth's surface erupted between 250 and 50 million years ago, and have been attributed to changes in plate velocity or mantle plumes. However, these mechanisms fail to explain the presence of strong subduction signatures observed in some Cretaceous kimberlites. This raises the question whether there is a subduction process that unifies our understanding of the timing of kimberlite eruptions. We develop a novel formulation for calculating subduction angle based on trench migration, convergence rate, slab thickness and density to connect the influx of slab material into the mantle with the timing of kimberlite eruptions. We find that subduction angles combined with peaks in slab flux predict pulses of kimberlite eruptions. High rates of subducting slab material trigger mantle return flow that stimulates fertile reservoirs in the mantle. These convective instabilities transport slab-influenced melt to the surface at a distance inbound from the trench corresponding to the subduction angle. Our deep-time slab dip formulation has numerous potential applications including modelling the deep carbon and water cycles, and an improved understanding of subduction-related mineral deposits.

3D geodynamic-geomorphologic modelling of deformation and exhumation at curved plate boundaries: Implications for the southern Alaskan plate corner

Plate corners with extreme exhumation rates are important because they offer a perspective for understanding the interactions between tectonics and surface processes. The southern Alaskan margin with its curved convergent plate boundary and associated zones of localized uplift is a prime location to study active orogeny. Here, we present the results of fully-coupled thermo-mechanical (geodynamic) and geomorphologic numerical modelling, the design of which captures the key features of the studied area: subduction of oceanic lithosphere (Pacific plate) is adjacent to a pronounced asymmetric indenter dipping at a shallow angle (Yakutat microplate), which in turn is bounded to the east by a dextral strike-slip shear zone (Fairweather fault). The resulting first-order deformation/rock uplift patterns show strong similarities with observations. In particular, relatively young thermochronological ages are reproduced along the plate-bounding (Fairweather) transform fault and in the area of its transition to convergence (the St. Elias syntaxis). The focused exhumation of the Chugach Core also finds its equivalent in model predicted zones of high rock uplift rates in an isolated region above the indenter. From these results, we suggest that the general exhumation patterns observed in southern Alaska are controlled by mutually reinforcing effects of tectonic deformation and surface erosion processes.

Topography of mountain belts controlled by rheology and surface processes

It is widely recognized that collisional mountain belt topography is generated by crustal thickening and lowered by river bedrock erosion, linking climate and tectonics^1-4. However, whether surface processes or lithospheric strength control mountain belt height, shape and longevity remains uncertain. Additionally, how to reconcile high erosion rates in some active orogens with long-term survival of mountain belts for hundreds of millions of years remains enigmatic. Here we investigate mountain belt growth and decay using a new coupled surface process^5,6 and mantle-scale tectonic model⁷. End-member models and the new non-dimensional Beaumont number, Bm, quantify how surface processes and tectonics control the topographic evolution of mountain belts, and enable the definition of three end-member types of growing orogens: type 1, non-steady state, strength controlled (Bm > 0.5); type 2, flux steady state⁸, strength controlled (Bm ≈ 0.4-0.5); and type 3, flux steady state, erosion controlled (Bm < 0.4). Our results indicate that tectonics dominate in Himalaya-Tibet and the Central Andes (both type 1), efficient surface processes balance high convergence rates in Taiwan (probably type 2) and surface processes dominate in the Southern Alps of New Zealand (type 3). Orogenic decay is determined by erosional efficiency and can be subdivided into two phases with variable isostatic rebound characteristics and associated timescales. The results presented here provide a unified framework explaining how surface processes and lithospheric strength control the height, shape, and longevity of mountain belts.

Episodic back-arc spreading centre jumps controlled by transform fault to overriding plate strength ratio

Spreading centre jumps are a common feature of oceanic back-arc basins. Jumps are conventionally suggested to be triggered by plate velocity changes, pre-existing weaknesses, or punctuated events such as the opening of slab windows. Here, we present 3D numerical models of back-arc spreading centre jumps evolving naturally in a homogeneous subduction system surrounded by continents without a trigger event. Spreading centres jump towards their subduction zone if the distance from trench to spreading centre becomes too long. In particular, jumps to a new spreading centre occur when the resistance on the boundary transform faults enabling relative motion of back-arc and neighbouring plates is larger than the resistance to break the overriding plate closer to trench. Time and distance of spreading centres jumps are, thus, controlled by the ratio between the transform fault and overriding plate strengths. Despite being less complex than natural systems, our models explain why narrow subducting plates (e.g. Calabrian slab), have more frequent and closely-spaced spreading jumps than wider subduction zones (e.g. Scotia). It also explains why wide back-arc basins undergo no spreading centre jumps in their life cycle.

Back-arc spreading centre jumps have been suggested to be controlled by a number of different drivers. Here, the authors, using 3D numerical models, show that transform faults can trigger back-arc spreading centre jumps, without the need of any ad hoc factors.

Southward expanding plate coupling due to variation in sediment subduction as a cause of Andean growth

Growth of the Andes has been attributed to Cenozoic subduction. Although climatic and tectonic processes have been proposed to be first-order mechanisms, their interaction and respective contributions remain largely unclear. Here, we apply three-dimensional, fully-dynamic subduction models to investigate the effect of trench-axial sediment transport and subduction on Andean growth, a mechanism that involves both climatic and tectonic processes. We find that the thickness of trench-fill sediments, a proxy of plate coupling (with less sediments causing stronger coupling), exerts an important influence on the pattern of crustal shortening along the Andes. The southward migrating Juan Fernandez Ridge acts as a barrier to the northward flowing trench sediments, thus expanding the zone of plate coupling southward through time. Consequently, the predicted history of Andean shortening is consistent with observations. Southward expanding crustal shortening matches the kinematic history of inferred compression. These results demonstrate the importance of climate-tectonic interaction on mountain building.

Basic Role of Extrusion Processes in the Late Cenozoic Evolution of the Western and Central Mediterranean Belts

Tectonic activity in the Mediterranean area (involving migrations of old orogenic belts, formation of basins and building of orogenic systems) has been determined by the convergence of the confining plates (Nubia, Arabia and Eurasia). Such convergence has been mainly accommodated by the consumption of oceanic and thinned continental domains, triggered by the lateral escapes of orogenic wedges. Here, we argue that the implications of the above basic concepts can allow plausible explanations for the very complex time-space distribution of tectonic processes in the study area, with particular regard to the development of Trench-Arc-Back Arc systems. In the late Oligocene and lower–middle Miocene, the consumption of the eastern Alpine Tethys oceanic domain was caused by the eastward to SE ward migration/bending of the Alpine–Iberian belt, driven by the Nubia–Eurasia convergence. The crustal stretching that developed in the wake of that migrating Arc led to formation of the Balearic basin, whereas accretionary activity along the trench zone formed the Apennine belt. Since the collision of the Anatolian–Aegean–Pelagonian system (extruding westward in response to the indentation of the Arabian promontory) with the Nubia-Adriatic continental domain, around the late Miocene–early Pliocene, the tectonic setting in the central Mediterranean area underwent a major reorganization, aimed at activating a less resisted shortening pattern, which led to the consumption of the remnant oceanic and thinned continental domains in the central Mediterranean area.

Subduction history of the Caribbean from upper-mantle seismic imaging and plate reconstruction

The margins of the Caribbean and associated hazards and resources have been shaped by a poorly understood history of subduction. Using new data, we improve teleseismic P-wave imaging of the eastern Caribbean upper mantle and compare identified subducted-plate fragments with trench locations predicted from plate reconstruction. This shows that material at 700-1200 km depth below South America derives from 90-115 Myr old westward subduction, initiated prior to Caribbean Large-Igneous-Province volcanism. At shallower depths, an accumulation of subducted material is attributed to Great Arc of the Caribbean subduction as it evolved over the past 70 Ma. We interpret gaps in these subducted-plate anomalies as: a plate window and tear along the subducted Proto-Caribbean ridge; tearing along subducted fracture zones, and subduction of a volatile-rich boundary between Proto-Caribbean and Atlantic domains. Phases of back-arc spreading and arc jumps correlate with changes in age, and hence buoyancy, of the subducting plate.

Seismic evidence for subduction-induced mantle flows underneath Middle America

Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.

The motions of subducted slabs are expected to drive mantle flow around slab edges, however, evidence of deep mantle flow has so far remained elusive. Here, the authors present a Full Waveform Inversion 3-D anisotropy model which allows them to infer deep subduction-induced mantle flows underneath the Mid-Americas and the Caribbean.

Oceanic crust recycling controlled by weakening at slab edges

Retreating subduction zones such as the Lesser Antilles, Gibraltar and Scotia have been migrating towards the Atlantic Ocean by cutting their way through the oceanic crust. This spontaneously retreating subduction is enabled by the development of faults at the edges of the slab, but the physical mechanisms controlling fault propagation and direction remain unknown. Here, using 3D numerical subduction models we show that oceanic lithosphere recycling is mainly controlled by the intensity of strain-induced weakening of fractures forming at the edges of the slab. Intense strain-induced weakening causes predominantly brittle fault propagation and slab narrowing until detachment. Without weakening, preponderantly ductile slab edge propagation occurs, which causes slab widening. This rheological control is not affected by the proximity of non-weakened passive continental margins. Natural examples suggest that slab edges follow convergent paths that could be controlled by fractures weakening due to deep water penetration into the oceanic lithosphere.

Influence of the Thickness of the Overriding Plate on Convergence Zone Dynamics

The important role played by the upper plate in convergence zones dynamics has long been underestimated but is now more and more emphasized. However, the influence of its thickness and/or strength on orogenic systems evolution remains largely unknown. Here we present results from 3D thermo-mechanical numerical simulations of convergence zones (including oceanic subduction followed by continental subduction/collision), in which we vary the rheological profile of the overriding plate (OP). For this, we systematically modify the crustal thickness of the overriding lithosphere and the temperature at the Moho to obtain a thermal thickness of the overriding lithosphere ranging from 80 to 180 km. While all models share a common global evolution (i.e., slab sinking, interaction between slab and the 660 km discontinuity, continental subduction/collision, and slab breakoff), they also highlight first-order differences arising from the variations in the OP strength (thermal thickness). With a thin/weak OP, slab rollback is favored, the slab dip is low, the mantle flow above the slab is vigorous, and the trench migrates at a high rate compared to a thick/strong OP. In addition, slab breakoff and back-arc basin formation events occur significantly earlier than in models involving a thick OP. Our models therefore highlight the major role played by the thickness/strength of the OP on convergence zone dynamics and illustrate its influence in a quantitative way.

Surface and mantle records reveal an ancient slab tear beneath Gondwana

Vertical slab-tearing has been widely reported in modern convergent settings profoundly influencing subduction and mantle dynamics. However, evaluating a similar impact in ancient convergent settings, where oceanic plates have been subducted and the geological record is limited, remains challenging. In this study, we correlate the lower mantle structure, which retained the past subduction configuration, with the upper-plate geological record to show a deep slab rupture interpreted as a large-scale tearing event in the early Mesozoic beneath southwestern Gondwana. For this purpose, we integrated geochronological and geological datasets with P-wave global seismic tomography and plate-kinematic reconstructions. The development of a Late Triassic-Early Jurassic slab-tearing episode supports (i) a slab gap at lower mantle depths, (ii) a contrasting spatiotemporal magmatic evolution, (iii) a lull in arc activity, and (iv) intraplate extension and magmatism in the Neuquén and Colorado basins. This finding not only has implications for identifying past examples of a fundamental process that shapes subduction zones, but also illustrates an additional mechanism to trigger slab-tearing in which plate rupture is caused by opposite rotation of slab segments.

Early exhumation of the Frontal Cordillera (Southern Central Andes) and implications for Andean mountain-building at ~33.5°S

The Andes are the modern active example of a Cordilleran-type orogen, with mountain-building and crustal thickening within the upper plate of a subduction zone. Despite numerous studies of this emblematic mountain range, several primary traits of this orogeny remain unresolved or poorly documented. The onset of uplift and deformation of the Frontal Cordillera basement culmination of the Southern Central Andes is such an example, even though this structural unit appears as a first-order topographic and geological feature. To solve for this, new (U-Th)/He ages on apatite and zircon from granitoids of the Frontal Cordillera at ~33.5°S are provided here. These data, interpreted as an age-elevation thermochronological profile, imply continuous exhumation initiating well before ~12–14 Ma, and at most by ~22 Ma when considering the youngest zircon grain from the lowermost sample. The age of exhumation onset is then refined to ~20 Ma by combining these results with data on sedimentary provenance from the nearby basins. Such continuous exhumation since ~20 Ma needs to have been sustained by tectonic uplift on an underlying crustal-scale thrust ramp. Such early exhumation and associated uplift of the Frontal Cordillera invalidate the classically proposed east-vergent models of the Andes at this latitude. Additionally, they provide further support to recent views on Andean mountain-building proposing that the Andes at ~33.5°S grew firstly over west-vergent basement structures.

Slab narrowing in the Central Mediterranean: the Calabro-Ionian subduction zone as imaged by high resolution seismic tomography

A detailed 3D image of the Calabro-Ionian subduction system in the central Mediterranean was obtained by means of a seismic tomography, exploiting a large dataset of local earthquakes and computing algorithms able to build a dense grid of measure nodes. Results show that the slab is continuous below the southern sector of the Calabro-Peloritan Arc, but the deformation processes developing at its edges are leading to its progressive narrowing, influencing tectonics and magmatism at the surface, and with possible stress concentration in the tip zones. In the southwest, the deformation occurring at a free slab edge lead to propagation of a vertical lithospheric tear in the overriding plate, which extends along a NW-SE fault system (Aeolian-Tindari-Letojanni) up to about 30 km into the Ionian Sea; further southeast, the lithosphere appears only flexed and not broken yet. In the northeast, the slab seems to break progressively, parallel to the trench. Finally, northwest of Mt. Etna, the tomography highlights low V_P that can be related to an upwelling of deep mantle material likely flowing laterally through a window opened by the complete slab detachment.

Lithospheric folding by flexural slip in subduction zones as source for reverse fault intraslab earthquakes

Subduction requires the permanent generation of a bend fold in the subducting slab which mechanics is not well understood. Lithospheric bending of subducting slabs was traditionally considered to be accommodated by orthogonal flexure, generating extensional outer rise earthquakes responsible of the external arc elongation during folding. Here we explore the possibility of lithospheric flexure being accommodated through simple shear deformation parallel to the slab (folding by flexural slip) and evaluate this process as source of earthquakes. The seismicity predicted by flexural slip dominated slab bending explains a significant amount of intermediate earthquakes observed in subduction zones with different degrees of coupling. This mechanism predicts the generation of intraslab thrust earthquakes with fault planes subparallel to the slab top. Being the orientations of the fault planes the same for the interface thrust earthquakes and the flexural-slip intraslab earthquakes, the amount of seismic moment liberated by the interface could be significantly lower than considered before. This proposed seismic source should be taken into account in models and hazard studies of subduction zones. Determining the seismic generating processes in subduction zones and their characteristics is a fundamental issue for the correct assessment of the associated seismic and tsunami risk.

Andean mountain building and magmatic arc migration driven by subduction-induced whole mantle flow

Subduction along the western margin of South America has been active since the Jurassic, but Andean orogeny started in the middle Cretaceous and was preceded by backarc extension in the Jurassic-Early Cretaceous. The timing and sequence of these events has remained unexplained. Here I present a four-dimensional buoyancy-driven whole-mantle subduction model implying that the ~200 Myr geological evolution can be attributed to sinking of a wide slab into a layered mantle, where upper-mantle wide-slab subduction causes backarc extension, while whole-mantle (upper+lower) wide-slab subduction drives Andean orogeny. The model reproduces the maximum shortening and crustal thickness observed in the Central Andes and their progressive northward and southward decrease. The subduction evolution coincides with a 29° decrease in slab dip angle, explaining ~200 km of Jurassic-present eastward migration of the Central Andean magmatic arc. Such arc migration negates proposed long-term subduction erosion and continental destruction, but is consistent with long-term crustal growth.

The Andean orogeny commenced in the Cretaceous, but was preceded by backarc extension starting in the Jurassic. Here, the author presents a 4D geodynamic subduction model reproducing the evolution of overriding plate deformation along with a progressive decrease in slab dip, explaining Jurassic-present arc migration.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

South-American plate advance and forced Andean trench retreat as drivers for transient flat subduction episodes

At two trench segments below the Andes, the Nazca Plate is subducting sub-horizontally over ∼200–300 km, thought to result from a combination of buoyant oceanic-plateau subduction and hydrodynamic mantle-wedge suction. Whether the actual conditions for both processes to work in concert existed is uncertain. Here we infer from a tectonic reconstruction of the Andes constructed in a mantle reference frame that the Nazca slab has retreated at ∼2 cm per year since ∼50 Ma. In the flat slab portions, no rollback has occurred since their formation at ∼12 Ma, generating ‘horse-shoe' slab geometries. We propose that, in concert with other drivers, an overpressured sub-slab mantle supporting the weight of the slab in an advancing upper plate-motion setting can locally impede rollback and maintain flat slabs until slab tearing releases the overpressure. Tear subduction re-establishes a continuous slab and allows the process to recur, providing a mechanism for the transient character of flat slabs.

How flat slabs at subduction zones are created remains unclear. Here, the authors show that the Nazca slab has retreated at ∼2 cm per year since 50 Ma but no rollback has occurred in the last ∼12 Myr in the flat slab, implying that an overpressured sub-slab mantle can impede rollback and maintain a flat slab.

Upper- and mid-mantle interaction between the Samoan plume and the Tonga-Kermadec slabs

Mantle plumes are thought to play a key role in transferring heat from the core-mantle boundary to the lithosphere, where it can significantly influence plate tectonics. On impinging on the lithosphere at spreading ridges or in intra-plate settings, mantle plumes may generate hotspots, large igneous provinces and hence considerable dynamic topography. However, the active role of mantle plumes on subducting slabs remains poorly understood. Here we show that the stagnation at 660 km and fastest trench retreat of the Tonga slab in Southwestern Pacific are consistent with an interaction with the Samoan plume and the Hikurangi plateau. Our findings are based on comparisons between 3D anisotropic tomography images and 3D petrological-thermo-mechanical models, which self-consistently explain several unique features of the Fiji-Tonga region. We identify four possible slip systems of bridgmanite in the lower mantle that reconcile the observed seismic anisotropy beneath the Tonga slab (V(SH)>V(SV)) with thermo-mechanical calculations.

Plate tectonic controls on atmospheric CO2 levels since the Triassic

Climate trends on timescales of 10s to 100s of millions of years are controlled by changes in solar luminosity, continent distribution, and atmosphere composition. Plate tectonics affect geography, but also atmosphere composition through volcanic degassing of CO2 at subduction zones and midocean ridges. So far, such degassing estimates were based on reconstructions of ocean floor production for the last 150 My and indirectly, through sea level inversion before 150 My. Here we quantitatively estimate CO2 degassing by reconstructing lithosphere subduction evolution, using recent advances in combining global plate reconstructions and present-day structure of the mantle. First, we estimate that since the Triassic (250-200 My) until the present, the total paleosubduction-zone length reached up to ∼200% of the present-day value. Comparing our subduction-zone lengths with previously reconstructed ocean-crust production rates over the past 140 My suggests average global subduction rates have been constant, ∼6 cm/y: Higher ocean-crust production is associated with longer total subduction length. We compute a strontium isotope record based on subduction-zone length, which agrees well with geological records supporting the validity of our approach: The total subduction-zone length is proportional to the summed arc and ridge volcanic CO2 production and thereby to global volcanic degassing at plate boundaries. We therefore use our degassing curve as input for the GEOCARBSULF model to estimate atmospheric CO2 levels since the Triassic. Our calculated CO2 levels for the mid Mesozoic differ from previous modeling results and are more consistent with available proxy data.

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.