Last Glacial Maximum Ice Research Articles

Late Holocene relative sea levels near Palmer Station, northern Antarctic Peninsula, strongly controlled by late Holocene ice-mass changes

Many studies of Holocene relative sea-level (RSL) changes across Antarctica assume that their reconstructions record uplift from glacial isostatic adjustment caused by the demise of the Last Glacial Maximum (LGM) ice sheets. However, recent analysis of GPS observations suggests that mantle viscosity beneath the Antarctic Peninsula is weaker than previously thought, which would imply that solid Earth motion is not controlled by post-LGM ice-sheet retreat but instead by late Holocene ice-mass changes. If this hypothesis is correct, one might expect to find Holocene RSL records that do not reflect a monotonic decrease in the rate of RSL fall but show variations in the rate of RSL change through the Holocene. We present a new record of late Holocene RSL change from Torgersen Island near Palmer Station in the western Antarctic Peninsula that shows an increase in the rate of relative sea-level fall from 3.0 ± 1.2 mm/yr to 5.1 ± 1.8 mm/yr during the late Holocene. Independent studies of the glacial history of the region provide evidence of ice-sheet changes over similar time scales that may be driving this change. When our RSL records are corrected for sea-surface height changes associated with glacial isostatic adjustment (GIA), the rate of post-0.79 ka land uplift at Torgersen Island, 5.3 ± 1.8 mm/yr, is much higher than the rate of uplift recorded at a nearby GPS site at Palmer Station prior to the Larsen B breakup in 2002 AD (1998-2002 AD; <0.1 mm/yr), but similar to the rates observed after 2002 AD (2002-2013 AD; 6–9 mm/yr). This substantial variation in uplift rates further supports the hypothesis that Holocene RSL rates of change are recording responses to late Holocene and recent changes in local ice loading rather than a post-LGM signal across portions of the Antarctic Peninsula. Thus middle-to-late Holocene RSL data may not be an effective tool for constraining the size of the LGM ice sheet across portions of the Antarctic Peninsula underlain by weaker mantle. In addition, current global-scale GIA models are unable to predict our observed changes in late Holocene RSL. Complexities in Earth structure and neoglacial history need to be taken into consideration in GIA models used for correcting modern satellite-based observations of ice-mass loss.

The influence of atmospheric grid resolution in a climate model-forced ice sheet simulation

Abstract. Coupled climate–ice sheet simulations have been growing in popularity in recent years. Experiments of this type are however challenging as ice sheets evolve over multi-millennial timescales, which is beyond the practical integration limit of most Earth system models. A common method to increase model throughput is to trade resolution for computational efficiency (compromise accuracy for speed). Here we analyze how the resolution of an atmospheric general circulation model (AGCM) influences the simulation quality in a stand-alone ice sheet model. Four identical AGCM simulations of the Last Glacial Maximum (LGM) were run at different horizontal resolutions: T85 (1.4∘), T42 (2.8∘), T31 (3.8∘), and T21 (5.6∘). These simulations were subsequently used as forcing of an ice sheet model. While the T85 climate forcing reproduces the LGM ice sheets to a high accuracy, the intermediate resolution cases (T42 and T31) fail to build the Eurasian ice sheet. The T21 case fails in both Eurasia and North America. Sensitivity experiments using different surface mass balance parameterizations improve the simulations of the Eurasian ice sheet in the T42 case, but the compromise is a substantial ice buildup in Siberia. The T31 and T21 cases do not improve in the same way in Eurasia, though the latter simulates the continent-wide Laurentide ice sheet in North America. The difficulty to reproduce the LGM ice sheets in the T21 case is in broad agreement with previous studies using low-resolution atmospheric models, and is caused by a substantial deterioration of the model climate between the T31 and T21 resolutions. It is speculated that this deficiency may demonstrate a fundamental problem with using low-resolution atmospheric models in these types of experiments.

Ice stream reorganization and glacial retreat on the northwest Greenland shelf

AbstractUnderstanding conditions at the grounding‐line of marine‐based ice sheets is essential for understanding ice sheet evolution. Offshore northwest Greenland, knowledge of the Last Glacial Maximum (LGM) ice sheet extent in Melville Bugt was previously based on sparse geological evidence. This study uses multibeam bathymetry, combined with 2‐D and 3‐D seismic reflection data, to present a detailed landform record from Melville Bugt. Seabed landforms include mega‐scale glacial lineations, grounding‐zone wedges, iceberg scours, and a lateral shear margin moraine, formed during the last glacial cycle. The geomorphology indicates that the LGM ice sheet reached the shelf edge before undergoing flow reorganization. After retreat of ~80 km across the outer shelf, the margin stabilized in a mid‐shelf position, possibly during the Younger Dryas (12.9–11.7 ka). The ice sheet then decoupled from the seafloor and retreated to a coast‐proximal position. This landform record provides an important constraint on deglaciation history offshore northwest Greenland.

Reconciling marine and terrestrial evidence for post LGM ice sheet retreat in southern McMurdo Sound, Antarctica

Abstract Retreat of the Antarctic ice sheets since the Last Glacial Maximum (LGM) contributed to sea-level rise, but the location, amount, and timing of ice mass loss has been controversial. This paper presents new 10Be exposure ages from glacially transported erratics which record post LGM retreat of grounded ice in the western Ross Sea. Ice elevation in southern McMurdo Sound was ≥520 m above present day sea level on the eastern side of Mount Discovery during the LGM, and the onset of major deglaciation in the region was after 14 ka. The ice surface lowered from ∼520 to 234 m above present day sea level between 14.0 ka and 10.3 ka and from 234 m to ∼30 m between 10.3 ka and 7.4 ka. This late-glacial and Holocene deglaciation chronology from southern McMurdo Sound is consistent with other records on the margins of the Ross Embayment, and implies that the western margins of the Ross Sea Ice Sheet (RSIS) experienced most mass loss during the early to middle Holocene. These 10Be exposure ages coupled with sediment provenance define a two-stage ice flow scenario for McMurdo Sound subdividing differing reconstructions into an early and late phase. Prior to Termination I, an expanded Koettlitz Glacier flowed north and northeast between Brown Peninsula and Mount Discovery and coalesced with northward flowing ice fed from the Skelton and Mulock Glaciers. Thinning and retreat of the Koettlitz Glacier and perhaps other outlet glaciers flowing through the Royal Society Range allowed ice grounded in the Ross Sea to flow westward and northward, north of Brown Peninsula. Grounding-line recession in the Ross Sea during the late-glacial and Holocene was likely driven by Southern Ocean warming and sea-level rise from the retreat of the Northern Hemisphere ice sheets and the outer margins of the Antarctic ice sheets.

Glacial isostatic uplift of the European Alps.

Following the last glacial maximum (LGM), the demise of continental ice sheets induced crustal rebound in tectonically stable regions of North America and Scandinavia that is still ongoing. Unlike the ice sheets, the Alpine ice cap developed in an orogen where the measured uplift is potentially attributed to tectonic shortening, lithospheric delamination and unloading due to deglaciation and erosion. Here we show that ∼90% of the geodetically measured rock uplift in the Alps can be explained by the Earth’s viscoelastic response to LGM deglaciation. We modelled rock uplift by reconstructing the Alpine ice cap, while accounting for postglacial erosion, sediment deposition and spatial variations in lithospheric rigidity. Clusters of excessive uplift in the Rhône Valley and in the Eastern Alps delineate regions potentially affected by mantle processes, crustal heterogeneity and active tectonics. Our study shows that even small LGM ice caps can dominate present-day rock uplift in tectonically active regions.

Biodiversity signature of the Last Glacial Maximum at South Georgia, Southern Ocean

AbstractAimHigh‐latitude biodiversity distributions can preserve signals of the timing and geography of past glaciations, and as such ground truth ice‐sheet models. Discrete polar archipelagos offer the fewest confounding factors for testing historic ice position records in extant biodiversity. At South Georgia, two competing geological hypotheses suggest that either the Last Glacial Maximum (LGM) ice was extensive, nearly covering the continental shelf (H1Big ice) or restricted to the inner fjords (H2Little ice). We examined the past configuration of the South Georgia ice cap using seabed biodiversity.LocationSouth Georgia, Southern Ocean.MethodsWe used a bespoke camera lander (SUCS) and Agassiz trawl deployments across ‘big ice’ and ‘little ice’ hypothesized positions ofLGMgrounded ice around the South Georgia continental shelf. We investigated faunal assemblage structure and richness, especially of brooders, and modelled low dispersal taxa, for example, those with limited pelagic larvae (bryozoans and sponges).ResultsWe found a striking ‘line’ of major richness discontinuity, with significantly higher richness, especially of brooders and low dispersal model taxa, mainly conforming to the ‘big ice’ hypothesized position. What few bryozoans and sponges occurred inside this line were a subset of those outside.Main conclusionsBenthic biodiversity is consistent with extensiveLGMgrounded ice advancing to near the shelf break in most, but not all locations around South Georgia's shelf, for example, the eastern shelf area. We suggest that most of the shelf is still undergoing recolonization from when grounded ice covered the shelf ~20 kyr ago. Our alternative hypothesis ofLGMice position, H3‘Limited‐Extensive ice’, best fitted our data and is easily further testable, but if verified, shows that shelf recolonization following glaciation is much slower than previously thought. This contrasts with surprisingly rapid colonization of continental shelves after ice‐shelf collapses, but these are not grounded, which may be crucial to polar recolonization rates.

Last Glacial Maximum precipitation pattern in the Alps inferred from glacier modelling

Abstract. During the Last Glacial Maximum (LGM), glaciers in the Alps reached a maximum extent, and broad sections of the foreland were covered by ice. In this study, we simulated the alpine ice cap using a glacier flow model to constrain the prevailing precipitation pattern with a geomorphological reconstruction of ice extent. For this purpose we forced the model using different temperature cooling and precipitation reduction factors. The use of the present-day precipitation pattern led to a systematic overestimation of the ice cover on the northern part of the Alps relative to the southern part. To reproduce the LGM ice cap, a more severe decrease in precipitation in the north than in the south was required. This result supports a southwesterly advection of atmospheric moisture to the Alps, sustained by a southward shift of the North Atlantic storm track during the LGM.

Cosmogenic-nuclide exposure ages from the Pensacola Mountains adjacent to the Foundation Ice Stream, Antarctica

We describe glacial-geological observations and cosmogenic-nuclide exposure ages from the Schmidt, Williams, and Thomas Hills in the Pensacola Mountains of Antarctica adjacent to the Foundation Ice Stream (FIS). Our aim is to learn about changes in the thickness and grounding line position of the Antarctic Ice Sheet in the Weddell Sea embayment between the Last Glacial Maximum (LGM) and the present. Glacial-geological observations from all three regions indicate that currently-ice-free areas were covered by ice during one or more past ice sheet expansions, and that this ice was typically frozen to its bed and thus non-erosive, permitting the accumulation of multiple generations of glacial drift. Cosmogenic-nuclide exposure-age data from glacially transported erratics are consistent with this interpretation in that we observe both (i) samples with Holocene exposure ages that display a systematic age-elevation relationship recording LGM-to-present deglaciation, and (ii) samples with older and highly scattered apparent exposure ages that were deposited in previous glacial-interglacial cycles and have experienced multiple periods of surface exposure and ice cover. Holocene exposure ages at the Thomas and Williams Hills, upstream of the present grounding line of the FIS, show that the FIS was at least 500 m thicker prior to 11 ka, and that 500 m of thinning took place between 11 and 4 ka. However, exposure-age data from the Schmidt Hills, downstream of the present grounding line of the FIS, show no evidence for LGM thickening of the FIS and, in fact, provide some evidence that the FIS could have been no more than 200 m thicker than present at the LGM. If all these observations are correct, they imply that the LGM and early Holocene ice surface slope in the vicinity of the present grounding line was steeper than present, which is inconsistent with glaciological model predictions of possible LGM ice sheet configurations. Specifically, scenarios in which the LGM grounding line of the FIS advanced to the outer continental shelf appear inconsistent with exposure-age data from the Schmidt Hills, whereas scenarios in which the FIS grounding line did not advance at the LGM appear inconsistent with exposure-age data from the Williams and Thomas Hills.

Dating the onset of LGM ice surface lowering in the High Alps

The general chronological outline of the Last Glacial Maximum (LGM) in the Alps is well known after more than two centuries of ice age research in Europe. Current studies focus on resolving the details of the deglaciation process after the LGM. Particularly few data in this context are available from the High Alps. Here we report chronological constraints on the onset of deglaciation in two different study areas from the Western and Eastern Alps: the Mont Blanc (Italy) and Zillertal Alps (Austria). We sampled 32 bedrock and boulder surfaces from high elevation ridges beneath the glacial trimline. The oldest 10Be exposure ages of ∼18.5 ka correspond to the initiation of lowering of the LGM ice surface. Identical ages from both study areas suggest synchronous decay of the LGM glaciers in the accumulation zones of the Western and Eastern Alps. Thus, the High Alpine ice surface lowered roughly synchronously to the downwasting of the glacier tongues in the forelands that was completed 19–18 ka. In both study areas, ages of numerous samples that are consistently 2–3 ka younger than 18.5 ka suggest the persistance of remnant ice patches at high elevations in the Lateglacial accumulation zones approximately until the Bølling-Allerød interstadial.

Rapid Last Glacial Maximum deglaciation in the Indian Himalaya coeval with midlatitude glaciers: New insights from 10Be‐dating of ice‐polished bedrock surfaces in the Chandra Valley, NW Himalaya

AbstractDespite a large number of dated glacial landforms in the Himalaya, the ice extent during the global Last Glacial Maximum (LGM) from 19 to 23 ka is only known to first order. New cosmogenic 10Be exposure ages from well‐preserved glacially polished surfaces, combined with published data, and an improved production rate scaling model allow reconstruction of the LGM ice extent and subsequent deglaciation in the Chandra Valley of NW India. We show that a &gt;1000 m thick valley glacier retreated &gt;150 km within a few thousand years after the onset of LGM deglaciation. By comparing the recession of the Chandra Valley Glacier and other Himalayan glaciers with those of Northern and Southern Hemisphere glaciers, we demonstrate that post‐LGM deglaciation was similar and nearly finished prior to the Bølling/Allerød interstadial. Our study supports the view that many Himalayan glaciers advanced during the LGM, likely in response to global variations in temperature.

Extensive glaciation in Transbaikalia, Siberia, at the Last Glacial Maximum

Successively smaller glacial extents have been proposed for continental Eurasia during the stadials of the last glacial period leading up to the Last Glacial Maximum (LGM). At the same time the large mountainous region east of Lake Baikal, Transbaikalia, has remained unexplored in terms of glacial chronology despite clear geomorphological evidence of substantial past glaciations. We have applied cosmogenic 10Be exposure dating and optically stimulated luminescence to establish the first quantitative glacial chronology for this region. Based on eighteen exposure ages from five moraine complexes, we propose that large mountain ice fields existed in the Kodar and Udokan mountains during Oxygen Isotope Stage 2, commensurate with the global LGM. These ice fields fed valley glaciers (>100 km in length) reaching down to the Chara Depression between the Kodar and Udokan mountains and to the valley of the Vitim River northwest of the Kodar Mountains. Two of the investigated moraines date to the Late Glacial, but indications of incomplete exposure among some of the sampled boulders obscure the specific details of the post-LGM glacial history. In addition to the LGM ice fields in the highest mountains of Transbaikalia, we report geomorphological evidence of a much more extensive, ice-cap type glaciation at a time that is yet to be firmly resolved.

Antarctic marine ice-sheet retreat in the Ross Sea during the early Holocene

Geological constraints on the timing of retreat of the Last Glacial Maximum (LGM) Antarctic Ice Sheets provide critical insights into the processes controlling marine-based ice-sheet retreat. The over-deepened, landward-sloping bathymetry of Antarctica's continental shelves is an ideal configuration for marine ice-sheet instability, with the potential for past and future ice-sheet collapse and accelerated sea-level rise. However, the chronology of retreat of the LGM ice sheet in the Ross Sea is largely constrained by imprecise radiocarbon chronology of bulk marine sediments or by coastal records that offer more reliable dating techniques but which may be influenced by local piedmont glaciers derived from East Antarctic outlet glaciers. Consequently, these coastal records may be ambiguous in the broader context of retreat in the central regions of the Ross Sea. Here, we present a sedimentary facies succession and foraminifera-based radiocarbon chronology from within the Ross Sea embayment that indicates glacial retreat and open-marine conditions to the east of Ross Island before 8.6 cal. (calibrated) kyr B.P., at least 1 k.y. earlier than indicated by terrestrial records in McMurdo Sound. Comparing these data to new modeling experiments, we hypothesize that marine-based ice-sheet retreat was triggered by oceanic forcings along most of the Pacific Ocean coastline of Antarctica, but continued Holocene retreat into the inner shelf region of the Ross Sea occurred primarily as a consequence of bathymetric controls on marine ice-sheet instability.

How might the North American ice sheet influence the northwestern Eurasian climate?

Abstract. It is now widely acknowledged that past Northern Hemisphere ice sheets covering Canada and northern Europe at the Last Glacial Maximum (LGM) exerted a strong influence on climate by causing changes in atmospheric and oceanic circulations. In turn, these changes may have impacted the development of the ice sheets themselves through a combination of different feedback mechanisms. The present study is designed to investigate the potential impact of the North American ice sheet on the surface mass balance (SMB) of the Eurasian ice sheet driven by simulated changes in the past glacial atmospheric circulation. Using the LMDZ5 atmospheric circulation model, we carried out 12 experiments under constant LGM conditions for insolation, greenhouse gases and ocean. In these experiments, the Eurasian ice sheet is removed. The 12 experiments differ in the North American ice-sheet topography, ranging from a white and flat (present-day topography) ice sheet to a full-size LGM ice sheet. This experimental design allows the albedo and the topographic impacts of the North American ice sheet onto the climate to be disentangled. The results are compared to our baseline experiment where both the North American and the Eurasian ice sheets have been removed. In summer, the sole albedo effect of the American ice sheet modifies the pattern of planetary waves with respect to the no-ice-sheet case, resulting in a cooling of the northwestern Eurasian region. By contrast, the atmospheric circulation changes induced by the topography of the North American ice sheet lead to a strong decrease of this cooling. In winter, the Scandinavian and the Barents–Kara regions respond differently to the American ice-sheet albedo effect: in response to atmospheric circulation changes, Scandinavia becomes warmer and total precipitation is more abundant, whereas the Barents–Kara area becomes cooler with a decrease of convective processes, causing a decrease of total precipitation. The gradual increase of the altitude of the American ice sheet leads to less total precipitation and snowfall and to colder temperatures over both the Scandinavian and the Barents and Kara sea sectors. We then compute the resulting annual surface mass balance over the Fennoscandian region from the simulated temperature and precipitation fields used to force an ice-sheet model. It clearly appears that the SMB is dominated by the ablation signal. In response to the summer cooling induced by the American ice-sheet albedo, high positive SMB values are obtained over the Eurasian region, leading thus to the growth of an ice sheet. On the contrary, the gradual increase of the American ice-sheet altitude induces more ablation over the Eurasian sector, hence limiting the growth of Fennoscandia. To test the robustness of our results with respect to the Eurasian ice sheet state, we carried out two additional LMDZ experiments with new boundary conditions involving both the American (flat or full LGM) and high Eurasian ice sheets. The most striking result is that the Eurasian ice sheet is maintained under full-LGM North American ice-sheet conditions, but loses ~ 10 % of its mass compared to the case in which the North American ice sheet is flat. These new findings qualitatively confirm the conclusions from our first series of experiments and suggest that the development of the Eurasian ice sheet may have been slowed down by the growth of the American ice sheet, offering thereby a new understanding of the evolution of Northern Hemisphere ice sheets throughout glacial–interglacial cycles.

Early break-up of the Norwegian Channel Ice Stream during the Last Glacial Maximum

We present 18 new cosmogenic 10Be exposure ages that constrain the breakup time of the Norwegian Channel Ice Stream (NCIS) and the initial retreat of the Scandinavian Ice Sheet from the Southwest coast of Norway following the Last Glacial Maximum (LGM). Seven samples from glacially transported erratics on the island Utsira, located in the path of the NCIS about 400 km up-flow from the LGM ice front position, yielded an average 10Be age of 22.0 ± 2.0 ka. The distribution of the ages is skewed with the 4 youngest all within the range 20.2–20.8 ka. We place most confidence on this cluster of ages to constrain the timing of ice sheet retreat as we suspect the 3 oldest ages have some inheritance from a previous ice free period. Three additional ages from the adjacent island Karmøy provided an average age of 20.9 ± 0.7 ka, further supporting the new timing of retreat for the NCIS. The 10Be ages from Utsira and Karmøy suggest that the ice stream broke up about 2000 years earlier than the age assignment based on 14C ages on foraminifera and molluscs from marine sediment cores. We postulate that the Scandinavian Ice Sheet flowed across the Norwegian Channel to Denmark and onto the North Sea plateau during early phases of the LGM. When the NCIS started to operate this ice supply to the North Sea was cut off and the fast flow of the NCIS also led to a lowering of the ice surface along the Norwegian Channel and thereby drawdown of the entire ice sheet. This facilitated rapid calving of the ice front in the North Sea and we reconstruct a large open bay across the entire northern North Sea by ∼20 ka based on our 10Be ages in the east and radiocarbon ages from marine cores in the west. Additional 10Be ages show that the mainland slightly east of the islands Utsira and Karmøy remained ice covered until about 16 ka, indicating almost no net ice-margin retreat for the 4000 years between 20 and 16 ka. After 16 ka the ice margin retreated quickly up-fjord.

A data-constrained large ensemble analysis of Antarctic evolution since the Eemian

Reconstructing the historical behaviour of the Antarctic Ice Sheet is essential for understanding and predicting climatic interactions, ice sheet sensitivities, and sea level contributions. Furthermore, to interpret such reconstructions with any degree of confidence, meaningful uncertainty estimates are required. Working toward this goal, this article presents results from a large-ensemble data-constrained study of Antarctic evolution over the Last Glacial cycle.Ice sheet chronologies have been generated using a 3-D glacial systems model that includes sheet/stream/shelf flow; a parametrized basal drag coefficient (accounting for sub-grid topographic roughness, sediment likelihood, and systematic model-to-observation thickness misfits); a sub-grid grounding-line flux parametrization; a visco-elastic bedrock response component; parametrized climate forcing; separate shelf and tidewater calving treatments; and a physically-motivated, empirically derived, sub-shelf melt component. There are 31 ensemble parameters to capture uncertainties in the glacial cycle climate, mass-balance processes, and ice dynamics. Once generated, the ensemble of model runs is constrained using a suite of observational data (including constraints on relative sea level, past ice thickness and grounding line retreat data-points, as well as present-day ice volumes and configuration) and an evaluation methodology that produces a misfit-to-observations score for each run. Assuming variants of a Gaussian error model, the scores are used to generate probability distributions for the past evolution of the ice sheet.In our ensemble-based analysis the Last Glacial Maximum (LGM) of Antarctic ice volume occurs at ∼25–24 ka. The LGM ice volume excess relative to present-day is likely between 5.6 and 14.3 m equivalent sea level (mESL), and with less confidence >10 mESL. There is little change in the grounded ice volume from 24 ka until 17–16 ka at which time widespread deglaciation commences. The Ross and Weddell Sea sectors are the largest sources of ice mass loss. The major period of grounding line retreat in the Weddell Sea sector, resulting in the present-day Ronne–Filchner shelf system, occurs after 10 ka; the grounding line in the Ross Sea sector start its major retreat phase after 12 ka. There is no significant contribution (<1.8 mESL) to meltwater pulse 1a.During the Eemian, the modelled Antarctic Ice Sheet has a minimal ice volume at 114 ka contributing ∼2 mESL relative to present-day. The mass loss is principally due to a major retreat of the grounding line in the Siple Coast sector of the Ross Sea.We present standard deviation plots of the ensemble results for present day and at the LGM. These plots highlight geographic areas that would most benefit from additional constraint and therefore offer priorities and guidance for future field campaigns.

The sensitivity of GNSS measurements in Fennoscandia to distinct three-dimensional upper-mantle structures

Abstract. The sensitivity of global navigation satellite system (GNSS) measurements in Fennoscandia to nearby viscosity variations in the upper mantle is investigated using a 3-D finite element model of glacial isostatic adjustment (GIA). Based on the lateral viscosity structure inferred from seismic tomography and the location of the ice margin at the last glacial maximum (LGM), the GIA earth model is subdivided into four layers, where each of them contains an amalgamation of about 20 blocks of different shapes in the central area. The sensitivity kernels of the three velocity components at 10 selected GNSS stations are then computed for all the blocks. We find that GNSS stations within the formerly glaciated area are most sensitive to mantle viscosities below and in its near proximity, i.e., within about 250 km in general. However, this can be as large as 1000 km if the stations lie near the center of uplift. The sensitivity of all stations to regions outside the ice margin during the LGM is generally negligible. In addition, it is shown that prominent structures in the second (250–450 km depth) and third layers (450–550 km depth) of the upper mantle may be readily detected by GNSS measurements, while the viscosity in the first mantle layer below the lithosphere (70–250 km depth) along the Norwegian coast, which is related to lateral lithospheric thickness variation there, can also be detected but with limited sensitivity. For future investigations on the lateral viscosity structure, preference should be on GNSS stations within the LGM ice margin. But these stations can be grouped into clusters to improve the inference of viscosity in a specific area. However, the GNSS measurements used in such inversion should be weighted according to their sensitivity. Such weighting should also be applied when they are used in combination with other GIA data (e.g., relative sea-level and gravity data) for the inference of mantle viscosity.

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.

Retreat history of the East Antarctic Ice Sheet since the Last Glacial Maximum

The East Antarctic Ice Sheet (EAIS) is the largest continental ice mass on Earth, and documenting its evolution since the Last Glacial Maximum (LGM) is important for understanding its present-day and future behaviour. As part of a community effort, we review geological evidence from East Antarctica that constrains the ice sheet history throughout this period (∼30,000 years ago to present). This includes terrestrial cosmogenic nuclide dates from previously glaciated regions, 14C chronologies from glacial and post-glacial deposits onshore and on the continental shelf, and ice sheet thickness changes inferred from ice cores and continental-scale ice sheet models. We also include new 14C dates from the George V Land – Terre Adélie Coast shelf. We show that the EAIS advanced to the continental shelf margin in some parts of East Antarctica, and that the ice sheet characteristically thickened by 300–400 m near the present-day coastline at these sites. This advance was associated with the formation of low-gradient ice streams that grounded at depths of >1 km below sea level on the inner continental shelf. The Lambert/Amery system thickened by a greater amount (800 m) near its present-day grounding zone, but did not advance beyond the inner continental shelf. At other sites in coastal East Antarctica (e.g. Bunger Hills, Larsemann Hills), very little change in the ice sheet margin occurred at the LGM, perhaps because ice streams accommodated any excess ice build up, leaving adjacent, ice-free areas relatively unaffected. Evidence from nunataks indicates that the amount of ice sheet thickening diminished inland at the LGM, an observation supported by ice cores, which suggest that interior ice sheet domes were ∼100 m lower than present at this time. Ice sheet recession may have started ∼18,000 years ago in the Lambert/Amery glacial system, and by ∼14,000 years ago in Mac.Robertson Land. These early pulses of deglaciation may have been responses to abrupt sea-level rise events such as Meltwater Pulse 1a, destabilising the margins of the ice sheet. It is unlikely, however, that East Antarctica contributed more than ∼1 m of eustatic sea-level equivalent to post-glacial meltwater pulses. The majority of ice sheet recession occurred after Meltwater Pulse 1a, between ∼12,000 and ∼6000 years ago, during a period when the adjacent ocean warmed significantly. Large tracts of East Antarctica remain poorly studied, and further work is required to develop a robust understanding of the LGM ice sheet expansion, and its subsequent contraction. Further work will also allow the contribution of the EAIS to post-glacial sea-level rise, and present-day estimates of glacio-isostatic adjustment to be refined.

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.

Influence of Last Glacial Maximum boundary conditions on the global water isotope distribution in an atmospheric general circulation model

Abstract. To understand the validity of δ18O 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 (δ18Oprecip) 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 δ18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and δ18Oprecip 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 δ18Oprecip 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 δ18Oprecip distribution among the simulations.