Ice Sheet Surface Research Articles

Direct measurements of meltwater runoff on the Greenland ice sheet surface

Meltwater runoff from the Greenland ice sheet surface influences surface mass balance (SMB), ice dynamics, and global sea level rise, but is estimated with climate models and thus difficult to validate. We present a way to measure ice surface runoff directly, from hourly in situ supraglacial river discharge measurements and simultaneous high-resolution satellite/drone remote sensing of upstream fluvial catchment area. A first 72-h trial for a 63.1-km2 moulin-terminating internally drained catchment (IDC) on Greenland's midelevation (1,207-1,381 m above sea level) ablation zone is compared with melt and runoff simulations from HIRHAM5, MAR3.6, RACMO2.3, MERRA-2, and SEB climate/SMB models. Current models cannot reproduce peak discharges or timing of runoff entering moulins but are improved using synthetic unit hydrograph (SUH) theory. Retroactive SUH applications to two older field studies reproduce their findings, signifying that remotely sensed IDC area, shape, and supraglacial river length are useful for predicting delays in peak runoff delivery to moulins. Applying SUH to HIRHAM5, MAR3.6, and RACMO2.3 gridded melt products for 799 surrounding IDCs suggests their terminal moulins receive lower peak discharges, less diurnal variability, and asynchronous runoff timing relative to climate/SMB model output alone. Conversely, large IDCs produce high moulin discharges, even at high elevations where melt rates are low. During this particular field experiment, models overestimated runoff by +21 to +58%, linked to overestimated surface ablation and possible meltwater retention in bare, porous, low-density ice. Direct measurements of ice surface runoff will improve climate/SMB models, and incorporating remotely sensed IDCs will aid coupling of SMB with ice dynamics and subglacial systems.

Blowing snow detection from ground-based ceilometers: application to East Antarctica

Abstract. Blowing snow impacts Antarctic ice sheet surface mass balance by snow redistribution and sublimation. However, numerical models poorly represent blowing snow processes, while direct observations are limited in space and time. Satellite retrieval of blowing snow is hindered by clouds and only the strongest events are considered. Here, we develop a blowing snow detection (BSD) algorithm for ground-based remote-sensing ceilometers in polar regions and apply it to ceilometers at Neumayer III and Princess Elisabeth (PE) stations, East Antarctica. The algorithm is able to detect (heavy) blowing snow layers reaching 30 m height. Results show that 78 % of the detected events are in agreement with visual observations at Neumayer III station. The BSD algorithm detects heavy blowing snow 36 % of the time at Neumayer (2011–2015) and 13 % at PE station (2010–2016). Blowing snow occurrence peaks during the austral winter and shows around 5 % interannual variability. The BSD algorithm is capable of detecting blowing snow both lifted from the ground and occurring during precipitation, which is an added value since results indicate that 92 % of the blowing snow is during synoptic events, often combined with precipitation. Analysis of atmospheric meteorological variables shows that blowing snow occurrence strongly depends on fresh snow availability in addition to wind speed. This finding challenges the commonly used parametrizations, where the threshold for snow particles to be lifted is a function of wind speed only. Blowing snow occurs predominantly during storms and overcast conditions, shortly after precipitation events, and can reach up to 1300 m a. g. l. in the case of heavy mixed events (precipitation and blowing snow together). These results suggest that synoptic conditions play an important role in generating blowing snow events and that fresh snow availability should be considered in determining the blowing snow onset.

Algae Drive Enhanced Darkening of Bare Ice on the Greenland Ice Sheet

AbstractSurface ablation of the Greenland ice sheet is amplified by surface darkening caused by light‐absorbing impurities such as mineral dust, black carbon, and pigmented microbial cells. We present the first quantitative assessment of the microbial contribution to the ice sheet surface darkening, based on field measurements of surface reflectance and concentrations of light‐absorbing impurities, including pigmented algae, during the 2014 melt season in the southwestern part of the ice sheet. The impact of algae on bare ice darkening in the study area was greater than that of nonalgal impurities and yielded a net albedo reduction of 0.038 ± 0.0035 for each algal population doubling. We argue that algal growth is a crucial control of bare ice darkening, and incorporating the algal darkening effect will improve mass balance and sea level projections of the Greenland ice sheet and ice masses elsewhere.

Basin-scale heterogeneity in Antarctic precipitation and its impact on surface mass variability

Abstract. Annually averaged precipitation in the form of snow, the dominant term of the Antarctic Ice Sheet surface mass balance, displays large spatial and temporal variability. Here we present an analysis of spatial patterns of regional Antarctic precipitation variability and their impact on integrated Antarctic surface mass balance variability simulated as part of a preindustrial 1800-year global, fully coupled Community Earth System Model simulation. Correlation and composite analyses based on this output allow for a robust exploration of Antarctic precipitation variability. We identify statistically significant relationships between precipitation patterns across Antarctica that are corroborated by climate reanalyses, regional modeling and ice core records. These patterns are driven by variability in large-scale atmospheric moisture transport, which itself is characterized by decadal- to centennial-scale oscillations around the long-term mean. We suggest that this heterogeneity in Antarctic precipitation variability has a dampening effect on overall Antarctic surface mass balance variability, with implications for regulation of Antarctic-sourced sea level variability, detection of an emergent anthropogenic signal in Antarctic mass trends and identification of Antarctic mass loss accelerations.

Evaluating Model Simulations of Twentieth-Century Sea-Level Rise. Part II: Regional Sea-Level Changes

Twentieth-century regional sea level changes are estimated from 12 climate models from phase 5 of the Climate Model Intercomparison Project (CMIP5). The output of the CMIP5 climate model simulations was used to calculate the global and regional sea level changes associated with dynamic sea level, atmospheric loading, glacier mass changes, and ice sheet surface mass balance contributions. The contribution from groundwater depletion, reservoir storage, and dynamic ice sheet mass changes are estimated from observations as they are not simulated by climate models. All contributions are summed, including the glacial isostatic adjustment (GIA) contribution, and compared to observational estimates from 27 tide gauge records over the twentieth century (1900–2015). A general agreement is found between the simulated sea level and tide gauge records in terms of interannual to multidecadal variability over 1900–2015. But climate models tend to systematically underestimate the observed sea level trends, particularly in the first half of the twentieth century. The corrections based on attributable biases between observations and models that have been identified in Part I of this two-part paper result in an improved explanation of the spatial variability in observed sea level trends by climate models. Climate models show that the spatial variability in sea level trends observed by tide gauge records is dominated by the GIA contribution and the steric contribution over 1900–2015. Climate models also show that it is important to include all contributions to sea level changes as they cause significant local deviations; note, for example, the groundwater depletion around India, which is responsible for the low twentieth-century sea level rise in the region.

Investigating the local-scale influence of sea ice on Greenland surface melt

Abstract. Rapid decline in Arctic sea ice cover in the 21st century may have wide-reaching effects on the Arctic climate system, including the Greenland ice sheet mass balance. Here, we investigate whether local changes in sea ice around the Greenland ice sheet have had an impact on Greenland surface melt. Specifically, we investigate the relationship between sea ice concentration, the timing of melt onset and open-water fraction surrounding Greenland with ice sheet surface melt using a combination of remote sensing observations, and outputs from a reanalysis model and a regional climate model for the period of 1979–2015. Statistical analysis points to covariability between Greenland ice sheet surface melt and sea ice within Baffin Bay and Davis Strait. While some of this covariance can be explained by simultaneous influence of atmospheric circulation anomalies on both the sea ice cover and Greenland melt, within Baffin Bay we find a modest correlation between detrended melt onset over sea ice and the adjacent ice sheet melt onset. This correlation appears to be related to increased transfer of sensible and latent heat fluxes from the ocean to the atmosphere in early sea ice melt years, increasing temperatures and humidity over the ice sheet that in turn initiate ice sheet melt.

Wave motion in a fluid under an inhomogeneous ice cover

This paper studies steady waves in fluid and in semi-infinite ice cover generated by a constant pressure distribution with a rectangular planform moving uniformly along the edge of ice cover at fixed distance. This load simulates the air-cushion vehicle (ACV). We consider two cases: (i) the surface of fluid is free outside of ice sheet, (ii) fluid is bounded by a solid vertical wall and the edge of ice cover can be either clamped or free. The fluid is assumed to be ideal incompressible and of finite depth. The ice sheet is modelled by elastic thin plate. The solution of linear hydroelastic problem is obtained by two methods: the Wiener-Hopf technique and matched eigenfunction expansions. The deflection of ice sheet and free surface elevation, as well as wave forces acting on ACV are investigated for different speeds of motion.

ANTARCTIC ICE SHEET SLOPE AND ASPECT BASED ON ICESAT’S REPEAT ORBIT MEASUREMENT

Abstract. Accurate information of ice sheet surface slope is essential for estimating elevation change by satellite altimetry measurement. A study is carried out to recover surface slope of Antarctic ice sheet from Ice, Cloud and land Elevation Satellite (ICESat) elevation measurements based on repeat orbits. ICESat provides repeat ground tracks within 200 meters in cross-track direction and 170 meters in along-track direction for most areas of Antarctic ice sheet. Both cross-track and along-track surface slopes could be obtained by adjacent repeat ground tracks. Combining those measurements yields a surface slope model with resolution of approximately 200 meters. An algorithm considering elevation change is developed to estimate the surface slope of Antarctic ice sheet. Three Antarctic Digital Elevation Models (DEMs) were used to calculate surface slopes. The surface slopes from DEMs are compared with estimates by using in situ GPS data in Dome A, the summit of Antarctic ice sheet. Our results reveal an average surface slope difference of 0.02 degree in Dome A. High resolution remote sensing images are also used in comparing the results derived from other DEMs and this paper. The comparison implies that our results have a slightly better coherence with GPS observation than results from DEMs, but our results provide more details and perform higher accuracy in coastal areas because of the higher resolution for ICESat measurements. Ice divides are estimated based on the aspect, and are weakly consistent with ice divides from other method in coastal regions.

On the importance of the albedo parameterization for the mass balance of the Greenland ice sheet in EC-Earth

Abstract. The albedo of the surface of ice sheets changes as a function of time due to the effects of deposition of new snow, ageing of dry snow, bare ice exposure, melting and run-off. Currently, the calculation of the albedo of ice sheets is highly parameterized within the earth system model EC-Earth by taking a constant value for areas with thick perennial snow cover. This is an important reason why the surface mass balance (SMB) of the Greenland ice sheet (GrIS) is poorly resolved in the model. The purpose of this study is to improve the SMB forcing of the GrIS by evaluating different parameter settings within a snow albedo scheme. By allowing ice-sheet albedo to vary as a function of wet and dry conditions, the spatial distribution of albedo and melt rate improves. Nevertheless, the spatial distribution of SMB in EC-Earth is not significantly improved. As a reason for this, we identify omissions in the current snow albedo scheme, such as separate treatment of snow and ice and the effect of refreezing. The resulting SMB is downscaled from the lower-resolution global climate model topography to the higher-resolution ice-sheet topography of the GrIS, such that the influence of these different SMB climatologies on the long-term evolution of the GrIS is tested by ice-sheet model simulations. From these ice-sheet simulations we conclude that an albedo scheme with a short response time of decaying albedo during wet conditions performs best with respect to long-term simulated ice-sheet volume. This results in an optimized albedo parameterization that can be used in future EC-Earth simulations with an interactive ice-sheet component.

Impact of MODIS sensor calibration updates on Greenland Ice Sheet surface reflectance and albedo trends

Abstract. We evaluate Greenland Ice Sheet (GrIS) surface reflectance and albedo trends using the newly released Collection 6 (C6) MODIS (Moderate Resolution Imaging Spectroradiometer) products over the period 2001–2016. We find that the correction of MODIS sensor degradation provided in the new C6 data products reduces the magnitude of the surface reflectance and albedo decline trends obtained from previous MODIS data (i.e., Collection 5, C5). Collection 5 and 6 data product analysis over GrIS is characterized by surface (i.e., wet vs. dry) and elevation (i.e., 500–2000 m, 2000 m and greater) conditions over the summer season from 1 June to 31 August. Notably, the visible-wavelength declining reflectance trends identified in several bands of MODIS C5 data from previous studies are only slightly detected at reduced magnitude in the C6 versions over the dry snow area. Declining albedo in the wet snow and ice area remains over the MODIS record in the C6 product, albeit at a lower magnitude than obtained using C5 data. Further analyses of C6 spectral reflectance trends show both reflectance increases and decreases in select bands and regions, suggesting that several competing processes are contributing to Greenland Ice Sheet albedo change. Investigators using MODIS data for other ocean, atmosphere and/or land analyses are urged to consider similar re-examinations of trends previously established using C5 data.

West Antarctic Ice Sheet Cloud Cover and Surface Radiation Budget from NASA A-Train Satellites

Clouds are an essential parameter of the surface energy budget influencing the West Antarctic Ice Sheet (WAIS) response to atmospheric warming and net contribution to global sea level rise. A 4-yr record of NASA A-Train cloud observations is combined with surface radiation measurements to quantify the WAIS radiation budget and constrain the three-dimensional occurrence frequency, thermodynamic phase partitioning, and surface radiative effect of clouds over West Antarctica (WA). The skill of satellite-modeled radiative fluxes is confirmed through evaluation against measurements at four Antarctic sites (WAIS Divide ice camp and Neumayer, Syowa, and Concordia stations). Owing to perennial high-albedo snow and ice cover, cloud infrared emission dominates over cloud solar reflection and absorption leading to a positive net all-wave cloud radiative effect (CRE) at the surface, with all monthly means and 99.15% of instantaneous CRE values exceeding zero. The annual-mean CRE at the WAIS surface is 34 W m−2, representing a significant cloud-induced warming of the ice sheet. Low-level liquid-containing clouds, including thin liquid water clouds implicated in radiative contributions to surface melting, are widespread and most frequent in WA during the austral summer. In summer, clouds warm the WAIS by 26 W m−2, on average, despite maximum offsetting shortwave CRE. Glaciated cloud systems are strongly linked to orographic forcing, with maximum incidence on the WAIS continuing downstream along the Transantarctic Mountains.

Antarctic climate and ice-sheet configuration during the early Pliocene interglacial at 4.23 Ma

Abstract. The geometry of Antarctic ice sheets during warm periods of the geological past is difficult to determine from geological evidence, but is important to know because such reconstructions enable a more complete understanding of how the ice-sheet system responds to changes in climate. Here we investigate how Antarctica evolved under orbital and greenhouse gas conditions representative of an interglacial in the early Pliocene at 4.23 Ma, when Southern Hemisphere insolation reached a maximum. Using offline-coupled climate and ice-sheet models, together with a new synthesis of high-latitude palaeoenvironmental proxy data to define a likely climate envelope, we simulate a range of ice-sheet geometries and calculate their likely contribution to sea level. In addition, we use these simulations to investigate the processes by which the West and East Antarctic ice sheets respond to environmental forcings and the timescales over which these behaviours manifest. We conclude that the Antarctic ice sheet contributed 8.6 ± 2.8 m to global sea level at this time, under an atmospheric CO2 concentration identical to present (400 ppm). Warmer-than-present ocean temperatures led to the collapse of West Antarctica over centuries, whereas higher air temperatures initiated surface melting in parts of East Antarctica that over one to two millennia led to lowering of the ice-sheet surface, flotation of grounded margins in some areas, and retreat of the ice sheet into the Wilkes Subglacial Basin. The results show that regional variations in climate, ice-sheet geometry, and topography produce long-term sea-level contributions that are non-linear with respect to the applied forcings, and which under certain conditions exhibit threshold behaviour associated with behavioural tipping points.

Ecological Modeling of the Supraglacial Ecosystem: A Process-based Perspective

Glacier and ice sheet surfaces are important microbe-dominated ecosystems that are changing rapidly due to climate change, with potentially significant impacts. A theoretical framework of the supraglacial (glacier surface) ecosystem is needed to enable its mathematical modelling, a necessary tool for understanding, quantifying and predicting present day and future ecosystem dynamics. Here we review key biological processes occurring on glacier and ice sheet surfaces and present three frameworks for constructing process-based models of the surface ecosystem, using the largest supraglacial ecosystem on Earth – the Greenland ice sheet surface – as an important example. The models are based on organic carbon transformations, but vary in numerical complexity and in the level of detail of biological processes. This perspective is intended to guide future supraglacial ecosystem model development, field data collection for parameterisation and validation purposes, and encourage inter-disciplinary collaboration between modellers and experimentalists.

Hypsometric amplification and routing moderation of Greenland ice sheet meltwater release

Abstract. Concurrent ice sheet surface runoff and proglacial discharge monitoring are essential for understanding Greenland ice sheet meltwater release. We use an updated, well-constrained river discharge time series from the Watson River in southwest Greenland, with an accurate, observation-based ice sheet surface mass balance model of the ∼ 12 000 km2 ice sheet area feeding the river. For the 2006–2015 decade, we find a large range of a factor of 3 in interannual variability in discharge. The amount of discharge is amplified ∼ 56 % by the ice sheet's hypsometry, i.e., area increase with elevation. A good match between river discharge and ice sheet surface meltwater production is found after introducing elevation-dependent transit delays that moderate diurnal variability in meltwater release by a factor of 10–20. The routing lag time increases with ice sheet elevation and attains values in excess of 1 week for the upper reaches of the runoff area at ∼ 1800 m above sea level. These multi-day routing delays ensure that the highest proglacial discharge levels and thus overbank flooding events are more likely to occur after multi-day melt episodes. Finally, for the Watson River ice sheet catchment, we find no evidence of meltwater storage in or release from the en- and subglacial environments in quantities exceeding our methodological uncertainty, based on the good match between ice sheet runoff and proglacial discharge.

Decoding ice sheet behavior using englacial layer slopes

AbstractThe complex flow fields of the Antarctic and Greenland ice sheets deform layers deposited as snow at the ice sheet surface, leaving a record of the regional flow history and/or local transitions in basal boundary conditions within the geometry of ice sheet layers. Ice‐penetrating radar reveals these layers, but radar data interpretations are limited by the challenges of quantitatively and reproducibly comparing observations with model output. We present a conceptual framework that relates along‐track reflector slope to gradients in the steady state velocity field of ice sheets. This method makes effective use of englacial reflectors in regions where it is challenging to image continuous layers and avoids the error propagation inherent to tracer‐transport methods, developing the potential for formal radar data assimilation in future modeling studies. We apply our method to radar data collected at the grounding line of Whillans Ice Stream, where enhanced bed friction produces characteristic reflector slopes reproducible using a higher‐order ice flow model.

A new approach to estimate ice dynamic rates using satellite observations in East Antarctica

Abstract. Mass balance changes of the Antarctic ice sheet are of significant interest due to its sensitivity to climatic changes and its contribution to changes in global sea level. While regional climate models successfully estimate mass input due to snowfall, it remains difficult to estimate the amount of mass loss due to ice dynamic processes. It has often been assumed that changes in ice dynamic rates only need to be considered when assessing long-term ice sheet mass balance; however, 2 decades of satellite altimetry observations reveal that the Antarctic ice sheet changes unexpectedly and much more dynamically than previously expected. Despite available estimates on ice dynamic rates obtained from radar altimetry, information about ice sheet changes due to changes in the ice dynamics are still limited, especially in East Antarctica. Without understanding ice dynamic rates, it is not possible to properly assess changes in ice sheet mass balance and surface elevation or to develop ice sheet models. In this study we investigate the possibility of estimating ice sheet changes due to ice dynamic rates by removing modelled rates of surface mass balance, firn compaction, and bedrock uplift from satellite altimetry and gravity observations. With similar rates of ice discharge acquired from two different satellite missions we show that it is possible to obtain an approximation of the rate of change due to ice dynamics by combining altimetry and gravity observations. Thus, surface elevation changes due to surface mass balance, firn compaction, and ice dynamic rates can be modelled and correlated with observed elevation changes from satellite altimetry.

Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model

Abstract. With the aim of studying the recent Greenland ice sheet (GrIS) surface mass balance (SMB) decrease relative to the last century, we have forced the regional climate MAR (Modèle Atmosphérique Régional; version 3.5.2) model with the ERA-Interim (ECMWF Interim Re-Analysis; 1979–2015), ERA-40 (1958–2001), NCEP–NCARv1 (National Centers for Environmental Prediction–National Center for Atmospheric Research Reanalysis version 1; 1948–2015), NCEP–NCARv2 (1979–2015), JRA-55 (Japanese 55-year Reanalysis; 1958–2014), 20CRv2(c) (Twentieth Century Reanalysis version 2; 1900–2014) and ERA-20C (1900–2010) reanalyses. While all these forcing products are reanalyses that are assumed to represent the same climate, they produce significant differences in the MAR-simulated SMB over their common period. A temperature adjustment of +1 °C (respectively −1 °C) was, for example, needed at the MAR boundaries with ERA-20C (20CRv2) reanalysis, given that ERA-20C (20CRv2) is ∼ 1 °C colder (warmer) than ERA-Interim over Greenland during the period 1980–2010. Comparisons with daily PROMICE (Programme for Monitoring of the Greenland Ice Sheet) near-surface observations support these adjustments. Comparisons with SMB measurements, ice cores and satellite-derived melt extent reveal the most accurate forcing datasets for the simulation of the GrIS SMB to be ERA-Interim and NCEP–NCARv1. However, some biases remain in MAR, suggesting that some improvements are still needed in its cloudiness and radiative schemes as well as in the representation of the bare ice albedo. Results from all MAR simulations indicate that (i) the period 1961–1990, commonly chosen as a stable reference period for Greenland SMB and ice dynamics, is actually a period of anomalously positive SMB (∼ +40 Gt yr−1) compared to 1900–2010; (ii) SMB has decreased significantly after this reference period due to increasing and unprecedented melt reaching the highest rates in the 120-year common period; (iii) before 1960, both ERA-20C and 20CRv2-forced MAR simulations suggest a significant precipitation increase over 1900–1950, but this increase could be the result of an artefact in the reanalyses that are not well-enough constrained by observations during this period and (iv) since the 1980s, snowfall is quite stable after having reached a maximum in the 1970s. These MAR-based SMB and accumulation reconstructions are, however, quite similar to those from Box (2013) after 1930 and confirm that SMB was quite stable from the 1940s to the 1990s. Finally, only the ERA-20C-forced simulation suggests that SMB during the 1920–1930 warm period over Greenland was comparable to the SMB of the 2000s, due to both higher melt and lower precipitation than normal.

The million-year evolution of the glacial trimline in the southernmost Ellsworth Mountains, Antarctica

An elevated erosional trimline in the heart of West Antarctica in the Ellsworth Mountains tells of thicker ice in the past and represents an important yet ambiguous stage in the evolution of the Antarctic Ice Sheet. Here we analyse the geomorphology of massifs in the southernmost Heritage Range where the surfaces associated with the trimline are overlain by surficial deposits that have the potential to be dated through cosmogenic nuclide analysis. Analysis of 100 rock samples reveals that some clasts have been exposed on glacially moulded surfaces for 1.4 Ma and perhaps more than 3.5 Ma, while others reflect fluctuations in thickness during Quaternary glacial cycles. Modelling the age of the glacially moulded bedrock surface based on cosmogenic 10Be, 26Al and 21Ne concentrations from a depth-profile indicates a minimum exposure age of 2.1–2.6 Ma. We conclude that the glacially eroded surfaces adjacent to the trimline predate the Last Glacial Maximum and indeed the Quaternary. Since erosion was by warm-based ice near an ice-sheet upper margin, we suggest it first occurred during the early glaciations of Antarctica before the stepped cooling of the mid-Miocene at ∼14 Ma. This was a time when the interior Antarctic continent had summers warm enough for tundra vegetation to grow and for mountain glaciers to consist of ice at the pressure melting point. During these milder conditions, and subsequently, erosion of glacial troughs is likely to have lowered the ice-sheet surface in relation to the mountains. This means that the range of orbitally induced cyclic fluctuations in ice thickness have progressively been confined to lower elevations.

Rapid thinning of the Laurentide Ice Sheet in coastal Maine, USA, during late Heinrich Stadial 1

Few data are available to infer the thinning rate of the Laurentide Ice Sheet (LIS) through the last deglaciation, despite its importance for constraining past ice sheet response to climate warming. We measured 31 cosmogenic 10Be exposure ages in samples collected on coastal mountainsides in Acadia National Park and from the slightly inland Pineo Ridge moraine complex, a ∼100-km-long glaciomarine delta, to constrain the timing and rate of LIS thinning and subsequent retreat in coastal Maine. Samples collected along vertical transects in Acadia National Park have indistinguishable exposure ages over a 300 m range of elevation, suggesting that rapid, century-scale thinning occurred at 15.2 ± 0.7 ka, similar to the timing of abrupt thinning inferred from cosmogenic exposure ages at Mt. Katahdin in central Maine (Davis et al., 2015). This rapid ice sheet surface lowering, which likely occurred during the latter part of the cold Heinrich Stadial 1 event (19–14.6 ka), may have been due to enhanced ice-shelf melt and calving in the Gulf of Maine, perhaps related to regional oceanic warming associated with a weakened Atlantic Meridional Overturning Circulation at this time. The ice margin subsequently stabilized at the Pineo Ridge moraine complex until 14.5 ± 0.7 ka, near the onset of Bølling Interstadial warming. Our 10Be ages are substantially younger than marine radiocarbon constraints on LIS retreat in the coastal lowlands, suggesting that the deglacial marine reservoir effect in this area was ∼1,200 14C years, perhaps also related to the sluggish Atlantic Meridional Overturning Circulation during Heinrich Stadial 1.

Mapping ice sheet grounding lines with CryoSat-2

The boundary between grounded and floating ice is an important glaciological parameter, because it delineates the lateral extent of an ice sheet and it marks the optimal location for computing ice discharge. We present a method for detecting the grounding line as the break in ice sheet surface slope, computed from CryoSat-2 elevation measurements using a plane-fitting solution. We apply this technique to map the break in surface slope in four topographically diverse sectors of Antarctica - Filchner-Ronne Ice Shelf, Ekström Ice Shelf, Amundsen Sea sector, and the Larsen-C Ice Shelf - using CryoSat-2 observations acquired between July 2010 and May 2014. An inter-comparison of the CryoSat-2 break in surface slope with independent measurements of the hinge line determined from quadruple-difference SAR interferometry (QDInSAR) shows good overall agreement between techniques, with a mean separation of 4.5km. In the Amundsen Sea sector, where in places over 35km of hinge line retreat has occurred since 1992, the CryoSat-2 break in surface slope coincides with the most recent hinge line position, recorded in 2011. The technique we have developed is automatic, is computationally-efficient, can be repeated given further data, and offers a complementary tool for monitoring changes in the lateral extent of grounded ice.