Modelling the effects of Hadean impacts

The Hadean Eon (~4.6 to 4.0 Ga) is short on rocks that represent it. In fact geologists only know of a single 20 km2 outcrop within that age span (~4.3 Ga): the Nuvvuagittuq Greenstone Belt (~4.3 Ga)on the eastern shore of Hudson Bay. Even that age remains disputed. But a few, tiny detrital zircon grains extracted from much younger sandstones yield an age range up to 4.4 Ga: barely enough to refute James Hutton’s ‘No vestige of a beginning’. So, the Hadean is long on speculation, most based on less than 3% of all the detrital zircon grains that have been dated. What zircons can tell us is based on their isotopic geochemistry, their trace-element content and even tinier granules of a few other minerals that they encapsulate. The data from them suggest the presence of some kind of felsic magma production that crystallised at low temperatures (~700° C) and was exposed to watery fluids. All very vague compared with what can be gleaned with confidence from post-4.0 Ga rock exposures. But there is a sound astronomical context and a theory based on geophysical and geochemical processes known from experiment and observation of later geology, that can shed a little light.

The planetary system began to form by gravitational accretion of material in a protoplanetary disc of pre-solar gas and dust. The first step would have been gravitational sticking together of dust particles. Fast when this cloud was dense but slowing as the available starting material was depleted by growing planetesimals. This early accretion would easily have radiated away the heat generated by the gravitational potential energy that was released. But that became less effective as the accreting bodies grew to sizes of tens to hundreds of kilometres. Studies of meteorites, formed by collisions of larger planetesimals, show that they became hot enough to melt their contents and even to undergo internal, geochemical differentiation. The current view of the next step is that gravitational perturbations associated with Jupiter drove bodies ranging from asteroidal to Mars size into chaotic motion through the Solar System. Assembly of protoplanets thereafter was dominated by collisions. In the case of the proto-Earth this involved its collision with another, Mars-sized body, to result in the formation of the Moon and the early Earth, each initially enveloped by magma oceans. This event can be considered to be the starting point for all subsequent geological processes on both bodies. But that did not ‘calm down’ planetary bombardment. Plenty of large asteroids were still around: their size range can be judged roughly from those that remain in the Asteroid Belt, that are up to 940 km across in the case of the dwarf planet Ceres. This repository of Hadean objects is what motivated Tim Johnson of Curtin University, Western Australia and three Australian colleagues to ponder on the influence on the Hadean Earth of far more bodies, large and small, hurtling around the early Solar System (Johnson, T.E. et al. 2026. Impact heating and the hidden Hadean. Science, v. 392, p. 1408-1412; DOI: 10.1126/science.aeb5402. PDF requests to tim.johnson@curtin.edu.au).

The impact history of the Earth has largely been expunged by tectonics, erosion and sedimentary burial. Johnson et al. assumed an early impact flux from the almost pristine ‘stratigraphy’ of lunar cratering scaled up to the roughly 13 times greater gravitational pull of the Earth. They calculated that energy being released by impacts and partly incorporated into the Earth during the Hadean outweighed that being generated by internal radioactive decay by several orders of magnitude. Hadean tectonics was thus thermally dominated by impact energy, whose supply probably fluctuated wildly because of different sizes of impacting bodies. By far the largest crater on the Moon – the South Pole-Aitken basin – is 2500 km across. It formed about 4.3 Ga ago when a body 200 km wide struck the lunar surface. Being larger and having a greater gravitational pull, Earth would have suffered up to ten collisions of this magnitude.

In Archaean and later times tectonics became the main means of shedding ‘smoothly’ generated internal radiogenic heating. Dated lunar rock samples strongly suggest that such awesome bombardment lasted until the early Archaean, around 3.8 Ga ago. Traces of this Late Heavy Bombardment are anomalous tungsten isotopes in gneisses of that age from West Greenland (see: Tungsten and Archaean heavy bombardment; July 2002). Internal heating now governs the physical behaviour of rock: whether it is ductile or brittle. Modern-style lithosphere is brittle, hence plate tectonics. The mantle beneath, in the long term, behaves in a ductile fashion, hence convection. As thermal energy built up with each massive impact neither thermal conduction nor bulk convection in the deeper mantle – i.e. the general state of Earth’s present thermal balance – would have been sufficient to check its effects. Rock would need to melt and magma move rapidly in vast amounts to the surface to dissipate energy by radiation into space: by far the most efficient planet-cooling process. The authors also modelled the geotherm – the variation of temperature with depth – established by conductive heat loss and radiation from the surface under Hadean conditions. This is shown in the figure and explained below.

Melting conditions in an early Hadean basaltic crust. Credit: Johnson et al., Fig 4

The thick white line is the modelled conductive geotherm for the ‘coolest’ impact-heating scenario; a usually safe scientific approach. The thin white line shows beginning of melting of hydrous basaltic crust: the ‘mafic solidus’ – the blue area to its left remains solid. The dark to light green shading towards the right marks increasing percentages of basalt melting in 10% steps (dashed white lines). The palest area at right represents a completely molten crust, beyond the ‘mafic liquidus’. The dashed purple line is the liquidus of mantle peridotite. Moving leftwards, the solid purple, pink and orange lines represent the beginning of melting (solidus) for peridotite, anhydrous basalt and sodium-rich granite respectively

The modelled Hadean geotherm shows very rapid temperature increase down to about 7.5 km. It passes across the solidi of granite, hydrous basalt, anhydrous basalt and mantle peridotite: everything begins to melt. Clearly, whatever its composition, the uppermost Hadean crust would have been in a partially molten condition below about 3.5 km. At depths of 10 km or more, between 40 to 70 % of basalt would be molten. The distinction between brittle and ductile becomes meaningless in the light of Johnson et al.’s analysis of Hadean impact heating. Not only does the modelling rule out any rigid lithosphere and plate tectonics during the Hadean, it also explains the almost complete absence today of tangible Hadean rock. In particular, continental crust dominated by granitic rocks was probably recycled continually and literally into the Hadean ‘melting pot’. Convection would have dominated Hadean tectonics, but rather than taking the modern form of isolated plumes it would have been chaotic.

Simulated convective patterns for a Hadean upper mantle subject only to radiogenic heating (A) compared with its dynamic behaviour when heated by continuous heavy bombardment The grey areas represent dense residues left by very high degrees of partial melting at more shallow depths (B). Credit: Johnson et al., Fig 3 A and B.

Suddenly, beginning about 3.9 Ga a rich record of albeit disputed tectonics emerges during the Palaeoarchaean and then evolves onwards to modern planetary behaviour. The heavy bombardment had stopped.

See also: Asteroid assault made ancient Earth too hot and chaotic for continents to form. EurekAlert; 25 June 2026. Why Earth Could Not Hold On to Its First Continents Until the Asteroids Stopped Falling. Science Blog; 25 June 2026.

Magmatism and water in the mantle

That Earth has always been such an active planet is largely due to water having continually being shifted into the mantle by subduction of oceanic lithosphere. Emplaced at temperatures around 1200°C the basaltic crust and ultramafic rocks of the lithospheric mantle become hydrothermally altered by interaction with the ocean, so that they contain a range of hydrous minerals. The mantle is estimated to contain between a quarter and four times the present volume of all ocean water. The vast bulk of the mantle is not undergoing partial melting at any one time. Most magmatic activity is linked to plate tectonics through linear belts such as those along oceanic rift systems and above subduction zones, with a small proportion at ocean islands above isolated mantle plumes and similarly, though more sparsely above hot spots below continents. Such plume-related magmatism has at intervals in the past been vastly bigger than now, forming flood-basalt provinces and ocean-floor plateaus, such as the Deccan Traps and the Ontong Java Plateau; some linked to extinction events. The particular role of water in mantle melting is its reduction in the temperature at which melting begins at depth: ‘dry’ mantle does not melt but remains solid, albeit ductile. Two recent studies have provided important insights into previously unsuspected roles that water can play deep in the mantle.

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

Jianfeng Yang of the Chinese Academy of Science and colleagues from China and University of Padua, Italy provide evidence that ancient subducted slabs that gather at the mantle transition zone (MTZ) may trigger ocean-island  and ocean-plateau volcanism  (Yang, J. et al. 2026.  Subduction legacies in the mantle transition zone modulate intraplate oceanic volcanism. Nature Communications, in press; DOI: 10.1038/s41467-026-73403-7). In fact there is a multiplicity of modes wherein subducted slabs interact with the MTZ, some are retained within it while, in one way or another, others eventually pass through it to the deeper mantle.  Long-dead relics of subduction zones trapped there form ‘reservoirs’ of water 410 to 660 km below the surface at concentrations far higher (1 to 3 %) than does pristine mantle (less than 0.1%). It is stored as OH ions in dense mafic minerals, such as ringwoodite a high-pressure form of olivine (Mg2SiO4) containing up to 2.6 % of OH ions, and bridgmanite (MgSiO3), which forms once subducted slabs pass into the mantle transition zone. If that transformed lithosphere rises above about 410 km, such minerals transform back into anhydrous olivine, thereby liberating their water. At such depths, where temperature in the surrounding dry mantle is about 1800°C the emergence of water triggers a decrease in the temperature at which the ancient slab and also the surrounding mantle can melt. The authors cite evidence that such a process has contributed to the Azores oceanic plateau where the crust is 10 to 20 km thick. It is conceivable that a similar process of deep water ‘recycling’ may have been associated with continental flood basalts. Yang et al.’s new insight may also help unravel hitherto puzzling geochemical anomalies in other kinds of basaltic igneous rocks, such as those which well-up at mid ocean ridges to form modern oceanic crust.

Slabs that descend deeper into the mantle retain their dense mafic minerals and thus the water trapped within them. That water may eventually be involved in transformations at much higher pressures and temperatures, as deep as the core-mantle boundary. One possibility is their retention in mantle plumes that rise from the CMB to facilitate partial melting once they pass through the MTZ

See also: Subduction legacies shape intraplate ocean volcanoes. Scienmag, 20 May 2026

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.

The ‘boring billion’ years of the Mesoproterozoic: plate tectonics and the eukaryotes

The emergence of the eukaryotes – of which we are a late-entry member – has been debated for quite a while. In 2023 Earth-logs reportedthat a study of ‘biomarker’ organic chemicals in Proterozoic sediments suggests that eukaryotes cannot be traced back further than about 900 Ma ago using such an approach. At about the same time another biomarker study showed signs of a eukaryote presence at around 1050 Ma. Both outcomes seriously contradicted a ‘molecular-clock’ approach based on the DNA of modern members of the Eukarya and estimates of the rate of genetic mutation. That method sought to deduce the time in the past when the last eukaryotic common ancestor (LECA) appeared. It pointed to about 2 Ga ago, i.e. a few hundred million years after the Great Oxygenation Event got underway. Since eukaryote metabolism depends on oxygen, the molecular-clock result seems reasonable. The biomarker evidence does not. But were the Palaeo- and Mesoproterozoic Eras truly ‘boring’? A recent paper by Dietmar Müller and colleagues from the Universities of Sydney and Adelaide, Australia definitely shows that geologically they were far from that (Müller, R.D. et al. 2025. Mid-Proterozoic expansion of passive margins and reduction in volcanic outgassing supported marine oxygenation and eukaryogenesis. Earth and Planetary Science Letters, v. 672; DOI: 10.1016/j.epsl.2025.119683).

Carbon influx (million tons per year) into tectonic plates and into the ocean-atmosphere system from 1800 Ma to present. The colour bands represent: total carbon influx into the atmosphere (mauve); sequestered in tectonic plates (green); net atmospheric influx i.e. total minus carbon sequestered into plates (orange). The widths of the bands show the uncertainties of the calculated masses shown as darker coloured lines.

From 1800 to 800 Ma two supercontinents– Nuna-Columbia and Rodinia – aggregated nearly all existing continental masses, and then broke apart. Continents had collided and then split asunder to drift. So plate tectonics was very active and encompassed the entire planet, as Müller et al’s palaeogeographic animation reveals dramatically. Tectonics behaved in much the same fashion through the succeeding Neoproterozoic and Phanerozoic to build-up then fragment the more familiar supercontinent of Pangaea. Such dynamic events emit magma to form new oceanic lithosphere at oceanic rift systems and arc volcanoes above subduction zones, interspersed with plume-related large igneous provinces and they wax and wane. Inevitably, such partial melting delivered carbon dioxide to the atmosphere. Reaction on land and in the rubbly flanks of spreading ridges between new lithosphere and dissolved CO2 drew down and sequestered some of that gas in the form of solid carbonate minerals. Continental collisions raised the land surface and the pace of weathering, which also acted as a carbon sink. But they also involved metamorphism that released carbon dioxide from limestones involved in the crustal transformation. This protracted and changing tectonic evolution is completely bound up through the rock cycle with geochemical change in the carbon cycle.

From the latest knowledge of the tectonic and other factors behind the accretion and break-up of Nuna and Rodinia, Müller et al. were able to model the changes in the carbon cycle during the ‘boring billion’ and their effects on climate and the chemistry of the oceans. For instance, about 1.46 Ga ago, the total length of continental margins doubled while Nuna broke apart. That would have hugely increased the area of shallow shelf seas where living processes would have been concentrated, including the photosynthetic emission of oxygen. In an evolutionary sense this increased, diversified and separated the ecological niches in which evolution could prosper. It also increased the sequestration of greenhouse gas through reactions on the flanks of a multiplicity of oceanic rift systems, thereby cooling the planet. Translating this into a geochemical model of the changing carbon cycle (see figure) suggests that the rate of carbon addition to the atmosphere (outgassing) halved during the Mesoproterozoic. The carbon cycle and probable global cooling bound up with Nuna’s breakup ended with the start of Rodinia’s aggregation about 1000 Ma ago and the time that biomarkers first indicate the presence of eukaryotes.

Simplified structures of (a) a prokaryote cell; (b) a simple eukaryote animal cell. Plants also contain organelles called chloroplasts

So, did tectonics play a major role in the rise of the Eukarya? Well, of course it did, as much as it was subsequently the changing background to the appearance of the Ediacaran animals and the evolutionary carnival of the Phanerozoic. But did it affect the billion-year delay of ‘eukaryogenesis’ during prolonged availability of the oxygen that such a biological revolution demanded? Possibly not. Lyn Margulis’s hypothesis of the origin of the basic eukaryote cell by a process of ‘endosymbiosis’ is still the best candidate 50 years on. She suggested that such cells were built from various forms of bacteria and archaea successively being engulfed within a cell wall to function together through symbiosis. Compared with prokaryote cells those of the eukaryotes are enormously complex. At each stage the symbionts had to be or become compatible to survive. It is highly unlikely that all components entered the relationship together. Each possible kind of cell assembly was also subject to evolutionary pressures. This clearly was a slow evolutionary process, probably only surviving from stage to stage because of the global presence of a little oxygen. But the eukaryote cell may also have been forced to restart again and again until a stable form emerged.

See also: New Clues Show Earth’s “Boring Billion” Sparked the Rise of Life. SciTechDaily, 3  November 2025

The final closure of the Iapetus Ocean

A symposium hosted by the Royal Society in 1965 aimed at resurrecting Alfred Wegener’s hypothesis of continental drift. During the half century since Wegener made his proposal in 1915, it had been studiously ignored by most geologists. The majority had bumbled along with the fixist ideology of their Victorian predecessors. The symposium launched what can only be regarded as a revolution in the Earth Sciences. In the three years following the symposium, the basic elements of plate tectonics had emerged from a flurry of papers, mainly centred on geophysical evidence. Geology itself became part of this cause célèbre through young scientists eager to make a name for themselves. The geological history of Britain, together with that of the eastern North America, became beneficiaries only four years after the Royal Society meeting (Dewey, J. 1969. Evolution of the Appalachian/Caledonian Orogen. Nature 222, 124–129; DOI: 10.1038/222124a0).

In Britain John Dewey, like a few other geologists, saw plate theory as key to understanding the many peculiarities revealed by geological structure, igneous activity and stratigraphy of the early Palaeozoic. These included very different Cambrian and Ordovician fossil assemblages in Scotland and Wales, now only a few hundred kilometres apart. The Cambro-Ordovician of NW Scotland was bounded to the SE by a belt of highly deformed and metamorphosed Proterozoic to Ordovician sediments and volcanics forming the Scottish Highlands. That was terminated to the SE by a gigantic fault zone containing slivers of possible oceanic lithosphere. The contorted and ‘shuffled’ Ordovician and Silurian sediments of the Southern Uplands of Scotland. The oldest strata seemed to have ocean-floor affinities, being deposited on another sliver of ophiolites.  A few tens of km south of that there was a very different Lower Palaeozoic stratigraphy in the Lake District of northern England. It included volcanic rocks with affinities to those of modern island arcs. A gap covered by only mildly deformed later Palaeozoic shelf and terrestrial sediments, dotted by inliers of Proterozoic sediments and volcanics separated the Lake District from yet another Lower Palaeozoic assembly of arc volcanics and marine sediments in Wales. Intervening in Anglesey was another Proterozoic block of deformed sediments that also included ophiolites.

Dewey’s tectonic assessment from this geological hodge-podge, which had made Britain irresistible to geologists through the 19th and early 20th centuries, was that it had resulted from blocks of crust (terranes), once separated by thousands of kilometres, being driven into each other. Britain was thus formed by the evolution and eventual destruction of an early Palaeozoic ocean, Iapetus: a product of plate tectonics. Scotland had a fundamentally different history from England and Wales; the unification of several terranes having taken over 150 Ma of diverse tectonic processes. Dewey concluded that the line of final convergence lay at a now dead, major subduction zone – the Iapetus Suture – roughly beneath the Solway Firth. During the 56 years since Dewey’s seminal paper on the Caledonian-Appalachian Orogeny details and modifications have been added at a rate of around one to two publications per year. The latest seeks to date when and where the accretion of 6 or 7 terranes was finally completed (Waldron, J.W.F. et al. 2025. Is Britain divided by an Acadian suture?  Geology, v. 53, p. 847–852; DOI: 10.1130/G53431.1).

Kernel density plots – smoothed versions of histograms – of detrital zircon ages in Silurian and Devonian sandstones from Wales. The bracketed words are stratigraphic epochs. Credit: Waldron et al. 2025, Fig 3A

John Waldron and colleagues from the University of Alberta and Acadia University in Canada and the British Geological Survey addressed this issue by extracting zircons from four late Silurian and early Devonian sandstones in North and South Wales. These sediments had been deposited between 433 and 393 Ma ago at the southernmost edge of the British Caledonide terrane assemblage towards the end of terrane assembly. The team dated roughly 250 zircons from each sandstone using the 207Pb/206Pb and 206Pb/238U methods. Each produced a range of ages, presumed to be those of igneous rocks from whose magma the zircon grains had crystallised. These data are expressed as plots of probable frequency against age.  Each pattern of ages is assumed to be a ‘fingerprint’ for the continental crust from which the zircons were eroded and transported to their resting place in their host sediment. In this case, the researchers were hoping to see signs of continental crust from the other side of the Caledonian orogen; i.e. from the Precambrian basement of the Laurentia continent.

The three late-Silurian sediments showed distinct zircon-age peaks around 600 Ma and a spread of smaller peaks extending to 2.2 Ga. This tallied with a sediment source in Africa, from which the southernmost Caledonian terrane was said to have split and moved northwards.  The Devonian sediment lacked signs of such an African ‘heritage’ but had a prominent age peak at about 1.0 Ga, absent from the Welsh Silurian sediments.  Not only is this a sign of different sediment provenance but closely follows the known age of a widespread magmatic pulse in the Laurentian continent. So, sediment transport from the opposite side of the Iapetus Ocean across the entire Caledonian orogenic belt was only possible after the end of the Silurian Period at around 410 Ma. There must have been an intervening barrier to sediment movement from Laurentia before that, such as deep ocean water further north. Previous studies from more northern Caledonian terranes show that Laurentian zircons arrived in the Southern Uplands of Scotland and the English Lake District around 432 Ma in the mid-Silurian. Waldron et al. suggest, on these grounds that the suture marking the final closure of the Iapetus Ocean lies between the English Lake District and Anglesey, rather than beneath the Solway. They hint that the late-Silurian to early Devonian granite magmatism that permeated the northern parts of the Caledonian-Appalachian orogen formed above northward subduction of the last relics of Iapetus, which presaged widespread crustal thickening known as the Acadian orogeny in North America.

Readers interested in this episode of Earth history should download Waldron et al.’s paper for its excellent graphics, which cannot be reproduced adequately here.

Gravity survey reveals signs of Archaean tectonics in Canadian Shield

Much of the Archaean Eon is represented by cratons, which occur at the core of continental parts of tectonic plates. Having low geothermal heat flow they are the most rigid parts of the continental crust.  The Superior Craton is an area that makes up much of the eastern part of the Canadian Shield, and formed during the Late Archaean from ~4.3 to 2.6 billion years (Ga) ago. Covering an area in excess of 1.5 million km2, it is the world’s largest craton. One of its most intensely studied components is the Abitibi Terrane, which hosts many mines. A granite-greenstone terrain, it consists of volcano-sedimentary supracrustal rocks in several typically linear greenstone belts separated by areas of mainly intrusive granitic bodies. Many Archaean terrains show much the same ‘stripey’ aspect on the grand scale. Greenstone belts are dominated by metamorphosed basaltic volcanic rock, together with lesser proportions of ultramafic lavas and intrusions, and overlying metasedimentary rocks, also of Archaean age. Various hypotheses have been suggested for the formation of granite-greenstone terrains, the latest turning to a process of ‘sagduction’. However the relative flat nature of cratonic areas tells geologists little about their deeper parts. They tend to have resisted large-scale later deformation by their very nature, so none have been tilted or wholly obducted onto other such stable crustal masses during later collisional tectonic processes. Geophysics does offer insights however, using seismic profiling, geomagnetic and gravity surveys.

The Geological Survey of Canada has produced masses of geophysical data as a means of coping with the vast size and logistical challenges of the Canadian Shield. Recently five Canadian geoscientists have used gravity data from the Canadian Geodetic Survey to model the deep crust beneath the huge Abitibi granite-greenstone terrain, specifically addressing variations in its density in three dimensions. They also used cross sections produced by seismic reflection and refraction data along 2-D survey lines (Galley, C. et al. 2025. Archean rifts and triple-junctions revealed by gravity modeling of the southern Superior Craton. Nature Communications, v. 16, article 8872; DOI: 10.1038/s41467-025-63931-z). The group found that entirely new insights emerge from the variation in crustal density down to its base at the Moho (Mohorovičić discontinuity). These data show large linear bulges in the Moho separated by broad zones of thicker crust.

Geology of the Abitibi Terrane (upper),; Depth to the Moho beneath the Abitibi Terrane with rifts and VMS deposits superimposed (lower). Credit: After Galley et al. Figs 1 and 5.

Galley et al. suggest that the zones are former sites of lithospheric extensional tectonics and crustal thinning: rifts from which ultramafic to mafic magmas emerged. They consider them to be akin to modern mid-ocean and continental rifts. Most of the rifts roughly parallel the trend of the greenstone belts and the large, long-lived faults that run west to east across the Abitibi Terrain. This suggests that rifts formed under the more ductile lithospheric condition of the Neoarchaean set the gross fabric of the granites and greenstones. Moreover, there are signs of two triple junctions where three rifts converge: fundamental features of modern plate tectonics. However, both rifts and junctions are on a smaller scale than those active at present. The rift patterns suggest plate tectonics in miniature, perhaps indicative of more vigorous mantle convection during the Archaean Eon.

There is an interesting spin-off. The Abitibi Terrane is rich in a variety of mineral resources, especially volcanic massive-sulfide deposits (VMS). Most of them are associated with the suggested rift zones. Such deposits form through sea-floor hydrothermal processes, which Archaean rifting and triple junctions would have focused to generate clusters of ‘black smokers’ precipitating large amounts of metal sulfides. Galley et al’s work is set to be applied to other large cratons, including those that formed earlier in the Archaean: the Pilbara and Kaapvaal cratons of Australia and South Africa. That could yield better insights into earlier tectonic processes and test some of the hypotheses proposed for them

See also: Archaean Rifts, Triple Junctions Mapped via Gravity Modeling. Scienmag, 6 October 2025

Sagduction of greenstone belts and formation of Archaean continental crust

Simplified geological map of the Archaean Yilgarn Craton in Western Australia. Credit: Geological Survey of Western Australia

Every ancient craton seen from space shows patterns that are unique to Archaean continental crust: elongated, ‘canoe-shaped’ greenstone belts enveloped by granitic gneisses, both of which are punctured by domes of younger, less deformed granites. The Yilgarn Craton of Western Australia is a typical granite-greenstone terrain. Greenstone belts contain lavas of ultramafic, basaltic and andesitic compositions, which in undeformed settings show the typical pillow structures formed by submarine volcanic extrusion. There are also layered mafic to ultramafic complexes, formed by fractional crystallisation, minor sedimentary sequences and occasionally more felsic lavas and ashes. The enveloping grey gneisses are dominantly highly deformed tonalite-trondhjemite-granodiorite (TTG) composition that suggest that they formed from large volumes of sodium-rich, silicic magmas, probably generated at depth by partial melting of hydrated basaltic rocks.

The heat producing radioactive isotopes of potassium, uranium and thorium in both the Archaean mantle and crust would have been more abundant before 2.5 Ga ago, because they decay over time. Consequently the Earth’s interior would have then generated more heat than now, gradually to escape by thermal conduction towards the cooler surface. The presence of pillow lavas and detrital sediments in greenstone belts indicate that surface temperatures during the Archaean Eon were below the boiling point of water; in fact probably much the same as in the tropics at present. Indeed there is evidence that Earth was then a water world. It may even have been so during the Hadean, as revealed by the oxygen-isotope data in 4.4 Ga zircon grains. The broad conclusion from such findings is that the Archaean geothermal gradient was much steeper; there would have been a greater temperature increase with depth than now and new crust would have cooled more slowly. Subduction of cool lithosphere would have been less likely than in later times, especially as higher mantle heat production would have generated new crust more quickly. Another likely possibility is that far more heat would have been moved by convection: there would have been more mantle-penetrating plumes and they would have been larger. Large mantle plumes of the Phanerozoic have generated vast ocean floor plateaus, such as the Kerguelen and Ontong Java Plateau.

A group of geoscience researchers at The University of Hong Kong and international colleagues recently completed a geological and geochemical study of the North China Craton, analysing their data in the light of recently emerging views on Archaean processes (Dingyi Zhao et al, A two-stage mantle plume-sagduction origin of Archean continental crust revealed by water and oxygen isotopes of TTGs, Science Advances, v. 11, article eadr9513  ; DOI: 10.1126/sciadv.adr9513).They found compelling evidence that ~2.5 Ga-old Neoarchaean TTG gneisses in the North China granite-greenstone terrain formed by partial melting of an earlier mafic-ultramafic greenstone crust with high water content. They consider this to support a two-stage model for the generation of the North China Craton’s crust above a vast mantle plume. The first stage at around 2.7 Ga was the arrival of the plume at the base of the lithosphere, which partially melted as a result of the decompression of the rising ultramafic plume. The resulting mafic magma created an oceanic plateau partly by underplating the older lithosphere, intruding it and erupting onto the older ocean floor. This created the precursors of the craton’s greenstones, the upper part of which interacted directly with seawater to become hydrothermally altered. They underwent minor partial melting to produce small TTG intrusions. A second plume arriving at ~2.5 Ga resulted in sinking of the greenstones under their own weight to mix or ‘hybridise’ with the re-heated lower crust. This caused the greenstones substantially to partially melt and so generate voluminous TTG magmas that rose as the greenstones subsided. . It seems likely that this dynamic, hot environment deformed the TTGs as they rose to create the grey gneisses so typical of Archaean granite-greenstone terranes. [Note: The key evidence for Dingyi Zhao et al.’s conclusions is that the two TTG pulses yielded the 2.7 and 2.5 Ga ages, and show significantly different oxygen isotope data (δ18O)].

Two stages of TTG gneiss formation in the North China Craton and the sinking (sagduction) of greenstone belts in the second phase. Credit: Dingyi Zhao et al., Fig 4)

Such a petrogenetic scenario, termed sagduction by Dingyi Zhao and colleagues, also helps explain the unique keel-like nature of greenstone belts, and abundant evidence of vertical tectonics in many Archaean terrains (see: Vertical tectonics and formation of Archaean crust; January 2002), Their model is not entirely new, but is better supported by data than earlier, more speculative ideas. That such processes have been recognised in the Neoarchaean – the North China Craton is one of the youngest granite-greenstone terrains – may well apply to far older Archaean continental crust generation. It is perhaps the last of a series of such events that began in the Hadean, as summarised in the previous Earth-logs post.

How the earliest continental crust may have formed

Detrital zircon grains extracted from sandstones deposited ~3 billion year (Ga) ago in Western Australia yield the ages at which these grains crystallised. The oldest formed at about 4.4 Ga; only 150 Ma after the origin of the Earth (4.55 Ga). Various lines of evidence suggest that they originally crystallized from magmas with roughly andesitic compositions, which some geochemists suggest to have formed the first continental crust (see: Zircons and early continents no longer to be sneezed at; February 2006). So far, no actual rocks of that age and composition have come to light. The oldest of these zircon grains also contain anomalously high levels of 18O, a sign that water played a role in the formation of these silicic magmas. Modern andesitic magmas – ultimately the source of most continental crust – typically form above steeply-dipping subduction zones where fluids expelled from descending oceanic crust encourage partial melting of the overriding lithospheric mantle. Higher radiogenic heat production in the Hadean and the early Archaean would probably have ensured that the increased density of later oceanic lithosphere needed for steep subduction could not have been achieved. If subduction occurred at all, it would have been at a shallow angle and unable to exert the slab-pull force that perpetuated plate tectonics in later times (see: Formation of continents without subduction, March, 2017).

Landsat image mosaic of the Palaeoarchaean granite-greenstone terrain of the Pilbara Craton, Western Australia. Granite bodies show as pale blobs, the volcanic and sedimentary greenstone belts in shades of grey.

Geoscientists have been trying to resolve this paradox for quite a while. Now a group from Australia, Germany and Austria have made what seems to be an important advance (Hartnady, M. I. H and 8 others 2025. Incipient continent formation by shallow melting of an altered mafic protocrust. Nature Communications, v. 16, article 4557; DOI: 10.1038/s41467-025-59075-9). It emerged from their geochemical studies of rocks in the Pilbara Craton of Western Australia that are about a billion years younger than the aforementioned ancient zircon grains. These are high-grade Palaeoarchaean metamorphic rocks known as migmatites that lie beneath lower-grade ‘granite-greenstone’ terrains that dominate the Craton, which Proterozoic deformation has forced to the surface. Their bulk composition is that of basalt which has been converted to amphibolite by high temperature, low pressure metamorphism (680 to 730°C at a depth of about 30 km). These metabasic rocks are laced with irregular streaks and patches of pale coloured rock made up mainly of sodium-rich feldspar and quartz, some of which cut across the foliation of the amphibolites. The authors interpret these as products of partial melting during metamorphism, and they show signs of having crystallised from a water-rich magma; i.e. their parental basaltic crust had been hydrothermally altered, probably by seawater soon after it formed. The composition of the melt rocks is that of trondhjemite, one of the most common types of granite found in Archaean continental crust. Interestingly, small amounts of trondhjemite are found in modern oceanic crust and ophiolites.

A typical migmatite from Antarctica showing dark amphibolites laced with quartzofeldspathic products of partial melting. Credit: Lunar and Planetary Laboratory, University of Arizona

The authors radiometrically dated zircon and titanite (CaTiSiO₅) – otherwise known as sphene – in the trondhjemites, to give an age of 3565 Ma. The metamorphism and partial melting took place around 30 Ma before the overlying granite-greenstone assemblages formed. They regard the amphibolites as the Palaeoarchaean equivalent of basaltic oceanic crust. Under the higher heat production of the time such primary crust would probably have approached the thickness of that at modern oceanic plateaux, such as Iceland and Ontong-Java, that formed above large mantle plumes. Michael Hartnady and colleagues surmise that this intracrustal partial melting formed a nucleus on which the Pilbara granite-greenstone terrain formed as the oldest substantial component of the Australian continent. The same nucleation may have occurred during the formation of similar early Archaean terrains that form the cores of most cratons that occur in all modern continents.

Modelling climate change since the Devonian

A consortium of geoscientists from Australia, Britain and France, led by Andrew Merdith of the University of Adelaide examines the likely climate cooling mechanisms that may have set off the two great ‘icehouse’ intervals in the last 541 Ma (Merdith, A.S. et al. 2025. Phanerozoic icehouse climates as the result of multiple solid-Earth cooling mechanisms. Science Advances, v. 11, article eadm9798: DOI: 10.1126/sciadv.adm9798). They consider the first to be the global cooling that began in the latter part of the Devonian culminating in the Carboniferous-Permian icehouse. The second is the Cenozoic global cooling to form the permanent Antarctic ice cap around 34 Ma and culminated in cyclical ice ages on the northern continents after 2.4 Ma during the Pleistocene. They dismiss the 40 Ma long, late Ordovician to early Silurian glaciation that left its imprint on North Africa and South America –  then combined in the Gondwana supercontinent. The data about two of the parameters used in their model – the degree of early colonisation of the continents by plants and their influence on terrestrial weathering are uncertain in that protracted event.  Yet the Hirnantian glaciation reached 20°S at its maximum extent in the Late Ordovician around 444 Ma to cover about a third of Gondwana: it was larger than the present Antarctic ice cap. For that reason, their study spans only Devonian and later times.

Fluctuation in evidence for the extent of glacial conditions since the Devonian: the ‘ice line’ is grey. The count of glacial proxy occurrences in each 10° of latitude through time is shown in the colour key. Credit: Merdith et al., Fig 2A.

Merdith et al. rely on four climatic proxies. The first of these comprises indicators of cold climates, such as glacial dropstones, tillites and evidence in sedimentary rocks of crystals of hydrated calcium carbonate (ikaite – CaCO3.6H2O) that bizarrely forms only at around 0°C . From such occurrences it is possible to define an ‘ice line’ linking different latitudes through geological time. Then there are estimates of global average surface temperature; low-latitude sea surface temperature; and estimates of atmospheric CO2. The ‘ice-line’ data records an additional, long period of glaciation in the Jurassic and early Cretaceous, but evidence does not extend to latitudes lower than 60°. It is regarded by Merdith et al. as an episode of ‘cooling’ rather than an ‘icehouse’. Their model assesses sources and sinks of COsince the Devonian Period.

The main natural source of the principal greenhouse gas CO2 is degassing through volcanism expelled from the mantle and breakdown of carbonate rock in subducted lithosphere. Natural sequestration of carbon involves weathering of exposed rock that releases dissolved CO2 and ions of calcium and magnesium.   A recently compiled set of plate reconstructions that chart the waxing and waning of tectonics since the Devonian Period allows them to model the tectonically driven release of carbon over time, with time scales on the order of tens to hundreds of Ma. The familiar Milanković forcing cycles on the order of tens to hundreds of ka are thus of no significance in Merdith et al.’s  broader conception of icehouse episodes  Their modelling shows high degassing during the Cretaceous, modern levels during the late Palaeozoic and early Mesozoic, and low emissions during the Devonian. The model also suggests that cooling stemmed from variations in the positions and configuration of continents over time.  Another crucial factor is the tempo of exposure of rocks that are most prone to weathering. The most important are rocks of the ocean lithosphere incorporated into the continents to form ophiolite masses. The release of soluble products of weathering into ocean basins through time acts as a fluctuating means of ‘fertilising’ so that more carbon can be sequestered in deep sediments in the form of organisms’ unoxidised tissue and hard parts made of calcium carbonates and phosphates. Less silicate weathering results in a boost to atmospheric CO2.

Only two long, true icehouse episodes emerge from the empirical proxy data, expressed by the ‘ice-line’ plots. Restricting the modelling to single global processes that might be expected to influence degassing or carbon sequestration produces no good fits to the climatic proxy data. Running the model with all the drivers “off” produces more or less continuous icehouse conditions since the Devonian. The model’s climate-related outputs thus imply that many complex processes working together in syncopation may have driven the gross climate vagaries over the last 400 Ma or so. A planet of Earth’s size without such complexity would throughout that period have had a high-CO2 warm climate. According to Andrew Merdith its fluctuation from greenhouse to icehouse conditions in the late Palaeozoic and the Cenozoic were probably due to “coincidental combination of very low rates of global volcanism, and highly dispersed continents with big mountains, which allow for lots of global rainfall and therefore amplify reactions that remove carbon from the atmosphere”.

Geological history is, almost by definition, somewhat rambling. So, despite despite the large investment in seeking a computed explanation of data drawn from the record, the outcome reflects that in a less than coherent account. To state that many complex processes working at once may have driven climate vagaries over the last 400 Ma or so, is hardly a major advance: palaeoclimatologists have said more or less the same for a couple of decades or more, but have mainly proposed single driving mechanisms. One aspect of Merdith et al.’s  results seems to be of particular interest. ‘Icehouse’ conditions seem to be rare events interspersed with broader ice-free periods. We evolved within the mammal-dominated ecosystems on the continents during the latest of these anomalous climatic episodes. And we and those ecosystems now rely on a cool world. As the supervisor of the project commented, ‘Over its long history, the Earth likes it hot, but our human society does not’.

Readers may like to venture into how some philosophers of science deal with a far bigger question; ‘Is intelligent life a rare, chance event throughout the universe?’ That is, might we be alone in the cosmos? In the same issue of Science Advances is a paper centred on just such questions (Mills, D.B. et al. 2025. A reassessment of the “hard-steps” model for the evolution of intelligent life. Science Advances, v. 11, article eads5698; DOI: 10.1126/sciadv.ads5698). It stems from cosmologist Brandon Carter’s ‘Anthropic Principle’ first developed at Nicolas Copernicus’s 500th birthday celebrations in 1973. This has since been much debated by scientists and philosophers – a gross understatement as it knocks the spots off the Drake Equation. To take the edge off what seems to be a daunting task, Mills et al. consider a corollary of the Anthropic Principle, the ‘hard steps model’. That, in a nutshell, postulates that the origin of humanity and its ability to ponder on observations of the universe required a successful evolutionary passage through a number of hard steps. It predicts that such intelligence is ‘exceedingly rare’ in the universe. Icehouse conditions are respectable candidates for evolutionary ‘hard steps’, and in the history of Earth there have been five of them.

A fully revised edition of Steve Drury’s book Stepping Stones: The Making of Our Home World can now be downloaded as a free eBook

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.

Drip tectonics beneath Türkiye

Tectonics and geomorphology of Turkey showing the main fault systems. The Konya basin is enclosed by the grey rectangle at centre. (Credit: Taymaz et al. Geological Society of London, Special Publication 291, p1-16, Fig 1)

The 1.5-2.0 km high Central Anatolian plateau in Türkiye has been rising since ~11 Ma ago: an uplift of about 1 km in the last 8 Ma. However, part of the southern Plateau shows signs of rapidly subsidence that has created the Konya Basin, marked by young lake sediments. Interferometric radar (InSAR) data from the European Space Agency’s Sentinel-1 satellite, which detects active movement of the Earth’s surface, reveal a crude, doughnut-shaped area of the surface that is subsiding at up to 50 mm per year. This ring of subsidence surrounds a core of active uplift that is about 50 km across (see the first figure). Expressed crudely, active subsidence suggests an excess of mass beneath the affected area, whereas uplift implies a mass deficit; in both cases within the lithosphere. So, when the InSAR data were published in 2020, it became clear that the lithosphere beneath Anatolia is doing something very strange.

Vertical velocities affecting the surface in the Konya Basin derived from InSAR data, velocities colour-coded cyan to blue show subsidence, yellow to red suggesting that the surface is rising. (Credit: Andersen et al., Fig 1c)

Canadian and Turkish geophysicists set out to find a tectonic reason for such aberrant behaviour (Andersen, A.J.  et al. 2024. Multistage lithospheric drips control active basin formation within an uplifting orogenic plateau. Nature Communications, v. 15, Article 7899; DOI: 10.1038/s41467-024-52126-7). They wondered if a process known as ‘drip tectonics’, first mooted as an explanation of anomalous features in some mountain belts in 2004 (see: Mantle dripping off mountain roots, October 2004; and A drop off the old block? May 2008) might be applicable to the Anatolian Plateau. The essence of this process is similar to the slab-pull force at the heart of subduction. Burial and cooling of basaltic material in oceanic lithosphere being driven beneath another tectonic plate converts its igneous mineralogy to the metamorphic rock eclogite, whose density exceeds that of mantle rocks. Gravity then acts to pull the changed material downwards. However, Anatolia shows little sign of subduction. But the mantle beneath shows seismic speed anomalies that hint at anomalously dense material.

Seismic tomography shows that in a large volume 100 to 200 km beneath the central part of the Plateau S-waves travel faster than in the surrounding mantle. The higher speed suggests a body that is denser and more rigid than its surroundings. This could be a sinking, detached block of ‘eclogitised’ lithosphere whose disconnection from the remaining continental lithosphere has been causing the uplift of the Plateau that began in the Late Miocene. A smaller high-speed anomaly lies directly under the Konya Basin, but at a shallower depth (50 to 80 km) just beneath the lithosphere-asthenosphere boundary. The authors suggest that this is another piece of the lower lithosphere that is beginning to sink and become a ‘drip’. Still mechanically attached to the lithosphere the sinking dense block is dragging the surface down.

Andersen et al. instead of relying on computer modelling created a laboratory analogue. This consisted of a tank full of a fluid polymer whose viscosity is a thousand times that of maple syrup that represents the Earth’s deep mantle beneath. They mimicked an overlying  plate by a layer of the same material with additional clay to render it more viscous – the model’s lithospheric mantle – with a ‘crust’ made of a sand of ceramic and silica spherules. A dense seed inserted into the model lithospheric mantle began to sink, dragging that material downwards in a ‘drip’. After that ‘drip’ had reached the bottom of the tank hours later, it became clear that another, smaller drip materialised along the track of the first and also began to sink. Monitoring of the surface of the ‘crust’ revealed that the initial drip did result in a basin. But the further down the drip fell the basin gradually became shallower: there was surface uplift. Once the initial drip had ‘bottomed-out’ the basin began to deepen again as the secondary drip formed and slowly moved downwards. The model seems to match the authors’ interpretation of the geophysics beneath the Anatolian Plateau. One drip created the potential for a lesser one, a bit like in inversion of the well-known slo-mo videos of a drop of milk falling into a glass of milk, when following the drop’s entry a smaller drop rebounds from the milky surface.

Cartoons of drip tectonics beneath the Anatolian Plateau. (a) Lower lithosphere detached from beneath Anatolia in the Late Miocene (10 to 8 Ma) descends into the mantle as it is ‘eclogitised’; (b) a smaller block beneath the Konya Basin beginning to ‘drip’, but still attached to the lithosphere. (Credit: Andersen et al., Fig 4)

In Anatolia the last 10 Ma has not been just ups and downs of the surface corresponding to drip tectonics. That was accompanied by volcanism, which can be explained by upwelling of mantle material displaced by lithospheric drips. When mantle rises and the pressure drops partial melting can occur, provided the mantle material rises faster than it can lose heat: adiabatic melting.

The onset of weathering in the late Archaean and stabilisation of the continents

Distribution of exposed Archaean cratons. The blue Proterozoic areas may, in part be underlain by cratons. (Credit: Groves, D.I. & Santosh, M. DOI:10.1016/j.gr.2020.06.008)

About 50% of continental crust is of Archaean age (2.5 to 4.0 Ga) in huge blocks above lithosphere more than 150 km thick. Younger continental lithosphere is significantly thinner – as low as 40 km. Since the end of the Archaean Eon these blocks have remained tectonically stable and only show signs of extensional, brittle fracture that have been exploited by basaltic dyke swarms. Such crystalline monstrosities have remained rigid for 2.5 billion years. They are termed cratons from the Greek word κράτο (kratos) for ‘might’ or ‘strength’. Numbers of cratons have been pushed together by later tectonics to form continental ‘cores’, separated from one another by highly deformed ‘mobile belts’ formed by younger collisional orogenies. Africa and South America have 4 cratons each, Eurasia 6 or 7, the other continents all have one

Considering how much cratons have been stressed by later tectonic forces, their implacable rigidity might seem surprising. This rigidity is thought to be due to cratons’ unusually low amounts of the main heat-producing elements (HPE) potassium, uranium and thorium, the decay of whose radioactive isotopes produces surface heat flow. Cratons have the lowest surface heat flow on the planet, so in bulk they must have low HPE content. This stems from the nature of cratons’ deepest parts: almost anhydrous, once igneous rocks of intermediate average composition known as granulites. They formed by metamorphism of earlier crustal rocks at depths of up to 70km, which drove out most of their original HPEs and water. The upper cratonic crust has much the same complement of HPEs as that of more recent continental crust. This bulk depletion of cratons has maintained unusually low temperatures in their deep continental crust. That has been immune from partial melting and thus ductile deformation since it formed.

Three billion year-old TTG gneiss in the Outer Hebrides, Scotland. (Credit: British Geological Survey)

Jesse Reimink and Andrew Smye of Pennsylvania State University, USA have considered the geochemistry and history of the world’s cratons to address the long-standing issue of their stability and longevity (Reimink, J.R. & Smye, A.J. 2024. Subaerial weathering drove stabilization of continents. Nature, v. 629, online article; DOI: 10.1038/s41586-024-07307-1). Their main focus is on how the Archaean lower crust lost most of it HPEs, and where they went. During much of the Archaean continental crust formed by partial melting of hydrated basaltic rocks at shallow depths. That generated sodium-rich silicic magmas from which the dominant grey tonalite-trondhjemite-granodiorite (TTG) gneisses of Archaean crust formed by extreme ductile deformation. Though TTGs originally contained sufficient heat-producing capacity to make them ductile during the early Archaean there is little evidence that they underwent extensive partial melting themselves. But they did after 3.0 Ga to produce swarms of granite plutons in the upper Archaean crust.

Complementing the late-Archaean granite ‘swarm’ are deep-crustal granulites with low HPE contents, which mainly formed around the same time. The granulites contain highly metamorphosed sedimentary rocks, which seem to have been sliced into the Archaean crust during its ductile deformation phase. Some of them have compositions that suggest that they are derived from clay-rich shales, their proportion reaching about 30% of all granulite-facies metasediments. Clay minerals are the products of chemical weathering of silicon- and aluminium-rich igneous rocks exposed to the atmosphere. When they form, they host K, U and Th. Also, their composition and high initial water contents are conducive to partial melting under high-temperature conditions, to become a source of granitic magmas. Crustal weathering is key to Reimink and Smye’s hypothesis for the development of cratons in the late Archaean.

There is growing evidence that high Archaean heat flow through oceanic lithosphere – the mantle contained more undecayed HPE isotopes than now – reduced its density. As a result Archaean oceanic basins were considerably shallower than they became in later times. Because of the lower volume of the basins during the Archaean, seawater extended across much of the continental surface. For most of the Archaean Eon Earth was a ‘waterworld’, with little subaerial weathering of its TTG upper crust. As the volume of exposed continental crust increased so did surface weathering to form clay minerals that selectively absorbed HPEs. Over time shales became tectonically incorporated deep into the thickening Archaean continental crust to form a zone with increased heat producing capacity and a higher water content. Once deep enough and heated by their own content of HPE they began partially melting to yield voluminous granitic magmas to which they contributed their load of HPEs. Being lower in density than the bulk of TTG crust the granite melts would have risen to reach the upper crust. They also took in HPEs from the deep TTG crust itself. According to Reimink and Smye this would have concentrated continental heat production in the upper crust, leaving the deeper crust drier, less able to melt and assume ductile properties, and thus to create the cratons.

The authors believe that such a redistribution of heat production in the ancient continental crust did not need any major change in global tectonics. All it required was decreasing oceanic heat flow to create deeper and more voluminous ocean basins, allowing more continental surface to emerge above sea level and dynamic burial of sedimentary products of subaerial weathering. They conclude: “The geological record can then be cast in terms of a pre-emergence (TTG-dominated) and post-emergence (granite-dominated) planet.” That seems very neat … but it seems unlikely that samples can be drilled from the depths where the ‘action’ took place. Geologists depend on exposures of Archaean middle to deep crust brought to the surface by fortuitous later tectonics.

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.

Aftershocks of ancient earthquakes

Any major earthquake is likely to be followed by aftershocks. Survivors of seismic devastation live in dread of them for weeks, even months. In reality the fault responsible for the initial event continues to move for longer than that. Commonly, aftershock activity dies down in magnitude and frequency over time, sometimes after a few weeks and in other cases much later to reach ‘normal background seismicity’ for the associated tectonic setting. Near a major plate boundary, such as the San Andreas Fault system in coastal California or the mid-Atlantic Ridge in Iceland, there is a continual risk of damaging seismic events, but the area around each major event becomes less risky a few tens of years afterwards. For instance, the Loma Prieta area on the San Andreas became quiescent sixteen years after the October 1989 Magnitude 6.9 earthquake that wrought havoc in San Francisco – and interrupted a Major League baseball match in the city. The December 1954, Magnitude 7.3 Dixie Valley earthquake in the active extensional zone of Nevada had a longer period of instability: 48 years. There is no fixed period for the aftermath, seismicity ‘stops when it stops’.

Earthquakes of greater than Magnitude 2.5 in eastern North America (see key to magnitudes at lower right). Those shown in blue date from 1568 to 1979, those in red between 1980 and 2016. (Credit: Chen & Liu, Fig 1)

Sometimes devastating earthquakes take place in what seem to be the least likely places: in tectonically ‘stable’ continental plate interiors. A Magnitude 7.9 earthquake in Sichuan Province, central China on 12 May 2008 left 86 thousand dead or missing, 374 thousand injured and 4.2 million homeless. It occurred in a region whose ancient fault systems had had little if any historic activity. One of the best studied records of seismic events in the middle of a continent is in the Mississippi River valley at the Missouri-Kentucky border, USA, near the town of New Madrid. This experienced three major earthquakes in 1811 and 1812 at Magnitudes estimated from 7.0 to 7.4. Seismicity there has continued ever since. Others that occurred long ago in the ‘stable’  North American continental crust were in South Carolina (1886) and southern Quebec, Canada (1663). They and the subsequent, lesser earthquakes that define clusters up to 250 km around them have been studied using spatial statistics (Chen, Y. & Liu, M. 2023. Long-Lived Aftershocks in the New Madrid seismic Zone and the Rest of Stable North America. Journal of Geophysics Research: Solid Earth, v. 128; DOI: 10.1029/2023JB026482). Yuxuan Chen and Mian Lui of Wuhan University, China and the University of Missouri, USA considered the dates of historic events, their estimated magnitudes and their proximity to other events in each cluster. The closer two events are the greater the chance that the later one is an aftershock of the first, although the relationship may also indicate a long-lived deformation process responsible for both. The authors suggest that this ‘nearest-neighbour’ approach may reveal that up to 65% of earthquakes in the New Madrid zone between 1980 and 2016 are aftershocks of the 1811-1812 major earthquake cluster, and a significant number of modern events in South Carolina could similarly relate to the 1886 Charleston earthquake. On the other hand, small modern earthquakes in Quebec are more likely to be part of the regional seismic background than to have any relationship to the large 17th century event.

Earthquakes are manifestations of deep-seated processes, most usually the build-up and release of strain in the lithosphere. If such processes persist they can result in long-lived earthquake swarms. So both delayed aftershocks and a high background of seismicity can contribute to the mapped clusters of historic events: a blend of relics of the past and modern deformation. They are yet to be detected in earthquake records associated with tectonic plate boundaries. A long history of movements within continents suggests that it is possible that long-delayed aftershocks may masquerade as foreshocks that presage greater events that are pending. Chen and Liu’s nearest-neighbour approach may therefore distinguish false alarms from real risk of major seismic motions.

See also: Some of today’s earthquakes may be aftershocks from quakes in the 1800s. Eureka|Alert, 13 November 2023

Plate tectonics loses another of its pioneers: W. Jason Morgan

The theory of plate tectonics had a long gestation. Continental drift, one of its central tenets, was first proposed by the meteorologist Alfred Wegener in 1912. Apart from a few enthusiasts of such a dynamic aspect of geology, such as Alex du Toit and Arthur Holmes, the majority of geoscientists remained with the non-revolutionary fixist ideology of their Victorian predecessors. Wegener’s stumbling block was his proposed driving mechanism – polflucht (flight from the poles) – which assumed that supercontinents had formed in polar regions to be subject to centrifugal force resulting from Earth’s rotation. This broke them apart to be driven towards the Equator. Such a mechanism being easily invalidated, most contemporary geologists preferred to ‘throw Wegener’s  baby out with the bathwater’. Yet every piece of his evidence that continents had moved around and most of his ideas about the nature of their movements were steadily verified and amplified over the next six decades, which attracted more curious and flexible scientists. What is now the central paradigm of the Earth Sciences had to wait for a set of major discoveries in the 1950s and ‘60s enabled by emerging technologies, such as the magnetometers used by Fred Vine and Drummond Matthews to discover sea-floor magnetic striping and thus sea-floor spreading. Their breakthrough presented a plausible mechanism for continental drift and launched a near frenzy of collaborative research among a global milieu of young geoscientists, one of whom being W. Jason Morgan.

W. Jason Morgan outside the Department of Earth Sciences, Princeton University. (Credit: Denise Applewhite, Princeton University)

His initial interest was in the great fracture zones on the floors of the Atlantic and Pacific Oceans. He grasped that each of them was very nearly a great circle. This was a central key to unifying seafloor spreading and continental drift – to move across a spherical surface every point on the seafloor had to follow such a path. Morgan recognised that the fracture zones could only result from rigid plates having to fracture to accommodate that motion. Using spherical geometry he was able to link together ridges, trenches and these huge transform faults with poles of rotation and triple junctions to predict plate motions in a quantitative manner. That insight provided a key to active earthquakes, mountain belts and volcanoes. His scientific unification was a result of genius: in just a few weeks Morgan established the fundamentals of what became known as plate tectonics.

W. Jason Morgan was one of the revolutionaries who made geology dynamic and launched its resurrection from the boring province of damp field workers in anoraks tramping across tracts of extremely puzzling rocks and structures, noses to the ground. He died at the age of 87 on 31 July 2023.

You can read an obituary by his former research student Richard Hey and his son Jason Phipps Morgan together with a fuller account of his career on Wikipedia.

Direct signs of what caused the Palaeocene-Eocene thermal maximum

Until about 56 Ma ago North America and Europe were connected: one of the last relics of the Pangaea supercontinent. Oxygen isotopes and magnesium/calcium ratios in the tests of both surface- and bottom-dwelling foraminifera suggest that around that time global mean surface temperature increased by about 5 to 6°C within 10 to 20 thousand years. The rate of global warming was comparable to that currently being induced by human activities. The Palaeocene-Eocene thermal maximum (PETM) is seen by climatologists as a dreadful warning of times to come in the not so distant future. The PETM event marks the most dramatic biological changes since the mass extinction at the Cretaceous-Palaeogene boundary 10 million years earlier. They included the rapid expansions of mammals and land plants and major extinction of deep-water foraminifera. The PETM also coincided with an equally profound excursion in the δ13C of carbon-rich strata of that age, whose extreme negative value marks the release of a huge mass of previously buried organic carbon into the atmosphere. It was probably methane, much more potent at delaying heat loss to space than carbon dioxide – methane has more than 80 times the warming effect of carbon dioxide. Since CH4 is soon oxidised to CO2 and H2O estimates of atmospheric greenhouse gas levels are generally expressed in terms of CO2. The PETM release was equivalent to about 4.4 x 1013metrictons over 50 ka; on average 0.24 gigatons per year compared with 0.51 Gt from energy-related sources in 2022.

During the Palaeocene, areas around the present North Atlantic were subject to basaltic continental volcanism before the rifting that opened the North Atlantic from 62 to 58 Ma. Magmatism, dominated by intrusions, began again at the Palaeocene-Eocene boundary from 56 to 54 Ma, linked to the start of continental rifting. Both episodes suggest a rising mantle plume. Once the rift had truly opened volcanism became restricted to the mid Atlantic ridge and a mantle plume remains active beneath Iceland. After geoscientists became aware of the PETM and its coincidence with North Atlantic igneous activity many palaeoclimatologists suggested methane release from organic-rich sediments heated by intrusion of basaltic sills below the opening seaway (but see 2022 post on alternative hypotheses). As with so many extreme geological events, choosing a most-likely scenario depends ultimately on tangible evidence. A convincing sign has been demonstrated dramatically in a recent study by a multinational team of geophysicists, oceanographers, geochemists, palaeontologists and sedimentologists (Berndt, C. and 35 others 2023. Shallow-water hydrothermal venting linked to the Palaeocene–Eocene Thermal Maximum. Nature Geoscience, v. 16, p. 803–809; DOI: 10.1038/s41561-023-01246-8).

Three-dimensional view of seismic reflection data off western Norway. The greytone lower part is a vertical ‘slice’. The coloured part shows the depth variation of sediments that fill hydrothermal vent systems beneath a horizontal unconformity. (Credit: Berndt et al, Fig 1b)

The breakthrough by Berndt et al. stemmed from a detailed 3-D seismic survey off the coast of Norway. It revealed an unconformity at the P-E boundary beneath which were clear signs of hundreds of large pockmarks, up to 80 m deep. Seismic reflection from older sediments beneath the unconformity showed the distinctive presence of intrusive sills of igneous rocks. The consortium drilled 20 boreholes into the seabed beneath the survey area. Five of them penetrated crater-like features to yield cores through the sediments that had filled them. The fills were muds, which were interleaved beds of volcanic ash in the sequences marking the P-E boundary suggesting an igneous influence. Organic remains in the muds established the depositional timing of several distinct layers and also gave clues to their depositional conditions. Those spanning the 50 ka of the PETM were dominated by plant debris, pollen and spores, together with abundant marine diatoms that live in very shallow water. Laminations in the muds dip radially inwards towards the deeper parts of some craters to define funnel-like structures. In others the sediments have been domed upwards. The sediments and their structures closely resemble those in blow-out craters formed during petroleum drilling accidents and in onshore maar volcanoes produced by sudden explosive eruptions on land. The pockmarks formed suddenly, to be filled by mobilised mud and volcanic ash.

The evidence points to explosive vents formed by massive degassing of deeper sediments induced by igneous intrusions. Such systems are common around active ocean-floor rifts: ‘black-‘ and ‘white smokers’, but those off Norway formed in shallow water. That has an important bearing on their potency during the PETM. Deep hydrothermal systems may emit methane, but it is oxidised to CO2 in seawater. Those very close to the surface vent their gas almost directly into the atmosphere before such oxidation can consume methane. Intrusive sills also underlie the eastern continental margin of Greenland, so such explosive hydrothermal vents may have been widespread during the initial rifting of the North Atlantic’.

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

New drill core penetrates the Mohorovičić Discontinuity (the ‘Moho’)

In 1909 Croatian geophysicist Andrija Mohorovičić examined seismograms of a shallow earthquake that shook the area around Zagreb. To his surprise the by-then familiar time sequence of P-waves followed by the slower S-waves appeared twice on seismic records up to 800 km away. The only explanation that he could come up with was that the first arrivals had travelled directly through the crust to the detector whereas the second set must have followed a longer path: it had travelled downwards to be refracted to reach the surface when it met rocks denser than those at the surface. His analysis revealed a sharp boundary between the Earth’s crust and its mantle at a depth of about 54 km below what was then Yugoslavia. Later workers confirmed this discovery and honoured its discoverer by naming it the Mohorovičić Discontinuity. Difficulty with pronouncing his name resulted in a geological nickname: ‘the Moho’. It can be detected everywhere: at 20 to 90 km beneath the continental surface and 5 to 10 km beneath the ocean floor, thus distinguishing between continental and oceanic crust.

In the late 1950s accelerating geological and oceanographic research that would culminate in the theory of plate tectonics turned its focus on drilling down to the Moho in much the same way as a lust for space travel spawned getting to the Moon. The difference was that the proposers of what became known as the Mohole Project were members of what amounted to a geoscientific glee club (The American Miscellaneous Society), which included a member of the well-financed US National Science Foundation’s Earth Science Panel. The idea emerged shortly after the Soviet Union had launched the Sputnik satellite and rumours emerged that it was proposing deep drilling into the continental crust beneath the Kola Peninsula.  The Mohole’s initial target was the 3.9 km deep floor of the Caribbean off Guadalupe in Mexico and required advanced methods of stabilisation for a new oceanographic ship that was to host the drilling rig.

Huge (tens of metres high) pillars or ‘chimneys’ of carbonates formed by the Lost City hydrothermal vent near the mid-Atlantic ridge (Credit: ETH Zurich)

The Mohole was spudded in 1961, but the deepest of five holes reached only 200 m beneath the sea floor. It recovered Miocene sediments and a few metres of basalt. Deep water drilling was somewhat more complicated than expected and about US$ 57 million was spent fruitlessly. The project was disbanded in 1966 with considerable acrimony and schadenfreude. Nonetheless, the Mohole fiasco made technical advances and did demonstrate the feasibility of offshore drilling. The petroleum industry benefitted and so did oceanography with the globe-spanning deep-sea drilling of ocean floor sediments. The sediment cores produced the 200 million-year exquisitely detailed record of climate change and vast amounts of geochemical data from the basaltic oceanic crust. In 2005 JOIDES (the Joint Oceanographic Institutions for Deep Earth Sampling) had another crack at the Moho. That venture centred on the intersection of the Mid-Atlantic Ridge and the Atlantis Fracture Zone close to the ‘Lost City’ hydrothermal vent. The area around the vent is the site of a huge low-angled extensional fault that has partly dragged the basaltic ocean crust off the mantle beneath causing it to bulge. This provided an excellent opportunity to drill through the Moho. All went well, but 54 days of drilling yielded 1.4 km of basalt but nothing resembling mantle rock. So, again, the Moho had thwarted Science (and research economics). But finally it is beginning to reveal it secrets (see: Voosen, P. 2023. Ocean drillers exhume a bounty of mantle rocks. Science, v. 380 (News) p. 876-877; DOI: 10.1126/science.adi9899

The area around the ‘Lost City’ vent was originally chosen for drilling to examine the chemical processes going on there. Hydrogen emitted by serpentinisation of mantle rocks can combine with carbon monoxide in hydrothermal fluids to create a wide variety of organic compounds, which could be the initial building blocks for the origin of life. As part of the International Ocean Drilling Programme JOIDES decided to launch IODP Expedition 399 to re-examine the area around ‘Lost City’ in more detail. The expedition first tried to continue drilling the 2005 hole, but failed yet again. Finally a new drill site aimed at penetrating the extensional detachment. Within a few days the drill bit punched into mantle rocks and over a 6-week period the expedition had recovered a kilometre of core. The technical accounts for each week of drilling give a flavour of what it must be like to be a part of such a ship-borne expedition as well as describing what emerged in the drill core. It seems like a bit of a jumble, dominated by the mineral olivine– the principal characteristic of the ultramafic mantle – almost pure in the rock dunite and mixed with pyroxenes in various kinds of peridotite. There are also coarse-grained rocks that contain plagioclase feldspar, which cut through the ultramafic materials – gabbros, troctolites and norites.  They are relics of intrusive basaltic magmas that did not make it to the seabed. The samples are variably altered by interaction with watery hydrothermal fluids, with lots of serpentine, talc and even asbestos: the drilling presented a health hazard for a few days. The rocks have been metamorphosed under pressure-temperature conditions of greenschist to amphibolite facies and subject to ductile deformation, probably because of the effect of extensional deformation. Whatever, there is plenty of material to be analysed, including for signs of microbial activity. So, the dreams of a 1950s academic drinking fraternity (they were all men!) have finally been realised. But since those pre-plate-tectonic times many geologists have seen and collected much the same, even putting their index fingers on the Moho itself in the time-honoured fashion. Intricate 3-D geology in ophiolite complexes such as that in Oman, provide such opportunities at the much lower cost of air travel, Land Cruiser hire and camping. Indeed what we know of the structure of the oceanic lithosphere – pillow lavas, sheeted dyke complexes, gabbro cumulates and serpentinised ultramafic mantle – has come from such bodies thrust onto continental crust at ancient plate margins. So, why the celebration in this case? They are the first samples of mantle from young oceanic lithosphere; the rocks of ophiolites may not have formed at mid-ocean ridges. These should give clues to the long-term magmatism that has created the vast abyssal basins that the mantle eventually reabsorbs by subduction. Then, of course, there is the link to biogenic processes at constructive margins that underpinned the return to the active hydrothermal venting at ‘Lost City’.

End-Ordovician mass extinction, faunal diversification, glaciation and true polar wander

Enormous events occurred between 460 and 435 Ma around the mid-point of the Palaeozoic Era and spanning the Ordovician-Silurian (O-S) boundary. At around 443 Ma the second-most severe mass extinction in Earth’s history occurred, which eliminated 50 to 60% of all marine genera and almost 85% of species: not much less than the Great Dying at the end of the Permian Period. The event was accompanied by one of the greatest biological diversifications known to palaeontology, which largely replaced the global biota initiated by the Cambrian Explosion. Centred on the Saharan region of northern Africa, Late Ordovician glacial deposits also occur in western South America and North America. At that time all the current southern continents and India were assembled in the Gondwana supercontinent, with continental masses that became North America, the Baltic region, Siberia and South China not far off: all the components that eventually collided to form Pangaea from the Late Silurian to the Carboniferous.

The mass extinction has troubled geologists for quite a while. There are few signs of major volcanism having been involved, although some geochemists have suggested that very high mercury concentrations in some Late Ordovician marine sediments bear witness to large, albeit invisible, igneous events. No large impact crater is known from those times, although there is a curious superabundance of extraterrestrial debris, including high helium-3, chromium and iridium concentrations, preserved in earlier Ordovician sedimentary rocks, around the Baltic Sea. Another suggestion, poorly supported by evidence, is destruction of the atmospheric ozone layer by a gamma-ray burst from some distant but stupendous supernova. A better supported idea is that the oceans around the time of the event lacked oxygen. Such anoxia can encourage solution of toxic metals and hydrogen sulfide gas. Unlike other mass extinctions, this one was long-drawn out with several pulses.

The glacial epoch also seems implicated somehow in the mass die-off, being the only one known to coincide with a mass extinction. It included spells of frigidity that exceeded those of the last Pleistocene glacial maximum, with the main ice cap having a volume of from 50 to 250 million cubic kilometres. The greatest of these, around 445 Ma, involved a 5°C fall in global sea-surface temperatures and a large negative spike in δ13C in carbon-rich sediments, both of which lasted for about a million years. The complex events around that time coincided with the highest ever extinction and speciation rates, the number of marine species being halved in a short space of time: a possible explanation for the δ13 C anomaly. Yet estimates of atmospheric CO2 concentration in the Late Ordovician suggests it was perhaps 8–16 times higher than today; Earth should have been a warm planet then. One probable contributor to extreme glacial conditions has been suggested to be that the South Pole at that time was well within Gondwana and thus isolated from the warming effect of the ocean. So, severe glaciation and a paradoxical combination of mass extinction with considerable biological diversification present quite an enigma.

A group of scientists based in Beijing, China set out to check the palaeogeographic position of South China between 460 and 435 Ma and evaluate those in  O-S sediments at locations on 6 present continents (Jing, X., Yang, Z., Mitchell, R.N. et al. 2022. Ordovician–Silurian true polar wander as a mechanism for severe glaciation and mass extinction. Nature Communications, v. 13, article 7941; DOI: 10.1038/s41467-022-35609-3). Their key tool is determining the position of the magnetic poles present at various times in the past from core samples drilled at different levels in these sedimentary sequences. The team aimed to test a hypothesis that in O-S times not only the entire lithosphere but the entire mantle moved relative to the Earth’s axis of rotation, the ‘slippage’ probably being at the Core-mantle boundary [thanks to Steve Rozario for pointing this out]. Such a ‘true polar wander’ spanning 20° over a mere  2 Ma has been detected during the Cretaceous, another case of a 90° shift over 15 Ma may have occurred at the time when Snowball Earth conditions first appeared in the Neoproterozoic around the time when the Rodinia supercontinent broke up and a similar event was proposed in 1994 for C-O times albeit based on sparse and roughly dated palaeomagnetic pole positions.

Xianqing Jing and colleagues report a wholesale 50° rotation of the lithosphere between 450 and 440 Ma that would have involved speeds of about 55 cm per year. It involved the Gondwana supercontinent and other continental masses still isolated from it moving synchronously in the same direction, as shown in the figure. From 460 to 450 Ma the geographic South Pole lay at the centre of the present Sahara. At 445 Ma its position had shifted to central Gondwana during the glacial period. By 440 Gondwana had moved further northwards so that the South Pole then lay at Gondwana’s southernmost extremity.

Palaeogeographic reconstructions charting true polar wander and the synchronised movement of all continental masses between 460 and 440 Ma. Note the changes in the trajectories of lines of latitude on the Mollweide projections. The grey band either side of the palaeo-Equator marks intense chemical weathering in the humid tropics. Credit Jing et al. Fig 5.

As well as a possible key to the brief but extreme glacial episode this astonishing journey by a vast area of lithosphere may help account for the mass extinction with rapid speciation and diversification associated with the O-S boundary. While the South Pole was traversing Gondwana as the supercontinent shifted the ‘satellite’ continental masses remained in or close to the humid tropics, exposed to silicate weathering and erosion. That is a means for extracting CO2 from the atmosphere and launching global cooling, eventually to result in glaciation over a huge tract of Gondwana around 445 Ma. Gondwana then moved rapidly into more clement climatic zones and was deglaciated a few million years later. The rapid movement of the most faunally diverse continental-shelf seas through different climate zones would have condemned earlier species to extinction simultaneous adaptation to changed conditions could have encouraged the appearance of new species and ecosystems. This does not require the catastrophic mechanisms largely established for the other mass extinction events. It seems that during the stupendous, en masse slippage of the Earth’s lithosphere plate tectonic processes still continued, yet it must have had a dynamic effect throughout the underlying mantle.

Yet the fascinating story does have a weak point. What if the position of the magnetic poles shifted during O-S times from their assumed rough coincidence with the geographic poles? In other words, did the self-exciting dynamo in the liquid outer core undergo a large and lengthy wobble? How the outer core’s circulation behaves depends on its depth to the solid core, yet the inner core seems only to have begun solidifying just before the onset of the Cambrian, about 100 Ma before the O-S events. It grew rapidly during the Palaeozoic, so the thickness of the outer core was continuously increasing. Fluid dynamic suggests that the form of its circulation may also have undergone changes, thereby affecting the shape and position of the geomagnetic field: perhaps even shifting its poles away from the geographic poles …

A Lower Jurassic environmental crisis

Curiously, one of the largest environmental disruptions during the Phanerozoic Eon (i.e. since 541 Ma ago) does not stand out in the way that the ‘Big Five’ mass extinctions do. Each of them killed off between 70 and 95% of all marine species. The Jurassic was a period of biological recovery from the End-Triassic extinction 201 Ma ago. Throughout its ~50 Ma duration extinction rates were below the average for the Phanerozoic, and they remained relatively low until the K-Pg mass extinction that drew the Mesozoic Era to a close at 66 Ma. Nevertheless, there were significant extinctions, such as the demise of several lineages of herbivorous dinosaurs towards the end of the Early Jurassic followed by the rise of the familiar, long-necked variety of eusauropods. Marine organisms that secreted hard parts made of calcium carbonate also experienced a collapse then. From time to time during the Jurassic and Cretaceous Periods the oceans lost a great deal of dissolved oxygen, increasing the chances of organic carbon being buried in marine sediments. Such oceanic anoxia resulted in the widespread deposition of hydrocarbon source rocks in the form of black bituminous muds. Overall, both the Jurassic and Cretaceous experienced  greenhouse climatic conditions, with  atmospheric CO2 levels rising to almost 3000 ppm and oxygen levels significantly lower than the modern 21%. Sea levels rose by up to 200 metres, thought to be due to fast sea-floor spreading and large areas of warm, buoyant oceanic lithosphere.

A notable ocean-anoxia event took place during the Lower Jurassic, around 183 Ma ago at the start of the Toarcian Age. This stratigraphic level was penetrated by a 1.5 km borehole sunk in 2015-2016 at Mochras in North Wales, UK, on the shore of Cardigan Bay. The core provided the thickest and most complete record ever recovered for this event, and has been analysed in exquisite detail using many techniques. The most revealing data have been published by a multinational team led by scientists from Trinity College, Dublin (Ruhl, M. et al. 2022. Reduced plate motion controlled timing of Early Jurassic Karoo-Ferrar large igneous province volcanism. Science Advances, v. 8, article eabo0866; DOI: 10.1126/sciadv.abo0866).

Plate boundaries around Gondwanaland and the Karoo-Ferrar large igneous province in the Early Jurassic (small yellow dots show dated localities) . Large pink dots: positions of Tristan de Cunha and Bouvet hotspots at the time (Credit: Ruhl et al. Fig 1A)

At the start of the Toarcian (183.7 Ma) the 187Os/186Os ratio of the samples begins to rise from 0.3 to almost 0.8 to fall back to 0.3 by 180.8 Ma. Osmium isotopes are a measure of continental weathering, and this ‘excursion’ surely signifies significant global warming and increases in atmospheric humidity and acidity that broke down rocks at the continental surface. Over the same period δ13C rises, decreases to by far the lowest value in the Lower Jurassic, rises again to gradually fall back. The start of the Toarcian seems to have experienced a major release of carbon then a profound sequestration of organic carbon, presumably through burial of dead organisms in the black mudstones that signify anoxic conditions. Remarkably, the 95 m thick Toarcian black-mudstone sequence also reveals a tenfold increase in its content of the element mercury, from 20 to 200 parts per billion (ppb), peaking at the same time (~182.8 Ma) as the most negative δ13C value was reached: the acme of carbon sequestration. A coincidence of massive organic carbon burial and increased mercury in marine sediments also happened at the time of the end-Permian mass extinction, although that does not necessarily imply exactly the same mechanism.

The early Toarcian geochemical trends, however, coincide with the initiation and duration of the Karoo-Ferrar large igneous province, which formed flood basalts, igneous dyke swarms and large volcanic centres in South Africa and Antarctica. That LIP may have emitted mercury, but so too may have increased chemical weathering of the land surface. Whichever, mercury forms an organic compound (methyl mercury) in water bodies. Readily incorporated into living organisms, that could explain the close parallel between the δ13C and Hg records in the Jurassic sediment core from Wales. The Karoo-Ferrar igneous activity itself presents a bit of a conundrum, as suggested by Ruhl et al. It happened at the very time that there was a 120° change in the direction of motion of the tectonic plate carrying along Africa and, indeed, the Gondwanaland supercontinent during the Jurassic. The directional change also involved local plate movement stopping for a while. According to the authors, it wasn’t a fortuitous coincidence of two mantle plumes from the core-mantle boundary hitting the bottom of the continental lithosphere below Africa and Antarctica at this tectonic ‘U-turn’. It is more likely that the pause gave existing plumes the opportunity and time to ‘erode’ the base of the continental lithosphere and rise. Decompression melting would then have produced the voluminous magmas. The two plumes were in place for a very long time and created seamount chains as plates moved over them. Both are still volcanically active: Tristan de Cunha on the mid-Atlantic Ridge, and Bouvet Island at a triple junction between South Africa and Antarctica.

So, a venture to unravel a period of profound environmental change during the Early Jurassic, which didn’t result in mass extinction, may well have spawned a new model for massive igneous events that did. Ruhl et al. suggest that the short-lived Siberian, North Atlantic and East African Rift LIPs each seem to have coincided with short episodes of tectonic slowing-down: LIPs may result in dramatic environmental change, but at the whim of plate tectonics.

See also: https://scitechdaily.com/surprising-discovery-shows-how-slowing-of-continental-plate-movement-controlled-earths-largest-volcanic-events/

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

Sun, sand and sangria on the Mediterranean Costas – and tsunamis?

You can easily spot a tourist returning from a few summer weeks on the coast of the western Mediterranean, especially during 2022’s record-breaking heat wave and wildfires: sunburnt and with a smoky aroma that expensive après-sun lotion can’t mask. Judging from the seismic records, they may have felt the odd minor earthquake too, perhaps putting it down to drink, lack of sleep and an overdose of trance music. Data from the last 100 years show that southern Spain and north-west Africa have a generally uniform distribution of seismic events, mostly less than Magnitude 5. Yet there is a distinct submarine zone running NNE to SSW from Almeria to the coast of western Algeria. It crosses the Alboran Basin, and reveals significantly more events greater than M 5. Most earthquakes in the region occurred at depths less than 30 km mainly in the crust. Five geophysicists from Spain and another two from Algeria and Italy have analysed the known seismicity of the region in the light of its tectonics and lithospheric structure (Gómez de la Peña, L., et al. 2022. Evidence for a developing plate boundary in the western Mediterranean. Nature Communications, v. 13, article 4786; DOI: 10.1038/s41467-022-31895-z).

Topography of the Alboran Basin beneath the western Mediterranean. The colours grey through blue to purple indicate increasing depth of seawater. Grey circles indicate historic earthquakes, the smallest being M 3 to 4, the largest greater than M 6. Green arrows show plate motions in the area measured using GPS. Active faults are marked in red (see key for types of motion). (Credit: based on Fig 1 of Gómez de la Peña et al.)

The West Alboran Basin is underlain by thinner continental crust (orange on the inset to the map) than beneath southern Spain and western Algeria. Normal crust underpins the Southern Alboran Basin. To the east are the deeper East Alboran and Algero-Balearic Basins, the floor of the latter being true oceanic crust and that of the former created in a now extinct island arc. Running ENE to WSW across the Alboran Basin are two ridges on the sea floor. Tectonic motions determined using the Global Positioning System reveal that the African plate is moving slowly westwards at up to 1 cm yr-1, about 2 to 3 times faster than the European plate. This reflected by the dextral strike-slip along the active ~E-W Yusuf Fault (YSF). This bends southwards to roughly parallel the Alboran Ridge, and becomes a large thrust fault that shows up on ship borne seismic reflection sections. The reflection seismic survey also shows that the shallow crust beneath the Alboran Ridge is being buckled under compression above the thrust. The thrust extends to the base of the African continental crust, which is beginning to override the arc crust of the East Alboran basin. Effectively, this system of major faults seems to have become a plate boundary between Africa and Europe in the last 5 million years and has taken up about 25 km of convergence between the two plates. An estimated 16 km of this has taken place across the Alboran Ridge Thrust which has detached the overriding African crust from the mantle beneath.

The authors estimate an 8.5 to 10 km depth beneath the Alboran fault system at which the overriding crust changes from ductile to brittle deformation – the threshold for strains being taken up by earthquakes. By comparison with other areas of seismic activity, they reckon that there is a distinct chance of much larger earthquakes (up to M 8) in the geologically near future. A great earthquake in this region, where the Mediterranean narrows towards the Strait of Gibraltar, may generate a devastating tsunami. An extension of the Africa-Europe plate boundary into the Atlantic is believed to have generated a major earthquake that launched a tsunami to destroy Lisbon and batter the Atlantic coasts of Portugal, Spain and NW Africa on 1st November 1755. The situation of the active plate boundary in the Alboran Basin may well present a similar, if not worse, risk of devastation.

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