Pollen and dinoflagellate cyst counts on sediment core M72/5_628-1 (MD72/5-25-GC1) from the Black Sea

High-resolution pollen and dinoflagellate cyst records from sediment core M72/5-25-GC1 were used to reconstruct vegetation dynamics in northern Anatolia and surface conditions of the Black Sea between 64 and 20 ka BP. During this period, the dominance of Artemisia in the pollen record indicates a steppe landscape and arid climate conditions. However, the concomitant presence of temperate arboreal pollen suggests the existence of glacial refugia in northern Anatolia. Long-term glacial vegetation dynamics reveal two major arid phases ~64-55 and 40-32 ka BP, and two major humid phases ~54-45 and 28-20 ka BP, correlating with higher and lower summer insolation, respectively. Dansgaard-Oeschger (D-O) cycles are clearly indicated by the 25-GC1 pollen record. Greenland interstadials are characterized by a marked increase in temperate tree pollen, indicating a spread of forests due to warm/wet conditions in northern Anatolia, whereas Greenland stadials reveal cold and arid conditions as indicated by spread of xerophytic biomes. There is evidence for a phase lag of ~500 to 1500 yr between initial warming and forest expansion, possibly due to successive changes in atmospheric circulation in the North Atlantic sector. The dominance of Pyxidinopsis psilata and Spiniferites cruciformis in the dinocyst record indicates brackish Black Sea conditions during the entire glacial period. The decrease of marine indicators (marine dinocysts, acritarchs) at ~54 ka BP and increase of freshwater algae (Pediastrum, Botryococcus) from 32 to 25 ka BP reveals freshening of the Black Sea surface water. This freshening is possibly related to humid phases in the region, to connection between Caspian Sea and Black Sea, to seasonal freshening by floating ice, and/or to closer position of river mouths due to low sea level. In the southern Black Sea, Greenland interstadials are clearly indicated by high dinocyst concentrations and calcium carbonate content, as a result of an increase in primary productivity. Heinrich events show a similar impact on the environment in the northern Anatolia/Black Sea region as Greenland stadials.

Database of the GLAMAP project

A uniform chronology for foraminifera-based sea surface temperature records has been established in more than 120 sediment cores obtained from the equatorial and eastern Atlantic up to the Arctic Ocean. The chronostratigraphy of the last 30,000 years is mainly based on published d18O records and 14C ages from accelerator mass spectrometry, converted into calendar-year ages. The high-precision age control provides the database necessary for the uniform reconstruction of the climate interval of the Last Glacial Maximum within the GLAMAP-2000 project.

(Table S1) Neodymium isotopic results from ODP Site 172-1063

The Atlantic Meridional Overturning Circulation (AMOC) plays an important role in the Northern Hemisphere climate system. Significant interest went into the question of how excessive freshwater input through melting of continental ice can affect its overturning vigor and, hence, heat supply, to higher northern latitudes. Such forcing can be tested by investigating its behavior during extreme iceberg discharge events into the open North Atlantic during the last glacial period, the so-called Heinrich events (HE). Here we present neodymium (Nd) isotope compositions of past seawater, a sensitive chemical water mass tag, extracted from sediments of Ocean Drilling Program Site 1063 in the western North Atlantic (Bermuda Rise), covering the period surrounding HE 2, the Last Glacial Maximum, and the early deglaciation. These data are compared with a record of the kinematic circulation tracer (231Pa/230Th)xs extracted from the same sediment core. Both tracers indicate significant circulation changes preceding intense ice rafting during HE 2 by almost 2 kyr. Moreover, the Nd isotope record suggests the presence of deeply ventilating North Atlantic Deep Water early during Marine Isotope Stage 2 until it was replaced by Southern Source Water at ~27 ka. The early switch to high (Pa/Th)xs and radiogenic epsilon-Nd in relation to intensified ice rafting during HE 2 suggests that ice rafting into the open North Atlantic during major HE 2 was preceded by an early change of the AMOC. This opens the possibility that variations in AMOC contributed to or even triggered the ice sheet instability rather than merely responding to it.

CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

We present new d13C measurements of atmospheric CO2 covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in d13C(atm) of 0.5 permil occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in d13C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO2]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of d13C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

Oceanographic investiations from the Sea of Ochotsk

The Sea of Okhotsk is a marginal sea of the Pacific Ocean, which is characterized by strong variations in the productivity and sediment supply due to sea ice transport and river input. Furthermore the variations in the hydrological cycle determine the formation of the SOIW (Sea of Okhotsk Intermediate Water) which plays an important role in the ventilation processes in the intermediate water of the N-Pacific. Isotope data measured on planktonic and benthic foraminifera, sedimentological and geochemical studies of sediment cores and surface samples from the Sea of Okhotsk are used to reconstruct the paleoceanography during the past 350.000 years. The dating and correlation of the sediments are based on oxygen isotope stratigraphy, absolute ages, magnetic susceptibility as well as a detailled tephrachronology of the entire basin. The sedimentation rates are characterized by temporal and spatial variations. The maximum sedimentation rate takes place at the continental slope off Sakhalin due to the input of the Amur River, the sea ice drift and the high productivity. The sedimentation rate in the eastern part of the Sea of Okhotsk is generelly high because of the influence of the nutrient-rich Kamchatka Current. In the central and northern parts of the Sea of Okhotsk, areas with low productivity and reduced terrestrial supply, the sedimentation rate is the lowest. The analyses of the surface sediment samples make it possible to characterize the (sub)- recent sediment supply and transportation processes. The bulk sediment measurements, isotope data and the accumulation rate of ice-rafted debris (IRD) show a dominant sea ice cover and a region with a high productivity as well as a high Amur River input in the western part of the sea. The eastern part of the Sea of Okhotsk, however, is marked by the predominance of warm and nutrient-rich water masses coming from the Kamchatka Current which restricts the sea ice cover. This is reflected in low content of ice-rafted debris and high productivity proxies as well as in isotope data. The deposits of the Sea of Okhotsk are characterized by terrestrial, biogenic and volcanogenic sediment input which varies temporally and spatially. Here, the sedimentation pattern is dominated by the terrestrial input. Bulk sediment measurements and sample analyses of the > 63 micron particle input make it possible to distinguish glacial and interglacial fluctuations. The sedimentation processes during glacial times are determined by a high content of ice-rafted debris, whereas the primary production is higher during interglacial periods. During the last glacial/interglacial cycle the IRD-distribution pattern indicates a strong sea ice transport in the western part and in large areas of the open sea in the eastern part of the Sea of Okhotsk with a relatively constant ice-drift system. The IRD flux in sediments of the oxygen isotope Stage 6 reflects a new sedimentation pattern in the eastern part of the sea. This high IRD accumulation rate indicates ice advances beyond the shelf margin and an iceberg transport from NE-E direction into the Sea of Okhotsk. The several large, brief, negative anomalies in d13C values of Neogloboquadrina pachyderma (s) show releases of methane from basin sediments which correspond to periods of relative sea level falls. The high sedimentation rates on the Sakhalin slope allow insights into the climatic history in Holocene and indicate shorter-scale variations oscillation in Stage 3, which correlate with the global climatic changes. These variations are described as Dansgaard-Oeschger cycles in Greenland ice cores and as Heinrich-Events in several marine sediment cores from the N-Atlantic.

Compilation of 220 sediment core d13C data from the glacial Atlantic Ocean

We compare a compilation of 220 sediment core d13C data from the glacial Atlantic Ocean with three-dimensional ocean circulation simulations including a marine carbon cycle model. The carbon cycle model employs circulation fields which were derived from previous climate simulations. All sediment data have been thoroughly quality controlled, focusing on epibenthic foraminiferal species (such as Cibicidoides wuellerstorfi or Planulina ariminensis) to improve the comparability of model and sediment core carbon isotopes. The model captures the general d13C pattern indicated by present-day water column data and Late Holocene sediment cores but underestimates intermediate and deep water values in the South Atlantic. The best agreement with glacial reconstructions is obtained for a model scenario with an altered freshwater balance in the Southern Ocean that mimics enhanced northward sea ice export and melting away from the zone of sea ice production. This results in a shoaled and weakened North Atlantic Deep Water flow and intensified Antarctic Bottom Water export, hence confirming previous reconstructions from paleoproxy records. Moreover, the modeled abyssal ocean is very cold and very saline, which is in line with other proxy data evidence.

(Table S1) Age determination of Atlantic Ocean sediment cores

Well-dated benthic foraminifer oxygen isotopic records (d18O) from different water depths and locations within the Atlantic Ocean exhibit distinct patterns and significant differences in timing over the last deglaciation. This has two implications: on the one hand, it confirms that benthic d18O cannot be used as a global correlation tool with millennial-scale precision, but on the other hand, the combination of benthic isotopic records with independent dating provides a wealth of information on past circulation changes. Comparing new South Atlantic benthic isotopic data with published benthic isotopic records, we show that (1) circulation changes first affected benthic d18O in the 1000-2200 m range, with marked decreases in benthic d18O taking place at ~17.5 cal. kyr B.P. (ka) due to the southward propagation of brine waters generated in the Nordic Seas during Heinrich Stadial 1 (HS1) cold period; (2) the arrival of d18O-depleted deglacial meltwater took place later at deeper North Atlantic sites; (3) hydrographic changes recorded in North Atlantic cores below 3000 m during HS1 do not correspond to simple alternations between northern- and southern-sourced water but likely reflect instead the incursion of brine-generated deep water of northern as well as southern origin; and (4) South Atlantic waters at ~44°S and ~3800 m depth remained isolated from better-ventilated northern-sourced water masses until after the resumption of North Atlantic Deep Water (NADW) formation at the onset of the Bølling-Allerod, which led to the propagation of NADW into the South Atlantic.

Temperature variabilty in the southeastern Black-Sea during 64-20ka BP

The Eurasian inland propagation of temperature anomalies during glacial millennial-scale climate variability is poorly understood but this knowledge is crucial to understanding hemisphere-wide atmospheric teleconnection patterns and climate mechanisms. Based on biomarkers and geochemical paleothermometers, a pronounced continental temperature variability between 64,000 and 20,000 years ago, coinciding with the Greenland Dansgaard-Oeschger cycles, was determined in a well-dated sediment record from the formerly enclosed Black Sea. Cooling during Heinrich events was not stronger than during other stadials in the Black Sea. This is corroborated by modeling results showing that regular Dansgaard-Oeschger cycles penetrated deeper into the Eurasian continent than Heinrich events. The pattern of coastal ice-rafted detritus suggests a strong dependence on the climate background state, with significantly milder winters during periods of reduced Eurasian ice sheets and an intensified meridional atmospheric circulation.

(Table S1) Age determination of sediment cores SK99_17 and MD76-131

Productivity in the Arabian Sea is one of the highest in the world. It is controlled by seasonally reversing monsoonal wind-driven upwelling of nutrient-rich deeper waters which fuel phytoplankton growth. The detailed history of upwelling-induced productivity in the eastern Arabian Sea is unknown. Here we present paleoproductivity records from a composite sediment core at the millennial scale during the last 80 kyr B.P. These records are based on relative abundance counts of planktonic foraminifera and organic carbon contents, which are shown to mainly vary in concert. The eastern Arabian Sea upwelling-induced productivity was higher in the glacial period than in the Holocene, but it fell repeatedly on millennial timescales. These productivity declines occurred during cold events in the North Atlantic region, with the most pronounced changes prevailing during the Heinrich events. Hence, seasonal monsoon winds that drive upwelling-induced productivity in the east were weak when the North Atlantic was cold. These weak winds resulted in stratification of the water column, comparable to today's Arabian Sea stratification in the intermonsoonal period. Combining the new eastern with published western Arabian Sea results shows that the entire biological factory was severely diminished during the North Atlantic Heinrich events, and the seasonal productivity change in the Arabian Sea monsoon system was reduced with year-round low productivity.

Dinocyst-derived sea-surface characteristics, palynomorph abundances, ice-rafted debris, and planktonic foraminifera abundances and d18O of the 35-41 ka BP interval of core MD99-2285

The last glacial millennial climatic events (i.e. Dansgaard-Oeschger and Heinrich events) constitute outstanding case studies of coupled atmosphere-ocean-cryosphere interactions. Here, we investigate the evolution of sea-surface and subsurface conditions, in terms of temperature, salinity and sea ice cover, at very high-resolution (mean resolution between 55 and 155 years depending on proxies) during the 35-41 ka cal BP interval covering three Dansgaard-Oeschger cycles and including Heinrich event 4, in a new unpublished marine record, i.e. the MD99-2285 core (62.69°N; -3.57s°E). We use a large panel of complementary tools, which notably includes dinocyst-derived sea-ice cover duration quantifications. The high temporal resolution and multiproxy approach of this work allows us to identify the sequence of processes and to assess ocean-cryosphere interactions occurring during these periodic ice-sheet collapse events. Our results evidence a paradoxical hydrological scheme where (i) Greenland interstadials are marked by a homogeneous and cold upper water column, with intensive winter sea ice formation and summer sea ice melting, and (ii) Greenland and Heinrich stadials are characterized by a very warm and low saline surface layer with iceberg calving and reduced sea ice formation, separated by a strong halocline from a less warm and saltier subsurface layer. Our work also suggests that this stadial surface/subsurface warming started before massive iceberg release, in relation with warm Atlantic water advection. These findings thus support the theory that upper ocean warming might have triggered European ice-sheet destabilization. Besides, previous paleoceanographic studies conducted along the Atlantic inflow pathways close to the edge of European ice-sheets suggest that such a feature might have occurred in this whole area. Nonetheless, additional high resolution paleoreconstructions are required to confirm such a regional scheme.

(Table 2) Age determination of Irminger Basin sediment cores

High-resolution records of coarse lithic content and oxygen isotope have been obtained in a piston core from the Irminger Basin. The last glacial period is characterized by numerous periods of increased iceberg discharges originating partly from Iceland and corresponding to millennial-scale instabilities of the coastal ice sheets and ice shelves in the Nordic area. A comparison with midlatitude sediment cores shows that ice-rafted material corresponding to the Heinrich events was deposited synchronously from 40° to 60°N. There are thus two oscillating systems: every 5-10 kyr massive iceberg armadas are released from large continental ice caps, whereas more frequent instabilities of the coastal ice sheets in the high latitude regions occur every 1.2-3.8 kyr. At the time of the Heinrich events the synchroneity of the response from all the northern hemisphere ice sheets attests the existence of strong interactions between the two systems.

(Table 1) Age determination of sediment cores EW9302-JPC1 and EW9302-JPC2

Results from two deep sea cores from northeast of Newfoundland at 1251 and 2527 m water depth, respectively, indicate that during the time period from 160,000 to 10,000 years BP, ice rafting events in the Labrador Sea were accompanied by rapid variations in deep and surface water circulation. Twelve ice-rafting events occurred, each coinciding with high concentrations of detrital carbonate and oxygen isotopic depletion of both surface and bottom waters. Eleven of these can be correlated with the North Atlantic Heinrich events H1-H11. The remaining very conspicuous ice-rafting event took place early in MIS substage 5e, at a time when the planktic faunal assemblage suggests marked warming of the sea surface. In the shallower core, benthic d13C values rise from a minimum during the deglaciation to peak substage 5e values following the last ice-rafting event, indicating that the ventilation of intermediate depths was renewed after the deglaciation was complete and continued throughout substage 5e. The benthic foraminifera suggest that this well-ventilated water mass was comparable to the modern Labrador Sea Water (LSW). The benthic faunas suggest that a relatively warm intermediate water mass entered the SE Labrador Sea during Heinrich events. Generally low benthic d13C values indicate that this water mass was poorly ventilated and rich in inorganic nutrients. Isotope data and benthic faunal distributions indicate that North Atlantic Deep Water (NADW) formed in the Norwegian-Greenland Sea reached the SE Labrador Sea between the Heinrich events.

Paleotemperature and productivity reconstruction of sediment core D11957P on the Iberian Margin

Foraminifera counts and climatic assemblages from the Tore Seamount are used to approach the glacial and interglacial changes in temperature and productivity on the Iberian Margin over the last 225 kyr. Chronostratigraphy is based on Globigerinoides ruber and Globigerina bulloides oxygen isotopes and supported by foraminifera and carbonate stadial fluctuations. Foraminifera indicate cooling from late interglacial stage 5 to the beginning of Termination I (TI). Neogloboquadnna pachyderma-s reflects cold conditions during glacial stages 4-2. In contrast, glacial stage 6 is dominated by warmer N. pachyderma-d and dutertrei and a restricted arctic assemblage. Past sea surface temperatures confirm the general cooling, reaching 4.3°C (SIMMAX.28) during stage 2. Multiple productivity proxies such as organic carbon, productivity-related foraminifera, and delta13C constrain the changes observed. A productivity increase occurs after interglacial stage 5, enhanced from late glacial stage 3 to TI Present-day satellite-detected phytoplankton plumes off Portugal would have accounted in the past glacial stages for the general productivity increase over the Tore. On top of this, welldefined peaks of organic carbon and productivity-related foraminifera correspond with Heinrich events 1-4.

(Table 1) Age determination of sediment core BC5, Golf of Lion

Hitherto unknown abundance peaks of left coiling (l.c.) Neogloboquadrina pachyderma from a Gulf of Lions piston core indicate that abrupt cold spells associated with Atlantic Heinrich events affected the Mediterranean. N. pachyderma (l.c.) is typical of (sub) polar waters in the open ocean. The southern edge of its glacial North Atlantic bioprovince reached south Portugal. Only trace abundances of N. pachyderma (l.c.) are known from Quaternary Mediterranean sediments, suggesting that no significant "invasions" occured via the Strait of Gibraltar. The Gulf of Lions abundance peaks therefore seem to reflect area-specific thriving of a normally rare but indigenous taxon in the western Mediterranean through local favorable habitat development. The general planktonic foraminiferal record suggests that the basic hydrographic regime in the Gulf of Lions, with wintertime deep convective overturn, was relatively stable over the past 60 kyr. Under these conditions, high abundances of N. pachyderma (l.c.) would essentially imply temperature reductions of the order of 5°-8° relative to the present.