Uvigerina spp. d18O and Mg/Ca measurements from sediment cores JR244-GC528 and MD07-3076

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 d18O 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 d13C and foraminifer/coral 14C. 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.

(Table S1) Age determination of sediment cores of the Gulf of Alaska

Oxygen isotope data from planktonic and benthic foraminifera, on a high-resolution age model (44 14C dates spanning 17,400 years), document deglacial environmental change on the southeast Alaska margin (59°33.32'N, 144°9.21'W, 682 m water depth). Surface freshening (i.e., d18O reduction of 0.8 per mil) began at 16,650 ± 170 cal years B.P. during an interval of ice proximal sedimentation, likely due to freshwater input from melting glaciers. A sharp transition to laminated hemipelagic sediments constrains retreat of regional outlet glaciers onto land circa 14,790 ± 380 cal years B.P. Abrupt warming and/or freshening of the surface ocean (i.e., additional d18O reduction of 0.9 per mil) coincides with the Bølling Interstade of northern Europe and Greenland. Cooling and/or higher salinities returned during the Allerød interval, coincident with the Antarctic Cold Reversal, and continue until 11,740 ± 200 cal years B.P., when onset of warming coincides with the end of the Younger Dryas. An abrupt 1 per mil reduction in benthic d18O at 14,250 ± 290 cal years B.P. likely reflects a decrease in bottom water salinity driven by deep mixing of glacial meltwater, a regional megaflood event, or brine formation associated with sea ice. Two laminated opal-rich intervals record discrete episodes of high productivity during the last deglaciation. These events, precisely dated here at 14,790 ± 380 to 12,990 ± 190 cal years B.P. and 11,160 ± 130 to 10,750 ± 220 cal years B.P., likely correlate to similar features observed elsewhere on the margins of the North Pacific and are coeval with episodes of rapid sea level rise. Remobilization of iron from newly inundated continental shelves may have helped to fuel these episodes of elevated primary productivity and sedimentary anoxia.

Geochemistry of sediment core PS2456-4

Sediment core PS2458 from the Laptev Sea continental margin (983-m water depth) stems from a position close to the paleoriver mouth of Lena and Yana rivers. It was dated by AMS-14C and analyzed in high resolution for oxygen isotopes of planktic foraminifers. Except the uppermost 100 cm and possibly the lowermost meter of the 8-m-long core, the sediments were deposited during the last deglaciation (14.5-8.0 cal-ka). According to 210Pb data, the uppermost 100 cm represents only the last 200 years. Planktic foraminifers are present throughout the dated deglacial interval, with the exception of a short time after ca. 13 cal-ka. Taking into account the global "ice volume effect" on the oxygen isotopic composition of the foraminifers, the isotopic record is considered to reflect salinity changes which were influenced by variable freshwater runoff and a growing marine influence during the postglacial transgression of the Laptev Sea shelf. The most conspicuous feature in the isotopic record is an outstanding peak dated to ca. 13 cal-ka. It is proposed that it represents a rapid outburst of large amounts of freshwater, possibly from an ice-dammed lake in the hinterland. Possible correlations to the onset of the cool Younger Dryas event in the northern hemisphere are discussed.

(Table 1) Age determination of sediment core JT96-09PC

he oxygen minimum zone (OMZ) off Vancouver Island was more oxygen depleted relative to modern conditions during the Allerød (~13.5 to 12.6 calendar kyr) and again from ~11 to 10 kyr. The timing of OMZ intensification is similar to that seen throughout the North Pacific, although the onset appears to have been delayed by ~1500 years off Vancouver Island. Radiocarbon dating of coeval benthic and planktonic foraminifera shows that between 16.0 and 12.6 kyr the age contrast between surface and intermediate waters (920 m depth) off Vancouver Island was similar to, or slightly less than, that today. There is no evidence of an increased age difference (i.e., decreased ventilation) during the deglaciation, particularly during the Allerød. However, sedimentary marine organic carbon concentration and mass accumulation rate increased substantially in the Allerød, suggesting that increased organic matter export was the principal cause of late Pleistocene OMZ intensification off Vancouver Island.

High-resolution FTIRS-inferred biogeochemical properties for MIS 11c in the sediment record of core ICDP5011-1 in Lake Elgygytgyn, NE Russia

Here we present a detailed multi-proxy record of the climate and environmental evolution at Lake El'gygytgyn, Far East Russian Arctic during the period 430-395 ka covering the marine isotope stage (MIS) 12/11 transition and the thermal maximum of super interglacial MIS 11c. The MIS 12/11 transition at Lake El'gygytgyn is characterized by initial warming followed by a cold reversal implying similarities to the last deglaciation. The thermal maximum of MIS 11c is characterized by full and remarkably stable interglacial conditions with mean temperatures of the warmest month (MTWM) ranging between ca. 10-15 °C; annual precipitation (PANN) ranging between ca. 300-600 mm; strong in-lake productivity coinciding with dark coniferous forests in the catchment; annual disintegration of the lake ice cover; and full mixis of the water column. Such conditions persisted, according to our age model, for ca. 27 ± 8 kyr between ca. 425-398 ka. The Lake El'gygytgyn record closely resembles the climate pattern recorded in Lake Baikal (SE Siberia) sediments and Antarctic ice cores, implying interhemispheric climate connectivity during MIS 11c.

(Appendix 1) Coral analysis and isostatic rebound effects from different Holes of IODP Expedition 310

Past sea-level records provide invaluable information about the response of ice sheets to climate forcing. Some such records suggest that the last deglaciation was punctuated by a dramatic period of sea-level rise, of about 20 metres, in less than 500 years. Controversy about the amplitude and timing of this meltwater pulse (MWP-1A) has, however, led to uncertainty about the source of the melt water and its temporal and causal relationships with the abrupt climate changes of the deglaciation. Here we show that MWP-1A started no earlier than 14,650 years ago and ended before 14,310 years ago, making it coeval with the Bølling warming. Our results, based on corals drilled offshore from Tahiti during Integrated Ocean Drilling Project Expedition 310, reveal that the increase in sea level at Tahiti was between 12 and 22 metres, with a most probable value between 14 and 18 metres, establishing a significant meltwater contribution from the Southern Hemisphere. This implies that the rate of eustatic sea-level rise exceeded 40 millimetres per year during MWP-1A.

(Table 2) Age determination of sediment core NIOP-C2_929

In this study we present a sea surface temperature (SST) record from the western Arabian Sea for the last 20,000 years. We produced centennial-scale d18O and Mg/Ca SST time series of core NIOP929 with focus on the glacial-interglacial transition. The western Arabian Sea is influenced by the seasonal NE and SW monsoon wind systems. Lowest SSTs occur during the SW monsoon season because of upwelling of cold water, and highest SSTs can be found in the low-productivity intermonsoon season. The Mg/Ca-based temperature record reflects the integrated SST of the SW and NE monsoon seasons. The results show a glacial-interglacial SST difference of ~2°C, which is corroborated by findings from other Arabian Sea cores. At 19 ka B.P. a yet undescribed warm event of several hundred years duration is found, which is also reflected in the d18O record. A second centennial-scale high SST/low d18O event is observed at 17 ka B.P. This event forms the onset of the stepwise yet persistent trend toward Holocene temperatures. Highest Mg/Ca-derived SSTs in the NIOP929 record occurred between 13 and 10 ka B.P. Interglacial SST is ~24°C, indicating influence of upwelling. The onset of Arabian Sea warming occurs when the North Atlantic is experiencing minimum temperatures. The rapid temperature variations at 19, 17, and 13 ka B.P. are difficult to explain with monsoon changes alone and are most likely also linked to regional hydrographic changes, such as trade wind induced variations in warm water advection.

Clay mineralogy of sediment cores from the Arctic Ocean

This study focuses on sedimentological investigations of sediment cores recovered during the international Arctic'91, expeditions with the German research ice breaker RV "Polarstern" to the European sector of the Arctic Ocean. Here, we deduce the last glacial/interglacial changes in transport mechanism and sedimentation from the clay mineral group smectite. We choose the smectites as an example of how sediment mineralogy can be linked with particular source regions (the Kara and Laptev seas), distinct transport mechanism (sea ice and surface currents) and sedimentation processes. Smectite contents in Arctic sediments discussed for two time slices, including the Last Glacial Maximum (LGM), and the last deglaciation (Termination I), reveal the highest variability subsequent to the retreat of the Eurasian ice sheets. Our results show that smectite anomalies in the Eurasian Basin are associated with distinct meltwater pulses and occurred around 13.5-13.0 14C ka B.P. Compelling evidence is provided that these anomalies are deduced from sea-ice entrained sediments from the eastern Kara Sea that entered the Arctic Ocean after ice-sheet break-up and eventually flooding of the Kara Sea. We propose that smectite anomalies in sediments of the eastern Arctic Ocean can be utilized to identify deglacial events and to help decipher configurations of the Eurasian ice sheets. The identification of smectite maxima along the modern sea-ice edge in the Eurasian Basin further indicates biologically enhanced sedimentation from melting sea ice allowing the reconstruction of seasonally open water in the region. Hence, considering the poor preservation conditions of primary paleoceanographic proxies in the Arctic Ocean, the clay mineral contents, particularly the smectite group, may be one alternative tool for paleoclimatic reconstruction in the Eurasian Basin.

(Table 1) Age determination of DSDP Hole 94-607 and sediment cores V30-101 and V23-81

Oxygen isotope measurements in Greenland ice demonstrate that a series of rapid warm-cold oscillations -called Dansgaard-Oeschger events- punctuated the last glaciation (Dansgard et al., 1993, doi:10.1038/364218a0). Here we present records of sea surface temperature from North Atlantic sediments spanning the past 90 kyr which contain a series of rapid temperature oscillations closely matching those in the ice-core record, confirming predictions that the ocean must bear the imprint of the Dansgaard-Oeschger events (Broecker et al., 1988, doi:10.1016/0033-5894(88)90082-8; 1990, doi:10.1029/PA005i004p00469). Moreover, we show that between 20 and 80 kyr ago, the shifts in ocean-atmosphere temperature are bundled into cooling cycles, lasting on average 10 to 15 kyr, with asymmetrical saw-tooth shapes. Each cycle culminated in an enormous discharge of icebergs into the North Atlantic (a 'Hein-rich event' (Bond et al., 1992, doi:10.1038/360245a0; Broecker et al., 1992, doi:10.1007/BF00193540), followed by an abrupt shift to a warmer climate. These cycles document a previously unrecognized link between ice sheet behaviour and ocean-atmosphere temperature changes. An important question that remains to be resolved is whether the cycles are driven by external factors, such as orbital forcing, or by inter-nal ice-sheet dynamics.

Sedimentation rates calculated on surface sediment samples from different site of the Atlantic and Pacific Oceans (Table 1)

Radiocarbon ages on CaCO3 from deep-sea cores offer constraints on the nature of the CaCO3 dissolution process. The idea is that the toll taken by dissolution on grains within the core top bioturbation zone should be in proportion to their time of residence in this zone. If so, dissolution would shift the mass distribution in favor of younger grains, thereby reducing the mean radiocarbon age for the grain ensemble. We have searched in vain for evidence supporting the existence of such an age reduction. Instead, we find that for water depths of more than 4 km in the tropical Pacific the radiocarbon age increases with the extent of dissolution. We can find no satisfactory steady state explanation and are forced to conclude that this increase must be the result of chemical erosion. The idea is that during the Holocene the rate of dissolution of CaCO3 has exceeded the rain rate of CaCO3. In this circumstance, bioturbation exhumes CaCO3 from the underlying glacial sediment and mixes it with CaCO3 raining from the sea surface.

(Table 1) Age determination of Arctic Ocean sediment cores

Accelerator mass spectrometer 14C dated stable isotope data from Neogloboquadrina pachyerma in cores raised from the Mendeleyev Ridge and slope provide evidence for significant influx of meltwater to the western Arctic Ocean during the early part of marine oxygen isotope stage 1 (OIS 1) and during several intervals within OIS 3. The strongest OIS 3 meltwater event occurred before ca. 45 ka (conventional radiocarbon age) and was probably related to the deglaciation at the beginning of OIS 3. Major meltwater input to the western Arctic Ocean during the last deglaciation coincides closely with the maximum rate of global sea-level rise as determined from the Barbados sea-level record, demonstrating a strong link between the global record and changes in the central Arctic Ocean. OIS 2, which includes the last glacial maximum, is very condensed or absent in the cores. Abundance and d13C values for N. pachyderma in the middle part of OIS 3 are similar to modern values, indicating high productivity and seasonal ice-free areas along the Arctic margin at that time. These records indicate that the Arctic Ocean was a source of heat and moisture to the northern polar atmosphere during parts of OIS 3.

Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.