Ice Surface Elevation Research Articles

Validation and Analysis of the ICESat-2 ATL11 Product: A Case Study of Lake Vostok

Elevation changes are crucial input data for mass balance assessment based on satellite altimetry. However, due to the accuracy limitations of altimetry data and mass conversion models, a significant uncertainty remains in estimating the mass balance of the Antarctic ice sheet, especially in East Antarctica. The ICESat-2 photon altimetry satellite enhances the precision of ice sheet surface elevation measurements to an accuracy of 2–4 cm. This improvement is particularly beneficial for detecting subtle elevation changes in East Antarctica. The ATL11 product from ICESat-2 offers a time series of ice surface elevations across Antarctica, enabling direct calculation of elevation change rates. To evaluate the capability of the ATL11 product in accurately depicting detailed elevation changes within the local terrain, we selected the Vostok Subglacial Lake as the validation region. Our research involves comparing fitted surface elevation change rates with in-situ data, while also considering surface mass balance, wind, and surface elevation to analyze the factors contributing to small differences in elevation changes within the local lake area. This analysis aims to identify the factors contributing to the minor variations in elevation changes within the local lake area. According to our analysis, the Vostok Lake surface elevation change rate is 2.00 ± 0.77 cm yr−1 from April 2019 to June 2023, with an average period of 356 ± 81 days. The results demonstrate that the ATL11 elevation product sequence has the potential to accurately characterize subtle elevation changes and seasonal variations in the Antarctic ice sheet.

Holocene thinning in central Greenland controlled by the Northeast Greenland Ice Stream

Ice-core records from the interior of the Greenland ice sheet suggest widespread thinning during the Holocene. However, the recurring underestimation of this thinning in numerical models raises concerns about both the veracity of such reconstructions and the reliability of glaciological models. Recent work suggests the 8000-year-old Northeast Greenland Ice Stream (NEGIS), including a now-extinct northern tributary, may have been an early influence on Greenland ice-sheet dynamics. Yet, the inaccurate reproduction of NEGIS-like dynamics in most models hampers investigation of whether this feature played a role in Holocene ice-sheet thinning. Here we show that grounding-line retreat in northeast Greenland triggers elevation changes at the northern summit via ice-dynamic effects modulated by the paleo NEGIS system. In our simulations, fast ice-stream flow caused by transiently imposed reduced basal shear stress following the northeast retreat explains 55% ( ± 18%) of the estimated ice thinning, showing that ice-stream dynamics is one of the main drivers of the NGRIP Holocene surface elevation drop. Our findings show that the ice-flow in northeast Greenland plays a large role in ice-surface elevation changes in central Greenland.

Early aerial expedition photos reveal 85 years of glacier growth and stability in East Antarctica

During the last few decades, several sectors in Antarctica have transitioned from glacial mass balance equilibrium to mass loss. In order to determine if recent trends exceed the scale of natural variability, long-term observations are vital. Here we explore the earliest, large-scale, aerial image archive of Antarctica to provide a unique record of 21 outlet glaciers along the coastline of East Antarctica since the 1930s. In Lützow-Holm Bay, our results reveal constant ice surface elevations since the 1930s, and indications of a weakening of local land-fast sea-ice conditions. Along the coastline of Kemp and Mac Robertson, and Ingrid Christensen Coast, we observe a long-term moderate thickening of the glaciers since 1937 and 1960 with periodic thinning and decadal variability. In all regions, the long-term changes in ice thickness correspond with the trends in snowfall since 1940. Our results demonstrate that the stability and growth in ice elevations observed in terrestrial basins over the past few decades are part of a trend spanning at least a century, and highlight the importance of understanding long-term changes when interpreting current dynamics.

A comparative study of mass balance in the Lambert Glacier and Amery Ice Shelf along the Chinese inland traverse during 2019–2023 using altimetry, gravity, and in-situ observations

Abstract. The new photon-counting laser altimetry satellite ICESat-2 was successfully launched on September 15, 2018 with an unprecedented ice surface elevation measurement accuracy of 2–4 cm. It is meaningful for accurately estimating volumetric changes in the Antarctic Ice Sheet. Cross-validation of different types of data, especially comparison with in-situ data, is important for the ice sheet mass balance results of a new satellite. This paper proposes an elevation-difference method of grid elevation change rate model, based on ICESat-2 ATL11 elevation time series data, to distinguish the linear change trend of ice sheet surface elevation from the elevation change resulting from periodic precipitation. In order to compare the estimated elevation-change rate with the flux computed from in-situ snow stake velocity measurements and GRACE-FO gravity survey data, we made corrections for firn air content, elastic, and glacier isostatic adjustment. Based on the ATL11 data from 2019 to 2023, our results show that the ice sheet change in Basin 11 along the CHINARE traverse is from ~0.019 m yr−1 to ~0.121 m yr−1, and the mass balance in the upstream of the traverse in Basin 11 is ~1.9±0.2 Gt yr−1. It is comparable to ~2.4±1.2 Gt yr−1 from GRACE-FO during the same time period. Furthermore, the flux across traverse-11.9±1.1 Gt yr−1 is comparable to that of ~-9.7±0.9 Gt yr−1 across the same flux gate during 1997–2009 which is calculated based on GNSS-derived ice velocity observations, considering the time period difference and uncertainties.

Terminus thinning drives recent acceleration of a Greenlandic lake-terminating outlet glacier

Abstract Ice-contact proglacial lakes affect ice dynamics and the transition of glacier termini from land- to lake-terminating has been shown to cause ice flow acceleration. In recent decades, the number and size of Greenlandic ice-marginal lakes has increased, highlighting the need to further understand these lake-terminating ice-margins as their influence on ice sheet mass balance increases. Here, time series of satellite-derived observations of ice velocity, surface elevation, and terminus position were generated at a lake-terminating outlet glacier, Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia in south-west Greenland. At Isortuarsuup Sermia, annual surface velocity at the terminus increased by a factor of 2.5 to 214 ± 4 m yr−1 (2013–2021), with the magnitude of this acceleration declining with distance up-glacier. Meanwhile, near-terminus surface elevation changed at a rate of −2.3 ± 1.1 m yr−1 (2012–2021). Conversely, velocity change at Kangaasarsuup Sermia was minimal, while surface elevation change was approximately half at comparable elevations (−1.2 ± 0.3 m yr−1). We attribute these dynamic differences to thinning at Isortuarsuup Sermia and subsequent retreat from a stabilising sublacustrine moraine, and emphasise the potential of proglacial lakes to enhance future rates of mass loss from the Greenland Ice Sheet.

The complex basal morphology and ice dynamics of the Nansen Ice Shelf, East Antarctica

Abstract. Ice shelf dynamics and morphology play an important role in the stability of floating bodies of ice by driving fracturing that can lead to calving, in turn impacting the ability of the ice shelf to buttress upstream grounded ice. Following a 2016 calving event at the Nansen Ice Shelf (NIS), East Antarctica, we collected airborne and ground-based radar data to map ice thickness across the shelf. We combine these data with published satellite-derived data to examine the spatial variations in ice shelf draft, the cause and effects of ice shelf strain rates, and the possibility that a suture zone may be channelizing ocean water and altering patterns of sub-ice-shelf melt and freeze-on. We also use our datasets to assess limitations that may arise from relying on hydrostatic-balance equations applied to ice surface elevation to determine ice draft morphology. We find that the Nansen Ice Shelf has a highly variable basal morphology driven primarily by the formation of basal fractures near the onset of the ice shelf suture zone. This morphology is reflected in the ice shelf strain rates but not in the calculated hydrostatic-balance thickness, which underestimates the scale of variability at the ice shelf base. Enhanced melt rates near the ice shelf terminus and in steep regions of the channelized suture zone, along with relatively thin ice in the suture zone, appear to represent vulnerable areas in the NIS. This morphology, combined with ice dynamics, induce strain that has led to the formation of transverse fractures within the suture zone, resulting in large-scale calving events. Similar transverse fractures at other Antarctic ice shelves may also be driven by highly variable morphology, and predicting their formation and evolution could aid projections of ice shelf stability.

Mass changes of the northern Antarctic Peninsula Ice Sheet derived from repeat bi-static synthetic aperture radar acquisitions for the period 2013–2017

Abstract. Some of the highest specific mass change rates in Antarctica are reported for the Antarctic Peninsula. However, the existing estimates for the northern Antarctic Peninsula (<70∘ S) are either spatially limited or are affected by considerable uncertainties. The complex topography, frequent cloud cover, limitations in ice thickness information, boundary effects, and uncertain glacial–isostatic adjustment estimates affect the ice sheet mass change estimates using altimetry, gravimetry, or the input-output method. Within this study, the first assessment of the geodetic mass balance throughout the ice sheet of the northern Antarctic Peninsula is carried out employing bi-static synthetic aperture radar (SAR) data from the TanDEM-X satellite mission. Repeat coverages from the austral winters of 2013 and 2017 are employed. Overall, coverage of 96.4 % of the study area by surface elevation change measurements and a total mass budget of -24.1±2.8 Gt a−1 are revealed. The spatial distribution of the surface elevation and mass changes points out that the former ice shelf tributary glaciers of the Prince Gustav Channel, Larsen A and B, and Wordie ice shelves are the hotspots of ice loss in the study area and highlights the long-lasting dynamic glacier adjustments after the ice shelf break-up events. The highest mass change rate is revealed for the Airy–Seller–Fleming glacier system at -4.9±0.6 Gt a−1, and the highest average surface elevation change rate of -2.30±0.03 m a−1 is observed at Drygalski Glacier. The comparison of the ice mass budget with anomalies in the climatic mass balance indicates, that for wide parts of the southern section of the study area, the mass changes can be partly attributed to changes in the climatic mass balance. However, imbalanced high ice discharge drives the overall ice loss. The previously reported connection between mid-ocean warming along the southern section of the west coast and increased frontal glacier recession does not repeat in the pattern of the observed glacier mass losses, excluding in Wordie Bay. The obtained results provide information on ice surface elevation and mass changes for the entire northern Antarctic Peninsula on unprecedented spatially detailed scales and with high precision and will be beneficial for subsequent analysis and modeling.

Extensive and anomalous grounding line retreat at Vanderford Glacier, Vincennes Bay, Wilkes Land, East Antarctica

Abstract. Wilkes Land, East Antarctica, has been losing mass at an accelerating rate over recent decades in response to enhanced oceanic forcing. Overlying the Aurora Subglacial Basin, it has been referred to as the “weak underbelly” of the East Antarctic Ice Sheet and is drained by several major outlet glaciers. Despite their potential importance, few of these glaciers have been studied in detail. This includes the six outlet glaciers which drain into Vincennes Bay, a region recently discovered to have the warmest intrusions of modified Circumpolar Deep Water (mCDW) ever recorded in East Antarctica. Here, we use satellite imagery; differential synthetic aperture radar interferometry (DInSAR); and remotely sensed datasets of ice-surface velocity, ice-surface elevation and grounding line position to investigate ice dynamics between 1963 and 2022. Our results support previous observations of extensive grounding line retreat at Vanderford Glacier, measured at 18.6 km between 1996 and 2020. The persistent grounding line retreat, averaging 0.8 km yr−1, places Vanderford Glacier as the fastest retreating glacier in East Antarctica, and the third fastest in Antarctica, across decadal timescales. Such rapid retreat is consistent with the hypothesis that warm mCDW is able to access deep cavities formed below the Vanderford Ice Shelf, driving high rates of basal melting close to the grounding line. With a retrograde slope observed inland along the Vanderford Trench, such oceanic forcing may have significant implications for the future stability of Vanderford Glacier.

Subglacial hydrology from high-resolution ice-flow simulations of the Rhine Glacier during the Last Glacial Maximum: a proxy for glacial erosion

Abstract. At the Last Glacial Maximum (LGM), the Rhine Glacier complex (Rhine and Linth glaciers) formed large piedmont lobes extending north into the Swiss and German Alpine forelands. Numerous overdeepened valleys there were formed by repeated glaciations. A characteristic of these overdeepened valleys is their location close to the LGM ice margin, away from the Alps. Numerical models of ice flow of the Rhine Glacier indicate a poor fit between the sliding distance, a proxy for glacial erosion, and the location of these overdeepenings. Calculations of the hydraulic potential based on the computed time-dependent ice surface elevations of the Rhine Glacier lobe obtained from a high-resolution thermo-mechanically coupled Stokes flow model are used to estimate the location of subglacial water drainage routes. Results indicate that the subglacial water discharge is high and focused along glacial valleys and overdeepenings when water pressure is equal to the ice overburden pressure. These conditions are necessary for subglacial water to remove basal sediments, expose fresh bedrock, and favor further erosion by quarrying and abrasion. Knowledge of the location of paleo-subglacial water drainage routes may be useful to understand patterns of subglacial erosion beneath paleo-ice masses that do not otherwise relate to the sliding of ice. Comparison of the erosion pattern from subglacial meltwater with those from quarrying and abrasion shows the importance of subglacial water flow in the formation of distal overdeepenings in the Swiss lowlands.

Subglacial lake activity beneath the ablation zone of the Greenland Ice Sheet

Abstract. Hydrologically active subglacial lakes can drain large volumes of water and sediment along subglacial pathways, affecting the motion and mass balance of ice masses and impacting downstream sediment dynamics. To date, only eight active lakes have been reported beneath the Greenland Ice Sheet (GrIS), and thus the understanding of their spatial distribution and dynamic processes is still lacking. Here, using ICESat-2 (Ice, Cloud, and land Elevation Satellite-2) ATL11 data, we identify 18 active subglacial lakes, 16 of which have not been previously reported. Multi-temporal ArcticDEM (digital elevation model of the Arctic) strip maps were used to extend the time series to verify lakes and determine their drainage history. The identification of active subglacial lakes beneath the GrIS is complicated by the occurrence of supraglacial lakes, which also fill and drain and are hypothesized to be almost co-located. We therefore used the temporal pattern of ice-surface elevation change to discriminate subglacial lakes and utilized the ability of ICESat-2 to penetrate through surface water to correct the elevation provided by the ATL11 data. A significant localized elevation anomaly (−16.03–10.30 m yr−1) was measured in all detected subglacial lakes after correction, revealing that six subglacial lakes are twinned with supraglacial lakes. The active subglacial lakes have large upstream hydrological catchments and are located near or below the equilibrium line. Lakes have a median area of 1.20 km2, and 12 lakes exhibited positive elevation-change rates during the ICESat-2 period. These observations illustrate the potential for combining ICESat-2 and the ArcticDEM to differentiate small subglacial lakes in the ablation zone and beneath supraglacial lakes.

Topography reconstruction and evolution analysis of outlet glacier using data from unmanned aerial vehicles in Antarctica

Reconstruction of glacial topography is important for assessing the ice dynamics of glaciers in the past and understanding how they may respond to climate change in the future. As an emerging strategy, unmanned aerial vehicles (UAVs) have been successfully used in glaciology applications to reconstruct surface topography and monitor the short-term dynamics of glaciers. However, none of these studies have focused on the ice dynamics of outlet glaciers in Antarctica. In this study, based on a combination of UAVs and a base station, we investigated Dalk Glacier, a typical marine-terminating glacier in East Antarctica, during two Chinese National Antarctic Research Expeditions from 2019 to 2020. By applying structure-from-motion and multi-view-stereo photogrammetry, high-resolution orthomosaics and digital elevation models of the glacial topography were reconstructed with centimeter-level accuracy via in situ validation, thus marking the first application of UAV observations in the monitoring of Antarctic outlet glaciers. Topographic evolution of the glacier was quantitatively analyzed from various aspects including ice velocity, surface elevation, crevasses, and ice front calving. A maximum ice velocity of ∼ 310 m a-1 at its terminus was observed, which was ∼ 90 m a-1 greater than that in satellite-based studies. We analyzed in detail the spatially heterogeneous changes in the ice velocity and surface elevation of the glacier under the influence of an ice rumple at its calving terminus. Combined with satellite images and existing datasets, we hypothesize that the large ice rumples at the glacier terminus could deform the glacier and potentially damage its structural integrity, thereby limiting the growth of its terminus.

Assessing UAV-based laser scanning for monitoring glacial processes and interactions at high spatial and temporal resolutions

Uncrewed Aerial Vehicles (UAVs), in combination with Structure from Motion (SfM) photogrammetry, have become an established tool for reconstructing glacial and ice-marginal topography, yet the method is highly dependent on several factors, all of which can be highly variable in glacial environments. However, recent technological advancements, related primarily to the miniaturisation of new payloads such as compact Laser Scanners (LS), has provided potential new opportunities for cryospheric investigation. Indeed, UAV-LS systems have shown promise in forestry, river, and snow depth research, but to date the method has yet to be deployed in glacial settings. As such, in this study we assessed the suitability of UAV-LS for glacial research by investigating short-term changes in ice surface elevation, calving front geometry and crevasse morphology over the near-terminus region of an actively calving glacier in southeast Iceland. We undertook repeat surveys over a 0.1 km2region of the glacier at sub-daily, daily, and weekly temporal intervals, producing directly georeferenced point clouds at very high spatial resolutions (average of >300 points per m−2at 40 m flying height). Our data has enabled us to: 1) Accurately map surface elevation changes (Median errors under 0.1 m), 2) Reconstruct the geometry and evolution of an active calving front, 3) Produce more accurate estimates of the volume of ice lost through calving, and 4) Better detect surface crevasse morphology, providing future scope to extract size, depth and improve the monitoring of their evolution through time. We also compared our results to data obtained in parallel using UAV-SfM, which further emphasised the relative advantages of our method and suitability in glaciology. Consequently, our study highlights the potential of UAV-LS in glacial research, particularly for investigating glacier mass balance, changing ice dynamics, and calving glacier behaviour, and thus we suggest it has a significant role in advancing our knowledge of, and ability to monitor, rapidly changing glacial environments in future.

Simulating surface height and terminus position for marine outlet glaciers using a level set method with data assimilation

We implement a data assimilation framework for integrating ice surface and terminus position observations into a numerical ice-flow model. The model uses the well-known shallow shelf approximation (SSA) coupled to a level set method to capture ice motion and changes in the glacier geometry. The level set method explicitly tracks the evolving ice–atmosphere and ice–ocean boundaries for a marine outlet glacier. We use an Ensemble Transform Kalman Filter to assimilate observations of ice surface elevation and lateral ice extent by updating the level set function that describes the ice interface. Numerical experiments on an idealized marine-terminating glacier demonstrate the effectiveness of our data assimilation approach for tracking seasonal and multi-year glacier advance and retreat cycles. The model is also applied to simulate Helheim Glacier, a major tidewater-terminating glacier of the Greenland Ice Sheet that has experienced a recent history of rapid retreat. By assimilating observations from remotely-sensed surface elevation profiles we are able to more accurately track the migrating glacier terminus and glacier surface changes. These results support the use of data assimilation methodologies for obtaining more accurate predictions of short-term ice sheet dynamics.

Inverting ice surface elevation and velocity for bed topography and slipperiness beneath Thwaites Glacier

Abstract. There is significant uncertainty over how ice sheets and glaciers will respond to rising global temperatures. Limited knowledge of the topography and rheology of the ice–bed interface is a key cause of this uncertainty as models show that small changes in the bed can have a large influence on predicted rates of ice loss. Most of our detailed knowledge of bed topography comes from airborne and ground-penetrating radar observations. However, these direct observations are not spaced closely enough to meet the requirements of ice-sheet models, so interpolation and inversion methods are used to fill in the gaps. Here we present the results of a new inversion of surface elevation and velocity data over Thwaites Glacier, West Antarctica, for bed topography and slipperiness (i.e. the degree of basal slip for a given level of drag). The inversion is based on a steady-state linear perturbation analysis of the shallow-ice-stream equations. The method works by identifying disturbances to surface flow which are caused by obstacles or sticky patches in the bed and can therefore be applied wherever the shallow-ice-stream equations hold and where surface data are available, even where the ice thickness is not well known. We assess the performance of the inversion for topography with the available radar data. Although the topographic output from the inversion is less successful where the bed slopes steeply, it compares well with radar data from the central trunk of the glacier for medium-wavelength features (5–50 km). This method could therefore be useful as an independent test of other interpolation methods such as mass conservation and kriging. We do not have data to allow us to assess the success of the slipperiness results from our inversions, but we provide maps that may guide future seismic data collection across Thwaites Glacier. The methods presented here show significant promise for using high-resolution satellite datasets, calibrated by sparser field datasets, to generate high-resolution bed topography products across the ice sheets and therefore contribute to reduced uncertainty in predictions of future sea-level rise.

Multi-decadal elevation changes of the land terminating sector of West Greenland

AbstractRegional assessments of ice elevation change provide insight into the processes controlling an ice sheet's geometric response to climate forcing. In Southwest Greenland's land terminating sector (SWLTS), it is presumed that ice surface elevation changes result solely from changing surface mass balance (SMB). Here we test this assumption by developing a multi-decadal (1985–2017) record of elevation change from digital elevation models (DEMs) and comparing it to regional climate model output and available records of ice speed. The SWLTS thinned by >12 m on average over the full 32-year period, but the change was highly variable in time and space. Thinning was amplified in the central region of the SWLTS, relative to the north and south. During 1985–2007, the north and south regions demonstrated net thickening while the central region thinned. Regional differences in elevation change are inconsistent with SMB anomalies, indicating that enhanced ice flow in the north and south contributed to thickening during this early time interval. While clear validation in the south is prevented by incomplete velocity data, historical surface speeds in the north were elevated. These findings support the interpretation that changing ice flow can influence ice surface elevation in the slow-moving SWLTS.

Recent hydrological response of glaciers in the Canadian Rockies to changing climate and glacier configuration

Abstract. Mountain snow and ice greatly influence the hydrological cycle of alpine regions by regulating both the quantity of and seasonal variations in water availability downstream. This study considers the combined impacts of climate and glacier changes due to recession on the hydrology and water balance of two high-elevation basins in the Canadian Rockies. A distributed, physically based, uncalibrated glacier hydrology model developed in the Cold Regions Hydrological Modelling platform (CRHM) was used to simulate the glacier mass balance and basin hydrology of the Peyto and Athabasca glacier basins in Alberta, Canada. Bias-corrected reanalysis data were used to drive the model. The model calculates the water balance of glacierized basins, influenced by the surface energy and mass balance, and considers the redistribution of snow by wind and avalanches. It was set up using hydrological response units based on elevation bands, surface slope, and aspect, as well as changing land cover. Aerial photos, satellite images and digital elevation models (DEMs) were assimilated to represent the changing configurations of glacier area and the exposure of ice and firn. Observations of glacier mass balance, snow, and glacier ice surface elevation changes at glacier and alpine tundra meteorological stations and streamflow discharge at the glacier outlets were used to evaluate the model performance. Basin hydrology was simulated over two periods, 1965–1975 and 2008–2018, using the observed glacier configurations for those time periods. Both basins have undergone continuous glacier loss over the last 3 to 5 decades, leading to a 6 %–31 % reduction in glacierized area, a 78 %–109 % increase in ice exposure, and changes to the elevation and slope of the glacier surfaces. Air temperatures are increasing, mainly due to increasing winter maximum and summer minimum daily temperatures. Annual precipitation has increased by less than 11 %, but rainfall ratios have increased by 29 %–44 %. The results show that changes in both climate and glacier configuration have influenced the melt rates and runoff and a shift of peak flows in the Peyto Glacier basin from August to July. Glacier melt contributions increased/decreased from 27 %–61 % to 43 %–59 % of the annual discharges. Recent discharges were 3 %–19 % higher than in the 1960s and 1970s. The results suggest that increased exposure of glacier ice and lower surface elevation due to glacier thinning were less influential than climate warming in increasing streamflow. Streamflow from these glaciers continues to increase.

Monitoring the Hydrological Activities of Antarctic Subglacial Lakes Using CryoSat-2 and ICESat-2 Altimetry Data

Monitoring the hydrological activities of subglacial lakes is critical to understanding the subglacial hydrological system and evaluating the internal mass changes of the Antarctic ice sheet. Drainage or filling events of active lakes lead to elevation changes in the ice surface. These changes can be observed by satellite altimetry, but the monitoring must be conducted continuously since the water movements in active subglacial lakes may occur frequently. We used CryoSat-2 Baseline-D and ICESat-2 data from 2010 to 2020 to obtain the time series of the ice surface elevation changes for 17 active lakes. We also evaluated the uncertainty of the time series derived from the CryoSat-2 data by cross-validation. The mean and RMS of the biases between the CryoSat-2-based and ICESat-2-based time series are generally less than 0.3 m and 1.0 m, respectively. However, the mean and RMS are greater over the lakes with rough ice surfaces, such as Whillans6, KT1, Mac3, and Slessor23. The drainage and filling events continue exhibiting in the extended period and the hydrological activities of SLW, L12, Whillans6, L78, and Mac1 occurred periodically. Furthermore, we inferred the hydrological connections between the lakes combining simulated water pathways.

PRETTY\xa0\u2013 passive GNSS-Reflectometry for CubeSats

The Passive Reflectometry and Dosimetry (PRETTY) is a small satellite mission developed by a consortium consisting of RUAG Space GmbH (RSA), Graz University of Technology (TUG), and Seibersdorf Laboratories (SL). The primary objective of the passive altimetry within the PRETTY mission is to measure the ice and sea ice surface elevation and roughness in the polar regions.

Quaternary ice thinning of David Glacier in the Terra Nova Bay region, Antarctica

Understanding the history of Antarctic glaciation is important for interpreting paleoclimatic changes and estimating the changes in climate, sea level, and ice volume in the future. Ice core studies of the East Antarctic Ice Sheet (EAIS) and marine sediment cores from the entire Ross Sea have employed numerous proxies to reconstruct the glacial history of the Antarctic region. However, the ice and marine core records can be biased because of their specific locations, such as the uppermost accumulation zone or the terminus of the ablation zone, thereby introducing significant uncertainties in ice modeling. In this study, we analyzed 34 new 10 Be and 26 Al samples from four benches that were glaciated in the past by David glacier and incorporate the present ice-free flat surfaces. We suggest that the David glacier experienced monotonic and stepwise vertical lowering along the flanks of Mt. Priestley since the early Pleistocene. The uppermost bedrock benches on Mt. Priestley were exposed at 1.77 ± 0.32 Ma, with no evidence of subsequent overriding by readvancing ice. At Mt. Priestley, the David glacier has been characterized by a cold-based regime since 1.77 Ma, with a denudation rate of only ∼16 cm/Ma, corresponding to the regional transition from warm to cold-based glaciation at 3.5 Ma. Simple exposure ages from two lower benches date to Marine Isotope Stage (MIS) 7 (234.1 ± 13.1 ka; 545 m asl) and MIS 4 (64.8 ± 13.7 ka; 222 m asl), suggesting that, since MIS 8, the overall lowering of glaciers has remained monotonic. The upper bench marks the lower limit of the MIS 8 glacial period and the upper limit of Penultimate Glacial Maximum (MIS 6), while the lower landform defines the upper limit of the last glacial period (MIS 4–2). The magnitude of Quaternary ice thinning at the David glacier was the highest (∼990 m) in the present terminal area (i.e., the most sensitive ablation zone), in contrast to the other outlet glaciers draining into the Terra Nova Bay, which experienced less ice lowering. The combination of the terrestrial (in situ 10 Be and 26 Al) and previous marine (authigenic 10 Be) cosmogenic data used in our study document the history of lowering of the David glacier driven by climatic changes during the Pleistocene. Both deglaciation and glaciation were limited during the mid-Pleistocene transition (MPT) and prior to the mid-Bruhnes event (MBE), due to the prevailing cold and arid climate, whereas deglaciation was dominant during other warm periods. • 10 Be & 26 Al, in combination with previous authigenic 10 Beelucidate the glacial history of David glacier during the entire Quaternary period. • Experienced monotonic terrestrial vertical lowering with complex and multiple cycles of horizontal retreat–advance at its marine margin. • The greatest ice lowering along David glacier during the Quaternary (∼990 m) occurred at its present terminal region. • The ice surface elevation of the East Antarctic Ice Sheet was greater during the Pliocene than the Pleistocene. • Ice thickness was greater in MIS 8 than MIS 6 and the Global LGM.

Ice-sheet expansion from the Ross Sea into McMurdo Sound, Antarctica, during the last two glaciations

An understanding of Antarctic Ice Sheet (AIS) behavior is important for future sea-level predictions. Here, we examine past ice-sheet history in the McMurdo Sound region of the western Ross Sea over the last two glacial-interglacial cycles in order to gain insight into the drivers of ice-sheet change. Surficial mapping, along with radiocarbon dates of lacustrine algae and uranium-thorium disequilibrium dates of lacustrine carbonates from ice-dammed lakes, allow reconstruction of the timing and origin of grounded ice in McMurdo Sound during the last glacial maximum (LGM) and penultimate glaciation. During the LGM, ice-surface elevation profiles and distribution of erratics indicate ice flow into southern McMurdo Sound from the Ross Sea, rather than seaward expansion of local glaciers from the Royal Society Range. The grounded ice in McMurdo Sound flowed westward to block the mouths of valleys in the Royal Society Range and to dam proglacial lakes. In Marshall Valley, maximum ice extent during what is termed the Ross Sea glaciation, was achieved by 18 ka and remained close to this position until after 14 ka. The pattern of surficial deposits suggests that ice during the penultimate glaciation, locally named the Marshall glaciation, was slightly more extensive than that of the LGM but had a similar Ross Sea origin. Maximum ice extent in Marshall Valley occurred at ∼145–150 ka; the grounded ice may have receded from the valley mouth shortly after 138 ka. Both ice expansions occurred broadly during times of low Antarctic air temperatures, which have been linked to insolation minima. However, the lack of widespread surface melting ablation zones on the AIS indicates that the link between ice expansion and orbital forcing is likely to be indirect and possibly driven through the ocean. Closer examination of the precise timing of the glacial maxima in Marshall Valley shows that the Marshall glaciation occurred synchronously with the penultimate global maximum; ice recession took place during Termination II. In contrast, maximum ice extent during the Ross Sea glaciation along the Royal Society Range occurred after the global LGM, during Termination I. Deglaciation was primarily an early Holocene event. We attribute this delayed maximum and deglaciation (relative to global events) to the effect of rising accumulation on ice-sheet mass balance.