Aftershocks of ancient earthquakes

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

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

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

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

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

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

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

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

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

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

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

Direct signs of what caused the Palaeocene-Eocene thermal maximum

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

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

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

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

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

News about when subduction began

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

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

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

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

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

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

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

Typical banded iron formation

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

A Lower Jurassic environmental crisis

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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