One-third of Southern Ocean productivity is supported by dust deposition.

Natural iron fertilization of the Southern Ocean by windblown dust has been suggested to enhance biological productivity and modulate the climate1-3. Yet, this process has never been quantified across the Southern Ocean and at annual timescales4,5. Here we combined 11 years of nitrate observations from autonomous biogeochemical ocean profiling floats with a Southern Hemisphere dust simulation to empirically derive the relationship between dust-iron deposition and annual net community production (ANCP) in the iron-limited Southern Ocean. Using this relationship, we determined the biological response to dust-iron in the pelagic perennially ice-free Southern Ocean at present and during the last glacial maximum (LGM). We estimate that dust-iron now supports 33% ± 15% of Southern Ocean ANCP. During the LGM, when dust deposition was 5-40-fold higher than today, the contribution of dust to Southern Ocean ANCP was much greater, estimated at 64% ± 13%. We provide quantitative evidence of basin-wide dust-iron fertilization of the Southern Ocean and the potential magnitude of its impact on glacial-interglacial timescales, supporting the idea of the important role of dust in the global carbon cycle and climate6-8.

Rhizobia-diatom symbiosis fixes missing nitrogen in the ocean.

Nitrogen (N2) fixation in oligotrophic surface waters is the main source of new nitrogen to the ocean1 and has a key role in fuelling the biological carbon pump2. Oceanic N2 fixation has been attributed almost exclusively to cyanobacteria, even though genes encoding nitrogenase, the enzyme that fixes N2 into ammonia, are widespread among marine bacteria and archaea3-5. Little is known about these non-cyanobacterial N2 fixers, and direct proof that they can fix nitrogen in the ocean has so far been lacking. Here we report the discovery of a non-cyanobacterial N2-fixing symbiont, 'Candidatus Tectiglobus diatomicola', which provides its diatom host with fixed nitrogen in return for photosynthetic carbon. The N2-fixing symbiont belongs to the order Rhizobiales and its association with a unicellular diatom expands the known hosts for this order beyond the well-known N2-fixing rhizobia-legume symbioses on land6. Our results show that the rhizobia-diatom symbioses can contribute as much fixed nitrogen as can cyanobacterial N2 fixers in the tropical North Atlantic, and that they might be responsible for N2 fixation in the vast regions of the ocean in which cyanobacteria are too rare to account for the measured rates.

Abyssal ocean overturning slowdown and warming driven by Antarctic meltwater.

The abyssal ocean circulation is a key component of the global meridional overturning circulation, cycling heat, carbon, oxygen and nutrients throughout the world ocean1,2. The strongest historical trend observed in the abyssal ocean is warming at high southern latitudes2-4, yet it is unclear what processes have driven this warming, and whether this warming is linked to a slowdown in the ocean's overturning circulation. Furthermore, attributing change to specific drivers is difficult owing to limited measurements, and because coupled climate models exhibit biases in the region5-7. In addition, future change remains uncertain, with the latest coordinated climate model projections not accounting for dynamic ice-sheet melt. Here we use a transient forced high-resolution coupled ocean-sea-ice model to show that under a high-emissions scenario, abyssal warming is set to accelerate over the next 30 years. We find that meltwater input around Antarctica drives a contraction of Antarctic Bottom Water (AABW), opening a pathway that allows warm Circumpolar Deep Water greater access to the continental shelf. The reduction in AABW formation results in warming and ageing of the abyssal ocean, consistent with recent measurements. In contrast, projected wind and thermal forcing has little impact on the properties, age and volume of AABW. These results highlight the critical importance of Antarctic meltwater in setting the abyssal ocean overturning, with implications for global ocean biogeochemistry and climate that could last for centuries.

Authigenic mineral phases as a driver of the upper-ocean iron cycle.

Iron is important in regulating the ocean carbon cycle1. Although several dissolved and particulate species participate in oceanic iron cycling, current understanding emphasizes the importance of complexation by organic ligands in stabilizing oceanic dissolved iron concentrations2-6. However, it is difficult to reconcile this view of ligands as a primary control on dissolved iron cycling with the observed size partitioning of dissolved iron species, inefficient dissolved iron regeneration at depth or the potential importance of authigenic iron phases in particulate iron observational datasets7-12. Here we present a new dissolved iron, ligand and particulate iron seasonal dataset from the Bermuda Atlantic Time-series Study (BATS) region. We find that upper-ocean dissolved iron dynamics were decoupled from those of ligands, which necessitates a process by which dissolved iron escapes ligand stabilization to generate a reservoir of authigenic iron particles that settle to depth. When this 'colloidal shunt' mechanism was implemented in a global-scale biogeochemical model, it reproduced both seasonal iron-cycle dynamics observations and independent global datasets when previous models failed13-15. Overall, we argue that the turnover of authigenic particulate iron phases must be considered alongside biological activity and ligands in controlling ocean-dissolved iron distributions and the coupling between dissolved and particulate iron pools.

Natural short-lived halogens exert an indirect cooling effect on climate.

Observational evidence shows the ubiquitous presence of ocean-emitted short-lived halogens in the global atmosphere1-3. Natural emissions of these chemical compounds have been anthropogenically amplified since pre-industrial times4-6, while, in addition, anthropogenic short-lived halocarbons are currently being emitted to the atmosphere7,8. Despite their widespread distribution in the atmosphere, the combined impact of these species on Earth's radiative balance remains unknown. Here we show that short-lived halogens exert a substantial indirect cooling effect at present (-0.13 ± 0.03 watts per square metre) that arises from halogen-mediated radiative perturbations of ozone (-0.24 ± 0.02 watts per square metre), compensated by those from methane (+0.09 ± 0.01 watts per square metre), aerosols (+0.03 ± 0.01 watts per square metre) and stratospheric water vapour (+0.011 ± 0.001 watts per square metre). Importantly, this substantial cooling effect has increased since 1750 by -0.05 ± 0.03 watts per square metre (61 per cent), driven by the anthropogenic amplification of natural halogen emissions, and is projected to change further (18-31 per cent by 2100) depending on climate warming projections and socioeconomic development. We conclude that the indirect radiative effect due to short-lived halogens should now be incorporated into climate models to provide a more realistic natural baseline of Earth's climate system.

A well-oxygenated eastern tropical Pacific during the warm Miocene.

The oxygen content of the oceans is susceptible to climate change and has declined in recent decades1, with the largest effect in oxygen-deficient zones (ODZs)2, that is, mid-depth ocean regions with oxygen concentrations <5 µmol kg-1 (ref. 3). Earth-system-model simulations of climate warming predict that ODZs will expand until at least 2100. The response on timescales of hundreds to thousands of years, however, remains uncertain3-5. Here we investigate changes in the response of ocean oxygenation during the warmer-than-present Miocene Climatic Optimum (MCO; 17.0-14.8 million years ago (Ma)). Our planktic foraminifera I/Ca and δ15N data, palaeoceanographic proxies sensitive to ODZ extent and intensity, indicate that dissolved-oxygen concentrations in the eastern tropical Pacific (ETP) exceeded 100 µmol kg-1 during the MCO. Paired Mg/Ca-derived temperature data suggest that an ODZ developed in response to an increased west-to-east temperature gradient and shoaling of the ETP thermocline. Our records align with model simulations of data from recent decades to centuries6,7, suggesting that weaker equatorial Pacific trade winds during warm periods may lead to decreased upwelling in the ETP, causing equatorial productivity and subsurface oxygen demand to be less concentrated in the east. These findings shed light on how warm-climate states such as during the MCO may affect ocean oxygenation. If the MCO is considered as a possible analogue for future warming, our findings seem to support models suggesting that the recent deoxygenation trend and expansion of the ETP ODZ may eventually reverse3,4.

Coral reefs benefit from reduced land-sea impacts under ocean warming.

Coral reef ecosystems are being fundamentally restructured by local human impacts and climate-driven marine heatwaves that trigger mass coral bleaching and mortality1. Reducing local impacts can increase reef resistance to and recovery from bleaching2. However, resource managers lack clear advice on targeted actions that best support coral reefs under climate change3 and sector-based governance means most land- and sea-based management efforts remain siloed4. Here we combine surveys of reef change with a unique 20-year time series of land-sea human impacts that encompassed an unprecedented marine heatwave in Hawai'i. Reefs with increased herbivorous fish populations and reduced land-based impacts, such as wastewater pollution and urban runoff, had positive coral cover trajectories predisturbance. These reefs also experienced a modest reduction in coral mortality following severe heat stress compared to reefs with reduced fish populations and enhanced land-based impacts. Scenario modelling indicated that simultaneously reducing land-sea human impacts results in a three- to sixfold greater probability of a reef having high reef-builder cover four years postdisturbance than if either occurred in isolation. International efforts to protect 30% of Earth's land and ocean ecosystems by 2030 are underway5. Our results reveal that integrated land-sea management could help achieve coastal ocean conservation goals and provide coral reefs with the best opportunity to persist in our changing climate.

Biological carbon pump estimate based on multidecadal hydrographic data.

The transfer of photosynthetically produced organic carbon from surface to mesopelagic waters draws carbon dioxide from the atmosphere1. However, current observation-based estimates disagree on the strength of this biological carbon pump (BCP)2. Earth system models (ESMs) also exhibit a large spread of BCP estimates, indicating limited representations of the known carbon export pathways3. Here we use several decades of hydrographic observations to produce a top-down estimate of the strength of the BCP with an inverse biogeochemical model that implicitly accounts for all known export pathways. Our estimate of total organic carbon (TOC) export at 73.4 m (model euphotic zone depth) is 15.00 ± 1.12 Pg C year-1, with only two-thirds reaching 100 m depth owing to rapid remineralization of organic matter in the upper water column. Partitioned by sequestration time below the euphotic zone, τ, the globally integrated organic carbon production rate with τ > 3 months is 11.09 ± 1.02 Pg C year-1, dropping to 8.25 ± 0.30 Pg C year-1 for τ > 1 year, with 81% contributed by the non-advective-diffusive vertical flux owing to sinking particles and vertically migrating zooplankton. Nevertheless, export of organic carbon by mixing and other fluid transport of dissolved matter and suspended particles remains regionally important for meeting the respiratory carbon demand. Furthermore, the temperature dependence of the sequestration efficiency inferred from our inversion suggests that future global warming may intensify the recycling of organic matter in the upper ocean, potentially weakening the BCP.

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.

Agricultural pesticide land budget and river discharge to oceans.

Pesticides are ubiquitous environmental pollutants negatively affecting ecosystem and human health1,2. About 3 Tg of pesticides are used annually in agriculture to protect crops3. How much of these pesticides remain on land and reach the aquifer or the ocean is uncertain. Monitoring their environmental fate is challenging, and a detailed picture of their mobility in time and space is largely missing4. Here, we develop a process-based model accounting for the hydrology and biogeochemistry of the 92 most used agricultural pesticide active substances to assess their pathways through the principal catchments of the world and draw a near-present picture of the global land and river budgets, including discharge to oceans. Of the 0.94 Tg net annual pesticide input in 2015 used in this study, 82% is biologically degraded, 10% remains as residue in soil and 7.2% leaches below the root zone. Rivers receive 0.73 Gg of pesticides from their drainage at a rate of 10 to more than 100 kg yr-1 km-1. By contrast to their fate in soil, only 1.1% of pesticides entering rivers are degraded along streams, exceeding safety levels (concentrations >1 µg l-1) in more than 13,000 km of river length, with 0.71 Gg of pesticide active ingredients released to oceans every year. Herbicides represent the prevalent pesticide residue on both land (72%) and river outlets (62%).

The evolution of the marine carbonate factory.

Calcium carbonate formation is the primary pathway by which carbon is returned from the ocean-atmosphere system to the solid Earth1,2. The removal of dissolved inorganic carbon from seawater by precipitation of carbonate minerals-the marine carbonate factory-plays a critical role in shaping marine biogeochemical cycling1,2. A paucity of empirical constraints has led to widely divergent views on how the marine carbonate factory has changed over time3-5. Here we use geochemical insights from stable strontium isotopes to provide a new perspective on the evolution of the marine carbonate factory and carbonate mineral saturation states. Although the production of carbonates in the surface ocean and in shallow seafloor settings have been widely considered the predominant carbonate sinks for most of the history of the Earth6, we propose that alternative processes-such as porewater production of authigenic carbonates-may have represented a major carbonate sink throughout the Precambrian. Our results also suggest that the rise of the skeletal carbonate factory decreased seawater carbonate saturation states.

Observed poleward freshwater transport since 1970.

Warming-induced global water cycle changes pose a significant challenge to global ecosystems and human society. However, quantifying historical water cycle change is difficult owing to a dearth of direct observations, particularly over the ocean, where 77% and 85% of global precipitation and evaporation occur, respectively1-3. Air-sea fluxes of freshwater imprint on ocean salinity such that mean salinity is lowest in the warmest and coldest parts of the ocean, and is highest at intermediate temperatures4. Here we track salinity trends in the warm, salty fraction of the ocean, and quantify the observed net poleward transport of freshwater in the Earth system from 1970 to 2014. Over this period, poleward freshwater transport from warm to cold ocean regions has occurred at a rate of 34-62 milli-sverdrups (mSv = 103 m3 s-1), a rate that is not replicated in the current generation of climate models (the Climate Model Intercomparison Project Phase 6 (CMIP6)). In CMIP6 models, surface freshwater flux intensification in warm ocean regions leads to an approximately equivalent change in ocean freshwater content, with little impact from ocean mixing and circulation. Should this partition of processes hold for the real world, the implication is that the historical surface flux amplification is weaker (0.3-4.6%) in CMIP6 compared with observations (3.0-7.4%). These results establish a historical constraint on poleward freshwater transport that will assist in addressing biases in climate models.

Arctic Ocean annual high in [Formula: see text] could shift from winter to summer.

Long-term stress on marine organisms from ocean acidification will differ between seasons. As atmospheric carbon dioxide (CO2) increases, so do seasonal variations of ocean CO2 partial pressure ([Formula: see text]), causing summer and winter long-term trends to diverge1-5. Trends may be further influenced by an unexplored factor-changes in the seasonal timing of [Formula: see text]. In Arctic Ocean surface waters, the observed timing is typified by a winter high and summer low6 because biological effects dominate thermal effects. Here we show that 27 Earth system models simulate similar timing under historical forcing but generally project that the summer low, relative to the annual mean, eventually becomes a high across much of the Arctic Ocean under mid-to-high-level CO2 emissions scenarios. Often the greater increase in summer [Formula: see text], although gradual, abruptly inverses the chronological order of the annual high and low, a phenomenon not previously seen in climate-related variables. The main cause is the large summer sea surface warming7 from earlier retreat of seasonal sea ice8. Warming and changes in other drivers enhance this century's increase in extreme summer [Formula: see text] by 29 ± 9 per cent compared with no change in driver seasonalities. Thus the timing change worsens summer ocean acidification, which in turn may lower the tolerance of endemic marine organisms to increasing summer temperatures.

Volcanic trigger of ocean deoxygenation during Cordilleran ice sheet retreat.

North Pacific deoxygenation events during the last deglaciation were sustained over millennia by high export productivity, but the triggering mechanisms and their links to deglacial warming remain uncertain1-3. Here we find that initial deoxygenation in the North Pacific immediately after the Cordilleran ice sheet (CIS) retreat4 was associated with increased volcanic ash in seafloor sediments. Timing of volcanic inputs relative to CIS retreat suggests that regional explosive volcanism was initiated by ice unloading5,6. We posit that iron fertilization by volcanic ash7-9 during CIS retreat fuelled ocean productivity in this otherwise iron-limited region, and tipped the marine system towards sustained deoxygenation. We also identify older deoxygenation events linked to CIS retreat over the past approximately 50,000 years (ref. 4). Our findings suggest that the apparent coupling between the atmosphere, ocean, cryosphere and solid-Earth systems occurs on relatively short timescales and can act as an important driver for ocean biogeochemical change.

Enhanced ocean oxygenation during Cenozoic warm periods.

Dissolved oxygen (O2) is essential for most ocean ecosystems, fuelling organisms' respiration and facilitating the cycling of carbon and nutrients. Oxygen measurements have been interpreted to indicate that the ocean's oxygen-deficient zones (ODZs) are expanding under global warming1,2. However, models provide an unclear picture of future ODZ change in both the near term and the long term3-6. The paleoclimate record can help explore the possible range of ODZ changes in warmer-than-modern periods. Here we use foraminifera-bound nitrogen (N) isotopes to show that water-column denitrification in the eastern tropical North Pacific was greatly reduced during the Middle Miocene Climatic Optimum (MMCO) and the Early Eocene Climatic Optimum (EECO). Because denitrification is restricted to oxygen-poor waters, our results indicate that, in these two Cenozoic periods of sustained warmth, ODZs were contracted, not expanded. ODZ contraction may have arisen from a decrease in upwelling-fuelled biological productivity in the tropical Pacific, which would have reduced oxygen demand in the subsurface. Alternatively, invigoration of deep-water ventilation by the Southern Ocean may have weakened the ocean's 'biological carbon pump', which would have increased deep-ocean oxygen. The mechanism at play would have determined whether the ODZ contractions occurred in step with the warming or took centuries or millennia to develop. Thus, although our results from the Cenozoic do not necessarily apply to the near-term future, they might imply that global warming may eventually cause ODZ contraction.

Pliocene decoupling of equatorial Pacific temperature and pH gradients.

Ocean dynamics in the equatorial Pacific drive tropical climate patterns that affect marine and terrestrial ecosystems worldwide. How this region will respond to global warming has profound implications for global climate, economic stability and ecosystem health. As a result, numerous studies have investigated equatorial Pacific dynamics during the Pliocene (5.3-2.6 million years ago) and late Miocene (around 6 million years ago) as an analogue for the future behaviour of the region under global warming1-12. Palaeoceanographic records from this time present an apparent paradox with proxy evidence of a reduced east-west sea surface temperature gradient along the equatorial Pacific1,3,7,8-indicative of reduced wind-driven upwelling-conflicting with evidence of enhanced biological productivity in the east Pacific13-15 that typically results from stronger upwelling. Here we reconcile these observations by providing new evidence for a radically different-from-modern circulation regime in the early Pliocene/late Miocene16 that results in older, more acidic and more nutrient-rich water reaching the equatorial Pacific. These results provide a mechanism for enhanced productivity in the early Pliocene/late Miocene east Pacific even in the presence of weaker wind-driven upwelling. Our findings shed new light on equatorial Pacific dynamics and help to constrain the potential changes they will undergo in the near future, given that the Earth is expected to reach Pliocene-like levels of warming in the next century.

Summertime increases in upper-ocean stratification and mixed-layer depth.

The surface mixed layer of the world ocean regulates global climate by controlling heat and carbon exchange between the atmosphere and the oceanic interior1-3. The mixed layer also shapes marine ecosystems by hosting most of the ocean's primary production4 and providing the conduit for oxygenation of deep oceanic layers. Despite these important climatic and life-supporting roles, possible changes in the mixed layer during an era of global climate change remain uncertain. Here we use oceanographic observations to show that from 1970 to 2018 the density contrast across the base of the mixed layer increased and that the mixed layer itself became deeper. Using a physically based definition of upper-ocean stability that follows different dynamical regimes across the global ocean, we find that the summertime density contrast increased by 8.9 ± 2.7 per cent per decade (10-6-10-5 per second squared per decade, depending on region), more than six times greater than previous estimates. Whereas prior work has suggested that a thinner mixed layer should accompany a more stratified upper ocean5-7, we find instead that the summertime mixed layer deepened by 2.9 ± 0.5 per cent per decade, or several metres per decade (typically 5-10 metres per decade, depending on region). A detailed mechanistic interpretation is challenging, but the concurrent stratification and deepening of the mixed layer are related to an increase in stability associated with surface warming and high-latitude surface freshening8,9, accompanied by a wind-driven intensification of upper-ocean turbulence10,11. Our findings are based on a complex dataset with incomplete coverage of a vast area. Although our results are robust within a wide range of sensitivity analyses, important uncertainties remain, such as those related to sparse coverage in the early years of the 1970-2018 period. Nonetheless, our work calls for reconsideration of the drivers of ongoing shifts in marine primary production, and reveals stark changes in the world's upper ocean over the past five decades.

A lithium-isotope perspective on the evolution of carbon and silicon cycles.

The evolution of the global carbon and silicon cycles is thought to have contributed to the long-term stability of Earth's climate1-3. Many questions remain, however, regarding the feedback mechanisms at play, and there are limited quantitative constraints on the sources and sinks of these elements in Earth's surface environments4-12. Here we argue that the lithium-isotope record can be used to track the processes controlling the long-term carbon and silicon cycles. By analysing more than 600 shallow-water marine carbonate samples from more than 100 stratigraphic units, we construct a new carbonate-based lithium-isotope record spanning the past 3 billion years. The data suggest an increase in the carbonate lithium-isotope values over time, which we propose was driven by long-term changes in the lithium-isotopic conditions of sea water, rather than by changes in the sedimentary alterations of older samples. Using a mass-balance modelling approach, we propose that the observed trend in lithium-isotope values reflects a transition from Precambrian carbon and silicon cycles to those characteristic of the modern. We speculate that this transition was linked to a gradual shift to a biologically controlled marine silicon cycle and the evolutionary radiation of land plants13,14.

Seasonal origin of the thermal maxima at the Holocene and the last interglacial.

Proxy reconstructions from marine sediment cores indicate peak temperatures in the first half of the last and current interglacial periods (the thermal maxima of the Holocene epoch, 10,000 to 6,000 years ago, and the last interglacial period, 128,000 to 123,000 years ago) that arguably exceed modern warmth1-3. By contrast, climate models simulate monotonic warming throughout both periods4-7. This substantial model-data discrepancy undermines confidence in both proxy reconstructions and climate models, and inhibits a mechanistic understanding of recent climate change. Here we show that previous global reconstructions of temperature in the Holocene1-3 and the last interglacial period8 reflect the evolution of seasonal, rather than annual, temperatures and we develop a method of transforming them to mean annual temperatures. We further demonstrate that global mean annual sea surface temperatures have been steadily increasing since the start of the Holocene (about 12,000 years ago), first in response to retreating ice sheets (12 to 6.5 thousand years ago), and then as a result of rising greenhouse gas concentrations (0.25 ± 0.21 degrees Celsius over the past 6,500 years or so). However, mean annual temperatures during the last interglacial period were stable and warmer than estimates of temperatures during the Holocene, and we attribute this to the near-constant greenhouse gas levels and the reduced extent of ice sheets. We therefore argue that the climate of the Holocene differed from that of the last interglacial period in two ways: first, larger remnant glacial ice sheets acted to cool the early Holocene, and second, rising greenhouse gas levels in the late Holocene warmed the planet. Furthermore, our reconstructions demonstrate that the modern global temperature has exceeded annual levels over the past 12,000 years and probably approaches the warmth of the last interglacial period (128,000 to 115,000 years ago).

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.