A new explanation for the Neoproterozoic Snowball Earth episodes

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

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

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

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

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

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

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

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

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

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

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

Kicking-off planetary Snowball conditions

Untitled-1
Artist’s impression of the glacial maximum of a Snowball Earth event (Source: NASA)

Twice in the Cryogenian Period of the Neoproterozoic, glacial- and sea ice extended from both poles to the Equator, giving ‘Snowball Earth’ conditions. Notable glacial climates in the Phanerozoic – Ordovician, Carboniferous-Permian and Pleistocene – were long-lived but restricted to areas around the poles, so do not qualify as Snowball Earth conditions. It is possible, but less certain, that Snowball Earth conditions also prevailed during the Palaeoproterozoic at around 2.4 to 2.1 billion years ago. This earlier episode roughly coincided with the ‘Great Oxidation Event’, and one explanation for it is that the rise of atmospheric oxygen removed methane, a more powerful greenhouse gas than carbon dioxide, by oxidizing it to CO2 and water. That may well have been a consequence of the evolution of the cyanobacteria, their photosynthesis releasing oxygen to the atmosphere. The Neoproterozoic ‘big freezes’ are associated with rapid changes in the biosphere, most importantly with the rise of metazoan life in the form of the Ediacaran fauna, the precursor to the explosion in animal diversity during the Cambrian. Indeed all major global coolings, restricted as well as global, find echoes in the course of biological evolution. Another interwoven factor is the rock cycle, particularly volcanism and the varying pace of chemical weathering. The first releases CO2 from the mantle, the second helps draw it down from the atmosphere when weak carbonic acid in rainwater rots silicate minerals (see: Can rock weathering halt global warming, July 2020). All such interplays between major and sometimes minor ‘actors’ in the Earth system influence climate and, in turn, climate inevitably affects all the rest. With such complexity it is hardly surprising that there is a plethora of theories about past climate shifts.

As well as a link with fluctuations in the greenhouse effect, climate is influenced by changes in the amount of solar heating, for which there are yet more options to consider. For instance, the increase in Earth’s albedo (reflectivity) that results from ice cover, may lead through a feedback effect to runaway cooling, particularly once ice extends beyond the poorly illuminated poles. Volcanic dust and sulfate aerosols in the stratosphere also increase albedo and the tendency to cooling, as would interplanetary dust. More complexity to befuddle would-be modellers of ancient climates. Yet it is safe to say that, within the maelstrom of contributory factors, the freeze-overs of Snowball conditions must have resulted from our planet passing through some kind of threshold in the Earth System. Two theoretical scientists from the Department of Earth, Atmospheric, and Planetary Sciences at the Massachusetts Institute of Technology have attempted to cut through the log-jam by modelling the dynamics of the interplay between the ice-albedo feedback and the carbon-silicate cycle of weathering (Arnscheidt, C.W. & Rothman, D.H. 2020. Routes to global glaciation. Proceedings of the Royal Society A, v. 476, article 0303 online; DOI: 10.1098/rspa.2020.0303). Their mathematical approach involves two relatively simple, if long-winded, equations based on parameters that express solar heating, albedo, surface temperature and pressure, and the rate of volcanic outgassing of CO2; a simplification that sets biological processes to one side.

Unlike previous models, theirs can simulate varying rates, particularly of changes in solar energy input. The key conclusion of the paper is that if solar heating decreases faster than a threshold rate the more a planet’s surface water is likely to freeze from pole to pole. The authors suggest that a Snowball Earth event would result from a 2% fall in received solar radiation over about ten thousand years: pretty quick in a geological sense. Such a trigger might stem from a volcanic ‘winter’ scenario, an increase in clouds seeded by spores of primitive marine algae or other factors. The real ‘tipping point’ would probably be the high albedo of ice. There is a warning in this for the present, when a variety of means of decreasing solar input have been proposed as a ‘solution’ to global warming.

Because the Earth orbits the Sun in the ‘Goldilocks Zone’ and is volcanically active even global glaciation would be temporary, albeit of the order of millions of years. The cold would have shut down weathering so that volcanic CO2 could slowly build up in the atmosphere: the greenhouse effect would rescue the planet. Further from the Sun, a planet would not have that escape route, regardless of its atmospheric concentration of greenhouse gases: a neat lead-in to another recent paper about the ancient climate of Mars (Grau Galofre, A. et al. 2020. Valley formation on early Mars by subglacial and fluvial erosion. Nature Geoscience, early online article; DOI: 10.1038/s41561-020-0618-x)

A Martian channel system: note later cratering (credit: European Space Agency)

There is a lot of evidence from both high-resolution orbital images of the Martian surface and surface ‘rovers’ that surface water was abundant over a long period in Mars’s early history. The most convincing are networks of channels, mainly in the southern hemisphere highlands. They are not the vast channelled scablands, such as those associated with Valles Marineris, which probably resulted from stupendous outburst floods connected to catastrophic melting of subsurface ice by some means. There are hundreds of channel networks, that resemble counterparts on Earth. Since rainfall and melting of ice and snow have carved most terrestrial channel networks, traditionally those on Mars have been attributed to similar processes during an early warm and wet phase. The warm-early Mars hypothesis extends even to interpreting the smooth low-lying plains of its northern hemisphere – about a third of Mars’s surface area – as the site of an ocean in those ancient times. Of course, a big question is, ‘Where did all that water go?’ Another relates to the fact that the early Sun emitted considerably less radiation 4.5 billion years ago than it does now: a warm-wet early Mars is counterintuitive.

Anna Grau Galofre of the University of British Columbia and co-authors found that many of the networks on Mars clearly differ in morphology from one another, even in small areas of its surface. Drainage networks on Earth conform to far fewer morphological types. By comparing the variability on Mars with channel-network shapes on Earth, the authors found a close match for many with those that formed beneath the ice sheet that covered high latitudes of North America during the last glaciation. Some match drainage patterns typical of surface-water erosion, but both types are present in low Martian latitudes: a suggestion of ‘Snowball Mars’ conditions? The authors reached their conclusions by analysing six mathematical measures that describe channel morphology for over ten thousand individual valley systems. Previous analyses of individual systems discovered on high-resolution images have qualitative comparisons with terrestrial geomorphology

See also: Chu, J. 2020. “Snowball Earths” May Have Been Triggered by a Plunge in Incoming Sunlight – “Be Wary of Speed” (SciTech Daily 29 July 2020); Early Mars was covered in ice sheets, not flowing rivers, researchers say (Science Daily, 3 August 2020)