Search

The presented study aims to apply a method for monitoring the retreat of mountain glaciers in the Alpine region of the Swiss Alps, which consists of using optical satellite imagery and their spectral capabilities to observe snow and ice objects on the earth's surface. Satellite imagery with its big legacy archive dating back to 1972 for Landsat imagery, can be of big help to track and monitor the alpine glaciers retreat for long periods and serve as a database for modelling and predicting the glaciers retreat in the future. By combining different satellite data processing approaches, results have been obtained on the spatial distribution and dynamics of the Morteratsch Glacier over a period of 51 years. The focus of the study is to track the changes that have occurred along the positions of the ice front (terminus) in the glacier's ablation zone to obtain information about the glacier's dynamics during the study period. [ABSTRACT FROM AUTHOR]

Climate change has had a significant impact on glacier recession, particularly in the Arctic, where glacier meltwater is an important contributor to global sea-level rise. Therefore, it is important to accurately quantify glacier recession within this sensitive region, using multiple observations of glacier extent. In this study, we mapped 480 glaciers in Novaya Zemlya, Russian Arctic, using object-based image analysis applied to multispectral Landsat satellite imagery in Google Earth Engine and quantify the area changes between 1986–89 and 2019–21. The results show that in 1986–89, the total glacierized area was 22 990 ± 301 km2, in 2000–01 the area was 22 525 ± 308 km2 and by 2019–21 the glacier area reduced to 21 670 ± 292 km2, representing a total of 5.8% reduction in glacier area between 1986–89 and 2019–21. Higher glacier area loss was observed on the Barents Sea coast (7.3%) compared to the Kara (4.2%), reflecting previously observed differences in warming trends. The accuracy of the automatically generated outlines of each layer (1986–89, 2000–01 and 2019–21) was evaluated by comparing with manually corrected outlines (reference data) using random sampling, resulting in an overall accuracy estimate of between 96 and 97% compared to the reference data. This automated approach in Google Earth Engine is a promising tool for rapidly mapping glacier change that reduces the amount of time required to generate accurate glacier outlines.

We developed automated ablation stakes to measure colocated in situ changes in relative glacier-surface elevation and climatological drivers of ablation. The designs, refined over 10 years of development and deployments, implement open-source hardware and common building materials. The ablation stakes record distance to the snow/ice surface, air temperature and relative humidity every 1–15 min. Using these high-frequency data, we demonstrate that melt factors calculated using standard melt-rate vs temperature regressions converge over averaging windows of approximately 12 h or greater. Furthermore, we evaluate an integral approach to estimating temperature-index melt factors for ablation. In a test case on Glaciar Perito Moreno, Argentina, this integral approach reveals an overall positive-degree-day melt factor of 7.5 mm w.e. $^\circ$C−1 d−1. We describe four deployments with iteratively improved designs and provide a list of materials required to construct an automated ablation stake.

The geodetic mass balance of a glacier corresponds to glacier-wide volume changes, converted to mass changes using density assumptions. It is typically calculated by differencing multi-temporal digital elevation models. In this study, we show how the annual geodetic mass balance of a glacier can be derived from uncrewed aerial vehicle (UAV) data. The presented workflow is applied to two small- to medium-sized glaciers in the Kyrgyz Tien Shan (Central Asia): Bordu glacier and Sary-Tor glacier. The obtained geodetic mass balance is compared with the glaciological mass balance derived from a network of ablation stakes and snow pits. A previously calibrated mass-balance model is used to correct for the difference in acquisition dates. The results show that the determined geodetic mass balance matches closely with the glaciological mass balance. Besides, for both glaciers the geodetic mass balance does not seem to be particularly sensitive to the assumptions regarding volume-to-mass conversion. Therefore, our results demonstrate that UAVs can serve as a valuable instrument to quantify the annual geodetic mass balance and to validate the glaciological mass balance. The conventional glaciological mass-balance estimation often relies on interpolation and extrapolation methods, whereas UAVs offer the potential for direct data acquisition over the entire glacier surface.

Aldegondabreen is a relatively small (5.3 km2) land-terminating glacier located in Nordeskiöld Land of Svalbard, ~10 km southwest of Barentsburg. Cumulative mass balance during 2002–20 equalled −21.79 m w.e., which corresponds to 37% of the total mass loss. The annual mass balance (Ba) varied from −0.24 to −2.19 m w.e., while the winter mass balance (Bw) ranged between 0.36 and 0.85 m w.e. Ba and Bw were strongly correlated with the positive degree-day sums and solid precipitation amounts, respectively, measured at the Barentsburg weather station. There was also a strong correlation (r = 0.76) between Ba and Bw, which shows that winter accumulation amplifies the consecutive summer ice melt by changing the surface albedo. The trend in both observational series is not detectable because the period from 2005 to 2013 was characterized by relatively high Bw and not very negative Ba values. This was also observed on the other Svalbard glaciers, and was related to prevailing north-westerly flows over Svalbard during the summer. Therefore, the decadal periodicity of the Aldegondabreen mass balance follows general archipelago patterns that are determined by regional-scale factors. Thus, the surface mass-balance time series, which is now the longest one in the central part of the Spitsbergen Island, is representative for the archipelago.

We use the apparent horizontal shift of an englacial tephra layer outcrop to calculate local glacier mass balance on Sléttjökull, a lobe of Mýrdalsjökull in Southern Iceland. For this approach, the dipping angle of the englacial tephra layer in the glacier upstream of the outcrop and the flow velocity of the ice need to be known. An earlier investigation was expanded by the application of ground-penetrating radar, detecting the depth of the tephra along tracks with a total length of 10 km. Interpolation between the tracks enables us to derive the dipping angle of the layer along several flow lines. Together with glacier surface velocities, determined from feature tracking, we are able to estimate the local surface mass balance from the horizontal displacement of the tephra outcrop using freely available satellite imagery without additional fieldwork. The earlier local balance series was extended to the period 2014/15 to 2019/20. Although the results for the individual profiles differ slightly from each other, they show the same temporal pattern and clear variations from year to year. The results are compared to traditional mass-balance data from Hofsjökull. The two series show a good agreement in their interannual variability.

Future mass loss projections of the Greenland ice sheet require understanding of the processes at a glacier terminus, especially of iceberg calving. We present detailed and high-rate terrestrial radar interferometer observations of Eqip Sermia and Bowdoin Glacier, two outlet glaciers in Greenland with comparable dimensions and investigate iceberg calving, surface elevation, velocity, strain rates and their links to air temperature, tides and topography. The results reveal that the two glaciers exhibit very different flow and calving behaviour on different timescales. Ice flow driven by a steep surface slope with several topographic steps leads to high velocities, areas of extension and intense crevassing, which triggers frequent but small calving events independent of local velocity gradients. In contrast, ice flow under smooth surface slopes leaves the ice relatively intact, such that sporadic large-scale calving events dominate, which initiate in areas with high shearing. Flow acceleration caused by enhanced meltwater input and tidal velocity variations were observed for terminus sections close to floatation. Firmly grounded terminus sections showed no tidal signal and a weak short-term reaction to air temperature. These results demonstrate reaction timescales to external forcings from hours to months, which are, however, strongly dependent on local terminus geometry.

Climate change has had a significant impact on glacier recession, particularly in the Arctic, where glacier meltwater is an important contributor to global sea-level rise. Therefore, it is important to accurately quantify glacier recession within this sensitive region, using multiple observations of glacier extent. In this study, we mapped 480 glaciers in Novaya Zemlya, Russian Arctic, using object-based image analysis applied to multispectral Landsat satellite imagery in Google Earth Engine and quantify the area changes between 1986–89 and 2019–21. The results show that in 1986–89, the total glacierized area was 22 990 ± 301 km2, in 2000–01 the area was 22 525 ± 308 km2 and by 2019–21 the glacier area reduced to 21 670 ± 292 km2, representing a total of 5.8% reduction in glacier area between 1986–89 and 2019–21. Higher glacier area loss was observed on the Barents Sea coast (7.3%) compared to the Kara (4.2%), reflecting previously observed differences in warming trends. The accuracy of the automatically generated outlines of each layer (1986–89, 2000–01 and 2019–21) was evaluated by comparing with manually corrected outlines (reference data) using random sampling, resulting in an overall accuracy estimate of between 96 and 97% compared to the reference data. This automated approach in Google Earth Engine is a promising tool for rapidly mapping glacier change that reduces the amount of time required to generate accurate glacier outlines. [ABSTRACT FROM AUTHOR]

Close-range sensors are employed to observe glaciological processes that operate over short timescales (e.g. iceberg calving, glacial lake outburst floods, diurnal surface melting). However, under poor weather conditions optical instruments fail while the operation of radar systems below 17 GHz do not have sufficient angular resolution to map glacier surfaces in detail. This letter reviews the potential of millimetre-wave radar at 94 GHz to obtain high-resolution 3-D measurements of glaciers under most weather conditions. We discuss the theory of 94 GHz radar for glaciology studies, demonstrate its potential to map a glacier calving front and summarise future research priorities. [ABSTRACT FROM AUTHOR]

Seismic Rayleigh wave ellipticity measurements are the horizontal-to-vertical ratio of the Rayleigh wave particle motion, and are sensitive to the subsurface structure beneath a seismic station. H/V ratios measured from the ambient vibrations of the Earth are being increasingly used in glaciological applications to determine glacier and ice sheet thickness, seismic velocities and firn properties. Using the newly developed degree-of-polarisation (DOP-E) method which exploits the polarisation properties of seismic noise, we identify and extract Rayleigh waves from seismic stations in Greenland, and relate them to sea ice processes and the geology of the upper crust. Finally, we provide some suggestions for future applications of DOP-E method to gain greater insight into seasonal and long-term variability of sea ice formation and breakup as well as the monitoring of ice sheet thickness, subglacial environment and firn layers in the poles. [ABSTRACT FROM AUTHOR]

Glacier monitoring has been internationally coordinated for more than 125 years. Despite this long history, there is no authoritative answer to the popular question: ‘Which glaciers are the largest in the world?’ Here, we present the first systematic assessment of this question and identify the largest glaciers in the world – distinct from the two ice sheets in Greenland and Antarctica but including the glaciers on the Antarctic Peninsula. We identify the largest glaciers in two domains: on each of the seven geographical continents and in the 19 first-order glacier regions defined by the Global Terrestrial Network for Glaciers. Ranking glaciers by area is non-trivial. It depends on how a glacier is defined and mapped and also requires differentiating between a glacier and a glacier complex, i.e. glaciers that meet at ice divides such as ice caps and icefields. It also depends on the availability of a homogenized global glacier inventory. Using separate rankings for glaciers and glacier complexes, we find that the largest glacier complexes have areas on the order of tens of thousands of square kilometers whereas the largest glaciers are several thousands of square kilometers. The world's largest glaciers and glacier complexes are located in the Antarctic, Arctic and Patagonia.

Glacier change in Kenai Fjords National Park in southcentral Alaska affects local terrestrial, fresh water and marine ecosystems and will likely impact ecotourism. We used Landsat 4–8 imagery from 1984 through 2021 to manually map lower glacier ice margins for 19 maritime glaciers in Kenai Fjords National Park. Of these glaciers, six are tidewater, three are lake-terminating, six are land-terminating and four terminated in more than one environment throughout the study period. We used the mapped ice margins to quantify seasonal terminus position and areal change, including distinguishing between ice loss at glacier termini and along glacier margins. Overall, 13 glaciers substantially retreated (more than 2σ), 14 lost substantial area and only two underwent both net advance and area gain. The glaciers that had insubstantial length and area changes were predominantly tidewater. Cumulatively, the lower reaches of these 19 glaciers lost 42 km2 of ice, which was nearly evenly distributed between the terminus and the lateral margins. The rapid rate of glacier change and subsequent land cover changes are highly visible to visitors and locals at Kenai Fjords National Park, and this study quantifies those changes in terms of glacier length and area.

We quantify the surface elevation changes along Yala Glacier in Langtang Valley, Nepal Himalaya, since 1981 using geodetic methods to understand the recent evolution and current state of small debris-free glaciers across the region. We analyse differential global positioning system measurements and aerial stereo imagery that were acquired along Yala Glacier in 2007, 2009, 2012 and 2015 to generate digital elevation models for each calculation period. Continuous surface lowering has mainly been observed across the down-glacier area during the calculation periods, although a large degree of variability exists, with this lowering trend propagating up-glacier in recent years. The area-weighted glacier mass balances range from −0.98 ± 0.27 to −0.26 ± 0.30 m w.e. a−1 for the five calculation periods (1981–2007, 2007–2009, 2009–2012, 2012–2015 and 2007–2015). These calculated mass-balance data reveal that Yala Glacier has undergone accelerated mass loss since the late 2000s, which is consistent with the results of previous in situ measurement and remote-sensing studies.

Future mass loss projections of the Greenland ice sheet require understanding of the processes at a glacier terminus, especially of iceberg calving. We present detailed and high-rate terrestrial radar interferometer observations of Eqip Sermia and Bowdoin Glacier, two outlet glaciers in Greenland with comparable dimensions and investigate iceberg calving, surface elevation, velocity, strain rates and their links to air temperature, tides and topography. The results reveal that the two glaciers exhibit very different flow and calving behaviour on different timescales. Ice flow driven by a steep surface slope with several topographic steps leads to high velocities, areas of extension and intense crevassing, which triggers frequent but small calving events independent of local velocity gradients. In contrast, ice flow under smooth surface slopes leaves the ice relatively intact, such that sporadic large-scale calving events dominate, which initiate in areas with high shearing. Flow acceleration caused by enhanced meltwater input and tidal velocity variations were observed for terminus sections close to floatation. Firmly grounded terminus sections showed no tidal signal and a weak short-term reaction to air temperature. These results demonstrate reaction timescales to external forcings from hours to months, which are, however, strongly dependent on local terminus geometry. [ABSTRACT FROM AUTHOR]

Regional 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.

Glacier motion responds dynamically to changing meltwater inputs, but the multi-decadal response of basal sliding to climate remains poorly constrained due to its sensitivity across multiple timescales. Observational records of glacier motion provide critical benchmarks to decode processes influencing glacier dynamics, but multi-decadal records that precede satellite observation and modern warming are rare. Here we present a record of motion in the ablation zone of Saskatchewan Glacier that spans seven decades. We combine in situ and remote-sensing observations to inform a first-order glacier flow model used to estimate the relative contributions of sliding and internal deformation on dynamics. We find a significant increase in basal sliding rates between melt-seasons in the 1950s and those in the 1990s and 2010s and explore three process-based explanations for this anomalous behavior: (i) the glacier surface steepened over seven decades, maintaining flow-driving stresses despite sustained thinning; (ii) the formation of a proglacial lake after 1955 may support elevated basal water pressures; and (iii) subglacial topography may cause dynamic responses specific to Saskatchewan Glacier. Although further constraints are necessary to ascertain which processes are of greatest importance for Saskatchewan Glacier's dynamic evolution, this record provides a benchmark for studies of multi-decadal glacier dynamics.

Antennas used in glacial environmental sensor networks and reported in the last two decades have been reviewed. A link budget framework for designing such antenna systems is presented and used to design an antenna system for deployment at the Thwaites glacier, Antarctica. Design details of two left hand circularly polarized cross dipole antennas, one for englacial sensor probes and the other for supraglacial surface receivers are presented. The probe antenna is a 3D bent cross dipole that fits within a borehole of 8 cm diameter while providing a 1 dBic gain at 433 MHz in ice. The surface receiver antenna is a planar printed antenna providing a gain of 6.1 dBic with a quarter wave reflector. Both antennas provide 3 dB beamwidths of at least 50° in the xz and yz vertical planes catering for transmitter-receiver antenna misalignments caused by extended deployments. The antennas displayed good circular polarization and polarization purity traits. The 3 dB axial ratio bandwidths of both the antennas remained 54.9 %. The total efficiencies of the bent cross dipole and the surface receiver antennas were noted as 69.7 % and 86.9 % respectively. Lastly, the 433 MHz band has been validated for achieving englacial communication ranges of up to 2300 metres.

An end of summer snowline (EOSS) photographic dataset for Aotearoa New Zealand contains over four decades of equilibrium line altitude (ELA) observations for more than 50 index glaciers. This dataset provides an opportunity to create a climatological ELA reference series that has several applications. Our work screened out EOSS sites that had low temporal coverage and also removed limited observations when the official survey did not take place. Snowline data from 41 of 50 glaciers in the EOSS dataset were retained and included in a normalised master snowline series that spans 1977–2020. Application of the regionally representative normalised master snowline series in monthly and seasonally resolved climate response function analyses showed consistently strong relationships with austral warm-season temperatures for land-based stations west of the Southern Alps and the central Tasman Sea. There is a trend towards higher regional snowlines since the 1990s that has been steepening in recent decades. If contemporary decadal normalised master snowline series trends are maintained, the average Southern Alps snowline elevation will be displaced at least 200 m higher than normal by the 2025–2034 decade. More frequent extremely high snowlines are expected to drive more extreme cumulative mass-balance losses that will reduce the glacierised area of Aotearoa New Zealand.

Close-range sensors are employed to observe glaciological processes that operate over short timescales (e.g. iceberg calving, glacial lake outburst floods, diurnal surface melting). However, under poor weather conditions optical instruments fail while the operation of radar systems below 17 GHz do not have sufficient angular resolution to map glacier surfaces in detail. This letter reviews the potential of millimetre-wave radar at 94 GHz to obtain high-resolution 3-D measurements of glaciers under most weather conditions. We discuss the theory of 94 GHz radar for glaciology studies, demonstrate its potential to map a glacier calving front and summarise future research priorities.

Seismic Rayleigh wave ellipticity measurements are the horizontal-to-vertical ratio of the Rayleigh wave particle motion, and are sensitive to the subsurface structure beneath a seismic station. H/V ratios measured from the ambient vibrations of the Earth are being increasingly used in glaciological applications to determine glacier and ice sheet thickness, seismic velocities and firn properties. Using the newly developed degree-of-polarisation (DOP-E) method which exploits the polarisation properties of seismic noise, we identify and extract Rayleigh waves from seismic stations in Greenland, and relate them to sea ice processes and the geology of the upper crust. Finally, we provide some suggestions for future applications of DOP-E method to gain greater insight into seasonal and long-term variability of sea ice formation and breakup as well as the monitoring of ice sheet thickness, subglacial environment and firn layers in the poles.

We quantify the surface elevation changes along Yala Glacier in Langtang Valley, Nepal Himalaya, since 1981 using geodetic methods to understand the recent evolution and current state of small debris-free glaciers across the region. We analyse differential global positioning system measurements and aerial stereo imagery that were acquired along Yala Glacier in 2007, 2009, 2012 and 2015 to generate digital elevation models for each calculation period. Continuous surface lowering has mainly been observed across the down-glacier area during the calculation periods, although a large degree of variability exists, with this lowering trend propagating up-glacier in recent years. The area-weighted glacier mass balances range from −0.98 ± 0.27 to −0.26 ± 0.30 m w.e. a−1 for the five calculation periods (1981–2007, 2007–2009, 2009–2012, 2012–2015 and 2007–2015). These calculated mass-balance data reveal that Yala Glacier has undergone accelerated mass loss since the late 2000s, which is consistent with the results of previous in situ measurement and remote-sensing studies. [ABSTRACT FROM AUTHOR]

Glacier monitoring has been internationally coordinated for more than 125 years. Despite this long history, there is no authoritative answer to the popular question: 'Which glaciers are the largest in the world?' Here, we present the first systematic assessment of this question and identify the largest glaciers in the world – distinct from the two ice sheets in Greenland and Antarctica but including the glaciers on the Antarctic Peninsula. We identify the largest glaciers in two domains: on each of the seven geographical continents and in the 19 first-order glacier regions defined by the Global Terrestrial Network for Glaciers. Ranking glaciers by area is non-trivial. It depends on how a glacier is defined and mapped and also requires differentiating between a glacier and a glacier complex, i.e. glaciers that meet at ice divides such as ice caps and icefields. It also depends on the availability of a homogenized global glacier inventory. Using separate rankings for glaciers and glacier complexes, we find that the largest glacier complexes have areas on the order of tens of thousands of square kilometers whereas the largest glaciers are several thousands of square kilometers. The world's largest glaciers and glacier complexes are located in the Antarctic, Arctic and Patagonia. [ABSTRACT FROM AUTHOR]

Glacier change in Kenai Fjords National Park in southcentral Alaska affects local terrestrial, fresh water and marine ecosystems and will likely impact ecotourism. We used Landsat 4–8 imagery from 1984 through 2021 to manually map lower glacier ice margins for 19 maritime glaciers in Kenai Fjords National Park. Of these glaciers, six are tidewater, three are lake-terminating, six are land-terminating and four terminated in more than one environment throughout the study period. We used the mapped ice margins to quantify seasonal terminus position and areal change, including distinguishing between ice loss at glacier termini and along glacier margins. Overall, 13 glaciers substantially retreated (more than 2 σ), 14 lost substantial area and only two underwent both net advance and area gain. The glaciers that had insubstantial length and area changes were predominantly tidewater. Cumulatively, the lower reaches of these 19 glaciers lost 42 km2 of ice, which was nearly evenly distributed between the terminus and the lateral margins. The rapid rate of glacier change and subsequent land cover changes are highly visible to visitors and locals at Kenai Fjords National Park, and this study quantifies those changes in terms of glacier length and area. [ABSTRACT FROM AUTHOR]

Regional 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. [ABSTRACT FROM AUTHOR]

Glacier motion responds dynamically to changing meltwater inputs, but the multi-decadal response of basal sliding to climate remains poorly constrained due to its sensitivity across multiple timescales. Observational records of glacier motion provide critical benchmarks to decode processes influencing glacier dynamics, but multi-decadal records that precede satellite observation and modern warming are rare. Here we present a record of motion in the ablation zone of Saskatchewan Glacier that spans seven decades. We combine in situ and remote-sensing observations to inform a first-order glacier flow model used to estimate the relative contributions of sliding and internal deformation on dynamics. We find a significant increase in basal sliding rates between melt-seasons in the 1950s and those in the 1990s and 2010s and explore three process-based explanations for this anomalous behavior: (i) the glacier surface steepened over seven decades, maintaining flow-driving stresses despite sustained thinning; (ii) the formation of a proglacial lake after 1955 may support elevated basal water pressures; and (iii) subglacial topography may cause dynamic responses specific to Saskatchewan Glacier. Although further constraints are necessary to ascertain which processes are of greatest importance for Saskatchewan Glacier's dynamic evolution, this record provides a benchmark for studies of multi-decadal glacier dynamics. [ABSTRACT FROM AUTHOR]

An end of summer snowline (EOSS) photographic dataset for Aotearoa New Zealand contains over four decades of equilibrium line altitude (ELA) observations for more than 50 index glaciers. This dataset provides an opportunity to create a climatological ELA reference series that has several applications. Our work screened out EOSS sites that had low temporal coverage and also removed limited observations when the official survey did not take place. Snowline data from 41 of 50 glaciers in the EOSS dataset were retained and included in a normalised master snowline series that spans 1977–2020. Application of the regionally representative normalised master snowline series in monthly and seasonally resolved climate response function analyses showed consistently strong relationships with austral warm-season temperatures for land-based stations west of the Southern Alps and the central Tasman Sea. There is a trend towards higher regional snowlines since the 1990s that has been steepening in recent decades. If contemporary decadal normalised master snowline series trends are maintained, the average Southern Alps snowline elevation will be displaced at least 200 m higher than normal by the 2025–2034 decade. More frequent extremely high snowlines are expected to drive more extreme cumulative mass-balance losses that will reduce the glacierised area of Aotearoa New Zealand. [ABSTRACT FROM AUTHOR]

The overall trend of rapid retreat of Alpine glaciers contains considerable variability of responses at the scale of individual glaciers. As a step towards a regional assessment of glacier state that allows a detailed differentiation of single glaciers, we explore the potential of a self-organizing maps (SOM) algorithm to identify and cluster recurring patterns of thickness change at glaciers in western Austria. Using digital elevation models and glacier inventories for three time periods, we compute the frequency distribution of surface elevation change over the area of each glacier in the data set, for each period. The results of the SOM clustering show a distinct pattern shift over time: From 1969 to 1997, surface elevation change occurred at relatively uniform rates across a given glacier. Since 1997, the distribution of surface elevation change at individual glaciers has been far less uniform, indicating accelerated processes of disintegration. Tracking the evolution of individual glaciers throughout the time periods via the clusters highlights both the broader regional trend as well as glaciers that deviate from this trend, e.g. some very small, high elevation glaciers that have reverted to reduced and more uniform volume loss patterns.

The systematic monitoring of glaciers is essential to both evaluate water resource availability and better understand the effects of climate change. The increased speed of glacier changes observed in the past few years requires a more frequent update of the glacier inventories than in the past; however, the high human supervision required by the state-of-the-art techniques is discouraging their systematic application over large areas. This article proposes a novel approach to exploit the large volume of data provided by Copernicus Sentinel missions for detecting glacier outlines, including debris-covered glaciers. In detail, our method exploits the Sentinel-1 and Sentinel-2 multitemporal images to build a composite image representing the glacier conditions during the yearly maximum ablation period. The Sentinel-2 multispectral images are classified with a support vector machine and composed to a mosaic that represents the information of the maximum glacier ablation. At the same time, the Sentinel-1 time series are exploited to build a multitemporal coherence composite that represents all the snow-covered and glaciated areas together with all the moving surfaces. This information is used together with the Sentinal-2 composite to detect the debris-covered part of the glaciers. The proposed method was tested in the Central East Alps and presented an overall accuracy of 92% with respect to a reference inventory over South Tyrol and an agreement of 90% with respect to the latest glacier inventory of the Alps from Sentinel-2. The proposed approach enables to assist glacier experts in identifying glacier outlines over large areas and in short time

Seismic Rayleigh wave ellipticity measurements are the horizontal-to-vertical ratio of the Rayleigh wave particle motion, and are sensitive to the subsurface structure beneath a seismic station. H/V ratios measured from the ambient vibrations of the Earth are being increasingly used in glaciological applications to determine glacier and ice sheet thickness, seismic velocities and firn properties. Using the newly developed degree-of-polarisation (DOP-E) method which exploits the polarisation properties of seismic noise, we identify and extract Rayleigh waves from seismic stations in Greenland, and relate them to sea ice processes and the geology of the upper crust. Finally, we provide some suggestions for future applications of DOP-E method to gain greater insight into seasonal and long-term variability of sea ice formation and breakup as well as the monitoring of ice sheet thickness, subglacial environment and firn layers in the poles. [ABSTRACT FROM AUTHOR]

Glacier mass-balance observations at seasonal resolution have been performed since 1914 at two sites on Claridenfirn, Switzerland. The measurements are the longest uninterrupted records of glacier mass balance worldwide. Here, we provide a complete re-analysis of the 106-year series (1914–2020), focusing on both point and glacier-wide mass balance. The approaches to evaluate and homogenize the direct observations are described in detail. Based on conservative assumptions, average uncertainties of $\pm$0.25 m w.e. are estimated for glacier-wide mass balances at the annual scale. It is demonstrated that long-term variations in mass balance are clearly driven by melting, whereas decadal changes in accumulation are uncorrelated with mass balance and can only be relevant in short periods. Mass change of Claridenfirn is impacted by dry calving at a frontal ice cliff. Considerations of ice volume flux at a cross-profile reveal long-term variations in frontal ice loss accounting for $\sim$9% of total annual ablation on average. The effect of changes in frontal ablation mostly explains $\lt$10% of the mass-balance difference relative to the period 1960–1990, but accounts for $\sim$20% in 2010–2020. Glacier mass changes are discussed in the context of observations throughout the European Alps indicating that Claridenfirn is regionally representative.

During the last few decades, the lake-terminating glaciers in the Himalaya have receded faster than the land-terminating glaciers as proglacial lakes have exacerbated the mass loss of their host glaciers. Monitoring the impacts of glacier recession and dynamics on lake extent and water volume provides an approach to assess the mass interplay between glaciers and proglacial lakes. We describe the recession of Longbasaba Glacier and estimate the mass wastage and its contribution to the water volume of its proglacial lake. The results show that the glacier area has decreased by 3% during 1988–2018, with a more variable recession prior to 2008 than in the last decade. Longbasaba Lake has expanded by 164% in area and 237% in water volume, primarily as a result of meltwater inflow produced from surface lowering of the glacier. Over the periods 1988–2000 and 2000–18, the mass loss contributed by glacier thinning has decreased from 81 to 61% of the total mass loss, accompanied by a nearly doubled contribution from terminus retreat. With the current rate of retreat, Longbasaba glacier is expected to terminate in its proglacial lake for another four decades. The hazard risk of this lake is expected to continue to increase in the near future because of the projected continued glacier mass loss and related lake expansion.

Due to climate change, the Qinghai Hoh Xil World Heritage Site is under a significant threat from glacier melting and related glaciological hazards. A comprehensive understanding of glacier change trends will be important for launching the policy of sustainable development and for taking precautionary actions when abnormal glacier movements are recognized. Current research neglects the phenomenon in glacier movement monitoring, where the glacier variation pattern responds abnormally over different time scales. In our work, the most advanced techniques, Google Earth Engine (GEE) and offset tracking (OT), were utilized to construct and compare long-term series of glacier changes over decadal and annual scales, respectively. In total, 264 glaciers were extracted, and similar trends, namely, general degradation except one glacier expressing a fast increase in XinQing Feng (XQF), slow recession in Malan (ML), and complex and seasonal degradation in Galadandong (GLDD), were confirmed by GEE and OT. By comparing the performance of these two methods, GEE shows its merits when applied to detect the large magnitude of geomorphic and glaciological change with an average uncertainty error of approximately 8 m/yr, while OT is more suitable for detecting regular and continuous displacement with high accuracy up to 1 m/yr. Our results highlight the asynchrony of glacier variation within different time scales and the significant complementary role of integrating GEE and OT approaches for world heritage site conservation and protection.

Unoccupied Aerial Vehicles (UAVs) equipped with optical instruments are increasingly deployed in high mountain environments to investigate and monitor glacial and periglacial processes. The comparison and fusion of UAV data with airborne and terrestrial data offers the opportunity to analyse spatio-temporal changes in the mountains and to upscale findings from local UAV surveys to larger areas. However, due to the lack of gridded high-resolution data in alpine terrain, the specific challenges and uncertainties associated with the comparison and fusion of multi-temporal data from different platforms in this environment are not well known. Here we make use of UAV, airborne, and terrestrial data from four (peri)glacial alpine study sites with different topographic settings. The aim is to assess the accuracy of UAV photogrammetric products in complex terrain, to point out differences to other products, and to discuss best practices regarding the fusion of multi-temporal data. The surface geometry and characteristic geomorphological features of the four alpine sites are well captured by the UAV data, but the positional accuracies vary greatly. They range from 15 cm (root-mean-square error) for the smallest survey area (0.2 km2) with a high ground control point (GCP) density (40 GCPs km−2) to 135 cm for the largest survey area (>2.5 km2) with a lower GCP density (

Close-range (CR) L-band radiometry and quasi-simultaneous in situ snow characterizations were conducted in May 2019 at the Swiss Camp research site in the ablation zone of the western Greenland ice sheet. Snow liquid-water and its melt/refreeze are retrieved from microwave antenna temperatures measured with the ground-based L-band radiometer ELBARA-III. The emission model (EM) used in the retrieval is a two-layer configuration of the ‘L-Band Specific Microwave Emission Model of Layered Snowpack’. Consistent snow wetness retrievals were achieved from both single- and multi-angle CR observations of L-band antenna temperatures. This suggests that multi-angle observation is not a pre-requisite for snow wetness retrieval. Therefore, in addition to soil moisture and ocean salinity (SMOS) multi-angle measurements, snow wetness can be estimated from spaceborne L-band brightness temperatures measured at a single observation angle, such as from NASA's SMAP satellite. Our results provide partial validation of a recently presented snow wetness retrieval approach based on the same EM and applied over Greenland using multi-angle SMOS brightness temperatures. Agreement between measured CR antenna temperatures and SMOS brightness temperatures is found to be within the 95% confidence intervals of ELBARA-III and SMOS measurement uncertainties. Our measurements confirm the modeled response of antenna temperatures to diurnal variations of snow wetness.

The 2007–19 glaciological mass-balance series of Mera Glacier in the Everest Region, East Nepal, is reanalysed using the geodetic mass balance assessed by differencing two DEMs obtained from Pléiades stereo-images acquired in November 2012 and in October 2018. The glaciological glacier-wide annual mass balance of Mera Glacier has to be systematically decreased by 0.11 m w.e. a−1 to match the geodetic mass balance. We attribute part of the positive bias of the glaciological mass balance to an over-estimation of the accumulation above 5520 m a.s.l., likely due to a measurement network unable to capture its spatial variability. Over the period 2007–19, Mera Glacier has lost mass at a rate of −0.41 ± 0.20 m w.e. a−1, in general agreement with regional averages for the central Himalaya. We observe a succession of negative mass-balance years since 2013.

In this article, we present a method for combining image-based displacements to build time series. This method takes advantage of the redundancy of these displacements, which comes from multiple possible combinations between images. The proposed method combines common master series with a different master date to construct a single time series of relative displacements. We were able to test this method on displacements computed using images from a stereo time-lapse device recording images of the Argentière glacier during the summer and fall of 2019. Our method was compared to two other displacement aggregation strategies: a simple common master approach and the classical inversion method. In order to perform this comparison, displacement data from four permanent GPS are used. The results show that our method provides a more accurate time series of relative displacements than those obtained with the other methods. The higher accuracy of the proposed method makes the detection of seasonal variations in glacier velocity possible.

Glaciers in the Alps and several other regions in the world have experienced strong negative mass balances over the past few decades. Some of them are disappearing, undergoing exceptionally negative mass balances that impact the mean regional value, and require replacement. In this study, we analyse the geomorphometric characteristics of 46 mass-balance glaciers in the Alps and the long-term mass-balance time series for a subset of nine reference glaciers. We identify regime shifts in the mass-balance time series (when non-climatic controls started impacting) and develop a glacier vulnerability index (GVI) as a proxy for their possible future development, based on criteria such as hypsometric index, breaks in slope, thickness distribution and elevation change pattern. We found that the subset of 46 mass-balance glaciers reflects the characteristics of the total glacier sample very well and identified a region-specific variability of the mass balance. As the GVI is strongly related to cumulative glacier mass balances, it can be used as a pre-selector of future mass-balance glaciers. We conclude that measurements on rapidly shrinking glaciers should be continued as long as possible to identify regime shifts in hind-cast and better understand the impacts of climatic variability on such glaciers.

Using 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−1 observed 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−1 for 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 m3 in 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.

Deglaciation of high mountain rockwalls alters slope stability as rockwalls become more sensitive to modifications in environmental factors (e.g. seasonal temperature variations). In the past decades, increasing efforts focused on studying deglaciated Alpine rockwalls. Yet, currently deglaciating rockfaces remain unstudied. Here, we quantify surface area variations of massive ice bodies lying on high mountain rockwalls (ice aprons) in the French sector of the Mont Blanc massif between the end of the Little Ice Age (LIA) and 2018. Surface area estimates are computed from terrestrial and aerial oblique photographs via photogrammetry. This technique allows using photographs taken without scientific intent, and to tap into diverse historical or recent photographic catalogs. We derive an ice apron surface area model from precipitation records and the positive degree-days. The studied ice aprons shrank from 1854 to the 1950s, before expanding until the end of the 1990s. The beginning of the 21st century shows a decrease in surface area, leading to the complete melt of one of the studied ice aprons in 2017. Observed variations correlate with modeled surface area, suggesting strong sensitivity of ice aprons to changes in climatic variables. By studying site-specific correlations, we explore the importance of local drivers over the balance of ice aprons.

We reanalyzed mass balance records at Taku and Lemon Creek Glaciers to better understand the relative roles of hypsometry, local climate and dynamics as mass balance drivers. Over the 1946–2018 period, the cumulative mass balances diverged. Tidewater Taku Glacier advanced and gained mass at an average rate of +0.25 ± 0.28 m w.e. a–1, contrasting with retreat and mass loss of −0.60 ± 0.15 m w.e. a−1 at land-terminating Lemon Creek Glacier. The uniform influence of regional climate is demonstrated by strong correlations among annual mass balance and climate data. Regional warming trends forced similar statistically significant decreases in surface mass balance after 1989: −0.83 m w.e. a–1 at Taku Glacier and −0.81 m w.e. a–1 at Lemon Creek Glacier. Divergence in cumulative mass balance arises from differences in glacier hypsometry and local climate. Since 2013 negative mass balance and glacier-wide thinning prevailed at Taku Glacier. These changes initiated terminus retreat, which could increase dramatically if calving begins. The future mass balance trajectory of Taku Glacier hinges on dynamics, likely ending the historic dichotomy between these glaciers.

Glacial retreat in response to warming climates in the arid Xinjiang region of northwestern China directly impacts downstream water resources available for local communities. We used high-resolution satellite imagery from 1969 to 2014 to delineate spatial changes in 54 active glaciers in the upper Kaidu River Basin in the Tian Shan as well as their past expanses during the Little Ice Age (LIA). We manually delineated their boundaries based on the interpretation of glacial, geomorphic and topographic features. From the total glacier surface area, we estimated glacier volume and mass. From 1969 to 2014, glacier area decreased by 10.1 ± 1.0 km2 (relative loss of 34.2 ± 3.5%) and mass by 1.025 ± 0.108 Gt (relative loss of 43 ± 4.6%). From the LIA maximum (est. 1586 CE) to 1969, relative losses were less (25.7 ± 4.3% area loss and 33.1 ± 5.7% mass loss). Our results indicate that glacier recession is accelerating over time and that the glaciers are currently losing over 1.5 times more relative area than elsewhere in the Tian Shan. Using linear and non-linear projections, we estimate that these glaciers may disappear between 2050 and 2150 CE if climatic warming continues at the same pace.

Individual glacier changes are still poorly documented in the Georgian Caucasus. In this paper, the change of Chalaati and Zopkhito glaciers in Georgian Caucasus has been studied between 1960 and 2014. Glacier geometries are reconstructed from archival topographic maps, Corona and Landsat images, along with modern field surveys. For the first time in the Georgian Caucasus aerial photogrammetric survey of both glacier termini was performed (2014) using a drone or Unmanned Aerial Vehicle, where high-resolution orthomosaics and digital elevation models were produced. We show that both glaciers have experienced area loss since 1960: 16.2±4.9 per cent for Chalaati Glacier and 14.6±5.1 per cent for Zopkhito Glacier with corresponding respective terminus retreat by ~675 m and ~720 m. These were accompanied by a rise in the equilibrium line altitudes of ~35 m and ~30 m, respectively. The glacier changes are a response to regional warming in surface air temperature over the last half century. We used a long-term temperature record from the town of Mestia and short-term meteorological observations at Chalaati and Zopkhito glaciers to estimate a longer-term air temperature record for both glaciers. This analysis suggests an increase in the duration of the melt season over the 54-year period, indicating the importance of summertime air temperature trends in controlling glacier loss in the Georgian Caucasus. We also observed supra-glacial debris cover increase for both glaciers over the last half century: from 6.16±6.9 per cent to 8.01±6.8 per cent for Chalaati Glacier and from 2.80±6.3 per cent to 8.53±5.7 per cent for Zopkhito Glacier.

Stationary ice-penetrating radar (sIPR) systems can be used to monitor temporal changes in electromagnetically sensitive properties of glaciers and ice sheets. We describe a system intended for autonomous operation in remote glacial environments, and document its performance during deployments in cold and temperate settings. The design is patterned after an existing impulse radar system, with the addition of a fibre-optic link and timing module to control transmitter pulses, a micro-UPS (uninterruptable power supply) to prevent uncontrolled system shutdown and a customized satellite telemetry scheme. Various implementations of the sIPR were deployed on the Kaskawulsh Glacier near an ice-marginal lake in Yukon, Canada, for 44–77 days in summers 2014, 2015 and 2017. Pronounced perturbations to englacial radiostratigraphy were observed commensurate with lake filling and drainage, and are interpreted as changes in englacial water storage. Another sIPR was deployed in 2015–2016 on ice island PII-A-1-f, which originated from the Petermann Glacier in northwest Greenland. This system operated autonomously for almost a year during which changes in thickness of the ice column were clearly detected.

Herein, we present results obtained from time-lapse imagery acquired by a digital single-lens reflex camera during 2014–18 used to monitor the Planpincieux Glacier on the Italian side of the Grandes Jorasses (Mont Blanc massif). We processed the images using image cross-correlation to measure the surface kinematics of the most active lobe of the glacier that presents a high velocity and frequent ice detachments. During the monitoring, we observed two or three periods of sharp acceleration per year that culminated with large break-offs followed by analogous decelerations. Overall, we registered more than 350 failures with a volume >100 m3, of which, 14 events had volumes larger than 10 000 m3. The study identified a monotonic positive relationship between the velocity and failure volume that may be used to estimate the volume of the collapses before an event. We identified the thresholds of velocity and acceleration that characterise the activation of the speed-up periods. The study allowed the characterisation of three different instability processes that lead to the break-off of ice chunks from the glacier terminus: (i) disaggregation, (ii) slab fracture and (iii) water tunnelling failure which can be differentiated based on the rheology, the volume involved and the trigger process.