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.

(Table 1) Thermoluminescence datings of glacial deposits from several sites in Schleswig-Holstein, north Germany

First thermoluminescence (TL) datings of glacial deposits from several well-known sites in northwest Germany (Schleswig-Holstein) indicate a possible Middle Weichselian Glaciation in this region. The TL dates obtained imply an ice sheet at about 70-40 ka B.P. here, whose presence in this part of Europe has not been previously considered. The reliability of the TL dates was cross-checked on different lithologies; the dating was performed in two independent laboratories. These dates and their interpretations, however, must still be confirmed both geologically and paleontologically on suitable sites where glacial sediments are sandwiched between Eemian deposits and interstadial deposits older than the Late Weichselian.

High-resolution multi-proxy record of the last glacial period from Lake Van, eastern Anatolian high plateau, Turkey

Detailed analyses of the Lake Van pollen, Ca/K ratio and stable oxygen isotope record allow the identification of millennial-scale vegetation and environmental changes in eastern Anatolia throughout the last glacial (~75-15 ka BP). The climate within the last glacial was cold and dry, with low arboreal pollen (AP) levels. The driest and coldest period corresponds to Marine Isotope Stage (MIS) 2 (~28-14.5 ka BP) dominated by the highest values of xerophytic steppe vegetation. Our high-resolution multi proxy record shows rapid expansions and contractions of tree populations that reflects variability in temperature and moisture availability. This rapid vegetation and environmental changes can be linked to the stadial-interstadial pattern of the Dansgaard-Oeschger (DO) events as recorded in the Greenland ice cores. Periods of reduced moisture availability were characterized by enhanced xerophytic species and high terrigenous input from the Lake Van catchment area. Furthermore, comparison with the marine realm reveals that the complex atmosphere-ocean interaction can be explained by the strength and position of the westerlies, which is responsible for the supply of humidity in eastern Anatolia. Influenced by diverse topography of the Lake Van catchment, larger DO interstadials (e.g. DO 19, 17-16, 14, 12 and 8) show the highest expansion of temperate species within the last glacial. However, Heinrich events (HE), characterized by highest concentrations of ice-rafted debris (IRD) in marine sediments, are identified in eastern Anatolia by AP values not lower and high steppe components not more abundant than during DO stadials. In addition, this work is a first attempt to establish a continuous microscopic charcoal record over the last glacial in the Near East, which documents an initial immediate response to millennial-scale climate and environmental variability and enables us to shed light on the history of fire activity during the last glacial.

Deep thermohaline Circulation in the low-latitude Atlantic during the Last Glacial

Present-day low-latitude eastern and western Atlantic basins are geochemically distinct below the sill depth of the Mid-Atlantic Ridge. While Antarctic Bottom Water (AABW) circulates freely in the western Atlantic, flow into the eastern Atlantic is restricted below 4 km which results in filling the abyssal depths of this basin with water of geochemical similarity to nutrient depleted North Atlantic Deep Water. Using carbon isotopes and Cd/Ca ratios in benthic foraminifera we reconstruct the geochemistry of these basins during the last glacial maximum. Results indicate that deep eastern and western Atlantic basins became geochemically identical during the last glacial. This was achieved by shoaling of the upper surface of AABW above the sill depth of the Mid-Atlantic Ridge, which allowed bottom waters in both basins to be filled with the same water mass. Although AABW became the dominant water mass in the deep eastern Atlantic basin during the glacial, Holocene-glacial delta13C-PO4 shifts in this basin are in Redfield proportions, unlike the disproportionate Holocene-glacial delta13C-PO4 shifts observed in the Southern Ocean. By examining the composition of deep and intermediate waters throughout the Atlantic, we show that this effect was induced by a change in gradient of the delta13C-PO4 deepwater mixing line during glacial times. Evidence from high-latitude planktonic data suggests that the change in gradient of the deepwater mixing line was brought about through a significant reduction in the thermodynamic effect on Southern Ocean surface waters. By using coupled delta13C-PO4 data to constrain the composition of end member water masses in the glacial Atlantic, we conclude that deep waters in the low-latitude glacial Atlantic were composed of a mixture of northern and southern source waters in a ratio of 1:3.

(Table 3) Sedimentation rates of stage 1 and last glacial maximum from different Holes during cruises of RRS Discovery in 1984-1985 and 1987-1988, and RRS James Clark Ross in 1992-1993

Within the Scotia Sea, the axis of the Antarctic Circumpolar Current (ACC) is geographically confined, and sediments therefore contain a record of palaeo-flow speed uncomplicated by ACC axis migration. We outline Holocene and Last Glacial Maximum (LGM) current-controlled sedimentation using data from 3.5-kHz profiles, cores and current meter moorings. Geophysical surveys show areas of erosion and deposition controlled by Neogene basement topography. Deposition occurs in mounded sediment drifts or flatter areas, where 500-1000 m of sediment overlies acoustic basement. 3.5-kHz profiles show parallel, continuous sub-bottom reflectors with highest sedimentation rates in the centre of the drifts, and reflectors converging towards marginal zones of non-deposition. Locally, on the flanks of continental blocks (e.g. South Georgia), downslope processes are dominant. The absence of mudwaves on the sediment drifts may result from the unsteadiness of ACC flow. A core transect from the ACC axis south to the boundary with the Weddell Gyre shows a southward decrease in biogenic content, controlled by the Polar Front and the spring sea-ice edge. Both these features lay farther north at LGM. The cores have been dated by relative abundance of the radiolarian Cycladophora davisiana, and by changes in the biogenic Ba content, a palaeoproductivity indicator. Sedimentation rates range from 3 to 17 cm/ka. The grain size of Holocene sediments shows a coarsening trend from south to north, consistent with strongest bottom-current flow near the ACC axis, though interpretation is complicated by the presence of biogenic grains. Year-long current meter records indicate mean speeds from 7 cm/s in the south to 12 cm/s in the north, with benthic storm frequency increasing northwards. LGM sediments are predominantly terrigenous and show a clearer northward-coarsening trend, with well-sorted silts in the northern Scotia Sea. Assuming a constant terrigenous source, this implies stronger ACC flow at the LGM, contrasting with weaker Weddell Gyre flow deduced from earlier work.

Sedimentology and geochemistry of late glacial and Holocene sediments from Laguna Potrok Aike, Patagonia

A high-resolution multiproxy geochemical approach was applied to the sediments of Laguna Potrok Aike in an attempt to reconstruct moist and dry periods during the past 16 000 years in southeastern Patagonia. The age-depth model is inferred from AMS 14C dates and tephrochronology, and suggests moist conditions during the Lateglacial and early Holocene (16 000-8700 cal. BP) interrupted by drier conditions before the beginning of the Holocene (13 200-11 400 cal. BP). Data also imply that this period was a major warm phase in southeastern Patagonia and was approximately contemporaneous with the Younger Dryas chronozone in the Northern Hemisphere (12 700-11 500 cal. BP). After 8650 cal. BP a major drought may have caused the lowest lake level of the record. Since 7300 cal. BP, the lake level rose and was variable until the 'Little Ice Age', which was the dominant humid period after 8650 cal. BP.

Stable oxygen isotope and Mg/Ca ratios of planktonic foraminifera from the Holocene and the Last Glacial Maximum

We measured oxygen isotopes and Mg/Ca ratios in the surface-dwelling planktonic foraminifer Globigerinoides ruber (white s.s.) and the thermocline dweller Pulleniatina obliquiloculata to investigate upper ocean spatial variability in the Indo-Pacific Warm Pool (IPWP). We focused on three critical time intervals: the Last Glacial Maximum (LGM; 18-21.5 ka), the early Holocene (8-9 ka), and the late Holocene (0-2 ka). Our records from 24 stations in the South China Sea, Timor Sea, Indonesian seas, and western Pacific indicate overall dry and cool conditions in the IPWP during the LGM with a low thermal gradient between surface and thermocline waters. During the early Holocene, sea surface temperatures increased by ~3°C over the entire region, indicating intensification of the IPWP. However, in the eastern Indian Ocean (Timor Sea), the thermocline gradually shoaled from the LGM to early Holocene, reflecting intensification of the subsurface Indonesian Throughflow (ITF). Increased surface salinity in the South China Sea during the Holocene appears related to northward displacement of the monsoonal rain belt over the Asian continent together with enhanced influx of saltier Pacific surface water through the Luzon Strait and freshwater export through the Java Sea. Opening of the freshwater portal through the Java Sea in the early Holocene led to a change in the vertical structure of the ITF from surface- to thermocline-dominated flow and to substantial freshening of Timor Sea thermocline waters.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (d13Catm ), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of d13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 years BP). The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4 permil shift to heavier values between the mean d13Catm level in the Penultimate (~ 140 000 years BP) and Last Glacial Maximum (~ 22 000 years BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 years, but with different phasing and magnitudes. Furthermore, a 5000 years lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS 5.5 (120 000 years BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Constraints on the magnitude and patterns of ocean cooling at the Last Glacial Maximum - Supplementary material

Observation-based reconstructions of sea surface temperature from relatively stable periods in the past, such as the Last Glacial Maximum, represent an important means of constraining climate sensitivity and evaluating model simulations. The first quantitative global reconstruction of sea surface temperatures during the Last Glacial Maximum was developed by the Climate Long-Range Investigation, Mapping and Prediction (CLIMAP) project in the 1970s and 1980s. Since that time, several shortcomings of that earlier effort have become apparent. Here we present an updated synthesis of sea surface temperatures during the Last Glacial Maximum, rigorously defined as the period between 23 and 19 thousand years before present, from the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project. We integrate microfossil and geochemical reconstructions of surface temperatures and include assessments of the reliability of individual records. Our reconstruction reveals the presence of large longitudinal gradients in sea surface temperature in all of the ocean basins, in contrast to the simulations of the Last Glacial Maximum climate available at present.

Plant-wax, alkenone and major element analyses of ODP Site 659 for the Last Glacial cycle and associated coretop data

We present a hydrologic reconstruction of the Sahara-Sahel transition, covering the complete last glacial cycle (130 ka), based on a combination of plant-wax-specific hydrogen (dD) and carbon isotopes (d13C). The dD and d13C signatures of long-chain n-alkanes from ODP Site 659 off NW Africa reveal a significant anti-correlation. Complementary to published pollen data, we infer that this plant-wax signal reflects sensitive responses of the vegetation cover to precipitation changes in the Sahel region, as well as varying contributions from biomes north of the Sahara (C3 domain) by North-East Trade Winds (NETW). During arid phases, especially the northern parts of the Sahel likely experienced crucial water stress, which resulted in a pronounced contraction of the vegetation cover, thus reducing the amount of C4 plant waxes from the region. The increase in NETW strength during dry periods further promoted a more pronounced C3-plant-wax signal derived from the North African C3 plant domain. During humid periods, the C4-dominated Sahelian environments spread northward into the Saharan realm, in association with lower NETW inputs of C3 plant waxes. Arid-humid cycles deduced from plant-wax dD are in accordance with concomitant changes in weathering intensity reflected in varying major element distributions. Environmental shifts are generally linked to periods with large fluctuations in Northern Hemisphere summer insolation. During Marine Isotope Stages 2 and 3, when insolation variability was low, coupling of the hydrologic regime to alkenone-based estimates of NE Atlantic sea-surface temperatures becomes apparent.