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.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Isotopic data of atmospheric methane from EDML and Vostok ice cores over a full glacial cycle

The response of natural CH4 sources to climate changes will be an important factor to consider as concentrations of this potent greenhouse gas continue to increase. Polar ice cores provide the means to assess this sensitivity in the past and have shown a close connection between CH4 levels and northern hemisphere temperature variability over the last glacial cycle. However, the contribution of the various CH4 sources and sinks to these changes is still a matter of debate. Contemporaneous stable CH4 isotope records in ice cores provide additional boundary conditions for assessing changes in the CH4 sources and sinks. Here we present new ice core CH4 isotope data covering the last 160,000 years, showing a clear decoupling between CH4 loading and carbon isotopic variations over most of the record. We suggest that d13CH4 variations were not dominated by a change in the source mix but rather by climate- and CO2-related ecosystem control on the isotopic composition of the methane precursor material, especially in seasonally inundated wetlands in the tropics. In contrast, relatively stable d13CH4 intervals occurred during large CH4 loading changes concurrently with past climate changes implying that most CH4 sources (most notably tropical wetlands) responded simultaneously.

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.

Glacial cold-water coral: ages, isotope concentrations and ratios

A set of 40 Uranium-series datings obtained on the reef-forming scleractinian cold-water corals Lophelia pertusa and Madrepora oculata revealed that during the past 400 kyr their occurrence in the Gulf of Cádiz (GoC) was almost exclusively restricted to glacial periods. This result strengthens the outcomes of former studies that coral growth in the temperate NE Atlantic encompassing the French, Iberian and Moroccan margins dominated during glacial periods, whereas in the higher latitudes (Irish and Norwegian margins) extended coral growth prevailed during interglacial periods. Thus it appears that the biogeographical limits for sustained cold-water coral growth along the NE Atlantic margin are strongly related to climate change. By focussing on the last glacial-interglacial cycle, this study shows that palaeo-productivity was increased during the last glacial. This was likely driven by the fertilisation effect of an increased input of aeolian dust and locally intensified upwelling. After the Younger Dryas cold event, the input of aeolian dust and productivity significantly decreased concurrent with an increase in water temperatures in the GoC. This primarily resulted in reduced food availability and caused a widespread demise of the formerly thriving coral ecosystems. Moreover, these climate induced changes most likely caused a latitudinal shift of areas withoptimum coral growth conditions towards the northern NE Atlantic where more suitable environmental conditions established with the onset of the Holocene.

Last glacial maximum sediment record of the southern Weddell Sea shelf

Abstract: The history of grounded ice-sheet extent on the southern Weddell Sea shelf during the Last Glacial Maximum (LGM) and the timing of post-LGM ice-sheet retreat are poorly constrained. Several glaciological models reconstructed widespread grounding and major thickening of the Antarctic Ice Sheet in the Weddell Sea sector at the LGM. In contrast, recently published onshore data and modelling results concluded only very limited LGM-thickening of glaciers and ice streams feeding into the modern Filchner and Ronne ice shelves. These studies concluded that during the LGM ice shelves rather than grounded ice covered the Filchner and Ronne troughs, two deep palaeo-ice stream troughs eroded into the southern Weddell Sea shelf. Here we review previously published and unpublished marine geophysical and geological data from the southern Weddell Sea shelf. The stratigraphy and geometry of reflectors in acoustic sub-bottom profiles are similar to those from other West Antarctic palaeo-ice stream troughs, where grounded ice had advanced to the shelf break at the LGM. Numerous cores from the southern Weddell Sea shelf recovered sequences with properties typical for subglacially deposited tills or subglacially compacted sediments. These data sets give evidence that grounded ice had advanced across the shelf during the past, thereby grounding in even the deepest parts of the Filchner and Ronne troughs. Radiocarbon dates from glaciomarine sediments overlying the subglacial deposits are limited, but indicate that the ice grounding occurred at the LGM and that ice retreat started before ~15.1 corrected 14C kyrs before present (BP) on the outer shelf and before ~7.7 corrected 14C kyrs BP on the inner shelf, which is broadly synchronous with ice retreat in other Antarctic sectors. The apparent mismatch between the ice-sheet reconstructions from marine and terrestrial data can be attributed to ice streams with very low surface profiles (similar to those of "ice plains") that had advanced through Filchner Trough and Ronne Trough at the LGM. Considering the global sea-level lowstand of ~130 metres below present, a low surface slope of the expanded LGM-ice sheet in the southern Weddell Sea can reconcile grounding-line advance to the shelf break with limited thickening of glaciers and ice streams in the hinterland. This scenario implies that ice-sheet growth in the Weddell Sea sector during the LGM and ice-sheet drawdown throughout the last deglaciation could only have made minor contributions to the major global sea-level fluctuations during these times.

Sedimentology of the northwestern Weddel Sea continental shelf and slope

Geologie cores on two profiles oriented normaly to the continental shelf and slope, have been investigated to reconstruct the Quaternary sedimentary history of the southeast continental border of South Orkney (NW Weddell Sea). The sediments were described macroscopically and their fabric investigated by use of X-radiographs. Laboratory work comprised detailed grain-size analysis, determination of the watercontent, carbonate, organic carbon and sand fraction.composition. Stable oxygen and carbon isotopes have been measured On planktonic foraminifera. Palaeomagnetism, analysis of 230Th-content and detailed comparison of the lithlogic Parameters with the oxygen isotope stages (Martinson curve) were used for stratigraphic classification of the sediments. The sediment cores from the continental slope comprise a maximum age of 300,000 years B. P.. Bottom currents, ice rafting and biogenic input are the main sources of sediment. Based on lithologic parameters a distinction between glacial and interglacial facies is possible. Silty clays without microfossils and few bioturbation characterise the sediments of the glacial facies. Only small amounts of icerafted debris can be recognized. This type of sediment was accumulated during times of lower sea-level and drastically reduced rate of bottom water production. Based on grain-size distribution, bottom current velocities of 0.01 cmls were calculated. Thick sea-ice coverage reduced biogenic production in the surface water, and as consequence benthic communities were depleted. Because of the reduced benthic life, sediments are only slithly bioturbated. At the beginning of the interglacial Stage, the sea-level rised rapidly, and calving rate of icebergs, combined with input of ice-rafted material, increased considerably. Sediments of this transition facies are silty cliiys with a high proportion of coarse ice-rafted debris, but without microfossils. With the onset of bottom water production in connection with shelf ice water, sediments of interglacial facies were formed. They consist of silty clays to clayey silts with considerable content of sand and gravel. Sediments are strongly bioturbated. Based On the sediment caracteristics, current velocities of the bottom water were calculated to be of 0.96 cmls for interglacials. At the southern slope of a NW/SE-striking ridge, bottom water current is channelized, resulting in a drastic increase of current velocities. Current velocities up to 7.5 cm/s lead to formation of residual sediments. While the continental slope has predominantly fine sediments, the South Orkney shelf are mainly sandy silts and silty sands with a high proportion of gravel. These sediments were formed dominantly by ice-rafting during Brunhes- and Matuyama-Epoch. Currents removed the fine fraction of the sediments. Based on microfossil contents it was not possible to differentiate sediments from glacial to interglacial. In the upper Parts of the cores graded sequences truncated by erosion were observed. These sequences were formed during Brunhes-Epoch by strong currents with velocities decreasing periodically from about 7.5 cm/s to about 1 cm/s. Sediments with a high proportion of siliceous microfossils but barren of foraminifera compose the lower part of the shelf cores. These sediments have formed during the warmer Matuyama-Epoch.

(Table S1) Chronostratigraphic constrains for sediment core AND1-1B

Thirty years after oxygen isotope records from microfossils deposited in ocean sediments confirmed the hypothesis that variations in the Earth's orbital geometry control the ice ages (Hays et al., 1976, doi:10.1126/science.194.4270.1121), fundamental questions remain over the response of the Antarctic ice sheets to orbital cycles (Raymo and Huybers, 2008, doi:10.1038/nature06589). Furthermore, an understanding of the behaviour of the marine-based West Antarctic ice sheet (WAIS) during the 'warmer-than-present' early-Pliocene epoch (~5-3 Myr ago) is needed to better constrain the possible range of ice-sheet behaviour in the context of future global warming (Solomon et al., 2007). Here we present a marine glacial record from the upper 600 m of the AND-1B sediment core recovered from beneath the northwest part of the Ross ice shelf by the ANDRILL programme and demonstrate well-dated, ~40-kyr cyclic variations in ice-sheet extent linked to cycles in insolation influenced by changes in the Earth's axial tilt (obliquity) during the Pliocene. Our data provide direct evidence for orbitally induced oscillations in the WAIS, which periodically collapsed, resulting in a switch from grounded ice, or ice shelves, to open waters in the Ross embayment when planetary temperatures were up to ~3° C warmer than today ( Kim and Crowley, 2000, doi:10.1029/1999PA000459) and atmospheric CO2 concentration was as high as ~400 p.p.m.v. (van der Burgh et al., 1993, doi:10.1126/science.260.5115.1788, Raymo et al., 1996, doi:10.1016/0377-8398(95)00048-8). The evidence is consistent with a new ice-sheet/ice-shelf model (Pollard and DeConto, 2009, doi:10.1038/nature07809) that simulates fluctuations in Antarctic ice volume of up to +7 m in equivalent sea level associated with the loss of the WAIS and up to +3 m in equivalent sea level from the East Antarctic ice sheet, in response to ocean-induced melting paced by obliquity. During interglacial times, diatomaceous sediments indicate high surface-water productivity, minimal summer sea ice and air temperatures above freezing, suggesting an additional influence of surface melt (Huybers, 2006, doi:10.1126/science.1125249) under conditions of elevated CO2.

Sedimentary record of the western Amundsen Sea Embayment after the last glacial maximum

Reliable dating of glaciomarine sediments deposited on the Antarctic shelf since the Last Glacial Maximum (LGM) is very challenging because of the general absence of calcareous (micro-) fossils and the recycling of fossil organic matter. As a consequence, radiocarbon (14C) ages of the acid-insoluble organic fraction (AIO) of the sediments bear uncertainties that are very difficult to quantify. In this paper we present the results of three different chronostratigraphic methods to date a sedimentary unit consisting of diatomaceous ooze and diatomaceous mud that was deposited following the last deglaciation at five core sites on the inner shelf in the western Amundsen Sea (West Antarctica). In three cores conventional 14C dating of the AIO in bulk sediment samples yielded age reversals down-core, but at all sites the AIO 14C ages obtained from diatomaceous ooze within the diatom-rich unit yielded similar uncorrected 14C ages ranging from 13,517±56 to 11,543±47 years before present (yr BP). Correction of these ages by subtracting the core-top ages, which are assumed to reflect present-day deposition (as indicated by 21044 Pb dating of the sediment surface at one core site), yielded ages between ca. 10,500 and 8,400 calibrated years before present (cal yr BP). Correction of the AIO ages of the diatomaceous ooze by only subtracting the marine reservoir effect (MRE) of 1,300 years indicated deposition of the diatom-rich sediments between 14,100 and 11,900 cal yr BP. Most of these ages are consistent with age constraints between 13.0 and 8.0 ka BP for the diatom-rich unit, which we obtained by correlating the relative palaeomagnetic intensity (RPI) records of three of the sediment cores with global and regional reference curves for palaeomagnetic intensity. As a third dating technique we applied conventional 53 radiocarbon dating of the AIO included in acid-cleaned diatom hard parts that were extracted from the diatomaceous ooze. This method yielded uncorrected 14C ages of only 5,111±38 and 5,106±38 yr BP, respectively. We reject these young ages, because they are likely to be overprinted by the adsorption of modern atmospheric carbon dioxide onto the surfaces of the extracted diatom hard parts prior to sample graphitisation and combustion for 14C dating. The deposition of the diatom-rich unit in the western Amundsen Sea suggests deglaciation of the inner shelf before ca. 13 ka BP. The deposition of diatomaceous oozes on other parts of the Antarctic shelf around the same time, however, seems to be coincidental rather than directly related.

Sea-surface temperature and sea ice distribution of the Southern Ocean at the EPILOG Last Glacial Maximum

Based on the quantitative study of diatoms and radiolarians, summer sea-surface temperature (SSST) and sea ice distribution were estimated from 122 sediment core localities in the Atlantic, Indian and Pacific sectors of the Southern Ocean to reconstruct the last glacial environment at the EPILOG (19.5-16.0 ka or 23 000-19 000 cal yr. B.P.) time-slice. The statistical methods applied include the Imbrie and Kipp Method, the Modern Analog Technique and the General Additive Model. Summer SSTs reveal greater surface-water cooling than reconstructed by CLIMAP (Geol. Soc. Am. Map Chart. Ser. MC-36 (1981) 1), reaching a maximum (4-5 °C) in the present Subantarctic Zone of the Atlantic and Indian sector. The reconstruction of maximum winter sea ice (WSI) extent is in accordance with CLIMAP, showing an expansion of the WSI field by around 100% compared to the present. Although only limited information is available, the data clearly show that CLIMAP strongly overestimated the glacial summer sea ice extent. As a result of the northward expansion of Antarctic cold waters by 5-10° in latitude and a relatively small displacement of the Subtropical Front, thermal gradients were steepened during the last glacial in the northern zone of the Southern Ocean. Such reconstruction may, however, be inapposite for the Pacific sector. The few data available indicate reduced cooling in the southern Pacific and give suggestion for a non-uniform cooling of the glacial Southern Ocean.

Age of bottom sediments from the Indian Ocean

It has been found that oxygen-isotope and paleotemperature curves based on types of planktonic foraminiferal thanatocenoses in three sediment cores, from the tropical, southern temperate, and southern glacial zones of the Indian Ocean can be readily correlated with each other. The sediment cores revealed three epochs of cold climate during the past 700 ky; these are probably connect with worldwide epochs of cooling during Pleistocene that led to advance of ice sheets during continental glaciations in the northern hemisphere.

Glacial landscape evolution in the Uummannaq region, West Greenland

Funded by Department of Geography (Durham University) Department of Geography and the Environment (University of Aberdeen) Royal Geographical Society-IBG Carnegie Trust for the Universities of Scotland

Dissolution index of Globigerina bulloides in recent and Last Glacial Maximum sediments

The modern Atlantic Ocean, dominated by the interactions of North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW), plays a key role in redistributing heat from the Southern to the Northern Hemisphere. In order to reconstruct the evolution of the relative importance of these two water masses, the NADW/AABW transition, reflected by the calcite lysocline, was investigated by the Globigerina bulloides dissolution index (BDX?). The depth level of the Late Glacial Maximum (LGM) calcite lysocline was elevated by several hundred metres, indicating a more corrosive water mass present at modern NADW level. Overall, the small range of BDX? data and the gradual decrease in preservation below the calcite lysocline point to a less stratified Atlantic Ocean during the LGM. Similar preservation patterns in the West and East Atlantic demonstrate that the modern west-east asymmetry did not exist due to an expansion of southern deep waters compensating for the decrease in NADW formation.