Description: The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) – induced by the closure of the gateway – was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling. Highlights: ► We study the effect of the Panama seaway on Pacific equatorial thermocline depth. ► Results from twelve model experiments are examined. ► Eastward net throughflow leads to a reduction in Atlantic overturning. ► We find a relationship between Panama closure and Pacific thermocline depth.

Description: Numerical models are important tools for understanding the processes and feedbacks in the Earth system, including those involving changes in atmospheric CO2 (CO2,atm) concentrations. Here, we compile 55 published model studies (consisting of 778 individual simulations) that assess the impact of six forcing mechanisms on millennial-scale CO2,atm variations: changes in freshwater supply to the North Atlantic and Southern Ocean, the strength and position of the southern-hemisphere westerlies, Antarctic sea ice extent, and aeolian dust fluxes. We generally find agreement on the direction of simulated CO2,atm change across simulations, but the amplitude of change is inconsistent, primarily due to the different complexities of the model representation of Earth system processes. When freshwater is added to the North Atlantic, a reduced Atlantic Meridional Overturning Circulation (AMOC) is generally accompanied by an increase in Southern Ocean- and Pacific overturning, reduced Antarctic sea ice extent, spatially varying export production, and changes in carbon storage in the Atlantic (rising), in other ocean basins (generally decreasing) and on land (more varied). Positive or negative CO2,atm changes are simulated during AMOC minima due to a spatially and temporally varying dominance of individual terrestrial and oceanic drivers (and compensating effects between them) across the different models. In contrast, AMOC recoveries are often accompanied by rising CO2,atm levels, which are mostly driven by ocean carbon release (albeit from different regions). The magnitude of simulated CO2,atm rise broadly scales with the duration of the AMOC perturbation (i.e., the stadial length). When freshwater is added to the Southern Ocean, reduced deep-ocean ventilation drives a CO2,atm drop via reduced carbon release from the Southern Ocean. Although the impacts of shifted southern-hemisphere westerlies are inconsistent across model simulations, their intensification raises CO2,atm via enhanced Southern Ocean Ekman pumping. Increased supply of aeolian dust to the ocean, and thus iron fertilisation of marine productivity, consistently lowers modelled CO2,atm concentrations via more efficient nutrient utilisation. The magnitude of CO2,atm change in response to dust flux variations, however, largely depends on the complexity of models' marine ecosystem and iron cycle. This especially applies to simulations forced by Antarctic sea ice changes, in which the direction of simulated CO2,atm change varies greatly across model hierarchies. Our compilation highlights that no single (forcing) mechanism can explain observed past millennial-scale CO2,atm variability, and identifies important future needs in coupled carbon cycle-climate modelling to better understand the mechanisms governing CO2,atm changes in the past.

Description: Deepwater circulation significantly changed during the last deglaciation from a shallow to a deep-reaching overturning cell. This change went along with a drawdown of isotopically light waters into the abyss and a deep ocean warming that changed deep ocean stratification from a salinity-to a temperature-controlled mode. Yet, the exact mechanisms causing these changes are still unknown. Furthermore, the long-standing idea of a complete shutdown of North Atlantic deepwater formation during Heinrich Stadial 1 (HS1) (17.5–14.6 kyr BP) remains prevalent. Here, we present a new compilation of benthic δ13C and δ18O data from the North Atlantic at high temporal resolution with consistent age models, established as part of the international PAGES working group OC3, to investigate deepwater properties in the North Atlantic. The extensive compilation, which includes 105 sediment cores, reveals different water masses during HS1. A water mass with heavy δ13C and δ18O signature occupies the Iceland Basin, whereas between 20 and 50°N, a distinct tongue of 18O depleted, 13C enriched water reaches down to 4000 m water depths. The heavy δ13C signature indicates active deepwater formation in the North Atlantic during HS1. Differences in its δ18O signature indicate either different sources or an alteration of the deepwater on its southward pathway. Based on these results, we discuss concepts of deepwater formation in the North Atlantic that help to explain the deglacial change from a salinity-driven to a temperature-driven circulation mode. Highlights • Spatial analyses of benthic δ13C and δ18O data from OC3 Atlantic compilation for HS1. • Heavy δ13C, light δ18O waters migrated into deep western North Atlantic basin during HS1. • Active deepwater formation between 30 and 60°N in the North Atlantic during HS1. • New concepts for transport of isotopically light δ18O into deep ocean. • Major contribution of North Atlantic waters to deglacial deep ocean stratification changes.

Description: Numerical models are important tools for understanding the processes and feedbacks in the Earth system, including those involving changes in atmospheric CO2 (CO2,atm) concentrations. Here, we compile 55 published model studies (consisting of 778 individual simulations) that assess the impact of six forcing mechanisms on millennial-scale CO2,atm variations: changes in freshwater supply to the North Atlantic and Southern Ocean, the strength and position of the southern-hemisphere westerlies, Antarctic sea ice extent, and aeolian dust fluxes. We generally find agreement on the direction of simulated CO2,atm change across simulations, but the amplitude of change is inconsistent, primarily due to the different complexities of the model representation of Earth system processes. When freshwater is added to the North Atlantic, a reduced Atlantic Meridional Overturning Circulation (AMOC) is generally accompanied by an increase in Southern Ocean- and Pacific overturning, reduced Antarctic sea ice extent, spatially varying export production, and changes in carbon storage in the Atlantic (rising), in other ocean basins (generally decreasing) and on land (more varied). Positive or negative CO2,atm changes are simulated during AMOC minima due to a spatially and temporally varying dominance of individual terrestrial and oceanic drivers (and compensating effects between them) across the different models. In contrast, AMOC recoveries are often accompanied by rising CO2,atm levels, which are mostly driven by ocean carbon release (albeit from different regions). The magnitude of simulated CO2,atm rise broadly scales with the duration of the AMOC perturbation (i.e., the stadial length). When freshwater is added to the Southern Ocean, reduced deep-ocean ventilation drives a CO2,atm drop via reduced carbon release from the Southern Ocean. Although the impacts of shifted southern-hemisphere westerlies are inconsistent across model simulations, their intensification raises CO2,atm via enhanced Southern Ocean Ekman pumping. Increased supply of aeolian dust to the ocean, and thus iron fertilisation of marine productivity, consistently lowers modelled CO2,atm concentrations via more efficient nutrient utilisation. The magni- tude of CO2,atm change in response to dust flux variations, however, largely depends on the complexity of models' marine ecosystem and iron cycle. This especially applies to simulations forced by Antarctic sea ice changes, in which the direction of simulated CO2,atm change varies greatly across model hierarchies. Our compilation highlights that no single (forcing) mechanism can explain observed past millennial-scale CO2,atm variability, and identifies important future needs in coupled carbon cycle-climate modelling to better understand the mechanisms governing CO2,atm changes in the past.

Description: We present the first version of the Ocean Circulation and Carbon Cycling (OC3) working group database, of oxygen and carbon stable isotope ratios from benthic foraminifera in deep ocean sediment cores from the Last Glacial Maximum (LGM, 23-19 ky) to the Holocene (〈10 ky) with a particular focus on the early last deglaciation (19-15 ky BP). It includes 287 globally distributed coring sites, with metadata, isotopic and chronostratigraphic information, and age models. A quality check was performed for all data and age models, and sites with at least millennial resolution were preferred. Deep water mass structure as well as differences between the early deglaciation and LGM are captured by the data, even though its coverage is still sparse in many regions. We find high correlations among time series calculated with different age models at sites that allow such analysis. The database provides a useful dynamical approach to map physical and biogeochemical changes of the ocean throughout the last deglaciation.