Terrestrial and submarine evidence for the extent and timing of the Last Glacial Maximum and the onset of deglaciation on the maritime-Antarctic and sub-Antarctic islands

This paper is the maritime and sub–Antarctic contribution to the Scientific Committee for Antarctic Research (SCAR) Past Antarctic Ice Sheet Dynamics (PAIS) community Antarctic Ice Sheet reconstruction. The overarching aim for all sectors of Antarctica was to reconstruct the Last Glacial Maximum (LGM) ice sheet extent and thickness, and map the subsequent deglaciation in a series of 5000 year time slices. However, our review of the literature found surprisingly few high quality chronological constraints on changing glacier extents on these timescales in the maritime and sub–Antarctic sector. Therefore, in this paper we focus on an assessment of the terrestrial and offshore evidence for the LGM ice extent, establishing minimum ages for the onset of deglaciation, and separating evidence of deglaciation from LGM limits from those associated with later Holocene glacier fluctuations. Evidence included geomorphological descriptions of glacial landscapes, radiocarbon dated basal peat and lake sediment deposits, cosmogenic isotope ages of glacial features and molecular biological data. We propose a classification of the glacial history of the maritime and sub–Antarctic islands based on this assembled evidence. These include: (Type I) islands which accumulated little or no LGM ice; (Type II) islands with a limited LGM ice extent but evidence of extensive earlier continental shelf glaciations; (Type III) seamounts and volcanoes unlikely to have accumulated significant LGM ice cover; (Type IV) islands on shallow shelves with both terrestrial and submarine evidence of LGM (and/or earlier) ice expansion; (Type V) Islands north of the Antarctic Polar Front with terrestrial evidence of LGM ice expansion; and (Type VI) islands with no data. Finally, we review the climatological and geomorphological settings that separate the glaciological history of the islands within this classification scheme.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.

A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

North Atlantic Eddy-Driven Jet in interglacial and glacial winter climates

The atmospheric westerly flow in the North Atlantic (NA) sector is dominated by atmospheric waves or eddies generating via momentum flux convergence, the so-called eddy-driven jet. The position of this jet is variable and shows for the present-day winter climate three preferred latitudinal states: a northern, central, and southernposition in the NA. Here, the authors analyze the behavior of the eddy-driven jet under different glacial and interglacial boundary conditions using atmosphere–land-only simulations with the CCSM4 climate model. As state-of-the-art climate models tend to underestimate the trimodality of the jet latitude, the authors apply a bias correction and successfully extract the trimodal behavior of the jet within CCSM4. The analysis shows that during interglacial times (i.e., the early Holocene and the Eemian) the preferred jet positions are rather stable and the observed multimodality is the typical interglacial character of the jet. During glacial times, the jet is strongly enhanced, its position is shifted southward, and the trimodal behavior vanishes. This is mainly due to the presence of the Laurentide ice sheet (LIS). The LIS enhances stationary waves downstream, thereby accelerating and displacing the NA eddy-driven jet by anomalous stationary momentum flux convergence. Additionally, changes in the transient eddy activity caused by topography changes as well as other glacial boundary conditions lead to an acceleration of the westerly winds over the southern NA at the expenseof more northernareas. Consequently, bothstationaryand transient eddiesfoster the southward shift of the NA eddy-driven jet during glacial winter times.

Glacier advances in northeastern Turkey before and during the global Last Glacial Maximum

Abstract Our study in the Başyayla Valley in northeastern Anatolia showed evidence of four glacier advances that built terminal and lateral moraines. Surface exposure dating of boulders on these moraines showed that the Maximum Ice Extent (MIE) was asynchronous with the global Last Glacial Maximum (LGM; 22.1 ± 4.3 thousand years; ka). The local {MIE} took place at least 57.0 ± 3.5 ka ago. The extent of the Başyayla Glacier during this advance is not known exactly because the boulders are only preserved on a lateral moraine. The next advance was prior to 41.5 ± 2.5 ka, and it descended down the valley to approximately 2320 m above sea level (m a.s.l.), with a glacier length of 5.3 km. During the early global LGM, the Başyayla Glacier extended for a distance of 4.9 km down to approx. 2430 m a.s.l. The last recorded advance occurred during the global LGM. This extension was 0.7 km smaller than the local {MIE} and its terminus reached 2490 m a.s.l. only. The exposure ages of boulders in a retreat position at an altitude of approx. 3045 m a.s.l. indicate that the valley has remained ice-free since the Lateglacial period. Therefore, the Lateglacial extent was limited to the cirque system in the uppermost part of the catchment. Furthermore, Holocene glacier oscillations seem to be either absent or restricted to solifluction in the whole catchment and to rock glacier movements in the southern tributary of the Başyayla Valley system.

Timing of retreat of the Reuss Glacier (Switzerland) at the end of the Last Glacial Maximum

We used cosmogenic 10Be and 36Cl to establish the timing of the onset of deglaciation after the Last Glacial Maximum of the Reuss Glacier, one of the piedmont lobes of the Alpine ice cap that reached the northern Alpine foreland in Switzerland. In this study, we sampled erratic boulders both at the frontal position in the foreland (Lenzburg and Wohlen, canton Aargau) and at the lateral Alpine border position (Seeboden moraine, Rigi, canton Schwyz). The minimum age for the beginning of retreat is 22.2 ± 1.0 ka at the frontal (terminal) position and 20.4 ± 1.0 ka at the lateral position. These ages are directly comparable with exposure ages from the other piedmont lobes in the northern Alpine foreland. Our data from the mountain called Rigi, do not support the hypothesis that boulders located external to the Seeboden moraine were deposited prior to the last glacial cycle. We present a first exposure age from an erratic boulder in a retreat position in the Alpine foreland. The Reuss Glacier was approximately 12 km behind the maximal extent no later than at 18.6 ± 0.9 ka.

Antarctic Zone nutrient conditions during the last two glacial cycles

In a sediment core from the Pacific sector of the Antarctic Zone (AZ) of the Southern Ocean, we report diatom-bound N isotope (δ15Ndb) records for total recoverable diatoms and two distinct diatom assemblages (pennate and centric rich). These data indicate tight coupling between the degree of nitrate consumption and Antarctic climate across the last two glacial cycles, with δ15Ndb (and thus the degree of nitrate consumption) increasing at each major Antarctic cooling event. Coupled with evidence from opal- and barium-based proxies for reduced export production during ice ages, the δ15Ndb increases point to ice age reductions in the supply of deep ocean-sourced nitrate to the AZ surface. The two diatom assemblages and species abundance data indicate that the δ15Ndb changes are not the result of changing species composition. The pennate and centric assemblage δ15Ndb records indicate similar changes but with a significant decline in their difference during peak ice ages. A tentative seasonality-based interpretation of the centric-to-pennate δ15Ndb difference suggests that late summer surface waters became nitrate free during the peak glacials.

Cool Glacial Temperatures and Changes in Moisture Source Recorded in Oman Groundwaters

Concentrations of atmospheric noble gases (neon, argon, krypton, and xenon) dissolved in groundwaters from northern Oman indicate that the average ground temperature during the Late Pleistocene (15,000 to 24,000 years before present) was 6.5° ± 0.6°C lower than that of today. Stable oxygen and hydrogen isotopic groundwater data show that the origin of atmospheric water vapor changed from a primarily southern, Indian Ocean source during the Late Pleistocene to a dominantly northern, Mediterranean source today. The reduced northern water vapor source is consistent with a drier Last Glacial Maximum through much of northern Africa and Arabia.

Global pulses of organic carbon burial in deep-sea sediments during glacial maxima

The burial of organic carbon in marine sediments removes carbon dioxide from the ocean–atmosphere pool, provides energy to the deep biosphere, and on geological timescales drives the oxygenation of the atmosphere. Here we quantify natural variations in the burial of organic carbon in deep-sea sediments over the last glacial cycle. Using a new data compilation of hundreds of sediment cores, we show that the accumulation rate of organic carbon in the deep sea was consistently higher (50%) during glacial maxima than during interglacials. The spatial pattern and temporal progression of the changes suggest that enhanced nutrient supply to parts of the surface ocean contributed to the glacial burial pulses, with likely additional contributions from more efficient transfer of organic matter to the deep sea and better preservation of organic matter due to reduced oxygen exposure. These results demonstrate a pronounced climate sensitivity for this global carbon cycle sink.

Deep water provenance and dynamics of the (de)glacial Atlantic meridional overturning circulation

Reconstructing past modes of ocean circulation is an essential task in paleoclimatology and paleoceanography. To this end, we combine two sedimentary proxies, Nd isotopes (εNd) and the 231Pa/230Th ratio, both of which are not directly involved in the global carbon cycle, but allow the reconstruction of water mass provenance and provide information about the past strength of overturning circulation, respectively. In this study, combined 231Pa/230Th and εNd down-core profiles from six Atlantic Ocean sediment cores are presented. The data set is complemented by the two available combined data sets from the literature. From this we derive a comprehensive picture of spatial and temporal patterns and the dynamic changes of the Atlantic Meridional Overturning Circulation over the past ∼25 ka. Our results provide evidence for a consistent pattern of glacial/stadial advances of Southern Sourced Water along with a northward circulation mode for all cores in the deeper (>3000 m) Atlantic. Results from shallower core sites support an active overturning cell of shoaled Northern Sourced Water during the LGM and the subsequent deglaciation. Furthermore, we report evidence for a short-lived period of intensified AMOC in the early Holocene.

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

Glacial-interglacial record of ODP Site 181-1119 in the Canterbury Bright, Southwest Pacific Ocean

Ocean Drilling Program (ODP) Site 1119 is located at water depth 395 m near the subtropical front (STF; here represented by the Southland Front), just downslope from the shelf edge of eastern South Island, New Zealand. The upper 86.19 metres composite depth (mcd) of Site 1119 sediment was deposited at an average sedimentation rate of 34 cm/kyr during Marine Isotope Stages (MIS) 1-8 (0-252 ka), and is underlain across a ~25 kyr intra-MIS 8 unconformity by MIS 8.5-11 (277-367 ka) and older sediment deposited at ~14 cm/kyr. A time scale is assigned to Site 1119 using radiocarbon dates for the period back to ~39 ka, and, prior to then, by matching its climatic record with that of the Vostok ice core, which it closely resembles. Four palaeoceanographic proxy measures for surface water masses vary together with the sandy-muddy, glacial-interglacial (G/I) cyclicity at the site. Interglacial intervals are characterised by heavy delta13C, high colour reflectance (a proxy for carbonate content), low Q-ray (a proxy for clay content) and light delta18O; conversely, glacial intervals exhibit light delta13C, low reflectance, high Q-ray and heavy delta18O signatures. Early interglacial intervals are represented by silty clays with 10-105-cm-thick beds of sharp-based (Chondrites-burrowed), shelly, graded, fine sand. The sands are rich in foraminifera, and were deposited distant from the shoreline under the influence of longitudinal flow in relatively deep water. Glacial intervals comprise mostly micaceous silty clay, though with some thin (2-10 cm thick) sands present also at peak cold periods, and contain the cold-water scallop Zygochlamys delicatula. Interglacial sandy intervals are characterised by relatively low sedimentation rates of 5-32 cm/kyr; cold climate intervals MIS 10, 6 and 2 have successively higher sedimentation rates of 45, 69 and 140 cm/kyr. Counter-intuitively,and forced by the bathymetric control of a laterally-moving shoreline during G/I and I/G transitions, the 1119 core records a southeasterly (seaward) movement of the STF during early glacial periods, accompanied by the incursion of subtropical water (STW) above the site, and northwesterly (landward) movement during late glacial and interglacial times, resulting in a dominant influence then of subantarctic surface water (SAW). The history of passage of these different water masses at the site is clearly delineated by their characteristic delta13C values. The intervals of thin, graded sands-muds which occur within MIS 2-3, 6, 7.4 and 10 indicate the onset at times of peak cold of intermittent bottom currents caused by strengthened and expanded frontal flows along the STF, which at such times lay near Site 1119 in close proximity to seaward-encroaching subantarctic waters within the Bounty gyre. In common with other nearby Southern Hemisphere records, the cold period which represents the last glacial maximum lasted between ~23-18 ka at Site 1119, during which time the STF and Subantarctic Front (SAF) probably merged into a single intense frontal zone around the head of the adjacent Bounty Trough.