Ice Elevation Changes Research Articles

Seawater Intrusion in the Observed Grounding Zone of Petermann Glacier Causes Extensive Retreat

AbstractUnderstanding grounding line dynamics is critical for projecting glacier evolution and sea level rise. Observations from satellite radar interferometry reveal rapid grounding line migration forced by oceanic tides that are several kilometers larger than predicted by hydrostatic equilibrium, indicating the transition from grounded to floating ice is more complex than previously thought. Recent studies suggest seawater intrusion beneath grounded ice may play a role in driving rapid ice loss. Here, we investigate its impact on the evolution of Petermann Glacier, Greenland, using an ice sheet model. We compare model results with observed changes in grounding line position, velocity, and ice elevation between 2010 and 2022. We match the observed retreat, speed up, and thinning using 3‐km‐long seawater intrusion that drive peak ice melt rates of 50 m/yr; but we cannot obtain the same agreement without seawater intrusion. Including seawater intrusion in glacier modeling will increase the sensitivity to ocean warming.

Post-little ice age glacial geomorphology of contrasting topographic settings at Skálafellsjökull, southeast Iceland

ABSTRACT Glacial geomorphological mapping from the southern margin of Skálafellsjökull, southeast Iceland, depicts a topographically diverse mountainside, influencing glacier dynamics, landform formation and glacier retreat since the Little Ice Age maximum in ∼1890. The glacial landforms present are typical of southeast Icelandic temperate glaciers, comprising recessional push moraines, including sawtooth moraines, and associated fluting. Study area A demonstrates an abandoned lobe confined by steep V-shaped topography, displaying moraines and minimal fluting, suggesting low preservation of landforms, and changes in glacier behaviour. At study area B, the sawtooth moraine morphology demonstrates changes in the glacier margin as the ice interacted with a series of topographic benches during active recession. The steep-sided valley at study area C illustrates densely spaced arcuate moraines, reflecting subtle changes in ice elevation. This mapping provides a framework for further investigations into glacier retreat rates and the influence of local topography and climate.

Rapid demise and committed loss of Bowman Glacier, northern Ellesmere Island, Arctic Canada

AbstractUsing historical and recent aerial photography and structure from motion (SfM) multiview stereo (MVS) techniques, we reconstruct the 1959 and 2018 ice surface topography and determine the geodetic mass balance of Bowman Glacier, a small mountain glacier on northern Ellesmere Island. This is combined with optical satellite imagery to reconstruct the evolution in extent of the glacier over six decades, and ground-penetrating radar measurements of ice thickness to estimate the remaining ice volume. Between 1959 and 2020, Bowman Glacier lost 78% of its extent (reducing from 2.75 to 0.61 km2), while average annual area loss rates have nearly tripled in the past two decades. Over the 1959–2018 period, glacier-wide ice-thickness change averaged −22.7 ± 4.7 m, corresponding to a mean specific annual mass balance of −347.0 ± 71.4 mm w.e. a−1. Projecting rates of area and volume change into the future indicates that the glacier will likely entirely disappear between 2030 and 2060. This study demonstrates the potential of SfM-MVS processing to generate elevation products from 1950/60s historical aerial photographs, and to extend observations of ice elevation and glacier volume change for the Canadian Arctic, prior to the satellite record.

Detection of Crustal Uplift Deformation in Response to Glacier Wastage in Southern Patagonia

The Southern Patagonian Icefield (SPI) is the largest continuous ice mass in the Southern Hemisphere outside Antarctica. It has been shrinking since the Little Ice Age (LIA) period, with increasing rates in recent years. An uplift of crustal deformation in response to this deglaciation process has been expected. The goal of this investigation is to analyze the crustal deformation caused by ice retreat using time-series data from continuous GPS stations (2015–2020) in the northern area of the SPI. For this purpose, we installed two continuous GPS stations on rocky nunataks of the SPI (the GRCS near Greve glacier and the GBCS close by Cerro Gorra Blanca). In addition, ice elevation changes (2000–2019) were analyzed by the co-registration of the SRTM digital elevation model and ICESat elevation data points. The results of the vertical components are positive (36.55 ± 2.58 mm a−1), with a maximum at GBCS, indicating the highest rate of crustal uplift ever continuously recorded in Patagonia; in addition, the mean horizontal velocities reached 11.7 mm a−1 with an azimuth of 43°. The negative ice elevation changes detected in the region have also accelerated in the recent two decades, with a median Δh (elevation change) of −3.36 ± 0.01 m a−1 in the ablation zone. The seasonality of the GPS signals was contrasted with the water levels of the main Patagonian lakes around the SPI, detecting a complex interplay between them. Hence, the study sheds light on the knowledge of the crustal uplift as evidence of the wastage experienced by the SPI glaciers.

Ice Elevation Change Based on GNSS Measurements along the Korth-Traverse in Southern Greenland

In 1912, a Swiss expedition led by meteorologist Alfred de Quervain crossed the Greenland ice sheet on a route from Disko Bay to Tasiilaq. Based on that, in 2002, a series of geodetic expeditions carried out by W. Korth and later by T. Hitziger began along the same traverse as in 1912, with the last measurements taken in May 2021. The statically collected GPS/GNSS data provide very accurate elevation changes at 36 points along the almost 700 km long crossing over a period of 19 years. According to this, there is a maximum increase of 2.1 m in the central area and a decrease of up to 38.7 m towards the coasts (influence Ilulissat Isbræ). By using kinematic GNSS measurements, there is a very dense profile with a spacing of a few meters. The comparison of those measurements is performed using crossing points or minimum distances and gives equivalent results for both methods. It is shown that local ice topography is preserved, and thus gaps in data sets can be caught. Areas of accumulation and ablation on the ice sheet can be identified, showing the widespread influence of outlet glaciers up to 200 km. The data can be used for direct verification of altimetry data, such as IceSat. Both IceSat elevations and their changes can be compared.

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 changes in the McMurdo Ice Shelf Transition Zone and Hut Point Peninsula, West Antarctica

McMurdo Ice Shelf (MIS) and Ross Island act as both lateral margins and pinning points for the Ross Ice Shelf, as well as logistic hubs for the United States and New Zealand Antarctic Programs. Recent thinning and retreat of MIS has motivated an evaluation of the future stability of this critical region and potential adaptation. If MIS were to collapse, it could initiate sizeable teleconnections across West Antarctica, as well as stifle Antarctic research that relies on MIS/Ross Island-based logistics. One logistics component already experiencing change is the Transition Zone (TZ) road, which connects MIS with the research stations on Hut Point Peninsula (HPP), Ross Island. Here we assess the vulnerability of the TZ road and a proposed rerouting location (TZ Hillside), by using a combination of in situ, remote sensing, and numerical modeling approaches. We evaluate overall elevation change by differencing late summer, WorldView-derived, DEMs (2011–2015), and evaluate the practicality of both a stable ice surface and bare earth road using a numerical ice volume flux model, informed by GPR and GPS measurements of ice thickness and velocity (2015–2016). Results show overall ice elevation change of up to 3 m, and 204 ± 24 m3 a−1 of ice flux across the TZ Hillside. Combined, these results indicate that the TZ Hillside is unsuitable for either type of road, and point to MIS terminus retreat and TZ Hillside melt as a continuing risk for TZ stability, which could prelude broader changes.

Satellite and ground based estimates for ice surface velocities in the part of central Dronning Maud Land, East Antarctica: Implications for ice flux calculations

The ice sheet and glaciers of Antarctica and Greenland represent the largest sources of freshwater on planet Earth. The understanding and quantification of their dynamic properties such as albedo, precipitation, ice mass movement, and ice elevation changes are critical for the improved climate and mass balance models. The present study utilizes space-borne optical and synthetic aperture radar (SAR) imagery to measure the ice surface velocity at high spatial resolution for a part of the central Dronning Maud Land (cDML), East Antarctica. The datasets from Landsat-8 and Sentinel-1 SAR satellite are used for ice stream velocity estimation using feature-offset tracking and differential interferometric SAR (DInSAR) methods. The derived velocity products are validated with ground based stakes network at annual time scale. The fundamental ice flow laws are used to estimate the ice outflux or discharge for selected ice stream drainage basins of cDML at fluxgate locations. The ice stream basin has been delineated using combination of elevation, slope and continental scale velocity maps. The ice influx for study area is estimated using ECMWF fifth generation reanalysis (ERA5) and Regional Atmospheric Climate Model (RACMO) v2.3 model outputs. The estimated influx and outflux are in the ranges of 0.18–4.167 Gt/y and 0.201 to 1.278 Gt/y respectively, indicating net positive mass balance for the selected area.

Increased Ice Thinning over Svalbard Measured by ICESat/ICESat-2 Laser Altimetry

A decade-long pronounced increase in temperatures in the Arctic, especially in the Barents Sea region, resulted in a global warming hotspot over Svalbard. Associated changes in the cryosphere are the consequence and lead to a demand for monitoring of the glacier changes. This study uses spaceborne laser altimetry data from the ICESat and ICESat-2 missions to obtain ice elevation and mass change rates between 2003–2008 and 2019. Elevation changes are derived at orbit crossover locations throughout the study area, and regional volume and mass changes are estimated using a hypsometric approach. A Svalbard-wide annual elevation change rate of −0.30 ± 0.15 m yr−1 was found, which corresponds to a mass loss rate of −12.40 ± 4.28 Gt yr−1. Compared to the ICESat period (2003–2009), thinning has increased over most regions, including the highest negative rates along the west coast and areas bordering the Barents Sea. The overall negative regime is expected to be linked to Arctic warming in the last decades and associated changes in glacier climatic mass balance. Further, observed increased thinning rates and pronounced changes at the eastern side of Svalbard since the ICESat period are found to correlate with atmospheric and oceanic warming in the respective regions.

A near 90-year record of the evolution of El Morado Glacier and its proglacial lake, Central Chilean Andes

AbstractUsing an ensemble of close- and long-range remote sensing, lake bathymetry and regional meteorological data, we present a detailed assessment of the geometric changes of El Morado Glacier in the Central Andes of Chile and its adjacent proglacial lake between 1932 and 2019. Overall, the results revealed a period of marked glacier down wasting, with a mean geodetic glacier mass balance of −0.39 ± 0.15 m w.e.a−1observed for the entire glacier between 1955 and 2015 with an area loss of 40% between 1955 and 2019. We estimate an ice elevation change of −1.00 ± 0.17 m a−1for the glacier tongue between 1932 and 2019. The increase in the ice thinning rates and area loss during the last decade is coincident with the severe drought in this region (2010–present), which our minimal surface mass-balance model is able to reproduce. As a result of the glacier changes observed, the proglacial lake increased in area substantially between 1955 and 2019, with bathymetry data suggesting a water volume of 3.6 million m3in 2017. This study highlights the need for further monitoring of glacierised areas in the Central Andes. Such efforts would facilitate a better understanding of the downstream impacts of glacier downwasting.

ICESat-2星载光子计数激光雷达数据处理与应用研究进展

The ICESat-2/ATLAS employed a micro-pulse multi-beam photon-counting laser altimeter system for the first time. ICESat-2/ALTAS was significantly different from ICESat-1/GLAS in terms of detection mechanism, data processing methods and the breadth and depth of data applications. Firstly, the key configuration parameters, data and products of ICESat-2/ATLAS were introduced. Secondly, the key technologies of noise removal and classification algorithms for ICESat-2/ATLAS data as well as the adaptability and problems of each algorithm were elaborated. Thirdly, the applications of ICESat-2/ATLAS data in ice sheet and sea ice elevation measurement and change monitoring, ground elevation and forest height extraction, forest biomass estimation, lake level and storage change monitoring were summarized. Finally, the development trend and prospect of photon-counting LiDAR data processing methods and applications were expected.

Monitoring changes of the Antarctic Ice sheet by GRACE, ICESat and GNSS

Abstract In this study, we estimate the ice mass changes, the ice elevation changes and the vertical displacements in Antarctica based on analysis of multi-geodetic datasets that involve the satellite gravimetry (GRACE), the satellite altimetry (ICESat) and the global navigation satellite systems (GNSS). According to our estimates, the total mass change of the Antarctic ice sheet from GRACE data is −162.91 Gt/yr over the investigated period between April 2002 and June 2017. This value was obtained after applying the GIA correction of −98.12 Gt/yr derived from the ICE-5G model of the glacial iso-static adjustment. A more detailed analysis of mass balance changes for three individual drainage regions in Antarctica reveal that the mass loss of the West Antarctic ice sheet was at a rate of −143.11 Gt/yr. The mass loss of the Antarctic Peninsula ice sheet was at a rate of −24.31 Gt/yr. The mass of the East Antarctic ice sheet increased at a rate of 5.29 Gt/yr during the investigated period. When integrated over the entire Antarctic ice sheet, average rates of ice elevation changes over the period from March 2003 to October 2009 derived from ICESat data represent the loss of total ice volume of −155.6 km3.The most prominent features in ice volume changes in Antarctica are characterized by a strong dynamic thinning and ice mass loss in the Amundsen Sea Embayment that is part of the West Antarctic ice sheet. In contrast, coastal regions between Dronning Maud Land and Enderby Land exhibit a minor ice increase, while a minor ice mass loss is observed in Wilkes Land. The vertical load displacement rates estimated from GRACE and GPS data relatively closely agree with the GIA model derived based on the ice-load history and the viscosity profile. For most sites, the GRACE signal appears to be in phase and has the same amplitude as that obtained from the GPS vertical motions while other sites exhibit some substantial differences possibly attributed to thermo-elastic deformations associated with surface temperature.

Net retreat of Antarctic glacier grounding lines

Grounding lines are a key indicator of ice-sheet instability, because changes in their position reflect imbalance with the surrounding ocean and affect the flow of inland ice. Although the grounding lines of several Antarctic glaciers have retreated rapidly due to ocean-driven melting, records are too scarce to assess the scale of the imbalance. Here, we combine satellite altimeter observations of ice-elevation change and measurements of ice geometry to track grounding-line movement around the entire continent, tripling the coverage of previous surveys. Between 2010 and 2016, 22%, 3% and 10% of surveyed grounding lines in West Antarctica, East Antarctica and at the Antarctic Peninsula retreated at rates faster than 25 m yr−1 (the typical pace since the Last Glacial Maximum) and the continent has lost 1,463 km2 ± 791 km2 of grounded-ice area. Although by far the fastest rates of retreat occurred in the Amundsen Sea sector, we show that the Pine Island Glacier grounding line has stabilized, probably as a consequence of abated ocean forcing. On average, Antarctica’s fast-flowing ice streams retreat by 110 metres per metre of ice thinning. Grounding lines in parts of West Antarctica, East Antarctica and the Antarctic Peninsula retreated faster than typical post-glacial pace, according to satellite observations and ice geometry measurements.

CryoSat-2 swath interferometric altimetry for mapping ice elevation and elevation change

For more than 25 years, satellite radar altimetry has provided continuous information on the state of the cryosphere and on its contribution to global sea-level rise. The technique typically delivers maps of ice-sheet elevation and elevation change with 3–10 km spatial resolution and seasonal to monthly temporal resolution. Here we show how the interferometric mode of CryoSat-2 can be used to map broad (5 km-wide) swaths of surface elevation with fine (500 m) spatial resolution from each satellite pass, providing a step-change in the capability of satellite altimetry for glaciology. These swaths of elevation data contain up to two orders of magnitude more surface elevation measurements than standard altimeter products, which provide single elevation measurements based on the range to the Point-Of-Closest-Approach (POCA) in the vicinity of the sub-satellite ground track. The swath elevations allow a more dense, statistically robust time series of elevation change to be formed with temporal resolution of a factor 5 higher than for POCA. The mean differences between airborne altimeter and CryoSat-2 derived ice sheet elevations and elevation rates range from −0.93 ± 1.17 m and 0.29 ± 1.25 m a−1, respectively, at the POCA, to −1.50 ± 1.73 m and 0.04 ± 1.04 m a−1, respectively, across the entire swath. We demonstrate the potential of these data by creating and evaluating elevation models of: (i) the Austfonna Ice Cap (Svalbard), (ii) western Greenland, and (iii) Law Dome (East Antarctica); and maps of ice elevation change of: (iv) the Amundsen Sea sector (West Antarctica), (v) Icelandic ice caps, and (vi) above an active subglacial lake system at Thwaites Glacier (Antarctica), each at 500 m spatial posting – around 10 times finer than possible using traditional approaches based on standard altimetry products.

Reply to “Basal buoyancy and fast-moving glaciers: in defense of analytic force balance” by C. J. van der Veen (2016)

Abstract. Two approaches to ice-sheet modeling are available. Analytical modeling is the traditional approach (Van der Veen, 2016). It solves the force (momentum), mass, and energy balances to obtain three-dimensional solutions over time, beginning with the Navier–Stokes equations for the force balance. Geometrical modeling employs simple geometry to solve the force and mass balance in one dimension along ice flow (Hughes, 2012a). It is useful primarily to provide the first-order physical basis of ice-sheet modeling for students with little background in mathematics. The geometric approach uses changes in ice-bed coupling along flow to calculate changes in ice elevation and thickness, using a floating fraction ϕ along a flow line or flow band, where ϕ = 0 for sheet flow, 0 < ϕ < 1 for stream flow, and ϕ = 1 for shelf flow. An attempt is made to reconcile the two approaches.

Recent Deceleration of the Ice Elevation Change of Ecology Glacier (King George Island, Antarctica)

Glacier change studies in the Antarctic Peninsula region, despite their importance for global sea level rise, are commonly restricted to the investigation of frontal position changes. Here we present a long term (37 years; 1979–2016) study of ice elevation changes of the Ecology Glacier, King George Island ( 62 ∘ 11 ′ S, 58 ∘ 29 ′ W). The glacier covers an area of 5.21 km 2 and is located close to the H. Arctowski Polish Antarctic Station, and therefore has been an object of various multidisciplinary studies with subject ranging from glaciology, meteorology to glacial microbiology. Hence, it is of great interest to assess its current state and put it in a broader context of recent glacial change. In order to achieve that goal, we conducted an analysis of archival cartographic material and combined it with field measurements of proglacial lagoon hydrography and state-of-art geodetic surveying of the glacier surface with terrestrial laser scanning and satellite imagery. Overall mass loss was largest in the beginning of 2000s, and the rate of elevation change substantially decreased between 2012–2016, with little ice front retreat and no significant surface lowering. Ice elevation change rate for the common ablation area over all analyzed periods (1979–2001–2012–2016) has decreased from −1.7 ± 0.4 m/year in 1979–2001 and −1.5 ± 0.5 m/year in 2001–2012 to −0.5 ± 0.6 m/year in 2012–2016. This reduction of ice mass loss is likely related to decreasing summer temperatures in this region of the Antarctic Peninsula.

An ice sheet model of reduced complexity for paleoclimate studies

Abstract. IceBern2D is a vertically integrated ice sheet model to investigate the ice distribution on long timescales under different climatic conditions. It is forced by simulated fields of surface temperature and precipitation of the Last Glacial Maximum and present-day climate from a comprehensive climate model. This constant forcing is adjusted to changes in ice elevation. Due to its reduced complexity and computational efficiency, the model is well suited for extensive sensitivity studies and ensemble simulations on extensive temporal and spatial scales. It shows good quantitative agreement with standardized benchmarks on an artificial domain (EISMINT). Present-day and Last Glacial Maximum ice distributions in the Northern Hemisphere are also simulated with good agreement. Glacial ice volume in Eurasia is underestimated due to the lack of ice shelves in our model. The efficiency of the model is utilized by running an ensemble of 400 simulations with perturbed model parameters and two different estimates of the climate at the Last Glacial Maximum. The sensitivity to the imposed climate boundary conditions and the positive degree-day factor β, i.e., the surface mass balance, outweighs the influence of parameters that disturb the flow of ice. This justifies the use of simplified dynamics as a means to achieve computational efficiency for simulations that cover several glacial cycles. Hysteresis simulations over 5 million years illustrate the stability of the simulated ice sheets to variations in surface air temperature.

Detailed ice loss pattern in the northern Antarctic Peninsula: widespread decline driven by ice front retreats

Abstract. The northern Antarctic Peninsula (nAP, < 66° S) is one of the most rapidly changing glaciated regions on earth, yet the spatial patterns of its ice mass loss at the glacier basin scale have to date been poorly documented. We use satellite laser altimetry and satellite stereo-image topography spanning 2001–2010, but primarily 2003–2008, to map ice elevation change and infer mass changes for 33 glacier basins covering the mainland and most large islands in the nAP. Rates of ice volume and ice mass change are 27.7± 8.6 km3 a−1 and 24.9± 7.8 Gt a−1, equal to −0.73 m a−1 w.e. for the study area. Mass loss is the highest for eastern glaciers affected by major ice shelf collapses in 1995 and 2002, where twelve glaciers account for 60% of the total imbalance. However, losses at smaller rates occur throughout the nAP, at both high and low elevation, despite increased snow accumulation along the western coast and ridge crest. We interpret the widespread mass loss to be driven by decades of ice front retreats on both sides of the nAP, and extended throughout the ice sheet due to the propagation of kinematic waves triggered at the fronts into the interior.

Rapid bedrock uplift in the Antarctic Peninsula explained by viscoelastic response to recent ice unloading

Since 1995 several ice shelves in the Northern Antarctic Peninsula have collapsed and triggered ice-mass unloading, invoking a solid Earth response that has been recorded at continuous GPS (cGPS) stations. A previous attempt to model the observation of rapid uplift following the 2002 breakup of Larsen B Ice Shelf was limited by incomplete knowledge of the pattern of ice unloading and possibly the assumption of an elastic-only mechanism. We make use of a new high resolution dataset of ice elevation change that captures ice-mass loss north of 66°S to first show that non-linear uplift of the Palmer cGPS station since 2002 cannot be explained by elastic deformation alone. We apply a viscoelastic model with linear Maxwell rheology to predict uplift since 1995 and test the fit to the Palmer cGPS time series, finding a well constrained upper mantle viscosity but less sensitivity to lithospheric thickness. We further constrain the best fitting Earth model by including six cGPS stations deployed after 2009 (the LARISSA network), with vertical velocities in the range 1.7 to 14.9 mm/yr. This results in a best fitting Earth model with lithospheric thickness of 100–140 km and upper mantle viscosity of 6×1017–2×1018 Pas – much lower than previously suggested for this region. Combining the LARISSA time series with the Palmer cGPS time series offers a rare opportunity to study the time-evolution of the low-viscosity solid Earth response to a well-captured ice unloading event.

航空機高度計によるアラスカ/ユーコンにおける氷河表面標高変化の測定 (2009-2012)

航空機高度計によるアラスカ/ユーコンにおける氷河表面標高変化の測定 (2009-2012)