Marine ice-rafted debris records and age determinations along the northern Eurasian margin

Ice-rafted debris (IRD) (>2 mm), input in eight sediment cores along the Eurasian continental margin (Arctic Ocean), have been studied over the last two glacial/interglacial cycles. Together with the revised chronologies and new micropaleontological data of two cores from the northern Barents Sea (PS2138) and northeastern Kara Sea (PS2741) spanning Marine Isotope Stages (MIS) 6 to 1, the IRD data give new insights into the glacial history of northern Eurasian ice-sheets over the last 150 ka. The chronologies of the cores are based on stable isotope records, AMS 14C datings, paleomagnetic and biostratigraphic data. Extensive episodes of northern Barents Sea ice-sheet growth, probably to the shelf edge, occurred during the late Weichselian (MIS 2) and the Saalian (MIS 6). Major IRD discharge at the MIS 4/3-transition hints to another severe glaciation, probably onto the outer shelf, during MIS 4. IRD-based instabilities of the marine-based ice margin along the northern Barents Sea between MIS 4 and 2 are similar in timing with North Atlantic Heinrich events and Nordic Seas IRD events, suggesting similar atmospheric cooling over a broad region or linkage of ice-sheet fluctuations through small sea-level events. In the relatively low-precipitation areas of eastern Eurasia, IRD peak values during Termination II and MIS 4/3-transition suggest a Kara Sea ice-sheet advance onto the outer shelf, probably to the shelf edge, during glacial MIS 6 and 4. This suggests that during the initial cooling following the interglacials MIS 5, and possibly MIS 7, the combined effect of sustained inflow of Atlantic water into the Arctic Ocean and penetration of moisture-bearing cyclones into easterly direction supported major ice build-up during Saalian (MIS 6) and Mid-Weichselian (MIS 4) glaciation. IRD peak values in MIS 5 indicate at least two advances of the Severnaya Semlya ice-sheet to the coast line during the Early Weichselian. In contrast, a distinct Kara Sea ice advance during the Late Weichselian (MIS 2) is not documented by the IRD records along the northeastern Kara Sea margin.

Colour reflectance and core descriptions from sediment cores from the Weddell Sea

The timing of the last maximum extent of the Antarctic ice sheets relative to those in the Northern Hemisphere remains poorly understood. We develop a chronology for the Weddell Sea sector of the East Antarctic Ice Sheet that, combined with ages from other Antarctic ice-sheet sectors, indicates that the advance to and retreat from their maximum extent was within dating uncertainties synchronous with most sectors of Northern Hemisphere ice sheets. Surface climate forcing of Antarctic mass balance would probably cause an opposite response, whereby a warming climate would increase accumulation but not surface melting. Our new data support teleconnections involving sea-level forcing from Northern Hemisphere ice sheets and changes in North Atlantic deep-water formation and attendant heat flux to Antarctic grounding lines to synchronize the hemispheric ice sheets.

Sedimentology of the southern Bellingshausen Sea

A major trough ('Belgica Trough') eroded by a palaeo-ice stream crosses the continental shelf of the southern Bellingshausen Sea (West Antarctica) and is associated with a trough mouth fan ('Belgica TMF') on the adjacent continental slope. Previous marine geophysical and geological studies investigated the bathymetry and geomorphology of Belgica Trough and Belgica TMF, erosional and depositional processes associated with bedform formation, and the temporal and spatial changes in clay mineral provenance of subglacial and glaciomarine sediments. Here, we present multi-proxy data from sediment cores recovered from the shelf and uppermost slope in the southern Bellingshausen Sea and reconstruct the ice-sheet history since the last glacial maximum (LGM) in this poorly studied area of West Antarctica. We combined new data (physical properties, sedimentary structures, geochemical and grain-size data) with published data (shear strength, clay mineral assemblages) to refine a previous facies classification for the sediments. The multi-proxy approach allowed us to distinguish four main facies types and to assign them to the following depositional settings: 1) subglacial, 2) proximal grounding-line, 3) distal sub-ice shelf/subsea ice, and 4) seasonal open-marine. In the seasonal open-marine facies we found evidence for episodic current-induced winnowing of near-seabed sediments on the middle to outer shelf and at the uppermost slope during the late Holocene. In addition, we obtained data on excess 210Pb activity at three core sites and 44 AMS 14C dates from the acid-insoluble fraction of organic matter (AIO) and calcareous (micro-)fossils, respectively, at 12 sites. These chronological data enabled us to reconstruct, for the first time, the timing of the last advance and retreat of the West Antarctic Ice Sheet (WAIS) and the Antarctic Peninsula Ice Sheet (APIS) in the southern Bellingshausen Sea. We used the down-core variability in sediment provenance inferred from clay mineral changes to identify the most reliable AIO 14C ages for ice-sheet retreat. The palaeo-ice stream advanced through Belgica Trough after ~36.0 corrected 14C ka before present (B.P.). It retreated from the outer shelf at ~25.5 ka B.P., the middle shelf at ~19.8 ka B.P., the inner shelf in Eltanin Bay at ~12.3 ka B.P., and the inner shelf in Ronne Entrance at ~6.3 ka B.P.. The retreat of the WAIS and APIS occurred slowly and stepwise, and may still be in progress. This dynamical ice-sheet behaviour has to be taken into account for the interpretation of recent and the prediction of future mass-balance changes in the study area. The glacial history of the southern Bellingshausen Sea is unique when compared to other regions in West Antarctica, but some open questions regarding its chronology need to be addressed by future work.

Clay mineralogy and shear strength of subglacial and glaciomarine sediments in the southern Bellingshausen Sea

The Belgica Trough and the adjacent Belgica Trough Mouth Fan in the southern Bellingshausen Sea (Pacific sector of the Southern Ocean) mark the location of a major outlet for the West Antarctic Ice Sheet during the Late Quaternary. The drainage basin of an ice stream that advanced through Belgica Trough across the shelf during the last glacial period comprised an area exceeding 200,000 km**2 in the West Antarctic hinterland. Previous studies, mainly based on marine-geophysical data from the continental shelf and slope, focused on the bathymetry and seafloor bedforms, and the reconstruction of associated depositional processes and ice- drainage patterns. In contrast, there was only sparse information from seabed sediments recovered by coring. In this paper, we present lithological and clay mineralogical data of 21 sediment cores collected from the shelf and slope of the southern Bellingshausen Sea. Most cores recovered three lithological units, which can be attributed to facies types deposited under glacial, transitional and seasonally open-marine conditions. The clay mineral assemblages document coinciding changes in provenance. The relationship between the clay mineral assemblages in the subglacial and proglacial sediments on the shelf and the glacial diamictons on the slope confirms that a grounded ice stream advanced through Belgica Trough to the shelf break during the past, thereby depositing detritus eroded in the West Antarctic hinterland as soft till on the shelf and as glaciogenic debris flows on the slope. The thinness of the transitional and seasonally open-marine sediments in the cores suggests that this ice advance occurred during the last glacial period. Clay mineralogical, acoustic sub-bottom and seismic data furthermore demonstrate that the palaeo-ice stream probably reworked old sedimentary strata, including older tills, on the shelf and incorporated this debris into its till bed. The geographical heterogeneity of the clay mineral assemblages in the sub- and proglacial diamictons and gravelly deposits indicates that they were eroded from underlying sedimentary strata of different ages. These strata may have been deposited during either different phases of the last glacial period or different glacial and interglacial periods. Additionally, the clay mineralogical heterogeneity of the soft tills recovered on the shelf suggests that the drainage area of the palaeo-ice stream flowing through Belgica Trough changed through time.

Grain size composition of sediment cores from the Weddell Sea, Antarctica

Sedimentary processes in the southeastern Weddell Sea are influenced by glacial-interglacial ice-shelf dynamics and the cyclonic circulation of the Weddell Gyre, which affects all water masses down to the sea floor. Significantly increased sedimentation rates occur during glacial stages, when ice sheets advance to the shelf edge and trigger gravitational sediment transport to the deep sea. Downslope transport on the Crary Fan and off Dronning Maud and Coats Land is channelized into three huge channel systems, which originate on the eastern-, the central and the western Crary Fan. They gradually turn from a northerly direction eastward until they follow a course parallel to the continental slope. All channels show strongly asymmetric cross sections with well-developed levees on their northwestern sides, forming wedge-shaped sediment bodies. They level off very gently. Levees on the southeastern sides are small, if present at all. This characteristic morphology likely results from the process of combined turbidite-contourite deposition. Strong thermohaline currents of the Weddell Gyre entrain particles from turbidity-current suspensions, which flow down the channels, and carry them westward out of the channel where they settle on a surface gently dipping away from the channel. These sediments are intercalated with overbank deposits of high-energy and high-volume turbidity currents, which preferentially flood the left of the channels (looking downchannel) as a result of Coriolis force. In the distal setting of the easternmost channel-levee complex, where thermohaline currents are directed northeastward as a result of a recirculation of water masses from the Enderby Basin, the setting and the internal structures of a wedge-shaped sediment body indicate a contourite drift rather than a channel levee. Dating of the sediments reveals that the levees in their present form started to develop with a late Miocene cooling event, which caused an expansion of the East Antarctic Ice Sheet and an invigoration of thermohaline current activity.

(Table 1) Age determination and isotopic signature of last glacial Patagonian caly/silt

Ice cores provide a record of changes in dust flux to Antarctica, which is thought to reflect changes in atmospheric circulation and environmental conditions in dust source areas (Forster et al., 2007; Diekmann et al. 2000, doi:10.1016/S0031-0182(00)00138-3; Winckler et al., 2008, doi:10.1126/science.1150595; Reader et al., 1999, doi:10.1029/1999JD900033; Mahowald et al., 1999, doi:10.1029/1999JD900084; Petit et al., 1999, doi:10.1038/20859; 1990, doi:10.1038/343056a0 Delmonte et al., 2009, doi:10.1029/2008GL033382; Lambert et al., 2008, doi:10.1038/nature06763). Isotopic tracers suggest that South America is the dominant source of the dust (Grousset et al., 1992, doi:10.1016/0012-821X(92)90177-W; Basile et al., 1997, doi:10.1016/S0012-821X(96)00255-5; Gaiero et al., 2007, doi:10.1016/j.chemgeo.2006.11.003), but it is unclear what led to the variable deposition of dust at concentrations 20-50 times higher than present in glacial-aged ice (Petit et al., 1990, doi:10.1038/343056a0; Lambert et al., 2008, doi:10.1038/nature06763). Here we characterize the age and composition of Patagonian glacial outwash sediments, to assess the relationship between the Antarctic dust record from Dome C (refs Lambert et al., 2008, doi:10.1038/nature06763; Wolff et al., 2006, doi:10.1038/nature04614) and Patagonian glacial fluctuations (Sugden et al., 2005; McCulloch et al., 2005, doi:10.1111/j.0435-3676.2005.00260.x; Kaplan et al., 2008, doi:10.1016/j.quascirev.2007.09.013) for the past 80,000 years. We show that dust peaks in Antarctica coincide with periods in Patagonia when rivers of glacial meltwater deposited sediment directly onto easily mobilized outwash plains. No dust peaks were noted when the glaciers instead terminated directly into pro-glacial lakes. We thus propose that the variable sediment supply resulting from Patagonian glacial fluctuations may have acted as an on/off switch for Antarctic dust deposition. At the last glacial termination, Patagonian glaciers quickly retreated into lakes, which may help explain why the deglacial decline in Antarctic dust concentrations preceded the main phase of warming, sea-level rise and reduction in Southern Hemisphere sea-ice extent (Wolff et al., 2006, doi:10.1038/nature04614).

Comparison of glacial/interglacial carbonate concentrations in Atlantic sediments

To assess the regional effects of glaciation on sedimentation in the Atlantic Ocean we compare sediment types, distributions, and rates between Recent (core top) and last glacial maximum (LGM: ~18,000 years B.P.) stratigraphic levels. Based upon smear slides and carbonate analyses in 178 cores we find that glacial age carbonate content is generally lower than Recent. During both the Recent and LGM, carbonate content shows an east/west asymmetry with western basins exhibiting lower carbonate values. Input of ice-rafted detritus into the North Atlantic during LGM time interrupts this topographic control on carbonate distribution considerably farther south than at present; in the South Atlantic this effect is minor. Comparison of LGM and Recent sediment distributions indicates that the LGM seafloor was dominated by biogenic oozes, calcareous clays, and clays, while the Recent is dominated by biogenic oozes and marls. Coarse-grained detritus is much more prevalent in LGM sediments, derived not only from glacial input but also from fluvial and aeolian sources. Sedimentation rates, calculated from LGM to Recent sediment thickness in cores, are <4 cm/1000 yr for most of the ocean. Higher rates are typical of the continental margin off the Amazon River, the North American Basin, and a small region off west equatorial Africa.

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.

Individual foraminifera variability from the western and eastern equatorial Pacific during the Late Holocene and Last Glacial Maximum

El Niño-Southern Oscillation (ENSO) is a major source of global interannual variability, but its response to climate change is uncertain. Paleoclimate records from the Last Glacial Maximum (LGM) provide insight into ENSO behavior when global boundary conditions (ice sheet extent, atmospheric partial pressure of CO2) were different from those today. In this work, we reconstruct LGM temperature variability at equatorial Pacific sites using measurements of individual planktonic foraminifera shells. A deep equatorial thermocline altered the dynamics in the eastern equatorial cold tongue, resulting in reduced ENSO variability during the LGM compared to the Late Holocene. These results suggest that ENSO was not tied directly to the east-west temperature gradient, as previously suggested. Rather, the thermocline of the eastern equatorial Pacific played a decisive role in the ENSO response to LGM climate.

Glacial North Atlantic: Sea surface conditions reconstructed by GLAMAP 2000

The response of the tropical ocean to global climate change and the extent of sea ice in the glacial nordic seas belong to the great controversies in paleoclimatology. Our new reconstruction of peak glacial sea surface temperatures (SSTs) in the Atlantic is based on census counts of planktic foraminifera, using the Maximum Similarity Technique Version 28 (SIMMAX-28) modern analog technique with 947 modern analog samples and 119 well-dated sediment cores. Our study compares two slightly different scenarios of the Last Glacial Maximum (LGM), the Environmental Processes of the Ice Age: Land, Oceans, Glaciers (EPILOG), and Glacial Atlantic Ocean Mapping (GLAMAP 2000) time slices. The comparison shows that the maximum LGM cooling in the Southern Hemisphere slightly preceeded that in the north. In both time slices sea ice was restricted to the north western margin of the nordic seas during glacial northern summer, while the central and eastern parts were ice-free. During northern glacial winter, sea ice advanced to the south of Iceland and Faeroe. In the central northern North Atlantic an anticyclonic gyre formed between 45° and 60°N, with a cool water mass centered west of Ireland, where glacial cooling reached a maximum of >12°C. In the subtropical ocean gyres the new reconstruction supports the glacial-to-interglacial stability of SST as shown by CLIMAP Project Members (CLIMAP) [1981]. The zonal belt of minimum SST seasonality between 2° and 6°N suggests that the LGM caloric equator occupied the same latitude as today. In contrast to the CLIMAP reconstruction, the glacial cooling of the tropical east Atlantic upwelling belt reached up to 6°-8°C during Northern Hemisphere summer. Differences between these SIMMAX-based and published U37[k]- and Mg/Ca-based equatorial SST records are ascribed to strong SST seasonalities and SST signals that were produced by different planktic species groups during different seasons.

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

Relative sea-level history of Marguerite Bay, Antarctic Peninsula derived from optically stimulated luminescence-dated beach cobbles

Calmette Bay within Marguerite Bay along the western side of the Antarctic Peninsula contains one of the most continuous flights of raised beaches described to date in Antarctica. Raised beaches extend to 40.8 m above sea level (masl) and are thought to reflect glacial isostatic adjustment due to the retreat of the Antarctic Peninsula Ice Sheet. Using optically stimulated luminescence (OSL), we dated quartz extracts from cobble surfaces buried in raised beaches at Calmette Bay. The beaches are separated into upper and lower beaches based on OSL ages, geomorphology, and sedimentary fabric. The two sets of beaches are separated by a prominent scarp. One of our OSL ages from the upper beaches dates to 9.3 thousand years ago (ka; as of 1950) consistent with previous extrapolation of sea-level data and the time of ice retreat from inner Marguerite Bay. However, four of the seven ages from the upper beaches date to the timing of glaciation. We interpret these ages to represent reworking of beaches deposited prior to the Last Glacial Maximum (LGM) by advancing and retreating LGM ice. Ages from the lower beaches record relative sea-level fall due to Holocene glacial-isostatic adjustment. We suggest a Holocene marine limit of 21.7 masl with an age of 5.5-7.3 ka based on OSL ages from Calmette Bay and other sea-level constraints in the area. A marine limit at 21.7 masl implies half as much relative sea-level change in Marguerite Bay during the Holocene as suggested by previous sea-level reconstructions. No evidence for a relative sea-level signature of neoglacial events, such as a decrease followed by an increase in RSL fall due to ice advance and retreat associated with the Little Ice Age, is found within Marguerite Bay indicating either: (1) no significant neoglacial advances occurred within Marguerite Bay; (2) rheological heterogeneity allows part of the Antarctic Peninsula (i.e. the South Shetland Islands) to respond to rapid ice mass changes while other regions are incapable of responding to short-lived ice advances; or (3) the magnitude of neoglacial events within Marguerite Bay is too small to resolve through relative sea-level reconstructions. Although the application of reconstructing sea-level histories using OSL-dated raised beach deposits provides a better understanding of the timing and nature of relative sea-level change in Marguerite Bay, we highlight possible problems associated with using raised beaches as sea-level indices due to post-depositional reworking by storm waves.