Changes In Ocean Circulation Research Articles

An optimal transformation method for inferring ocean tracer sources and sinks

Abstract. The geography of changes in the fluxes of heat, carbon, freshwater and other tracers at the sea surface is highly uncertain and is critical to our understanding of climate change and its impacts. We present a state estimation framework wherein prior estimates of boundary fluxes can be adjusted to make them consistent with the evolving ocean state. In this framework, we define a discrete set of ocean water masses distinguished by their geographical, thermodynamic and chemical properties for specific time periods. Ocean circulation then moves these water masses in geographic space. In phase space, geographically adjacent water masses are able to mix together, representing a convergence, and air–sea property fluxes move the water masses over time. We define an optimisation problem whose solution is constrained by the physically permissible bounds of changes in ocean circulation, air–sea fluxes and mixing. As a proof-of-concept implementation, we use data from a historical numerical climate model simulation with a closed heat and salinity budget. An inverse model solution is found for the evolution of temperature and salinity that is consistent with “true” air–sea heat and freshwater fluxes which are introduced as model priors. When biases are introduced into the prior fluxes, the inverse model finds a solution closer to the true fluxes. This framework, which we call the optimal transformation method, represents a modular, relatively computationally cost-effective, open-source and transparent state estimation tool that complements existing approaches.

Revisiting the mid-Pleistocene transition ocean circulation crisis.

The mid-Pleistocene transition (MPT) [~1.25 to 0.85 million years ago (Ma)] marks a shift in the character of glacial-interglacial climate (1, 2). One prevailing hypothesis for the origin of the MPT is that glacial deep ocean circulation fundamentally changed, marked by a circulation "crisis" at ~0.90 Ma (marine isotope stages 24 to 22) (3). Using high-resolution paired neodymium, carbon, and oxygen isotope data from the South Atlantic Ocean (Cape Basin) across the MPT, we find no evidence of a substantial change in deep ocean circulation. Before and during the early MPT (~1.30 to 1.12 Ma), the glacial deep ocean variability closely resembled that of the most recent glacial cycle. The carbon storage facilitated by developing deep ocean stratification across the MPT required only modest circulation adjustments.

The role of ocean circulation and regolith removal in triggering the Mid-Pleistocene Transition: Insights from authigenic Nd isotopes

Approximately 1,250,000 to 700,000 years ago, the pacing of glacial-interglacial cycles changed from 41,000 years to ∼100,000 years, a shift known as the ‘Mid-Pleistocene Transition’ (MPT). The cause – or causes – of this shift remain uncertain. However, changes in ocean circulation and removal of northern hemisphere regolith have both been proposed as potential triggers. Here, we present continuous, orbitally resolved reconstructions of deep ocean neodymium isotopes from three locations in the equatorial Atlantic, Indian and southwest Pacific oceans, spanning 1.7 million years from the Holocene to before the MPT, to test these two hypotheses. We find that global seawater neodymium isotope variability over glacial-interglacial cycles is controlled by changes in both neodymium input to the North Atlantic and deep ocean mixing. Using this neodymium isotope data, we show that enhanced northern hemisphere regolith removal began approximately 1.45 million years ago, ∼200,000 years prior to the onset of the MPT and ∼500,000 years prior to a major expansion in northern hemisphere ice sheets between ∼900,000 and 870,000 years ago. This ice sheet expansion was immediately preceded by an interval of reduced mixing of Atlantic-sourced waters into the deep southwest Pacific Ocean. We hypothesize that this circulation reorganization reflected increased stratification of the deep Southern Ocean interior, possibly in response to cooling and Antarctic sea ice expansion at this time. Taken together, these data suggest an expansion and/or thickening of northern hemisphere ice sheets during the MPT was facilitated by a combination of reduced northern hemisphere regolith cover alongside Southern Ocean circulation changes. Together, these shifts allowed the build up of larger northern hemisphere ice sheets that were more resistant to deglaciation, facilitating the longer glacial cycles of the post-MPT world.

Uncovering the oxic-suboxic microenvironment change in seamount flank through authigenic clay minerals in basaltic substrate of ferromanganese crust, Magellan seamount

In low-temperature ocean environments, basalt alteration by seawater precipitates authigenic clay minerals that serve as proxies for reconstructing paleo-ocean conditions because they reflect surrounding oxic-suboxic conditions. In this study, alteration rims on basaltic substrate associated with ferromanganese (Fe-Mn) crust from the Magellan seamount KC-7 were identified by microscopic analyses. Mineralogical and geochemical analyses indicate that the alteration rims contain K-enriched Fe-smectite and glauconite which suggest that seawater-basalt interaction occurred under oxic conditions and in the presence of organic-rich suboxic conditions, respectively. These disparate environmental conditions suggest that the environment changed before and after Fe-Mn crusts formed. During the Cenozoic hyperthermal events, oxygen-rich bottom water was supplied by upwelling driven by the geomorphological influence of the seamounts, which may have led to basalt alteration. The K-enriched Fe-smectites, which indicate oxic condition, formed via seawater-basalt interactions before the Fe-Mn crust incrustation. Later, during the Eocene, the opening of the Drake Passage enhanced the supply of oxygen-rich seawater to the Magellan Seamounts, thereby enabling the formation of hydrogenetic Fe-Mn crust. After the incrustation of seamount flanks with Fe-Mn crusts, the carbonate fluorapatite (CFA), a product of the global phosphatization event, filled the pores in the Fe-Mn crusts during Oi-1 glaciation. As a result, seawater-basalt interactions decreased and led to suboxic conditions, in which glauconite formed as organic matter was remineralized under the organic-rich conditions in the basaltic substrate. This authigenic clay mineral formation sequence suggests that changes in ocean circulation and subsequent changes in the oxic-suboxic conditions in the basaltic substrate occurred on the western Pacific seamount KC-7.

Calibrating lower-middle Miocene mammal faunas and unravelling climate change during the Miocene Climate Optimum: the Bardenas Reales de Navarra record (Ebro basin, NE Iberian Peninsula)

AbstractThe chronology of lower Miocene Iberian small mammal faunas is still poorly constrained given the scarcity of well dated sedimentary successions including small mammal fossil localities. Such scarcity has prevented also an accurate understanding of the response of European terrestrial ecosystem to global changes across the Miocene climate optimum (MCO), one of the best analogues of present-day global warming. Here we present an updated fossil small mammal record of the Bardenas Reales de Navarra (western Ebro basin, Spain), where an expanded lower to middle Miocene continental succession is superbly exposed. Previous and new magnetostratigraphic results from this succession have enabled us to propose, along with additional magnetostratigraphically-dated Iberian faunas, a new chronology for local zones Y to D (Mammal Neogene zones MN2 to MN5). In addition to that, the studied small mammal faunas point to a gradual increase in temperature and humidity conditions in SW Europe between 20 and 15.5 Ma, which appears to be coupled with the progressive shift towards warmer regional (Atlantic) and global conditions across the MCO, thereby pointing to gradual changes in oceanic circulation as the main driver of this period of global warmth. The evolution of sedimentary facies appears to indicate a threshold response of the Ebro basin hydrological balance to the MCO, whereas pedogenic formation of magnetic minerals seems to be linked to periods of enhanced climate variability. These results highlight the need of combining different paleoenvironmental indicators in order to obtain a reliable view of the response of continental ecosystems to global warming.

Towards the construction of regional marine radiocarbon calibration curves: an unsupervised machine learning approach

Abstract. Radiocarbon may serve as a powerful dating tool in palaeoceanography, but its accuracy is limited by the need to calibrate radiocarbon dates to calendar ages. A key problem is that marine radiocarbon dates must be corrected for past offsets from either the contemporary atmosphere (i.e. “reservoir age” offsets) or a modelled estimate of the global average surface ocean (i.e. delta-R offsets). This presents a challenge because the spatial distribution of reservoir ages and delta-R offsets can vary significantly, particularly over periods of major marine hydrographic and/or carbon cycle change such as the last deglaciation. Modern reservoir age and delta-R estimates therefore have limited applicability. While forward modelling of past R-age variability has been proposed as a means of resolving this problem, this requires accurate a priori knowledge of past global radiocarbon budget closure (i.e. production, and cycling), which we currently lack. In this context, the construction of empirical regional marine calibration curves could provide a way forward. However, the spatial reach of such calibrations and their robustness subject to (uncertain) temporal changes in climate and ocean circulation would need to be tested. Here, we use unsupervised machine learning techniques to define distinct regions of the surface ocean that exhibit coherent behaviour in terms of their radiocarbon age offsets from the contemporary atmosphere (R ages), regardless of the causes of R-age variability. We apply multiple algorithms (k-means, k-medoids, and hierarchical clustering) to outputs from two different numerical models spanning a range of climate states, forcings, and timescales of adjustment. Comparisons between the cluster assignments across model runs confirm some robust regional patterns that likely stem from constraints imposed by large-scale ocean and atmospheric physics. At the coarsest scale, regions of coherent R-age variability correspond to the major ocean basins. By further dividing basin-scale shape-based clusters into amplitude-based subclusters, we recover regional associations, such as increased high-latitude R ages, or the propagation of R-age anomalies from regions of deep mixing in the Southern Ocean to upwelling sites in the eastern equatorial Pacific, which cohere with known modern oceanographic processes. We show that the medoids (i.e. the most representative locations) for these regional sub-clusters provide significantly better approximations of simulated local R-age variability than constant offsets from the global surface average. This remains true when cluster assignments obtained from one model simulation are applied to simulated R-age time series from another. Further, model-based clusters are found to be broadly consistent with existing reservoir age reconstructions that span the last ∼30 kyr. We therefore propose that machine learning provides a promising approach to the problem of defining regions for which empirical marine radiocarbon calibration curves may eventually be generated.

Investigating ocean circulation dynamics through data assimilation: A mathematical study using the Stommel box model with rapid oscillatory forcings.

We investigate ocean circulation changes through the lens of data assimilation using a reduced-order model. Our primary interest lies in the Stommel box model, which reveals itself to be one of the most practicable models that has the ability of reproducing the meridional overturning circulation. The Stommel box model has at most two regimes: TH (temperature driven circulation with sinking near the north pole) and SA (salinity driven with sinking near the equator). Currently, the meridional overturning is in the TH regime. Using box-averaged Met Office EN4 ocean temperature and salinity data, our goal is to provide a probability that a future regime change occurs and establish how this probability depends on the uncertainties in initial conditions, parameters, and forcings. We will achieve this using data assimilation tools and DAPPER within the Stommel box model with fast oscillatory regimes.

Is the Regime Shift in Gulf Stream Warm Core Rings Detected by Satellite Altimetry? An Inter‐Comparison of Eddy Identification and Tracking Products

AbstractDownstream of Cape Hatteras, the vigorously meandering Gulf Stream forms anticyclonic warm core rings (WCRs) that carry warm Gulf Stream and Sargasso Sea waters into the cooler, fresher Slope Sea, and forms cyclonic cold core rings (CCRs) that carry Slope Sea waters into the Sargasso Sea. The Northwest Atlantic shelf and open ocean off the U.S. East Coast have experienced dramatic changes in ocean circulation and water properties in recent years, with significant consequences for marine ecosystems and coastal communities. Some of these changes may be related to a reported regime shift in the number of WCRs formed annually, with a doubling of WCRs shed after 2000. Since the regime shift was detected using a regional eddy‐tracking product, primarily based on sea surface temperatures and relies on analyst skill, we examine three global eddy‐tracking products as an automated and potentially more objective way to detect changes in Gulf Stream rings. Currently, global products rely on altimeter‐measured sea surface height (SSH), with WCRs registering as sea surface highs and CCRs as lows. To identify eddies, these products use either SSH contours or a Lagrangian approach, with particles seeded in satellite‐based surface geostrophic velocity fields. This study confirms the three global products are not well suited for statistical analysis of Gulf Stream rings and suggests that automated WCR identification and tracking comes at the price of accurate identification and tracking. Furthermore, a shift to a higher energy state is detected in the Northwest Atlantic, which coincides with the regime shift in WCRs.

Enhanced ocean heat storage efficiency during the last deglaciation

Proxy reconstructions suggest that increasing global mean sea surface temperature (GMSST) during the last deglaciation was accompanied by a comparable or greater increase in global mean ocean temperature (GMOT), corresponding to a large heat storage efficiency (HSE; ∆GMOT/∆GMSST). An increased GMOT is commonly attributed to surface warming at sites of deepwater formation, but winter sea ice covered much of these source areas during the last deglaciation, which would imply an HSE much less than 1. Here, we use climate model simulations and proxy-based reconstructions of ocean temperature changes to show that an increased deglacial HSE is achieved by warming of intermediate-depth waters forced by mid-latitude surface warming in response to greenhouse gas and ice sheet forcing as well as by reduced Atlantic meridional overturning circulation associated with meltwater forcing. These results, which highlight the role of surface warming and oceanic circulation changes, have implications for our understanding of long-term ocean heat storage change.

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.

Role of land-ocean interactions in stepwise Northern Hemisphere Glaciation

The investigation of triggers causing the onset and intensification of Northern Hemisphere Glaciation (NHG) during the late Pliocene is essential for understanding the global climate system, with important implications for projecting future climate changes. Despite their critical roles in the global climate system, influences of land-ocean interactions on high-latitude ice sheets remain largely unexplored. Here, we present a high-resolution Asian dust record from Ocean Drilling Program Site 1208 in the North Pacific, which lies along the main route of the westerlies. Our data indicate that atmosphere-land-ocean interactions affected aeolian dust emissions through modulating moisture and vegetation in dust source regions, highlighting a critical role of terrestrial systems in initiating the NHG as early as 3.6 Myr ago. Combined with additional multi-proxy and model results, we further show that westerly wind strength was enhanced, mainly at low-to-middle tropospheric levels, during major glacial events at about 3.3 and 2.7 Myr ago. We suggest that coupled responses of Earth’s surface dynamics and atmospheric circulation in the Plio-Pleistocene likely involved feedbacks related to changes in paleogeography, ocean circulation, and global climate.

Revising chronological uncertainties in marine archives using global anthropogenic signals: a case study on the oceanic 13C Suess effect

Abstract. Marine sediments are excellent archives for reconstructing past changes in climate and ocean circulation. Overlapping with instrumental records, they hold the potential to elucidate natural variability and contextualize current changes. Yet, dating uncertainties of traditional approaches (e.g., up to ± 30–50 years for the last 2 centuries) pose major challenges for integrating the shorter instrumental records with these extended marine archives. Hence, robust sediment chronologies are crucial, and most existing age model constraints do not provide sufficient age control, particularly for the 20th century, which is the most critical period for comparing proxy records to historical changes. Here we propose a novel chronostratigraphic approach that uses anthropogenic signals such as the oceanic 13C Suess effect and spheroidal carbonaceous fly-ash particles to reduce age model uncertainties in high-resolution marine archives. As a test, we apply this new approach to a marine sediment core located at the Gardar Drift, in the subpolar North Atlantic, and revise the previously published age model for this site. We further provide a refined estimate of regional reservoir corrections and uncertainties for Gardar Drift.

Controls on Early Cretaceous South Atlantic Ocean circulation and carbon burial – a climate model–proxy synthesis

Abstract. Black shale sediments from the Barremian to Aptian South Atlantic document the intense and widespread burial of marine organic carbon during the initial stages of seafloor spreading between Africa and South America. The enhanced sequestration of atmospheric CO2 makes these young ocean basins potential drivers of the Early Cretaceous carbon cycle and climate perturbations. The opening of marine gateways between initially restricted basins and related circulation and ventilation changes are a commonly invoked explanation for the transient formation and disappearance of these regional carbon sinks. However, large uncertainties in palaeogeographic reconstructions limit the interpretation of available palaeoceanographic data and prevent any robust model-based quantifications of the proposed circulation and carbon burial changes. Here, we present a new approach to assess the principal controls on the Early Cretaceous South Atlantic and Southern Ocean circulation changes under full consideration of the uncertainties in available boundary conditions. Specifically, we use a large ensemble of 36 climate model experiments to simulate the Barremian to Albian progressive opening of the Falkland Plateau and Georgia Basin gateways with different configurations of the proto-Drake Passage, the Walvis Ridge, and atmospheric CO2 concentrations. The experiments are designed to complement available geochemical data across the regions and to test circulation scenarios derived from them. All simulations show increased evaporation and intermediate water formation at subtropical latitudes that drive a meridional overturning circulation whose vertical extent is determined by the sill depth of the Falkland Plateau. The densest water masses formed in the southern Angola Basin and potentially reached the deep Cape Basin as Walvis Ridge Overflow Water. Palaeogeographic uncertainties are as important as the lack of precise knowledge of atmospheric CO2 levels for the simulated temperature and salinity spread in large parts of the South Atlantic. Overall temperature uncertainties reach up to 15 °C and increase significantly with water depth. The ensemble approach reveals temporal changes in the relative importance of geographic and radiative forcings for the simulated oceanographic conditions and, importantly, nonlinear interactions between them. The progressive northward opening of the highly restricted Angola Basin increased the sensitivity of local overturning and upper-ocean stratification to atmospheric CO2 concentrations due to large-scale changes in the hydrological cycle, while the chosen proto-Drake Passage depth is critical for the ocean dynamics and CO2 response in the southern South Atlantic. Finally, the simulated processes are integrated into a recent carbon burial framework to document the principal control of the regional gateway evolution on the progressive shift from the prevailing saline and oxygen-depleted subtropical water masses to the dominance of ventilated high-latitude deep waters.

An Assessment of Changes in the Thermal Environment during the COVID-19 Lockdown: Case Studies from the Greenland and Norwegian Seas

The COVID-19 lockdown had a significant impact on human activities, reducing anthropogenic heat and CO2 emissions. To effectively assess the impact of the lockdown on the thermal environment, we used the sliding paired t-test, which we improved from the traditional sliding t-test by introducing the paired t-test for sliding statistical tests, to test the abrupt change in the thermal environment. Furthermore, an additive decomposition model and wavelet analysis method were used to analyze the characteristics of trend and irregular change, coherence, and phase difference of the time series data with respect to the thermal environment. We chose the Greenland Sea and the Norwegian Sea, regions highly sensitive to changes in climate and ocean circulation, as case studies and used remote sensing data of the sea surface temperature (SST) and the atmospheric CO2 concentration data obtained from the Goddard Earth Sciences Data and Information Services Center from January 2015 to December 2021 for the analysis. The results show that although the annual spatial mean SST in 2020 is lower than the mean of all 7 years in most areas of the two seas, there is no evidence of a significant mutation in the decrease in the SST during the lockdown in 2020 compared with the temperatures before, according to the sliding paired t-test. The analysis of the irregular components of the monthly mean SST decomposed by an additive decomposition model also does not show the anomalously low SST during the lockdown in 2020. In addition, the lockdown had almost no impact on the increasing trend of CO2 concentration. The wavelet analysis also shows that there is no obvious anomaly in coherence or phase difference between the periodic variation of the SST and the CO2 concentrations in 2020 compared with other years. These results suggest that the direct effect of the COVID-19 lockdown on the thermal environment of the study area could be negligible.

Eccentricity Paces Late Pleistocene Glaciations

AbstractLate Pleistocene glacial terminations are caused by rising atmospheric CO2 occurring in response to atmospheric and ocean circulation changes induced by increased discharge from Northern Hemisphere ice sheets. While climate records place glacial terminations coincident with decreasing orbital precession, it remains unclear why a specific precession minimum causes a termination. We compare the orbital and ice volume configuration at each precession minima over the last million years to demonstrate that eccentricity, through its control on precession amplitude, period and coherence with obliquity, along with ice sheet size, determine whether a given precession minimum will cause a termination. We also demonstrate how eccentricity controls obliquity maxima and precession minima coherence, varying the duration of glaciations. Glaciations lasting ∼100 thousand years are controlled by Earth's eccentricity cycle of the same period, while the shortest (20–40 ka) and longest (155 ka) occupy the maxima and minimums of the 400 thousand year eccentricity cycle.

Primary Drivers of the Recent Changes Observed in Global Ocean Circulation Patterns in Mozambique

Purpose: The aim of the study was to examine the primary drivers of the recent changes observed in global ocean circulation patterns in Mozambique Methodology: This study adopted a desk methodology. A desk study research design is commonly known as secondary data collection. This is basically collecting data from existing resources preferably because of its low cost advantage as compared to a field research. Our current study looked into already published studies and reports as the data was easily accessed through online journals and libraries. Findings: The study revealed that primary drivers of the recent changes observed in global ocean circulation patterns are involving a combination of natural variability and anthropogenic influences. Key factors identified include global warming, alterations in wind patterns, sea surface temperature anomalies, and increased freshwater input from melting ice caps and glaciers. These elements interact in intricate ways to influence the dynamics of ocean currents, which are critical for regulating Earth's climate and supporting marine ecosystems. Unique Contribution to Theory, Practice and Policy: Bronfenbrenner's Ecological Systems Theory & Abraham Maslow's Hierarchy of Needs may be used to anchor future studies on primary drivers of the recent changes observed in global ocean circulation patterns in Mozambique. Increase investment in global ocean monitoring systems, such as Argo floats, satellite missions, and autonomous underwater vehicles, to collect high-resolution data on temperature, salinity, and currents. This data is crucial for validating models and understanding the current state and trends of ocean circulation. Strengthen international cooperation to address global warming, a primary driver of changes in ocean circulation. This includes adhering to and expanding commitments under the Paris Agreement to reduce greenhouse gas emissions.

Causes and multiyear predictability of the rapid acceleration of U.S. Southeast Sea level rise after 2010

The rate of sea level rise (SLR) along the Southeast Coast of the U.S. increased significantly after 2010. While anthropogenic radiative forcing causes an acceleration of global mean SLR, regional changes in the rate of SLR are strongly influenced by internal variability. Here we use observations and climate models to show that the rapid increase in the rate of SLR along the U.S. Southeast Coast after 2010 is due in part to multidecadal buoyancy-driven Atlantic meridional overturning circulation (AMOC) variations, along with heat transport convergence from wind-driven ocean circulation changes. We show that an initialized decadal prediction system can provide skillful regional SLR predictions induced by AMOC variations 5 years in advance, while wind-driven sea level variations are predictable 2 years in advance. Our results suggest that the rate of coastal SLR and its associated flooding risk along the U.S. southeastern seaboard are potentially predictable on multiyear timescales.

Leeuwin current dynamics in the SE Indian Ocean and implications for regional surface hydrography since the latest Miocene: Results from ODP Site 763A

This study examines the relationship between the ocean circulation dynamics in the eastern Indian Ocean and climate variability in Northwest (NW) Australia over the last 6.1 million years. The vertical water column of the southeastern Indian Ocean was analyzed using planktic foraminifera assemblages, and oxygen and carbon isotope ratios of surface and subsurface dwelling foraminifera from Ocean Drilling Program (ODP) Site 763A. The X-ray fluorescence data of bulk sediment from ODP Site 763A was also used to understand variability in relative humidity and aridity in NW Australia associated with surface ocean warming. The findings suggest a pronounced change in surface oceanography and climate in the region from 5.4 to 5 Ma, indicating a weak Leeuwin Current (LC) and reduction in mixed layer thickness. Subsequently, between 5.0 and 3.6 Ma, a strongly stratified water column with a thick mixed layer and deep thermocline emerged owing to the intensified Indonesian Throughflow (ITF). This led to a strengthened LC and diminished surface productivity, corresponding to a wet phase in NW Australia. The closure of the Indonesian Seaway during the Middle Pliocene (3.6–3.2 Ma) restricted the flow of West Pacific Warm Water into the Indian Ocean, weakened the LC and intensified the dry phase in NW Australia. Significant shifts in planktic foraminifera abundances are observed in the early Pleistocene (2.6–1.2 Ma) perhaps linked to major changes in surface ocean circulation. During the Mid-Pleistocene Transition (1.2–0.6 Ma), the West Australian Current exerted a significant influence over Site 763A.

Bacterial magnetofossil evidence for enhanced Pacific Ocean respired carbon storage during buildup of Antarctic glaciation

Global cooling with the onset of Antarctic glaciation ca. 34 Ma across the Eocene–Oligocene transition (EOT) terminated the early Cenozoic greenhouse climate state and marked the beginning of icehouse conditions. Although a pCO2 decline is considered to have been a major cause of this climate shift, the associated carbon-sequestration mechanism remains unclear. Here, we assessed ocean production and circulation changes across the EOT using numerical simulations combined with a novel proxy, namely, bacterial magnetofossils, the abundance and morphology of which are sensitive to sedimentary organic matter accumulation and oxygenation. We detected production and oxygenation declines in the equatorial Pacific Ocean coeval with increased biological production in the Southern Ocean after the EOT. Corroborated by simulation results and evidence from the Subantarctic region, we interpret this counterintuitive combination as a result of enhanced bottom-water formation and biological pump efficiency in the Southern Ocean due to Antarctic glacial buildup across the EOT. These results provide key evidence for deep Pacific Ocean deoxygenation and increased respired carbon concentrations, which amplified CO2 decline across the EOT.

Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene

Recent geologic and modeled evidence suggests that the grounding line of the Siple Coast of the West Antarctic Ice Sheet (WAIS) retreated hundreds of kilometers beyond its present position in the middle to late Holocene and readvanced within the past 1.7 ka. This grounding line reversal has been attributed to both changing rates of isostatic rebound and regional climate change. Here, we test these two hypotheses using a proxy-informed ensemble of ice sheet model simulations with varying ocean thermal forcing, global glacioisostatic adjustment (GIA) model simulations, and coupled ice sheet-GIA simulations that consider the interactions between these processes. Our results indicate that a warm to cold ocean cavity regime shift is the most likely cause of this grounding line reversal, but that GIA influences the rate of ice sheet response to oceanic changes. This implies that the grounding line here is sensitive to future changes in sub-ice shelf ocean circulation.