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

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.

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.

The world’s oldest crust in the Nuvvuagittuq Greenstone Belt, Quebec

Since 1999, the rocks generally acknowledged to be the oldest on Earth were part of the Acasta gneisses in the Slave Craton in Canada’s Northwest Territories; specifically the Idiwhaa tonalitic gneisses. Zircons extracted from that unit yielded an age of 4.02 billion years (Ga) using U-Pb radimetric dating, revealing the time of their crystallisation from granitic magma. But nine years later some metabasaltic rocks from the tiny (20 km2) Nuvvuagittuq Greenstone Belt on the eastern shore of Hudson Bay were dated using the Sm-Nd method at almost 4.3 Ga (see: At last, 4.0 Ga barrier broken; November 2008). Taken at face value the metabasaltic rocks seemed to be well within the Hadean Eon (4.6 to 4.0 Ga) and could thus represent primary crust of that antiquity. However, U-Pb dating of zircons from thin sodium-rich granitic rocks (trondhjemites) that intrude them yielded ages no older than about 3.8 Ga. Similar ages emerged from zircons found in metasediments interleaved in the dominant mafic unit. Discrepancies between the two completely different dating methods resulted in the Hadean antiquity of the mafic rocks having been disputed since 2008. It was possible that the Sm-Nd results from the metabasalts may have resulted from the original mafic magmas having inherited a Hadean Sm-Nd isotopic ‘signature’ from their mantle source. That is, they may have been contaminated and could have formed in the early Archaean.

Glacially smoothed outcrops near Inukjuak, Quebec that reveals rocks of the Nuvvuagittuq Greenstone Belt. Credit: Jonathan O’Neil, University of Ottawa

Jonathan O’Neil, now at Ottawa University in Canada, led the first isotopic investigation of the Nuvvuagittuq Greenstone Belt and has engaged in research there ever since. Further field and laboratory studies revealed that the previously dated mafic rocks had been intruded by large, chemically differentiated gabbro sills. A team of geochemists from the University of Ottawa and Carleton University, including O’Neil, has now published isotopic evidence from the intrusions that suggests a Hadean age for their parent magma (C. Sole et al. 2025. Evidence for Hadean mafic intrusions in the Nuvvuagittuq Greenstone Belt, CanadaScience, v. 388, p. 1431-1435. DOI: 10.1126/science.ads8461). The authors used the decay schemes of two radioactive samarium isotopes 147Sm and 146Sm; a significant advance in radiometric dating. The first decays to 143Nd with a half-life of about 1011 years, the second to 142Nd with a much shorter half life of about 108 years. Due to its more rapid decay, in geological terms,146Sm is now much rarer than 147Sm. Consequently, using the short-lived 146Sm-142Nd decay system is technically more difficult than that of the 147Sm-143Nd system. But the team managed to get good results from both the ‘fast’ and the ‘slow’ decay schemes. They tally nicely, yielding ages of 4157 and 4196 Ma.  The gabbros provide a minimum age for the metabasalts that they cut through. The original 4.3 Ga Sm-Nd date for the metabasalts is thus plausible. Sole and colleagues consider the dominant metabasaltic rocks to have formed a primary crust in late Hadean times that was invaded by later mantle-derived mafic magma about 100 Ma later. The granitic rocks that constitute about one third of the Nuvvuagittuq terrain seem to have been generated by partial melting more than 300 Ma later still, during the Palaeoarchaean.

Perhaps similar techniques will now be deployed in granite-greenstone terrains in other cratons. Many of the older ones, generally designated as Palaeoarchaean in age, also contain abundant metamorphosed mafic and ultramafic igneous rocks. Perhaps their origin was akin to those of Nuvvuagittuq; i.e. more Hadean crust may await unmasking. Meanwhile, there seems to be more to discover from Nuvvuagittuq. For instance, some of the rocks suggested to be metasediments interleaved in the metabasalts show intricate banding that resembles products of bacterial mat accumulation in younger terrains. Signs of Hadean life?

Since the first reliable radiometric dating of Archaean rocks in 1971, there has been an element of competition to date the oldest rocks on Earth: to push history back towards the initial formation of the Earth. It is one of the most disputatious branches of Earth history. Rivalry may play a significant part in driving the science, as well as the development of novel dating techniques and the continuing discovery of clearly old relationships using ‘old-fashioned’ relative dating, such as signs of intrusion, unconformities etcetera. But in some cases there is a darker side: the potential for profit. Recently, samples from Nuvvuagittuq appeared for sale on the Internet, priced at $10,000. They may have been collected under the guise of supplying museums by a group that shipped-in mechanical excavators in 2016. Unsurprisingly this angered the local Innuit community of Inukjuak. They were also worried about bona fide collection for scientific research that had left parts of the small, once pristine area somewhat battered, including cultural features such as an inukshuk navigational monument. Their fury at commercial exploitation of their homeland resulted in the community council closing the area to collecting in 2024. I emphasise that this violation of basic geological ethics was by commercial rock collectors and dealers, not academic geologists. The local people are now considering careful issue of research permits so that important research can continue. But further rock collecting may remain banned.

See also: New Research Verifies Northern Canada Hosts Earth’s Oldest Rocks. Scienmag, 26 June 2025; Gramling, C. 2025. Earth’s oldest rocks may be at least 4.16 billion years old. ScienceNews.

PS With many thanks to ‘Piso Mojado’ for alerting me to this paper

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.

The prospect of climate chaos following major volcano eruptions

It hardly needs saying that volcanoes present a major hazard to people living in close proximity. The inhabitants of the Roman cities of Herculaneum and Pompeii in the shadow of Vesuvius were snuffed out by an incandescent pyroclastic during the 79 CE eruption of the volcano. Since December 2023 long-lasting eruptions from the Sundhnúksgígar crater row on the Reykjanes Penisula of Iceland have driven the inhabitants of nearby Grindavík from their homes, but no injuries or fatalities have been reported. Far worse was the 1815 eruption of Tambora on Sumbawa, Indonesia, when at least 71,000 people perished. But that event had much wider consequences, which lasted into 1817 at least. As well as an ash cloud the huge plume from Tambora injected 28 million tons of sulfur dioxide into the stratosphere. In the form of sulfuric acid aerosols, this reflected so much solar energy back into space that the Northern Hemisphere cooled by 1° C, making 1816 ‘the year without a summer’. Crop failures in Europe and North America doubled grain prices, leading to widespread social unrest and economic depression. That year also saw unusual weather in India culminate in a cholera outbreak, which spread to unleash the 1817 global pandemic. Tambora is implicated in a global death toll in the tens of millions. Thanks to the record of sulfur in Greenland ice cores it has proved possible to link past volcanic action to historic famines and epidemics, such as the Plague of Justinian in 541 CE. If they emit large amounts of sulfur gases volcanic eruptions can result in sudden global climatic downturns.

The ash plume towering above Pinatubo volcano in the Philippines on 12 June 1991, which rose to 40 km (Credit: Karin Jackson U.S. Air Force)

With this in mind Markus Stoffel, Christophe Corona and Scott St. George of the University of Geneva, Switzerland, CNRS, Grenoble France and global insurance brokers WTW, London, respectively, have published a Comment in Nature warning of this kind of global hazard (Stoffel, M., Corona, C. & St. George, S. 2024.  The next massive volcano eruption will cause climate chaos — we are unprepared. Nature v. 635, p. 286-289; DOI: 10.1038/d41586-024-03680-z). The crux of their argument is that there has been nothing approaching the scale of Tambora for the last two centuries. The 1991 eruption of Pinatubo fed the stratosphere with just over a quarter of Tambora’s complement of SO2, and decreased global temperatures by around 0.6°C during 1991-2. Should one so-called Decade Volcanoes – those located in densely populated areas, such as Vesuvius – erupt within the next five years actuaries at Lloyd’s of London estimate economic impacts of US$ 3 trillion in the first year and US$1.5 trillion over the following years. But that is based on just the local risk of ash falls, lava and pyroclastic flows, mud slides and lateral collapse, not global climatic effects. So, a Tambora-sized or larger event is not countenanced by the world’s most famous insurance underwriter: probably because its economic impact is incalculable. Yet the chances of such a repeat certainly are conceivable. A 60 ka record of sulfate in the Greenland ice cores allows the probability of eruptions on the scale of Tambora to be estimated. The data suggest that there is a one-in-six chance that one will occur somewhere during the 21st century, but not necessarily at a site judged by volcanologists to be precarious . Nobody expected the eruption from the Pacific Ocean floor of the Hunga Tonga-Hunga Ha’apai volcano on January 15, 2022: the largest in the last 30 years.

The authors insist that climate-changing eruptions now need to be viewed in the context of anthropogenic global warming. Superficially, it might seem that a few volcanic winters and years without a summer could be a welcome, albeit short-term, solution. However, Stoffel, Corona and St. George suggest that the interaction of a volcano-induced global cooling with climatic processes would probably be very complex. Global warming heats the lower atmosphere and cools the stratosphere. Such steady changes will affect the height to which explosive volcanic plumes may reach. Atmospheric circulation patterns are changing dramatically as the weather of 2024 seems to show. The same may be said for ocean currents that are changing as sea-surface temperatures increase. Superimposing volcano-induced cooling of the sea surface adds an element of chaos to what is already worrying. What if a volcanic winter coincided with an el Niño event? The Intergovernmental Panel on Climate Change that projects climate changes is ‘flying blind’ as regards volcanic cooling. Another issue is that our knowledge of the effects in 1815 of Tambora concerned a very different world from ours: a global population then that was eight times smaller than now; very different patterns of agriculture and habitation; a world with industrial production on a tiny proportion of the continental surface. Stoffel, Corona and St. George urge the IPCC to shed light on this major blind spot. Climate modellers need to explore the truly worst-case scenarios since a massive volcanic eruption is bound to happen one day. Unlike global warming from greenhouse-gas emission, there is absolutely nothing that can be done to avert another Tambora.

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.

Repeated climate and ecological stress during the run-up to the K-Pg extinction

The Cretaceous-Palaeogene mass extinction is no longer an event that polarises geologists’ views between a slow volcanic driver (The Deccan large igneous province) and a near instantaneous asteroid impact (Chicxulub). There is now a broad consensus that both processes were involved in weakening the Late Cretaceous biosphere and snuffing out much of it around 66 Ma ago. Yet is still no closure as regards the details. From a palaeontologist’s standpoint the die-off varied dramatically between major groups of animals. For instance, the non-avian dinosaurs disappeared completely while those that evolved to modern birds did not. Crocodiles came through it largely unscathed unlike aquatic dinosaurs. In the seas those animals that lived in the water column, such as ammonites, were far more affected than were denizens of the seafloor. But much the same final devastation was visited on every continent and ocean. However, lesser and more restricted extinctions occurred before the Chicxulub impact.

Scientists from Norway, Canada, the US, Italy, the UK and Sweden have now thrown light on the possibility that climate change during the last half-million years of the Cretaceous may have been eroding biodiversity and disrupting ecosystems (Callegaro, S. et al. 2023. Recurring volcanic winters during the latest Cretaceous: Sulfur and fluorine budgets of Deccan Traps lavas. Science Advances, v. 9, article eadg8284; DOI: 10.1126/sciadv.adg8284). Almost inevitably, they turned to the record of Deccan volcanism that overlapped the K-Pg event, specifically the likely composition of the gases that the magmas may have belched into the atmosphere. Instead of choosing the usual suspect carbon dioxide and its greenhouse effect, their focus was on sulfur and fluorine dissolved in pyroxene grains from 15 basalts erupted in the 10 Formations of the Deccan flood-basalt sequence. From these analyses they were able to estimate the amounts of the two elements in the magma erupted in each of these 10 phases.

Exposed section through a small part of the Deccan Traps in the Western Ghats of Maharashtra, India. (Credit: Gerta Keller, Princeton University)

The accompanying image of a famous section through the Deccan Traps SE of Mumbai clearly shows that 15 sampled flows could reveal only a fraction of the magmas’ variability: there are 12 flows in the foreground alone. The mountain beyond shows that the pale-coloured sequence is underlain by many more flows, and the full Deccan sequence is about 3.5 km thick. Clearly, flood-basalt volcanism is in no way continuous, but builds up from repeated lava flows that can be as much as 50 m thick. Each of them is capped by a red, clay-rich soil or bole – from the Greek word bolos (βόλος) meaning ‘clod of earth’. Weathering of basalt would have taken a few centuries to form each bole. Individual Deccan flows extend over enormous areas: one can be traced for 1500 km. At the end of volcanism the pile extended over roughly 1.5 million km2 to reach a volume of half a million km3.

Fluorine is a particularly toxic gas with horrific effects on organisms that ingest it. In the form of hydrofluoric acid (HF) – routinely used to dissolve rock – it penetrates tissue very rapidly to react with calcium in the blood to form calcium fluoride. This causes very severe pain, bone damage and other symptoms of skeletal fluorosis. The 1783-4 eruption of the Laki volcanic fissure in Iceland emitted an estimated 8,000 t of HF gas that wiped out more than half the domestic animals as a result of their eating contaminated grass. The famine that followed the eruption killed 20 to 25% of Iceland’s people: exhumed human skeletons buried in the aftermath show the distinctive signs of endemic skeletal fluorosis. This small flood-basalt event had global repercussions, as the Wikipedia entry for Laki documents. Volcanic sulfur emissions in the form of SO2 gas react with water vapour to form sulphuric acid aerosols in a reflective haze. If this takes place in the stratosphere as a result of powerful eruptions, as was the case with the 1991 Pinatubo eruption in the Philippines, the high-altitude haze lingers and spreads. This results in reduced solar warming: a so-called ‘volcanic winter’. In the Pinatubo aftermath global temperatures fell by about 0.5°C during 1991-3. Unsurprisingly, volcanic sulfur emissions also result in acid rainfall. Moreover, inhaling the sulphur-rich haze at low altitudes causes victims to choke as their respiratory tissues swell: an estimated 23,000 people in Britain died in this way when the 1783-4 Laki eruption haze spread southwards Sara Calegaro and colleagues found that the fluorine and sulfur contents of Deccan magmas fluctuated significantly during the eruptive phases. They suggest that fluorine emissions were far above those from Laki, perhaps leading to regional fluorine toxicity around the site of the Deccan flood volcanism but not extinctions. Global cooling due to sulphuric acid aerosols in the stratosphere is suggested to have happened repeatedly, albeit briefly, as eruption waxed and waned during each phase. Magmas rich in volatiles would have been more likely to erupt explosively to inject SO2 to stratospheric altitudes (above 10 to 20 km). The authors do not attempt to model when such cooling episodes may have occurred: data from only 15 levels in the Deccan Traps do not have the time-resolution to achieve that. They do, however, show that this large igneous province definitely had the potential to generate ‘volcanic winters’ and toxic episodes. Time and time again ecosystems globally and regionally would have experienced severe stress, the most important perhaps being disruption of the terrestrial and marine food chains.

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

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

Did giant impacts trigger formation of the bulk of continental crust?

Earth is the only one of the rocky Inner Planets that has substantial continental crust, the rest being largely basaltic worlds. That explains a lot. For a start, it means that almost 30 percent of its surface area stands well above the average level of the basaltic ocean basins – more than 5 km – because of the difference in density between continental and oceanic lithosphere. Without continents and the inability of subduction to draw them back  into the mantle  Earth would remain a water-world as it is thought to have been during the Hadean and early Archaean Eons. The complex processes involved in geochemical differentiation and the repeated reworking of the continents through continual tectonic and sedimentary processes has further enriched parts of them in all manner of useful elements and chemical compounds. And, of course, the land has had a huge biosphere since the Devonian period that subsequently helped to draw down CO­2 well as evolving us.

It has been estimated that during the Archaean (4.0 to 2.5 Ga) around 75% of continental crust formed. Much of this Archaean crust is made up of sodium-rich granitoids: grey tonalite-trondhjemite-granodiorite (TTG) gneisses in the main. Their patterns of trace elements strongly suggest that their parent magmas formed by partial melting at shallow depths (25 to 50 km). Their source was probably basalts altered by hydrothermal fluids to amphibolites, unlike the post-Archaean dominance of melting associated with subducted slabs of lithosphere. Yet most of the discourse on early continents has centred on when plate tectonics began and when they became strong enough to avoid disruption into subductible ‘chunks’. Yet 10 years ago geochemists at the University of St Andrews in Scotland used hafnium and oxygen isotopes in Archaean zircons to suggest that the first continents grew very quickly in the Hadean and early Archaean at around 3.0 km3 yr-1, slowing to an average of 0.8 km3 yr-1 after 3.5 Ga. In 2017 Geochemists working on one of the oldest cratons in the Pilbara region of Western Australia developed a new, multistage model for early crust formation that did not have a subduction component. They proposed that high degrees of mantle melting first produced a mafic-ultramafic crust of komatiites, which became the source for a 3.5 Ga mafic magma with a geochemistry similar to those of modern island-arc basalts. If a crust of that composition attained a thickness greater than 25 km and was itself partially melted at its base, theoretically it could have generated TTG magma and Archaean continental crust. Three members of that team from Curtin University, Western Australia, and others have now contributed to formulating a new possibility for early continent formation (Johnson, T.E. et al. 2022.  Giant impacts and the origin and evolution of continents. Nature, v. 608, p. 330–335; DOI: 10.1038/s41586-022-04956-y).

The distinctive Archaean granite-greenstone terrain of the Pilbara craton of Western Australia. TTG granites are shown in reds in the form of domes, which are enveloped by metamorphosed sediments and mafic-ultramafic volcanics in khaki and emerald green. Other colours signify post Archaean rocks. (Credit: Warren B. Hamilton; Earth’s first two billion years. GSA, 2007)

Tim Johnson and colleagues base their views on oxygen isotopes in Archaean zircon grains from the Pilbara. The zircons’ O-isotopes fall into three kinds of cluster: low 18O that indicate a hydrothermally altered source; intermediate 18O suggesting a mantle source; high 18O signifying contamination by metasedimentary and volcanic rocks. The first two alternate in the 3.6 to 3.4 Ga period; 4 clusters with mantle connotations occupy the 3.4 to 3.0 Ga range; a cluster with supracrustal contamination follows 3.0 Ga. This record can be reconciled agreeably with the geological and broad geochemical history of the Pilbara craton. But there is another connection: the Late Heavy Bombardment (LHB) recognised on most rocky bodies in the Solar System.

Bodies with much more sluggish internal processes than the Earth have preserved much of their earliest surfaces and the damage they have suffered since the Hadean. The Moon is the best example. Its earliest rocks in the lunar Highlands record a vast number of impact craters. Their relative ages, deduced from older ones being affected by later ones, backed up by radiometric ages of materials produced by impacts, such as melt spherules and basaltic magmas that flooded the lunar maria, revealed the time span of the LHB. The maria formed between 4.2 and 3.2 billion years ago and the damage done then is shown starkly by the dark maria that make up the ‘face’ of the Man in the Moon. The lunar bombardment was at a maximum between 4.1 and 3.8 Ga but continued until 3.5 Ga, dropping off sharply from its maximum effects. Earth preserves no tangible sign of the LHB, but because it is larger and more massive than the Moon, and both have always been in much the same orbit around the Sun, it must have been subject to impacts on a far grander scale. Projectiles carry kinetic energy that enables them to do geological work when they impact: 1/2 x mass x speed2. The minimum speed of an impact is the same as the target’s escape velocity – 2.4 km s-1 for the Moon and 11.2 km s-1 for the Earth. So the energy of an object hitting the Earth would be 20 times more than if it struck the lunar surface. Taking into account the Earth’s larger cross sectional area, the amount of geological work done here by the LHB would have been as much as 300 times greater than that on Earth’s battered satellite.

The Earth’s early geological history was rarely seen in that context before the 21st century, but that is the framework plausibly adopted by Johnson and colleagues. Archaean  sediments in South Africa contain several beds of impact spherules older than 3.2 Ga, as do those of the Pilbara. The LHB also left a geochemical imprint on Earth in the form of anomalous isotope proportions of tungsten in 3.8 Ga gneisses from West Greenland (See: Tungsten and Archaean heavy bombardment and Evidence builds for major impacts in Early Archaean; respectively, July and August 2002). Johnson et al. suggest a 3-stage process for the evolution of the Pilbara craton: First a giant impact akin to the lunar Maria that formed a nucleus of mafic-ultramafic crust from shallow melting of the mantle; its chemical fractionation to produce low-magnesium basalts; and in turn their melting to form TTG magmas and thus a continental nucleus. They conclude:

‘The search for evidence of the Late Heavy Bombardment on Earth has been a long one. However, all along it seems that the evidence was right beneath our feet.’

I agree wholeheartedly, but would add that, until quite recently, many scientists who referred to extraterrestrial influences over Earth history were either pilloried or lampooned by their peers as purveyors of ‘whizz-bang’ science. So, many ‘kept their powder dry’. The weight of evidence and a reversal of wider opinion over the last couple of decades has made such hypotheses acceptable. But it has also opened the door to less plausible notions, such as an impact cause for sudden climate change and even for mythological catastrophes such as the destruction of Sodom and Gomorrah!

See also: Timmer, J. 2022. Did giant impacts start plate tectonics? arsTechnica 11 August 2022.

Signs of massive hydrocarbon burning at the end of the Triassic

One of the ‘Big Five’ mass extinctions occurred at the end of the Triassic Period (~201 Ma), whose magnitude matches that of the more famous end-Cretaceous (K-Pg) event. It roughly coincided with the beginning of break-up of the Pangaea supercontinent that was accompanied by a major episode of volcanism preserved in the Central Atlantic Magmatic Province (CAMP). Eastern North America, West Africa and northern South America reveal scattered patches of CAMP flood basalts, swarms of dykes and large intrusive sills. Like all mass extinctions, that at the Triassic-Jurassic boundary left a huge selection of vacant or depleted ecological niches ready for evolution to fill by later adaptive radiation of surviving organisms. Because it coincided with continental break-up and drift, unlike other such events, evolution proceeded in different ways on the various wandering land masses and in newly formed seas (see  an excellent animation of the formation and break-up of Pangaea – move the slider to 3 minutes for the start of break-up). The Jurassic was a period of explosive evolution among all groups of organisms. The most notable changes were among marine cephalopods, to give rise to a bewildering variety of ammonite species, and on land with the appearance and subsequent diversification of dinosaurs.

Pangaea at the end of the Triassic (top) and in Middle Cretaceous times (Credit: screen shots from animation by Christopher Scotese)

Many scientists have ascribed the origin of these events to the CAMP magmatic activity and the release of huge amounts of methane to trigger rapid global warming. In October 2021 one group focused on a special role for the high percentages of magma that never reached the surface and formed huge intrusions that spread laterally in thick sedimentary sequences to ‘crack’ hydrocarbons to their simplest form, CH4 or methane. A sedimentary origin of the methane, rather than its escape from the mantle, is indicated by the carbon-isotope ‘signature’ of sediments deposited shortly after the Tr-J event. The lighter isotope 12C rose significantly relative to 13C, suggesting an organic source – photosynthesis selectively takes up the lighter isotope.

By examining the element mercury (Hg) in deep ocean sediments from a Tr-J sedimentary section now exposed in Japan, scientists from China, the US and Norway have added detail to the methane-release hypothesis (Shen, J et al. 2022. Mercury evidence for combustion of organic-rich sediments during the end-Triassic crisis. Nature Communications, v. 13, article 1307; DOI:10.1038/s41467-022-28891-8). The relative proportions of Hg isotopes strongly suggest that the mercury had been released, as was the methane, from organic-rich sediments rather than from the CAMP magmas (i.e. ultimately from the mantle) through gasification and then burning at the surface.

The hypothesis is enlivened by a separate study (Fox C.P. et al. 2022. Flame out! End-Triassic mass extinction polycyclic aromatic hydrocarbons reflect more than just fire. Earth and Planetary Science Letters, v. 584, article 117418; DOI: 10.1016/j.epsl.2022.117418) that sees magmatic heating as being not so important. Calum Fox and colleagues at Curtin University, Western Australia analysed sediments from a Triassic-Jurassic sedimentary sequence near the Severn Bridge in SW England, focusing on polycyclic hydrocarbons in them. Their results show little sign of the kinds of organic chemical remnants of modern wildfires. Instead they suggest a greater contribution from soil erosion by acid rain that increased input of plant debris to a late Triassic marine basin

See also: How a major volcanic eruption paved the way for the rise of the dinosaurs Eureka Alert 23 March 2022;  Soil erosion and wildfire: another nail in coffin for Triassic era. Science Daily, 21 March 2022

Pinpointing the source of Martian meteorites and a stab at magmatism on Mars

Most meteorites found on the Earth’s surface are fragments of small bodies left over from the accretion of the planets around 4.5 billion years ago, thanks largely to collisions among larger, asteroid-sized bodies. A minority have other origins: some as debris from otherwise icy comets and a few that have been flung off other rocky planets or large moons by crater-forming impacts. Meteorites suspected to have originated through impact are ‘rocky’ – i.e.  made of silicates – and have textures and mineral contents suggesting they formed late in planetary evolution. Most are igneous with basaltic or ultramafic composition: respectively lavas and cumulates formed in magma chambers. Some are breccias, hinting at a pyroclastic origin. The radiometric ages of such planetary fragments are generally far younger than the times when the solar system and planets formed.  Almost 300 have been classified as coming from Mars, only two of which are older than 1400 Ma. The most numerous group of Martian meteorites, known as shergottites, crystallised between 575 and 150 Ma ago to form crust of igneous origin. During the journey from their source to Earth meteorites are exposed to high-energy cosmic rays that generate a variety of new isotopes, from whose relative proportions their travel time can be estimated. The shergottites all seem to have been blasted from Mars a mere 1.1 Ma ago, suggesting that a single impact launched them. So, identifying their source crater on Mars would enable the shergottites to be treated in the same way as samples collected by geologists from a small locality on Earth. Their geochemistry should give important clues to processes within Mars over a time period that spans the late-Precambrian to early Cretaceous on Earth.

Kuiper crater on the Moon, with rays and secondary craters. (Credit: NASA/Johns Hopkins University, USA)

There are many craters on Mars, so homing-in on a single source for shergottite meteorites might seem a tall order. A strategy for doing that depends on recognising craters formed by impacts with sufficient energy to eject debris at the escape velocity from Martian gravity: about 5 km s-1 compared with 11 km s-1 for Earth. Calculations suggest that such impacts would produce craters larger than 3 km across. Large ejecta travelling at slower speeds from them would fall back to produce smaller craters arranged radially from the main crater, forming distinctive rays. Anthony Lagain and colleagues from Curtin University, Western Australia and other institutions in Australia, USA, France and Côte d’ Ivoire adapted a detection algorithm to locate craters less than 1 km across that formed in rays around larger craters (Lagain, A. and 10 others 2021. The Tharsis mantle source of depleted shergottites revealed by 90 million impact craters. Nature Communications, v. 12, article6352; DOI: 10.1038/s41467-021-26648-3). They used 100 m resolution images of thermal emission from the Martian surface that most clearly distinguish large craters that have ejecta deposits around them. Then they turned to images with 0.25 m resolution covering the visible spectrum that can spot very small craters. The authors’ analysis compiled around 90 million impact craters smaller than 300 metres across (a quarter the size of the celebrated Meteor Crater in Arizona).

Laser-altimetry data that show two large impact craters and their ejecta aprons on the Tharsis Plateau of Mars and two of its huge volcanoes: grey-brown-red-orange-yellow-green = high-to-low elevations. (Credit: NASA / JPL-Caltech / Arizona State University)

Dust storms on Mars gradually fill and obscure small craters and ejecta rays, so the younger the impact event, the more visible are rays and secondary small craters. Luckily, just two large craters on Mars have well-preserved rays that contain high densities of small secondary craters. Both of them lie on the Tharsis Plateau near the Martian Equator. This is a vast bulge on the planet’s surface – 5000 km across and rising to 7 km – characterised by three enormous shield volcanoes that rise to 18 km above the average elevation of Mars. The authors judge that one or the other crater is the source for shergottite meteorites, and that this meteorite class collectively samples the most recent igneous rocks that form the Tharsis Plateau. So vast is its mass, that the plateau has probably built-up over most of Mars’s history. One hypothesis is that the bulging has progressively developed over a huge thermal anomaly that has supported a mantle superplume for billions of years from which basaltic magma has steadily moved to the surface.

This model of a perpetual hot spot beneath Tharsis implies that the magmas that it has generated in the past have progressively depleted the underlying mantle in the incompatible trace elements that preferentially enter magma rather than remaining in solid minerals during partial melting. Having been able to suggest that the 575 to 150 Ma-old shergottites represent the upper crust of Tharsis that formed at that late stage in its history, Lagain et al. use those meteorites’ well-established trace-element geochemistry to test that hypothesis. They do indeed suggest their derivation by partial melting of mantle rocks that had in earlier times been strongly depleted in incompatible elements. One of the greatest mysteries about Mars’ evolution may have been resolved without the need for a crewed mission.

Influence of massive igneous intrusions on end-Triassic mass extinction

About 200 Ma ago, the break-up of the Pangaea supercontinent was imminent. The signs of impending events are spread through the eastern seaboard of North America, West Africa and central and northern South America. Today, they take the form of isolated patches of continental flood basalts, dyke swarms – probably the feeders for much more extensive flood volcanism – and large intrusive sills. Break-up began with the separation of North America from Africa and the start of sea-floor spreading that began to form the Central Atlantic Ocean: hence the name Central Atlantic Magmatic Province (CAMP) for the igneous activity. It all kicked off at the time of the Triassic-Jurassic stratigraphic boundary, and a mass extinction with a similar magnitude to that at the end of the Cretaceous. Disappearances of animals in the oceans and on continents were selective rather than general, as were extinctions of land plants. The mass extinction is estimated to have taken about ten thousand years. It left a great variety of ecological niches ready for re-occupation. On land a small group of reptiles with a substantial destiny entered some of these vacant niches. They evolved explosively to the plethora of later dinosaurs as their descendants became separated as a result of continental drift and adaptive radiation.

Flood basalts of the Central Atlantic Magmatic Province in Morocco (Credit: Andrea Marzoli)

The end-Triassic mass extinction, like three others of the Big Five, was thus closely associated in time with massive continental flood volcanism: indeed one of the largest such events. Within at most 10 ka large theropod dinosaurs entered the early Jurassic scene of eastern North America. The Jurassic was a greenhouse world whose atmosphere had about five times more CO2, a mean global surface temperature between 5 and 10°C higher and deep ocean temperatures 8°C above those at present. Was mantle carbon transported by CAMP magmas the main source (widely assumed until recently) or, as during the end-Permian mass extinction, was buried organic carbon responsible? A multinational group of geoscientists have closely examined samples from a one million cubic kilometre stack of intrusive basaltic sills, dated at 201 Ma, in the Amazon basin of Brazil that amount to about a third of all CAMP magmatism (Capriolo, M. and 11 others 2021. Massive methane fluxing from magma–sediment interaction in the end-Triassic Central Atlantic Magmatic ProvinceNature Communications, v. 12, article 5534; DOI: 10.1038/s41467-021-25510-w).

The team focussed on fluid inclusions in quartz within the basaltic sills that formed during the late stages of their crystallisation. The tiny inclusions contain methane gas and tiny crystals of halite (NaCl) as well as liquid water. Such was the bulk composition of the intrusive magma that the presence of around 5% of quartz in the basalts would be impossible without their magma having assimilated large volumes of silica-rich sedimentary rocks such as shales. The host rocks for the huge slab of igneous sills are sediments of Palaeozoic age: a ready source for contamination by both organic carbon and salt. The presence of methane in the inclusions suggests that more complex hydrocarbons had been ‘cracked’ by thermal metamorphism. Moreover, it is highly unlikely to have been derived from the mantle, partly because methane has been experimentally shown not to be soluble in basaltic magmas whereas CO2 is. The authors conclude that both quartz and methane entered the sills in hydrothermal fluids generated in adjacent sediments. Thermal metamorphism of the sediments would also have driven such fluids to the surface to inject methane directly to the atmosphere. Methane is 25 times as potent as carbon dioxide at trapping heat in the atmosphere, yet it combines with the hydroxyl (OH) radical to form CO2 and water vapour within about 12 years. Nevertheless during continuous emission methane traps 84 times more heat in the atmosphere than would an equivalent mass of carbon dioxide.

Calculations suggest about seven trillion tonnes of methane were generated by the CAMP intrusions in Brazil. Had the magmas mainly been extruded as flood basalts then perhaps global warming at the close of the Triassic would have been far less. Extinctions and subsequent biological evolution would have taken very different paths; dinosaurs may not have exploded onto the terrestrial scene so dramatically during the remaining 185 Ma of the Mesozoic. So it seems important to attempt an explanation of why CAMP magmas in Brazil did not rise to the surface but stayed buried as such stupendous igneous intrusions. Work on smaller intrusive sills suggests that magmas that are denser than the rocks that they pass through – as in a large, thick sedimentary basin – are forced by gravity to take a lateral ‘line of least resistance’ to intrude along sedimentary bedding. That would be aided by the enormous pressure of steam boiled from wet sedimentary rocks forcing beds apart. In areas where only thin sedimentary cover rests on crystalline, more dense igneous and metamorphic rocks, basaltic magma has a greater likelihood of rising through vertical dyke swarms to reach the surface and form lava floods.

Weak lithosphere delayed the formation of continents

There are very few tangible signs that the Earth had continents at the surface before about 4 billion years (Ga) ago. The most cited evidence that they may have existed in the Hadean Eon are zircon grains with radiometric ages up to 4.4 Ga that were recovered from much younger sedimentary rocks in Western Australia. These tiny grains also show isotopic anomalies that support the existence of continental material, i.e. rocks of broadly granitic composition, only 100 Ma after the Earth formed (see: Zircons and early continents no longer to be sneezed at; February 2006). So, how come relics of such early continents have yet to be discovered in the geological record? After all granitic rocks – in the broad sense – which form continents are so less dense than the mantle that modern subduction is incapable of recycling them en masse. Indeed, mantle convection of any type in the hotter Earth of the Hadean seems unlikely to have swallowed continents once they had formed. Perhaps they are hiding in another guise among younger rocks of the continental crust. But, believe me; geologists have been hunting for them, to no avail, in every scrap of existing continental crust since 1971 when gneisses found in West Greenland by New Zealander Vic McGregor turned out to be almost 3.8 Ga old. This set off a grail-quest, which still continues, to negate James Hutton’s ‘No vestige of a beginning …’ concept of geological time.

There is another view. Early continental lithosphere may have returned to the mantle piece by piece by other means. One that has been happening since the Archaean is as debris from surface erosion and its transportation to the ocean floor, thence to be subducted along with denser material of the oceanic lithosphere. Another possibility is that before 4 Ga continental lithosphere had far less strength than characterised it in later times; it may have been continually torn into fragments small enough for viscous drag to defy buoyancy and consign them into the mantle by convective processes. Two things seem to confer strength on continental lithosphere younger than 4 billion years: its depleted surface heat flow and heat-production that stem from low concentrations of radioactive isotopes of uranium, thorium and potassium in the lower crust and sub-continental mantle; bolstering by cratons that form the cores of all major continents. Three geoscientists at Monash University in Victoria, Australia have examined how parts of early convecting mantle may have undergone chemical and thermal differentiation (Capitanio, F.A. et al. 2020. Thermochemical lithosphere differentiation and the origin of cratonic mantle.  Nature, v. 588, p. 89-94; DOI: 10.1038/s41586-020-2976-3). These processes are an inevitable outcome of the tendency for mantle melting to begin as it becomes decompressed when pressure decreases when it rises during convection. Continual removal of the magmas produced in this way would remove not only much of the residue’s heat-producing capacity – U, Th and K preferentially enter silicate melts – but also its content of volatiles, especially water. Even if granitic magmas were completely recycled back to the mantle by the greater vigour of the hot, early Earth, at least some of the residue of partial melting would remain. Its dehydration would increase its viscosity (strength). Over time this would build what eventually became the highly viscous thick mantle roots (tectosphere) on which increasing amounts of the granitic magmas could stabilise to establish the oldest cratons. Over time more and more such cratonised crust would accumulate, becoming increasingly unlikely to be resorbed into the mantle. Although cratons are not zoned in terms of the age of their constituent rocks, they do jumble together several billion years’ worth of continental crust in what used to be called ‘the Basement Complex’.

Development of depleted and viscous sub-continental mantle on the early Earth – a precedes b – TTG signifies tonalite-trondhjemite-granodiorite rocks typical of Archaean cratons (Credit, Capitanio et al.; Fig 5)

Early in this process, heat would have made much of the lithosphere too weak to form rigid plates and the tectonics with which geologists are so familiar from the later parts of Earth’s history. The evolution that Capitanio et al. propose suggests that the earliest rigid plates were capped by Archaean continental crust. That implies subduction of oceanic lithosphere starting at their margins, with intra-oceanic destructive plate margins and island arcs being a later feature of tectonics. It is in the later, Proterozoic Eon that evidence for accretion of arc terranes becomes obvious, plastering their magmatic products onto cratons, further enlarging the continents.

How continental keels and cratons may have formed

There is Byzantine ring to the word craton: hardly surprising as it stems from the Greek kratos meaning ‘might’ or ‘strength’. Yes, the ancient cores of the continents were well named, for they are mighty. Some continents, such as Africa, have several of them: probably relics of very ancient supercontinents that have split and spread again and again. Cratons overlie what are almost literally the ‘keels’ of continents. Unlike other mantle lithosphere beneath continental crust (150 km on average) cratonic lithosphere extends down to 350 km and is rigid. Upper mantle rocks at that depth elsewhere are mechanically weaker and constitute the asthenosphere. Geologists only have evidence from the near-surface on which to base ideas of how cratons formed. Their exposed rocks are always Precambrian in age, from 1.5 to 3.5 billion years old, though in some cases they are covered by a thin veneer of later sedimentary rocks that show little sign of deformation. No cratons formed after the Palaeoproterozoic and they are the main repositories of Archaean rock. Their crust is thicker than elsewhere and dominated at the surface by crystalline rocks of roughly granitic composition. Cratons have the lowest amount of heat flowing out from the Earth’s interior; i.e. heat produced by the decay of long-lived radioactive isotopes of uranium, thorium and potassium. This relative coolness provides an explanation for the rigidity of cratons relative to younger continental lithosphere. Because granitic rocks are well-endowed with heat-producing isotopes, the implication of low heat flow is that the deeper parts of the crust are strongly depleted in them. As a result the deep mantle in cratonic keels is at higher pressure and lower temperature than elsewhere beneath the continental surface. Ideal conditions for the formation of diamonds in mantle rock, so that cratonic keels are their main source – they get to the surface in magma pipes when small amounts of partial melting take place in the lithospheric mantle.

The low heat flow through cratons beckons the idea that the heat-producing elements U, Th and K were at some stage driven from depth. An attractive hypothesis is that they were carried in low-density granitic magmas formed by partial melting of mantle lithosphere during the Precambrian that rose to form continental crust. Yet there is an abundance of younger granite plutons that are associated with thinner continental lithosphere. This seeming paradox suggests different kinds of magmagenesis and tectonics during the early Precambrian. Russian and Australian geoscientists have proposed an ingenious explanation (Perchuk, A.L. et al. 2020. Building cratonic keels in Precambrian plate tectonics. Nature, v. 586, p. 395-401; DOI: 10.1038/s41586-020-2806-7). The key to their hypothesis lies in the 2-layered nature of mantle keels beneath cratons, as revealed by seismic studies. Modelling of the data suggests that the layering resulted from different degrees of partial melting in the upper mantle during Precambrian subduction.

Development of a cratonic keel from melt-depleted lithospheric mantle during early Precambrian subduction. Mantle temperature is 250°C higher than it is today. The oceanic lithosphere being subducted in (a) has become a series of stagnant slabs in (b) (credit: Perchuk et al.; Fig. 2)

Perchuk et al. suggest that high degrees of partial melting of mantle associated with subduction zones produced the bulk of magma that formed the Archaean and Palaeoproterozoic crust. This helps explain large differences between the bulk compositions of ancient and more recent continental crust, which involves less melting. The residue left by high degrees of melting of mantle rock in the early Precambrian would have had a lower density than the rest of the mantle. While older oceanic crust at ancient subduction zones would be transformed to a state denser than the mantle as a whole and thus able to sink, this depleted lithospheric mantle would not. In its hot ductile state following partial melting, this mantle would be ‘peeled’ from the associated oceanic crust to be emplaced below. The figure shows one of several outcomes of a complex magmatic-thermomechanical model ‘driven’ by assumed Archaean conditions in the upper mantle and lithosphere An excellent summary of modern ideas on the start of plate tectonics and evolution of the continents is given by:Hawkesworth, C.J., Cawood, P.A. & Dhuime, B. 2020. The evolution of the continental crust and the onset of plate tectonics. In Topic: The early Earth crust and its formation, Frontiers in Earth Sciences; DOI: 10.3389/feart.2020.00326

The Younger Dryas and volcanic eruptions

The issue of the Younger Dryas (YD) cold ‘hiccup’  between 12.9 to 11.7 thousand years (ka) ago during deglaciation and general warming has been the subject of at least 10 Earth-logs commentaries in the last 15 years (you can check them via the Palaeoclimatology logs). I make no apologies for what might seem to be verging on a personal obsession, because it isn’t. That 1200-year episode is bound up with major human migrations on all the northern continents: it may be more accurate to say ‘retreats’. Cooling to near-glacial climates was astonishingly rapid, on the order of a few decades at most. The YD was a shock, and without it the major human transition from foraging to agriculture might, arguably, have happened more than a millennium before it did. There is ample evidence that at 12.9 ka ocean water in the North Atlantic was freshened by a substantial input of meltwater from the decaying ice sheet on northern North America, which shut down the Gulf Stream (see: Tracking ocean circulation during the last glacial period, April 2005; The Younger Dryas and the Flood, June 2006). Such an event has many supporters. Less popular is that it was caused by some kind of extraterrestrial impact, based on various lines of evidence assembled by what amounts to a single consortium of enthusiasts. Even more ‘outlandish’ is a hypothesis that it all kicked off with radiation from a coincident supernova in the constellation Vela in the Southern sky, which is alleged to have resulted in cosmogenic 14C and 10Be anomalies at 12.9 ka. Another coincidence has been revealed by 12.9 ka-old volcanic ash in a sediment core from a circular volcanogenic lake or maar in Germany (see: Did the Younger Dryas start and end at the same times across Europe? January 2014). Being in a paper that sought to chart climate variations during the YD in a precisely calibrated and continuous core, the implications of that coincidence have not been explored fully, until now.

The Laacher See caldera lake in the recently active Eifel volcanic province in western Germany

A consortium of geochemists from three universities in Texas, USA has worked for some time on cave-floor sediments in Hall’s Cave, Texas as they span the YD. In particular, they sought an independent test of evidence for the highly publicised and controversial causal impact in the form of anomalous concentrations of the highly siderophile elements (HSE) osmium, iridium, platinum, palladium and rhenium (Sun, N. et al. 2020. Volcanic origin for Younger Dryas geochemical anomalies ca. 12,900 cal B.P.. Science Advances, v. 6, article eaax8587; DOI: 10.1126/sciadv.aax8587). There is a small HSE ‘spike’ at the 12.9 ka level, but there are three larger ones that precede it and one at about 11 ka. Two isotopes of the element osmium are often used to check the ultimate source of that element through the 187Os/188Os ratio, as can the relative proportions of the HSE elements compared with those in chondritic meteorites. The presence of spikes other than at the base of the YD does not disprove the extraterrestrial causal hypothesis, but the nature of those that bracket the mini-glacial time span not only casts doubt on it, they suggest a more plausible alternative. The 187Os/188Os data from each spike are ambiguous: they could either have arisen from partial melting of the mantle or from an extraterrestrial impact. But the relative HSE proportions point unerringly to the enriched layers having been inherited from volcanic gas aerosols. Two fit dated major eruptions of  the active volcanoes Mount Saint Helens (13.75 to 13.45 ka) and Glacier Peak (13.71 to 13.41 ka) in the Cascades province of western North America. Two others in the Aleutian and Kuril Arcs are also likely sources. The spike at the base of the YD exactly matches the catastrophic volcanic blast that excavated the Laacher See caldera in the Eifel region of western Germany, which ejected 6.3 km3 of sulfur-rich magma (containing 2 to 150 Mt of sulfur). Volcanic aerosols blasted into the stratosphere then may have dispersed throughout the Northern Hemisphere: a plausible mechanism for climatic cooling.

Sun et al. have not established the Laacher See explosion as the sole cause of the Younger Dryas. However, its coincidence with the shutdown of the Gulf Stream would have added a sudden cooling that may have amplified climatic effects of the disappearance of the North Atlantic’s main source of warm surface water. Effects of the Laacher See explosion may have been a tipping point, but it was one of several potential volcanic injections of highly reflective sulfate aerosols that closely precede and span the YD.

See also: Cooling of Earth caused by eruptions, not meteors (Science Daily, 31 July 2020)

Fossil fuel, mercury and the end-Palaeozoic catastrophe

Siberian flood-basalt flows in the Putorana Plateau, Taymyr Peninsula, Russia. (Credit: Paul Wignall)

The end of the Permian Period (~252 Ma ago) saw the loss of 90% of marine fossil species and 70% of those known from terrestrial sediments: the greatest known extinction in Earth’s history. In their naming of newly discovered life forms, palaeontologists can become quite lyrical. Extinctions, however, really stretch their imagination. They call the Permo-Triassic boundary event ‘The Great Dying’. Why not ‘Permageddon’? Sadly, that was snaffled in the 1980s by an astonishingly short-haired heavy-metal tribute band. Enough bathos … The close of the Palaeozoic left a great many ecological niches to be filled by adaptive radiation during the Triassic and later Mesozoic times. Coinciding with the largest known flood-basalt outpouring – the three million cubic kilometres of Siberian Traps – the P-Tr event seemed to be ‘done and dusted’ after that possible connection was discovered in the mid 1990s. Notwithstanding, the quest for a gigantic, causative impact crater continues (see: Palaeobiology Earth-logs, May, September and October 2004), albeit among a dwindling circle of enthusiasts. The Siberian Traps are suitably vast to snuff the fossil record, for their eruption must have belched all manner of climate-changing gases and dusts into the atmosphere; CO2 to encourage global warming; SO2 and dusts as cooling agents. There is also evidence of a role for geochemical toxicity (see: Nickel, life and the end-Permian extinction, June 2014). The extinctions accompanied not only climate change but also a catastrophic fall in atmospheric oxygen content (see: Homing in on the great end-Permian extinction, April 2003; When rain kick-started evolution, December 2019). Recovery of the biosphere during the early Triassic was exceedingly slow.

Research focussed on the P-Tr boundary eventually uncovered an element of pure chance. Shales in Canada that span the boundary show major, negative δ13C excursions in the carbon-isotope record that coincide with fly ash in the analysed layers. This material is similar in all respects to that emitted from coal-fired power stations (see: Coal and the end-Permian mass extinction, March 2011). The part of Siberia onto which the flood basalts were erupted is rich in Permian coal measures and oil shales that lay close to the surface 252 Ma ago. The coal ash and massive emissions of CO2 may have resulted from their burning by the flood basalt event. Now evidence has emerged that this did indeed happen (Elkins-Tanton, L.T. et al. 2020. Field evidence for coal combustion links the 252 Ma Siberian Traps with global carbon disruption. Geology, v. 48, early publication; DOI: 10.1130/G47365.1).

The US, Canadian and Russian team found large quantities of burnt coal and woody material, and bituminous blobs in 600 m thick volcanic ashes at the base of the Siberian traps themselves. They concluded that the magma chamber from which the flood basalts emerged had incorporated sizeable volumes of the coal measures, leading to their combustion and distillation. This would have released CO2 enriched in light 12C due to isotopic fractionation by biological means, i.e. its δ13C would have been sufficiently negative to affect the carbon locked up in the Canadian P-Tr boundary-layer shales that show the sharp isotopic anomalies. The magnitude of the anomalies suggest that between six to ten thousand billion tons of carbon released as CO2 or methane by interaction of the Siberian Traps with sediments through which their magma passed could have created the global δ13C anomalies. That is about one tenth of the organic carbon originally locked in the Permian coal measures beneath the flood basalts

Another paper whose publication coincided with that by Elkins-Tanton et al. suggests that environmental mercury appears to have followed the same geochemical course as did carbon at the end of the Palaeozoic Era (Dal Corso, J. and 9 others 2020. Permo–Triassic boundary carbon and mercury cycling linked to terrestrial ecosystem collapse. Nature Communications, v. 11, paper 2962; DOI: 10.1038/s41467-020-16725-4). This group, based at Leeds and Oxford Universities, UK and the University of Geosciences in Wuhan, China, base their findings on biogeochemical modelling of the global carbon and mercury cycles at the end of the Permian. Their view is that the coincidence in marine sediments at the P-Tr boundary of a short-lived spike in mercury and an anomaly in its isotopic composition with the depletion in 13C, described earlier, shows an intimate link between mercury and the biological carbon cycle in the oceans at the time. They suggest that this synergy marks ecosystem collapse and derives ‘from a massive oxidation of terrestrial biomass’; i.e. burning of organic material on the land surface. Their modelling hints at huge wildfires in equatorial peatlands but also a role for the Siberian flood-basalt volcanism and the incorporation of coal measures into the Siberian Trap magma chamber.

Earliest plate tectonics tied down?

Papers that ponder the question of when plate tectonics first powered the engine of internal geological processes are sure to get read: tectonics lies at the heart of Earth science. Opinion has swung back and forth from ‘sometime in the Proterozoic’ to ‘since the very birth of the Earth’, which is no surprise. There are simply no rocks that formed during the Hadean Eon of any greater extent than 20 km2. Those occur in the 4.2 billion year (Ga) old Nuvvuagittuq greenstone belt on Hudson Bay, which have been grossly mangled by later events. But there are grains of the sturdy mineral zircon ZrSiO4)  that occur in much younger sedimentary rocks, famously from the Jack Hills of Western Australia, whose ages range back to 4.4 Ga, based on uranium-lead radiometric dating. You can buy zircons from Jack Hills on eBay as a result of a cottage industry that sprang up following news of their great antiquity: that is, if you do a lot of mineral separation from the dust and rock chips that are on offer, and they are very small. Given a laser-fuelled SHRIMP mass spectrometer and a lot of other preparation kit, you could date them. Having gone to that expense, you might as well analyse them chemically using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) to check out their trace-element contents. Geochemist Simon Turner of Macquarie University in Sydney, Australia, and colleagues from Curtin University in Western Australia and Geowissenschaftliches Zentrum Göttingen in Germany, have done all this for 32 newly extracted Jack Hills zircons, whose ages range from 4.3 to 3.3 Ga (Turner, S. et al. 2020. An andesitic source for Jack Hills zircon supports onset of plate tectonics in the HadeanNature Communications, v. 11, article 1241; DOI: 10.1038/s41467-020-14857-1). Then they applied sophisticated geochemical modelling to tease out what kinds of Hadean rock once hosted these grains that were eventually eroded out and transported to come to rest in a much younger sedimentary rock.

Artist’s impression of the old-style hellish Hadean (Credit : Dan Durday, Southwest Research Institute)

Zircons only form duuring the crystallisation of igneous magmas, at around 700°C, the original magma having formed under somewhat hotter conditions – up to 1200°C for mafic compositions. In the course of their crystallising, minerals take in not only the elements of which they are mainly composed, zirconium, silicon and oxygen in the case of zircon , but many other elements that the magma contains in low concentrations. The relative proportions of these trace elements that are partitioned from the magma into the growing mineral grains are more or less constant and unique to that mineral, depending on the particular composition of the magma itself. Using the proportions of these trace elements in the mineral gives a clue to the original bulk composition of the parent magma. The Jack Hills zircons  mainly  reflect an origin in magmas of andesitic composition, intermediate in composition between high-silica granites and basalts that have lower silica contents. Andesitic magmas only form today by partial melting of more mafic rocks under the influence of water-rich fluid driven upwards from subducting oceanic lithosphere. The proportions of trace elements in the zircons could only have formed in this way, according to the authors.

Interestingly, the 4.2 Ga Nuvvuagittuq greenstone belt contains metamorphosed mafic andesites, though any zircons in them have yet to be analysed in the manner used by Turner et al., although they were used to date those late-Hadean rocks. The deep post-Archaean continental crust, broadly speaking, has an andesitic composition, strongly suggesting its generation above subduction zones. Yet that portion of Archaean age is not andesitic on average, but a mixture of three geochemically different rocks. It is referred to as TTG crust from those three rock types (trondhjemite, tonalite and granodiorite). That TTG nature of the most ancient continental crust has encouraged most geochemists to reject the idea of magmatic activity controlled by plate tectonics during the Archaean and, by extension, during the preceding Hadean. What is truly remarkable is that if mafic andesites – such as those implied by the Jack Hills zircons and found in the Nuvvuagittuq greenstone belt – partially melted under high pressures that formed garnet in them, they would have yielded magmas of TTG composition. This, it seems, puts plate tectonics in the frame for the whole of Earth’s evolution since it stabilised several million years after the catastrophic collision that flung off the Moon and completely melted the outer layers of our planet. Up to now, controversy about what kind of planet-wide processes operated then have swung this way and that, often into quite strange scenarios. Turner and colleagues may have opened a new, hopefully more unified, episode of geochemical studies that revisit the early Earth . It could complement the work described in An Early Archaean Waterworld published on Earth-logs earlier in March 2020.

Better dating of Deccan Traps, and the K-Pg event

Predictably, the dialogue between the supporters of the Deccan Trap flood basalts and the Chicxulub impact as triggers that were responsible for the mass extinction at the end of the Mesozoic Era (the K-Pg event) continues. A recent issue of Science contains two new approaches focussing on the timing of flood basalt eruptions in western India relative to the age of the Chicxulub impact. One is based on dating the lavas using zircon U-Pb geochronology (Schoene, B. et al. 2019. U-Pb constraints on pulsed eruption of the Deccan Traps across the end-Cretaceous mass extinction. Science, v. 363, p. 862-866; DOI: 10.1126/science.aau2422), the other using 40Ar/39Ar dating of plagioclase feldspars (Sprain, C.G. et al. 2019. The eruptive tempo of Deccan volcanism in relation to the Cretaceous-Paleogene boundary. Science, v. 363, p. 866-870; DOI: 10.1126/science.aav1446). Both studies were initiated for the same reason: previous dating of the sequence of flows in the Deccan Traps was limited by inadequate sampling of the flow sequence and/or high analytical uncertainties. All that could be said with confidence was that the outpouring of more than a million cubic kilometres of plume-related basaltic magma lasted around a million years (65.5 to 66.5 Ma) that encompassed the sudden extinction event and the possibly implicated Chicxulub impact. The age of the impact, as recorded by its iridium-rich ejecta found in sediments of the Denver Basin in Colorado, has been estimated from zircon U-Pb data at 66.016 ± 0.050 Ma; i.e. with a precision of around 50 thousand years.

407458aa.2
The Deccan Traps in the Western Ghats of India (Credit: Wikipedia)

Because basalts rarely contain sufficient zircons to estimate a U-Pb age of their eruption, Blair Schoene and colleagues collected them from palaeosols or boles that commonly occur between flows and sometimes incorporate volcanic ash. Their data cover 23 boles and a single zircon-bearing basalt. Sprain et al. obtained 40Ar/39Ar ages from 19 flows, which they used to supplement 5 ages obtained by their team in previous studies that used the same analytical methods and 4 palaeosol ages from an earlier paper by Schoene’s group.

The zircon U-Pb data from palaeosols, combined with estimates of magma volumes that contributed to the lava sequence between each dated stratigraphic level, provide a record of the varying rates at which lavas accumulated. The results suggest four distinct periods of high-volume eruption separated by long. periods of relative quiescence. The second such pulse precedes the K-Pg event by up to 100 ka, the extinction and impact occurring in a period of quiescence. A few tens of thousand years after the event Deccan magmatism rose to its maximum intensity. Schoene’s group consider that this supports the notion that both magmatism and bolide impact drove environmental deterioration that culminated in mass extinction.

The Ar-Ar data derived from the basalt flows themselves, seem to tell a significantly different story. A plot of basalt accumulation, similarly derived from dating and stratigraphy, shows little if any sign of major magmatic pulses and periods of quiescence. Instead, Courtney Sprain’s team distinguish an average eruption rate of around 0.4 km3 per year before the K-Pg event and 0.6 km3 per year following it. Yet they observe from climate proxy data that there seems to have been only minor climatic change (about 2 to 3 °C warming) during the period around and after the K-Pg event when some 75% of the lavas flooded out. Yet during the pre-extinction period of slower effusion global temperature rose by 4°C then fell back to pre-eruption levels immediately before the K-Pg event. This odd mismatch between magma production and climate, based on their data, prompts Sprain et al. to speculate on possible shifts in the emission of climate-changing gases during the period Deccan volcanism: warming by carbon dioxide – either from the magma or older carbon-rich sediments heated by it; cooling induced by stratospheric sulfate aerosols formed by volcanogenic SO2 emissions. That would imply a complex scenario of changes in the composition of gas emissions of either type. They suggest that one conceivable trigger for the post-extinction climate shift may have been exhaustion of the magma source’s sulfur-rich volatile content before the Chicxulub impact added enough energy to the Earth system to generate the massive extrusions that followed it. But their view peters out in a demand for ‘better understanding of [the Deccan Traps’] volatile release’.

A curious case of empiricism seeming to resolve the K-Pg conundrum, on the one hand, yet pushing the resolution further off, on the other …

Read more on Palaeobiology and Magmatism