Ice Sheet Changes Research Articles

Direct observations of submarine melt and subsurface geometry at a tidewater glacier.

Ice loss from the world's glaciers and ice sheets contributes to sea level rise, influences ocean circulation, and affects ecosystem productivity. Ongoing changes in glaciers and ice sheets are driven by submarine melting and iceberg calving from tidewater glacier margins. However, predictions of glacier change largely rest on unconstrained theory for submarine melting. Here, we use repeat multibeam sonar surveys to image a subsurface tidewater glacier face and document a time-variable, three-dimensional geometry linked to melting and calving patterns. Submarine melt rates are high across the entire ice face over both seasons surveyed and increase from spring to summer. The observed melt rates are up to two orders of magnitude greater than predicted by theory, challenging current simulations of ice loss from tidewater glaciers.

Modelling the Antarctic Ice Sheet across the mid-Pleistocene transition – implications for Oldest Ice

Abstract. The international endeavour to retrieve a continuous ice core, which spans the middle Pleistocene climate transition ca. 1.2–0.9 Myr ago, encompasses a multitude of field and model-based pre-site surveys. We expand on the current efforts to locate a suitable drilling site for the oldest Antarctic ice core by means of 3-D continental ice-sheet modelling. To this end, we present an ensemble of ice-sheet simulations spanning the last 2 Myr, employing transient boundary conditions derived from climate modelling and climate proxy records. We discuss the imprint of changing climate conditions, sea level and geothermal heat flux on the ice thickness, and basal conditions around previously identified sites with continuous records of old ice. Our modelling results show a range of configurational ice-sheet changes across the middle Pleistocene transition, suggesting a potential shift of the West Antarctic Ice Sheet to a marine-based configuration. Despite the middle Pleistocene climate reorganisation and associated ice-dynamic changes, we identify several regions conducive to conditions maintaining 1.5 Myr (million years) old ice, particularly around Dome Fuji, Dome C and Ridge B, which is in agreement with previous studies. This finding strengthens the notion that continuous records with such old ice do exist in previously identified regions, while we are also providing a dynamic continental ice-sheet context.

ICE SHEET ELEVATION MAPPING AND CHANGE DETECTION WITH THE ICE, CLOUD AND LAND ELEVATION SATELLITE-2

Abstract. On September 15, 2018, ICESat-2 (Ice, Cloud, and land Elevation satellite) was successfully launched to measure ice sheet and glacier elevation change, sea ice freeboard, and vegetation. This paper describes the computation of surface elevation change rates obtained with SERAC (Surface Elevation Reconstruction And Change detection) from ICESat-2 observations. After summarizing some relevant aspects of ICESat-2 and its sole instrument ATLAS (Advanced Topographic Laser Altimetry System) the paper focuses on how we calculate time series of elevation change rates from ICESat-2’s data product ATL03. Since real ICESat-2 data suitable for generating time series of several time epochs are not yet available, we used simulated data for this study. We will start generating time series from real ICESat-2 data after the conclusion of the ongoing calibration and validation phase and we expect to present real-world examples at the WG III/9 meeting in June, 2019 in Enschede, The Netherlands.

Deglacial Thinning of the Laurentide Ice Sheet in the Adirondack Mountains, New York, USA, Revealed by 36Cl Exposure Dating

AbstractFuture changes in sea level will largely be dictated by changes in the world's ice sheets. Yet the magnitude and rate at which these ice sheets will respond to climate change remain uncertain, necessitating a deeper investigation into past ice sheet‐climate interactions. Numerous studies have documented the timing and pattern of Laurentide Ice Sheet margin retreat since the Last Glacial Maximum, but few studies have provided vertical constraints necessary for accurate sea level contribution estimates. Here we present 21 36Cl ages from boulder and bedrock samples along vertical transects spanning ~1,000 m of relief from multiple peaks in the Adirondack Mountains of northeastern New York, USA. Our exposure ages span the Last Glacial Maximum through the last deglaciation, with the highest‐elevation sites (~1,500 m) ranging between 25 and 19 ka, and the lower elevation sites (≤1,300 m) between 16 and 13 ka. These data suggest gradual ice sheet thinning of 200 m initiated at ~20 ka followed by more rapid surface lowering of 1,000 m approximately coincident with Bølling‐Allerød warming.

Downscaling GRACE Predictions of the Crustal Response to the Present‐Day Mass Changes in Greenland

AbstractThe GRACE mission has had a revolutionary impact on the study of Earth system processes, but it provides a band‐limited representation of mass changes. This is particularly problematic when studying mass changes that tend to be concentrated in fairly narrow zones near the edges of the Greenland Ice Sheet. In this study, coarse‐resolution estimates of the mass change derived from GRACE have been enhanced by the introduction of heuristic scaling factors applied to model surface mass balance and observed surface elevation change. Corresponding results indicate large spatial heterogeneity in the gridded scaling factors at the 0.5° × 0.5° scale, reflecting significant mass losses concentrated along the ice sheet margin and relatively small internal ice sheet changes at higher elevations. The scaled GRACE‐derived vertical displacements are in the range from −2 to 14 mm/year from 2003 to 2015. The Greenland GPS network was used to examine the downscaling GRACE predictions of the crustal displacements. The results show consistency of scaled GRACE‐predicted and GPS‐observed seasonal and long‐term uplift in major drainage basins of Greenland. Our results indicate that GRACE predictions underestimate vertical displacements at sites located in regions characterized by concentrated loads, but perform well in other regions. Differences between predicted and observed uplift rates are mainly caused by the sensitivity kernels, because GPS and GRACE estimates are based on weighted averages of mass loss in different sensitivity ranges. Moreover, a large uncertainty in the glacial isostatic adjustment correction may also cause errors in the GPS‐to‐GRACE ratio.

On recovering Last Interglacial changes in the Antarctic ice sheet

On recovering Last Interglacial changes in the Antarctic ice sheet

Holocene Indian Ocean sea level, Antarctic melting history and past Tsunami deposits inferred using sea level reconstructions from the Sri Lankan, Southeastern Indian and Maldivian coasts

Holocene sea level change in the northern Indian Ocean was studied using geochemical and geophysical approaches. Molluscs were sub-sampled for radiocarbon dating from sediment cores retrieved from a south Sri Lankan coastal lagoon. They were then combined with previously published sediment core radiocarbon ages from the same lagoon. We observe a ∼20-fold reduction in sedimentation rates at around 4 ka, attributed to a decrease in the rate of sea level rise. Previously reported radiocarbon ages based on total organic matter in the lagoon show remarkable agreement with our new dates suggesting stable and homogeneous sedimentation during the last 8 ka. Comparison of down core age trends of both data sets reveal that sediments are older than coeval molluscs, except for the horizons adjacent to the paleo Tsunami layers, indicating fresh and younger sediments are provided during Tsunami events, which could become a new proxy for identifying paleo Tsunamis. Compilation of previously reported sea level indicators from Sri Lanka, Southeastern India and the Maldives, together with predicted sea level obtained from a glacio-hydro-isostatic adjustment model (GIA), suggest that 3–4 m of global sea level equivalent ice sheet melting occurred during the Mid Holocene due to the retreat of the Antarctic and/or Greenland ice sheets. Previous works suggests late Holocene (ca. 4 ka) climate anomalies in both the low and high latitudes. We suggest the low latitude climate anomaly, transmitted via atmosphere to the high latitude during the late Holocene, seems to have induced changes in polar ice sheets.

The GRISLI ice sheet model (version 2.0): calibration and validation for multi-millennial changes of the Antarctic ice sheet

Abstract. In this paper, we present the GRISLI (Grenoble ice sheet and land ice) model in its newest revision (version 2.0). Whilst GRISLI is applicable to any given ice sheet, we focus here on the Antarctic ice sheet because it highlights the importance of grounding line dynamics. Important improvements have been implemented in the model since its original version (Ritz et al., 2001). Notably, GRISLI now includes a basal hydrology model and an explicit flux computation at the grounding line based on the analytical formulations of Schoof (2007) or Tsai et al. (2015). We perform a full calibration of the model based on an ensemble of 300 simulations sampling mechanical parameter space using a Latin hypercube method. Performance of individual members is assessed relative to the deviation from present-day observed Antarctic ice thickness. To assess the ability of the model to simulate grounding line migration, we also present glacial–interglacial ice sheet changes throughout the last 400 kyr using the best ensemble members taking advantage of the capacity of the model to perform multi-millennial long-term integrations. To achieve this goal, we construct a simple climatic perturbation of present-day climate forcing fields based on two climate proxies: atmospheric and oceanic. The model is able to reproduce expected grounding line advances during glacial periods and subsequent retreats during terminations with reasonable glacial–interglacial ice volume changes.

Comparison of four calving laws to model Greenland outlet glaciers

Abstract. Calving is an important mechanism that controls the dynamics of marine terminating glaciers of Greenland. Iceberg calving at the terminus affects the entire stress regime of outlet glaciers, which may lead to further retreat and ice flow acceleration. It is therefore critical to accurately parameterize calving in ice sheet models in order to improve the projections of ice sheet change over the coming decades and reduce the uncertainty in their contribution to sea-level rise. Several calving laws have been proposed, but most of them have been applied only to a specific region and have not been tested on other glaciers, while some others have only been implemented in 1-D flowline or vertical flowband models. Here, we test and compare several calving laws recently proposed in the literature using the Ice Sheet System Model (ISSM). We test these calving laws on nine tidewater glaciers of Greenland. We compare the modeled ice front evolution to the observed retreat from Landsat data collected over the past 10 years, and assess which calving law has better predictive abilities for each glacier. Overall, the von Mises tensile stress calving law is more satisfactory than other laws for simulating observed ice front retreat, but new parameterizations that better capture the different modes of calving should be developed. Although the final positions of ice fronts are different for forecast simulations with different calving laws, our results confirm that ice front retreat highly depends on bed topography, irrespective of the calving law employed. This study also confirms that calving dynamics needs to be 3-D or in plan view in ice sheet models to account for complex bed topography and narrow fjords along the coast of Greenland.

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.

Simulating the Last Interglacial Greenland stable water isotope peak: The role of Arctic sea ice changes

Last Interglacial (LIG), stable water isotope values (δ18O) measured in Greenland deep ice cores are at least 2.5‰ higher compared to the present day. Previous isotopic climate simulations of the LIG do not capture the observed Greenland δ18O increases. Here, we use the isotope-enabled HadCM3 (UK Met Office coupled atmosphere-ocean general circulation model) to investigate whether a retreat of Northern Hemisphere sea ice was responsible for this model-data disagreement. Our results highlight the potential significance of sea ice changes on the LIG Greenland isotopic maximum. Sea ice loss in combination with increased sea surface temperatures, over the Arctic, affect δ18O: water vapour enriched in heavy isotopes and a shorter distillation path may both increase δ18O values over Greenland. We show, for the first time, that simulations of the response to Arctic sea ice reduction are capable of producing the likely magnitude of LIG δ18O increases at NEEM, NGRIP, GIPS2 and Camp Century ice core sites. However, we may underestimate δ18O changes at the Renland, DYE3 and GRIP ice core locations. Accounting for possible ice sheet changes is likely to be required to produce a better fit to the LIG ice core δ18O values.

Is the Last Glacial Maximum a reverse analog for future hydroclimate changes in the Americas?

Future hydroclimate change is expected to generally follow a wet-get-wetter, dry-get-drier (WWDD) pattern, yet key uncertainties remain regionally and over land. It has been previously hypothesized that lake levels of the Last Glacial Maximum (LGM) could map a reverse analog to future hydroclimate changes due to reduction of CO2 levels at this time. Potential complications to this approach include, however, the confounding effects of factors such as the Laurentide Ice Sheet and lake evaporation changes. Using the ensemble output of six coupled climate models, lake energy and water balance models, an atmospheric moisture budget analysis, and additional CO2 sensitivity experiments, we assess the effectiveness of the LGM as a reverse analog for future hydroclimate changes for a transect from the drylands of North America to southern South America. The model ensemble successfully simulates the general pattern of lower tropical lake levels and higher extratropical lake levels at LGM, matching 82% of the lake proxy records. The greatest model-data mismatch occurs in tropical and extratropical South America, potentially as a result of underestimated changes in temperature and surface evaporation. Thermodynamic processes of the mean circulation best explain the direction of lake changes observed in the proxy record, particularly in the tropics and Pacific coasts of the extratropics, and produce a WWDD pattern. CO2 forcing alone cannot account for LGM lake level changes, however, as the enhanced cooling from the Laurentide ice sheet appears necessary to generate LGM dry anomalies in the tropics and to deepen anomalies in the extratropics. LGM performance as a reverse analog is regionally dependent as anti-correlation between LGM and future P − E is not uniformly observed across the study domain.

A new global GPS data set for testing and improving modelled GIA uplift rates

Glacial isostatic adjustment (GIA) is the response of the solid Earth to past ice loading, primarily, since the Last Glacial Maximum, about 20 K yr BP. Modelling GIA is challenging because of large uncertainties in ice loading history and also the viscosity of the upper and lower mantle. GPS data contain the signature of GIA in their uplift rates but these also contain other sources of vertical land motion (VLM) such as tectonics, human and natural influences on water storage that can mask the underlying GIA signal. In this study, we use about 4000 GPS vertical velocities as observational estimates of global GIA uplift rates, after correcting for major elastic deformation effects. A novel fully automatic strategy is developed to postprocess the GPS time-series and to correct for non-GIA artefacts. Before estimating vertical velocities and uncertainties, we detect outliers and jumps and correct for atmospheric mass loading displacements. We correct the resulting velocities for the elastic response of the solid Earth to global changes in ice sheets, glaciers and ocean loading, as well as for changes in the Earth’s rotational pole relative to the 20th century average. We then apply a spatial median filter to remove sites where local effects are dominant to leave approximately 4000 GPS sites. The resulting novel global GPS data set shows a clean GIA signal at all post-processed stations and is therefore suitable to investigate the behavior of global GIA forward models. The results are transformed from a frame with its origin in the centre of mass of the total Earth’s system (CM) into a frame with its origin in the centre of mass of the solid Earth (CE) before comparison with 13 global GIA forward model solutions, with best fits with Pur-6- VM5 and ICE-6G predictions. The largest discrepancies for all models were identified for Antarctica and Greenland, which may be due to either uncertain mantle rheology, ice loading history/magnitude and/or GPS errors.

The devil's in the disequilibrium: multi-component analysis of dissolved carbon and oxygen changes under a broad range of forcings in a general circulation model

Abstract. The complexity of dissolved gas cycling in the ocean presents a challenge for mechanistic understanding and can hinder model intercomparison. One helpful approach is the conceptualization of dissolved gases as the sum of multiple, strictly defined components. Here we decompose dissolved inorganic carbon (DIC) into four components: saturation (DICsat), disequilibrium (DICdis), carbonate (DICcarb), and soft tissue (DICsoft). The cycling of dissolved oxygen is simpler, but can still be aided by considering O2, O2sat, and O2dis. We explore changes in these components within a large suite of simulations with a complex coupled climate–biogeochemical model, driven by changes in astronomical parameters, ice sheets, and radiative forcing, in order to explore the potential importance of the different components to ocean carbon storage on long timescales. We find that both DICsoft and DICdis vary over a range of 40 µmol kg−1 in response to the climate forcing, equivalent to changes in atmospheric pCO2 on the order of 50 ppm for each. The most extreme values occur at the coldest and intermediate climate states. We also find significant changes in O2 disequilibrium, with large increases under cold climate states. We find that, despite the broad range of climate states represented, changes in global DICsoft can be quantitatively approximated by the product of deep ocean ideal age and the global export production flux. In contrast, global DICdis is dominantly controlled by the fraction of the ocean filled by Antarctic Bottom Water (AABW). Because the AABW fraction and ideal age are inversely correlated among the simulations, DICdis and DICsoft are also inversely correlated, dampening the overall changes in DIC. This inverse correlation could be decoupled if changes in deep ocean mixing were to alter ideal age independently of AABW fraction, or if independent ecosystem changes were to alter export and remineralization, thereby modifying DICsoft. As an example of the latter, we show that iron fertilization causes both DICsoft and DICdis to increase and that the relationship between these two components depends on the climate state. We propose a simple framework to consider the global contribution of DICsoft+DICdis to ocean carbon storage as a function of the surface preformed nitrate and DICdis of dense water formation regions, the global volume fractions ventilated by these regions, and the global nitrate inventory.

The glacial isostatic adjustment signal at present day in northern Europe and the British Isles estimated from geodetic observations and geophysical models

Abstract. The glacial isostatic adjustment (GIA) signal at present day is constrained via the joint inversion of geodetic observations and GIA models for a region encompassing northern Europe, the British Isles, and the Barents Sea. The constraining data are Global Positioning System (GPS) vertical crustal velocities and GRACE (Gravity Recovery and Climate Experiment) gravity data. When the data are inverted with a set of GIA models, the best-fit model for the vertical motion signal has a χ2 value of approximately 1 and a maximum a posteriori uncertainty of 0.3–0.4 mm yr−1. An elastic correction is applied to the vertical land motion rates that accounts for present-day changes to terrestrial hydrology as well as recent mass changes of ice sheets and glaciered regions. Throughout the study area, mass losses from Greenland dominate the elastic vertical signal and combine to give an elastic correction of up to +0.5 mm yr−1 in central Scandinavia. Neglecting to use an elastic correction may thus introduce a small but persistent bias in model predictions of GIA vertical motion even in central Scandinavia where vertical motion is dominated by GIA due to past glaciations. The predicted gravity signal is generally less well-constrained than the vertical signal, in part due to uncertainties associated with the correction for contemporary ice mass loss in Svalbard and the Russian Arctic. The GRACE-derived gravity trend is corrected for present-day ice mass loss using estimates derived from the ICESat and CryoSat missions, although a difference in magnitude between GRACE-inferred and altimetry-inferred regional mass loss rates suggests the possibility of a non-negligible GIA response here either from millennial-scale or Little Ice Age GIA.

An Overview of Interactions and Feedbacks Between Ice Sheets and the Earth System

AbstractIce sheet response to forced changes—such as that from anthropogenic climate forcing—is closely regulated by two‐way interactions with other components of the Earth system. These interactions encompass the ice sheet response to Earth system forcing, the Earth system response to ice sheet change, and feedbacks resulting from coupled ice sheet/Earth system evolution. Motivated by the impact of Antarctic and Greenland ice sheet change on future sea level rise, here we review the state of knowledge of ice sheet/Earth system interactions and feedbacks. We also describe emerging observation and model‐based methods that can improve understanding of ice sheet/Earth system interactions and feedbacks. We particularly focus on the development of Earth system models that incorporate current understanding of Earth system processes, ice dynamics, and ice sheet/Earth system couplings. Such models will be critical tools for projecting future sea level rise from anthropogenically forced ice sheet mass loss.

Past ice stream and ice sheet changes on the continental shelf off the Sabrina Coast, East Antarctica

Our understanding of the response of the Antarctic ice sheet to climate and ocean changes requires the improvement of current ice-atmosphere-ocean models and the accurate determination of boundary conditions such as ice thickness and extent at key time intervals, so that satellite gravity observations and isostatic models can be adjusted. However, large portions of the Antarctic margin remain understudied or lack suitable data. One key area where data are lacking, is the Sabrina Coast portion of the East Antarctic Ice Sheet (EAIS) margin where the Totten Glacier, which has the largest ice discharge in East Antarctica, is accelerating, thinning and loosing mass at high rates. In this work, we present the results of the first geological and geophysical marine survey to the continental shelf offshore of the Dalton Ice Tongue and Moscow University ice shelf, east of the Totten Glacier. The data presented include multibeam swath bathymetry and multichannel seismic, focusing on the sea floor morphology and sedimentary section above a regional angular unconformity separating pre- and post-Miocene glacial strata. Sea floor scouring and iceberg keel marks on the outer shelf, associated with gullies on the upper slope indicate that ice expanded in the past and grounded ~5 km from the shelf edge at ~450–500 mbsl, extending ~155 km north of the current Moscow University Ice Shelf. A nearly 1000 m deep area in the inner-middle shelf, oriented NW with paleo-ice flow direction indicated by mega scale glacial lineations (MSGL) and drumlins, is interpreted as a cross shelf glacial trough. A series of geomorphic associations on the north-eastern side of the glacial trough includes glacial lobes, grounding zone wedges (GZW), glacial lineations and transverse ridges, which indicates slower ice, grounding line stabilization and collapse. These geomorphic associations are organized in 4 four sets representing different past ice-flow configurations reflecting changes in ice flow direction, grounding line position, location of fast and slow ice areas, and retreat pattern. Some of the geomorphic features identified are compatible with the presence of an organized subglacial drainage, and others are with rapid grounding line collapse. A well-preserved series of GZWs occurring at different water depths implies they were formed during different glacial stages or cycles. The inferred diminishing ice thickness for consecutives GZWs indicates that the margin of the Antarctic ice sheet evolved to a less extensive coverage of the continental shelf through successive glacial stages or cycles.The identification of different ice flow configurations, evidence of subglacial water and past ice margin collapse indicates a dynamic ice sheet margin with varying glacial conditions and retreat modes. We observe that some of the described morphological associations are similar to those found in the Amundsen sea sector of the West Antarctic Ice Sheet (WAIS) where they are associated with ice sheet and ice stream collapse. Although further studies are needed to assess the precise timing and rates of the glacial processes involved, we conclude that there is enough evidence to support the hypothesis that the EAIS margin can behave as dynamically as the WAIS margin, especially during glacial retreat and icesheet margin collapse.

Spatial variation in Late Ordovician glacioeustatic sea-level change

Mass extinction of Late Ordovician marine fauna closely coincided with southern hemisphere glaciation. The sequence stratigraphic architecture of shallow marine deposits informs estimates of glacioeustatic sea-level change at sites both proximal and distal to the reconstructed Ordovician ice sheet(s) and contemporaneous changes in ice volume. A recent correlation framework for the stratigraphic architectures of one near and one far field Late Ordovician margin concluded that the Late Ordovician glaciation encompassed multiple long-term cycles of ice volume growth and retreat with superimposed higher frequency cycles. Here we posit that—similar to Cenozoic glacial cycles—glacial isostatic adjustment can preclude synchronous and similar magnitude (or directional) changes in Late Ordovician sea level between ice proximal and ice distal locations and, hence, distort a globally correlative sequence stratigraphy. We explored whether long-duration (i.e., million year) Late Ordovician glacial cycles should produce a globally coherent, eustatic record of sea-level change between ice proximal and ice distal margins using a gravitationally self-consistent theory that accounts for the deformational, gravitational and rotational perturbations to sea level on a viscoelastic Earth model. We adopted a Late Ordovician paleogeography and a synthetic continental ice-sheet distribution and volume informed by the areal extent of glaciogenic deposits and geochemical records, respectively. We demonstrate that modeled million year Late Ordovician glacial cycles produce sea-level histories on near and far field margins that differ from eustasy, and from one another, due primarily to elastic flexure and associated gravitational effects. While predicted far-field sea-level histories faithfully preserve the temporal structure of modeled glacioeustasy, their amplitude may differ from eustasy by as much as 30–40%. The impact of glacial isostatic adjustment is largest at the margins of glaciated continents, and these effects can be of the same order of magnitude as the eustatic, and even induce a local sea-level rise during an episode of ice growth and eustatic sea-level fall, and vice versa. In this regard, stratal surfaces of maximum regression and flooding expressed at near-field margins need not reflect global (‘eustatic’) trends in ice sheet growth and decay, respectively, and thus may not provide chronostratigraphic horizons for correlation with far-field sequence stratigraphic architectures.

Glacier extent in sub-Antarctic Kerguelen archipelago from MIS 3 period: Evidence from 36Cl dating

Documenting sub-Antarctic glacier variations during the local last glacial maximum is of major interest to better understand their sensitivity to atmospheric and oceanic temperature changes in conjunction with Antarctic ice sheet changes. However, data are sparse because evidence of earlier glacier extents is for most sub-Antarctic islands located offshore making their observation complex. Here, we present 22 cosmogenic 36Cl surface exposure ages obtained from five sites at Kerguelen to document the glacial history. The 36Cl ages from roche moutonnee surfaces, erratics and boulders collected on moraines span from 41.9 ± 4.4 ka to 14.3 ± 1.1 ka. Ice began to retreat on the eastern part of the main island before 41.4 ± 4.4 ka. Slow deglaciation occurred from ∼41 to ∼29 ka. There is no evidence of advances between 29 ka and the Antarctic Cold Reversal (ACR) period (∼14.5–12.9 ka) period. During the ACR, however, the Bontemps and possibly Belvedere moraines were formed by the advance of a Cook Ice Cap outlet glacier and a local glacier on the Presque Ile Jeanne d’Arc, respectively. This glacier evolution differs partly from that of glaciers in New Zealand and in Patagonia. These asynchronous glacier changes in the sub-Antarctic region are however in agreement with sea surface temperature changes recorded around Antarctica, which suggest differences in the climate evolution of the Indo-Pacific and Atlantic sectors of Antarctica.

Middle to late Holocene chronology of the western margin of the Greenland Ice Sheet: A comparison with Holocene temperature and precipitation records

ABSTRACT The response of the Greenland Ice Sheet to future climate change is relatively unconstrained. Determining the extents and rates of ice-margin fluctuations during the Holocene provides a longer-term perspective on ice-sheet changes and offers an analogue of how the ice-sheet may respond to future changes. Here, we present sixteen new 10Be ages of boulders on moraines, boulders perched on bedrock, and bedrock surfaces that mark the timing of ice-margin fluctuations during the Holocene in the Kangerlussuaq region of southern west Greenland. We show that the Keglen moraines date to 8.0 ± 0.3 ka (n = 6) and that the average ice-margin retreat rate slowed from about 49 to 13 m yr−1 after about 8.0 ka, likely in response to the ice margin retreating onto land at the head of the fjord Kangerlussuaq at this time. The average retreat rate further slowed to less than 1 m yr−1 between 6.8 ka and 4.2 cal kyr BP, a time when nearby paleoclimate records document warm summers and increased precipitation. Finally, we show that the historical advances of the ice margin occurred during the past 200 years, likely in response to cooler summer temperatures.