Onset of coupled atmosphere-ocean oxygenation 2.3 billion years ago.

The initial rise of molecular oxygen (O2) shortly after the Archaean-Proterozoic transition 2.5 billion years ago was more complex than the single step-change once envisioned. Sulfur mass-independent fractionation records suggest that the rise of atmospheric O2 was oscillatory, with multiple returns to an anoxic state until perhaps 2.2 billion years ago1-3. Yet few constraints exist for contemporaneous marine oxygenation dynamics, precluding a holistic understanding of planetary oxygenation. Here we report thallium (Tl) isotope ratio and redox-sensitive element data for marine shales from the Transvaal Supergroup, South Africa. Synchronous with sulfur isotope evidence of atmospheric oxygenation in the same shales3, we found lower authigenic 205Tl/203Tl ratios indicative of widespread manganese oxide burial on an oxygenated seafloor and higher redox-sensitive element abundances consistent with expanded oxygenated waters. Both signatures disappear when the sulfur isotope data indicate a brief return to an anoxic atmospheric state. Our data connect recently identified atmospheric O2 dynamics on early Earth with the marine realm, marking an important turning point in Earth's redox history away from heterogeneous and highly localized 'oasis'-style oxygenation.

Uncovering the Ediacaran phosphorus cycle.

Phosphorus is a limiting nutrient that is thought to control oceanic oxygen levels to a large extent1-3. A possible increase in marine phosphorus concentrations during the Ediacaran Period (about 635-539 million years ago) has been proposed as a driver for increasing oxygen levels4-6. However, little is known about the nature and evolution of phosphorus cycling during this time4. Here we use carbonate-associated phosphate (CAP) from six globally distributed sections to reconstruct oceanic phosphorus concentrations during a large negative carbon-isotope excursion-the Shuram excursion (SE)-which co-occurred with global oceanic oxygenation7-9. Our data suggest pulsed increases in oceanic phosphorus concentrations during the falling and rising limbs of the SE. Using a quantitative biogeochemical model, we propose that this observation could be explained by carbon dioxide and phosphorus release from marine organic-matter oxidation primarily by sulfate, with further phosphorus release from carbon-dioxide-driven weathering on land. Collectively, this may have resulted in elevated organic-pyrite burial and ocean oxygenation. Our CAP data also seem to suggest equivalent oceanic phosphorus concentrations under maximum and minimum extents of ocean anoxia across the SE. This observation may reflect decoupled phosphorus and ocean anoxia cycles, as opposed to their coupled nature in the modern ocean. Our findings point to external stimuli such as sulfate weathering rather than internal oceanic phosphorus-oxygen cycling alone as a possible control on oceanic oxygenation in the Ediacaran. In turn, this may help explain the prolonged rise of atmospheric oxygen levels.

A 200-million-year delay in permanent atmospheric oxygenation.

The rise of atmospheric oxygen fundamentally changed the chemistry of surficial environments and the nature of Earth's habitability1. Early atmospheric oxygenation occurred over a protracted period of extreme climatic instability marked by multiple global glaciations2,3, with the initial rise of oxygen concentration to above 10-5 of the present atmospheric level constrained to about 2.43 billion years ago4,5. Subsequent fluctuations in atmospheric oxygen levels have, however, been reported to have occurred until about 2.32 billion years ago4, which represents the estimated timing of irreversible oxygenation of the atmosphere6,7. Here we report a high-resolution reconstruction of atmospheric and local oceanic redox conditions across the final two glaciations of the early Palaeoproterozoic era, as documented by marine sediments from the Transvaal Supergroup, South Africa. Using multiple sulfur isotope and iron-sulfur-carbon systematics, we demonstrate continued oscillations in atmospheric oxygen levels after about 2.32 billion years ago that are linked to major perturbations in ocean redox chemistry and climate. Oxygen levels thus fluctuated across the threshold of 10-5 of the present atmospheric level for about 200 million years, with permanent atmospheric oxygenation finally arriving with the Lomagundi carbon isotope excursion at about 2.22 billion years ago, some 100 million years later than currently estimated.

Onset of the aerobic nitrogen cycle during the Great Oxidation Event.

The rise of oxygen on the early Earth (about 2.4 billion years ago) caused a reorganization of marine nutrient cycles, including that of nitrogen, which is important for controlling global primary productivity. However, current geochemical records lack the temporal resolution to address the nature and timing of the biogeochemical response to oxygenation directly. Here we couple records of ocean redox chemistry with nitrogen isotope (15N/14N) values from approximately 2.31-billion-year-old shales of the Rooihoogte and Timeball Hill formations in South Africa, deposited during the early stages of the first rise in atmospheric oxygen on the Earth (the Great Oxidation Event). Our data fill a gap of about 400 million years in the temporal 15N/14N record and provide evidence for the emergence of a pervasive aerobic marine nitrogen cycle. The interpretation of our nitrogen isotope data in the context of iron speciation and carbon isotope data suggests biogeochemical cycling across a dynamic redox boundary, with primary productivity fuelled by chemoautotrophic production and a nitrogen cycle dominated by nitrogen loss processes using newly available marine oxidants. This chemostratigraphic trend constrains the onset of widespread nitrate availability associated with ocean oxygenation. The rise of marine nitrate could have allowed for the rapid diversification and proliferation of nitrate-using cyanobacteria and, potentially, eukaryotic phytoplankton.

Limited oxygen production in the Mesoarchean ocean.

The Archean Eon was a time of predominantly anoxic Earth surface conditions, where anaerobic processes controlled bioessential element cycles. In contrast to "oxygen oases" well documented for the Neoarchean [2.8 to 2.5 billion years ago (Ga)], the magnitude, spatial extent, and underlying causes of possible Mesoarchean (3.2 to 2.8 Ga) surface-ocean oxygenation remain controversial. Here, we report δ15N and δ13C values coupled with local seawater redox data for Mesoarchean shales of the Mozaan Group (Pongola Supergroup, South Africa) that were deposited during an episode of enhanced Mn (oxyhydr)oxide precipitation between ∼2.95 and 2.85 Ga. Iron and Mn redox systematics are consistent with an oxygen oasis in the Mesoarchean anoxic ocean, but δ15N data indicate a Mo-based diazotrophic biosphere with no compelling evidence for a significant aerobic nitrogen cycle. We propose that in contrast to the Neoarchean, dissolved O2 levels were either too low or too limited in extent to develop a large and stable nitrate reservoir in the Mesoarchean ocean. Since biological N2 fixation was evidently active in this environment, the growth and proliferation of O2-producing organisms were likely suppressed by nutrients other than nitrogen (e.g., phosphorus), which would have limited the expansion of oxygenated conditions during the Mesoarchean.

Biological regulation of atmospheric chemistry en route to planetary oxygenation.

Emerging evidence suggests that atmospheric oxygen may have varied before rising irreversibly ∼2.4 billion years ago, during the Great Oxidation Event (GOE). Significantly, however, pre-GOE atmospheric aberrations toward more reducing conditions-featuring a methane-derived organic-haze-have recently been suggested, yet their occurrence, causes, and significance remain underexplored. To examine the role of haze formation in Earth's history, we targeted an episode of inferred haze development. Our redox-controlled (Fe-speciation) carbon- and sulfur-isotope record reveals sustained systematic stratigraphic covariance, precluding nonatmospheric explanations. Photochemical models corroborate this inference, showing Δ36S/Δ33S ratios are sensitive to the presence of haze. Exploiting existing age constraints, we estimate that organic haze developed rapidly, stabilizing within ∼0.3 ± 0.1 million years (Myr), and persisted for upward of ∼1.4 ± 0.4 Myr. Given these temporal constraints, and the elevated atmospheric CO2 concentrations in the Archean, the sustained methane fluxes necessary for haze formation can only be reconciled with a biological source. Correlative δ13COrg and total organic carbon measurements support the interpretation that atmospheric haze was a transient response of the biosphere to increased nutrient availability, with methane fluxes controlled by the relative availability of organic carbon and sulfate. Elevated atmospheric methane concentrations during haze episodes would have expedited planetary hydrogen loss, with a single episode of haze development providing up to 2.6-18 × 1018 moles of O2 equivalents to the Earth system. Our findings suggest the Neoarchean likely represented a unique state of the Earth system where haze development played a pivotal role in planetary oxidation, hastening the contingent biological innovations that followed.

Pathways for Neoarchean pyrite formation constrained by mass-independent sulfur isotopes.

It is generally thought that the sulfate reduction metabolism is ancient and would have been established well before the Neoarchean. It is puzzling, therefore, that the sulfur isotope record of the Neoarchean is characterized by a signal of atmospheric mass-independent chemistry rather than a strong overprint by sulfate reducers. Here, we present a study of the four sulfur isotopes obtained using secondary ion MS that seeks to reconcile a number of features seen in the Neoarchean sulfur isotope record. We suggest that Neoarchean ocean basins had two coexisting, significantly sized sulfur pools and that the pathways forming pyrite precursors played an important role in establishing how the isotopic characteristics of each of these pools was transferred to the sedimentary rock record. One of these pools is suggested to be a soluble (sulfate) pool, and the other pool (atmospherically derived elemental sulfur) is suggested to be largely insoluble and unreactive until it reacts with hydrogen sulfide. We suggest that the relative contributions of these pools to the formation of pyrite depend on both the accumulation of the insoluble pool and the rate of sulfide production in the pyrite-forming environments. We also suggest that the existence of a significant nonsulfate pool of reactive sulfur has masked isotopic evidence for the widespread activity of sulfate reducers in the rock record.

Fluctuations in Precambrian atmospheric oxygenation recorded by chromium isotopes.

Geochemical data suggest that oxygenation of the Earth's atmosphere occurred in two broad steps. The first rise in atmospheric oxygen is thought to have occurred between approximately 2.45 and 2.2 Gyr ago, leading to a significant increase in atmospheric oxygen concentrations and concomitant oxygenation of the shallow surface ocean. The second increase in atmospheric oxygen appears to have taken place in distinct stages during the late Neoproterozoic era ( approximately 800-542 Myr ago), ultimately leading to oxygenation of the deep ocean approximately 580 Myr ago, but details of the evolution of atmospheric oxygenation remain uncertain. Here we use chromium (Cr) stable isotopes from banded iron formations (BIFs) to track the presence of Cr(VI) in Precambrian oceans, providing a time-resolved picture of the oxygenation history of the Earth's atmosphere-hydrosphere system. The geochemical behaviour of Cr is highly sensitive to the redox state of the surface environment because oxidative weathering processes produce the oxidized hexavalent [Cr(VI)] form. Oxidation of reduced trivalent [Cr(III)] chromium on land is accompanied by an isotopic fractionation, leading to enrichment of the mobile hexavalent form in the heavier isotope. Our fractionated Cr isotope data indicate the accumulation of Cr(VI) in ocean surface waters approximately 2.8 to 2.6 Gyr ago and a likely transient elevation in atmospheric and surface ocean oxygenation before the first great rise of oxygen 2.45-2.2 Gyr ago (the Great Oxidation Event). In approximately 1.88-Gyr-old BIFs we find that Cr isotopes are not fractionated, indicating a decline in atmospheric oxygen. Our findings suggest that the Great Oxidation Event did not lead to a unidirectional stepwise increase in atmospheric oxygen. In the late Neoproterozoic, we observe strong positive fractionations in Cr isotopes (delta(53)Cr up to +4.9 per thousand), providing independent support for increased surface oxygenation at that time, which may have stimulated rapid evolution of macroscopic multicellular life.

Episodic intensification of marine phosphorus burial over the last 80 million years.

Marine phosphatization events cause episodic carbonate fluorapatite (CFA) precipitation on seamounts, and are commonly linked to growth hiatuses in ferromanganese (Fe-Mn) crusts. However, the complete record of these events and their paleoenvironmental significance remains poorly understood, in large part due to poor age constraints. Here, we apply U-Pb dating to CFA in Fe-Mn crusts from Western Pacific seamounts. These data exhibit good alignment with Sr isotope ages, revealing six potential phosphatization events. This established CFA chronology tightens the timespan of phosphatization events and refines the age framework of Fe-Mn crusts. We subsequently utilize a multiproxy approach to demonstrate that the phosphatization events occurred coeval with the expansion of oceanic oxygen minimum zones. The Western Pacific Fe-Mn crusts thus document major perturbations in global oceanic phosphorus cycling, which appear to have been driven by climate-induced increases in primary productivity linked to changes in global ocean circulation.

Geographic range of plants drives long-term climate change.

Long computation times in vegetation and climate models hamper our ability to evaluate the potentially powerful role of plants on weathering and carbon sequestration over the Phanerozoic Eon. Simulated vegetation over deep time is often homogenous, and disregards the spatial distribution of plants and the impact of local climatic variables on plant function. Here we couple a fast vegetation model (FLORA) to a spatially-resolved long-term climate-biogeochemical model (SCION), to assess links between plant geographical range, the long-term carbon cycle and climate. Model results show lower rates of carbon fixation and up to double the previously predicted atmospheric CO2 concentration due to a limited plant geographical range over the arid Pangea supercontinent. The Mesozoic dispersion of the continents increases modelled plant geographical range from 65% to > 90%, amplifying global CO2 removal, consistent with geological data. We demonstrate that plant geographical range likely exerted a major, under-explored control on long-term climate change.

Dynamic redox and nutrient cycling response to climate forcing in the Mesoproterozoic ocean.

Controls on Mesoproterozoic ocean redox heterogeneity, and links to nutrient cycling and oxygenation feedbacks, remain poorly resolved. Here, we report ocean redox and phosphorus cycling across two high-resolution sections from the ~1.4 Ga Xiamaling Formation, North China Craton. In the lower section, fluctuations in trade wind intensity regulated the spatial extent of a ferruginous oxygen minimum zone, promoting phosphorus drawdown and persistent oligotrophic conditions. In the upper section, high but variable continental chemical weathering rates led to periodic fluctuations between highly and weakly euxinic conditions, promoting phosphorus recycling and persistent eutrophication. Biogeochemical modeling demonstrates how changes in geographical location relative to global atmospheric circulation cells could have driven these temporal changes in regional ocean biogeochemistry. Our approach suggests that much of the ocean redox heterogeneity apparent in the Mesoproterozoic record can be explained by climate forcing at individual locations, rather than specific events or step-changes in global oceanic redox conditions.

Biological diversification linked to environmental stabilization following the Sturtian Snowball glaciation.

The body fossil and biomarker records hint at an increase in biotic complexity between the two Cryogenian Snowball Earth episodes (ca. 661 million to ≤650 million years ago). Oxygen and nutrient availability can promote biotic complexity, but nutrient (particularly phosphorus) and redox dynamics across this interval remain poorly understood. Here, we present high-resolution paleoredox and phosphorus phase association data from multiple globally distributed drill core records through the non-glacial interval. These data are first correlated regionally by litho- and chemostratigraphy, and then calibrated within a series of global chronostratigraphic frameworks. The combined data show that regional differences in postglacial redox stabilization were partly controlled by the intensity of phosphorus recycling from marine sediments. The apparent increase in biotic complexity followed a global transition to more stable and less reducing conditions in shallow to mid-depth marine environments and occurred within a tolerable climatic window during progressive cooling after post-Snowball super-greenhouse conditions.

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.

Pyrite mega-analysis reveals modes of anoxia through geological time.

The redox structure of the water column in anoxic basins through geological time remains poorly resolved despite its importance to biological evolution/extinction and biogeochemical cycling. Here, we provide a temporal record of bottom and pore water redox conditions by analyzing the temporal distribution and chemistry of sedimentary pyrite. We combine machine-reading techniques, applied over a large library of published literature, with statistical analysis of element concentrations in databases of sedimentary pyrite and bulk sedimentary rocks to generate a scaled analysis spanning the majority of Earth's history. This analysis delineates the prevalent anoxic basin states from the Archaean to present day, which are associated with diagnostic combinations of five types of syngenetic pyrite. The underlying driver(s) for the pyrite types are unresolved but plausibly includes the ambient seawater inventory, precipitation kinetics, and the (co)location of organic matter degradation coupled to sulfate reduction, iron (oxyhydr)oxide dissolution, and pyrite precipitation.

The transition to a sulphidic ocean approximately 1.84 billion years ago.

The Proterozoic aeon (2.5 to 0.54 billion years (Gyr) ago) marks the time between the largely anoxic world of the Archean (> 2.5 Gyr ago) and the dominantly oxic world of the Phanerozoic (< 0.54 Gyr ago). The course of ocean chemistry through the Proterozoic has traditionally been explained by progressive oxygenation of the deep ocean in response to an increase in atmospheric oxygen around 2.3 Gyr ago. This postulated rise in the oxygen content of the ocean is in turn thought to have led to the oxidation of dissolved iron, Fe(II), thus ending the deposition of banded iron formations (BIF) around 1.8 Gyr ago. An alternative interpretation suggests that the increasing atmospheric oxygen levels enhanced sulphide weathering on land and the flux of sulphate to the oceans. This increased rates of sulphate reduction, resulting in Fe(II) removal in the form of pyrite as the oceans became sulphidic. Here we investigate sediments from the approximately 1.8-Gyr-old Animikie group, Canada, which were deposited during the final stages of the main global period of BIF deposition. This allows us to evaluate the two competing hypotheses for the termination of BIF deposition. We use iron-sulphur-carbon (Fe-S-C) systematics to demonstrate continued ocean anoxia after the final global deposition of BIF and show that a transition to sulphidic bottom waters was ultimately responsible for the termination of BIF deposition. Sulphidic conditions may have persisted until a second major rise in oxygen between 0.8 to 0.58 Gyr ago, possibly reducing global rates of primary production and arresting the pace of algal evolution.

Stepwise Earth oxygenation is an inherent property of global biogeochemical cycling.

Oxygenation of Earth's atmosphere and oceans occurred across three major steps during the Paleoproterozoic, Neoproterozoic, and Paleozoic eras, with each increase having profound consequences for the biosphere. Biological or tectonic revolutions have been proposed to explain each of these stepwise increases in oxygen, but the principal driver of each event remains unclear. Here we show, using a theoretical model, that the observed oxygenation steps are a simple consequence of internal feedbacks in the long-term biogeochemical cycles of carbon, oxygen, and phosphorus, and that there is no requirement for a specific stepwise external forcing to explain the course of Earth surface oxygenation. We conclude that Earth's oxygenation events are entirely consistent with gradual oxygenation of the planetary surface after the evolution of oxygenic photosynthesis.

Possible links between extreme oxygen perturbations and the Cambrian radiation of animals.

The role of oxygen as a driver for early animal evolution is widely debated. During the Cambrian explosion, episodic radiations of major animal phyla occurred coincident with repeated carbon isotope fluctuations. However, the driver of these isotope fluctuations and potential links to environmental oxygenation are unclear. Here, we report high-resolution carbon and sulphur isotope data for marine carbonates from the southeastern Siberian Platform that document the canonical explosive phase of the Cambrian radiation from ~524 to ~514 Myr ago. These analyses demonstrate a strong positive covariation between carbonate δ13C and carbonate-associated sulphate δ34S through five isotope cycles. Biogeochemical modelling suggests that this isotopic coupling reflects periodic oscillations in atmospheric O2 and the extent of shallow ocean oxygenation. Episodic maxima in the biodiversity of animal phyla directly coincided with these extreme oxygen perturbations. Conversely, the subsequent Botoman-Toyonian animal extinction events (~514 to ~512 Myr ago) coincided with decoupled isotope records that suggest a shrinking marine sulphate reservoir and expanded shallow marine anoxia. We suggest that fluctuations in oxygen availability in the shallow marine realm exerted a primary control on the timing and tempo of biodiversity radiations at a crucial phase in the early history of animal life.

Early Palaeozoic ocean anoxia and global warming driven by the evolution of shallow burrowing.

The evolution of burrowing animals forms a defining event in the history of the Earth. It has been hypothesised that the expansion of seafloor burrowing during the Palaeozoic altered the biogeochemistry of the oceans and atmosphere. However, whilst potential impacts of bioturbation on the individual phosphorus, oxygen and sulphur cycles have been considered, combined effects have not been investigated, leading to major uncertainty over the timing and magnitude of the Earth system response to the evolution of bioturbation. Here we integrate the evolution of bioturbation into the COPSE model of global biogeochemical cycling, and compare quantitative model predictions to multiple geochemical proxies. Our results suggest that the advent of shallow burrowing in the early Cambrian contributed to a global low-oxygen state, which prevailed for ~100 million years. This impact of bioturbation on global biogeochemistry likely affected animal evolution through expanded ocean anoxia, high atmospheric CO2 levels and global warming.

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

Marine oxygen production and open water supported an active nitrogen cycle during the Marinoan Snowball Earth.

The Neoproterozoic Earth was punctuated by two low-latitude Snowball Earth glaciations. Models permit oceans with either total ice cover or substantial areas of open water. Total ice cover would make an anoxic ocean likely, and would be a formidable barrier to biologic survival. However, there are no direct data constraining either the redox state of the ocean or marine biological productivity during the glacials. Here we present iron-speciation, redox-sensitive trace element, and nitrogen isotope data from a Neoproterozoic (Marinoan) glacial episode. Iron-speciation indicates deeper waters were anoxic and Fe-rich, while trace element concentrations indicate surface waters were in contact with an oxygenated atmosphere. Furthermore, synglacial sedimentary nitrogen is isotopically heavier than the modern atmosphere, requiring a biologic cycle with nitrogen fixation, nitrification and denitrification. Our results indicate significant regions of open marine water and active biologic productivity throughout one of the harshest glaciations in Earth history.