Late inception of a resiliently oxygenated upper ocean

Episodic body size variations of early Paleozoic trilobites associated with marine redox changes

Body size greatly affects how organisms interact with their environments. However, the macroevolutionary patterns of body size across many major metazoan clades and their constraining mechanisms remain elusive. A new high-resolution body size dataset covering 2435 species from 1091 genera of Cambrian and Ordovician trilobites reveals that body size evolution changes episodically, with three marked reductions in size. Such a pattern rules out a persistent Cope's rule dynamic. Rather, we find a strong temporal link between body size changes and major fluctuations in marine redox, supporting the hypothesis that marine oxygen levels exerted a primary control on the tempo and mode of trilobite body size evolution. These further imply a dominant role for marine oxygen in early animal evolution.

Marine sulphate captures a Paleozoic transition to a modern terrestrial weathering environment

The triple oxygen isotope composition of sulphate minerals has been used to constrain the evolution of Earth's surface environment (e.g., pO2, pCO2 and gross primary productivity) throughout the Proterozoic Eon. This approach presumes the incorporation of atmospheric O2 atoms into riverine sulphate via the oxidative weathering of pyrite. However, this is not borne out in recent geological or modern sulphate records, where an atmospheric signal is imperceptible and where terrestrial pyrite weathering occurs predominantly in bedrock fractures that are physically more removed from atmospheric O2. To better define the transition from a Proterozoic to a modern-like weathering regime, here we present new measurements from twelve marine evaporite basins spanning the Phanerozoic. These data display a step-like transition in the triple oxygen isotope composition of evaporite sulphate during the mid-Paleozoic (420 to 387.7 million years ago). We propose that the evolution of early root systems deepened the locus of pyrite oxidation and reduced the incorporation of O2 into sulphate. Further, the early Devonian proliferation of land plants increased terrestrial organic carbon burial, releasing free oxygen that fueled increased redox recycling of soil-bound iron and resulted in the final rise in pO2 to modern-like levels.

Evolution of the iodine cycle and the late stabilization of the Earth's ozone layer

The origin of complex life and the evolution of terrestrial ecosystems are fundamental aspects of the natural history on Earth. Here, we present evidence for a protracted stabilization of the Earth's ozone layer. The destruction of atmospheric ozone today is inherently linked to the cycling of marine and atmospheric iodine. Supported by multiple independent lines of geological evidence and examined through an iodine mass balance model, we find that elevated marine iodide content prevailed through most of Earth's history. Since the rise of oxygen ~2.4 billion years ago, high marine iodide concentrations would have led to significant inorganic iodine emissions to the atmosphere, facilitating catalytic ozone destruction and resulting in atmospheric ozone instability with periodic or persistently lower ozone levels. At a global scale, unstable and low ozone levels likely persisted for about two billion years until the early Phanerozoic, roughly 0.5 billion years ago. The delayed stabilization of the Earth's ozone layer holds significant implications for the tempo and direction of the evolution of life, in particular life on land.

Opportunities for Earth Observation to Inform Risk Management for Ocean Tipping Points

As climate change continues, the likelihood of passing critical thresholds or tipping points increases. Hence, there is a need to advance the science for detecting such thresholds. In this paper, we assess the needs and opportunities for Earth Observation (EO, here understood to refer to satellite observations) to inform society in responding to the risks associated with ten potential large-scale ocean tipping elements: Atlantic Meridional Overturning Circulation; Atlantic Subpolar Gyre; Beaufort Gyre; Arctic halocline; Kuroshio Large Meander; deoxygenation; phytoplankton; zooplankton; higher level ecosystems (including fisheries); and marine biodiversity. We review current scientific understanding and identify specific EO and related modelling needs for each of these tipping elements. We draw out some generic points that apply across several of the elements. These common points include the importance of maintaining long-term, consistent time series; the need to combine EO data consistently with in situ data types (including subsurface), for example through data assimilation; and the need to reduce or work with current mismatches in resolution (in both directions) between climate models and EO datasets. Our analysis shows that developing EO, modelling and prediction systems together, with understanding of the strengths and limitations of each, provides many promising paths towards monitoring and early warning systems for tipping, and towards the development of the next generation of climate models.

Phanerozoic co-evolution of O2-CO2 and ocean habitability

This perspective reviews how atmospheric compositions, animals and marine algae evolved together to determine global ocean habitability during the past 500 million years.

Deep-water first occurrences of Ediacara biota prior to the Shuram carbon isotope excursion in the Wernecke Mountains, Yukon, Canada

Ediacara-type macrofossils appear as early as ~575 Ma in deep-water facies of the Drook Formation of the Avalon Peninsula, Newfoundland, and the Nadaleen Formation of Yukon and Northwest Territories, Canada. Our ability to assess whether a deep-water origination of the Ediacara biota is a genuine reflection of evolutionary succession, an artifact of an incomplete stratigraphic record, or a bathymetrically controlled biotope is limited by a lack of geochronological constraints and detailed shelf-to-slope transects of Ediacaran continental margins. The Ediacaran Rackla Group of the Wernecke Mountains, NW Canada, represents an ideal shelf-to-slope depositional system to understand the spatiotemporal and environmental context of Ediacara-type organisms' stratigraphic occurrence. New sedimentological and paleontological data presented herein from the Wernecke Mountains establish a stratigraphic framework relating shelfal strata in the Goz/Corn Creek area to lower slope deposits in the Nadaleen River area. We report new discoveries of numerous Aspidella hold-fast discs, indicative of frondose Ediacara organisms, from deep-water slope deposits of the Nadaleen Formation stratigraphically below the Shuram carbon isotope excursion (CIE) in the Nadaleen River area. Such fossils are notably absent in coeval shallow-water strata in the Goz/Corn Creek region despite appropriate facies for potential preservation. The presence of pre-Shuram CIE Ediacara-type fossils occurring only in deep-water facies within a basin that has equivalent well-preserved shallow-water facies provides the first stratigraphic paleobiological support for a deep-water origination of the Ediacara biota. In contrast, new occurrences of Ediacara-type fossils (including juvenile fronds, Beltanelliformis, Aspidella, annulated tubes, and multiple ichnotaxa) are found above the Shuram CIE in both deep- and shallow-water deposits of the Blueflower Formation. Given existing age constraints on the Shuram CIE, it appears that Ediacaran organisms may have originated in the deeper ocean and lived there for up to ~15 million years before migrating into shelfal environments in the terminal Ediacaran. This indicates unique ecophysiological constraints likely shaped the initial habitat preference and later environmental expansion of the Ediacara biota.

Late Paleozoic oxygenation of marine environments supported by dolomite U-Pb dating

Understanding causal relationships between evolution and ocean oxygenation hinges on reliable reconstructions of marine oxygen levels, typically from redox-sensitive geochemical proxies. Here, we develop a proxy, using dolomite U-Pb geochronology, to reconstruct seawater U/Pb ratios. Dolomite samples consistently give U-Pb dates and initial 207Pb/206Pb ratios lower than expected from their stratigraphic ages. These observations are explained by resetting of the U-Pb system long after deposition; the magnitude of deviations from expected initial 207Pb/206Pb are a function of the redox-sensitive U/Pb ratios during deposition. Reconstructed initial U/Pb ratios increased notably in the late-Paleozoic, reflecting an increase in oxygenation of marine environments at that time. This timeline is consistent with documented shifts in some other redox proxies and supports evolution-driven mechanisms for the oxygenation of late-Paleozoic marine environments, as well as suggestions that early animals thrived in oceans that on long time scales were oxygen-limited compared to today.

The role of iron in the formation of Ediacaran 'death masks'

The Ediacara biota are an enigmatic group of Neoproterozoic soft-bodied fossils that mark the first major radiation of complex eukaryotic and macroscopic life. These fossils are thought to have been preserved via pyritic "death masks" mediated by seafloor microbial mats, though little about the chemical constraints of this preservational pathway is known, in particular surrounding the role of bioavailable iron in death mask formation and preservational fidelity. In this study, we perform decay experiments on both diploblastic and triploblastic animals under a range of simulated sedimentary iron concentrations, in order to characterize the role of iron in the preservation of Ediacaran organisms. After 28 days of decay, we demonstrate the first convincing "death masks" produced under experimental laboratory conditions composed of iron sulfide and probable oxide veneers. Moreover, our results demonstrate that the abundance of iron in experiments is not the sole control on death mask formation, but also tissue histology and the availability of nucleation sites. This illustrates that Ediacaran preservation via microbial death masks need not be a "perfect storm" of paleoenvironmental porewater and sediment chemistry, but instead can occur under a range of conditions.

Evolution of the crustal phosphorus reservoir

The release of phosphorus (P) from crustal rocks during weathering plays a key role in determining the size of Earth's biosphere, yet the concentration of P in crustal rocks over time remains controversial. Here, we combine spatial, temporal, and chemical measurements of preserved rocks to reconstruct the lithological and chemical evolution of Earth's continental crust. We identify a threefold increase in average crustal P concentrations across the Neoproterozoic-Phanerozoic boundary (600 to 400 million years), showing that preferential biomass burial on shelves acted to progressively concentrate P within continental crust. Rapid compositional change was made possible by massive removal of ancient P-poor rock and deposition of young P-rich sediment during an episode of enhanced global erosion. Subsequent weathering of newly P-rich crust led to increased riverine P fluxes to the ocean. Our results suggest that global erosion coupled to sedimentary P-enrichment forged a markedly nutrient-rich crust at the dawn of the Phanerozoic.

Reduction in animal abundance and oxygen availability during and after the end-Triassic mass extinction

The end-Triassic biodiversity crisis was one of the most severe mass extinctions in the history of animal life. However, the extent to which the loss of taxonomic diversity was coupled with a reduction in organismal abundance remains to be quantified. Further, the temporal relationship between organismal abundance and local marine redox conditions is lacking in carbonate sections. To address these questions, we measured skeletal grain abundance in shallow-marine limestones by point counting 293 thin sections from four stratigraphic sections across the Triassic/Jurassic boundary in the Lombardy Basin and Apennine Platform of western Tethys. Skeletal abundance decreased abruptly across the Triassic/Jurassic boundary in all stratigraphic sections. The abundance of skeletal organisms remained low throughout the lower-middle Hettangian strata and began to rebound during the late Hettangian and early Sinemurian. A two-way ANOVA indicates that sample age (p < .01, η²  = 0.30) explains more of the variation in skeletal abundance than the depositional environment or paleobathymetry (p < .01, η²  = 0.15). Measured I/Ca ratios, a proxy for local shallow-marine redox conditions, show this same pattern with the lowest I/Ca ratios occurring in the early Hettangian. The close correspondence between oceanic water column oxygen levels and skeletal abundance indicates a connection between redox conditions and benthic organismal abundance across the Triassic/Jurassic boundary. These findings indicate that the end-Triassic mass extinction reduced not only the biodiversity but also the carrying capacity for skeletal organisms in early Hettangian ecosystems, adding to evidence that mass extinction of species generally leads to mass rarity among survivors.

Extreme variability in atmospheric oxygen levels in the late Precambrian

Mapping the history of atmospheric O2 during the late Precambrian is vital for evaluating potential links to animal evolution. Ancient O2 levels are often inferred from geochemical analyses of marine sediments, leading to the assumption that the Earth experienced a stepwise increase in atmospheric O2 during the Neoproterozoic. However, the nature of this hypothesized oxygenation event remains unknown, with suggestions of a more dynamic O2 history in the oceans and major uncertainty over any direct connection between the marine realm and atmospheric O2. Here, we present a continuous quantitative reconstruction of atmospheric O2 over the past 1.5 billion years using an isotope mass balance approach that combines bulk geochemistry and tectonic recycling rate calculations. We predict that atmospheric O2 levels during the Neoproterozoic oscillated between ~1 and ~50% of the present atmospheric level. We conclude that there was no simple unidirectional rise in atmospheric O2 during the Neoproterozoic, and the first animals evolved against a backdrop of extreme O2 variability.

Continental configuration controls ocean oxygenation during the Phanerozoic

The early evolutionary and much of the extinction history of marine animals is thought to be driven by changes in dissolved oxygen concentrations ([O₂]) in the ocean^1-3. In turn, [O₂] is widely assumed to be dominated by the geological history of atmospheric oxygen (pO₂)^4,5. Here, by contrast, we show by means of a series of Earth system model experiments how continental rearrangement during the Phanerozoic Eon drives profound variations in ocean oxygenation and induces a fundamental decoupling in time between upper-ocean and benthic [O₂]. We further identify the presence of state transitions in the global ocean circulation, which lead to extensive deep-ocean anoxia developing in the early Phanerozoic even under modern pO₂. Our finding that ocean oxygenation oscillates over stable thousand-year (kyr) periods also provides a causal mechanism that might explain elevated rates of metazoan radiation and extinction during the early Palaeozoic Era⁶. The absence, in our modelling, of any simple correlation between global climate and ocean ventilation, and the occurrence of profound variations in ocean oxygenation independent of atmospheric pO₂, presents a challenge to the interpretation of marine redox proxies, but also points to a hitherto unrecognized role for continental configuration in the evolution of the biosphere.

Decreasing Phanerozoic extinction intensity as a consequence of Earth surface oxygenation and metazoan ecophysiology

The decline in background extinction rates of marine animals through geologic time is an established but unexplained feature of the Phanerozoic fossil record. There is also growing consensus that the ocean and atmosphere did not become oxygenated to near-modern levels until the mid-Paleozoic, coinciding with the onset of generally lower extinction rates. Physiological theory provides us with a possible causal link between these two observations-predicting that the synergistic impacts of oxygen and temperature on aerobic respiration would have made marine animals more vulnerable to ocean warming events during periods of limited surface oxygenation. Here, we evaluate the hypothesis that changes in surface oxygenation exerted a first-order control on extinction rates through the Phanerozoic using a combined Earth system and ecophysiological modeling approach. We find that although continental configuration, the efficiency of the biological carbon pump in the ocean, and initial climate state all impact the magnitude of modeled biodiversity loss across simulated warming events, atmospheric oxygen is the dominant predictor of extinction vulnerability, with metabolic habitat viability and global ecophysiotype extinction exhibiting inflection points around 40% of present atmospheric oxygen. Given this is the broad upper limit for estimates of early Paleozoic oxygen levels, our results are consistent with the relative frequency of high-magnitude extinction events (particularly those not included in the canonical big five mass extinctions) early in the Phanerozoic being a direct consequence of limited early Paleozoic oxygenation and temperature-dependent hypoxia responses.

A long-term record of early to mid-Paleozoic marine redox change

The extent to which Paleozoic oceans differed from Neoproterozoic oceans and the causal relationship between biological evolution and changing environmental conditions are heavily debated. Here, we report a nearly continuous record of seafloor redox change from the deep-water upper Cambrian to Middle Devonian Road River Group of Yukon, Canada. Bottom waters were largely anoxic in the Richardson trough during the entirety of Road River Group deposition, while independent evidence from iron speciation and Mo/U ratios show that the biogeochemical nature of anoxia changed through time. Both in Yukon and globally, Ordovician through Early Devonian anoxic waters were broadly ferruginous (nonsulfidic), with a transition toward more euxinic (sulfidic) conditions in the mid-Early Devonian (Pragian), coincident with the early diversification of vascular plants and disappearance of graptolites. This ~80-million-year interval of the Paleozoic characterized by widespread ferruginous bottom waters represents a persistence of Neoproterozoic-like marine redox conditions well into the Phanerozoic.

Marine iodine emissions in a changing world

Iodine is a critical trace element involved in many diverse and important processes in the Earth system. The importance of iodine for human health has been known for over a century, with low iodine in the diet being linked to goitre, cretinism and neonatal death. Research over the last few decades has shown that iodine has significant impacts on tropospheric photochemistry, ultimately impacting climate by reducing the radiative forcing of ozone (O₃) and air quality by reducing extreme O₃ concentrations in polluted regions. Iodine is naturally present in the ocean, predominantly as aqueous iodide and iodate. The rapid reaction of sea-surface iodide with O₃ is believed to be the largest single source of gaseous iodine to the atmosphere. Due to increased anthropogenic O₃, this release of iodine is believed to have increased dramatically over the twentieth century, by as much as a factor of 3. Uncertainties in the marine iodine distribution and global cycle are, however, major constraints in the effective prediction of how the emissions of iodine and its biogeochemical cycle may change in the future or have changed in the past. Here, we present a synthesis of recent results by our team and others which bring a fresh perspective to understanding the global iodine biogeochemical cycle. In particular, we suggest that future climate-induced oceanographic changes could result in a significant change in aqueous iodide concentrations in the surface ocean, with implications for atmospheric air quality and climate.

A persistently low level of atmospheric oxygen in Earth's middle age

Resolving how Earth surface redox conditions evolved through the Proterozoic Eon is fundamental to understanding how biogeochemical cycles have changed through time. The redox sensitivity of cerium relative to other rare earth elements and its uptake in carbonate minerals make the Ce anomaly (Ce/Ce*) a particularly useful proxy for capturing redox conditions in the local marine environment. Here, we report Ce/Ce* data in marine carbonate rocks through 3.5 billion years of Earth's history, focusing in particular on the mid-Proterozoic Eon (i.e., 1.8 - 0.8 Ga). To better understand the role of atmospheric oxygenation, we use Ce/Ce* data to estimate the partial pressure of atmospheric oxygen (pO₂) through this time. Our thermodynamics-based modeling supports a major rise in atmospheric oxygen level in the aftermath of the Great Oxidation Event (~ 2.4 Ga), followed by invariant pO₂ of about 1% of present atmospheric level through most of the Proterozoic Eon (2.4 to 0.65 Ga).

The evolution of complex life and the stabilization of the Earth system

The half-billion-year history of animal evolution is characterized by decreasing rates of background extinction. Earth's increasing habitability for animals could result from several processes: (i) a decrease in the intensity of interactions among species that lead to extinctions; (ii) a decrease in the prevalence or intensity of geological triggers such as flood basalt eruptions and bolide impacts; (iii) a decrease in the sensitivity of animals to environmental disturbance; or (iv) an increase in the strength of stabilizing feedbacks within the climate system and biogeochemical cycles. There is no evidence that the prevalence or intensity of interactions among species or geological extinction triggers have decreased over time. There is, however, evidence from palaeontology, geochemistry and comparative physiology that animals have become more resilient to an environmental change and that the evolution of complex life has, on the whole, strengthened stabilizing feedbacks in the climate system. The differential success of certain phyla and classes appears to result, at least in part, from the anatomical solutions to the evolution of macroscopic size that were arrived at largely during Ediacaran and Cambrian time. Larger-bodied animals, enabled by increased anatomical complexity, were increasingly able to mix the marine sediment and water columns, thus promoting stability in biogeochemical cycles. In addition, body plans that also facilitated ecological differentiation have tended to be associated with lower rates of extinction. In this sense, Cambrian solutions to Cambrian problems have had a lasting impact on the trajectory of complex life and, in turn, fundamental properties of the Earth system.

Reconciling proxy records and models of Earth's oxygenation during the Neoproterozoic and Palaeozoic

A hypothesized rise in oxygen levels in the Neoproterozoic, dubbed the Neoproterozoic Oxygenation Event, has been repeatedly linked to the origin and rise of animal life. However, a new body of work has emerged over the past decade that questions this narrative. We explore available proxy records of atmospheric and marine oxygenation and, considering the unique systematics of each geochemical system, attempt to reconcile the data. We also present new results from a comprehensive COPSE biogeochemical model that combines several recent additions, to create a continuous model record from 850 to 250 Ma. We conclude that oxygen levels were intermediate across the Ediacaran and early Palaeozoic, and highly dynamic. Stable, modern-like conditions were not reached until the Late Palaeozoic. We therefore propose that the terms Neoproterozoic Oxygenation Window and Palaeozoic Oxygenation Event are more appropriate descriptors of the rise of oxygen in Earth's atmosphere and oceans.

Ediacaran reorganization of the marine phosphorus cycle

The evolution of macroscopic animals in the latest Proterozoic Eon is associated with many changes in the geochemical environment, but the sequence of cause and effect remains a topic of intense research and debate. In this study, we use two apparently paradoxical observations—that massively phosphorus-rich rocks first appear at this time, and that the median P content of rocks does not change—to argue for a change in internal marine P cycling associated with rising sulfate levels. We argue that this change was self-sustaining, setting in motion a cascade of biogeochemical transformations that led to conditions favorable for major ecological and evolutionary change.

The Ediacaran Period (635 to 541 Ma) marks the global transition to a more productive biosphere, evidenced by increased availability of food and oxidants, the appearance of macroscopic animals, significant populations of eukaryotic phytoplankton, and the onset of massive phosphorite deposition. We propose this entire suite of changes results from an increase in the size of the deep-water marine phosphorus reservoir, associated with rising sulfate concentrations and increased remineralization of organic P by sulfate-reducing bacteria. Simple mass balance calculations, constrained by modern anoxic basins, suggest that deep-water phosphate concentrations may have increased by an order of magnitude without any increase in the rate of P input from the continents. Strikingly, despite a major shift in phosphorite deposition, a new compilation of the phosphorus content of Neoproterozoic and early Paleozoic shows little secular change in median values, supporting the view that changes in remineralization and not erosional P fluxes were the principal drivers of observed shifts in phosphorite accumulation. The trigger for these changes may have been transient Neoproterozoic weathering events whose biogeochemical consequences were sustained by a set of positive feedbacks, mediated by the oxygen and sulfur cycles, that led to permanent state change in biogeochemical cycling, primary production, and biological diversity by the end of the Ediacaran Period.

Persistent global marine euxinia in the early Silurian

The second pulse of the Late Ordovician mass extinction occurred around the Hirnantian-Rhuddanian boundary (~444 Ma) and has been correlated with expanded marine anoxia lasting into the earliest Silurian. Characterization of the Hirnantian ocean anoxic event has focused on the onset of anoxia, with global reconstructions based on carbonate δ²³⁸U modeling. However, there have been limited attempts to quantify uncertainty in metal isotope mass balance approaches. Here, we probabilistically evaluate coupled metal isotopes and sedimentary archives to increase constraint. We present iron speciation, metal concentration, δ⁹⁸Mo and δ²³⁸U measurements of Rhuddanian black shales from the Murzuq Basin, Libya. We evaluate these data (and published carbonate δ²³⁸U data) with a coupled stochastic mass balance model. Combined statistical analysis of metal isotopes and sedimentary sinks provides uncertainty-bounded constraints on the intensity of Hirnantian-Rhuddanian euxinia. This work extends the duration of anoxia to >3 Myrs – notably longer than well-studied Mesozoic ocean anoxic events.

The Late Ordovician mass extinction has been attributed to extended marine anoxia. Here, the authors use a metal isotope mass balance model and find the marine anoxic event lasted over 3 million years, notably longer than the anoxic event associated with the Permian-Triassic extinction and Cretaceous ocean anoxic events.

Global sea-surface iodide observations, 1967-2018

The marine iodine cycle has significant impacts on air quality and atmospheric chemistry. Specifically, the reaction of iodide with ozone in the top few micrometres of the surface ocean is an important sink for tropospheric ozone (a pollutant gas) and the dominant source of reactive iodine to the atmosphere. Sea surface iodide parameterisations are now being implemented in air quality models, but these are currently a major source of uncertainty. Relatively little observational data is available to estimate the global surface iodide concentrations, and this data has not hitherto been openly available in a collated, digital form. Here we present all available sea surface (<20 m depth) iodide observations. The dataset includes values digitised from published manuscripts, published and unpublished data supplied directly by the originators, and data obtained from repositories. It contains 1342 data points, and spans latitudes from 70°S to 68°N, representing all major basins. The data may be used to model sea surface iodide concentrations or as a reference for future observations.

Susceptibility of algae to Cr toxicity reveals contrasting metal management strategies

At the Paleozoic-Mesozoic boundary, the dominance of marine eukaryotic algae shifted from the green (chlorophyll b) to the red (chlorophyll c) superfamily. Selection pressures caused by the bioavailability of trace metals associated with increasing oxygenation of the ocean may have played a key role in this algal revolution. From a scan of elemental compositions, a significant difference in the cellular Cr/P quota was found between the two superfamilies. Here, the different responses to high levels of Cr exposure reveal contrasting strategies for metal uptake and homeostasis in these algal lineages. At high Cr(VI) concentrations, red lineages experience growth inhibition through reduced photosynthetic capability, while green lineages are completely unaffected. Moreover, Cr(VI) has a more significant impact on the metallomes of red lineage algae, in which metal/P ratios increased with increasing Cr(VI) concentration for many trace elements. Green algae have higher specificity transporters to prevent Cr(VI) from entering the cell, and more specific intracellular stores of Cr within the membrane fraction than the red algae, which accumulate more Cr mistakenly in the cytosol fraction via lower affinity transport mechanisms. Green algal approaches require greater nutrient investments in the more numerous transport proteins required and management of specific metals, a strategy better adapted to the resource-rich coastal waters. By contrast, the red algae are nutrient-efficient with fewer and less discriminate metal transporters, which can be fast and better adapted in the oligotrophic, oxygenated open ocean, which has prevailed since the deepening of the oxygen minimum zones at the start of the Mesozoic era.

Energy metabolism in anaerobic eukaryotes and Earth's late oxygenation

Eukaryotes arose about 1.6 billion years ago, at a time when oxygen levels were still very low on Earth, both in the atmosphere and in the ocean. According to newer geochemical data, oxygen rose to approximately its present atmospheric levels very late in evolution, perhaps as late as the origin of land plants (only about 450 million years ago). It is therefore natural that many lineages of eukaryotes harbor, and use, enzymes for oxygen-independent energy metabolism. This paper provides a concise overview of anaerobic energy metabolism in eukaryotes with a focus on anaerobic energy metabolism in mitochondria. We also address the widespread assumption that oxygen improves the overall energetic state of a cell. While it is true that ATP yield from glucose or amino acids is increased in the presence of oxygen, it is also true that the synthesis of biomass costs thirteen times more energy per cell in the presence of oxygen than in anoxic conditions. This is because in the reaction of cellular biomass with O2, the equilibrium lies very far on the side of CO2. The absence of oxygen offers energetic benefits of the same magnitude as the presence of oxygen. Anaerobic and low oxygen environments are ancient. During evolution, some eukaryotes have specialized to life in permanently oxic environments (life on land), other eukaryotes have remained specialized to low oxygen habitats. We suggest that the Km of mitochondrial cytochrome c oxidase of 0.1-10 μM for O2, which corresponds to about 0.04%-4% (avg. 0.4%) of present atmospheric O2 levels, reflects environmental O2 concentrations that existed at the time that the eukaryotes arose.

Upper ocean oxygenation, evolution of RuBisCO and the Phanerozoic succession of phytoplankton

Evidence is compiled to demonstrate a redox scale within Earth's photosynthesisers that correlates the specificity of their RuBisCO with organismal metabolic tolerance to anoxia, and ecological selection by dissolved O2/CO2 and nutrients. The Form 1B RuBisCO found in the chlorophyte green algae, has a poor selectivity between the two dissolved substrates, O2 and CO2, at the active site. This enzyme appears adapted to lower O2/CO2 ratios, or more "anoxic" conditions and therefore requires additional energetic or nutrient investment in a carbon concentrating mechanism (CCM) to boost the intracellular CO2/O2 ratio and maintain competitive carboxylation rates under increasingly high O2/CO2 conditions in the environment. By contrast the coccolithophores and diatoms evolved containing the more selective Rhodophyte Form 1D RuBisCO, better adapted to a higher O2/CO2 ratio, or more oxic conditions. This Form 1D RuBisCO requires lesser energetic or nutrient investment in a CCM to attain high carboxylation rates under environmentally high O2/CO2 ratios. Such a physiological relationship may underpin the succession of phytoplankton in the Phanerozoic oceans: the coccolithophores and diatoms took over the oceanic realm from the incumbent cyanobacteria and green algae when the upper ocean became persistently oxygenated, alkaline and more oligotrophic. The facultatively anaerobic green algae, able to tolerate the anoxic conditions of the water column and a periodically inundated soil, were better poised to adapt to the fluctuating anoxia associated with periods of submergence and emergence and transition onto the land. The induction of a CCM may exert a natural limit to the improvement of RuBisCO efficiency over Earth history. Rubisco specificity appears to adapt on the timescale of ∼100 Myrs. So persistent elevation of CO2/O2 ratios in the intracellular environment around the enzyme, may induce a relaxation in RuBisCO selectivity for CO2 relative to O2. The most efficient RuBisCO for net carboxylation is likely to be found in CCM-lacking algae that have been exposed to hyperoxic conditions for at least 100 Myrs, such as intertidal brown seaweeds.

Stepwise oxygenation of the Paleozoic atmosphere

Oxygen is essential for animal life, and while geochemical proxies have been instrumental in determining the broad evolutionary history of oxygen on Earth, much of our insight into Phanerozoic oxygen comes from biogeochemical modelling. The GEOCARBSULF model utilizes carbon and sulphur isotope records to produce the most detailed history of Phanerozoic atmospheric O₂ currently available. However, its predictions for the Paleozoic disagree with geochemical proxies, and with non-isotope modelling. Here we show that GEOCARBSULF oversimplifies the geochemistry of sulphur isotope fractionation, returning unrealistic values for the O₂ sourced from pyrite burial when oxygen is low. We rebuild the model from first principles, utilizing an improved numerical scheme, the latest carbon isotope data, and we replace the sulphur cycle equations in line with forwards modelling approaches. Our new model, GEOCARBSULFOR, produces a revised, highly-detailed prediction for Phanerozoic O₂ that is consistent with available proxy data, and independently supports a Paleozoic Oxygenation Event, which likely contributed to the observed radiation of complex, diverse fauna at this time.

The GEOCARBSULF model provides the most detailed reconstructions of Phanerozoic O₂, but its predictions are not supported by geochemical data. Here, a GEOCARBSULF model rebuilt from first principles, with the addition of an amended sulphur cycle and the latest isotope records, supports a Paleozoic Oxygenation Event.