Late glacial and deglacial age model and biomarker analyses of sediment core MSM05/5_712-2 from the western continental slope of Svalbard

The transition from last glacial to deglacial and subsequently to modern interglacial climate conditions was accompanied by abrupt shifts in the palaeoceanographic setting in the subpolar North Atlantic. Knowledge about the role that sea ice coverage played during these rapid climate reversals is limited since most marine sediment cores from the higher latitudes provide only a coarse temporal resolution and often poorly preserved microfossils. Here we present a highly resolved reconstruction of sea ice conditions that characterised the eastern Fram Strait - a key area for water mass exchange between the Arctic Ocean and the North Atlantic - for the past 30 ka BP. This reconstruction is based on the distribution of the sea ice biomarker IP25 and phytoplankton derived biomarkers in a sediment core from the continental slope of western Svalbard. During the late glacial (30 ka to 19 ka BP), recurrent advances and retreats of sea ice characterised the study area and point to a hitherto less considered oceanic (and/or atmospheric) variability. A long-lasting perennial sea ice coverage in eastern Fram Strait persisted only at the very end of the Last Glacial Maximum (i.e. from 19.2 to 17.6 ka BP) and was abruptly reduced at the onset of Heinrich Event 1 - coincident with or possibly even inducing the collapse of the Atlantic Meridional Overturning Circulation (AMOC). Further maximum sea ice conditions prevailed during the Younger Dryas cooling event and support the assumption of an AMOC reduction due to increased formation and export of Arctic sea ice through Fram Strait. A significant retreat of sea ice and sea surface warming are observed for the Early Holocene.

Water mass evolution of the Greenland Sea since late glacial times

Four sediment cores from the central and northern Greenland Sea basin, a crucial area for the renewal of North Atlantic deep water, were analyzed for planktic foraminiferal fauna, planktic and benthic stable oxygen and carbon iso- topes as well as ice-rafted debris to reconstruct the environ- mental variability in the last 23 kyr. During the Last Glacial Maximum, the Greenland Sea was dominated by cold and sea-ice bearing surface water masses. Meltwater discharges from the surrounding ice sheets affected the area during the deglaciation, influencing the water mass circulation. During the Younger Dryas interval the last major freshwater event occurred in the region. The onset of the Holocene interglacial was marked by an increase in the advection of Atlantic Wa- ter and a rise in sea surface temperatures (SST). Although the thermal maximum was not reached simultaneously across the basin, benthic isotope data indicate that the rate of overturn- ing circulation reached a maximum in the central Greenland Sea around 7ka. After 6-5ka a SST cooling and increas- ing sea-ice cover is noted. Conditions during this so-called "Neoglacial" cooling, however, changed after 3 ka, probably due to enhanced sea-ice expansion, which limited the deep convection. As a result, a well stratified upper water column amplified the warming of the subsurface waters in the central Greenland Sea, which were fed by increased inflow of At- lantic Water from the eastern Nordic Seas. Our data reveal that the Holocene oceanographic conditions in the Green- land Sea did not develop uniformly. These variations were a response to a complex interplay between the Atlantic and Polar water masses, the rate of sea-ice formation and melting and its effect on vertical convection intensity during times of Northern Hemisphere insolation changes.

Late Glacial and Holocene data from sediment cores MD99-2236 and MD99-2256

The Laurentide Ice Sheet (LIS) was a large, dynamic ice sheet in the early Holocene. The glacial events through Hudson Strait leading to its eventual demise are recorded in the well-dated Labrador shelf core, MD99-2236 from the Cartwright Saddle. We develop a detailed history of the timing of ice-sheet discharge events from the Hudson Strait outlet of the LIS during the Holocene using high-resolution detrital carbonate, ice rafted detritus (IRD), d18O, and sediment color data. Eight detrital carbonate peaks (DCPs) associated with IRD peaks and light oxygen isotope events punctuate the MD99-2236 record between 11.5 and 8.0 ka. We use the stratigraphy of the DCPs developed from MD99-2236 to select the appropriate DeltaR to calibrate the ages of recorded glacial events in Hudson Bay and Hudson Strait such that they match the DCPs in MD99-2236. We associate the eight DCPs with H0, Gold Cove advance, Noble Inlet advance, initial retreat of the Hudson Strait ice stream (HSIS) from Hudson Strait, opening of the Tyrrell Sea, and drainage of glacial lakes Agassiz and Ojibway. The opening of Foxe Channel and retreat of glacial ice from Foxe Basin are represented by a shoulder in the carbonate data. DeltaR of 350 years applied to the radiocarbon ages constraining glacial events H0 through the opening of the Tyrell Sea provided the best match with the MD99-2236 DCPs; DeltaR values and ages from the literature are used for the younger events. A very close age match was achieved between the 8.2 ka cold event in the Greenland ice cores, DCP7 (8.15 ka BP), and the drainage of glacial lakes Agassiz and Ojibway. Our stratigraphic comparison between the DCPs in MD99-2236 and the calibrated ages of Hudson Strait/Bay deglacial events shows that the retreat of the HSIS, the opening of the Tyrell Sea, and the catastrophic drainage of glacial lakes Agassiz and Ojibway at 8.2 ka are separate events that have been combined in previous estimates of the timing of the 8.2 ka event from marine records. SW Iceland shelf core MD99-2256 documents freshwater entrainment into the subpolar gyre from the Hudson Strait outlet via the Labrador, North Atlantic, and Irminger currents. The timing of freshwater release from the LIS Hudson Strait outlet in MD99-2236 matches evidence for freshwater forcing and LIS icebergs carrying foreign minerals to the SW Iceland shelf between 11.5 and 8.2 ka. The congruency of these records supports the conclusion of the entrainment of freshwater from the retreat of the LIS through Hudson Strait into the subpolar gyre and provides specific time periods when pulses of LIS freshwater were present to influence climate.

(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.

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.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Last Glacial Maximum carbon isotopic composition

We present Holocene and last glacial maximum (LGM) oxygen and carbon isotope measurements on Planulina wuellerstorfi in six southeast Pacific cores. Sedimentation rates are low in this part of the ocean, and measurements were made on individual foraminiferal shells in order to identify the Holocene and glacial individuals on the basis of their extreme d18O. The new d13C data were combined with previous P. wuellerstorfi data for interpretation of global thermohaline circulation. Data from the Southern Ocean were examined closely for regional coherency and a few anomalous d13C values suspected of having productivity overprint were removed. The resulting global d13C distributions and gradients indicate that the deep water circulation was similar during the Holocene and LGM. This interpretation brings d13C data to a better agreement with Cd/Ca data and marks a sharp contrast with a widely held view based on d13C measurements that the glacial Southern Ocean was the terminus of the thermohaline circulation. The proposed presence of glacial North Atlantic Deep Water does not necessarily contradict the postulated presence of Glacial North Atlantic Intermediate Water.

Global compilation of last glacial maximum oxygen isotopic ratios for planktonic and benthic foraminifera

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 yr). By using a model that is able to explicitly simulate the sensor (d18O), 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 d18O 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 d18O 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 d18O 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 trong mean annual cooling 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 d18O 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 d18Osw anomaly for the north Indian Ocean that contrasts with a large negative d18Osw anomaly in the China Sea between the LGM and 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 d18Osw between the LGM and present are not spatially homogenous. 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 d18O data. Nonetheless, our data-model comparison support a heterogeneous cooling of few degrees (2-4°C) in the LGM Ocean.

Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates

Mineral dust aerosols play a major role in present and past climates. To date, we rely on climate models for estimates of dust fluxes to calculate the impact of airborne micronutrients on biogeochemical cycles. Here we provide a new global dust flux data set for Holocene and Last Glacial Maximum (LGM) conditions based on observational data. A comparison with dust flux simulations highlights regional differences between observations and models. By forcing a biogeochemical model with our new data set and using this model's results to guide a millennial-scale Earth System Model simulation, we calculate the impact of enhanced glacial oceanic iron deposition on the LGM-Holocene carbon cycle. On centennial timescales, the higher LGM dust deposition results in a weak reduction of <10?ppm in atmospheric CO2 due to enhanced efficiency of the biological pump. This is followed by a further ~10?ppm reduction over millennial timescales due to greater carbon burial and carbonate compensation.

Gridded climate data from 5 Global Climate Models (GCM) of the Last Glacial Maximum (LGM) downscaled to 30 arc seconds for Europe, with links to NetCDF files

Studies on the impact of historical, current and future global change require very high-resolution climate data (less or equal 1km) as a basis for modelled responses, meaning that data from digital climate models generally require substantial rescaling. Another shortcoming of available datasets on past climate is that the effects of sea level rise and fall are not considered. Without such information, the study of glacial refugia or early Holocene plant and animal migration are incomplete if not impossible. Sea level at the last glacial maximum (LGM) was approximately 125m lower, creating substantial additional terrestrial area for which no current baseline data exist. Here, we introduce the development of a novel, gridded climate dataset for LGM that is both very high resolution (1km) and extends to the LGM sea and land mask. We developed two methods to extend current terrestrial precipitation and temperature data to areas between the current and LGM coastlines. The absolute interpolation error is less than 1°C and 0.5 °C for 98.9% and 87.8% of all pixels for the first two 1 arc degree distance zones. We use the change factor method with these newly assembled baseline data to downscale five global circulation models of LGM climate to a resolution of 1km for Europe. As additional variables we calculate 19 'bioclimatic' variables, which are often used in climate change impact studies on biological diversity. The new LGM climate maps are well suited for analysing refugia and migration during Holocene warming following the LGM.

Paleoenvironmental changes in the western Mediterranean since the last glacial maximum: High resolution multiproxy record from the Algero-Balearic basin

The present study uses a multiproxy approach in order to further understand the evolution of climate responses in the western Mediterranean as of the Last Glacial Maximum. Sediments from ODP Site 975 in the Algero-Balearic basin have been analysed at high resolution, both geochemically andmineralogicallly. The resulting data have been used as proxies to establish a sedimentary regime, primary marine productivity, the preservation of the proxies and oxygen conditions. Fluctuations in detrital element concentrations were mainly the consequence of wet/arid oscillations. Productivity has been established using Ba excess, according to which marine productivity appears to have been greatest during cold events Heinrich 1 and Younger Dryas. The S1 time interval was not as marked by increases in productivity as was the eastern Mediterranean. In contrast, the S1 interval was first characterized by a decreasing trend and then by a fall in productivity after the 8.2 ky BP dry-cold event. Since then productivity has remained low. Here we report that there was an important redox event in this basin, probably a consequence of the major oceanographic circulation change occurring in the western Mediterranean at 7.7 ky BP. This circulation change led to reventilation as well as to diagenetic remobilization of redox-sensitive elements and organic matter oxidation. Comparisons between our paleoceanographic reconstruction for this basin and those regarding other Mediterranean basins support the hypothesis that across the Mediterranean there were different types of responses to climate forcing mechanism. The Algero-Balearic basin is likely to be a key area for further understanding of the relationships between the North Atlantic and the eastern Mediterranean basins.