Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Boron, carbonate carbon and oxygen isotope analyses of samples from the Yangtze Platform (China), Congo craton (Namibia), and Karatau microcontinent (Kazakhstan) of the Marinoan glaciation

Boron isotope patterns preserved in cap carbonates deposited in the aftermath of the younger Cryogenian (Marinoan, ca. 635 Ma) glaciation confirm a temporary ocean acidification event on the continental margin of the southern Congo craton, Namibia. To test the significance of this acidification event and reconstruct Earth's global seawater pH states at the Cryogenian-Ediacaran transition, we present a new boron isotope data set recorded in cap carbonates deposited on the Yangtze Platform in south China and on the Karatau microcontinent in Kazakhstan. Our compiled d11B data reveal similar ocean pH patterns for all investigated cratons and confirm the presence of a global and synchronous ocean acidification event during the Marinoan deglacial period, compatible with elevated postglacial pCO2 concentrations. Differences in the details of the ocean acidification event point to regional distinctions in the buffering capacity of Ediacaran seawater.

(Fig. 1D, 3E) Sea surface temperature and alkenone C37:4 for sediment core MD01-2444 during the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport.

Diatom oxygen and silicon isotopes data from IODP Hole 323-U1341B in the Bering Sea through the late Pliocene and over the onset of major Northern Hemisphere Glaciation

The Bering Sea is one of the most biologically productive regions in the marine system and plays a key role in regulating the flow of waters to the Arctic Ocean and into the subarctic North Pacific Ocean. Cores from Integrated Ocean Drilling Program (IODP) Expedition 323 to the Bering Sea provide the first opportunity to obtain reconstructions from the region that extend back to the Pliocene. Previous research at Bowers Ridge, south Bering Sea, has revealed stable levels of siliceous productivity over the onset of major Northern Hemisphere Glaciation (NHG) (circa 2.85-2.73 Ma). However, diatom silica isotope records of oxygen (d18Odiatom) and silicon (d30Sidiatom) presented here demonstrate that this interval was associated with a progressive increase in the supply of silicic acid to the region, superimposed on shift to a more dynamic environment characterized by colder temperatures and increased sea ice. This concluded at 2.58 Ma with a sharp increase in diatom productivity, further increases in photic zone nutrient availability and a permanent shift to colder sea surface conditions. These transitions are suggested to reflect a gradually more intense nutrient leakage from the subarctic northwest Pacific Ocean, with increases in productivity further aided by increased sea ice- and wind-driven mixing in the Bering Sea. In suggesting a linkage in biogeochemical cycling between the south Bering Sea and subarctic Northwest Pacific Ocean, mainly via the Kamchatka Strait, this work highlights the need to consider the interconnectivity of these two systems when future reconstructions are carried out in the region.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Accumulation rates of volcanic glass and ice-rafted debris in sediments of the North Pacific Ocean

Mass accumulation rates (MAR) of different components of North Pacific deep-sea sediment provide detailed information about the timing of the onset of major Northern Hemisphere glaciation that occurred at 2.65 Ma. An increase in explosive volcanism in the Kamchatka-Kurile and Aleutian arcs occured at this same time, suggesting a link between volcanism and glaciation. Sediments recovered by piston-coring techniques during ODP Leg 145 provide a unique opportunity to undertake a detailed test of this possibility. Here we use volcanic glass as a proxy for explosive volcanism and ice-rafted debris (IRD) as a proxy for glaciation. The MAR of both glass and IRD increase markedly at 2.65 Ma. Further, the flux of the volcanic glass increased just prior the flix of ice-radted material, suggesting that the cooling resulting from explosive volcanic eruptions may have been the ultimate trigger for the mid-Pliocene glacial intensification.

Dinoflagellate cyst assemblages, oxygen isotope records of foraminifera, and alkenone and Mg/Ca-based SST estimates for the time interval 3.40-3.18 Ma (Late Pliocene) from five North Atlantic and Caribbean sediment cores

The early Late Pliocene (3.6 to ~3.0 million years ago) is the last extended interval in Earth's history when atmospheric CO2 concentrations were comparable to today's and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ~3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Our geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current (NAC) and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the Greenland ice sheet during MIS M2, despite near-modern atmospheric CO2 concentrations. Before and after MIS M2, heat transport via the NAC was crucial in maintaining warm climates comparable to those predicted for the end of this century.

Stable isotope analysis and temperature reconstruction data from DSDP Hole 94-609B and ODP Hole 162-984B

The Greenland ice sheet is accepted as a key factor controlling the Quaternary glacial scenario. However, the origin and mechanisms of major Arctic glaciation starting at 3.15 Ma and culminating at 2.74 Ma are still controversial. For this phase of intense cooling Ravelo et al. proposed a complex gradual forcing mechanism. In contrast, our new submillennial-scale paleoceanographic records from the Pliocene North Atlantic suggest a far more precise timing and forcing for the initiation of northern hemisphere glaciation (NHG), since it was linked to a 2-3 °C surface water warming during warm stages from 2.95 to 2.82 Ma. These records support previous models, claiming that the final closure of the Panama Isthmus (3.0- ~2.5 Ma induced an increased poleward salt and heat transport. Associated strengthening of North Atlantic Thermohaline Circulation and in turn, an intensified moisture supply to northern high latitudes resulted in the build-up of NHG, finally culminating in the great, irreversible climate crash at marine isotope stage G6 (2.74 Ma). In summary, there was a two-step threshold mechanism that marked the onset of NHG with glacial-to-interglacial cycles quasi-persistent until today.

Sedimentological and clay mineralogical analysis of ODP sites in the central Fram Strait

This study presents the results of high-resolution sedimentological and clay mineralogical investigations on sediments from ODP Sites 908A and 909AlC located in the central Fram Strait. The objective was to reconstruct the paleoclimate and paleoceanography of the high northern latitudes since the middle Miocene. The sediments are characterised in particular by a distinctive input of ice-rafted material, which most probably occurs since 6 Ma and very likely since 15 Ma. A change in the source area at 1 1.2 Ma is clearly marked by variations within clay mineral composition and increasing accumulation rates. This is interpreted as a result of an increase in water mass exchange through the Fram Strait. A further period of increasing exchange between 4-3 Ma is identified by granulometric investigations and points to a synchronous intensification of deep water production in the North Atlantic during this time interval. A comparison of the components of coarse and clay fraction clearly shows that both are not delivered by the Same transport process. The input of the clay fraction can be related to transport mechanisms through sea ice and glaciers and very likely also through oceanic currents. A reconstruction of source areas for clay minerals is possible only with some restrictions. High smectite contents in middle and late Miocene sediments indicate a background signal produced by soil formation together with sediment input, possibly originating from the Greenland- Scotland Ridge. The applicability of clay mineral distribution as a climate proxy for the high northern latitudes can be confirmed. Based on a comparison of sediments from Site 909C, characterised by the smectite/illite and chlorite ratio, with regional and global climatic records (oxygen isotopes), a middle Miocene cooling phase between 14.8-14.6 Ma can be proposed. A further cooling phase between 10-9 Ma clearly shows similarities in its Progress toward drastic decrease in carbonate sedimentation and preservation in the eastern equatorial Pacific. The modification in sea water and atmosphere chemistry may represent a possible link due to the built-up of equatorial carbonate platforms. Between 4.8-4.6 Ma clay mineral distribution indicates a distinct cooling trend in the Fram Strait region. This is not accompanied by relevant glaciation, which would otherwise be indicated by the coarse fraction. The intensification of glaciation in the northern hemisphere is distinctly documented by a rapid increase of illite and chlorite starting from 3.3 Ma, which corresponds to oxygen isotope data trends from North Atlantic.

Facies analysis on sediment core CRP-3 from the Ross Sea, Antarctica

Cape Roberts Project drill core 3 (CRP-3) was obtained from Roberts ridge, a sea-floor high located at 77°S, 12 km offshore from Cape Roberts in western McMurdo Sound, Antarctica. The recovered core is about 939 m long and comprises strata dated as being early Oligocene (possibly latest Eocene) in age, resting unconformably on ~116 m of basement rocks consisting of Palaeozoic Beacon Supergroup sediments. The core includes ten facies commonly occurring in five major associations that are repeated in particular sequences throughout the core and which are interpreted as representing different depositional environments through time. Depositional systems inferred to be represented in the succession include: outer shelf, inner shelf, nearshore to shoreface each under iceberg influence, deltaic and/or grounding-line fan, and ice proximal-ice marginal-subglacial (mass flow/rainout diamictite/subglacial till) singly or in combination. The record is taken to represent the initial talus/alluvial fan setting of a glaciated rift margin adjacent to the block-uplifted Transantarctic Mountains. Development of a deltaic succession upcore was probably associated with the formation of palaeo-Mackay valley with temperate glaciers in its headwaters. At that stage glaciation was intense enough to support glaciers ending in the sea elsewhere along the coast, but a local glacier was fluctuating down to the sea by the time the youngest part of CRP-3 was being deposited. Changes in palaeoenvironmental interpretations in this youngest part of the core are used to estimate relative glacial proximity to the drillsite through time. These inferred glacial fluctuations are compared with the global d180 and Mg/Ca curves to evaluate the potential of glacial fluctuations on Antarctica for influencing these records of global change. Although the comparisons are tentative at present, the records do have similarities, but there are also some differences that require further evaluation.

Neogene planktonic foraminifera of ODP Site 130-806

The Neogene biostratigraphy presented here is based on the study of 230 samples through 737 m of pelagic sediment in Hole 806B. Sediment accumulation is interrupted only once in the uppermost lower Miocene (Zone N6), apparently coincident with a widespread deep-sea hiatus. Preservation of planktonic foraminifers through the section ranges from good to moderately poor. One hundred and ten species of planktonic foraminifers were identified; taxonomic notes on most species are included. All of the standard low-latitude Neogene foraminiferal zones are delineated, with the exceptions of Zones N8 and N9 because of a high first occurrence of Orbulina, and Zones N18 and N19 because of a high first occurrence of Sphaeroidinella dehiscens. Good agreement exists between the published account of the variation in planktonic foraminiferal species richness and the rates of diversification and turnover, and measurements of these evolutionary indexes in the record of Hole 806B. The global pattern of change in tropical/transitional species richness is paralleled in Hole 806B, with departures caused by either ecological conditions peculiar to the western equatorial Pacific or by inexactness in the estimation of million-year intervals in Hole 806B. Temporal changes in the relative abundance of taxa in the sediment assemblages, considered in light of their depth habitats, reveal a detailed picture of historical change in the structure of the upper water column over the Ontong Java Plateau. The dominance of surface dwellers (Paragloborotalia kugleri, P. mayeri, Dentoglobigerina altispira, Globigerinita glutinata, and Globigerinoides spp.) throughout the lower and middle Miocene is replaced by a more equitable distribution of surface (D. altispira and Globigerinoides spp.), intermediate (Globorotalia menardii plexus), and deep (Streptochilus spp.) dwellers in the late Miocene, following the closing of the Indo-Pacific Seaway and the initiation of large-scale glaciation in the Antarctic. The shoaling of the thermocline along the equator engendered by these climatic and tectonic events persisted through the Pliocene, when initial increases in the abundance of a new set of shallow, intermediate, and deep dwelling species of planktonic foraminifers coincide with the closing of the Panamanian Seaway.

(Table 1) Maceral composition and organic carbon chemistry of ODP Leg 104 samples

The Cenozoic sediments sampled in ODP Leg 104 on the Vøring Plateau show a distinct variability of the total organic carbon content (TOC) and the accumulation rates of TOC. Based on the geochemical and organic-petrographic characterization of the sedimentary organic matter (OM), the allochthonous and autochthonous proportion of the OM could be quantified. The results clearly demonstrate that high TOC percentages and TOC accumulation rates in Cenozoic sediment sections display a generally high input of allochthonous organic matter. Oxidized and partly well-rounded organic particles built up the main portion of OM within the Miocene, TOC-rich sediments. The most probable source of this oxidized OM are reworked sediments from the Scandinavian shelf. Changes in the input of these organic particles are to some degree correlative with sea-level changes. The Cenozoic accumulation of autochthonous OM is low and does not reveal a clear variation during the Miocene and early Pliocene. In spite of a high accumulation rate of biogenic opal during the Early Miocene, the accumulation rate of autochthonous TOC is low. The autochthonous particle assemblage is dominated by relatively inert OM, like dinoflagellate cysts. This points to an intensive biological and/or early diagenetic degradation of the marine OM under well oxidized bottom water conditions during the last 23 Myr. Nevertheless, a continuation of marine OM degradation during later stages of diagenesis cannot be excluded. A prominent dominance of allochthonous OM over autochthonous is documented with the beginning of the Pliocene. At 2.45 Ma the episodic occurrence of ice-rafted, thermally mature OM reflects the onset of the glacial erosion of Mesozoic, coal and black shale bearing sediments on the Scandinavian and Barents Sea shelves. The first occurrence of these, in view of the actual burial depth, thermally overmature OM particles is, therefore, a marker for the beginning of the strong Scandinavian glaciation and the advance of the glacial front toward the shelves.

Elevation lines of modelled ice sheet extension during last glacial maximum

A numerical ice-sheet model was used to reconstruct the Late Weichselian glaciation of the Eurasian High Arctic, between Franz Josef Land and Severnaya Zemlya. An ice sheet was developed over the entire Eurasian High Arctic so that ice flow from the central Barents and Kara seas toward the northern Russian Arctic could be accounted for. An inverse approach to modeling was utilized, where ice-sheet results were forced to be compatible with geological information indicating ice-free conditions over the Taymyr Peninsula during the Late Weichselian. The model indicates complete glaciation of the Barents and Kara seas and predicts a "maximum-sized" ice sheet for the Late Weichselian Russian High Arctic. In this scenario, full-glacial conditions are characterized by a 1500-m-thick ice mass over the Barents Sea, from which ice flowed to the north and west within several bathymetric troughs as large ice streams. In contrast to this reconstruction, a "minimum" model of glaciation involves restricted glaciation in the Kara Sea, where the ice thickness is only 300 m in the south and which is free of ice in the north across Severnaya Zemlya. Our maximum reconstruction is compatible with geological information that indicates complete glaciation of the Barents Sea. However, geological data from Severnaya Zemlya suggest our minimum model is more relevant further east. This, in turn, implies a strong paleoclimatic gradient to colder and drier conditions eastward across the Eurasian Arctic during the Late Weichselian.