Three different glacier surges at a spot: what satellites observe and what not

In the Karakoram, dozens of glacier surges occurred in the past 2 decades, making the region a global hotspot. Detailed analyses of dense time series from optical and radar satellite images revealed a wide range of surge behaviour in this region: from slow advances longer than a decade at low flow velocities to short, pulse-like advances over 1 or 2 years with high velocities. In this study, we present an analysis of three currently surging glaciers in the central Karakoram: North and South Chongtar Glaciers and an unnamed glacier referred to as NN9. All three glaciers flow towards the same small region but differ strongly in surge behaviour. A full suite of satellites (e.g. Landsat, Sentinel-1 and 2, Planet, TerraSAR-X, ICESat-2) and digital elevation models (DEMs) from different sources (e.g. Shuttle Radar Topography Mission, SRTM; Satellite Pour l’Observation de la Terre, SPOT; High Mountain Asia DEM, HMA DEM) are used to (a) obtain comprehensive information about the evolution of the surges from 2000 to 2021 and (b) to compare and evaluate capabilities and limitations of the different satellite sensors for monitoring surges of relatively small glaciers in steep terrain. A strongly contrasting evolution of advance rates and flow velocities is found, though the elevation change pattern is more similar. For example, South Chongtar Glacier had short-lived advance rates above 10 km yr−1, velocities up to 30 m d−1, and surface elevations increasing by 170 m. In contrast, the neighbouring and 3-times-smaller North Chongtar Glacier had a slow and near-linear increase in advance rates (up to 500 m yr−1), flow velocities below 1 m d−1 and elevation increases up to 100 m. The even smaller glacier NN9 changed from a slow advance to a full surge within a year, reaching advance rates higher than 1 km yr−1. It seems that, despite a similar climatic setting, different surge mechanisms are at play, and a transition from one mechanism to another can occur during a single surge. The sensor inter-comparison revealed a high agreement across sensors for deriving flow velocities, but limitations are found on small and narrow glaciers in steep terrain, in particular for Sentinel-1. All investigated DEMs have the required accuracy to clearly show the volume changes during the surges, and elevations from ICESat-2 ATL03 data fit neatly to the other DEMs. We conclude that the available satellite data allow for a comprehensive observation of glacier surges from space when combining different sensors to determine the temporal evolution of length, elevation and velocity changes

Three different glacier surges at a spot: What satellites observe and what not

In the Karakoram, dozens of glacier surges occurred in the past 2 decades, making the region a global hotspot. Detailed analyses of dense time series from optical and radar satellite images revealed a wide range of surge behaviour in this region: from slow advances longer than a decade at low flow velocities to short, pulse-like advances over 1 or 2 years with high velocities. In this study, we present an analysis of three currently surging glaciers in the central Karakoram: North and South Chongtar Glaciers and an unnamed glacier referred to as NN9. All three glaciers flow towards the same small region but differ strongly in surge behaviour. A full suite of satellites (e.g. Landsat, Sentinel-1 and 2, Planet, TerraSAR-X, ICESat-2) and digital elevation models (DEMs) from different sources (e.g. Shuttle Radar Topography Mission, SRTM; Satellite Pour l'Observation de la Terre, SPOT; High Mountain Asia DEM, HMA DEM) are used to (a) obtain comprehensive information about the evolution of the surges from 2000 to 2021 and (b) to compare and evaluate capabilities and limitations of the different satellite sensors for monitoring surges of relatively small glaciers in steep terrain. A strongly contrasting evolution of advance rates and flow velocities is found, though the elevation change pattern is more similar. For example, South Chongtar Glacier had short-lived advance rates above 10 yr-1, velocities up to 30 d-1, and surface elevations increasing by 170 m. In contrast, the neighbouring and 3-times-smaller North Chongtar Glacier had a slow and near-linear increase in advance rates (up to 500 yr-1), flow velocities below 1 d-1 and elevation increases up to 100 m. The even smaller glacier NN9 changed from a slow advance to a full surge within a year, reaching advance rates higher than 1 yr-1. It seems that, despite a similar climatic setting, different surge mechanisms are at play, and a transition from one mechanism to another can occur during a single surge. The sensor inter-comparison revealed a high agreement across sensors for deriving flow velocities, but limitations are found on small and narrow glaciers in steep terrain, in particular for Sentinel-1. All investigated DEMs have the required accuracy to clearly show the volume changes during the surges, and elevations from ICESat-2 ATL03 data fit neatly to the other DEMs. We conclude that the available satellite data allow for a comprehensive observation of glacier surges from space when combining different sensors to determine the temporal evolution of length, elevation and velocity changes

Deriving glacier surface velocities from repeat optical images

The velocity of glaciers is important for many aspects in glaciology. Mass accumulated in the accumulation area is transported down to the ablation area by deformation and sliding due to the gravitational force, and hence gla­cier velocity is connected to the mass balance of glaciers. It also contributes directly to the mass balance of calving glaciers because it is an important control of the ice discharge rate for such glaciers. Changing glacier velocities is an indicator of instable glaciers, and monitoring velocity over time can make people aware of possible hazards that may arise from instable glaciers. The movement of glaciers is also important for transporting material and for eroding the landscape. The focus of this thesis is to further develop image matching within glaciology. In image matching, images from two di.erent times are compared us­ing correlation techniques to derive glacier displacement over the time period. Most studies have concentrated on using image matching to derive glacier velocities instead of developing this method further. To be able to derive the densest possible velocity grids for all glaciers in the world, image matching methods over glacier surfaces have to be explored further. So far all images that have been used to derive velocity in glaciology have been high or medium spatial resolution images. Low resolution images cover large sections in one image, and this makes them suited for investigating the velocity of large areas such as Antarctic ice shelves. We derive velocities for Antarctic ice shelves using MODIS images with a spatial resolution of 250 m to test whether these images are suited for deriving ice shelf velocity. Because the accuracy is about one fourth of a pixel, and it is possible to use images acquired several years apart due to the low surface transformation, MODIS images are well suited for deriving velocity of Antarctic ice shelves and also to monitor their changes over time. We found when comparing di.erent image matching methods over dif­ferent glacier surfaces that the most commonly used method, normalized cross-correlation, generally performs worse compared to orientation correla­tion and the matching part of the program COSI-Corr. The only situation where normalized cross-correlation outperforms the two other methods are on narrow glaciers where small window sizes are needed. COSI-Corr per­forms best overall, but orientation correlation performs almost as well. In addition orientation correlation is the only method that manages to match striped Landsat images after the failure of the Scan Line Corrector. Both orientation correlation and COSI-Corr are considered to be methods well suited for global glacier velocity mapping. Normalized cross-correlation can supplement these two methods on narrow glaciers. The effort that has been put into developing image matching in glaciology since the start of this study, both in this study and in other studies, makes it possible to derive glacier velocities over large regions, and only computer processing time hinders automatic matching of glacier velocities worldwide. Global glacier velocities can give valuable insights. We show in this thesis that it can give information about how glaciers respond to climate change. Glacier velocity of .ve regions of the world with negative mass balance is derived, and in all regions the general glacier speed is decreasing over the last decades. In addition global glacier velocities can be used to understand glacier dynamics, and predict glacier hazards. It can be tested against gla­cier inventory parameters, and it can be used to estimate erosion rates and transport times

One Decade of Glacier Mass Changes on the Tibetan Plateau Derived from Multisensoral Remote Sensing Data

The Tibetan Plateau (TP) with an average altitude of 4,500 meters above sea level is characterized by many glaciers and ice caps. Glaciers are a natural indicator for climate variability in this high mountain environment where meteorological stations are rare or non-existent. In addition, the melt water released from the Tibetan glaciers is feeding the headwaters of the major Asian river systems and contributes to the rising levels of endorheic lakes on the plateau. As many people directly rely on the glacier melt water a continuous glacier monitoring program is necessary in this region. In situ measurements of glaciers are important, but are spatial limited due to large logistical efforts, physical constrains and high costs. Remote sensing techniques can overcome this gap and are suitable to complement in situ measurements on a larger scale. In the last decade several remote sensing studies dealt with areal changes of glaciers on the TP. However, glacier area changes only provide a delayed signal to a changing climate and the amount of melt water released from the glaciers cannot be quantified. Therefore it is important to measure the glacier mass balance. In order to estimate glacier mass balances and their spatial differences on the TP, several remote sensing techniques and sensors were synthesized in this thesis. In a first study data from the Ice Cloud and Elevation Satellite (ICESat) mission were employed. ICESat was in orbit between 2003 and 2009 and carried a laser altimeter which recorded highly accurate surface elevation measurements. As in mid-latitudes these measurements are rather sparse glaciers on the TP were grouped into eight climatological homogeneous sub-regions in order to perform a statistical sound analysis of glacier elevation changes. To assess surface elevation changes of a single mountain glacier from ICESat data, an adequate spatial sampling of ICESat measurements need to be present. This is the case for the Grosser Aletschgletscher, located in the Swiss Alps which served as a test site in this thesis. In another study data from the current TanDEM-X satellite mission and from the Shuttle Radar Topography Mission (SRTM) conducted in February 2000 were employed to calculate glacier elevation changes. In a co-authored study, these estimates could be compared with glacier elevation changes obtained from the current French Pléiades satellite mission. In order to calculate glacier mass balances, the derived elevation changes were combined with assumptions about glacier area and ice density in all studies. In this thesis contrasting patterns of glacier mass changes were found on the TP. With an ICESat derived estimate of -15.6±10.1 Gt/a between 2003 and 2009 the average glacier mass balance on the TP was clearly negative. However, some glaciers in the central and north-western part of the TP showed a neutral mass balance or a slightly positive anomaly which was also confirmed by data from the current TanDEM-X satellite mission. A possible explanation of this anomaly in mass balance could be a compensation of the temperature driven glacier melt due to an increase in precipitation

Mapping Samudra Tapu glacier:A holistic approach utilizing radar and optical remote sensing data for glacier radar facies mapping and velocity estimation

Himalayan glaciers have shown more sensitivity and visible changes to the climate change and global warming in the last 150 years. The highly rugged topography and inaccessible remote areas makes satellite images as the most appropriate source of information retrieval. We performed remote sensing based glacier change study for Samudra Tapu glacier, located in the Chandra basin of North-West Himalaya. In the present study, the capabilities of both optical and microwave remote sensing data was analysed in glacier change study in terms of its coverage, shift in equilibrium line altitude (ELA) and surface velocity over a period from 2000 to 2021. Multi Sensor (RISAT-1, Sentinel-1) time series of C-band SAR data along with a object oriented classification technique were used to identify different glacier facies such as percolation facies, icefalls, bare ice facies, refreeze snow and supraglacial debris. These classified maps were also used to detect the snow line and firn line along with ELA, aided with elevation information from digital elevation model (DEM). It was identified that more than 50 % of the total glacier area still lies into accumulation region. Further, we estimated the glacier surface velocity using Differential Interferometric Synthetic Aperture Radar (DInSAR) technique using European Remote Sensing Satellite (ERS-1/2) tandem data of 1996. High value of coherence was observed from the SAR return signal for one-day temporal difference. A mean velocity of 17–24 cm/day was found for the months of March and May 1996, highest flow rates were seen in the high accumulation area located in the Eastern and Southern Aspect of glacier. Spatial analysis of velocity patterns with respect to slope and aspect show that high rates of flow was found in southern slopes and movement rates generally increase with increase in slope. Feature tracking approach was used to estimate the glacier flow for long term and seasonal basis using optical and SAR datasets (IRS-1C, 1D PAN, Landsat-7, 8 PAN, and TANDEM-x) during 1999–2020 period. The results suggest that glacier flow varies with season, i.e., high velocity during spring-summer season, as compared to late summer or winter and, the rate of ice flow changes over the years. The mean glacier velocity reduced to 49.5 m/year during 2013–2020 time, as compared to 67.67 m/year during 1999–2003 time. These results of reducing glacier velocity and changing snow line altitude indicates enhanced glacier's melt rate and overall negative mass balance for Smudra tapu glacier.</p

Glacier area change (1993–2019) and its relationship to debris cover, proglacial lakes, and morphological parameters in the Chandra-Bhaga Basin, Western Himalaya, India

Peer reviewe

Overall recession and mass budget of Gangotri Glacier, Garhwal Himalayas, from 1965 to 2015 using remote sensing data

Thinning rates for the debris-covered Gangotri Glacier and its tributary glaciers during the period 1968-2014, length variation and area vacated at the snout from 1965 to 2015, and seasonal variation of ice-surface velocity for the last two decades have been investigated in this study. It was found that the mass loss of Gangotri and its tributary glaciers was slightly less than those reported for other debris-covered glaciers in the Himalayan regions. The average velocity during 2006-14 decreased by ~6.7% as compared with that during 1993-2006. The debris-covered area of the main trunk of Gangotri Glacier increased significantly from 1965 until 2015 with the maximum rate of increase (0.8 ± 0.2 km2 a-1) during 2006-15. The retreat (~9.0 ± 3.5 m a-1) was less in recent years (2006-2015) but the down-wasting (0.34 ± 0.2 m a-1) in the same period (2006-2014) was higher than that (0.20 ± 0.1 m a-1) during 1968-2006. The study reinforced the established fact that the glacier length change is a delayed response to climate change and, in addition, is affected by debris cover, whereas glacier mass balance is a more direct and immediate response. Therefore, it is recommended to study the glacier mass balance and not only the glacier extent, to conclude about a glacier's response to climate change.Publisher PDFPeer reviewe

Evolution of Glacial and High-Altitude Lakes in the Sikkim, Eastern Himalaya Over the Past Four Decades (1975–2017)

Global climate change is significantly triggering the dynamic evolution of high-mountain lakes which may pose a serious threat to downstream areas, warranting their systematic and regular monitoring. This study presents the first temporal inventory of glacial and high-altitude lakes in the Sikkim, Eastern Himalaya for four points in time i.e., 1975, 1991, 2000, and 2017 using Hexagon, TM, ETM+, and OLI images, respectively. First, a baseline data was generated for the year 2000 and then the multi-temporal lake changes were assessed. The annual mapping of SGLs was also performed for four consecutive years (2014–2017) to analyze their nature and occurrence pattern. The results show an existence of 463 glacial and high-altitude lakes (&gt; 0.003 km2) in 2000 which were grouped into four classes: supraglacial (SGL; 50) pro/peri glacial lake in contact with glacier (PGLC; 35), pro/peri glacial lake away from glacier (PGLA; 112) and other lakes (OL; 266). The mean size of lakes is 0.06 km2 and about 87% lakes have area &lt; 0.1 km2. The number of lakes increased (by 9%) from 425 in 1975 to 466 in 2017, accompanied by a rapid areal expansion from 25.17 ± 1.90 to 31.24 ± 2.36 km2 (24%). The maximum expansion in number (106%) and area (138%) was observed in SGLs, followed by PGLCs (number: 34%; area: 90%). Contrarily no significant change was found in other lakes. The annual SGL mapping reveals that their number (area) increased from 81 (543,153 m2) to 96 (840,973 m2) between 2014 and 2017. Occurrence pattern of SGLs shows that maximum number of lakes (&gt; 80%) are persistent in nature, followed by drain-out (15–20%) and recurring type lakes (7–8%). The new-formed lakes (9–17%) were consistently noticed in all the years (2014–2017). The results of this study underline that regional climate is accelerating the cryosphere thawing and if the current trend continues, further glacier melting will likely occur. Therefore, formation of new lakes and expansion of existing lakes is expected in the study area leading to increase in potential of glacial lake outburst floods. Thereby, persistent attention should be paid to the influences of climatic change in the region

Four decades of glacier variations at Muztagh Ata (eastern Pamir) : A multi-sensor study including Hexagon KH-9 and Pléiades data

Previous in situ measurements have indicated a slight mass gain at Muztagh Ata in the eastern Pamir, contrary to the global trend. We extend these measurements both in space and time by using remote sensing data and present four decades of glacier variations in the entire mountain massif. Geodetic mass balances and area changes were determined at glacier scale from stereo satellite imagery and derived digital elevation models (DEMs). This includes Hexagon KH-9 (year 1973), ALOS-PRISM (2009), Pléiades (2013) and Landsat 7 ETMC data in conjunction with the SRTM-3 DEM (2000). In addition, surface velocities of Kekesayi Glacier, the largest glacier at Muztagh Ata, were derived from amplitude tracking of TerraSARX images (2011). Locally, we observed strong spatial and temporal glacier variations during the last four decades, which were, however, on average not significant for the entire massif. Some south-west-exposed glaciers fluctuated or advanced, while glaciers with other aspects rather experienced continuous shrinkage. Several glaciers such as Kekesayi indicate no measurable change at their frontal position, but clear down-wasting despite mostly thick debris coverage at low altitudes. The surface velocity of this debriscovered glacier reach up to 20 cm per day, but the lowest part of the tongue appears to be stagnant. The low velocity or even stagnancy at the tongue is likely one reason for the down-wasting. On average, the glaciers showed a small, insignificant shrinkage from 274.3 ± 10.6 km2 in 1973 to 272.7 ± 1.0 km2 in 2013 (-0.02 ± 0.1%a-1). Average mass changes in the range of -0.03 ± 0.33mw.e. a-1 (1973-2009) to -0.01 ± 0.30mw.e. a-1 (1973-2013) reveal nearly balanced budgets for the last 40 years. Indications of slightly positive rates after 1999 (+0.04 ± 0.27 mw.e. a-1) are not significant, but confirmed by measurements in the field.Publisher PDFPeer reviewe