Biogeochemistry of Sea Ice

The monitoring of glaciers in Switzerland has a long tradition, yet glacier changes during the 20th century are only known through sparse observations. Here, we estimate a halving of Swiss glacier volumes between 1931 and 2016 by mapping historical glacier elevation changes at high resolution. Our analysis relies on a terrestrial image archive known as TerrA, which covers about 86 % of the Swiss glacierised area with 21,703 images acquired during the period 1916–1947 (1931 on average). We developed a semi-automated workflow to generate digital elevation models (DEMs) from these images, resulting in a 45 % total glacier coverage. Using the geodetic method, we estimate a Swiss-wide glacier mass balance of –0.52 ± 0.09 m w.e. a−1 between 1931 and 2016. This equates to a 51.5 ± 6.1 % loss in glacier volume. We find that low elevation, high debris cover, and gently sloping glacier termini are conductive to particularly high mass losses. In addition to these glacier-specific, quasi- centennial elevation changes, we present a new inventory of glacier outlines with known timestamps and complete attributes from around 1931. The fragmented spatial coverage and temporal heterogeneity of the TerrA archive are the largest sources of uncertainty in our glacier-specific estimates, reaching up to 0.50 m w.e. a−1. We suggest that the high-resolution mapping of historic surface elevations could unlock great potentials also for research fields other than glaciology.

Abstract. Increased rates of glacier retreat and thinning need accurate local estimates of glacier elevation change to predict future changes in glacier runoff and their contribution to sea level rise. Glacier elevation change is typically derived from digital elevation models (DEMs) tied to surface change analysis from satellite imagery. Yet, the rugged topography in mountain regions can cast shadows onto glacier surfaces, making it difficult to detect local glacier elevation changes in remote areas. A rather untapped resource comprises precise, time-stamped metadata on the solar position and angle in satellite images. These data are useful for simulating shadows from a given DEM. Accordingly, any differences in shadow length between simulated and mapped shadows in satellite images could indicate a change in glacier elevation relative to the acquisition date of the DEM. We tested this hypothesis at five selected glaciers with long-term monitoring programmes. For each glacier, we projected cast shadows onto the glacier surface from freely available DEMs and compared simulated shadows to cast shadows mapped from ∼40 years of Landsat images. We validated the relative differences with geodetic measurements of glacier elevation change where these shadows occurred. We find that shadow-derived glacier elevation changes are consistent with independent photogrammetric and geodetic surveys in shaded areas. Accordingly, a shadow cast on Baltoro Glacier (the Karakoram, Pakistan) suggests no changes in elevation between 1987 and 2020, while shadows on Great Aletsch Glacier (Switzerland) point to negative thinning rates of about 1 m yr−1 in our sample. Our estimates of glacier elevation change are tied to occurrence of mountain shadows and may help complement field campaigns in regions that are difficult to access. This information can be vital to quantify possibly varying elevation-dependent changes in the accumulation or ablation zone of a given glacier. Shadow-based retrieval of glacier elevation changes hinges on the precision of the DEM as the geometry of ridges and peaks constrains the shadow that we cast on the glacier surface. Future generations of DEMs with higher resolution and accuracy will improve our method, enriching the toolbox for tracking historical glacier mass balances from satellite and aerial images.

Accurate measurements of glacier surface topography and their changes play an essential role in various glaciological studies related to glacier dynamics and mass balance. The focus of this study is on mapping glacier digital elevation model (DEM) and elevation changes in the western Qilian Mountains, northern Tibetan Plateau, by synergistically using the TanDEM-X (TDX) bistatic Interferometric Synthetic Aperture Radar (InSAR) data in 2013 and Shuttle Radar Topography Mission (SRTM) DEM in 2000. The first high-resolution and high-precision glacier DEM is derived in this region by a TDX InSAR procedure with a nonlocal (NL) filter. Validated against the Ice, Cloud, and land Elevation Satellite (ICESat) height references, the absolute height error of the TanDEM-X DEM derived with the NL filter and the Goldstein filter with the parameters investigated is, respectively, 1.493 ± 0.747 and 1.857 ± 1.709 m. Further, four combinations of differential phase method (DiffPha) and DEM differencing method (DiffDem) with Goldstein filter and NL filter are applied to estimate glacier elevation changes between 2000 and 2013. The synergistic use of the DiffPha method and the NL filter is superior to other three combinations in terms of uncertainty and noise reduction. Generally, a clear surface thinning can be found in most glacier tongue regions, the maximum value of elevation lowering up to approximately −40 m, whereas a slight thickening is detected in accumulation areas, which are in agreement with the height difference results between GPS measurements and SRTM DEM over Laohugou Glacier No.12. This study demonstrates the potential of the TanDEM-X bistatic InSAR in mapping surface topography and elevation changes of valley glaciers in the Tibetan Plateau.

The Antarctic Peninsula (AP), encompassing the ice sheet and its peripheral glaciers, is a highly dynamic component of the cryosphere, disproportionately contributing to sea level rise. However, a large spread remains between the mass changes estimated using gravimetry, altimetry and the input/output method. Among these techniques, the satellite (radar or laser) altimetry method has a resolution of, at best, 1 km, which is too coarse to resolve the complex pattern of changes in the Peninsula. Therefore, we use digital elevation models (DEMs; 30x30 m) to map elevation changes for the entire Peninsula, combining 476 DEMs derived from SPOT5-HRS satellite images (2006-2008) and 2525 strips of the Reference Elevation Model of Antarctica (2020-2022) to provide a comprehensive 14-year record. We bias-corrected each DEM using near-synchronous ICESat/-2 laser altimetry measurements.Our observations cover 70% of the AP ice sheet and 60% of its peripheral glaciers, including for regions of the Peninsula poorly studied to date and decipher a spatially complex pattern of elevation changes. After correction with different models of firn air content and solid-earth response, we find that between 2007 and 2021, the AP ice sheet lost -27 ± 9 Gt/yr while its peripheral glaciers lost -14 ± 2 Gt/yr. For the AP ice sheet, our new estimate is 4 to 5 times more negative than the one obtained in IMBIE using purely altimetry data (-6 ± 6 Gt/yr from 2006 to 2018) and in better agreement with gravimetry and the input/output method. Our study highlights the importance of resolving fine scale elevation changes of glaciers and ice sheets.

Changes in glaciers and ice caps provide some of the clearest evidence of climate change, and as such they constitute key variables for early detection strategies in global climate-related observations. These changes have impacts on global sea level fluctuations, the regional to local natural hazard situation, as well as on societies dependent on glacier meltwater. Internationally coordinated collection and publication of standardised information about ongoing glacier changes was initiated back in 1894. The compiled data sets on the global distribution and changes in glaciers and ice caps provide the backbone of the numerous scientific publications on the latest findings about surface ice on land. Since the very beginning, the compiled data has been published by the World Glacier Monitoring Service and its predecessor organisations. However, the corresponding data tables, formats and meta-data are mainly of use to specialists. It is in order to fill the gaps in access to glacier data and related background information that this publication aims to provide an illustrated global view of the available data sets related to glaciers and ice caps, their distribution around the globe, and the changes that have occurred since the maximum extents of the so-called Little Ice Age (LIA).
\nInternational glacier monitoring has produced a range of unprecedented data compilations including some 36 000 length change observations and roughly 3 400 mass balance measurements for approximately 1 800 and 230 glaciers, respectively. The observation series are drawn from around the globe; however, there is a strong bias towards the Northern Hemisphere and Europe. A first attempt to compile a world glacier inventory was made in the 1970s based mainly on aerial photographs and maps. It has resulted to date in a detailed inventory of more than 100 000 glaciers covering an area of about 240 000 km2 and in preliminary estimates, for the remaining ice cover of some 445 000 km2 for the second half of the 20th century. This inventory task continues through the present day, based mainly on satellite images.
\nThe moraines formed towards the end of the Little Ice Age, between the 17th and the second half of the 19th century, are prominent features of the landscape, and mark Holocene glacier maximum extents in many mountain ranges around the globe. From these positions, glaciers worldwide have been shrinking significantly, with strong glacier retreats in the 1940s, stable or growing conditions around the 1920s and 1970s, and again increasing rates of ice loss since the mid 1980s. However, on a time scale of decades, glaciers in various mountain ranges have shown intermittent re-advances. When looking at individual fluctuation series, one finds a high rate of variability and sometimes widely contrasting behaviour of neighbouring ice bodies. In the current scenarios of climate change, the ongoing trend of worldwide and rapid, if not accelerating, glacier shrinkage on the century time scale is most likely of a non-periodic nature, and may lead to the deglaciation of large parts of many mountain ranges in the coming decades. Such rapid environmental changes require that the international glacier monitoring efforts make use of the swiftly developing new technologies, such as remote sensing and geo-informatics, and relate them to the more traditional field observations, in order to better face the challenges of the 21st century.

Elevation changes of glaciers revealed by SfM-digital elevation models in the Khumbu region, Nepal Himalaya, from the 1970s to 2016

Greenlandic glaciers distinct from the ice sheet make up 12% of the global glacierized area and store about 10% of the global glacier ice volume (Farinotti et al., 2019). However, knowledge about recent climate change-induced volume changes of these 19,000 individual glaciers is limited. The small number of available glaciological and geodetic mass-balance observations have a limited spatial coverage, and the representativeness of these measurements for the region is largely unknown, factors which make a regional assessment of mass balance challenging. Here we use two recently released digital elevation models (DEMs) to assess glacier-wide elevation changes of 1,526 glaciers covering 3,785 km2 in west-central Greenland: The historical AeroDEM representing the surface in 1985 and a TanDEM-X composite representing 2010–2014. The results show that on average glacier surfaces lowered by about 14.0 ± 4.6 m from 1985 until 2012 or 0.5 ± 0.2 m yr−1, which is equivalent to a sample mass loss of ~45.1 ± 14.9 Gt in total or 1.7 ± 0.6 Gt yr−1. Challenges arise from the nature of the DEMs, such as large areas of data voids, fuzzy acquisition dates, and potential radar penetration. We compared several different interpolation methods to assess the best method to fill data voids and constrain unknown survey dates and the associated uncertainties with each method. The potential radar penetration is considered negligible for this assessment in view of the overall glacier changes, the length of the observation period, and the overall uncertainties. A comparison with earlier studies indicates that for glacier change assessments based on ICESat, data selection and averaging methodology strongly influences the results from these spatially limited measurements. This study promotes improved assessments of the contribution of glaciers to sea-level rise and encourages to extend geodetic glacier mass balances to all glaciers on Greenland.

Recent elevation changes of Svalbard glaciers derived from ICESat laser altimetry

Glacier observation data from major mountain regions of the world are key to improving our understanding of glacier changes: they deliver fundamental baseline information for climatological, hydrological, and hazard assessments. In many mountain ecosystems, as well as in the adjacent lowlands, glaciers play a crucial role in freshwater provision and regulation. This article first presents the state of the art on glacier monitoring and related strategies within the framework of the Global Terrestrial Network for Glaciers (GTN-G). Both in situ measurements of changes in glacier mass, volume, and length as well as remotely sensed data on glacier extents and changes over entire mountain ranges provide clear indications of climate change. Based on experiences from capacity-building activities undertaken in the Tropical Andes and Central Asia over the past years, we also review the state of the art on institutional capacity in these regions and make further recommendations for sustainable mountain development. The examples from Peru, Ecuador, Colombia, and Kyrgyzstan demonstrate that a sound understanding of measurement techniques and of the purpose of measurements is necessary for successful glacier monitoring. In addition, establishing durable institutions, capacity-building programs, and related funding is necessary to ensure that glacier monitoring is sustainable and maintained in the long term. Therefore, strengthening regional cooperation, collaborating with local scientists and institutions, and enhancing knowledge sharing and dialogue are envisaged within the GTN-G. Finally, glacier monitoring enhances the resilience of the populations that depend on water resources from glacierized mountains or that are affected by hazards related to glacier changes. We therefore suggest that glacier monitoring be included in the development of sustainable adaptation strategies in regions with glaciated mountains.

Due to climate warming, the glaciers of the Tibetan Plateau have experienced rapid mass loss and the patterns of glacier changes have exhibited high spatiotemporal heterogeneity, especially in interior areas. As the largest ice field within the Tibetan Plateau, the Puruogangri Ice Field has attracted a lot of attention from the scientific community. However, relevant studies that are based on satellite data have mainly focused on a few periods from 2000–2016. Long-term and multiperiod observations remain to be conducted. To this end, we estimated the changes in the glacier area and mass balance of the Puruogangri Ice Field over five subperiods between 1975 and 2021, based on multisource remote sensing data. Specifically, we employed KH-9 and Landsat images to estimate the area change from 1975 to 2021 using the band ratio method. Subsequently, based on KH-9 DEM, SRTM DEM, Copernicus DEM and ZY-3 DEM data, we evaluated the glacier elevation changes and mass balance over five subperiods during 1975–2021. The results showed that the total glacier area decreased from 427.44 ± 12.43 km2 to 387.87 ± 11.02 km2 from 1975 to 2021, with a decrease rate of 0.86 km2 a−1. The rate of area change at a decade timescale was −0.74 km2 a−1 (2000–2012) and −1.00 km2 a−1 (2012–2021). Furthermore, the rates at a multiyear timescale were −1.23 km2 a−1, −1.83 km2 a−1 and −0.42 km2 a−1 for 2012–2015, 2015–2017 and 2017–2021, respectively. In terms of the glacier mass balance, the region-wide results at a two-decade timescale were −0.23 ± 0.02 m w.e. a−1 for 1975–2000 and −0.29 ± 0.02 m w.e. a−1 for 2000–2021, indicating a sustained and relatively stable mass loss over the past nearly five decades. After 2000, the loss rate at a decade timescale was −0.04 ± 0.01 m w.e. a−1 for 2000–2012 and −0.17 ± 0.01 m w.e. a−1 for 2012–2021, indicating an increasing loss rate over recent decades. It was further found that the mass loss rate was −0.12 ± 0.02 m w.e. a−1 for 2012–2015, −0.03 ± 0.01 m w.e. a−1 for 2015–2017 and −0.40 ± 0.03 m w.e. a−1 for 2017–2021. These results indicated that a significant portion of the glacier mass loss mainly occurred after 2017. According to our analysis of the meteorological measurements in nearby regions, the trends of precipitation and the average annual air temperature both increased. Combining these findings with the results of the glacier changes implied that the glacier changes seemed to be more sensitive to temperature increase in this region. Overall, our results improved our understanding of the status of glacier changes and their reaction to climate change in the central Tibetan Plateau.

Ice elevation and areal changes of glaciers from the Northern Patagonia Icefield, Chile

Glaciers in the central Andes of Chile are fundamental freshwater sources for ecosystems and communities. Overall, glaciers in this region have shown continuous recession and down-wasting, but long-term glacier mass balance studies providing precise estimates of these changes are scarce. Here, we present the first long-term (1955–2013/2015), region-specific glacier elevation and mass change estimates for the Maipo River Basin, from which the densely populated metropolitan region of Chile obtains most of its freshwater supply. We calculated glacier elevation and mass changes using historical topographic maps, Shuttle Radar Topography Mission (SRTM), TerraSAR-X add-on for Digital Elevation Measurements (TanDEM-X), and airborne Light Detection and Ranging (LiDAR) digital elevation models. The results indicated a mean regional glacier mass balance of −0.12 ± 0.06 m w.e.a−1, with a total mass loss of 2.43 ± 0.26 Gt for the Maipo River Basin between 1955–2013. The most negative glacier mass balance was the Olivares sub-basin, with a mean value of −0.29 ± 0.07 m w.e.a−1. We observed spatially heterogeneous glacier elevation and mass changes between 1955 and 2000, and more negative values between 2000 and 2013, with an acceleration in ice thinning rates starting in 2010, which coincides with the severe drought. Our results provide key information to improve glaciological and hydrological projections in a region where water resources are under pressure.

The study evaluates the global glacier inventories available for the study area viz., RGI, GAMDAM and ICIMOD, with the newly generated Kashmir University Glacier Inventory (KUGI) for three Himalaya basins; Jhelum, Suru and Chenab in the north-western Himalaya, comprising of 2096 glaciers spread over an area of 3300 km2. The KUGI was prepared from the Landsat data supplemented by Digital Elevation Model, Google Earth images and limited field surveys. The KUGI comprises of 154 glaciers in the Jhelum, 328 in the Suru and 1614 in the Chenab basin, corresponding to the glacier area of 85.9 ± 11.4 km2, 487 ± 16.2 km2 and 2727 ± 90.2 km2 respectively. The investigation revealed that most of the glaciers in the study area are < 1 km2 in size, however, the glaciers in 1–5 km2 size class cover most (55.8 %) of the glacier area. Majority of the glaciers, both in terms of number and area, are at 4500–5500 m asl except in the Jhelum where the glaciers are mostly situated between 4000–5000 m asl altitude. The glaciers in the three basins mainly harbor slopes ranging from 10–30°. It was also observed that the southern aspects host more number of large-sized glaciers than the northern aspects. Comparative analyses of the inventories revealed that the GAMDAM (RAB = 0.75) and RGI (RAB = 0.73) inventories are consistent with the KUGI. However, discrepancies were observed in the debris-covered and shadowed glaciers particularly in the ICIMOD inventory. The glacier elevation changes were also estimated for glaciers in the three basins using the Tandem-X and SRTM-C DEMs from 2000 to 2012. The investigation revealed a strong control of glacial morphology, topography, and debris cover on glacier thinning. Glacier elevation change of −1.33 ± 0.8 m a−1 was observed in the Jhelum basin but a similar glacier elevation changes of −1.08 ± 0.7 m a−1 and −1.09 ± 0.8 m a−1 was observed in the Suru and Chenab basins respectively. Evaluation of the glacier inventories and assessment of glacier elevation change in the data-scarce Himalaya, reported in this article, would constitute a reliable database for research particularly in hydrology, glaciology, and climate change. The dataset is freely available at http://doi.org/10.5281/zenodo.4461799 (Romshoo et al., 2021).

Glaciers play an important role in the climate system. The elevation change of a glacier is an important parameter in studies of glacier dynamics. Only a few ground-based measurements of high mountain glaciers are available due to their remoteness, high elevation, and complex topography. The acquisition from the German TanDEM-X (TerraSAR-X add-on for Digital Elevation Measurement) SAR imaging configuration provides a reliable data sources for studying the elevation change of glaciers. In this study, the bistatic TanDEM-X data that cover the Geladandong Mountain on the Tibetan Plateau were processed with SAR interferometry technique and the elevation changes of the mountain’s glaciers during 2000–2014 were obtained. The results indicated that although distinct positive and negative elevation changes were found for different glacier tongues, the mean elevation change was about -0.14±0.26 m a-1. Geoscience Laser Altimeter System (GLAS) data were obtained for comparison and verification. The investigation using GLAS data demonstrated the efficacy of the proposed method in determining glacier elevation change. Thus, the presented approach is appropriate for monitoring glacier elevation change and it constitutes a valuable tool for studies of glacier dynamics.

Quantifying glacier elevation and mass changes is essential for understanding glacier dynamics as well as the interaction between volcanic activity and ice cover. This study investigates glacier elevation and mass changes within the Klyuchevskaya Volcanic Group (KVG) on the Kamchatka Peninsula using TanDEM-X and SRTM C-band SAR data combined with a differential SAR interferometric approach. Elevation and mass changes are assessed for the period 2000–2020, demonstrating the suitability of TanDEM-X digital elevation models for geodetic glacier analysis in volcanically active environments. Cumulative mass loss 2000-2020 amounts to −0.782 ± 0.058 Gt. For the total glacierized area of 204.15 km², an average elevation change rate of −0.347 ± 0.011 m a⁻¹ is derived, corresponding to a specific mass balance of −295 ± 23 kg m⁻² a⁻¹ for the period 2012-2020, with locally much higher losses. Marked temporal variability is observed, with strongly increased mass loss after 2015/16 (-0.528±0.014 m a-1) coinciding with intensified volcanic activity. Enhanced supraglacial debris cover following frequent and larger eruptions significantly influences glacier mass budgets, as supported by Landsat 8 Normalized Difference Snow Index analyses. Despite the absence of field-based debris thickness measurements, spatial patterns across individual glaciers highlight the critical role of volcanic debris in modulating glacier response.