Glacier Volume Change Research Articles

Assessment of methods for reconstructing Little Ice Age glacier surfaces on the examples of Novaya Zemlya and the Swiss Alps

Climate change since the end of the Little Ice Age (LIA) has driven observed glacier volume loss, significantly contributing to the rise of global sea level. The related calculation of volume change requires knowledge of glacier surfaces from at least two points in time, usually represented by two digital elevation models (DEMs). These are typically derived from photogrammetric techniques using stereo images, but such images do not go back to the LIA. Accordingly, several techniques have been developed to reconstruct LIA glacier surfaces from historic outlines and modern DEMs. Here we first evaluate various surface interpolation methods by replicating modern glacier surfaces from outline elevation points and analyse elevation differences and uncertainties. Secondly, we investigate different GIS-based methods for LIA surface reconstruction including a new method that is based on up-scaling of recent glacier-specific elevation change data and works also for ice caps without lateral moraines. The methods were tested on 90 glaciers (covering 643 km2) in southern Novaya Zemlya and 266 glaciers (524 km2) in the Bernese Alps of Switzerland. As in previous studies, we also found that the Natural Neighbor and Topo to Raster interpolation methods in ESRI's ArcGIS performed best for glacier surface reconstruction and that all methods are challenged by replicating the variable surface curvature of glaciers. The new reconstruction method shows the smallest mean difference to the reference dataset (RMSE of 26.7 vs. 39.8 m). The often neglected small surface lowering in the accumulation area can increase the derived glacier volume changes by 30–50 % and should thus be considered whenever possible. Applying the up-scaling method to both test regions revealed that elevation change rates over the last two decades (−0.8 m a−1 for Novaya Zemlya and -1.14 m a−1 for the Bernese Alps) were much higher than from 1850 until today (−0.13 m a−1 and -0.22 m a−1, respectively).

Rapid mass losses of Urumqi River Basin glaciers, eastern Tianshan Mountains revealed from multi-temporal DEMs, 1964–2021

ABSTRACT Glacier volume changes have remarkable impacts on regional water resources; however, their estimation uncertainties remain large. In this study, the surface elevations of seven glaciers at the headwaters of the Urumqi River exhibited obvious thinning at the rate of 0.32 m a−1 from 1964 to 2021. Correlations between the ice volume change (dV) derived from surface elevation conversion and the area change (dS) were established for valley, cirque, and hanging glaciers. The ice volumes of valley, cirque, and hanging glaciers in the Urumqi River Basin decreased by 0.99 km3, 0.18 km3, and 0.59 km3, respectively, and the total ice volume decreased by 1.76 km3 during 1964–2021. Ice volume losses are accelerating. The correlation between dV and dS for valley glaciers produced a large estimation difference for the ice volume change of Urumqi Glacier No. 1 (UG1) compared with the glaciological mass balance (15%), while the estimation differences of ice volume changes for the three types of glaciers in the Urumqi River Basin were within acceptable limits (10%) based on the geodetic results. This new method has vast potential in estimating ice volume changes at the basin or regional scale, or even across western China.

Late Holocene glacier and climate fluctuations in the Mackenzie and Selwyn mountain ranges, northwestern Canada

Abstract. Over the last century, northwestern Canada experienced some of the highest rates of tropospheric warming globally, which caused glaciers in the region to rapidly retreat. Our study seeks to extend the record of glacier fluctuations and assess climate drivers prior to the instrumental record in the Mackenzie and Selwyn mountains of northwestern Canada. We collected 27 10Be surface exposure ages across nine cirque and valley glacier moraines to constrain the timing of their emplacement. Cirque and valley glaciers in this region reached their greatest Holocene extents in the latter half of the Little Ice Age (1600–1850 CE). Four erratic boulders, 10–250 m distal from late Holocene moraines, yielded 10Be exposure ages of 10.9–11.6 ka, demonstrating that by ca. 11 ka, alpine glaciers were no more extensive than during the last several hundred years. Estimated temperature change obtained through reconstruction of equilibrium line altitudes shows that since ca. 1850 CE, mean annual temperatures have risen 0.2–2.3 ∘C. We use our glacier chronology and the Open Global Glacier Model (OGGM) to estimate that from 1000 CE, glaciers in this region reached a maximum total volume of 34–38 km3 between 1765 and 1855 CE and had lost nearly half their ice volume by 2019 CE. OGGM was unable to produce modeled glacier lengths that match the timing or magnitude of the maximum glacier extent indicated by the 10Be chronology. However, when applied to the entire Mackenzie and Selwyn mountain region, past millennium OGGM simulations using the Max Planck Institute Earth System Model (MPI-ESM) and the Community Climate System Model 4 (CCSM4) yield late Holocene glacier volume change temporally consistent with our moraine and remote sensing record, while the Meteorological Research Institute Earth System Model 2 (MRI-ESM2) and the Model for Interdisciplinary Research on Climate (MIROC) fail to produce modeled glacier change consistent with our glacier chronology. Finally, OGGM forced by future climate projections under varying greenhouse gas emission scenarios predicts 85 % to over 97 % glacier volume loss by the end of the 21st century. The loss of glaciers from this region will have profound impacts on local ecosystems and communities that rely on meltwater from glacierized catchments.

Climate Variability and Glacier Evolution at Selected Sites Across the World: Past Trends and Future Projections

AbstractThe availability of freshwater from glaciers and snowmelt is of vital importance for people and ecosystems in the context of global climate change. Here, we focus on 25 glaciers located in different climates and latitudes and investigate their recent (1958–2020) and future projected trends (2020–2050 and 2070–2100) in monthly precipitation (Pr), maximum and minimum temperatures, ice mass loss, and their relationships with cloud properties. The study sites are located in Temperate Europe (France), the Inner (Ecuador, Venezuela, and Colombia) and Outer Tropics (Bolivia and Peru), Central America (Mexico), tropical Southeast Asia (Indonesia), Equatorial Africa (Uganda), and the Southern dry and Patagonian Andes (Chile and Argentina). The climate analyses are based on TerraClimate data (Monthly Climate and Climatic Water Balance for Global Terrestrial Surfaces) and 28 CORDEX (Coordinated Regional Climate Downscaling Experiment) climate simulations. Our findings reveal that, extrapolating current glacier volume change trends, almost half of the studied glaciers are likely to vanish (95%–100% volume loss) by 2050, with widespread warming and drying trends since 1958. A shift toward wetter conditions at Pico Humboldt (Venezuela) and Martial Este (Argentina) identifiable in the CORDEX simulation will very likely not have a limiting impact on glacier mass loss owing to increasing temperatures, which will raise the elevation of the rain/snow limit. Our results provide useful new information to better understand glacier‐climate relationships and future scenarios dominated by ice mass loss trends across the globe. These findings suggest serious consequences for future water availability, which exacerbate the vulnerability of local populations and ecosystems.

Applying UAV-Based Remote Sensing Observation Products in High Arctic Catchments in SW Spitsbergen

In the age of remote sensing, particularly with new generation Uncrewed Aerial Vehicles (UAVs), there is a broad spectrum of applications, especially in remote and rapidly changing areas such as the Arctic. Due to challenging conditions in this region, there is a scarcity of detailed spatial studies with data that may be used to estimate changes in glacier volume and geomorphological changes caused by permafrost freeze–thaw cycles. Drone-based Digital Elevation Models (DEM) offer a finer spatial resolution with higher accuracy than airborne and satellite-based products that can be used for acquiring, interpreting, and precisely representing spatial data in broad studies. In this study, we evaluate a UAV-based DEM of two High Arctic catchments, Fuglebekken and Ariebekken, located on Spitsbergen Island. The surveys were carried out in July 2022 using a DJI Matrice 300 RTK drone equipped with a photogrammetric Zenmuse P1 camera. A total of 371 images were taken, covering an area of 7.81 km2. The DEM was created by the Structure-from-Motion technique and achieved a centimetre-level accuracy by overlapping very high-resolution images. The final resolution of the DEM was found to be 0.06 m in Fuglebekken and 0.07 m in Ariebekken, with a horizontal and vertical RMSE of 0.09 m and 0.20 m, respectively. The DJI Matrice 300 RTK drone-based DEM is compared and correlated with the aerial mission of the Svalbard Integrated Arctic Earth Observing System (SIOS) conducted in July 2020 and the satellite-based ArcticDEM acquired in July 2018. This allowed the detection of elevation changes and identification of landscape evolution, such as moraine breaches and coastal erosion. We also highlight the usage of DEM in providing detailed morphometric characteristics and hydrological parameters, such as the delineation of catchments and stream channels. The final products are available at the IG PAS Data Portal.

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.

Monitoring glacier characteristics and their mass balance using a multidimensional approach over the glaciers of the Chandra basin, western Himalaya

ABSTRACT Glacier change monitoring provides a significant understanding of glacier health and its influence on freshwater availability. Therefore, we have utilized remotely sensed data to estimate the glacier characteristics, mass balance, and volume change during 2013–2019 over the Sutri Dhaka and Batal glaciers, Chandra basin, western Himalayas. The glacier feature classification and shifting in isoline were mapped on a spatio-temporal scale. Furthermore, the surface velocity was measured using Coregistration of Optically Sensed Images and Correlation (COSI-Corr) and then used for estimating the ice thickness. The estimated velocity was compared with the stake observations and National Aeronautics and Space Administration (NASA)- Inter-Mission Time Series of Land Ice Velocity and Elevation (ITS_LIVE) data. Further, the glaciers’ geodetic mass balance and volume change were calculated, indicating an average value of −0.03 ± 0.5, −0.34 ± 0.40 and −0.13 ± 0.14 m water equivalent (w.e.) year−1 during 2000–2009, 2009–2020, and 2000–2020, respectively. This multidimensional analysis of the glaciers is critical, as the meltwater from both glaciers contributes to the Chandra basin and also influences the local ecosystem.

Uncertainty Analysis of Digital Elevation Models by Spatial Inference From Stable Terrain

The monitoring of Earth’s and planetary surface elevations at larger and finer scales is rapidly progressing through the increasing availability and resolution of digital elevation models (DEMs). Surface elevation observations are being used across an expanding range of fields to study topographical attributes and their changes over time, notably in glaciology, hydrology, volcanology, seismology, forestry and geomorphology. However, DEMs frequently contain large-scale instrument noise and varying vertical precision that lead to complex patterns of errors. Here, we present a validated statistical workflow to estimate, model, and propagate uncertainties in DEMs. We review the state-of-the-art of DEM accuracy and precision analyses, and define a conceptual framework to consistently address those. We show how to characterize DEM precision by quantifying the heteroscedasticity of elevation measurements, i.e. varying vertical precision with terrain- or sensor-dependent variables, and the spatial correlation of errors that can occur across multiple spatial scales. With the increasing availability of high-precision observations, our workflow based on independent elevation data acquired on stable terrain can be applied almost anywhere on Earth. We illustrate how to propagate uncertainties for both pixel-scale and spatial elevation derivatives, using terrain slope and glacier volume changes as examples. We find that uncertainties in DEMs are largely underestimated in the literature, and advocate that new metrics of DEM precision are essential to ensure the reliability of future Earth and planetary surface elevation assessments.

Identificación de clústeres glaciares a lo largo de los Andes chilenos usando variables topoclimáticas

Utilizando variables topográficas y climáticas, presentamos clústeres glaciares en los Andes chilenos (17.6-55.4°S), donde se ejecutó el algoritmo de aprendizaje automático no supervisado Partitioning Around Medoids (PAM). Los resultados clasificaron 23,974 glaciares dentro de trece clústeres, que muestran condiciones específicas en términos de cantidades anuales y mensuales de precipitación, temperatura y radiación solar. En los Andes secos, los valores medios anuales de cinco clústeres glaciares (C1-C5) muestran una diferencia de precipitación y temperatura de hasta 400 mm (29 y 33°S) y 8°C (33°S), con una resta de elevación promedio de 1800 m entre glaciares clústeres C1 y C5 (18 a 34°S). Mientras que en los Andes húmedos las mayores diferencias se observaron en la latitud del Campo de Hielo Patagónico Sur (50°S), donde los valores medios anuales de precipitación y temperatura muestran una precipitación marítima por encima de 3700 mm/año (C12), donde el aire húmedo occidental juega un papel importante, y por debajo de 1000 mm/año al este del Campo de Hielo Patagónico Sur (C10), con diferencias de temperatura cercanas a 4°C y una resta de elevación promedio de 500 m. Esta clasificación confirma que los glaciares chilenos no pueden agruparse solo latitudinalmente, contribuyendo a una mejor comprensión de los cambios recientes en el volumen de los glaciares a escala regional.

Novel Techniques for Void Filling in Glacier Elevation Change Data Sets

The increasing availability of digital elevation models (DEMs) facilitates the monitoring of glacier mass balances on local and regional scales. Geodetic glacier mass balances are obtained by differentiating DEMs. However, these computations are usually affected by voids in the derived elevation change data sets. Different approaches, using spatial statistics or interpolation techniques, were developed to account for these voids in glacier mass balance estimations. In this study, we apply novel void filling techniques, which are typically used for the reconstruction and retouche of images and photos, for the first time on elevation change maps. We selected 6210 km2 of glacier area in southeast Alaska, USA, covered by two void-free DEMs as the study site to test different inpainting methods. Different artificially voided setups were generated using manually defined voids and a correlation mask based on stereoscopic processing of Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) acquisition. Three “novel” (Telea, Navier–Stokes and shearlet) as well as three “classical” (bilinear interpolation, local and global hypsometric methods) void filling approaches for glacier elevation data sets were implemented and evaluated. The hypsometric approaches showed, in general, the worst performance, leading to high average and local offsets. Telea and Navier–Stokes void filling showed an overall stable and reasonable quality. The best results are obtained for shearlet and bilinear void filling, if certain criteria are met. Considering also computational costs and feasibility, we recommend using the bilinear void filling method in glacier volume change analyses. Moreover, we propose and validate a formula to estimate the uncertainties caused by void filling in glacier volume change computations. The formula is transferable to other study sites, where no ground truth data on the void areas exist, and leads to higher accuracy of the error estimates on void-filled areas. In the spirit of reproducible research, we publish a software repository with the implementation of the novel void filling algorithms and the code reproducing the statistical analysis of the data, along with the data sets themselves.

Future projection of cryospheric and hydrologic regimes in Koshi River basin, Central Himalaya, using coupled glacier dynamics and glacio-hydrological models

AbstractClimate-induced cryospheric changes can have a significant impact on the downstream water availability. In this study, the Open Global Glacier Model (OGGM) and the Glacio-hydrological Degree-day Model (GDM) are integrated to project the response of cryospheric and hydrological systems to climate change until 2100. The study area comprises six sub-basins of glacierized Koshi River basin covering Nepalese and Chinese territories. The output from OGGM is provided as input to GDM along with the spatial and hydro-meteorological data. The average glacier area change in all the sub-basins from 2021 to 2100 is estimated as 65 and 85% decrease and the average glacier volume change is estimated as 76 and 86% decrease for RCP 4.5 and 8.5 scenarios, respectively. The future simulated discharge shows an increasing trend in pre-monsoon and monsoon seasons and a decreasing trend in post-monsoon and winter seasons after 2060 in all the sub-basins, which can lead to wetter wet seasons and drier dry seasons in the far future. A shift in peak flow is observed from August to July in most of the sub-basins. The coupled modelling technique used in this study can largely improve our understanding of glacio-hydrological dynamics in the Himalayan region.

Volume-area scaling for debris-covered glaciers

AbstractA volume-area scaling relation is commonly used to estimate glacier volume or its future changes on a global scale. The presence of an insulating supraglacial debris cover alters the mass-balance profile of a glacier, potentially modifying the scaling relation. Here, the nature of scaling relations for extensively debris-covered glaciers is investigated. Theoretical arguments suggest that the volume-area scaling exponent for these glaciers is ~7% smaller than that for clean glaciers. This is consistent with the results from flowline-model simulations of idealised glaciers, and the available data from the Himalaya. The best-fit scale factor for debris-covered Himalayan glaciers is ~60% larger compared to that for the clean ones, implying a significantly larger stored ice volume in a debris-covered glacier compared to a clean one having the same area. These results may help improve scaling-based estimates of glacier volume and future glacier changes in regions where debris-covered glaciers are abundant.

Deformation Time Series and Driving-Force Analysis of Glaciers in the Eastern Tienshan Mountains Using the SBAS InSAR Method.

Glacier melting is one of the important causes of glacier morphology change and can provide basic parameters for calculating glacier volume change and glacier mass balance, which, in turn, is important for evaluating water resources. However, it is difficult to obtain large-scale time series of glacier changes due to the cloudy and foggy conditions which are typical of mountain areas. Gravity-measuring satellites and laser altimetry satellites can monitor changes in glacier volume over a wide area, while synthetic-aperture radar satellites can monitoring glacier morphology with a high spatial and temporal resolution. In this article, an interferometric method using a short temporal baseline and a short spatial baseline, called the Small Baseline Subsets (SBAS) Interferometric Synthetic-Aperture Radar (InSAR) method, was used to study the average rate of glacier deformation on Karlik Mountain, in the Eastern Tienshan Mountains, China, by using 19 Sentinel-1A images from November 2017 to December 2018. Thus, a time series analysis of glacier deformation was conducted. It was found that the average glacier deformation in the study region was −11.77 ± 9.73 mm/year, with the observation sites generally moving away from the satellite along the Line of Sight (LOS). Taking the ridge line as the dividing line, it was found that the melting rate of southern slopes was higher than that of northern slopes. According to the perpendicular of the mountain direction, the mountain was divided into an area in the northwest with large glaciers (Area I) and an area in the southeast with small glaciers (Area II). It was found that the melting rate in the southeast area was larger than that in the northwest area. Additionally, through the analysis of temperature and precipitation data, it was found that precipitation played a leading role in glacier deformation in the study region. Through the statistical analysis of the deformation, it was concluded that the absolute value of deformation is large at elevations below 4200 m while the absolute value of the deformation is very small at elevations above 4500 m; the direction of deformation is always away from the satellite along the LOS and the absolute value of glacier deformation decreases with increasing elevation.

Change in the Extent of Glaciers and Glacier Runoff in the Chinese Sector of the Ile River Basin between 1962 and 2012

Change in glacier area in the Kuksu and Kunes river basins, which are tributaries to the internationally important Ile River, were assessed at two different time steps between 1962/63, 1990/93, and 2010/12. Overall, glaciers lost 191.3 ± 16.8 km2 or 36.9 ± 6.5% of the initial area. Glacier wastage intensified in the latter period: While in 1962/63–1990/93 glaciers were losing 0.5% a−1, in 1990/93–2010/12, they were losing 1.2% a−1. Streamflow of the Ile River and its tributaries do not exhibit statistically significant change during the vegetative period between May and September. Positive trends were observed in the Ile flow in autumn, winter, and early spring. By contrast, the calculation of the total runoff from the glacier surface (including snow and ice melt) using temperature-index method and runoff forming due to melting of multiyear ice estimated from changes in glacier volume at different time steps between the 1960s and 2010s, showed that their absolute values and their contribution to total river runoff declined since the 1980s. This change is attributed to a strong reduction in glacier area.

Changes in glacier volume on Mt. Gongga, southeastern Tibetan Plateau, based on the analysis of multi-temporal DEMs from 1966 to 2015

ABSTRACTThe accelerated retreat of glaciers and the reduction of glacier ice reserves caused by climate change can significantly affect regional water resources and hydrological cycles. Changes in glacier thickness are among the key indicators that reflect this process. We analyzed changes observed in the elevation of glacier surfaces in the Gongga Mountains (GGM) using multi-temporal Digital Elevation Models (DEMs) derived from topographic maps, SRTM, ICESat and ZY-3 data. The results showed that the mean rate of change in glacier surface altitude in the GGM was ~−26.7 ± 2.03 m (0.54 ± 0.04 m a−1) between 1966 and 2015. The mean melt rates differed over different time periods, latterly showing an accelerating trend. As a general rule, glaciers appear to be losing more volume at lower than at higher elevations. Further analysis of these results suggests that supraglacial debris coverage in the GGM promotes higher rates of mass loss.

Evaluating Glacier Volume Changes since the Little Ice Age Maximum and Consequences for Stream Flow by Integrating Models of Glacier Flow and Hydrology in the Cordillera Blanca, Peruvian Andes

Evaluating the historical contribution of the volume loss of ice to stream flow based on reconstructed volume changes through the Little Ice Age (LIA) can be directly related to the understanding of glacier-hydrology in the current epoch of rapid glacier loss that has disquieting implications for a water resource in the Cordillera Blanca in the Peruvian Andes. However, the accurate prediction of the future glacial meltwater availability for the developing regional Andean society needs more extensive quantitative estimation from long-term glacial meltwater of reconstructed glacial volume. Modeling the LIA paleoglaciers through the mid-19th century (with the most extensive recent period of mountain glacier expansion having occurred around 1850 AD) in different catchments of the Cordillera Blanca allows us to reconstruct glacier volume and its change from likely combinations of climatic control variables and time. We computed the rate and magnitude of centennial-scale glacier volume changes for glacier surfaces between the LIA and the modern era, as defined by 2011 Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Digital Elevation Model Version 2 (GDEM V2) and 2008 Light Detection and Range (LiDAR) data. The model simulation showed good agreement with the observed geomorphic data and the volume and surface area (V-S) scaling remained within the 25% error range in the reconstructed simulation. Also, we employed a recently demonstrated approach (Baraër, M. et al.) to calculate meltwater contribution to glacierized catchment runoff. The results revealed multiple peaks of both mean annual and dry season discharge that have never been shown in previous research on the same mountain range.

Estimación del Área de la Superficie y el Cambio de Volumen del Glaciar del Nevado Champará (Cordillera Blanca, Perú) a partir de las Imágenes y los Modelos de Elevación Digital del Sensor ASTER / Terra (2000-2010)

El objetivo de este trabajo es la estimación de la variación del área y volumen del glaciar del nevado Champará a partir de las imágenes y los Modelos de Elevación Digital (MED) del sensor ASTER a bordo del satélite TERRA. Se han utilizado siete imágenes ASTER (2000-05, 2003-06, 2003-06, 2006-07, 2007-08, 2009-07, 2010-05) que corresponden al periodo 2000-2010. Para estimar el área del glaciar, se ha utilizado el Índice de Nieve de Diferencia Normalizada (NDSI) y el Índice del Agua de Diferencia Normalizada (NDWI). La variación del volumen del glaciar para el periodo 2003-2010 se determinó a partir de los MED (con resolución espacial 30 m) optenidos que se generaron de las bandas 3B y 3N del sensor ASTER, teniendo en consideración las sombras y las nubes sobre el glaciar que generan incertidumbre y eliminando 1533 pixeles (que representan resoluciones espaciales de altitud mayor de 30 m), quedando con 73,035 de un total de 74,568 pixeles. En cuanto a la desglaciación episódica en el periodo 2000-2010, se ha observado una disminución del área glaciar de aproximadamente 50% con relación al año 2000, y la tendencia del ritmo de cambio episódico es de una disminución exponencial. El cambio de volumen glaciar estimado se ajusta aun comportamiento lineal con una tasa aproximada de -0.085 m3 /año para todos los pixeles del MED, y con el enmascaramiento que permite eliminar pixeles anómalos es de -0.0159 m3 /año con un coeficiente de correlación de aproximadamente R2=0.9303 y R2=0.9954, respectivamente. La disminución de altitud acumulada para el periodo 2003-2010 tiende a un comportamiento lineal de aumento negativo de hielo y nieve glaciar del nevado Champará.

Brief communication: Unabated wastage of the Juneau and Stikine icefields (southeast Alaska) in the early 21st century

Abstract. The large Juneau and Stikine icefields (Alaska) lost mass rapidly in the second part of the 20th century. Laser altimetry, gravimetry and field measurements suggest continuing mass loss in the early 21st century. However, two recent studies based on time series of Shuttle Radar Topographic Mission (SRTM) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) digital elevation models (DEMs) indicate a slowdown in mass loss after 2000. Here, the ASTER-based geodetic mass balances are recalculated carefully avoiding the use of the SRTM DEM because of the unknown penetration depth of the C-band radar signal. We find strongly negative mass balances from 2000 to 2016 (−0.68 ± 0.15 m w.e. a−1 for the Juneau Icefield and −0.83 ± 0.12 m w.e. a−1 for the Stikine Icefield), in agreement with laser altimetry, confirming that mass losses are continuing at unabated rates for both icefields. The SRTM DEM should be avoided or used very cautiously to estimate glacier volume change, especially in the North Hemisphere and over timescales of less than ∼ 20 years.

Elevation and Mass Changes of the Southern Patagonia Icefield Derived from TanDEM-X and SRTM Data

The contribution to sea level rise from Patagonian icefields is one of the largest mass losses outside the large ice sheets of Antarctica and Greenland. However, only a few studies have provided large-scale assessments in a spatially detailed way to address the reaction of individual glaciers in Patagonia and hence to better understand and explain the underlying processes. In this work, we use repeat radar interferometric measurements of the German TerraSAR-X-Add-on for Digital Elevation Measurements (TanDEM-X) satellite constellation between 2011/12 and 2016 together with the digital elevation model from the Shuttle Radar Topography Mission (SRTM) in 2000 in order to derive surface elevation and mass changes of the Southern Patagonia Icefield (SPI). Our results reveal a mass loss rate of −11.84 ± 3.3 Gt·a−1 (corresponding to 0.033 ± 0.009 mm·a−1 sea level rise) for an area of 12573 km2 in the period 2000–2015/16. This equals a specific glacier mass balance of −0.941 ± 0.19 m w.e.·a−1 for the whole SPI. These values are comparable with previous estimates since the 1970s, but a magnitude larger than mass change rates reported since the Little Ice Age. The spatial pattern reveals that not all glaciers respond similarly to changes and that various factors need to be considered in order to explain the observed changes. Our multi-temporal coverage of the southern part of the SPI (south of 50.3° S) shows that the mean elevation change rates do not vary significantly over time below the equilibrium line. However, we see indications for more positive mass balances due to possible precipitation increase in 2014 and 2015. We conclude that bi-static radar interferometry is a suitable tool to accurately measure glacier volume and mass changes in frequently cloudy regions. We recommend regular repeat TanDEM-X acquisitions to be scheduled for the maximum summer melt extent in order to minimize the effects of radar signal penetration and to increase product quality.

Effect of Topography on Subglacial Discharge and Submarine Melting During Tidewater Glacier Retreat

To first order, subglacial discharge depends on climate, which determines precipitation fluxes and glacier mass balance, and the rate of glacier volume change. For tidewater glaciers, large and rapid changes in glacier volume can occur independent of climate change due to strong glacier dynamic feedbacks. Using an idealized tidewater glacier model, we show that these feedbacks produce secular variations in subglacial discharge that are influenced by subglacial topography. Retreat along retrograde bed slopes (into deep water) results in rapid surface lowering and coincident increases in subglacial discharge. Consequently, submarine melting of glacier termini, which depends on subglacial discharge and ocean thermal forcing, also increases during retreat into deep water. Both subglacial discharge and submarine melting subsequently decrease as glacier termini retreat out of deep water and approach new steady state equilibria. In our simulations, subglacial discharge reached peaks that were 6–17% higher than preretreat values, with the highest values occurring during retreat from narrow sills, and submarine melting increased by 14% for unstratified fjords and 51% for highly stratified fjords. Our results therefore indicate that submarine melting acts in concert with iceberg calving to cause tidewater glacier termini to be unstable on retrograde beds. The full impact of submarine melting on tidewater glacier stability remains uncertain, however, due to poor understanding of the coupling between submarine melting and iceberg calving.