Geodetic Mass Balance of Haxilegen Glacier No. 51, Eastern Tien Shan, from 1964 to 2018

The eastern Tien Shan hosts substantial mid-latitude glaciers, but in situ glacier mass balance records are extremely sparse. Haxilegen Glacier No. 51 (eastern Tien Shan, China) is one of the very few well-measured glaciers, and comprehensive glaciological measurements were implemented from 1999 to 2011 and re-established in 2017. Mass balance of Haxilegen Glacier No. 51 (1999–2015) has recently been reported, but the mass balance record has not extended to the period before 1999. Here, we used a 1:50,000-scale topographic map and long-range terrestrial laser scanning (TLS) data to calculate the area, volume, and mass changes for Haxilegen Glacier No. 51 from 1964 to 2018. Haxilegen Glacier No. 51 lost 0.34 km2 (at a rate of 0.006 km2 a−1 or 0.42% a−1) of its area during the period 1964–2018. The glacier experienced clearly negative surface elevation changes and geodetic mass balance. Thinning occurred almost across the entire glacier surface, with a mean value of −0.43 ± 0.12 m a−1. The calculated average geodetic mass balance was −0.36 ± 0.12 m w.e. a−1. Without considering the error bounds of mass balance estimates, glacier mass loss over the past 50 years was in line with the observed and modeled mass balance (−0.37 ± 0.22 m w.e. a−1) that was published for short time intervals since 1999 but was slightly less negative than glacier mass loss in the entire eastern Tien Shan. Our results indicate that Riegl VZ®-6000 TLS can be widely used for mass balance measurements of unmonitored individual glaciers.

The Response of Brewster Glacier to Five Decades of Climate

<p>Small perturbations in climate can produce measurable changes to the size of a glacier. Documenting such changes is important for quantifying water storage changes, and understanding glacier-climate interactions. By using all available geodetic data, such as Landsat imagery, Shuttle Radar Topography Mission, GNSS and photogrammetric techniques, as well as ground penetrating radar for the construction of a bed DEM, it is found that Brewster Glacier decreased in volume from 1967 to 2017, losing ∼56% of its volume, with a period of volume increase of ∼10% from 1986 to 1997. The overall pattern of geodetic mass balance is similar to the glaciological mass balance record, however, the geodetic method tends to show more negative values by an average of ∼0.6 m w.e. Contrary to many other New Zealand glaciers, which experienced an advance from 1983 to 2008, Brewster Glacier continued to retreat by 390 m during the study period, at an average rate of 7.8 m a⁻¹, but at a significantly reduced rate of ∼2 m a⁻¹ from 1997 until 2005. By comparing the records of Brewster Glacier and Fox and Franz Josef glaciers, we explore the differences in response and reaction times resulting from glacier area-altitude distribution, and climatic setting. Furthermore, DEMs produced by this study are now available for use by a New Zealand wide glacier monitoring programme.</p>

The Response of Brewster Glacier to Five Decades of Climate

<p>Small perturbations in climate can produce measurable changes to the size of a glacier. Documenting such changes is important for quantifying water storage changes, and understanding glacier-climate interactions. By using all available geodetic data, such as Landsat imagery, Shuttle Radar Topography Mission, GNSS and photogrammetric techniques, as well as ground penetrating radar for the construction of a bed DEM, it is found that Brewster Glacier decreased in volume from 1967 to 2017, losing ∼56% of its volume, with a period of volume increase of ∼10% from 1986 to 1997. The overall pattern of geodetic mass balance is similar to the glaciological mass balance record, however, the geodetic method tends to show more negative values by an average of ∼0.6 m w.e. Contrary to many other New Zealand glaciers, which experienced an advance from 1983 to 2008, Brewster Glacier continued to retreat by 390 m during the study period, at an average rate of 7.8 m a⁻¹, but at a significantly reduced rate of ∼2 m a⁻¹ from 1997 until 2005. By comparing the records of Brewster Glacier and Fox and Franz Josef glaciers, we explore the differences in response and reaction times resulting from glacier area-altitude distribution, and climatic setting. Furthermore, DEMs produced by this study are now available for use by a New Zealand wide glacier monitoring programme.</p>

High-resolution inventory to capture glacier disintegration in the Austrian Silvretta

A new high-resolution glacier inventory captures the rapid decay of the glaciers in the Austrian Silvretta for the years 2017 and 2018. Identifying the glacier outlines offers a wide range of possible interpretations of glaciers that have evolved into small and now totally debris-covered cryogenic structures. In previous inventories, a high proportion of active bare ice allowed a clear delineation of the glacier margins even by optical imagery. In contrast, in the current state of the glacier only the patterns and amounts of volume change allow us to estimate the area of the buried glacier remnants. We mapped the glacier outlines manually based on lidar elevation models and patterns of volume change at 1 to 0.5 m spatial resolution. The vertical accuracy of the digital elevation models (DEMs) generated from six to eight lidar points per square metre is of the order of centimetres. Between 2004/2006 and 2017/2018, the 46 glaciers of the Austrian Silvretta lost −29 ± 4 % of their area and now cover 13.1 ± 0.4 km2. This is only 32 ± 2 % of their Little Ice Age (LIA) extent of 40.9 ± 4.1 km2. The area change rate increased from 0.6 %/yr (1969–2002) to −2.4 %/yr (2004/2006–2017/2018). The Sentinel-2-based glacier inventory of 2018 deviates by just 1 % of the area. The annual geodetic mass balance referring to the area at the beginning of the period showed a loss increasing from −0.2 ± 0.1 m w.e./yr (1969–2002) to −0.8 ± 0.1 m w.e./yr (2004/2006–2017/2018) with an interim peak in 2002–2004/2006 of −1.5 ± 0.7 m w.e./yr. To keep track of the buried ice and its fate and to distinguish increasing debris cover from ice loss, we recommend inventory repeat frequencies of 3 to 5 years and surface elevation data with a spatial resolution of 1 m.

Geodetic Mass Balance of the South Shetland Islands Ice Caps, Antarctica, from Differencing TanDEM-X DEMs

Although the glaciers in the Antarctic periphery currently modestly contribute to sea level rise, their contribution is projected to increase substantially until the end of the 21st century. The South Shetland Islands (SSI), located to the north of the Antarctic Peninsula, are lacking a geodetic mass balance calculation for the entire archipelago. We estimated its geodetic mass balance over a 3–4-year period within 2013–2017. Our estimation is based on remotely sensed multispectral and interferometric SAR data covering 96% of the glacierized areas of the islands considered in our study and 73% of the total glacierized area of the SSI archipelago (Elephant, Clarence, and Smith Islands were excluded due to data limitations). Our results show a close to balance, slightly negative average specific mass balance for the whole area of −0.106 ± 0.007 m w.e. a−1, representing a mass change of −238 ± 12 Mt a−1. These results are consistent with a wider scale geodetic mass balance estimation and with glaciological mass balance measurements at SSI locations for the same study period. They are also consistent with the cooling trend observed in the region between 1998 and the mid-2010s.

Mass balances of Yala and Rikha Samba glaciers, Nepal, from 2000 to 2017

The glacier mass balance is an important variable to describe the climate system and is used for various applications like water resource management or runoff modelling. The direct or glaciological method and the geodetic method are the standard methods to quantify glacier mass changes, and both methods are an integral part of international glacier monitoring strategies. In 2011, we established two glacier mass-balance programmes on Yala and Rikha Samba glaciers in the Nepal Himalaya. Here we present the methods and data of the directly measured annual mass balances for the first six mass-balance years for both glaciers from 2011/2012 to 2016/2017. For Yala Glacier we additionally present the directly measured seasonal mass balance from 2011 to 2017, as well as the mass balance from 2000 to 2012 obtained with the geodetic method. In addition, we analysed glacier length changes for both glaciers. The directly measured average annual mass-balance rates of Yala and Rikha Samba glaciers are −0.80 ± 0.28 and −0.39 ± 0.32 m w.e. a−1, respectively, from 2011 to 2017. The geodetically measured annual mass-balance rate of Yala Glacier based on digital elevation models from 2000 and 2012 is −0.74 ± 0.53 m w.e. The cumulative mass loss for the period 2011 to 2017 for Yala and Rikha Samba glaciers is −4.80 ± 0.69 and −2.34 ± 0.79 m w.e., respectively. The mass loss on Yala Glacier from 2000 to 2012 is −8.92 ± 6.33 m w.e. The winter balance of Yala Glacier is positive, and the summer balance is negative in every investigated year. The summer balance determines the annual balance. Compared to regional mean geodetic mass-balance rates in the Nepalese Himalaya, the mean mass-balance rate of Rikha Samba Glacier is in a similar range, and the mean mass-balance rate of Yala Glacier is more negative because of the small and low-lying accumulation area. During the study period, a change of Yala Glacier's surface topography has been observed with glacier thinning and downwasting. The retreat rates of Rikha Samba Glacier are higher than for Yala Glacier. From 1989 to 2013, Rikha Samba Glacier retreated 431 m (−18.0 m a−1), and from 1974 to 2016 Yala Glacier retreated 346 m (−8.2 m a−1). The data of the annual and seasonal mass balances, point mass balance, geodetic mass balance, and length changes are accessible from the World Glacier Monitoring Service (WGMS, 2021), https://doi.org/10.5904/wgms-fog-2021-05.

Evolution of the firn pack of Kaskawulsh Glacier, Yukon: meltwater effects, densification, and the development of a perennial firn aquifer

In spring 2018, two firn cores (21 and 36 m in length) were extracted from the accumulation zone of Kaskawulsh Glacier, St. Elias Mountains, Yukon. The cores were analyzed for ice layer stratigraphy and density and compared against historical measurements made in 1964 and 2006. Deep meltwater percolation and refreezing events were evident in the cores, with a total ice content of 2.33±0.26 m in the 36 m core and liquid water discovered below a depth of 34.5 m. Together with the observed ice content, surface energy balance and firn modelling indicate that Kaskawulsh Glacier firn retained about 86 % of its meltwater in the years 2005–2017. For an average surface ablation of 0.38 m w.e. yr−1 over this period, an estimated 0.28 m w.e. yr−1 refroze in the firn, 0.05 m w.e. yr−1 was retained as liquid water, and 0.05 m w.e. yr−1 drained or ran off. The refrozen meltwater is associated with a surface lowering of 0.73±0.23 m between 2005 and 2017 (i.e., surface drawdown that has no associated mass loss). The firn has become denser and more ice-rich since the 1960s and contains a perennial firn aquifer (PFA), which may have developed over the past decade. This illustrates how firn may be evolving in response to climate change in the St. Elias Mountains, provides firn density information required for geodetic mass balance calculations, and is the first documented PFA in the Yukon–Alaska region.

First geodetic mass balance estimate of the bulk of the South Shetland Islands ice caps

&lt;p&gt;The Antarctic Peninsula ice sheet is an important contributor to sea-level rise and the glaciers in its peripheral islands have a large potential to increase their contribution under a warming climate. This region has undergone a complex history of climate change during recent decades, which justifies a close monitoring of their glaciers. The South Shetland Islands (SSI) is one of the northernmost archipelagos in this region, but it is lacking a geodetic mass balance (GMB) calculation for the entire archipelago. We have estimated the GMB of the SSI over a 3-4 years period within 2013-2017 (depending on the data availability for each island). Our estimation is based on remotely-sensed multispectral and interferometric SAR data covering 96% of the glacierized areas of the islands considered in our study, and 73% of the total glacierized area of the SSI archipelago (Elephant, Clarence and Smith Islands were excluded due to overly large slopes for SAR or limited input data). Our Results show a close-to-balance overall status during the analyzed period, with specific mass balances ranging from -0.680&amp;#177;0.071 to 0.209&amp;#177;0.025 m w.e. a&lt;sup&gt;-1&lt;/sup&gt; on Low and Livingston islands, respectively. The average specific mass balance for the whole area is -0.064&amp;#177;0.015 m w.e. a&lt;sup&gt;-1&lt;/sup&gt;, representing an ice mass loss of 0.144&amp;#177;0.035 Gt a&lt;sup&gt;-1&lt;/sup&gt;. This result is consistent with the cooling trend observed in the region between 1998 and 2017, and with the mass balance estimates by the glaciological method performed in various glaciers in the AP region (and the SSI in particular).&lt;/p&gt;

Investigating the bias of TanDEM-X digital elevation models of glaciers on the Tibetan Plateau: impacting factors and potential effects on geodetic mass-balance measurements

The TanDEM-X DEM is a valuable data source for estimating glacier mass balance. However, the accuracy of TanDEM-X elevation over glaciers can be affected by microwave penetration and phase decorrelation. To investigate the bias of TanDEM-X DEMs of glaciers on the Tibetan Plateau, these DEMs were subtracted from SPOT-6 DEMs obtained around the same time at two study sites. The average bias over the studied glacier areas in West Kunlun (175.0 km2) was 2.106 ± 0.012 m in April 2014, and it was 1.523 ± 0.011 m in Geladandong (228.8 km2) in October 2013. By combining backscatter coefficients and interferometric coherence maps, we found surface decorrelation and baseline decorrelation can cause obvious bias in addition to microwave penetration. If the optical/laser data and winter TanDEM-X data were used as new and historic elevation sources for mass-balance measurements over an arbitrary observation period of 10 years, the glacier mass loss rates in West Kunlun and Geladandong would be potentially underestimated by 0.218 ± 0.016 and 0.158 ± 0.011 m w.e. a−1, respectively. The impact is therefore significant, and users should carefully treat the bias of TanDEM-X DEMs when retrieving a geodetic glacier mass balance.

Recent Increases in Winter Snowfall Provide Resilience to Very Small Glaciers in the Julian Alps, Europe

Very small glaciers (<0.5 km2) account for more than 80% of the total number of glaciers and more than 15% of the total glacier area in the European Alps. This study seeks to better understand the impact of extreme snowfall events on the resilience of very small glaciers and ice patches in the southeastern European Alps, an area with the highest mean annual precipitation in the entire Alpine chain. Mean annual precipitation here is up to 3300 mm water equivalent, and the winter snow accumulation is approximately 6.80 m at 1800 m asl averaged over the period 1979–2018. As a consequence, very small glaciers and ice/firn patches are still present in this area at rather low altitudes (1830–2340 m). We performed repeated geodetic mass balance measurements on 14 ice bodies during the period 2006–2018 and the results show an increase greater than 10% increase in ice volume over this period. This is in accordance with several extreme winter snow accumulations in the 2000s, promoting a positive mass balance in the following years. The long-term evolution of these very small glaciers and ice bodies matches well with changes in mean temperature of the ablation season linked to variability of Atlantic Multidecadal Oscillation. Nevertheless, the recent behaviour of such residual ice masses in this area where orographic precipitation represents an important component of weather amplification is somehow different to most of the Alps. We analysed synoptic meteorological conditions leading to the exceptional snowy winters in the 2000s, which appear to be related to the influence and modification of atmospheric planetary waves and Arctic Amplification, with further positive feedbacks due to change in local sea surface temperature and its interactions with low level flows and the orography. Although further summer warming is expected in the next decades, we conclude that modification of storm tracks and more frequent occurrence of extreme snowfall events during winter are crucial in ensuring the resilience of small glacial remnants in maritime alpine sectors.