How do subducted slabs accumulate at different mantle depths?

Seismic tomography provides no evidence that slabs of oceanic lithosphere descend intact through the whole mantle to the core-mantle boundary. It might once have happened when they were capped by abundant high-density rocks, such as Precambrian banded-iron formations. A great many actively descending slabs have been shown to cease sinking, slide sideways and accumulate at depths around 660 and 1000 km. Until recently these discontinuities were been generally ascribed to transitions in the structure of the dominant mafic mineral olivine (Mg2SiO4) in mantle peridotite induced by increasing pressure and temperature. The resulting increases in mantle density supposedly form barriers to further slab descent. Pressure-induced mineral transitions in the slabs themselves that increase their density, such as pyroxene to garnet, may somehow be inhibited thereby leading to stagnation in slab descent. That may be true for the 660 km discontinuity, but for stagnation at 1000 km deep no such density-changing mineral transitions have shown up in high-P high-T mineralogical experiments. Some other process must therefore be responsible for slab descent to that depth. Recent work by geoscientists at several universities in China gives insights into what may be going on (Li, J., Li, K., Li, J. et al 2026. Dual slab stagnation depths controlled by grain-size-induced sporadic low-viscosity zones at around 1000 km depth. Nature Communications DOI: 10.1038/s41467-026-69987-9).

Four different modes of subduction at island arcs. Credit: Li et al. Fig 6

Jing Li and colleagues have focussed on the possibility that changes in the bulk viscosity of the mantle may play an important role. Their approach is twofold: experimental mineral physics and geodynamic modelling. Results suggest that recrystallization in the mantle when deeply penetrating slabs pass through it may patchily reduce the mantle’s grain size and thus its viscosity; the more so with larger volumes of subducted slab material. In turn, the resulting physical heterogeneity probably disrupts the steady downward passage of the slabs; fine-grained, less viscous zones ‘lubricating’ slab penetration, unchanged zones hindering it. The authors link such hypothetical micro-structural processes to modes of subduction that are currently active. They consider four modes of active subduction beneath island arcs with either a slow or a fast rate of trench retreat (see Figure). A slowly retreating trench system combined with low-viscosity patches at depth (Mode 1) results in penetration below 660 km and slab stagnation at 1000 km. Slow trench retreat with a homogenous lower mantle (Mode 2) gives rise to penetration and buckling of the descending slab between 660 and 1000 km. Fast trench retreat with a deeper low-viscosity zone (Mode 3), or with a homogeneous lower mantle (Mode 4) both result in slab stagnation at 660 km.

The models developed by Jing Li et al convincingly simulate various results of seismic tomography beneath island arcs. Interestingly, they suggest that the eventual assimilation of older slab materials into the deeper mantle (‘fossil’ slabs) may play a major role in mineral comminution and reduced mantle strength. That may leave behind low viscosity zones that later subduction may exploit. In fact, there are signs of possible fossil slabs in seismic tomograms more than 1000 km below the present Pacific Ocean floor in the form of zones of high P-wave velocity.

This work shows that plate tectonics is far from ‘done-and-dusted’, the mantle being far from uniform in its properties. Li et al’s results potentially open up new insights into whole-mantle convection, in which older tectonic events influence plate motions that are currently operating and the triggering of plumes rising from the deepest mantle. It also hints that such complex physical mixing of subducted material into the mantle may have resulted in the geochemical heterogeneities that increasingly emerge from analysis of magmas with ultimate origins in the mantle.

See also:Grain Size Creates Dual Slab Stagnation Zones at 1000 km. Scienmag 3 March 2026

Vanished continents of the Hadean Eon: the zircon key

Over the last few decades improved analytical techniques have made it possible to analyse tiny mineral grains for a variety of trace elements and several isotopes. Zircons obtained directly from crushed granitic igneous rocks vary in chemistry according to the magmatic processes that generated them and their tectonic context. Elevated ratios between uranium and niobium (U/Nb) and scandium and ytterbium (Sc/Yb) are characteristic of zircons in intermediate granites. These contain 52 to 63 % SiO2 – between mafic and felsic magmas – which formed by melting of hydrated mafic crust in settings akin to modern continental arcs; i.e. in subduction zones. But such partial melting can also take place where the base of continental crust delaminates and ‘drips’ into the mantle. That process is part of what is known as stagnant lid tectonics, believed by many to have been important in the Palaeoarchaean and Hadean. Such a process would have involved nearly anhydrous conditions and thus different geochemical partitioning of elements in the magmas and minerals that crystallised from them. Exposures of crystalline continental crust become increasingly rare further back in geological time, and there are none older than 4.0 Ga – i.e. of Hadean age – with a granitic component. Consequently studying the generation of continental crust in the Hadean and the early Archaean is almost entirely dependent on ancient zircons that found their way into much younger sedimentary rocks. The most famous of these occur as detrital grains in the 3.6 Ga Jack Hills conglomerate of Western Australia. Others have been extracted from similar ~3.3 Ga sedimentary rocks in the Barberton Greenstone Belt of South Africa and Eswatini.

Cartoon of possible Hadean stagnant lid tectonics, dominated by mantle plumes. (Credit: Bédard, J.H. 2018, Fig 3B, DOI: 10.1016/j.gsf.2017.01.005)

John Valley of the University of Wisconsin-Madison, USA, and co-workers from the US, Germany, Australia and France have worked on a large number of zircons newly extracted from Jack Hills. They have radiometrically dated them, and analysed Nb, Sc, U and Yb trace elements and hafnium (Hf) and oxygen isotopes Together with data from earlier studies, including Barberton zircons, they have teased out some remarkable insights into  ‘continent-forming’ magmatism as far back in time as 4.4 billion years ago (Valley, J.W. and 11 others 2026. Contemporaneous mobile- and stagnant-lid tectonics on the Hadean Earth. Nature, Open access; DOI: 10.1038/s41586-025-10066-2). More than 70% of the >4.0 Ga Jack Hills zircons have elevated U/Nb and Sc/Yb ratios, which suggest that they formed in a setting akin to continental-arc subduction (CAS) zones, to produce now-vanished Hadean continental crust. The remainder seem to represent processes at mid-ocean ridge (MOR) and oceanic island (OI) settings. In contrast, the bulk of Barberton zircons of Hadean age show OI affinities, with only around 22% showing Nb–Sc–U–Yb signatures of probable CAS origins. From about 4.4 to 3.8 Ga two distinct forms of continental crust generation seem to have operated on Earth. In the erosional source region for the Barberton zircons their host granites seem to have formed during the Hadean and Eoarchaean by remelting of foundered lower crust, i.e. probably in a stagnant-lid-like tectonic setting. But at around 3.6 Ga they ‘flip’ to a subduction-like setting. The zircons yielded by Jack Hills conglomerates suggest substantially different conditions: alternating CAS and OI settings during the Hadean and a fall-off in crust generation during the Eoarchaean (4.0 to 3.8 Ga).

Plots of Sc/Yb and U/Nb against ages of zircons (vertical scale logarithmic). Black points are from Jack Hills, red from Barberton. The yellow field represents zircons formed in subduction zones; green suggests stagnant lid tectonics; grey the overlap between the two settings. Credit: Valley et al. Fig 3 a and b.

The mixed Hadean zircon signatures from Jack Hills possibly indicate that they were derived by erosion and transport from several distinct terranes that had been generated by two different processes: some kind of upper crustal recycling and stagnant lid tectonics. Meanwhile, that part of the Hadean Earth represented by the Barberton zircons may have been a long-lived regime of stagnant lid tectonics, replaced by dominant subduction at the end of the Eoarchaean.  Yet the data suggest that into the Palaeoarchaean (3.6 to 3.2 Ga) and perhaps later, lid tectonics continued to operate somewhere, but at no time after 4.4 Ga was the Earth entirely subject to lid tectonics. Likewise, the authors insist that subduction was not of the plate-tectonic style, referring to some form of recycling of hydrated upper crustal mafic and ultramafic rocks into the mantle to undergo partial melting. Plate tectonics as we know it probably developed later in the Archaean. The early Earth had much higher heat flow than in later times, and thus the lithosphere was more ductile rather than brittle. The essence of modern tectonics is a series of rigid plates that extend down to the asthenosphere. When they deform it is largely through brittle failure of the entire lithosphere.

How India accelerated towards Eurasia at the end of the Cretaceous

About 70 Ma ago the magnetic striping of the Indian Ocean floor suggests that the Indian subcontinent was then moving towards the huge, almost stationary Eurasian continent at about 8 cm per year. Over the next 5 Ma this convergence rate underwent a tectonically startling acceleration to reach 18 cm yr-1 by around the time of the Cretaceous-Palaeogene boundary (65 Ma): more than doubling the approach rate. Thereafter it slowed, eventually to a few centimetres per year once collision and building of the Himalayan mountain belt were more or less complete about 30 Ma ago. This cannot easily be explained by a speeding up of the sea-floor spreading rate at an Indian Ocean ridge to the south, 18 cm yr-1 being as fast as tectonic forces can manage at present. At that time ocean floor to the north of India was being subducted beneath Eurasia, and basaltic volcanism was flooding what is now the Deccan Plateau on western India. A couple of suggestions have been made: two northward subduction zones may have developed or the mantle plume feeding the Deccan flood basalts may have driven the tectonic acceleration. A third possibility is that the subduction was somehow lubricated. That approach has recently been considered by geoscientists from China and Singapore  (Zhou, H. et al. 2024. India–Eurasia convergence speed-up by passive-margin sediment subduction. Nature, v. 635, p. 114-120; DOI: 10.1038/s41586-024-08069-6).

Hao Zhou and colleagues studied the isotopic and trace-element geochemistry of volcanic and plutonic igneous complexes to the north of the Himalaya. They were emplaced in arc environments in three stages: from 98 to 89; 65 to 60; and 57 to 50 Ma. In this tectonic setting fluids rise from the subducted slab to induce the mantle part of the overriding lithosphere to partially melt. That yields magmas which penetrate the crust above. The first and last magmatic events produced similar isotopic and trace-element ‘signatures’, which suggest fluids rose from subducted ocean lithosphere.  But those in the latest Cretaceous to earliest Palaeocene are markedly different. Instead of showing signs of their magmas being entirely mantle derived like the earlier and later groups, the 65 to 60 Ma rocks exhibit clear evidence of partial melting having incorporated materials that had originated in older continental crust. The authors suggest that this crustal contamination stemmed from sediments that had been deposited at the northern margin of the Indian subcontinent during the Mesozoic. These sediments had formed by weathering of the ancient rocks that underpin India, transport of the debris by rivers and deposition on the seafloor as water-saturated sands, silts and clays. Once those sediments were subducted beneath what is now Tibet they would yield fluids with a geochemical ‘fingerprint’ inherited from old continental crust. Moreover, far more fluids than subducted oceanic crust could ever release would rise into the overriding lithosphere than.

The fluids rising from a subducted wedge of sediments may have reduced friction between the overriding Eurasian lithosphere and the subducted slab derived from the Indian tectonic plate. That scenario would not only have lubricated subduction, but allowed compressive forces in the overriding lithosphere to relax. Both would have allowed convergence of the two plates to move significantly faster as the sediments were progressively consumed. Once completed, convergence would have slowed without such ‘lubrication’.Earlier continent-continent collision zones, such as those that united Pangaea and older supercontinents may well have involved such tectonic surges. And the same kind of process may eventually speed up the reassembly of the latest distribution of continents.

Watch an animation of the India-Eurasia convergence (just over 3 minutes long)compiled by Christopher Scotese of Northwestern University in Evanston, Illinois, USA, which is a component of his Paleomap Project. It starts by following India from its current position to its origin in the break-up of Gondwanaland ~100 Ma ago. The last half reverses the motions to show India’s slow collision with Eurasia.

A new explanation for the Neoproterozoic Snowball Earth episodes

The Cryogenian Period that lasted from 860 to 635 million years ago is aptly named, for it encompassed two maybe three episodes of glaciation. Each left a mark on every modern continent and extended from the poles to the Equator. In some way, this series of long, frigid catastrophes seems to have been instrumental in a decisive change in Earth’s biology that emerged as fossils during the following Ediacaran Period (635 to 541 Ma). That saw the sudden appearance of multicelled organisms whose macrofossil remains – enigmatic bag-like, quilted and ribbed animals – are found in sedimentary rocks in Australia, eastern Canada and NW Europe. Their type locality is in the Ediacara Hills of South Australia, and there can be little doubt that they were the ultimate ancestors of all succeeding animal phyla. Indeed one of them Helminthoidichnites, a stubby worm-like animal, is a candidate for the first bilaterian animal and thus our own ultimate ancestor. Using the index for Palaeobiology or the Search Earth-logs pane you can discover more about them in 12 posts from 2006 to 2023. The issue here concerns the question: Why did Snowball Earth conditions develop? Again, refresh your knowledge of them, if you wish, using the index for Palaeoclimatology or Search Earth-logs. From 2000 onwards you will find 18 posts: the most for any specific topic covered by Earth-logs. The most recent are Kicking-off planetary Snowball conditions (August 2020) and Signs of Milankovich Effect during Snowball Earth episodes (July 2021): see also: Chapter 17 in Stepping Stones.

One reason why Snowball Earths are so enigmatic is that CO2 concentrations in the Neoproterozoic atmospheric were far higher than they are at present. In fact since the Hadean Earth has largely been prevented from being perpetually frozen over by a powerful atmospheric greenhouse effect. Four Ga ago solar heating was about 70 % less intense than today, because of the ‘Faint Young Sun’ paradox. There was a long episode of glaciation (from 2.5 to 2.2 Ga) at the start of the Palaeoproterozoic Era during which the Great Oxygenation Event (GOE) occurred once photosynthesis by oxygenic bacteria became far more common than those that produced methane. This resulted in wholesale oxidation to carbon dioxide of atmospheric methane whose loss drove down the early greenhouse effect – perhaps a narrow escape from the fate of Venus. There followed the ‘boring billion years’ of the Mesoproterozoic during which tectonic processes seem to have been less active. in that geologically tedious episode important proxies (carbon and sulfur isotopes) that relate to the surface part of the Earth System ‘flat-lined’.  The plethora of research centred on the Cryogenian glacial events seems to have stemmed from the by-then greater complexity of the Precambrian Earth System.

Since the GOE the main drivers of Earth’s climate have been the emission of CO2 and SO2 by volcanism, the sedimentary burial of carbonates and organic carbon in the deep oceans, and weathering. Volcanism in the context of climate is a two-edged sword: CO2 emission results in greenhouse warming, and SO2 that enters the stratosphere helps reflect solar radiation away leading to cooling. Silicate minerals in rocks are attacked by hydrogen ions (H+) produced by the solution of CO2 in rain water to form a weak acid (H2CO3: carbonic acid). A very simple example of such chemical weathering is the breakdown of calcium silicate:

CaSiO3  +  2CO2  + 3H2O  =  Ca2+  +  2HCO3  +  H4SiO4  

The reaction results in calcium and bicarbonate ions being dissolved in water, eventually to enter the oceans where they are recombined in the shells of planktonic organisms as calcium carbonate. On death, their shells sink and end up in ocean-floor sediments along with unoxidised organic carbon compounds. The net result of this part of the carbon cycle is reduction in atmospheric CO2 and a decreased greenhouse effect: increased silicate weathering cools down the climate. Overall, internal processes – particularly volcanism – and surface processes – weathering and carbonate burial – interact. During the ‘boring billion’ they seem to have been in balance. The two processes lie at the core of attempts to model global climate behaviour in the past, along with what is known about developments in plate tectonics – continental break-up, seafloor spreading and orogenies – and large igneous events resulting from mantle plumes. A group of geoscientists from the Universities of Sydney and Adelaide, Australia have evaluated the tectonic factors that may have contributed to the first and longest Snowball Earth of the Neoproterozoic: the Sturtian glaciation (717 to 661 Ma) (Dutkiewicz, A. et al. 2024. Duration of Sturtian “Snowball Earth” glaciation linked to exceptionally low mid-ocean ridge outgassing. Geology, v. 52, online early publication; DOI: 10.1130/G51669.1).

Palaeogeographic reconstructions (Robinson projection) during the early part of the Sturtian global glaciation: LEFT based on geological data from Neoproterozoic terrains on modern continents; RIGHT based on palaeomagnetic pole positions from those terrains. Acronyms refer to each terrains, e.g. Am is Amazonia, WAC is the West African Craton. Orange lines are ocean ridges, those with teeth are subduction zone. (Credit: Dutkiewicz et al., parts of Fig. 1)

Shortly before the Sturtian began there was a major flood volcanism event, forming the Franklin large igneous province, remains of which are in Arctic Canada. The Franklin LIP is a subject of interest for triggering the Sturtian, by way of a ‘volcanic winter’ effect from SO2 emissions or as a sink for CO through its weathering. But both can be ruled out as no subsequent LIP is associated with global cooling and the later, equally intense Marinoan global glaciation (655 to 632 Ma) was bereft of a preceding LIP. Moreover, a world of growing frigidity probably could not sustain the degree of chemical weathering to launch a massive depletion in atmospheric CO2. In search of an alternative, Adriana Dutkiewicz and colleagues turned to the plate movements of the early Neoproterozoic. Since 2020 there have been two notable developments in modelling global tectonics of that time, which was dominated by the evolution of the Rodinia supercontinent. One is based largely on geological data from the surviving remnants of Rodinia (download animation), the other uses palaeomagnetic pole positions to fix their relative positions: the results are very different (download animation).

Variations in ocean ridge lengths, spreading rates and oceanic crust production during the Neoproterozoic estimated from the geological (orange) and palaeomagnetic (blue) models. Credit: Dutkiewicz et al., parts of Fig. 2)

The geology-based model has Rodinia beginning to break up around 800 Ma ago with a lengthening of global constructive plate margins during disassembly. The resulting continental drift involved an increase in the rate of oceanic crust formation from 3.5 to 5.0 km2 yr-1. Around 760 Ma new crust production more than halved and continued at a much slowed rate throughout the Cryogenian and the early part of the Ediacaran Period.  The palaeomagnetic model delays breakup of the Rodinia supercontinent until 750 Ma, and instead of the rate of crust production declining through the Cryogenian it more than doubles and remains higher than in the geological model until the late Ediacaran. The production of new oceanic crust is likely to govern the rate at which CO2 is out-gassed from the mantle to the atmosphere. The geology-based model suggests that from 750 to 580 Ma annual CO2 additions could have been significantly below what occurred during the Pleistocene ice ages since 2.5 Ma ago. Taking into account the lower solar heat emission, such a drop is a plausible explanation for the recurrent Snowball Earths of the Neoproterozoic. On the other hand, the model based on palaeomagnetic data suggests significant warming during the Cryogenian contrary to a mass of geological evidence for the opposite.

A prolonged decrease in tectonic activity thus seems to be a plausible trigger for global glaciation. Moreover, reconstruction of Precambrian global tectonics using available palaeomagnetic data seems to be flawed, perhaps fatally. One may ask, given the trends in tectonic data: How did the Earth repeatedly emerge from Snowball episodes? The authors suggest that the slowing or shut-down of silicate weathering during glaciations allowed atmospheric CO2 to gradually build up as a result of on-land volcanism associated with subduction zones that are a quintessential part of any tectonic scenario.

This kind of explanation for recovery of a planet and its biosphere locked in glaciation is in fact not new. From the outset of the Snowball Earth hypothesis much the same escape mechanisms were speculated and endlessly discussed. Adriana Dutkiewicz and colleagues have fleshed out such ideas quite nicely, stressing a central role for tectonics. But the glaring disparities between the two models show that geoscientists remain ‘not quite there’. For one thing, carbon isotope data from the Cryogenian and Ediacaran Periods went haywire: living processes almost certainly played a major role in the Neoproterozoic climatic dialectic.

News about when subduction began

Tangible signs of past subduction take the form of rocks whose mineralogy shows that they have been metamorphosed under conditions of high pressure and low temperature, and then returned to the surface somehow. Ocean-crust basaltic rocks become blueschist and eclogite. The latter is denser than mantle peridotite so that oceanic lithosphere can sink and be recycled. That provides the slab-pull force, which is the major driver of plate tectonics. Unfortunately, neither blueschists nor eclogites are found in metamorphic complexes older than about 800 Ma. This absence of direct proof of subduction and thus modern style plate tectonics has resulted in lively discussion and research seeking indirect evidence for when it did begin, the progress of which since 2000 you can follow through the index for annual logs about tectonics. An interesting new approach emerged in 2017 that sought a general theory for the evolution of silicate planets, which involves the concept of ‘lid tectonics’. A planet in a stagnant-lid phase has a lithosphere that is weak as a result of high temperatures: indeed so weak and warm that subduction was impossible. Stagnant-lid tectonics does not recycle crustal material back to its source in the mantle and it simply builds up the lithosphere. Once planetary heat production wanes below a threshold level that permits a rigid lithosphere, parts of the lid can be driven into the mantle. The beginnings of this mobile-lid phase and thus plate tectonics of some kind involves surface materials in mantle convection: the may be recycled.

Cartoon of possible Hadean stagnant lid tectonics, dominated by mantle plumes. (Credit: Bédard, J.H. 2018, Fig 3B, DOI: 10.1016/j.gsf.2017.01.005)

A group of geochemists from China, Canada and Australia have sought evidence for recycled crustal rocks from silicon and oxygen isotopes in the oldest large Archaean terrane, the  4.0 Ga old Acasta Gneiss Complex in northern Canada (Zhang, Q. and 10 others 2023. No evidence of supracrustal recycling in Si-O isotopes of Earth’s oldest rocks 4 Ga ago. Science Advances, v.9, article eadf0693; DOI: 10.1126/sciadv.adf0693). Silicon has three stable isotopes 28Si, 29Si, and 30Si. As happens with a number of elements, various geochemical processes are able to selectively change the relative proportions of such isotopes: a process known as isotope fractionation. As regards silicon isotopes used to chart lithosphere recycling, the basic steps are as follows: Organisms that now remove silicon from solution in seawater to form their hard parts and accumulate in death as fine sediments like flint had not evolved in the Archaean. Because of that reasonable supposition it has been suggested that seawater during the Archaean contained far more dissolved silicon than it does now. Such a rich source of Si would have entered Archaean oceanic crust and ocean-floor sediments to precipitate silica ‘cement’. The heaviest isotope 30Si would have left solution more easily than the lighter two. Should such silicified lithosphere have descended to depths in the mantle where it could partially melt the anomalously high 30Si would be transferred to the resulting magmas.

Proportions of 30Si in zircons, quartz and whole rock for Acasta gneisses (coloured), other Archaean areas (grey) and Jack Hills zircons (open circles. Vertical lines are error bars. (Credit: simplified from Zhang et al. Fig 1)

Stable-isotope analyses by Zhang et al. revealed that zircon and quartz grains and bulk rock samples from the Acasta gneisses, with undisturbed U-Pb ages, contain 30Si in about the same proportions relative to silicon’s other stable isotopes as do samples of the mantle. So it seems that the dominant trondhjemite-tonalite-granodiorite (TTG) rocks that make up the oldest Acasta gneisses were formed by partial melting of a source that did not contain rocks from the ocean crust. Yet the Acasta Gneiss Complex also contains younger granitic rocks (3.75 to 3.50 Ga) and they are significantly more enriched in 30Si, as expected from a deep source that contained formerly oceanic rocks. A similar ‘heavy’ silicon-isotope signature is also found in samples from other Archaean terranes that are less than 3.8 Ga old. Thus a major shift from stagnant-lid tectonics to the mobile-lid form may have occurred at the end of the Hadean. But apart from the Acasta Gneiss Complex only one other, much smaller Hadean terrane has been discovered, the 4.2 Ga Nuvvuagittuq Greenstone Belt. It occupies a mere 20 km2 on the eastern shore of Hudson Bay in Canada, and appears to be a sample of Hadean oceanic crust. It does include TTG gneisses, but they are about 3.8 Ga old and contain isotopically heavy silicon. So it seems unlikely that testing this hypothesis with silicon-isotope data from other Hadean gneissic terranes will be possible for quite a while, if at all.

Did Precambrian BIFs ‘fall’ into the mantle to trigger mantle plumes?

How the Earth has been shaped has depended to a large extent on a very simple variable among rocks: their density. Contrasts in density between vast rock masses are expressed when gravity attempts to maintain a balance of forces. The abrupt difference in elevation of the solid surface at the boundaries of oceans and continents – the Earth’s hypsometry – stems from the contrasted densities of continental and oceanic crust: the one dominated by granitic rocks (~2.8 t m-3) the other by those of basaltic composition (~ 3.0 t m-3). Astronomers have estimated that Earth’s overall density is about 5.5 t m-3 – it is the densest planet in the Solar System. The underlying mantle makes up 68% of Earth’s mass, with a density that increases with depth from 3.3 to 5.4 t m-3 in a stepwise fashion, at a number of discontinuities, because mantle minerals undergo changes induced by pressure. The remaining one third of Earth’s mass resides in the iron-nickel core at densities between 9.5 to 14.5 t m-3. Such density layering is by no means completely stable. Locally increased temperatures in mantle rocks reduce their density sufficiently for masses to rise convectively to be replaced by cooler ones, albeit slowly. By far the most important form of convection affecting the lithosphere involves the resorption of oceanic lithosphere plates at destructive margins, which results in subduction. This is thought to be due to old, cold oceanic basalts undergoing metamorphism as pressure increases during subduction. They are transformed at depth to a mineral assemblage (eclogite) that is denser (3.4 to 3.5 t m-3) than the enveloping upper mantle. That density contrast is sufficient for gravity to pull slabs of oceanic lithosphere downwards. This slab-pull force is transmitted through oceanic lithosphere that remains at the surface to become the dominant driver of modern plate tectonics. As a result, extension of the surface oceanic lithosphere at constructive margins draws mantle upwards to partially melt at reduced pressure, thus adding new basaltic crust at mid-ocean rift systems to maintain a form of mantle convection. Seismic tomography shows that active subducted slabs become ductile about 660 km beneath the surface and below that no earthquakes are detected. Quite possibly, the density of the reconstituted lithospheric slab becomes less than that of the mantle below the 660 km discontinuity. So the subducted slab continues by moving sideways and buckling in response to the ‘push’ from its rigid upper parts above. But it has been suggested that some subducted slabs do finally sink to the core-mantle boundary, but that is somewhat conjectural.

Typical banded iron formation

There are sedimentary rocks whose density at the surface exceeds that of the upper mantle: banded iron formations (BIFs) that contain up to 60% iron oxides (mainly Fe2O3) and have an average density at the surface of around 3.5 t m-3. BIFs formed mainly in the late Archaean and early Proterozoic Eons  (3.2 to 1.0 Ga) and none are known from the last 400 Ma. They formed when soluble iron-2 (Fe2+) – being added to ocean water by submarine hydrothermal activity –was precipitated as Fe3+ in the form of iron oxide (Fe2O3) where oxygen was present in ocean water. With little doubt this happened only in shallow marine basins where cyanobacteria that appeared about 3.5 Ga ago had sufficient sunlight to photosynthesise. Until about 2.4 Ga the atmosphere and thus the bulk of ocean water contained very little oxygen so the oceans were pervaded by soluble iron so that BIFs were able to form wherever such biological activity was going on. Conceivably (but not proven), that BIF-forming biochemical reaction may even have operated far from land in ocean surface water, slowly to deposit Fe2O3 on the deep ocean floor. After 2.4 Ga oxygen began to build in the atmosphere after the Great Oxidation Event had begon. That time was also when the greatest production of BIFs took place. Strangely, the amount of BIF in the geological record fell during the next 600 Ma to rise again to a very high peak at 1.8 Ga. Since there must have been sufficient soluble iron and an increasing amount of available oxygen for BIFs to form throughout that ‘lean’ period the drop in BIF formation is paradoxical. After 1.0 Ga BIFs more or less disappear. By then so much oxygen was present in the atmosphere and from top to bottom in ocean water that soluble iron was mostly precipitated at its hydrothermal source on the ocean floor. Incidentally, modern ocean surface water far from land contains so little dissolved iron that little microbiological activity goes on there: iron is an essential nutrient so the surface waters of remote oceans are effectively ‘wet deserts’.

Plots of probability of LIPs and BIFs forming at the Earth’s surface during Precambrian times, based on actual occurrences (Credit: Keller, et al., modified Fig 1A)

Spurred by the fact that if a sea-floor slab dominated by BIFs was subducted it wouldn’t need eclogite formation to sink into the mantle, Duncan Keller of Rice University in Texas and other US and Canadian colleagues have published a ‘thought experiment’ using time-series data on LIPs and BIFs compiled by other geoscientists (Keller, D.S. et al. 2023. Links between large igneous province volcanism and subducted iron formations. Nature Geoscience, v. 16, article; DOI: 10.1038/s41561-023-01188-1.). Their approach involves comparing the occurrences of 54 BIFs through time with signs of activity in the mantle during the Palaeo- and Mesoproterozoic Eras, as marked by large igneous provinces (LIPs) during that time span. To do this they calculated the degree of correlation in time between BIFs and LIPs. The authors chose a minimum area for LIPs of 400 thousand km2 – giving a total of 66 well-dated examples. Because the bulk of Precambrian flood-basalt provinces, such as occurred during the Phanerozoic, have been eroded away, most of their examples are huge, well-dated dyke swarms that almost certainly fed such plateau basalts. Rather than a direct time-correlation, what emerged was a match-up that covered 74% of the LIPs with BIFs that had formed about 241 Ma earlier. They also found a less precise correlation between LIPs associated with 241 Ma older BIFs and protracted periods of stable geomagnetic field, known as ‘superchrons’. These are thought by geophysicists to be influenced by heat flow through the core-mantle boundary (CMB).

The high bulk density of BIFs at the surface would be likely to remain about 15 % greater than that of peridotite as pressure increased with depth in the mantle. Such slabs could therefore penetrate the 660 mantle discontinuity. Their subduction would probably result in their eventually ‘piling up’ in the vicinity of the CMB. The high iron content of BIFs may also have changed the way that the core loses heat, thereby triggering mantle plumes. Certainly, there is a complex zone of ultra-low seismic velocities (ULVZ) that signifies hot, ductile material extending above the CMB. Because BIFs’ high iron-content makes them thermally highly conductive compared with basalts and other sediments, they may be responsible. Clearly, Keller et al’s hypothesis is likely to be controversial and they hope that other geoscientists will test it with new or re-analysed geophysical data. But the possibility of BIFs falling to the base of the mantle spectacularly extends the influence of surface biological processes to the entire planet. And, indeed, it may have shaped the later part of its tectonic history having changed the composition of the deep mantle. The interconnectedness of the Earth system also demands that the consequences – plumes and large igneous provinces – would have fed back to the Precambrian biosphere. See also: Iron-rich rocks unlock new insights into Earth’s planetary history, Science Daily, 2 June 2023

The Earth System in action: land plants affected composition of continental crust

The essence of the Earth System is that all processes upon, above and beneath the surface interact in a bewildering set of connections. Matter and energy in all their forms are continually being exchanged, deployed and moved through complex cycles: involving rocks and sediments; water in its various forms; gases in the atmosphere; magmas; moving tectonic plates and much else besides. The central and massively dominant role of plate tectonics connects surface processes with those of our planet’s interior: the lithosphere, mantle and, arguably, the core. Interactions between the Earth System’s components impose changes in the dynamics and chemical processes through which it operates. Living processes have been a part of this for at least 3.5 billion years ago, in part through their role in the carbon cycle and thus the Earth’s climatic evolution. During the Silurian Period life became a pervasive component of the continental surface, first in the form of plants, to be followed by animals during the Devonian Period. Those novel changes have remained in place since about 430 Ma ago, plants being the dominant base of continental ecosystems and food chains.

Schematic diagram showing changes in river systems and their alluvium before and after the development of land plants. (Credit: Based on Spencer et al. 2022, Fig 4)

Land plants exude a variety of chemicals from their roots that break down rock to yield nutrient elements. So they play a dominant role in the formation of soil and are an important means of rock weathering and the production of clay minerals from igneous and metamorphic minerals. Plant root systems bind near-surface sediments thus increasing their resistance to erosion by wind and water, and to mass movement under gravity. This binding and plant canopies efficiently reduce dust transport, slow water flow on slopes and decrease the sediment load of flowing water. Plants and their roots also stabilise channels systems. There is much evidence that before the Devonian most rivers comprised continually migrating braided channels in which mainly coarse sands and gravels were rapidly deposited while silts and muds in suspension were shifted to the sea. Thereafter flow became dominated by larger and fewer channels meandering across wide tracts on which fine sediment could accumulate as alluvium on flood plains when channels broke their banks. Land plants more efficiently extract CO2 from the atmosphere through photosynthesis and the new regime of floodplains could store dead plant debris in the muds and also in thick peat deposits. As a result, greenhouse warming had dwindled by the Carboniferous, encouraging global cooling and glaciation. 

Judging the wider influence of the ‘greening of the land’ on other parts of the Earth system, particularly those that depend on internal  magmatic processes, relies on detecting geochemical changes in minerals formed as direct outcomes of plate tectonics. Christopher Spencer of Queen’s University in Kingston, Canada and co-workers at the Universities of Southampton, Cambridge and Aberdeen in the UK, and the China University of Geosciences in Wuhan set out to find and assess such a geochemical signal (Spencer, C., Davies, N., Gernon, T. et al. 2022. Composition of continental crust altered by the emergence of land plants. Nature Geoscience, v. 15 online publication; DOI: 10.1038/s41561-022-00995-2). Achieving that required analyses of a common mineral formed when magmas crystallise: one that can be precisely dated, contains diverse trace elements and whose chemistry remains little changed by later geological events. Readers of Earth-logs might have guessed that would be zircon (ZrSiO). Being chemically unreactive and hard, small zircon grains resist weathering and the abrasion of transport to become common minor minerals in sediments. Thousands of detrital zircon grains teased out from sediments have been dated and analysed in the last few decades. They span almost the entirety of geological history. Spencer et al. compiled a database of over 5,000 zircon analyses from igneous rocks formed at subduction zones over the last 720 Ma, from 183 publications by a variety of laboratories.

The approach considered two measures: the varying percentages of mudrocks in continental sedimentary sequences since 600 Ma ago; aspects of the hafnium- (Hf) and oxygen-isotope proportions measured in the zircons using mass spectrometry and their changes over the same time. Before ~430 Ma the proportion of mudrocks in continental sedimentary sequences is consistently much lower than it is in post post-Silurian, suggesting a link with the rise of continental plant cover (see second paragraph). The deviation of the 176Hf/177Hf ratio in an igneous mineral from that of chondritic meteorites (the mineral’s εHf value) is a guide to the source of the magma, negative values indicating a crustal source, whereas positive values suggest a mantle origin. The relative proportions of two oxygen isotopes 18O and 16O  in zircons, expressed as δ18O, indicates the proportion of products of weathering, such as clay minerals, involved in magma production – 18O selectively moves from groundwater to clay minerals when they form, increasing their δ18O.

While the two geochemical parameters express very different geological processes, the authors noticed that before ~430 Ma the two showed low correlation between their values in zircons. Yet, surprisingly, the parameters showed a considerable and consistent increase in their correlation in younger zircons, directly paralleling the ‘step change’ in the proportions of mudstones after the Silurian. Complex as their arguments are, based on several statistical tests, Spencer et al. conclude that the geologically sudden change in zircon geochemistry ultimately stems from land plants’ stabilisation of river systems. As a result more clay minerals formed by protracted weathering, increasing the δ18O in soils when they were eroded and transported. When the resulting marine mudrocks were subducted they transferred their oxygen-isotope proportions to magmas when they were partially melted.

That bolsters the case for dramatic geological consequences of the ‘greening of the land’. But did its effect on arc magmatism fundamentally change the bulk composition of post-Silurian additions to the continental crust? To be convinced of that I would like to see if other geochemical parameters in subduction-related magmas changed after 430 Ma. Many other elements and isotopes in broadly granitic rocks have been monitored since the emergence of high-precision rock-analysing technologies around 50 years ago. There has been no mention, to my knowledge, that the late-Silurian involved a magmatic game-changer to match that which occurred in the Archaean, also revealed by hafnium and oxygen isotopes in much more ancient zircons.   

See also: https://www.sci.news/othersciences/geoscience/land-plants-continental-crust-composition-11151.htmlhttps://www.eurekalert.org/news-releases/963296

Evidence for an early Archaean transition to subduction

Modern plate tectonics is largely driven by slab-pull: a consequence of high-pressure, low-temperature metamorphism of the oceanic crust far from its origin at an oceanic ridge. As it ages, basaltic crust cools, become increasingly hydrated by hydrothermal circulation of seawater through it and its density increases. That is why the abyssal plains of the ocean floor are so deep relative to the shallower oceanic ridges where it formed. Due to the decrease in the Earth’s internal heat production by decay of radioactive isotopes, once oceanic lithosphere breaks and begins to descend high-P low-T metamorphism transforms the basaltic crust to a denser form: eclogite, in which the dense, anhydrous minerals garnet and sodium-rich pyroxene (omphacite) form. Depending on local heat flow, the entire oceanic slab may then exceed the density of the upper mantle to drag the plate downwards under gravity. Metamorphic reactions of any P-T regime creates minerals less capable of holding water and drive H2O-rich fluids upwards into the overriding lithosphere, thus inducing it to partially melt. Magmas produced by this create volcanism at the surface, either at oceanic island arcs or near to continental margins, depending on the initial position of the plate subduction.

A direct proof of active subduction in the geological record is the presence of eclogite and related blueschists. Such rocks are unknown before 2100 Ma ago (mid-Palaeoproterozoic of the Democratic Republic of Congo) but there are geochemical means of ‘sensing’ plate tectonic control over arc magmatism (See: So, when did plate tectonics start up? February 2016).  The relative proportions of rare-earth elements in ancient magmatic rocks that make up the bulk of continental crust once seemed to suggest that plate tectonics started at the end of the Archaean Eon (~2500 Ma). That method, however, was quite crude and has been superseded by looking in great detail at the geochemistry of the Earth’s most durable mineral: zircon (ZrSiO4), which began more than two decades ago. Minute grains of that mineral most famously have pushed back the geological record into what was long believed to be half a billion years with no suggestion of a history: the Hadean. Zircon grains extracted from a variety of ancient sediments have yielded U-Pb ages of their crystallisation from igneous magma that extend back 4.4 billion years (Ga) (see: Pushing back the “vestige of a beginning”;January 2001).  

Though simple in their basic chemical formula, zircons sponge-up a large range of other trace elements from their parent magma. So, in a sense, each tiny grain is a capsule of their geochemical environment at the time they crystallised. In 2020 Australian geochemists presented the trace-element geochemistry of 32 zircons extracted from a 3.3 Ga old sedimentary conglomerate in the Jack Hills of Western Australia, which lie within an ancient continental nucleus or craton. They concluded that those zircons mainly reveal that they formed in andesitic magmas, little different from the volcanic rocks that are erupted today above subduction zones. From those data it might seem that some form of plate tectonics has been present since shortly after the Earth’s formation. Oxygen-isotope data from zircons are useful in checking whether zircons had formed in magmas derived directly from partial melting of mantle rocks or by recycling of crustal magmatic rocks through subduction. Such a study in 2012 (see: Charting the growth of continental crust; March 2012) that used a very much larger number of detrital zircon grains from Australia, Eurasia, North America, and South America seemed, in retrospect, to contradict a subduction-since-the-start view of Earth dynamics and crust formation. Instead it suggested that recycling of crust, and thus plate-tectonic subduction, first showed itself in zircon geochemistry at about 3 Ga ago.

Detailed chemical and isotopic analysis of zircons using a variety of instruments has steadily become faster and cheaper. Actually finding the grains is much easier than doing interesting things with them. It is a matter of crushing the host rock to ‘liberate’ the grains. Sedimentary hosts that have not been strongly metamorphosed are much more tractable than igneous rocks. Being denser than quartz, the dominant sedimentary mineral, zircon can be separated from it along with other dense, trace minerals, and from them in turn by various methods based on magnetic and electrical properties. Zircons can then be picked out manually because of their distinctive colours and shapes. A tedious process, but there are now several thousand fully analysed zircons aged between 3.0 to 4.4 Ga, from eleven cratons that underpin Australia, North America, India, Greenland and southern Africa. The latest come from a sandstone bed laid down about 3.31 Ga ago in the Barberton area of South Africa (Drabon, N. et al. 2022. Destabilization of Long‐Lived Hadean Protocrust and the Onset of Pervasive Hydrous Melting at 3.8 GaAGU Advances, v. 3, article e2021AV000520; DOI: 10.1029/2021AV000520). The authors measured lutetium (Lu), hafnium (Hf) and oxygen isotopes, and concentrations of a suite of trace element in 329 zircons from Barberton dated between 3.3 to 4.15 Ga.

A schematic model of transition from Hadean-Eoarchaean lid tectonics to a type of plate tectonics that subsequently evolved to its current form, based on hafnium isotope data in ancient zircons (credit: Bauer et al. 2020; Fig 3)

The Hf isotopes show two main groups relative to the values for chondritic meteorites (assumed to reflect the composition of the bulk Earth). Zircons dated between 3.8 and 4.15 Ga all show values below that expected for the whole Earth. Those between 3.3 and 3.8 Ga show a broader range of values that extend above chondritic levels. The transition in data at around 3.8 Ga is also present in age plots of uranium relative to niobium and scandium relative to ytterbium, and to a lesser extent in the oxygen isotope data. On the basis of these data, something fundamentally changed in the way the Earth worked at around 3.8 Ga. Nadja Drabon and colleagues ascribe the chemical features of Hadean and Eoarchaean zircons to an early protocrust formed by melting of chemically undepleted mantle. This gradually built up and remained more or less stable for more than 600 Ma, without being substantially remelted through recycling back to mantle depths. After 3.8 billion years ago, geochemical signatures of the zircons start showing similarities to those of zircons derived from modern subduction zones. Hf isotopes and trace-element geochemistry in 3.6 to 3.8 Ga-old  detrital zircons from other cratons are consistent with a 200 Ma transition from ‘lid’ tectonics (see: Lid tectonics on Earth; December 2017) to the familiar tectonics of rigid plates whose basalt-capped lithosphere ultimately returns to the mantle to be involved in formation of new magmas from which continental crust stems. Parts of plates bolstered by this new, low density crust largely remain at the surface.

While Drabon et al. do provide new data from South Africa’s Kaapvaal craton, their conclusions are similar to earlier work by other geochemists based on data from other area (e.g. Bauer, A.M. et al. 2020. Hafnium isotopes in zircons document the gradual onset of mobile-lid tectonicsGeochemical Perspectives Letters, v. 14; DOI: 10.7185/geochemlet.2015), which the accompanying figure illustrates.

See also: Earliest geochemical evidence of plate tectonics found in 3.8-billion-year-old crystal. Science Daily, 21 April 2022. 3.8-Billion-Year-Old Zircons Offer Clues to When Earth’s Plate Tectonics Began. SciNews, 26 April 2022

New ideas on how subduction works

Nowadays, plate tectonics is thought mainly to be driven by the sinking of old, relatively cold and dense oceanic lithosphere at subduction zones: slab-pull force dominates the current behaviour of the outermost Earth. At the eastern edge of Eurasia subduction beneath Japan has yet to consume Pacific Ocean lithosphere younger than 180 Ma (Middle Jurassic). The Pacific Plate extends eastwards from there for over 7000 km to its source at the East Pacific Rise. That spreading axis has disappeared quite recently beneath the North American Plate between Baha California and northern California. It has been subducted. Since, to a first approximation, sea-floor spreading is at the same pace either side of mid-ocean constructive plate margins, subduction at the western edge of the North America has consumed at least 7000 km of old ocean lithosphere. Slab-pull force there has been sustained for probably more than 250 Ma. As a result several former island arcs have been plastered onto the leading edge of the North American Plate to create the geological complexity of its western states. If at any time the weight of the subducting slab had caused it leading edge literally to snap and fall independently wouldn’t that have decreased slab-pull force or shut it off, and spreading at the East Pacific Rise, altogether? No, says the vast expanse of the West Pacific plate

That dichotomy once encouraged scientists of the plate-tectonic era to assume that a subducted slab remains as strong as rigid plates at the surface. They believed that subduction merely bends a plate so that it can slide into the mantle. The use of seismic waves (seismic tomography) to peer into the mantle has revealed a far more complex situation. Beneath North America traces of subducted slabs are highly deformed and must have lost their rigidity, yet they still maintain slab-pull force. Three geoscientists from the Swiss Federal Institute of Technology Zurich, Switzerland, and the University of Texas at Austin, USA (Gerya T. V., Becovici, D. & Becker, T.W. 2021. Dynamic slab segmentation due to brittle–ductile damage in the outer rise. Nature, v. 599, p 245-250; DOI: 10.1038/s41586-021-03937-x) used computer-generated models of how various forces and temperature conditions at small and large scales bear on the behaviour of slabs being subducted. Where a plate bends into a subduction zone its rigidity results in cracking and faulting of its no convex upper surface, while the base is compressed. Seismic anomalies in the descending slab reflect the formation of pulled-apart segments, similar to those in a bar of chocolate (for a possible example from an exhumed subduction zone see: A drop off the old block? May 2008). Thermo-mechanical modelling suggests that the slab becomes distinctly weakened through brittle damage and by reduction in grain size because of ductile deformation, yet each segment maintains a high viscosity relative to the surrounding mantle rocks. Under present conditions and those extrapolated back into the Proterozoic, where the slab is thinned between segments it remains sufficiently viscous to avoid segments detaching to sink independently of one another. Such delamination would reduce slab-pull force. Another process operates in the surrounding mantle. The occurrence of earthquakes in a subducted slab down to a depth of about 660 km – the level of a major discontinuity in the mantle where pressure induces a change in its mineralogy and density – confirms that a modern slab maintains some rigidity and deforms in a brittle fashion. But at this depth it cannot continue to descend steeply and travels horizontally along the discontinuity, pushed by the more shallow subduction. It can now become buckled as the mantle resists its lateral motion.

Left: the subduction zone beneath Japan defined by seismic tomography (yellow to red = lower seismic wave speeds – more ductile; yellow to blue = higher speeds – more rigid). Right: modelled evolution of viscosity in a similar subduction zone under modern conditions showing slab segmentation (blue to brown = increasing viscosity). (Credit: Gerya et al., Figs 4c & 1a-e)

Rather than trying to mimic the chaos beneath North America the authors compared their results with seismic tomography of the younger system of westward subduction beneath Japan. This allowed them to ‘calibrate’ their modelling against actual deep structure well-defined by seismic tomography. The tectonic jumble beneath North America probably resulted from a much longer history of eastwards subduction. The complexity there may be explained by successive foundering of deformed slabs into the deeper mantle looking a bit like a sheet of still viscous pie pastry dropped on its edge. This happened, perhaps, as island arcs that had formed in the eastern Pacific sporadically accreted to the continent as the intervening oceanic lithosphere was subducted.   

There is ample evidence that modern-style subduction was widespread back as far as the Palaeoproterozoic. But in the Archaean the evidence is fitful: some hints of subduction, but plenty of contrary evidence.  Gerya and co-workers suggest that higher heat production from radioactive decay mantle earlier in Earth’s history would have reduced plate strength and mantle resistance to slab penetration. Subduction may have occurred but was interrupted repeatedly by foundering/delamination of individual detached segments at much shallower depths. That implies weaker as well as intermittent slab pull, or even further back its complete absence, so that planetary recycling would then have required other mechanisms, such as ‘drip tectonics’.

See also: Crushed resistance: Tectonic plate sinking into a subduction zone and Fate of sinking tectonic plates is revealed, Science Daily, 11 November 2021

Nappe tectonics at the end of the Archaean

The beginning of modern-style plate tectonics is still debated in the absence of definite evidence. Because Earth’s mantle generates heat through radioactive decay and still contains heat left over from planetary accretion and core formation it must always have maintained some kind of heat transfer through some kind of circulatory motion involving the mantle and lithosphere. That must always too have involved partial melting and chemical differentiation that created materials whose density was lower than that of the mantle; e.g. continental crust. Since continental materials date back to more than 4 billion years ago and some may have been generated earlier in the Hadean, only to be lrgely resorbed, a generalised circulation and chemical differentiation have been Earth’s main characteristics from the start. One view is that early circulation was a form of vertical tectonics without subduction via a sort of ‘dripping’ or delamination of particularly dense crustal materials back into the mantle. A sophisticated model of how the hotter early Earth worked in this way has been called ‘lid tectonics’, from which plate tectonics evolved as the Earth cooled and developed a thicker, more rigid lithosphere. Such an outer layer would be capable of self-generating the slab pull that largely drives lateral motions of lithospheric plates. That process occurs once a slab of oceanic lithosphere becomes cool and dense enough to be subducted (see: How does subduction start?; August 2018).

The most convincing evidence for early plate tectonics would therefore be tangible signs of both subduction and large horizontal movements of lithospheric plates: common enough in the Neoproterozoic and Phanerozoic records, but not glaringly obvious in the earlier Archaean Eon. These unequivocal hallmarks have now emerged from studies of Archaean rocks in the Precambrian basement that underpins northern China and North Korea. The North China Craton has two main Archaean components: an Eastern Block of gneisses dated between 3.8 and 3.0 Ga and a Western Block of younger (2.6 to 2.5 Ga) gneisses, metavolcanics and metasediments. They are separated by a zone of high deformation. A key area for understanding the nature of the deformed Central Orogenic Belt is the Zanhuan Complex near the city of Kingtai (Zhong, YL. et al. 2021. Alpine-style nappes thrust over ancient North China continental margin demonstrate large Archean horizontal plate motions. Nature  Communications, v. 12, article6172, DOI: 10.1038/s41467-021-26474-7).

Schematic cross sections through the Zanhuan Complex of northern China, showing early and final development of the Central Orogenic Belt in the North China Block . (Credit: Zhong, YL. et al.;Figs 10b and c)

This small, complex area reveals that the older Eastern Block is unconformably overlain by Neoarchaean sediments, above which has been thrust a stacked series of nappes similar in size and form to those of the much younger Alpine orogenic belt of southern Europe. Though highly complex, the rocks involved having been folded and stretched by ductile processes, they are still recognisable as having originally been at the surface. Metavolcanics in the nappes can be assigned from their geochemistry to a late-Archaean fore-arc, through comparison with that of modern igneous rocks formed at such a setting in the Western Pacific. Thrust over the nappe complex is a jumble or mélange of highly deformed metasediments containing blocks of metabasalts and occasional ultramafic igneous rocks that geochemically resemble oceanic crust formed at a mid-ocean ridge. Some of them contain high-pressure minerals formed at depth in the mantle, indicating that they had once been subducted. The whole complex is cut by undeformed dykes of granitic composition dated at 2.5 Ga, confirming that the older rocks and the structures within them are Archaean in age. Thrust over the melange and tectonically underlying nappe complex are less-deformed volcanic rocks and granitic intrusions that closely resemble what is generally found in modern island arcs.

Orogenic belts bear witness to enormous crustal shortening caused by horizontal compressive forces. Assuming the average rate of modern subduction (2 cm yr-1) the 178 Ma history of the Zanhuan Complex implies more than 3,500 km of lateral transport. 2.5 billion years ago, higher radioactive heat production in the mantle would have made tectonic overturning considerably faster  The unconformity at the base of the complex suggests that it was driven over the equivalent of a modern passive, continental margin. So the complex provides direct evidence of horizontal plate tectonics and associated subduction during the latter stages of the Archaean that ranks in scale with that of many Phanerozoic orogenic belts, such as that of the European Alps. The Zanhuan Complex is a result of arc accretion that played a major role in many later orogens. The North China craton itself is reminiscent of continent-continent collision, as required in the formation of supercontinents.

Subduction and continental collision in the Himalaya

The Indian subcontinent after it separated from Madagascar in the Late Cretaceous to move northwards to its destined collision with Eurasia and the formation of the Himalaya. (Credit: Frame from an animation ©Christopher Scotese)

During the Early Cretaceous (~140 Ma ago) India, Madagascar, Antarctica and Australia parted company with Africa after 400 Ma of unity as components of the Gondwana supercontinent. By 120 Ma Antarctica and Australia split from India and Madagascar, and the Indian Ocean began to form. India moved northwards , leaving Madagascar in its wake after about 70 Ma ago. By 50 Ma the subcontinent began to collide with Eurasia, its northward motion driving before it crustal materials that eventually formed the Himalaya. This highly complex process is wonderfully documented in an animation made in 2015 by Christopher Scotese, Emeritus Professor in the Department of Earth and Environmental Sciences, Northwestern University, USA. At the start of its journey India moved northwards at a slow rate of about 5 cm per year. After 80 Ma it speeded up dramatically to 15 cm per year, about twice as fast as any modern continental drift and a pace that lasted for over 30 Ma until collision began. How could that, in a geological sense, sudden and sustained acceleration have been induced? It would have required a change in the slab-pull force that is the primary driver of plate tectonics, suggesting an increase in the amount of subduction in the Tethys Ocean that formerly lay between India and Eurasia, probably at two, now hidden destructive plate margins.

A group of geoscientists from Canada, the US and Pakistan has documented that collision in terms of the record of metamorphism experienced beneath the Himalaya as slab after slab of once near-surface rocks were driven beneath the rising orogen (Soret, M. et al. 2021. How Himalayan collision stems from subduction. Geology, v. 49, p. 894-898; DOI: 10.1130/G48803.1). The Western Himalaya has trapped a deformed and tilted magmatic rock sequence of an island arc – the Kohistan Arc – between  the Eurasian plate and a zone of crustal thickening and shortening that was thrust southward over the ancient metamorphic basement of India itself. That crust was mantled by a variety of younger sediments deposited on the Tethyan continental shelf of the northern Indian plate which became involved in the process of crustal thickening. The Kohistan Arc probably formed above one of the destructive margins that consumed the oceanic lithosphere of the now vanished Tethys Ocean. Two distinct types of rock make up the slabs stacked-up by thrusting.

The uppermost, which also forms the highest part of the Western Himalaya in the form of Nanga Parbat (at 8,126 metres the world’s ninth highest mountain) comprises rocks thought to represent Tethyan oceanic lithosphere subducted perhaps at the second destructive margin. Their mineral assemblages, especially those of eclogites, indicate that they have been metamorphosed under pressures corresponding to depths of up to 100 km, but at low temperatures along a geothermal gradient of about 7°C km-1, i.e. in a low heat-flow environment. These ultra-high pressure (UHP) metamorphic rocks formed at the start of the India-Eurasia collision. The sequence of sedimentary slabs now overridden by the UHP slab were metamorphosed at around the same time, but under very different conditions. Their burial reached only about 35 km – the normal thickness of the continental crust – and a temperature of about 600°C on a 30°C km-1 geothermal gradient. Detailed mineralogy of the UHP slab reveals that as it was driven over the metasediments it evolved to the same geothermal conditions.

Matthew Soret and his colleagues explain how this marked metamorphic duality may have arisen in rocks that are now part of the same huge thrust complex. Their results are consistent with slicing together of oceanic lithosphere in a subduction zone to form a tectonic wedge of UHP mineral assemblages at the same time as continental shelf sediments were metamorphosed under more normal geothermal conditions. This was happening just as India came into contact with Eurasia. When crustal thickening began in earnest through the inter-slicing of the two assemblages, pressure on the UHP rocks fell rapidly as a result of their being thrust over the dominantly metasedimentary shelf sequence. It also moved into a zone of normal heat flow, first heating up equally quickly and then following a path of decreasing pressure and temperature as erosion pared away the newly thickened crust. Both assemblages now became part of the same metamorphic regime. In this way a subduction system evolved to become incorporated in an orogenic zone as two continents collided; a complex process that finds parallels in other orogens such as the Alps.

The subduction pulley: a new feature of plate tectonics

Geological map of part of the Italian Alps. The Sesia-Lanzo Zone is 6 in the Key: a – highly deformed gneisses; b – metasedimentary schists with granite intrusions; c – mafic rocks; d – mixed mantle and crystalline basement rocks. (Credit: M. Assanelli, Universita degli Studi di Milano)

To a first approximation, as they say, the basis of plate tectonics is that the lithosphere is divided up into discrete, rigid plates that are bounded by lines of divergent, convergent and sideways relative motions: constructive, destructive and conservative plate margins. These are characterised by zones of earthquakes whose senses of motion roughly correspond to the nature of each boundary: normal, reverse and strike-slip, respectively. The seismicity is mainly confined to the lithosphere in the cases of constructive and conservative boundaries (i.e. shallow) but extends as deep as 700 km into the mantle at destructive margins, thereby defining the subduction of lithosphere that remains cool enough to retain its rigidity. Although the definition assumes that there is no deformation within plates, in practice that does occur for a wide variety of reasons in the form of intra-plate seismicity, mainly within continental lithosphere. Oceanic plate interiors are much stronger and largely ‘follow the rules’; they are generally seismically quiet.

One important feature of plate tectonics is the creation of new subduction zones when an earlier one eventually ceases to function. Where these form in an oceanic setting volcanism in the overriding plate creates island arcs. They create precursors of new continental crust because the density of magmas forming the new lithosphere confers sufficient buoyancy for them to be more difficult to subduct. Eventually island arcs become accreted onto continental margins through subduction of the intervening oceanic lithosphere. Joining them in such ‘docking’ are microcontinents, small fragments spalled from much older continents because of the formation of new constructive plate margins within them. It might seem that arcs and microcontinents behave like passive rafts to form the complex assemblages of terranes that characterise continental mountain belts, such as those of western North America, the Himalaya and the Alps. Yet evidence has emerged that such docking is much more complicated (Gün, E. et al. 2021. Pre-collisional extension of microcontinental terranes by a subduction pulleyNature Geoscience, v. 14, online publication; DOI: 10.1038/s41561-021-00746-9).

Erkan Gün and colleagues from the University of Toronto and Istanbul Technical University examined one of the terranes in the Italian Alps – the Sesia-Lanzo Zone (SLZ) – thought to have been a late-Carboniferous microcontinental fragment in the ocean that once separated Africa from Europe. When it accreted the SLZ was forced downwards to depths of up to 70 km and then popped up in the latter stages of the Alpine orogeny. It is now a high-pressure, low-temperature metamorphic complex, having reached eclogite facies during its evolution. Yet its original components, including granites that contain the high-pressure mineral jadeite instead of feldspar, are still recognisable. Decades of geological mapping have revealed that the SLZ sequence shows signs of large-scale extensional tectonics. Clearly that cannot have occurred after its incorporation into southern Europe, and must therefore have taken place prior to its docking. Similar features are present within the accreted microcontinental and island-arc terranes of Eastern Anatolia in Turkey. In fact, most large orogenic belts comprise hosts of accreted terranes that have been amalgamated into older continents.

An ‘engineering’ simplification of the subduction pulley. Different elements represent slab weight (slab pull force) transmitted through a pulley at the trench to a weak microcontinent and a strong oceanic lithosphere. (Credit: Gün et al., Fig. 4)

Lithospheric extension associated with convergent plate margins has been deduced widely in the form of back-arc basins. But these form in the plate being underidden by a subduction zone. Extension of the SLZ, however, must have taken place in the plate destined to be subducted. Gün et al. modelled the forces, lithospheric structure, deformation and tectonic consequences that may have operated to form the SLZ, for a variety of microcontinent sizes. The pull exerted by the subduction of oceanic lithosphere (slab pull) would exert extensional forces on the lithosphere as it approached the destructive plate boundary. Oceanic lithosphere is very strong and would remain intact, simply transmitting slab-pull force to the weaker continental lithosphere, which ultimately would be extended. This is what the authors call a subduction ‘pulley’ system. At some stage the microcontinent fails mechanically, part of it being detached to continue with the now broken slab down the subduction zone. The rest would become a terrane accreted to the overriding plate. Subduction at this site would stop because the linkage to the plate has broken. It may continue by being transferred to a new destructive margin ‘behind’ the accreted microcontinent. This would allow other weak continental and island-arc ‘passengers’ further out on the oceanic plate eventually to undergo much the same process.

The observed complexity of tectonic terranes in other vast assemblies of them, such as the northern Pacific coast of North America and in many more ancient orogenic belts, is probably as much a result of extension before accretion as the compressional deformation suffered afterwards. The theoretical work by Erkan Gün and colleagues will surely spur tectonicians to re-evaluate earlier models of orogenesis.

Note: Figure 2 in the paper by Gün et al. shows how the width (perpendicular to the subduction zone) affects the outcomes of the subduction pulley. View an animation of a subduction pulley

Diamonds and the deep carbon cycle

When considering the fate of the element carbon and CO2, together with all their climatic connotations, it is easy to forget that they may end up back in the Earth’s mantle from which they once escaped to the surface. In fact all geochemical cycles involve rock, so that elements may find their way into the deep Earth through subduction, and they could eventually come out again: the ‘logic’ of plate tectonics. Teasing out the various routes by which carbon might get to mantle is not so easily achieved. Yet one of the ways it escapes is through the strange magma that once produced kimberlite intrusions, in the form of pure-carbon crystals of diamond that kimberlites contain. A variety of petrological and geochemical techniques, some hinging on other minerals that occur as inclusions, has allowed mineralogists to figure out that diamonds may form at depths greater than about 150 km. Most diamonds of gem quality formed in unusually thick lithosphere beneath the stable, and relatively cool blocks of ancient continental crust known as cratons, which extends to about 250 km. But there are a few that reflect formation depths as great as 800 km that span two major discontinuities in the mantle (at 410 and 660 km depth). These transition zones are marked by sudden changes in seismic speed due to pressure-induced transformations in the structure and density of the main mantle mineral, olivine.

Diamond crystal containing a garnet and other inclusions (Credit: Stephen Richardson, University of Cape Town, South Africa)

Carbon-rich rocks that may be subducted are not restricted to limestones and carbon-rich mudstones. Far greater in mass are the basalts of oceanic crust. Not especially rich in carbon when they crystallised as igneous rocks, their progress away from oceanic spreading centres exposes them to infiltration by ocean water. Once heated, aqueous fluids cause basalts to be hydrothermally altered. Anhydrous feldspars, pyroxenes and olivines react with the fluids to break down to hydrated-silicate clays and dissolved metals. Dissolved carbon dioxide combines with released calcium and magnesium to form pervasive carbonate minerals, often occupying networks of veins. So there has been considerable dispute as to whether subducted sediments or igneous rocks of the oceanic crust are the main source of diamonds. Diamonds with gem potential form only a small proportion of recovered diamonds. Most are only saleable for industrial uses as the ultimate natural abrasive and so are cheaply available for research. This now centres on the isotopic chemistry of carbon and nitrogen in the diamonds themselves and the various depth-indicating silicate minerals that occur in them as minute inclusions, most useful being various types of garnet.

The depletion of diamonds in ‘heavy’ 13C once seemed to match that of carbonaceous shales and the carbonates in fossil shells, but recent data from carbonates in oceanic basalts reveals similar carbon, giving three possibilities. Yet, when their nitrogen-isotope characteristics are taken into account, even diamonds that formed at lithospheric depths do not support a sedimentary source (Regier, M.E. et al. 2020. The lithospheric-to-lower-mantle carbon cycle recorded in superdeep diamonds. Nature, v. 585, p. 234–238; DOI: 10.1038/s41586-020-2676-z). That leaves secondary carbonates in subducted oceanic basalts as the most likely option, the nitrogen isotopes more reminiscent of clays formed from igneous minerals by hydrothermal processes than those created by weathering and sedimentary deposition. However, diamonds with the deepest origins – below the 660 km mantle transition zone – suggest yet another possibility, from the oxygen isotopes of their inclusions combined with those of C and N in the diamonds. All three have tightly constrained values that most resemble those from pristine mantle that has had no interaction with crustal rocks. At such depths, unaltered mantle probably contains carbon in the form of metal alloys and carbides. Regier and colleagues suggest that subducted slabs reaching this environment – the lower mantle – may release watery fluids that mobilise carbon from such alloys to form diamonds. So, I suppose, such ultra-deep diamonds may be formed from the original stellar stuff that accreted to form the Earth and never since saw the ‘light of day’.

Recycling of continental crust through time

Because continental crust is so light – an average density of 2700 kg m-3 compared with the mantles’ value of 3300 – it has been widely believed that continents cannot be subducted en masse. Yet it is conceivable that sial can be ‘shaved’ from below during subduction and from above by erosion and added to subductable sediment on the ocean floor. Certainly, there is overwhelming evidence for the net growth of continents through time and plenty for periods of increased and dwindling growth in the past. In some ancient orogens there are substantial slabs of continental composition whose mineralogy bears witness to ultra-high pressure metamorphism at depths greater than that of the base of continents. These slabs had been caught-up in subduction but never reached sufficiently high density to be retained by the mantle; they eventually ‘bobbed up’ again. On the other hand, if early continents were less silica rich through incorporation of substantial proportions of rock with basaltic composition parts of them could founder if subjected to high-pressure, low-temperature metamorphism. But not all crustal recycling to the mantle is through subduction. Some abnormally highly elevated parts of the continents that rose quickly in geological terms, such as the Tibetan Plateau, may have formed by lower crustal slabs becoming detached or delaminated from their base. Again modelling can help assess the past magnitude of continental recycling (Chowdhury, P. et al. 2017. Emergence of silicic continents as the lower crust peels off on a hot plate-tectonic Earth. Nature Geoscience, v. 10, p. 698-703; DOI: 10.1038/NGEO3010).

Various lines of evidence suggest that between 65 to 70% of the present continental volume existed by 3 billion years ago, yet that does not manifest itself in the rock record; perhaps a sign that some has returned to the mantle. It is also widely suggested that plate tectonics in the modern style began at about that time. Pryadarshi Chowdhury and colleagues simulate what may happen at depth in continent-continent collision zones – the classic site of orogenies –at different times in the past. Under the hotter conditions in the early Archaean mantle delamination would have been more likely than it has been during the Phanerozoic; i.e. the peeling off and sinking of the denser, more mafic lower crust and the attached upper mantle. The authors show that increased mantle temperature further back in time increases the likelihood and extent of such delamination. It also encourages partial melting of the descending continental material so creating rising bodies of more silicic magma that add to the remaining continent at the surface. Together with the lower crust’s attachment of to a mantle slab, this ensures that the peeled off material is able to descend under its own load. Once below a depth of 250 km felsic rocks are doomed to further descent. Waning of radiogenic mantle heat production encourages descending slabs to fail and break from the connection with lithosphere at higher levels so that a smaller proportion of the lower crust becomes detached and recycled. This evolution suggests that less and less continental crust is recycled with time. This broadly fits with current geochemical ideas based on the record of radiogenic Nd-, Sr- and Pb-isotopes in rocks ranging from early Archaean to Phanerozoic age.

Plate tectonic graveyard

Where do old plates go to die? For the most part, down subduction zones to mix with their original source, the mantle. Earth-Pages has covered evidence for quite a few of the dead plates, which emerges from a geophysical technique known as seismic tomography – analogous to X-ray or magnetic resonance scans of the whole human body. For 20 years geophysicists have been analysing seismograms from many stations across the globe for every digitally recorded earthquake, i.e. virtually all of those since the 1970s. This form of depth sounding goes far beyond early deep-Earth seismometry that discovered the inner and outer core, various transition zones in the mantle and measured the average variation with depth of mantle properties. Tomography relies on complex models of the paths taken by seismic body waves and very powerful computing to assess variations in the speed of P- and S-waves as they travelled through the Earth: the more rigid/cool the mantle is the faster waves travel through it and vice versa. The result is images of deep structure in 2-D slices, but the quality of such sections depends, ironically, on plate tectonics. Most earthquakes occur at plate boundaries. Such linearly distributed, one-dimensional sources inevitably leave the bulk of the mantle as a blur. Around 20 different methodologies have been developed by the many teams working on seismic tomography. So sometimes conflicting images of the deep Earth have been produced.

Results of seismic tomography across Central America showing anomalously fast (in blue) P- (top) and S-wave (bottom) speeds in map view at a fixed mantle depth (1290 km, left) and as vertical sections (right). The blue zones at right are interpreted to show a steeply dipping slab that represents subduction of the eastern Pacific Cocos plate since about 175 Ma ago (credit: van der Meer, D.G et al. ‘Atlas of the Underworld)

The technique has come of age now that superfast computing and use of multiple models have begun to resolve some of tomography’s early problems. The latest outcome is astonishing: ‘The Atlas of the Underworld’ catalogues 94 2-D sections from surface to the core-mantle boundary each of which spans 40° or arc – about a ninth of the Earth’s circumference (see: van der Meer, D.G., van Hinsbergen, D.J.J., and Spakman, W., 2017, Atlas of the Underworld: slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity, Tectonophysics online; doi.org/10.1016/j.tecto.2017.10.004). Specifically, the Atlas locates remnants of relatively cold slabs in the mantle that are suspected to be remnants of former subduction zones, or those that connect to active subduction. The upper parts of active slabs are revealed by the earthquakes generate along them. At deeper levels they are too ductile to have seismicity, so what form they take has long been a mystery. Once subduction stops, so do the telltale earthquakes and the slabs ‘disappear’.

The slabs covered by the ‘Atlas’ only go back as far as the end of the Permian, when the current round of plate tectonics began as Pangaea started to break-up. It takes 250 Ma for slabs to reach the base of the mantle and beyond that time they will have heated up and begun to be mixed into the lower mantle and invisible. Nevertheless, the rich resource allows models of vanished Mesozoic to Recent plates and the tectonics in which they participated, based on geological information, to be evaluated and enriched. Just as important, the project opens up the possibility of finding out how the mantle ‘worked’ since Pangaea broke up, in 3-D; a key to more than plate tectonics, including the mantle’s chemical heterogeneity. Already it has been used to estimate changes in the total length of subduction zones since 250 Ma ago, and thus arc volcanism and CO2 emissions, which correlates with estimates of past atmospheric CO2 levels, climate and even sea levels.

See also:  Voosen, P. 2016. ‘Atlas of the Underworld’ reveals oceans and mountains lost to Earth’s history. Science; doi:10.1126/science.aal0411.

Lee, H. 2017. The Earth’s interior is teeming with dead plates. Ars Technica UK, 18 October 2017.

So, when did plate tectonics start up?

Tiny, 4.4 billion year old zircon grains extracted from much younger sandstones in Western Australia are the oldest known relics of the Earth system. But they don’t say much about early tectonic processes. For that, substantial exposures of rock are needed, of which the undisputedly oldest are the Acasta gneisses 300 km north of Yellowknife in Canada’s North West Territories, which have an age of slightly more than 4 Ga. The ‘world’s oldest rock’ has been something of a grail for geologists and isotope geochemists who have combed the ancient Archaean cratons for 5 decades. But since the discovery of metasediments with an age of 3.8 Ga in West Greenland during the 1970s they haven’t made much headway into the huge time gap between Earth’s accretion at 4.54 Ga and the oldest known rocks (the Hadean Eon).

The Deccan Traps shown as dark purple spot on ...
Continental cratons (orange) where very-old rocks are likely to lurk. (credit: Wikipedia)

There have been more vibrant research themes about the Archaean Earth system, specifically the issue of when our planet settled into its modern plate tectonic phase A sprinkling of work on reconstructing the deep structural framework of Archaean relics has convinced some that opposed motion of rigid, brittle plates was responsible for their geological architecture, whereas others have claimed signs of a more plastic and chaotic kind of deformation of the outer Earth. More effort has been devoted to using the geochemistry of all the dominant rocks found in the ancient cratons, seeking similarities with and differences from those of more recent vintage. There can be little doubt that the earliest processes did form crust whose density prevented or delayed it from being absorbed into the mantle. Even the 4.4 Ga zircons probably crystallized from magma that was felsic in composition. Once trapped by buoyancy at the surface and subsequently wrapped around by similarly low density materials continental crust formed as a more or less permanent rider on the Earth’s deeper dynamics. But did it all form by the same kinds of process that we know to be operating today?

Plate tectonics involves the perpetual creation of rigid slabs of basalt-capped oceanic lithosphere at oceanic rift systems and their motions and interactions, including those with continental crust. Ocean floor cools as it ages and becomes hydrated by seawater that enters it. The bulk of it is destined eventually to oppose, head-to-head, the motions of other such plates and to deform in some way. The main driving force for global tectonics begins when an old, cold plate does deform, breaks, bends and drives downwards. Increasing pressure on its cold, wet basaltic top transforms it into a denser form: from a wet basaltic mineralogy (feldspar+pyroxene+amphibole) to one consisting of anhydrous pyroxene and garnet (eclogite) from which watery fluid is expelled upwards. Eclogite’s density exceeds that of mantle peridotite and compels the whole slab of oceanic lithosphere to sink or subduct into the mantle, dragging the younger parts with it. This gravity-induced ‘slab pull’ sustains the sum total of all tectonic motion. The water rising from it induces the wedge of upper mantle above to melt partially, the resulting magma evolves to produce new felsic crust in island arcs whose destiny is to be plastered on to and enlarge older continental masses.

Relics of eclogites and other high-pressure, low-temperature versions of hydrated basalts incorporated into continents bear direct and unchallengeable witness to plate tectonics having operated back to about 800 Ma ago. Before that, evidence for plate tectonics is circumstantial and in need of special pleading. Adversarial to-ing and fro-ing seems to be perpetual, between geoscientists who see no reason to doubt that Earth has always behaved in this general fashion and others who see room for very different scenarios in the distant past. The non-Huttonian tendency suggests an early, more ductile phase when greater radioactive heat production in the mantle produced oceanic crust so fast that when it interacted with other slabs it was hot enough to resist metamorphic densification wherever it was forced down. Faster production of magma by the mantle without slab-pull could have produced a variety of ‘recycling’ turnover mechanisms that were not plate-tectonic.

One thing that geochemists have discovered is that the composition of Archaean continental crust is very different from that produced in later times. In 1985 Ross Taylor and Scott McLennan, then of the Australian National University, hit on the idea of using shales of different ages as proxies for the preceding continental crust from which they had been derived by long erosion. Archaean and younger shales differed in such a way that suggests that after 2.5 Ga (the end of the Archaean) vast amounts of feldspar were extracted from the continent-forming magmas. This left the later Precambrian and Phanerozoic upper crust depleted in the rare-earth element europium, which ended up in a mafic, feldspar-rich lower crust. On the other hand, no such mass fractionation had left such a signature before 2.5 Ga. Another ANU geochemist, now at the University of Maryland, Roberta Rudnick has subsequently carried this approach further, culminating in a recent paper (Tang, M., Chen, K and Rudnick, R.L. 2016. Archean upper crust transition from mafic to felsic marks the onset of plate tectonics. Science, v. 351, p. 372-375). This uses nickel, chromium and zinc concentrations in ancient igneous and sedimentary rocks to track the contribution of magnesium (the ‘ma’ in ‘mafic’) to the early continents. The authors found that between 3.0 to 2.5 Ga continental additions shifted from a dominant more mafic composition to one similar to that of later times by the end of the Archaean. Moreover, this accompanied a fivefold increase in the pace of continental growth. Such a spurt has long been suspected and widely suggested to mark to start of true plate tectonics: but an hypothesis bereft of evidence.

A better clue, in my opinion, came 30 years ago from a study of the geochemistry of actual crustal rocks that formed before and after 2.5 Ga (Martin, H. 1986. Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas. Geology, v. 14, p. 753-756). Martin showed that plutonic Archaean and post-Archaean felsic rocks of the continental crust lie in distinctly different fields on plots of their rare-earth element (REE) abundances. Archaean felsic plutonic rocks show a distinct trend of enrichment in light REE relative to heavy REE as measures of the degree of partial melting decreases, whereas the younger crustal rocks show almost constant, low values of heavy REE/light REE whatever the degree of melting. The conclusion he reached was that while in the post Archaean the source was consistent with modern subduction processes – i.e. partial melting of hydrated peridotite in the mantle wedge above subduction zones – but during the Archaean the source was hydrated, garnet-bearing amphibolite of basaltic composition, in the descending slab of subducted oceanic crust. Together with Taylor and McLennan’s lack of evidence for any fractional crystallization in Archaean continental growth, in contrast to that implicated in Post-Archaean times.

The geochemistry forces geologists to accept that a fundamental change took place in the generation and speed of continental growth at the end of the Archaean, marking a shift from a dominance of melting of oceanic, mafic crust to one where the upper mantle was the main source of felsic, low-density magmas. Yet, no matter how much we might speculate on indirect evidence, whether or not subduction, slab-pull and therefore plate tectonics dominated the Archaean remains an open question.

More on continental growth and plate tectonics

Seismic menace of the Sumatra plate boundary

More than a decade after the 26 December 2004 Great Aceh Earthquake and the Indian Ocean tsunamis that devastating experience and four more lesser seismic events (> 7.8 Magnitude) have show a stepwise shift in activity to the SE along the Sumatran plate boundary. It seems that stresses along the huge thrust system associated with subduction of the Indo-Australian Plate that had built up over 200 years of little seismicity are becoming unlocked from sector to sector along the Sumatran coast. Areas further to the SE are therefore at risk from both major earthquakes and tsunamis. A seismic warning system now operates in the Indian Ocean, but the effectiveness of communications to potential victims has been questioned since its installation. However, increasing sophistication of geophysical data and modelling allows likely zones at high risk to be assessed.

Recent Great Earthquakes in different segments of the Sumatra plate margin (credit: Tectonics Observatory, California Institute of Technology http://www.tectonics.caltech.edu/outreach/highlights/sumatra/why.html
Recent Great Earthquakes in different segments of the Sumatra plate margin (credit: Tectonics Observatory, California Institute of Technology http://www.tectonics.caltech.edu/outreach/highlights/sumatra/why.html

One segment is known to have experienced giant earthquakes in 1797 and 1833 but none since then. What is known as the Mentawai seismic gap lies between two other segments in which large earthquakes have occurred in the 21st century: it is feared that gap will eventually be filled by another devastating event. Geophysicists from the Institut de Physique du Globe de Paris and Nanyang Technological University in Singapore have published a high-resolution seismic reflection survey showing the subduction zone beneath the Mentawai seismic gap (Kuncoro, A.K. et al. 2015. Tsunamigenic potential due to frontal rupturing in the Sumatra locked zone. Earth and Planetary Science Letters, v. 432, p. 311-322). It shows that that the upper part of the zone, the accretionary wedge, is laced with small thrust-bounded ‘pop-ups’. The base of the accretionary wedge shows a series of small seaward thrusts above the subduction surface itself forming ‘piggyback’ or duplex structures.

Seismic reflection profile across part of the Sumatra plate boundary, showing structures produced by past seismicity. (credit: Kuncoro et al. 2015, Figure 3b)
Seismic reflection profile across part of the Sumatra plate boundary, showing structures produced by past seismicity. (credit: Kuncoro et al. 2015, Figure 3b)

The authors model the mechanisms that probably produced these intricate structures. This shows that the inactive parts of the plate margin have probably locked in stresses equivalent to of the order of 10 m of horizontal displacement formed by the average 5 to 6 cm of annual subduction of the Indo-Australian Plate over the two centuries since the last major earthquakes. Reactivation of the local structures by release of this strain would distribute it by horizontal movements of between 5.5 to 9.2 m and related 2 to 6.6 m vertical displacement in the pop-ups. That may suddenly push up the seafloor substantially during a major earthquake, thereby producing tsunamis. Whether or not this is a special feature of the Sumatra plate boundary that makes it unusually prone to tsunami production is not certain: such highly resolving seismic profiles need to be conducted over all major subduction zones to resolve that issue. What does emerge from the study is that a repeat of the 2004 Indian Ocean tsunamis is a distinct possibility, sooner rather than later.

Subduction and the water cycle

Note: Earth-Pages will be closing as of early July, but will continue in another form at Earth-logs

For many geoscientists and lay people the water cycle is considered to be part of the Earth’s surface system. That is, the cycle of evapotranspiration, precipitation and infiltration involving atmosphere, oceans, cryosphere, terrestrial hydrology and groundwater. Yet it links to the mantle through subduction of hydrated oceanic lithosphere and volcanism. The rate at which water vapour re-enters the surface part of the water cycle through volcanoes is reasonably well understood, but the same cannot be said about ‘recharge’ of the mantle through subduction.

Water cycle http://ga.water.usgs.gov/edu/water...
The water cycle as visualised by the US Geological Survey (credit: Wikipedia)

Subducted oceanic crust is old, cold and wet: fundamentals of plate theory. The slab-pull that largely drives plate tectonics results from phase transitions in oceanic crust that are part and parcel of low-temperature – high-pressure metamorphism. They involve the growth of the anhydrous minerals garnet and high-pressure pyroxene that constitute eclogite, the dense form taken by basalt that causes the density of subducted lithosphere to exceed that of mantle peridotite and so to sink. This transformation drives water out of subducted lithosphere into the mantle wedge overlying a subduction zone, where it encourages partial melting to produce volatile-rich andesitic basalt magma – the primary magma of island- and continental-arc igneous activity. Thus, most water that does reenter the mantle probably resides in the ultramafic lithospheric mantle in the form of hydrated olivine, i.e. the mineral serpentine, and that is hard to judge.

Water probably gets into the mantle lithosphere when the lithosphere bends to begin its descent. That is believed to involve faults – cold lithosphere is brittle – down which water can diffuse to hydrate ultramafic rocks. So the amount of water probably depends on the number of such bend-related faults. A way of assessing the degree of such faulting and thus the proportion of serpentinite is analysis of seismic records from subduction zones. This has been done from earthquake records from the West Pacific subduction zone descending beneath northern Japan (Garth, T. & Rietbrock, A. 2014. Order of magnitude increase in subducted H2O due to hydrated normal faults within the Wadati-Bennioff zone. Geology, on-line publication doi:10.1130/G34730.1). The results suggest that between 17 to 31% of the subducted mantle there has been serpentinised.

In a million years each kilometre along the length of this subduction zone would therefore transfer between 170 to 318 billion tonnes of water into the mantle; an estimate more than ten times previous estimates. The authors observe that at such a rate a subduction zone equivalent to the existing, 3400 km long Kuril and Izu-Bonin arcs that affect Japan would have transferred sufficient water to fill the present world oceans 3.5 times over the history of the Earth. Had the entire rate of modern subduction along a length of 55 thousand kilometres been maintained over 4.5 billion years, the world’s oceans would have been recycled through the mantle once every 80 million years. To put that in perspective, since the Cretaceous Chalk of southern England began to be deposited, the entire mass of ocean water has been renewed. Moreover, subduction has probably slowed considerably through time, so the transfer of water would have been at a greater pace in the more distant past.

Enhanced by Zemanta

How the great Tohoku-Sendai earthquake and tsunami happened

The great Tohoku earthquake (moment magnitude 9.0) of 11 March 2011 beneath the Pacific Ocean off the largest Japanese island of Honshu resulted in the devastating tsunami that tore many kilometres inland along its northern coast line and affected the entire Pacific Basin (see NOAA animation of the tsunami’s propagation) .

English: Sendai Rinkai Railway locomotive(SD55...
Railway locomotive thrown aside by the 11 March 2011 Tsunami in Japan. (credit: Wikipedia)

This article can now be read in full at Earth-logs in the Geohazards archive for 2017

Brittle-ductile deformation in subduction zones

Almenning, Norway. The red-brown mineral is ga...
Eclogite: the red-brown mineral is garnet, omphacite is green and there is some white quartz.(credit: Kevin Walsh via Wikipedia)

The ultra-dense form of basalt, eclogite made from mainly garnet and a strange high-pressure, low-temperature pyroxene (omphacite) that forms from plagioclase and some of the basalt’s ferromagnesian minerals, is possibly the most important rock there is. Without the basalt to eclogite transition that takes place when ocean-floor is subducted the density of the lithosphere would be insufficient to pull more ocean floor to destruction and maintain the planetary circulation otherwise known as plate tectonics. Since the transition involves the formation of anhydrous eclogite from old, cold and wet basalt water is driven upwards into the mantle wedge that lies over subduction zones. The encourages partial melting which creates andesite magmas and island arcs, the ultimate source of the Earth’s continental crust.

Despite being cold and rigid, subducted oceanic lithosphere somehow manages to be moved en masse, showing its track by earthquakes down to almost 700 km below the Earth’s surface.  A major ophiolite in the Western Alps on the Franco-Italian border escaped complete loss to the mantle by rebounding upwards after being subducted and metamorphosed under high-P, Low-T condition when the Alps began to form. So the basaltic crustal unit is eclogite and that preserves a petrographic  record of what actually happened as it descended (Angiboust, S. et al. 2012. Eclogite breccia in a subducted ophiolite: A record of intermediate depth earthquakes? Geology, v.  40, p. 707-710). The French geologists found breccias consisting of gabbroic eclogite blocks set in a matrix of serpentinite and talc. The blocks themselves are breccias too, with clasts of eclogite mylonite set in fine-grained lawsonite-bearing eclogite. The relationships in the breccias point to possibly earthquake-related processes, grinding and fracturing basalt as it was metamorphosed: an essentially brittle process, yet the shearing that forms mylonites does seem reminiscent of ductile deformation too.

The deformation seems to have been at the middle level of oceanic crust where oceanic basalt lavas formed above cumulate gabbro, their plutonic equivalents. Yet much deformation was also at the gabbro-serpentinite or crust-mantle boundary, where water loss from serpentine may have helped lubricate some of the processes. Clearly the Monviso ophiolite will soon become a place to visit for geophysicists as well as metamorphic petrologists.

The ultra-deep carbon cycle

A scattering of "brilliant" cut diam...
Image via Wikipedia

The presence of diamonds in the strange, potassium-rich, mafic to ultramafic igneous rocks known as kimberlites clearly demonstrates that there is carbon in the mantle, but it could have come from either biogenic carbon having moved down subduction zones or the original meteoritic matter that accreted to form the Earth. Both are distinct possibilities for which evidence can only be found within diamonds themselves as inclusions. There is a steady flow of publications focussed on diamond inclusions subsidised to some extent by companies that mine them (see Plate tectonics monitored by diamonds in EPN, 2 August 2011). The latest centres on the original source rocks of kimberlites and the depths that they reached (Walter, M.J. and 8 others 2011. Deep mantle cycling of oceanic crust: evidence from diamonds and their mineral inclusions. Science, v. 334, p. 54-57). The British, Brazilian and US team analysed inclusions in diamonds from Brazil, finding assemblages that are consistent with original minerals having formed below the 660 km upper- to lower-mantle seismic boundary and then adjusting to decreasing pressure as the kimberlite’s precursor rose to melt at shallower levels. The minerals – various forms of perovskite stable at deep-mantle pressures – from which the intricate composites of several lower-pressure phases exsolved suggest the diamonds originated around 1000 km below the surface; far deeper than did more common diamonds. Moreover, their geochemistry suggests that the inclusions formed from deeply subducted basalts of former oceanic crust.

Previous work on the carbon isotopes in ‘super-deep’ diamonds seemed to rule out a biogenic origin for the carbon, suggesting that surface carbon does not survive subduction into the lower mantle. In this case, however, the diamonds are made of carbon strongly enriched in light 12C relative to 13C, with δ13C values of around -20 ‰ (per thousand), which is far lower than that found in mantle peridotite and may have been subducted organic carbon. If that proves to be the case it extends the global carbon cycle far deeper than had been imagined, even by the most enthusiastic supporters of the Gaia hypothesis.

Plate tectonics monitored by diamonds

eclogite
Norwegian Eclogite. Image by kevinzim via Flickr

For more than 30 years a debate has raged about the antiquity of plate tectonics: some claim it has always operated since the Earth first acquired a rigid carapace not long after a molten state following formation of the Moon; others look to the earliest occurrences of island-arc volcanism, oceanic crust thrust onto continents as ophiolite complexes, and to high-pressure, low-temperature metamorphic rocks. The earliest evidence of this kind has been cited from as far apart in time as the oldest Archaean rocks of Greenland (3.9 Ga) and the Neoproterozoic (1 Ga to 542 Ma). A key feature produced by plate interactions that can be preserved are high-P, low-T rocks formed where old, cool oceanic lithosphere is pulled by its own increasing density into the mantle at subduction zones to form eclogites and blueschists. In the accessible crust, both rock types are unstable as well as rare and can be retrogressed to different metamorphic mineral assemblages by high-temperature events at lower pressures than those at which they formed. Relics dating back to the earliest subduction may be in the mantle, but that seems inaccessible. Yet, from time to time explosive magmatism from very deep sources brings mantle-depth materials to the surface in kimberlite pipes that are most commonly found in stabilised blocks of ancient continental crust or cratons. Again there is the problem of mineral stability when solids enter different physical conditions, but there is one mineral that preserves characteristics of its deep origins – diamond. Steven Shirer and Stephen Richardson of the Carnegie Institution of Washington and the University of Cape Town have shed light on early subduction by exploiting the relative ease of dating diamonds and their capacity for preserving other minerals captured within them (Shirey, S.B. & Richardson, S.H. 2011. Start of the Wilson cycle at 3 Ga shown by diamonds from the subcontinental mantle. Science, v. 333, p. 434-436). Their study used data from over four thousand silicate inclusions in previously dated large diamonds, made almost worthless as gemstones by their contaminants. It is these inclusions that are amenable to dating, principally by the Sm-Nd method. Adrift in the mantle high temperature would result in daughter isotopes diffusing from the minerals. Once locked within diamond that isotopic loss would be stopped by the strength of the diamond structure, so building up with time to yield an age of entrapment when sampled.  The collection spans five cratons in Australia, Africa, Asia and North America, and has an age spectrum from 1.0 to 3.5 Ga. Note that diamonds are not formed by subduction but grow as a result of reduction of carbonates or oxidation of methane in the mantle at depths between 125 to 175 km. In growing they may envelop fragments of their surroundings that formed by other processes.

A notable feature of the inclusions is that before 3.2 Ga only mantle peridotites (olivine and pyroxene) are trapped, whereas in diamonds younger than 3.0 Ga the inclusions are dominated by eclogite minerals (garnet and Na-, Al-rich omphacite pyroxenes). This dichotomy is paralleled by the rhenium and osmium isotope composition of sulfide mineral inclusions. To the authors these consistent features point to an absence of steep-angled subduction, characteristic of modern plate tectonics, from the Earth system before 3 Ga. But does that rule out plate tectonics in earlier times and cast doubt on structural and other evidence for it? Not entirely, because consumption of spreading oceanic lithosphere by the mantle can take place if basaltic rock is not converted to eclogite by high-P, low-T metamorphism when the consumed lithosphere is warmer than it generally is nowadays – this happens beneath a large stretch of the Central Andes where subduction is at a shallow angle. What Shirey and Richardson have conveyed is a sense that the dominant force of modern plate tectonics – slab-pull that is driven by increased density of eclogitised basalt – did not operate in the first 1.5 Ga of Earth history. Eclogite can also form, under the right physical conditions, when chunks of basaltic material (perhaps underplated magmatically to the base of continents) founder and fall into the mantle. The absence of eclogite inclusions seems also to rule out such delamination from the early Earth system. So whatever tectonic activity and mantle convection did take place upon and within the pre-3 Ga Earth it was probably simpler than modern geodynamics. The other matter is that the shift to dominant eclogite inclusions appears quite abrupt from the data, perhaps suggesting major upheavals around 3 Ga. The Archaean cratons do provide some evidence for a major transformation in the rate of growth of continental crust around 3 Ga; about 30-40 percent of modern continental material was generated in the following 500 Ma to reach a total of 60% of the current amount, the remaining 40% taking 2.5 Ga to form through modern plate tectonics

Atlantic subduction due soon!

Rio de Janeiro
Rio de Janeiro, a threatened city? Image by Alcindo Correa Filho via Flickr

Earthquake prediction has not had a good record, but it seems that vastly larger tectonic processes are now becoming the subject of risk analysis (Nikolaeva, K. et al. 2011. Numerical analysis of subduction initiation risk along the Atlantic American passive margins. Geology, v. 39, p. 463-466). The Swiss, Russian and Portuguese authors focus on the old (Jurassic ~170 Ma) and presumably cold oceanic lithosphere on the western flank of the Atlantic, against both the North and South American continents. Increased density with ageing imparts a potential downwards force, but that has to overcome resistance to plate failure at passive margins. The dominance of upper continental lithosphere by rheologically weak quartz tends to make it more likely to fail than more or less quartz-free oceanic lithosphere. So, if subduction at a passive continental margin is to take place, then where and when it begins depends on the nature of the abutting continental lithosphere. That on the Atlantic’s western flank varies a lot, ranging from 75-150 km thick. Consequently the temperature at the Moho, the junction between continental lithosphere and weaker asthenosphere, varies too. The loading by marginal sedimentation also plays a role, as do continent-wide forces associated with far-distant mountain ranges, such as the Western Cordillera and Andes, and the forces from opposed sea-floor spreading from the Juan de Fuca and East Pacific systems that affect the whole of western South America, most of Central America and the far NW of North America.

Analysing all pertinent forces acting along 9 lines of section through both North and South America, the authors’ focus fell on the relatively thin continental lithosphere of the Atlantic margin of South America. It is at its thinnest along the southernmost part of the margin adjacent to Brazil, where the Moho temperature reaches as high as 735°C: the weakest link in the American continental lithosphere, where there is seismicity and also indications of igneous activity. The modelling suggests that incipient deformation may begin off southern Brazil within 4 Ma to form a zone of overthrusting, eventually evolving towards failure of the ocean-continent interface and the start of proper subduction in the succeeding 20 Ma. Other stretches of the eastern Americas are deemed safe from subduction for considerably longer by virtue of their greater thickness, lower Moho temperatures and thus higher strength. It is an interesting situation because, insofar as I understand plate tectonics, extensional or compressional failure needed to generate plate boundaries must also propagate from the weak spots that first fail; plate boundaries are lines not points. If that does not happen, then the very strength of the overwhelming longer continent-ocean interface will surely prevent subduction at a single, albeit weak link.

Paper PDF at http://xa.yimg.com/kq/groups/13231164/1842350625/name/Geology-2011-Nikolaeva-463-6.pdf

Bouncing back from the deep

eclogite
Eclogite from Norway. Image by kevinzim via Flickr

Because the average density of the rocks making up the continental crust is about 2.7 t m-3 while that of the mantle is greater than 3.0 t m-3 it might seem as though continents cannot be subducted. Indeed, that was one of the first principles of plate tectonics, which would account for continental crust dating back to 4000 Ma, whereas there is no oceanic crust older than about 150 Ma. In the southern foothills of the Alps in Piemonte, Italy is a site which refutes the hypothesis in a stunning fashion. The minor ski resort of Monte Mucrone is backed by cliffs in what to all appearances is a common-or-garden granite: it even seems to contain phenocrysts of plagioclase feldspar. Microscopic examination of the megacrysts reveals them to be made up of a complex intergrowth between jadeite, a high-pressure sodic pyroxene, and quartz. This is exactly what should form if albite, the sodium-rich kind of plagioclase feldspar, if it descended to depths over 70 km below the surface, i.e. to mantle depths.

Monte Mucrone proves that continental materials can be subducted, but also reveals that these granites popped back up again when the forces of subduction were relieved at the end of the Alpine orogeny. Other examples have since turned up, but few so spectacular as continental rocks from Switzerland (Herwartz, D. et al. 2011. Tracing two orogenic cycles in one eclogite sample by Lu-Hf garnet chronometry. Nature Geoscience, v. 4, p. 178-183). The Adula nappe of the Swiss Lepontine Alps consists of granitoid gneisses and metasediments of continental affinities, associated with mafic and ultramafic metamorphic rocks. The mafic rocks include eclogites typical of high-pressure, low-temperature metamorphism characteristic of subduction. Their minerals record formation temperatures around 680°C at a depth of more than more than 80 km. Eclogites are beautiful green and red rocks containing high-pressure omphacite pyroxene and pyrope garnet. Garnets generally contain abundant rare-earth elements especially those with the highest atomic numbers. One of these is lutetium (Lu) that has a radioactive isotope 176Lu with a half-life of 3.78×1010 years to yield a daughter isotope of hafnium 176Hf; garnets can be dated using this method. Garnets are frequently zoned, and the Adula eclogites clearly show several generations of zonation. Zoning can form as metamorphic conditions change, so in itself is not unusual, but dating different generations is. The German team from the Universities of Bonn, Cologne and Münster found that the garnets defined two distinct isochrons, one of Variscan age of just over 330 Ma, the other Alpine around 38 Ma. Clearly the pre-Variscan crust (probably once part of the African continent) had been subducted twice but had wrested itself clear of the mantle’s clutches on both occasions, each time remaining more or less intact. One idea that stems from this coincidence is that the Variscan mountain belt that formed at the earlier subduction zone subsequently split at its high P – low T core, so that the eclogites lay at a new continental margin and could suffer the same extreme compression when new subduction began there.

It also turns out that the region in which  Monte Mucrone lies, the Sesia zone of the Western Alps, also records a double whammy of continental subduction, but a repetition that occurred during the early events of the  Alpine orogeny (Rubatto, D. et al. 2011. Yo-yo subduction recorded by accessory minerals in the Italian Western Alps. Nature Geoscience, v. 4, p. 338-342). The team of Australian, Swiss and Italian geologists focused on the P-T record preserved in zoned garnets, allanites and zircons and evidence for two generation of white micas in eclogites and blueschists. Backed by U-Pb dating of zircon and allanite zones, the authots uncovered two episodes of deep subduction separated by period of rapid exhumation over the period between 79 to 65 Ma ago. The double subduction took place while the African plate converged obliquely with Eurasia; a strike-slip configuration that probably resulted in large-scale switches from compression to extension.

See also: Bruekner, H.K. 2011. Double-dunk tectonics. Nature Geoscience, v. 4, p. 136-138