The Himalayan glaciers constitute the largest mass of ice found outside of the polar regions, with an area and volume of ~ 40,000 km2 (Azam et al., 2021) and 3500 km3(Azam et al., 2021), respectively. The Himalayan glaciers, which are also aptly known as the "water towers of the world," are the principal sources of perennial rivers like the Ganges, the Indus, and the Brahmaputra, and these rivers provide water for millions of people (Bajracharya et al., 2011; Bolch et al., 2012; Bhattacharya et al., 2016; Pradhan et al., 2021). Currently, the Himalayan glaciers are in a stage of depletion as a result of climate change (Anthwal et al., 2006; Bajracharya et al., 2008; Shrestha et al., 2011; Kaushik et al., 2022), manifested through negative mass balance, dimensional shrinkage, snowline up-shift and velocity reduction (Boltch et al., 2012; Farinotti et al., 2015; Zhou et al., 2018; King et al., 2019; Dehecq et al., 2019; Shukla and Garg, 2020; Hugonnet et al., 2021).
The velocity of a glacier is a crucial factor to comprehend its mass (Dehecq et al., 2019), ice volume (Sattar et al., 2019), surge activity (Quincey et al., 2011), topography (Sam et al., 2018), and its response to climate change (Garg et al., 2019; Zhao et al., 2020; Kaushik et al., 2022). To be more specific, decreasing glacier velocity denotes a reduction in ice mass, whereas increasing glacier velocity denotes an increase in ice mass within the glacier system (Heid and Kääb, 2012a). Reduced glacier velocity causes a noticeable shift in the glacier system because it can encourage the build-up of supraglacial debris-cover and the emergence of ice cliffs and supraglacial ponds (Sakai et al., 2002; Scherler et al., 2011; Miles et al., 2016; Kneib et al., 2021).
Velocity of the glaciers can be estimated in the field by utilizing Differential Global Navigation Satellite System (DGNSS) measurements to follow the position of stakes over time. However, this strategy can be time-consuming and expensive due to the difficult topography, harsh atmosphere, and remoteness of glaciated terrain (Willis, 1995; Watson and Quincey 2015; Dehecq et al 2015). Additionally, the entire glacier is rarely represented in field-based velocity studies since they typically measure point-based surface displacements with restricted geographical and temporal recordings (Taylor et al., 2021). Stake-based surface velocity measurements in the western and central Himalayas have only been the subject of a small number of investigations (Azam et al., 2014; Vincent et al., 2016). Given this, remote sensing technology provides a suitable alternative for determining surface velocities and offers the chance to obtain perfect glacier-wide spatial coverage, even in the most remote regions of the Himalayas (Scherler et al. 2008; Copland et al. 2011; Dehecq et al. 2015; 2019; Kumar et al. 2011; Saraswat et al. 2013; Paul et al. 2017).
A technology that shows promise for determining an object's movement over time is optical satellite data-based image correlation (Singh et al., 2018; Villarroel et al., 2018). The glacier surface velocity has been calculated using a variety of correlation methods, including cross-correlation in the spatial and frequency domains, sub-pixel phase correlation, and orientation correlation (Scambos et al., 1992; Rolstad et al., 1997; Herman et al., 2011). The scientific community has successfully used numerous satellite datasets ranging from low to medium resolution images like AWiFS, MODIS, and Landsat (MSS, TM, and OLI), to determine the displacement fields of glaciers (Scherler et al., 2008; Derkacheva et al., 2020; Kumar et al., 2020; Patel et al., 2021).
Recent studies have examined the glacier velocity in the Himalayas on a regional basis and revealed a consistent slowdown (Scherler et al., 2011; Sam et al., 2018; Dehecq et al., 2019). With the aid of 220 scenes from Landsat-7 panchromatic images taken between 1999 and 2000 and Sentinel-2 panchromatic images taken between 2017 and 2018, Zhou et al. (2021) examined 36,722 glaciers in the Himalayas and analysed the flow patterns of glaciers. According to the investigation, 32% of glaciers have accelerated, 24.5% have slowed down, and 43.5% have remained unchanged (Zhou et al., 2021). Literature suggests that, in comparison to the Western Himalayas, glaciers in the Eastern and Central Himalayas move more slowly (Bolch et al., 2012; Kaushik et al., 2022). Interestingly, western Himalayan glaciers accelerated during the final decade of the 20th century followed by a steady slow down (Scherler et al., 2011; Azam et al., 2014; Sam et al., 2018; Sahu and Gupta 2019; Garg et al., 2022), however, with intra-regional heterogeneity. The distribution of the Western Himalayas over a wide range of longitudes may be responsible for this notable divergence. Nevertheless, the long-term velocity of the glaciers in the Himalayan region has been analysed in restricted number of studies (Kaushik et al., 2022; Shukla and Garg 2019; Das et al., 2022; Patel et al., 2021; Singh et al., 2018; Shukla and Garg 2020; Yellala et al., 2019).
In the western Himalaya, Das et al. (2021) estimated glaciers velocity in the Chandra Bhaga basin over four-time frames (19992-94, 2000–2002, 2009–2011 and 2013–2019). They reported that the area of crevasses and slope shift was where clean glaciers displayed their highest velocity (> 80 m/y). Additionally, the study found that surface velocity near the front of the lake-terminating glaciers (32.3 m/y) was almost twice as high as the land-terminating glaciers (12.5 m/y). Earlier, Tiwari et al. (2014) analysed the Chhota Shigri glacier's surface velocity fluctuations between 2003 and 2009 and reported that the velocity fluctuates from about 20 m/y to about 40 m/y. Likewise, Garg et al. (2017) investigated the velocity of the Sakchum, Bara Shigri, and Chhota Shigri glaciers between 2003 and 2014 and found high heterogeneity in flow patterns. Yellala et al. (2019) estimated the annual and seasonal surface velocity of the Chhota Shigri and Bara Shigri glacier. Their findings revealed that the average SIV of the glaciers was 19.25m/y and 25.62m/y, respectively between 2016 and 2018. Patel et al. (2019) employed the particle image velocimetry (PIV) technology to calculate the Chhota Shigri glacier's surface velocity. According to their research, the SIV from 2006 to 2016 varied between 22.5 and 27 m/y. Based on data from more than two time periods, Sahu and Gupta (2019) assessed the velocities of four glaciers (Bara Shigri, Chhota Shigri, Gepang Gath, and Hamta) from 2000 to 2018 and they found that the annual surface ice velocity of the glaciers decreased by approximately 34% during the 18-year period.
Above discussion demonstrates that glacier velocity studies have always been conducted at the basin scale and mostly on the quinquennial or decadal scale. The knowledge about long-term inter-annual velocity on individual glacier and its drivers remains unconfirmed. In light of this, the current study concentrates on (a) estimating annual glacier velocity changes of 21st century (2000 to 2020) on one of the debris-covered glaciers of the western Himalaya and (b) comprehends the potential controls that affect the spatial and inter-annual velocity variations of the glacier.