Milankovich precession and the Palaeocene-Eocene Thermal Maximum

About 56 Ma ago there occurred some of the most dramatic biological changes since the mass extinction at the Cretaceous-Palaeogene boundary. They included rapid expansion and diversification of mammals and land plants, and a plunge in the number of deep-water foraminifera. Global cooling from the Cretaceous hothouse was rudely reversed by sudden global warming of about 5 to 10°C. Some climatologists have ascribed bugbear status to the Palaeocene-Eocene Thermal Maximum (PETM) as a possible scenario for future anthropogenic global warming. The widely accepted cause is a massive blurt into the Palaeocene atmosphere of greenhouse gases, but what caused it is enthusiastically debated. The climate shift is associated with a sudden decrease in the proportion of 13C in marine sediments: a negative spike in δ13C. Because photosynthesis favours the lighter 12C, organic matter has a low δ13C, so a great deal of buried organic carbon may have escaped from the ocean floor, most likely in the form of methane gas. However, massive burning of living terrestrial biomass would produce the same carbon-isotope signal, but absence of evidence for mass conflagration supports methane release. Methane is temporarily held in marine sediments in the form of gas hydrate (clathrate), an ice-like solid that forms at low temperatures on the deep seafloor. Warming of deep sea water or a decrease in pressure, if sea level falls, destabilise clathrates thereby releasing methane gas: the ‘clathrate gun hypothesis’. The main issue is what mechanism may have pulled the trigger for a monstrous methane release.

Massive leak of natural gas – mainly methane – off Sweden in the Baltic Sea, from the probably sabotaged Nord Stream pipeline. (Source: Swedish coastguard agency)

Many have favoured a major igneous event. Between 55.0 and 55.8 Ma basaltic magmatism– continuing today in Iceland – formed the North Atlantic Igneous Province. It involved large-scale intrusion of sills as well as outpourings of flood basalts and coincided with the initial rifting of Greenland from northern Europe (see: Smoking gun for end-Palaeocene warming: an igneous connection; July/August 2004). The occurrence of impact ejecta in end-Palaeocene sediments off the east coast of the US has spawned an extraterrestrial hypothesis for the warming, which could account for the negative spike in δ13C as the product of a burning terrestrial biosphere (see: Impact linked to the Palaeocene-Eocene boundary event; October 2016). Less headline-grabbing is the possibility that the event was part and parcel of the Milankovich effect: an inevitability in the complex interplay between the three astronomical components that affect Earth’s orbital and rotational behaviour: eccentricity, axial tilt and precession. A group of geoscientists from China and the US, led by Mingsong Li of Peking University, have investigated in minute detail the ups and downs of δ13C around 56 Ma in drill cores recovered from a sequence of Palaeocene and Eocene continental-shelf sediments in Maryland, USA (Li, M., Bralower, T.J. et al. 2022. Astrochronology of the Paleocene-Eocene Thermal Maximum on the Atlantic Coastal Plain. Nature Communications, v. 13, Article 5618; DOI: 10.1038/s41467-022-33390-x).

The study involved sampling sediment for carbon- and oxygen-isotope analysis at depth intervals between 3 and 10 cm over a 35 m section through the lower Eocene and uppermost Palaeocene. Calcium abundances in the core were logged at a resolution of 5 mm using an X-ray fluorescence instrument. The results link to variations in CaCO3 in the sediments across the PETM event. Another dataset involves semi-continuous measurements of magnetic susceptibility (MS) along the core. These measurements are able to indicate variations in delivery to the ocean of dissolved calcium and detrital magnetic minerals as climate and continental weathering vary through time. They are widely known to be good recorders of Milankovich cycles. After processing, the Ca and MS data sets show cyclical fluctuations relative to depth within the cores. ‘Tuning’ their frequencies to the familiar time series of Milankovich astronomical climate forcing reveals a close match to what would be expected if the climate fluctuations were paced by the 26 ka axial precession signal. My post of 17 June 2022 about the influence of precession over ‘iceberg armadas’ during the Pleistocene might be useful to re-read in this context. This correlation enabled the researchers to convert depth in the cores to time, so that the timing of fluctuations in carbon- and oxygen-isotope data that the PETM had created could be considered against various hypotheses for its cause. The ‘excursions’ of both began at the same time and reached the maxima of their changes from Palaeocene values over about 6,000 years. The authors consider that is far too long to countenance the release of methane as a result of asteroidal impact, or by massive burning of terrestrial vegetation. The other option that the beginning of the North Atlantic Igneous Province had been the trigger may also be ruled out on two grounds: the magmatism began earlier, and it continued for far longer. The onset of the PETM coincides with an extreme in precession-related climatic forcing. So Li et al. consider that a quirk in the Milankovich Effect could have played a role in triggering massive methane release. This might also explain features of the global calcium record in seafloor sediments as results of a brief period of ocean acidification during the PETM. Such an event would play havoc with carbonate-secreting organisms, such as foraminifera, by lowering the dissolved carbonate ion content on which they depend for their shells: hence their suffering considerable extinction. Of course, the other elements of astronomical forcing – eccentricity and axial tilt – would also have been operating on global climate at the time.  The long-term 100 and 405 ka eccentricity cycles may have played a role in amplifying warming, which may have resulted in increased burial of organic carbon and thus the amount of methane buried beneath the seabed.

A Lower Jurassic environmental crisis

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

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

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

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

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

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

See also:

Climate out of control after the Permian-Triassic mass extinction

The snuffing out of up to 90 percent of all terrestrial and marine species at the end of the Permian (252 Ma) was the outcome of lethal climatic warming. It probably stemmed from a stupendous episode of flood basalt volcanism and intrusions in what is now Siberia that burned vast amounts of peat or coal in the basin that the flows filled (see: Coal and the end-Permian mass extinction; March 2011). The carbon dioxide so released created planetary hyperthermia and toxic acid rain. For at least five million years Earth was an almost sterile world, a notable absence being dense vegetation on the land surface – the Early Triassic is devoid of coal, whereas there is plenty of Late Permian age. Much the same slow recovery of life is found in meagre collections of land and marine animal fossils of that age. Yet, other mass extinctions were followed by recovery and species diversification at a much faster pace.

One conceivable explanation could be the near absence of vegetation whose photosynthesis and burial would otherwise draw down CO2 and the same goes for its marine equivalent phytoplankton. But there is a powerful inorganic means of carbon sequestration: silicate weathering. The chemistry depends on carbon dioxide dissolved in water. For simple silicates it can be expressed as:

2CO2 + H2O + CaSiO3 → Ca2+ + 2HCO3 + SiO2.

The higher the ambient temperature, the faster such reactions proceed. Most silicates are more complex and many common ones, such as feldspars, include aluminium, so that another product of weathering is insoluble, fine-grained clay minerals. So various soluble metal ions (Ca, Mg, K, Na etc), dissolved bicarbonate ions, silica in various guises and clays eventually end up in the sea. Once there, it is possible for them to recombine, as for instance calcium and bicarbonate ions:

Ca2+ + 2HCO3→ CaCO3 + CO2 + H2O

Despite some CO2 gas being released, this reaction results in a net sequestration of carbon in calcium carbonate. Incidentally, the same kind of chemical reaction occurs in the soils produced by weathering. The carbonate may cement soils to form a hard crust of caliche or ‘calcrete’. Chemical weathering enhanced by a hot climate, it might seem, should reduce the greenhouse effect quickly: a feedback mechanism that normally stabilises climate. But that did not happen after the P-Tr extinction event, thereby stressing all remaining life forms. A group of scientists at the University of Waikato in New Zealand have developed a possible explanation for this potentially fatal hazard for life on Earth (Isson, T.T. et al. 2022. Marine siliceous ecosystem decline led to sustained anomalous Early Triassic warmth. Nature Communications, v. 13, article 3509; DOI: 10.1038/s41467-022-31128-3). It focuses on the silica (SiO2) released by chemical weathering, which enters the ocean in the form of a colloid: Si(OH)4, a form of silicic acid known as ‘reactive silica’. Under ‘normal’ conditions, this is removed by organisms, such as diatoms and radiolaria, and is constantly recycled on a time scale of about 400 years, some contributing to deep-ocean oozes in the form of chert. But, like all other marine organisms, they too were victims of the P-Tr mass extinction.

Examples of marine radiolaria (top)

Reactive silica colloids in seawater also participate in inorganic chemical reactions, combining with dissolved metal ions to form complex hydrated aluminosilicates, i.e. more clay minerals. The reactions change the alkalinity of seawater. As a result dissolved HCO3ions transform to CO2 gas and water. Despite the complexity of the chemistry that interweaves the carbon and silicon cycles, there is a simple conclusion. If the abundance of silica-secreting marine organisms falls drastically while continental weathering continues to deliver silica, clay-mineral formation on the ocean floor results in release of CO2 that reverses the effect of enhanced weathering and thus maintains hyperthermal conditions. The other outcome is that less chert and flint granules form Terry Isson and colleagues examined the varying proportion of chert in cores through Lower Triassic marine sediments. A ‘chert gap’characterises the 4 to 6 Ma following the P-Tr boundary event. This can be explained in part by extinction of silica-secreting organisms and by inorganic reactions converting the reactive silica that enhanced weathering delivered to the oceans to clay minerals. This supports the idea that the inorganic part of the silica cycle maintained greenhouse conditions in the absence of organic ‘competition’ for reactive silica. Many other biogeochemical cycles link biological and chemical processes that combine to affect climate: involving phosphorus, nitrogen and iron, to name but three.

A new twist to Pleistocene climate cycles

The combined gravitational pulls of the sun and moon modulate variations in local tidal range. High spring tides occur when the two bodies are opposed at full moon or in roughly the same direction at new Moon. When the positions of sun and moon are at right angles (1st quarter and 3rd quarter) their gravitational pulls partly cancel each other to give neap tides. Consequently, there are two tidal cycles every lunar month.  In a similar way, the varying gravitational pulls of the planets during their orbital cycles impart a repetitive harmony to Earths astronomical behaviour. But their combined effects are on the order of tens of thousand years. Milutin Milankovich (1879-1958), a Serbian engineer, pondered on the possible causes of Earth’s climatic variations, particularly the repetition of ice ages. He was inspired by 19th century astronomers’ suggestion that maybe the gravitational effects of other planets might be a fruitful line of research. Milankovich focussed on how the shape of Earth’s orbit, the tilt of its rotational axis and the way the axis wobbles like that of a spinning top affect the amount of solar heating at all points on the surface: the effects of varying eccentricity, obliquity and precession, respectively.

 Earlier astronomers had calculated cycles of gravitational effects on Earth of the orbits of Jupiter and Saturn of the three attributes of Earth’s astronomical behaviour and found periods of about 100, 41 and 23 thousand years (ka) respectively. The other 3 inner planets and the much more distant giants Uranus and Neptune also have gravitational effects on Earth, but they are negligible compared with those of the two nearest giant planets, because gravitation force varies with mass and inversely with the square of distance. Sadly, Milankovich was long dead when his hypothesis of astronomical climate forcing was verified in 1976 by frequency analysis of the record of oxygen isotopes in foraminifera found in two ocean sediment core from the Southern Indian Ocean. It revealed that all three periods interfered in complex ways during the Late Pleistocene, to dominate variations in sea-surface temperatures and the fluctuating volume of continental ice sheets for which δ18O is a proxy (see: Odds and ends about Milankovich and climate change; February 2017).

Precession of the axis of a spinning top and that of the Earth. At present the northern end of Earth’s axis points to what we now call the Pole Star. Around 11.5 ka from now it will point to the star Vega

This was as revolutionary for climatology as plate tectonics was for geology. We now know that in the early Pleistocene glacial-interglacial cycles were in lockstep with the 41 ka period of axial obliquity, and since 700 ka followed closely – but not perfectly – the 100 ka orbital eccentricity forcing. The transitional period between 1.25 and 0.7 Ma (the Mid-Pleistocene Transition or MPT) suggested neither one nor the other. Milankovich established that axial tilt variations have the greatest influence on solar heating, so the early 41 ka cycles were no surprise. But the dominance of orbital eccentricity on the last 700 ka certainly presented a puzzle, for it has by far the weakest influence on solar heating: 10 times less than those of axial obliquity and precession. The other oddity concerns the actual effect of axial precession on climate change. There are no obvious 23 ka cycles in the climate record, despite the precession signal being clear in frequency analysis and its effect on solar heating being almost as powerful as obliquity and ten times greater than that of orbital eccentricity. Precessional wobbling of the axis controls the time of year when one hemisphere or the other is closest to the Sun. At one extreme it will be the Northern and 11.5 ka later it will be the Southern. The times of solstices and equinoxes also change relative to the calendar that we use today.

There is an important, if obvious, point about astronomical forcing of climate. It is always there, with much the same complicated interactions between the factors: human activities have absolutely no bearing on them. Climatic ‘surprises’ are likely to continue!

Changes in ice-rafted debris (IRD) since 1.7 Ma in a sediment core from the North Atlantic (orange fill) compared with its oxygen-isotope (δ18O) record of changes in continental ice cover (blue fill). At the top are the modelled variations in 23 ka axial precession (lilac) and 41 ka obliquity (green). The red circles mark major interglacial episodes, blue diamonds show the onset of significant ice rafting and orange diamonds are terminations of ice-rafting (TIR). (Credit: Barker et al., Fig. 2)

Sea temperature and ice-sheet volume are not the only things that changed during the Pleistocene. Another kind of record from oceanic sediments concerns the varying proportion in the muddy layers of abnormally coarse sand grains and even small pebbles that have been carried by icebergs; they are known as ice-rafted debris (IRD). The North Atlantic Ocean floor has plenty of evidence for them appearing and disappearing on a layer-by-layer basis. They were first recognised in 1988 by an oceanographer called Helmut Heinrich, who proposed that six major layers rich in IRD in North Atlantic cores bear witness to iceberg ‘armadas’ launched by collapse, or ablation, at the front of surging ice sheets on Scandinavia, Greenland and eastern Canada. Heinrich events, along with Dansgaard-Oeschger events (rapid climatic warming followed by slower cooling) in the progression to the last glacial maximum have been ascribed to a variety of processes  operating on a ‘millennial’ scale. However, ocean-floor sediment cores are full of lesser fluctuations in IRD, back to at least 1.7 Ma ago. That record offers a better chance of explaining fluctuations in ice-sheet ablation. A joint European-US group has investigated their potential over the last decade or so (Barker, S. et al. 2022. Persistent influence of precession on northern ice sheet variability since the early Pleistocene. Science, v. 376, p. 961-967; DOI: 10.1126/science.abm4033). The authors noted that in each glacial cycle since 1.7 Ma the start of ice rafting consistently occurred during a time of decreasing axial obliquity. Yet the largest ablation events were linked to minima in the precession cycles. In the last 700 ka, such extreme events are associated with the terminations of each ice age.

In the earlier part of the record, the 41 ka obliquity ‘signal’ was sufficient to drive glacial-interglacial cycles, hence their much greater regularity and symmetry than those that followed the Mid-Pleistocene Transition. The earlier ice sheets in the Northern Hemisphere also had consistently smaller extents than those after the MPT. Although the records show a role for precession in pre-MPT times in the form of ice-rafting events, the lesser effect of precession on summer warming at higher latitudes, compared with that of axial obliquity, gave it no decisive influence. After 700 ka the northern ice sheets extended much further south – as far as 40°N in North America – where summer warming would always have been commensurately greater than at high northern latitudes. So they were more susceptible to melting during the increased summer warming driven by the precession cycles. When maximum summer heating induced by axial precession in the Northern Hemisphere coincided with that of obliquity the ice sheets as a whole would have become prone to catastrophic collapse.

It is hard to say whether these revelations have a bearing on future climate. Of course, astronomical forcing will continue relentlessly, irrespective of anthropogenic greenhouse gas emissions. Earth has been in an interglacial for the last 11.5 ka, since the Younger Dryas; i.e. about half a precession cycle ago. The combination of obliquity- and precession-driven influences suggest that climate should be cooling and has been since 6,000 years ago, until the Industrial Revolution intervened. Can the gravitational pull of the giant planets prevent a runaway greenhouse effect, or will human effects defy astronomical forces that continually distort Earth’s astronomical behaviour?

Climate and tectonics since 250 Ma

A central feature of the Earth’s climate system is the way that carbon bound in two gases – carbon dioxide (CO2) and methane (CH4) – controls the amount of incoming solar energy that is retained by the atmosphere. Indeed, without one or the other our home world would have been locked in frigidity since shortly after its formation: a sterile, ice-covered planet. The ‘greenhouse effect’ has been ever-present because the material from which the Earth accreted contained carbon as well as every other chemical element from hydrogen to uranium. Naturally reactive, it readily combines with hydrogen and oxygen to form methane and carbon dioxide, which would have escaped the inner Earth as gases to enter the earliest atmosphere as a ‘comfort blanket’, along with water vapour, another greenhouse gas.  Their combined effects have remained crudely balanced so that neither inescapable frigidity nor surface temperatures high enough to boil-off the oceans have ever occurred in the last 4.5 billion years. Earth has remained like the wee bear’s porridge in the Goldilocks story! Even so, global climate has fluctuated again and again from that akin to a steamy greenhouse, through long periods of moderation to extensive glacial conditions, including three that extended from pole-to-pole – ‘Snowball’ Earths –  during in the Precambrian. During the Phanerozoic the Earth has entered three long periods of generally low global temperatures, in the Ordovician, the Carboniferous and during the last 2.5 Ma  that allowed polar ice caps and sea-ice to extend a third of the way to the Equator. These were forced back and forth repeatedly by cyclical influences apparently triggered by astronomically controlled changes to Earth’s orbital and rotational parameters – the Milankovich Effect. Anthropogenic emissions of greenhouse gases in vast and increasing amounts now threaten to disrupt natural climate variation, effectively overthrowing the gravitational influences of distant giant planets that have controlled climate changes that shaped our own evolution since the genus Homo first emerged.

Bubbles of air trapped in cores through the ice sheets of Antarctica and Greenland record decreased volumes of land ice as CO2 content increased and the opposite during glacial episodes. Somehow in step with the astronomical forcing the Earth released greenhouse gas to warm the climate and drew it down to bring on cooling. Since all life forms are built from carbon-rich compounds and some extract it from the environment to build carbonate hard parts, climate and life on land and in the oceans are interlinked. In fact life and death are involved, because once dead organisms and their hard parts are buried before being oxidised in sediments on land, as in peat and ultimately coal, and on the ocean floors as limestones or carbonaceous mudstones, atmospheric carbon is sequestered. Exposed to acid water containing dissolved CO2 from the atmosphere or to oxygen, respectively, the two forms of carbon in solid form are released as greenhouse gas once more. Both take place when sedimentary deposits are exhumed as a result of erosion and tectonics. Another factor is the abundance of available nutrients, themselves released and distributed by erosion and agents of transportation. At present surface waters of the most distant parts of the oceans contains plenty of such nutrients, except for a vital one, dissolved iron. So they are wet ‘deserts’. It seems that during the much dustier times of glacial episodes iron in fine form reached far out into the world’s oceans so that phytoplankton at the base of the food chain ‘bloomed ‘and so did planktonic animals. Dead organisms ‘rained’ to the ocean floor so drawing down CO2 from the atmosphere and decreasing the greenhouse effect. The surface parts of the carbon and rock cycles are extremely complex and climatologists have yet to come to grips with modelling its future climates convincingly. Yet the carbon cycle and much deeper parts of the rock cycle are interwoven too.

Carbon in sedimentary rock can be heated by burial, and some can be subducted to great depths at destructive plate margins together. The same applies to in ocean-floor basalts that have been permeated by circulating sea water through hydrothermal circulation to form carbonates in the altered volcanic rock. In both cases carbon stored for hundreds of million years can be released by metamorphism in orogenic belts at zones of continental collision and deep below island arcs. Carbon from mantle depths that has never ‘seen the light of day’ is also added to the atmosphere when magmas form below oceanic constructive margins, hot spots and subduction zones, and where magmas flood the continental surface. Consequently, plate tectonics and deep mantle convection have surely played a long-term role in the evolution of our planet’s climate system. Geoscientists based in Australia and the UK have used geochemical data to reconstruct the stores of carbon in oceanic plates and thermodynamic modelling to track what may have happened to it and the climate through the last 250 Ma (Müller, R.D. et al. 2022. Evolution of Earth’s tectonic carbon conveyor belt. Nature, v. 605, p. 629-639; DOI: 10.1038/s41586-022-04420-x). Their review is an important step in understanding what underpins climate on a geological time scale, onto which much shorter-term surface influences are superimposed.

The amount of carbon being outgassed as CO2 each year along plate boundaries in the early Jurassic (185 Ma) shown in dark purple (low) to yellow (high). Also shown in shades of blue is the accumulation of carbon stored in each square metre of the ocean plates. Plate motions are shown as grey arrows (credit: Müller, R.D. et al. Clip from video in Supplementary Information)

At mid-ocean ridges basaltic magma wells up from mantle depths and loses much of its content of dissolved CO2. The annual outgassing at ridges, which depends on the global rate of plate formation, has varied from 13 to 30 million tonnes of carbon  (MtC yr-1) since the start of the Mesozoic Era 250 Ma ago. Similarly, there is greenhouse-gas escape from volcanic arcs above subduction zones, estimated to have ranged from 0 to 18 MtC yr-1. As an oceanic plate moves away from its source various processes sequester CO2 into the oceanic crust and upper mantle through accumulation of deep-sea sediments and hydrothermal alteration of basaltic crust and peridotite mantle (ranging from 30 to 311 MtC yr-1). Of this influx of carbon into oceanic plates between 36 to 103 MtC yr-1 has gone down subduction zones in descending slabs. Between 0 to 49 MtC yr-1 of that has been outgassed by arc volcanic activity or absorbed into the overriding plate. The rest continues down into the deep mantle, perhaps to form diamonds. Overall, when the rate at which oceanic plates grow is rapid and plate motion speeds up, outgassing should be high. When plate growth slows, so does the rate of CO2 release. Variations in plate growth can be estimated from the magnetic reversal stripes above the ocean floors.  The authors have released an animation of the break-up of Pangaea (well worth watching at full screen – you can skip the ad at the start), with the rate of carbon emission at ridges and volcanic arcs being colour-coded. Also shown is the storage of carbon within oceanic plats plates as time passes.

Length of mid-ocean ridges (orange) and subduction zones (blue) through the last 250 Ma (top). The areas of oceanic crust produced at ridges and consumed by subduction (bottom) (credit: Müller, R.D. et al., Figs 1a, 1c)

Before Pangaea began to break up at the end of the Triassic (200 Ma) the total length of mid-ocean ridges was at a minimum of about 40 thousand km. Through the Jurassic it never exceeded 50,000 km, but rose to a maximum of 80,000 km during the Cretaceous then declined slowly to the current length of 60,000 km. Throughout the last 250 Ma the length of subduction zones stayed roughly the same at about 65 thousand km – not always in the same places – although the overall rate of subduction changed in line with the rate of oceanic plate growth  (the volume that is added must be balanced roughly by the amount that returns to the mantle).  Between the end of the Jurassic and the mid-Cretaceous crustal production and destruction doubled, shown by the bottom plot in the figure above. The very fast  movement of plates and an increase in the global length of ridges during Jurassic to mid-Cretaceous times led to a dramatic increase in CO2 outgassing from ridges so that its content in the atmosphere rose as high as 1200 ppm – more than four times that before the Industrial Revolution. That level resulted in global ‘hothouse’ conditions during the Cretaceous. Another factor behind the Cretaceous climate was a decrease in the global complement of mountains. That led to decreases in erosion and the weathering of silicates by acid rain, thus reducing natural sequestration of carbon.

During the Cenozoic (after 65 Ma) declining ridge outgassing was actually outpaced by that associated with subduction, according to the modelling. That is strange, for by around 35 Ma glaciation had begun  on Antarctica as the Earth was cooling, which implies a major, unexpected sink for excess CO2. The most likely way this might have arisen is through increased erosion and silicate weathering on the exposed continents that consumed CO2 faster than tectonics was releasing the gas. The length of continental arcs shows no sign of a major increase during the Cenozoic, which might have accelerated that kind of sequestration, but a variety of proxies for signs of weathering definitely suggests that there was an upsurge. Also there was increased storage of carbon on the deep ocean floor, shown by the video. Increased calcium released by weathering to enter ocean water in solution would allow more planktonic organisms to secrete calcite (CaCO3) skeletons that would then fall to the ocean floor when they died.

There may be more to be discovered in this hugely complex interplay between tectonics and climate. For instance, when the bottom waters of the oceans are oxygenated by deep currents of cold dense seawater sinking from polar regions, carbon in tissues of sunken dead organism is oxidised to release CO2. If bottom waters are anoxic, this organic carbon is preserved in sediments. The authors mention this as something to be considered in their future work on  the ‘tectonic carbon conveyor belt’.

The end of the Carboniferous ‘icehouse’ world

From about 340 to 290 Ma the Earth experienced the longest episode of repeated ice ages of the Phanerozoic. The climate then was similar in many ways to that of the Pleistocene. The South Polar region was then within the Pangaea supercontinent and thus isolated from any warming effect from the surrounding ocean: much the same as modern Antarctica but on a much larger scale. Glaciation extended as far across what became the southern continents and India as did the continental ice sheets of the Northern Hemisphere during Pleistocene glacial maxima. Tropical sedimentary rocks of the time, display evidence for repeated alternations of high and low sea levels that mark cycles of glacial maxima and interglacial episodes akin to those of the Pleistocene. In fact they probably reflect the influence of changes in the Earth’s orbit and geometry of its axis of rotation very similar to those predicted by Milankovich from astronomical factors to explain Pleistocene climatic cycles. At the end of the Carboniferous what was an ‘ice-house’ world changed suddenly to its opposite – ‘greenhouse’ conditions – that persisted through the Mesozoic Era until the later part of the Cenozoic, when Antarctica developed is ice cap and global climate slowly cooled to become extremely cyclical once again.

Sedimentary evidence for global climates 320 Ma ago. As well as the large tracts of glaciogenic sediments, smaller occurrences and examples of polished rock surfaces over which ice had passed show the probable full extent (blue line) of ice sheets across the southern, Gondwana sector of Pangaea (Credit: after Fig 7.3, S104, Earth and Space, ©Open University 2007)

The end of the Carboniferous witnessed the collapse of the vast Equatorial rainforests, which formed the coal deposits that put ‘Carbon’ into the name of the Period. By its end this ecosystem had vanished to result in a minor mass extinction of both flora and fauna. Temperatures rose and aridity set in, to the extent that the latest Carboniferous in the British coalfields is marked by redbeds that presage the spread of desert conditions across the Equatorial parts of Pangaea during the succeeding Permian. A team of researchers based at the University of California at Davis have been studying data pertaining to this sudden change have now published their findings (Chan J. and 17 others 2022. Marine anoxia linked to abrupt global warming during Earth’s penultimate icehouse. Proceedings of the National Academy of Sciences, v. 119, article e2115231119; DOI: 10.1073/pnas.2115231119). They used carbon-, oxygen- and uranium isotopes, together with proxies for changes in atmospheric CO2 concentrations, to model changes in the carbon cycle in the Late Carboniferous of China.

Changes in uranium isotopes within marine carbonates are useful indicators of the amount of oxygen available in ocean water at the sea floor. Between 304 and 303.5 Ma ago oxygen content declined by around 30%, the peak of this anoxia being at 303.7 Ma. This occurred about 100 ka after atmospheric CO2 had risen to ~700 parts per million (ppm) from around 350 ppm in the preceding 300 ka, as marked by several proxies.  The authors suggest that the lower ‘baseline’ for the main greenhouse gas marked an extreme glacial maximum. Changes in the proportions of 18O relative to ‘lighter’ 16O in fossil shells suggest that sea-surface temperatures increased in step with the doubling of the greenhouse effect. At the same time there was a major marine transgression as sea level rose. This would have been accompanied by a massive increase in low density freshwater in surface ocean water derived from melting of Pangaea’s ice cap. The team suggests that the freshened surface layer could not sink to carry oxygen to deeper levels, thereby creating anoxic conditions across an estimated 23% of the global seafloor, and thus toxic ‘death zones’ for marine organisms.

One possibility for this sudden rise of atmospheric CO2 is a massive episode of volcanism, perhaps a large igneous province, but there is scanty evidence for that at the end of the Carboniferous. A coinciding sharp decrease in δ13C  in carbonate shells suggests that the excess carbon dioxide probably had an organic origin. So a more plausible hypothesis is massive burning on the continental surface. In the tropics, the huge coals swamps would have contained vast amounts of peat-like decayed vegetable matter as well as living green vegetation. How might that have caught fire? The peat precursor to Carboniferous coal deposits derived from photosynthesis on an unprecedented, and never repeated, scale during tens of million years of thriving tropical rain forest during that Period. This built up atmospheric oxygen levels to about 35%, compared with about 21% today. Insects, whose maximum size is governed by their ability to take in oxygen through spiracles in their bodies, and by the atmospheric concentration of oxygen, became truly huge during the earlier Carboniferous. The more oxygen in the air, the greater the chance that organic matter will catch fire. In fact wet vegetation can burn if oxygen levels rise above 25%. At the levels reached in the Carboniferous huge wildfires in forests and peatlands would have been inevitable. Evidence that huge fires did occur comes from the amount of charcoal found in Carboniferous coal seams, which reach 70% compared with the 4 to 8 % in more recent coals. They may have been ignited by lightning strikes or even spontaneous combustion if decay of vegetation generated sufficient heat, as sometimes happens today in wet haystacks or garden compost heaps.  But how in a short period around 304 Ma could 9 trillion tons of carbon dioxide be released in this way. The preceding  glacial super-maximum, like glacial maxima of the Pleistocene, may have been accompanied by decreased atmospheric humidity: this would dry out the vast surface peat deposits.

The succeeding Permian is famous for its extensive continental redbeds, and so too those of the Triassic. They are red because sediment grains are coated in the iron oxide hematite (Fe2O3). As on Mars, the redbeds are a vast repository for oxygen sequestered from the atmosphere by the oxidation of dissolved Fe2+ to insoluble Fe3+. This had been going on throughout the Permian, the nett result being that by 250 Ma atmospheric oxygen content has slumped to 16% and remained so low for another 50 million years. Photosynthesis failed to resupply oxygen against this inorganic depletion, and there are few coal deposits of Permian or Triassic age: for about 100 Ma Earth ceased to have green continents.

See also: Carbon, climate change and ocean anoxia in an ancient icehouse world. Science Daily, 2 May 2022. 

‘Smoking gun’ for Younger Dryas trigger refuted

In 2018 airborne ice-penetrating radar over the far northwest of the Greenland revealed an impact crater as large as the extent of Washington DC, USA beneath the Hiawatha Glacier. The ice surrounding it was estimated to be younger than 100 ka. This seemed to offer a measure of support for the controversial hypothesis that an impact may have triggered the start of the millennium-long Younger Dryas episode of frigidity (12.9 to 11.7 ka). This notion had been proposed by a group of scientists who claimed to have found mineralogical and geochemical signs of an asteroid impact at a variety of archaeological sites of roughly this age in North America, Chile and Syria. A new study of the Hiawatha crater by a multinational team, including the original discoverers of the impact structure, has focussed on sediments deposited beyond the edge of the Greenland ice cap by meltwater streams flowing along its base. (Kenny, G.G. et al. 2022. A Late Paleocene age for Greenland’s Hiawatha impact structure. Science Advances, v.8, article eabm2434; DOI: 10.1126/science.eabm2434).

Colour-coded subglacial topography from airborne radar sounding over the Hiawatha Glacier of NW Greenland (Credit: Kjaer et al. 2018; Fig. 1D)

Where meltwater emerges from the Hiawatha Glacier downstream of the crater there are glaciofluvial sands and gravels that began to build up after 2010 when rapid summer melting began, probably due to global warming. As luck would have it, the team found quartz grains that contained distinctive planar features that are characteristic of impact shock. They also found pebbles of glassy impact melts that contain clasts of bedrock, further grains of shocked quartz and tiny needles of plagioclase feldspar that crystallised from the melt. Also present were small grains of the mineral zircon (ZrSiO4), both as pristine crystals in the bedrock clasts and porous, grainy-textured grains showing signs of deformation in the feldspathic melt rock. So, two materials that can be radiometrically dated are available: feldspars suitable for the 40Ar/39Ar method and zircons for uranium-lead (U-Pb) dating. The feldspars proved to be about 58 million years old; i.e. of Late Palaeocene age. The pristine zircon grains from bedrock clasts yielded Palaeoproterozoic U-Pb ages (~1915 Ma), which is the general age of the Precambrian metamorphic basement that underpins northern Greenland. The deformed zircon samples have a very precise U-Pb age of 57.99±0.54 Ma. There seems little doubt that the impact structure beneath the Hiawatha Glacier formed towards the beginning of the Cenozoic Era.

During the Palaeocene, Northern Greenland was experiencing warm conditions and sediments of that age show that it was covered with dense forest. The group that since 2007 has been advocating the influence of an impact over the rapid onset of the Younger Dryas acknowledges that the Hiawatha crater cannot support their view. But they have an alternative: an airburst of an incoming projectile. Although scientists know such phenomena do occur, as one did over the Tunguska area in Siberia on the morning of 30 June 1908. Research on the Tunguska Event has discovered  geochemical traces that may implicate an extraterrestrial object, but coincidentally the area affected is underlain by the giant SIberian Traps large igneous province that arguably might account for geochemical anomalies. Airbursts need to have been observed to have irrefutable recognition. Two posts from October 2021 – A Bronze Age catastrophe: the destruction of Sodom and Gomorrah? and Wide criticism of Sodom airburst hypothesis emerges – suggest that some scientists question the data used repeatedly to infer extraterrestrial events by the team that first suggested an impact origin for the Younger Dryas.

See also: Voosen, P, 2022. Controversial impact crater under Greenland’s ice is surprisingly ancient. Science, v. 375, article adb1944;DOI: 10.1126/science.adb1944

The Mid-Pleistocene Transition: when glacial cycles changed to 100 ka

Before about a million years ago the Earth’s overall climate repeatedly swung from warm to cool roughly every 41 thousand years. This cyclicity is best shown by the variation of oxygen isotopes in sea-floor sediments. That evidence stems from the tendency during evaporation at the ocean surface for isotopically light  oxygen (16O) in seawater to preferentially enter atmospheric water vapour relative to 18O.  During cool episodes more water vapour that falls as snow at high latitudes fails to melt, so that glaciers grow. Continental ice sheets therefore extract and store 16O so that the proportion of the heavier 18O increases in the oceans. This shift shows up in the calcium carbonate (CaCO­3) shells of surface-dwelling organisms whose shells are preserved in sea-floor sediment. When the climate warms, the ice sheets melt and return the excess of 16O back to ocean-surface water, again marked by changed oxygen isotope proportions in plankton shells. The first systematic study of sea-floor oxygen isotopes over time revolutionised ideas about ancient climates in much the same way as sea-floor magnetic stripes revealed the existence of plate tectonics. Both provided incontrovertible explanations for changes observed in the geological record. In the case of oxygen isotopes climatic cyclicity could be linked to changes in the Earth’s orbital and rotational behaviour: the Milankovich Effect.

Glacial-interglacial cycles during the Pleistocene

The 41 ka cycles reflect periodic changes in the angle of the Earth’s rotational axis (obliquity), which have the greatest effect on how much solar heating occurs at high latitudes. However, between about 1200 and 600 ka the fairly regular, moderately intense 41 ka climate cycles shifted to more extreme, complex and longer 100 ka cycles at the ‘Mid-Pleistocene Transition’ (MPT). They crudely match cyclical variations in the shape of Earth’s orbit (eccentricity), but that has by far the least influence over seasonal solar heating. Moreover, modelling of the combined astronomical climate influences through the transition show little, if any, sign of any significant change in external climatic forcing. Thirty years of pondering on this climatic enigma has forced climatologists to wonder if the MPT was due to some sort of change in the surface part of the Earth system itself.

There are means of addressing the general processes at the Earth’s surface and how they may have changed by using other aspects of sea-floor geochemistry (Yehudai, M. and 8 others 2021. Evidence for a Northern Hemispheric trigger of the 100,000-y glacial cyclicity. Proceedings of the National Academy of Sciences, v. 118, article e2020260118; DOI: 10.1073/pnas.2020260118). For instance the ratio between the abundance of the strontium isotope 87Sr to that of 86Sr in marine sediments tells us about the progress of continental weathering around a particular ocean basin. The 87Sr/ 86Sr ratio is higher in rocks making up the bulk of the crystalline continental crust than that in basalts of the oceanic crust. That ratio is currently uniform throughout all ocean water. During the Cenozoic Era the ratio steadily increased in sea-floor sediments, reflecting the continual weathering and erosion of the continents. In the warm Pliocene (5.3 to 2.8 Ma) 87Sr/ 86Sr remained more or less constant, but began increasing again at the start of the Pleistocene with the onset of glaciation in the Northern Hemisphere. At about 1450 ka it began to increase more rapidly suggesting increased weathering, and then settled back to its earlier Pleistocene rate after 1100 ka. Another geochemical contrast between the continental and oceanic crust lies in the degree to which the ratio of two isotopes of neodymium (143Nd/144Nd) in rocks deviates from that in the Earth’s mantle – modelled from meteorite geochemistry – a measure signified by ЄNd. Magmatic rocks and young continental rocks have positive ЄNd values, but going back in time continental crust has increasingly negative ЄNd.

Yehudai et al analysed cores from deep-sea sediments that had been drilled between 41°N and 43°S in the Atlantic Ocean floor. They targeted layers designated as glacial and interglacial from their oxygen isotope geochemistry at different levels in the cores to check how ЄNd varied with time. The broad variations within each core look much the same, although at increasingly negative values from south to north, except in one case. The data from the most northerly Atlantic core show far more negative values of ЄNd, in both glacial and interglacial layers at around 950 ka ago, than do cores further to the south. The authors interpret this anomaly as showing a sudden increase in the amount of very old continental rocks – with highly negative ЄNd – that had become exposed at and ground from the base of the great northern ice sheets of North America, Greenland and Scandinavia. At present, the shield areas where the great ice sheets occurred until about 11 ka are almost entirely crystalline Precambrian basement, including the most ancient rocks that are known. Although broadly speaking the shields now have low relief, they are extremely rugged terrains of knobbly basement outcrops and depressions filled with millions of lakes. In the earlier Cenozoic they were covered by younger sedimentary rocks and soils formed by deep weathering, with less-negative ЄNd values. The authors conclude that around 950 ka that younger cover had largely been removed by glacial every every 41 ka or so since about 2.6 Ma ago, when glaciation of the Northern Hemisphere began.

The surface on which the North American ice sheet moved – typical Canadian Shield.

So what follows from that ЄNd anomaly? Yehudai et al suggest that in earlier Pleistocene times each successive ice sheet rested on soft rock; i.e. their bases were well lubricated. As a result, glaciers quickly reached the coast to break up and melt as icebergs drifted south. Exposure of the deeper, very resistant crystalline basement resulted in much more rugged base, as can be seen in northern Canada and Scandinavia today. Friction at their bases suddenly increased, so that much more ice was able to build up on the great shields surrounding the Arctic Ocean than had previously been possible. Shortly after 950 ka the sea-floor cores also reveal that deep ocean circulation weakened significantly in the following 100 ka. The influence on climate of regular, 41 ka changes in the tilt of the Earth’s rotational axis could therefore not be sustained in the later Pleistocene. The ice sheets could neither melt nor slide into the sea sufficiently quickly; indeed, bigger and more durable ice sheets would reflect away more solar heating than was previously possible as glacial gave way to interglacial. The 41 ka astronomical ‘pacemaker’ still operated, but ineffectually. A new and much more complex climate cyclicity set in. Insofar as climate change became stabilised, an overall ~100 ka pulsation emerged. Whether or not this fortuitously had the same pace as the weak influence of Earth’s changing orbital eccentricity remains to be addressed. The climate system just might be too complicated and sensitive for us ever to tell: it may even have little relevance in a climatically uncertain future.

See also: Why did glacial cycles intensify a million years ago? Science Daily, 8 November 2021.

Multiple impacts set back oxygen build-up in the Archaean

Earth’s present atmosphere contains oxygen because of one form of photosynthesis that processes water and carbon dioxide to make plant carbohydrates, leaving oxygen at a waste product. The photochemical trick that underpins oxygenic photosynthesis seems only to have evolved once. It was incorporated in a simple, single-celled organism or prokaryote, which lacks a cell nucleus but contains the necessary catalyst chlorophyll. Such an organism gave rise to cyanobacteria or blue-green bacteria, which still make a major contribution to replenishing atmospheric oxygen. Chloroplasts that perform the same function in plant cells are so like cyanobacteria that they were almost certainly co-opted during the evolution of a section of nucleus-bearing eukaryotes that became the ancestors of plants. A range of evidence suggests that oxygenic photosynthesis appeared during the Archaean Eon, the most tangible being the presence of stromatolites, which cyanobacteria mats or biofilms form today. These knobbly structures in carbonate sediments extend as far back as 3.5 billion years ago (see: Signs of life in some of the oldest rocks; September 2016). Yet it took a billion years before the first inklings of biogenic oxygen production culminated in the Great Oxygenation Event or GOE (see: Massive event in the Precambrian carbon cycle; January, 2012) at around 2400 Ma. Then, for the first time, oxidised iron in ancient soils turned them red. If oxygen was being produced, albeit in small amounts, in shallow, sunlit Archaean seas, why didn’t it build up in the atmosphere of those times? Geochemical analyses of Archaean sediments do point to trace amounts, with a few ‘whiffs’ of more substantial amounts. But they fall well below those of Meso- and Neoproterozoic and Phanerozoic times. One hypothesis is that Archaean oceans contained dissolved, ferrous iron (Fe2+) – a powerful reducing agent – with which available oxygen reacted to form insoluble ferric iron (Fe3+) oxides and hydroxides that formed banded iron formations (BIFS). The Fe2+ in this hypothesis is attributed to hydrothermal activity in basaltic oceanic crust. There is, however, another possibility for suppression of atmospheric oxygen accumulation in the Archaean and early-Palaeoproterozoic.

Summary of the evolution of atmospheric oxygen and related geological features. The percentage scale is logarithmic with the modern level being100%. Credit Alex Glass, Duke University

Simone Marchi of the Southwest Research Institute of Boulder, CO, USA and colleagues from the US, Austria and Germany suggest that planetary bombardment offers a plausible explanation (Marchi, S. et al 2021. Delayed and variable late Archaean atmospheric oxidation due to high collision rates on Earth. Nature Geoscience, v. 14 advance publication; DOI: 10.1038/s41561-021-00835-9). Over the last 20 years evidence of extraterrestrial impacts has emerged, in the form of thin spherule-bearing layers in Archaean sedimentary strata, probably formed by impacts of objects around 10 km across. So far 35 such layers have been identified from several locations in South Africa and Western Australia. They span the last billion years of the Archaean and the earliest Palaeoproterozoic, although they are not evenly spaced in time. The spherules represent droplets of mainly crustal but some meteoritic rocks that were vaporised by impacts and then condensed as liquid. Meteorites in particular contain reduced elements and compounds, including iron, whose oxidation by would remove free oxygen.

The evidence from spherule beds is supplemented by the team’s new calculations of the likely flux of impactors during the Archaean. These stem from re-evaluation of the lunar cratering record that is used to estimate the number and size of impacts on Earth up to 2.5 Ga ago. This flux amounts to the ‘leftovers’ of the catastrophic period around 4.1 Ga when the giant planets Jupiter and Saturn ran amok before they settled into their present orbits. Their perturbation of gravitational fields in the solar system injected a long-lived supply of potential impactors into the inner solar system, which is recorded by craters on the post-4.1 Ga lunar maria. The calculations suggest that the known spherule layers underestimate the true number of such collisions on Earth. Modelling by Marchi et al., based on the meteorite flux and the oxidation of vaporised materials produced by impacts, plausibly accounts for the delay in atmospheric oxygen build-up.

It is worth bearing in mind, however, that large impacts and their geochemical aftermath are, in a geological sense, instantaneous events widely spaced in time. They may have chemically ‘sucked’ oxygen out of the Archaean and early-Palaeoproterozoic atmosphere. Yet photosynthesising bacteria would have been generating oxygen continuously between such sudden events. The same goes for the supply of reduced ferrous iron and its circulation in the oceans of those times, capable of scavenging available oxygen through simple chemical reactions. In fact we can still observe that in action around ocean-floor hydrothermal vents where a host of reduced elements and compounds are oxidised by dissolved oxygen. The difference is that oxygen is now produced more efficiently on land and in the upper oceans and a less vigorous mantle is adding less iron-rich basalt magma to the crust: the balance has changed. Another issue is that the Great Oxygenation Event terminated the oxygen-starved conditions of the Archaean and Palaeoproterozoic in about 200 million years, despite the vast production of BIFs before and after it happened. The Wikipedia entry for the GOE provides a number of hypotheses for how that termination came about. Interestingly, one idea looks to a shortage of dissolved nickel that is vital for methane generating bacteria: a nickel ‘famine’. A geochemical setback for methanogens would have been a boost for oxygenic photosynthesisers and especially their waste product oxygen: methane quickly reacts with oxygen in the atmosphere to produce CO2 and water. Anomalously high nickel is a ‘signature element’ for meteorite bombardment, though it can be released by hydrothermal alteration of basalt. Had meteoritic nickel been fertilising methane-generating bacteria in the oceans prior to the GOE?

See also: A new Earth bombardment model. Science Daily, 21 October 2021.

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.