Search

We combine consistently dated benthic carbon isotopic records distributed over the entire Atlantic Ocean with numerical simulations performed by a glacial configuration of the Norwegian Earth System Model with active ocean biogeochemistry in order to interpret the observed Cibicides δ13C changes at the stadial–interstadial transition corresponding to the end of Heinrich Stadial 4 (HS4) in terms of ocean circulation and remineralization changes. We show that the marked increase in Cibicides δ13C observed at the end of HS4 between ∼2000 and 4200 m in the Atlantic can be explained by changes in nutrient concentrations as simulated by the model in response to the halting of freshwater input in the high-latitude glacial North Atlantic. Our model results show that this Cibicides δ13C signal is associated with changes in the ratio of southern-sourced (SSW) versus northern-sourced (NSW) water masses at the core sites, whereby SSW is replaced by NSW as a consequence of the resumption of deep-water formation in the northern North Atlantic and Nordic Seas after the freshwater input is halted. Our results further suggest that the contribution of ocean circulation changes to this signal increases from ∼40 % at 2000 m to ∼80 % at 4000 m. Below ∼4200 m, the model shows little ocean circulation change but an increase in remineralization across the transition marking the end of HS4. The simulated lower remineralization during stadials compared to during interstadials is particularly pronounced in deep subantarctic sites, in agreement with the decrease in the export production of carbon to the deep Southern Ocean during stadials found in previous studies.

The characteristics of the CISE-LOCEAN seawater isotope dataset (δ18O, δ2H, referred to as δD) are presented (https://doi.org/10.17882/71186; Waterisotopes-CISE-LOCEAN, 2021). This dataset covers the time period from 1998 to 2021 and currently includes close to 8000 data entries, all with δ18O, three-quarters of them also with δD, associated with a date stamp, space stamp, and usually a salinity measurement. Until 2010, samples were analyzed by isotopic ratio mass spectrometry and since then mostly by cavity ring-down spectroscopy (CRDS). Instrumental uncertainty in this dataset is usually as low as 0.03 ‰ for δ18O and 0.15 ‰ for δD. An additional uncertainty is related to the isotopic composition of the in-house standards that are used to convert data to the Vienna Standard Mean Ocean Water (VSMOW) scale. Different comparisons suggest that since 2010 the latter have remained within at most 0.03 ‰ for δ18O and 0.20 ‰ for δD. Therefore, combining the two uncertainties suggests a standard deviation of at most 0.05 ‰ for δ18O and 0.25 ‰ for δD. For some samples, we find that there has been evaporation during collection and storage, requiring adjustment of the isotopic data produced by CRDS, based on d-excess (δD − 8×δ18O). This adjustment adds an uncertainty in the respective data of roughly 0.05 ‰ for δ18O and 0.10 ‰ for δD. This issue of conservation of samples is certainly a strong source of quality loss for parts of the database, and “small” effects may have remained undetected. The internal consistency of the database can be tested for subsets of the dataset when time series can be obtained (such as in the southern Indian Ocean or North Atlantic subpolar gyre). These comparisons suggest that the overall uncertainty of the spatially (for a cruise) or temporally (over a year) averaged data is less than 0.03 ‰ for δ18O and 0.15 ‰ for δD. However, 18 comparisons with duplicate seawater data analyzed in other laboratories or with other datasets in the intermediate and deep ocean suggest a larger scatter. When averaging the 18 comparisons done for δ18O, we find a difference of 0.082 ‰ with a standard error of 0.016 ‰. Such an average difference is expected due to the adjustments applied at LOCEAN to saline water data produced either by CRDS or isotope ratio mass spectrometry (IRMS), but the scatter found suggests that care is needed when merging datasets from different laboratories. Examples of time series in the surface North Atlantic subpolar gyre illustrate the temporal changes in water isotope composition that can be detected with a carefully validated dataset.

Changes in water mass distribution are considered to be a significant contributor to the atmospheric CO2 concentration drop to around 186 ppm recorded during the Last Glacial Maximum (LGM). Yet simulating a glacial Atlantic Meridional Overturning Circulation (AMOC) in agreement with paleotracer data remains a challenge, with most models from previous Paleoclimate Modelling Intercomparison Project (PMIP) phases showing a tendency to simulate a strong and deep North Atlantic Deep Water (NADW) instead of the shoaling inferred from proxy records of water mass distribution. Conversely, the simulated Antarctic Bottom Water (AABW) is often reduced compared to its pre-industrial volume, and the Atlantic Ocean stratification is underestimated with respect to paleoproxy data. Inadequate representation of surface conditions, driving deep convection around Antarctica, may explain inaccurately simulated bottom water properties in the Southern Ocean. We investigate here the impact of a range of surface conditions in the Southern Ocean in the iLOVECLIM model using nine simulations obtained with different LGM boundary conditions associated with the ice sheet reconstruction (e.g., changes of elevation, bathymetry, and land–sea mask) and/or modeling choices related to sea-ice export, formation of salty brines, and freshwater input. Based on model–data comparison of sea-surface temperatures and sea ice, we find that only simulations with a cold Southern Ocean and a quite extensive sea-ice cover show an improved agreement with proxy records of sea ice, despite systematic model biases in the seasonal and regional patterns. We then show that the only simulation which does not display a much deeper NADW is obtained by parameterizing the sinking of brines along Antarctica, a modeling choice reducing the open-ocean convection in the Southern Ocean. These results highlight the importance of the representation of convection processes, which have a large impact on the water mass properties, while the choice of boundary conditions appears secondary for the model resolution and variables considered in this study.

Understanding the ocean circulation changes associated with abrupt climate events is key to better assessing climate variability and understanding its different natural modes. Sedimentary Pa∕Th, benthic δ13C and Δ14C are common proxies used to reconstruct past circulation flow rate and ventilation. To overcome the limitations of each proxy taken separately, a better approach is to produce multiproxy measurements on a single sediment core. Yet, different proxies can provide conflicting information about past ocean circulation. Thus, modelling them in a consistent physical framework has become necessary to assess the geographical pattern and the timing and sequence of the multiproxy response to abrupt circulation changes. We have implemented a representation of the 231Pa and 230Th tracers into the model of intermediate complexity iLOVECLIM, which already included δ13C and Δ14C. We have further evaluated the response of these three ocean circulation proxies to a classical abrupt circulation reduction obtained by freshwater addition in the Nordic Seas under preindustrial boundary conditions. The proxy response is shown to cluster in modes that resemble the modern Atlantic water masses. The clearest and most coherent response is obtained in the deep (> 2000 m) northwest Atlantic, where δ13C and Δ14C significantly decrease, while Pa∕Th increases. This is consistent with observational data across millennial-scale events of the last glacial. Interestingly, while in marine records, except in rare instances, the phase relationship between these proxies remains unclear due to large dating uncertainties, in the model the bottom water carbon isotope (δ13C and Δ14C) response lags behind the sedimentary Pa∕Th response by a few hundred years.

A complete understanding of past El Niño–Southern Oscillation (ENSO) fluctuations is important for the future predictions of regional climate using climate models. One approach to reconstructing past ENSO dynamics uses planktonic foraminifera as recorders of past climate to assess past spatio-temporal changes in upper ocean conditions. In this paper, we utilise a model of planktonic foraminifera populations, Foraminifera as Modelled Entities (FAME), to forward model the potential monthly average δ18Oc and temperature signal proxy values for Globigerinoides ruber, Globigerinoides sacculifer, and Neogloboquadrina dutertrei from input variables covering the period of the instrumental record. We test whether the modelled foraminifera population δ18Oc and Tc associated with El Niño events statistically differ from the values associated with other climate states. Provided the assumptions of the model are correct, our results indicate that the values of El Niño events can be differentiated from other climate states using these species. Our model computes the proxy values of foraminifera in the water, suggesting that, in theory, water locations for a large portion of the tropical Pacific should be suitable for differentiating El Niño events from other climate states. However, in practice it may not be possible to differentiate climate states in the sediment record. Specifically, comparison of our model results with the sedimentological features of the Pacific Ocean shows that a large portion of the hydrographically/ecologically suitable water regions coincide with low sediment accumulation rate at the sea floor and/or of sea floor that lie below threshold water depths for calcite preservation.

The oxygen-18 to oxygen-16 ratio recorded in fossil planktonic foraminifer shells has been used for over 50 years in many geoscience applications. However, different planktonic foraminifer species generally yield distinct signals, as a consequence of their specific living habitats in the water column and along the year. This complexity is usually not taken into account in model–data integration studies. To overcome this shortcoming, we developed the Foraminifers As Modeled Entities (FAME) module. The module predicts the presence or absence of commonly used planktonic foraminifers and their oxygen-18 values. It is only forced by hydrographic data and uses a very limited number of parameters, almost all derived from culture experiments. FAME performance is evaluated using the Multiproxy Approach for the Reconstruction of the Glacial Ocean surface (MARGO) Late Holocene planktonic foraminifer calcite oxygen-18 and abundance datasets. The application of FAME to a simple cooling scenario demonstrates its utility to predict changes in planktonic foraminifer oxygen-18 to oxygen-16 ratio in response to changing climatic conditions.

Thanks to its optimal location on the northern Brazilian margin, core MD09-3257 records both ocean circulation and atmospheric changes. The latter occur locally in the form of increased rainfall on the adjacent continent during the cold intervals recorded in Greenland ice and northern North Atlantic sediment cores (i.e., Greenland stadials). These rainfall events are recorded in MD09-3257 as peaks in ln(Ti ∕ Ca). New sedimentary Pa ∕ Th data indicate that mid-depth western equatorial water mass transport decreased during all of the Greenland stadials of the last 40 kyr. Using cross-wavelet transforms and spectrogram analysis, we assess the relative phase between the MD09-3257 sedimentary Pa ∕ Th and ln(Ti ∕ Ca) signals. We show that decreased water mass transport between a depth of ∼ 1300 and 2300 m in the western equatorial Atlantic preceded increased rainfall over the adjacent continent by 120 to 400 yr at Dansgaard–Oeschger (D–O) frequencies, and by 280 to 980 yr at Heinrich-like frequencies.We suggest that the large lead of ocean circulation changes with respect to changes in tropical South American precipitation at Heinrich-like frequencies is related to the effect of a positive feedback involving iceberg discharges in the North Atlantic. In contrast, the absence of widespread ice rafted detrital layers in North Atlantic cores during D–O stadials supports the hypothesis that a feedback such as this was not triggered in the case of D–O stadials, with circulation slowdowns and subsequent changes remaining more limited during D–O stadials than Heinrich stadials.

Abstract Excess 231Pa and 230Th (231Paxs and 230Thxs) can be used to reconstruct past oceanic sedimentation (230Th‐normalized flux) and circulation changes ((231Pa/230Th)xs,0, hereafter Pa/Th). These quantities are determined by computing the detrital and authigenic contributions from bulk sediment measurement. The method relies on the use of a chosen constant value of the detrital (238U/232Th) activity ratio (hereafter (U/Th)det). In this study, we have extracted the detrital fraction of the sediments from North Atlantic deep‐sea core SU90‐08 (43°03′1N, 30°02′5W, 3,080m) and determined its (U/Th)det value over the last 40 ky. We find that (U/Th)det varied significantly through time with a minimum value of 0.4 during the Holocene and a maximum value of 0.7 during the Last Glacial Maximum (LGM). The sensitivity of sedimentary 230Th‐normalized flux and Pa/Th is tested for our study site and for other North Atlantic sites. We show that the sensitivity is highly dependent on the core location and its terrigenous material supply. The 230Th‐normalized flux and Pa/Th signals are very robust in cores with low detrital contributions, whereas they are very sensitive to (U/Th)det changes in cores with higher detrital contribution. In the latter case, changes in 230Th‐normalized flux and Pa/Th due to the choice of a constant (U/Th)det can largely exceed the uncertainty on the 230Th‐normalized flux and Pa/Th, inducing potential biases in the amplitude and temporal variability of reconstructed sedimentation and ocean circulation changes.

Past climate is an important benchmark to assess the ability of climate models to simulate key processes and feedbacks. Numerous proxy records exist for stable isotopes of water and/or carbon, which are also implemented inside the components of a growing number of Earth system model. Model–data comparisons can help to constrain the uncertainties associated with transfer functions. This motivates the need of producing a comprehensive compilation of different proxy sources. We have put together a global database of proxy records of oxygen (δ18O), hydrogen (δD) and carbon (δ13C) stable isotopes from different archives: ocean and lake sediments, corals, ice cores, speleothems and tree-ring cellulose. Source records were obtained from the georeferenced open access PANGAEA and NOAA libraries, complemented by additional data obtained from a literature survey. About 3000 source records were screened for chronological information and temporal resolution of proxy records. Altogether, this database consists of hundreds of dated δ18O, δ13C and δD records in a standardized simple text format, complemented with a metadata Excel catalog. A quality control flag was implemented to describe age markers and inform on chronological uncertainty. This compilation effort highlights the need to homogenize and structure the format of datasets and chronological information as well as enhance the distribution of published datasets that are currently highly fragmented and scattered. We also provide an online portal based on the records included in this database with an intuitive and interactive platform (http://climateproxiesfinder.ipsl.fr/), allowing one to easily select, visualize and download subsets of the homogeneously formatted records that constitute this database, following a choice of search criteria, and to upload new datasets. In the last part, we illustrate the type of application allowed by our database by comparing several key periods highly investigated by the paleoclimate community. For coherency with the Paleoclimate Modelling Intercomparison Project (PMIP), we focus on records spanning the past 200 years, the mid-Holocene (MH, 5.5–6.5 ka; calendar kiloyears before 1950), the Last Glacial Maximum (LGM, 19–23 ka), and those spanning the last interglacial period (LIG, 115–130 ka). Basic statistics have been applied to characterize anomalies between these different periods. Most changes from the MH to present day and from LIG to MH appear statistically insignificant. Significant global differences are reported from LGM to MH with regional discrepancies in signals from different archives and complex patterns.

The Atlantic Ocean overturning circulation and its impact on respired carbon storage during glacial periods remains a mystery. Here, the authors combine radiocarbon and stable carbon isotope data from the Brazil and Iberian Margins to investigate ocean circulation and carbon storage at the Last Glacial Maximum.

We present and examine a multi-sensor global compilation of mid-Holocene (MH) sea surface temperatures (SST), based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst census counts. We assess the uncertainties originating from using different methodologies and evaluate the potential of MH SST reconstructions as a benchmark for climate-model simulations. The comparison between different analytical approaches (time frame, baseline climate) shows the choice of time window for the MH has a negligible effect on the reconstructed SST pattern, but the choice of baseline climate affects both the magnitude and spatial pattern of the reconstructed SSTs. Comparison of the SST reconstructions made using different sensors shows significant discrepancies at a regional scale, with uncertainties often exceeding the reconstructed SST anomaly. Apparent patterns in SST may largely be a reflection of the use of different sensors in different regions. Overall, the uncertainties associated with the SST reconstructions are generally larger than the MH anomalies. Thus, the SST data currently available cannot serve as a target for benchmarking model simulations. Further evaluations of potential subsurface and/or seasonal artifacts that may contribute to obscure the MH SST reconstructions are urgently needed to provide reliable benchmarks for model evaluations.

We use the fully coupled atmosphere–ocean three-dimensional model of intermediate complexity iLOVECLIM to simulate the climate and oxygen stable isotopic signal during the Last Glacial Maximum (LGM, 21 000 years). By using a model that is able to explicitly simulate the sensor (δ18O), results can be directly compared with data from climatic archives in the different realms. Our results indicate that iLOVECLIM reproduces well the main feature of the LGM climate in the atmospheric and oceanic components. The annual mean δ18O in precipitation shows more depleted values in the northern and southern high latitudes during the LGM. The model reproduces very well the spatial gradient observed in ice core records over the Greenland ice sheet. We observe a general pattern toward more enriched values for continental calcite δ18O in the model at the LGM, in agreement with speleothem data. This can be explained by both a general atmospheric cooling in the tropical and subtropical regions and a reduction in precipitation as confirmed by reconstruction derived from pollens and plant macrofossils. Data–model comparison for sea surface temperature indicates that iLOVECLIM is capable to satisfyingly simulate the change in oceanic surface conditions between the LGM and present. Our data–model comparison for calcite δ18O allows investigating the large discrepancies with respect to glacial temperatures recorded by different microfossil proxies in the North Atlantic region. The results argue for a strong mean annual cooling in the area south of Iceland and Greenland between the LGM and present (> 6 °C), supporting the foraminifera transfer function reconstruction but in disagreement with alkenones and dinocyst reconstructions. The data–model comparison also reveals that large positive calcite δ18O anomaly in the Southern Ocean may be explained by an important cooling, although the driver of this pattern is unclear. We deduce a large positive δ18Osw anomaly for the north Indian Ocean that contrasts with a large negative δ18Osw anomaly in the China Sea between the LGM and the present. This pattern may be linked to changes in the hydrological cycle over these regions. Our simulation of the deep ocean suggests that changes in δ18Osw between the LGM and the present are not spatially homogeneous. This is supported by reconstructions derived from pore fluids in deep-sea sediments. The model underestimates the deep ocean cooling thus biasing the comparison with benthic calcite δ18O data. Nonetheless, our data–model comparison supports a heterogeneous cooling of a few degrees (2–4 °C) in the LGM Ocean.

Although the Last Interglacial (LIG) is often considered as a possible analogue for future climate in high latitudes, its precise climate evolution and associated causes remain uncertain. Here we compile high-resolution marine sediment records from the North Atlantic, Labrador Sea, Norwegian Sea and the Southern Ocean. We document a delay in the establishment of peak interglacial conditions in the North Atlantic, Labrador and Norwegian Seas as compared to the Southern Ocean. In particular, we observe a persistent iceberg melting at high northern latitudes at the beginning of the LIG. It is associated with (1) colder and fresher surface-water conditions in the North Atlantic, Labrador and Norwegian Seas, and (2) a weaker ventilation of North Atlantic deep waters during the early LIG (129–125 ka) compared to the late LIG. Results from an ocean-atmosphere coupled model with insolation as a sole forcing for three key periods of the LIG show warmer North Atlantic surface waters and stronger Atlantic overturning during the early LIG (126 ka) than the late LIG (122 ka). Hence, insolation variations alone do not explain the delay in peak interglacial conditions observed at high northern latitudes. Additionally, we consider an idealized meltwater scenario at 126 ka where the freshwater input is interactively computed in response to the high boreal summer insolation. The model simulates colder, fresher North Atlantic surface waters and weaker Atlantic overturning during the early LIG (126 ka) compared to the late LIG (122 ka). This result suggests that both insolation and ice sheet melting have to be considered to reproduce the climatic pattern that we identify during the early LIG. Our model-data comparison also reveals a number of limitations and reinforces the need for further detailed investigations using coupled climate-ice sheet models and transient simulations.

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

The EPICA (European Project for Ice Coring in Antarctica) Dome C drilling in East Antarctica has now been completed to a depth of 3260 m, at only a few meters above bedrock. Here we present the new EDC3 chronology, which is based on the use of 1) a snow accumulation and mechanical flow model, and 2) a set of independent age markers along the core. These are obtained by pattern matching of recorded parameters to either absolutely dated paleoclimatic records, or to insolation variations. We show that this new time scale is in excellent agreement with the Dome Fuji and Vostok ice core time scales back to 100 kyr within 1 kyr. Discrepancies larger than 3 kyr arise during MIS 5.4, 5.5 and 6, which points to anomalies in either snow accumulation or mechanical flow during these time periods. We estimate that EDC3 gives accurate event durations within 20% (2σ) back to MIS11 and accurate absolute ages with a maximum uncertainty of 6 kyr back to 800 kyr.

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., (© 2023. The Author(s).)

Rapid changes in ocean circulation and climate have been observed in marine-sediment and ice cores over the last glacial period and deglaciation, highlighting the non-linear character of the climate system and underlining the possibility of rapid climate shifts in response to anthropogenic greenhouse gas forcing. To date, these rapid changes in climate and ocean circulation are still not fully explained. One obstacle hindering progress in our understanding of the interactions between past ocean circulation and climate changes is the difficulty of accurately dating marine cores. Here, we present a set of 92 marine sediment cores from the Atlantic Ocean for which we have established age-depth models that are consistent with the Greenland GICC05 ice core chronology, and computed the associated dating uncertainties, using a new deposition modeling technique. This is the first set of consistently dated marine sediment cores enabling paleoclimate scientists to evaluate leads/lags between circulation and climate changes over vast regions of the Atlantic Ocean. Moreover, this data set is of direct use in paleoclimate modeling studies.

Producing independent and accurate chronologies for marine sediments is a prerequisite to understand the sequence of millennial-scale events and reveal potential temporal offsets between marine and continental records, or between different marine records, possibly from different regions. The last 40 ky is a generally well-constrained period since radiocarbon ( 14 C) can be used as an absolute dating tool. However, in the northern North Atlantic, calendar ages cannot be directly derived from 14 C ages, due to temporal and spatial variations of surface reservoir ages. Alternatively, chronologies can be derived by aligning Greenland ice-core time series with marine surface records. Yet this approach suffers from the lack of clearly defined climatic events between 14.7 and 23.3 cal ky BP (hereafter ka), a crucial period encompassing Heinrich Stadial 1 and the onset of the last deglaciation. In this study, (i) we assess the benefits of 230 Th normalization to refine the sedimentation history between surface temperature alignment tie points and (ii) revisit the chronologies of three North Atlantic marine records. Our study supports the contention that the marked increase in the Greenland Ca 2+ record at 17.48 ka ± 0.21 ky (1 σ ) occurred within dating uncertainty of sea surface temperature cooling in the North Atlantic at the onset of Heinrich Stadial 1. This sharp feature might be useful for future chronostratigraphic alignments to remedy the lack of chronological constraint between 14.7 and 23.3 ka for North Atlantic marine records that are subject to large changes in 14 C surface reservoir age., (© 2019. The Authors.)

The influence of ocean circulation changes on atmospheric CO2 hinges primarily on the ability to alter the ocean interior's respired nutrient inventory. Here we investigate the Atlantic overturning circulation at the Last Glacial Maximum and its impact on respired carbon storage using radiocarbon and stable carbon isotope data from the Brazil and Iberian Margins. The data demonstrate the existence of a shallow well-ventilated northern-sourced cell overlying a poorly ventilated, predominantly southern-sourced cell at the Last Glacial Maximum. We also find that organic carbon remineralization rates in the deep Atlantic remained broadly similar to modern, but that ventilation ages in the southern-sourced overturning cell were significantly increased. Respired carbon storage in the deep Atlantic was therefore enhanced during the last glacial period, primarily due to an increase in the residence time of carbon in the deep ocean, rather than an increase in biological carbon export.

Millennial-scale climate changes during the last glacial period and deglaciation were accompanied by rapid changes in atmospheric CO2 that remain unexplained. While the role of the Southern Ocean as a 'control valve' on ocean-atmosphere CO2 exchange has been emphasized, the exact nature of this role, in particular the relative contributions of physical (for example, ocean dynamics and air-sea gas exchange) versus biological processes (for example, export productivity), remains poorly constrained. Here we combine reconstructions of bottom-water [O2], export production and (14)C ventilation ages in the sub-Antarctic Atlantic, and show that atmospheric CO2 pulses during the last glacial- and deglacial periods were consistently accompanied by decreases in the biological export of carbon and increases in deep-ocean ventilation via southern-sourced water masses. These findings demonstrate how the Southern Ocean's 'organic carbon pump' has exerted a tight control on atmospheric CO2, and thus global climate, specifically via a synergy of both physical and biological processes.

Explanations of the glacial-interglacial variations in atmospheric pCO2 invoke a significant role for the deep ocean in the storage of CO2. Deep-ocean density stratification has been proposed as a mechanism to promote the storage of CO2 in the deep ocean during glacial times. A wealth of proxy data supports the presence of a "chemical divide" between intermediate and deep water in the glacial Atlantic Ocean, which indirectly points to an increase in deep-ocean density stratification. However, direct observational evidence of changes in the primary controls of ocean density stratification, i.e., temperature and salinity, remain scarce. Here, we use Mg/Ca-derived seawater temperature and salinity estimates determined from temperature-corrected δ(18)O measurements on the benthic foraminifer Uvigerina spp. from deep and intermediate water-depth marine sediment cores to reconstruct the changes in density of sub-Antarctic South Atlantic water masses over the last deglaciation (i.e., 22-2 ka before present). We find that a major breakdown in the physical density stratification significantly lags the breakdown of the deep-intermediate chemical divide, as indicated by the chemical tracers of benthic foraminifer δ(13)C and foraminifer/coral (14)C. Our results indicate that chemical destratification likely resulted in the first rise in atmospheric pCO2, whereas the density destratification of the deep South Atlantic lags the second rise in atmospheric pCO2 during the late deglacial period. Our findings emphasize that the physical and chemical destratification of the ocean are not as tightly coupled as generally assumed.

Recent theories for glacial-interglacial climate transitions call on millennial climate perturbations that purged the deep sea of sequestered carbon dioxide via a "bipolar ventilation seesaw." However, the viability of this hypothesis has been contested, and robust evidence in its support is lacking. Here we present a record of North Atlantic deep-water radiocarbon ventilation, which we compare with similar data from the Southern Ocean. A striking coherence in ventilation changes is found, with extremely high ventilation ages prevailing across the deep Atlantic during the last glacial period. The data also reveal two reversals in the ventilation gradient between the deep North Atlantic and Southern Ocean during Heinrich Stadial 1 and the Younger Dryas. These coincided with periods of sustained atmospheric CO2 rise and appear to have been driven by enhanced ocean-atmosphere exchange, primarily in the Southern Ocean. These results confirm the operation of a bipolar ventilation seesaw during deglaciation and underline the contribution of abrupt regional climate anomalies to longer-term global climate transitions.

The deuterium excess of polar ice cores documents past changes in evaporation conditions and moisture origin. New data obtained from the European Project for Ice Coring in Antarctica Dome C East Antarctic ice core provide new insights on the sequence of events involved in Termination II, the transition between the penultimate glacial and interglacial periods. This termination is marked by a north-south seesaw behavior, with first a slow methane concentration rise associated with a strong Antarctic temperature warming and a slow deuterium excess rise. This first step is followed by an abrupt north Atlantic warming, an abrupt resumption of the East Asian summer monsoon, a sharp methane rise, and a CO(2) overshoot, which coincide within dating uncertainties with the end of Antarctic optimum. Here, we show that this second phase is marked by a very sharp Dome C centennial deuterium excess rise, revealing abrupt reorganization of atmospheric circulation in the southern Indian Ocean sector.

Past glacial-interglacial increases in the concentration of atmospheric carbon dioxide (CO2) are thought to arise from the rapid release of CO2 sequestered in the deep sea, primarily via the Southern Ocean. Here, we present radiocarbon evidence from the Atlantic sector of the Southern Ocean that strongly supports this hypothesis. We show that during the last glacial period, deep water circulating around Antarctica was more than two times older than today relative to the atmosphere. During deglaciation, the dissipation of this old and presumably CO2-enriched deep water played an important role in the pulsed rise of atmospheric CO2 through its variable influence on the upwelling branch of the Antarctic overturning circulation.

An 60Fe peak in a deep-sea FeMn crust has been interpreted as due to the signature left by the ejecta of a supernova explosion close to the solar system 2.8+/-0.4 Myr ago [Knie, Phys. Rev. Lett. 93, 171103 (2004)10.1103/PhysRevLett.93.171103]. In an attempt to confirm this interpretation with better time resolution and obtain a more direct flux estimate, we measured 60Fe concentrations along a dated marine sediment. We find no 60Fe peak at the expected level from 1.7 to 3.2 Myr ago. Possible causes for the discrepancy are discussed.

To determine the mechanisms governing the last deglaciation and the sequence of events that lead to deglaciation, it is important to obtain a temporal framework that applies to both continental and marine climate records. Radiocarbon dating has been widely used to derive calendar dates for marine sediments, but it rests on the assumption that the 'apparent age' of surface water (the age of surface water relative to the atmosphere) has remained constant over time. Here we present new evidence for variation in the apparent age of surface water (or reservoir age) in the North Atlantic ocean north of 40 degrees N over the past 20,000 years. In two cores we found apparent surface-water ages to be larger than those of today by 1,230 +/- 600 and 1,940 +/- 750 years at the end of the Heinrich 1 surge event (15,000 years BP) and by 820 +/- 430 to 1,010 +/- 340 years at the end of the Younger Dryas cold episode. During the warm Bølling-Allerød period, between these two periods of large reservoir ages, apparent surface-water ages were comparable to present values. Our results allow us to reconcile the chronologies from ice cores and the North Atlantic marine records over the entire deglaciation period. Moreover, the data imply that marine carbon dates from the North Atlantic north of 40 degrees N will need to be corrected for these highly variable effects.