Spatiotemporal patterns of High Mountain Asia's snowmelt season identified with an automated snowmelt detection algorithm, 1987–2016

Abstract

Abstract. High Mountain Asia (HMA) -encompassing the Tibetan Plateau and surrounding mountain ranges -is the primary water source for much of Asia, serving more than a billion downstream users. Many catchments receive the majority of their yearly water budget in the form of snow, which is poorly monitored by sparse in situ weather networks. Both the timing and volume of snowmelt play critical roles in downstream water provision, as many applications -such as agriculture, drinking-water generation, and hydropower -rely on consistent and predictable snowmelt runoff. Here, we examine passive microwave data across HMA with five sensors (SSMI, SSMIS, AMSR-E, AMSR2, and GPM) from 1987 to 2016 to track the timing of the snowmelt season -defined here as the time between maximum passive microwave signal separation and snow clearance. We validated our method against climate model surface temperatures, optical remote-sensing snow-cover data, and a manual control dataset (n = 2100, 3 variables at 25 locations over 28 years); our algorithm is generally accurate within 3-5 days. Using the algorithm-generated snowmelt dates, we examine the spatiotemporal patterns of the snowmelt season across HMA. The climatically short (29-year) time series, along with complex interannual snowfall variations, makes determining trends in snowmelt dates at a single point difficult.We instead identify trends in snowmelt timing by using hierarchical clustering of the passive microwave data to determine trends in self-similar regions. We make the following four key observations. (1) The end of the snowmelt season is trending almost universally earlier in HMA (negative trends). Changes in the end of the snowmelt season are generally between 2 and 8 days decade −1 over the 29-year study period (5-25 days total). The length of the snowmelt season is thus shrinking in many, though not all, regions of HMA. Some areas exhibit later peak signal separation (positive trends), but with generally smaller magnitudes than trends in snowmelt end. (2) Areas with long snowmelt periods, such as the Tibetan Plateau, show the strongest compression of the snowmelt season (negative trends). These trends are apparent regardless of the time period over which the regression is performed. (3) While trends averaged over 3 decades indicate generally earlier snowmelt seasons, data from the last 14 years (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) exhibit positive trends in many regions, such as parts of the Pamir and Kunlun Shan. Due to the short nature of the time series, it is not clear whether this change is a reversal of a long-term trend or simply interannual variability. (4) Some regions with stable or growing glaciers -such as the Karakoram and Kunlun Shan -see slightly later snowmelt seasons and longer snowmelt periods. It is likely that changes in the snowmelt regime of HMA account for some of the observed heterogeneity in glacier response to climate change. While the decadal increases in regional temperature have in ge...

The effect of Indian Summer Monsoon rainfall on surface water δD values in the central Himalaya

“…The time series data (Figure ) indicate the largest snow melt fraction of Arun main stem waters in April, where δD values drop after winter but d‐excess remains on a comparably high level. This corresponds well to other studies reporting major snow melt impact in spring as direct surface run‐off (Smith et al, ) and protracted run‐off during the ISM onset in May to June after temporal storage in the vadose zone (Buttle, ).…”

“…Our analysis revealed that the interannual ISM strength variability in 2011 and 2012 was reflected by surface water δD values of both tributary and main stem waters. Although surface water represent a complex mixture of surface run‐off, ground‐water fed base flow, and other potential water sources, such as glacial and snow melt (Andermann et al, ; Oshun, Dietrich, Dawson, & Fung, ; Racoviteanu et al, ; Smith et al, ; Wilson et al, ), our results imply that the Arun Valley and especially its southern tributaries are very sensitive to ISM δD values. The observation of surface water δD values reflecting interannual variable ISM trajectories and amount has an important consequence for the temporal resolution of potential stable isotope records being used for palaeohydrology: Changes of δD precip values in the Arun Valley due to ISM variability may theoretically be recorded annually and are rather limited by temporal recording capacities of the proxy.…”

“…For the Shillong Plateau in eastern India, Breitenbach et al () showed that local ISM δD precip and δ 18 O precip values mostly reflect the variability of Bay of Bengal surface water and the moisture transport history, for example, convective processes over the SE Arabian Sea (Lekshmy, Midhun, Ramesh, & Jani, ), with minor influence of actual rainfall amount. The southern Himalayan front is influenced by varying quantities of ISM (Bookhagen & Burbank, ; Hren, Bookhagen, Blisniuk, Booth, & Chamberlain, ; Wulf, Bookhagen, & Scherler, ; Yu et al, ), Winter Westerly Disturbances (WWD) snowfall (Sanwal et al, ; Smith, Bookhagen, & Rheinwalt, ), and glacial melt (Racoviteanu, Armstrong, & Williams, ; Wilson, Williams, Kayastha, & Racoviteanu, ). Additionally, mesoscale weather anomalies occurring along mountain ranges and plateau margins might affect δD precip and δ 18 O precip values (Blisniuk & Stern, ; Cannon, Carvalho, Jones, & Bookhagen, ; Rohrmann et al, ).…”

Clustering River Profiles to Classify Geomorphic Domains

“…Clustering algorithms have been used to group one-dimensional data in many diverse fields, including economics (e.g., Abido, 2003), computational science (e.g., March, 1983), biological science (e.g., Eisen et al, 1998;Girvan & Newman, 2002), and environmental science (e.g., Maschler et al, 2018;Rheinwalt et al, 2015;Smith et al, 2017). Many applications of one-dimensional clustering algorithms deal with the analysis of time series data, for example, where a metric such as air temperature is measured at the same time intervals at a series of different spatial locations.…”

An improved Terra–Aqua MODIS snow cover and Randolph Glacier Inventory 6.0 combined product (MOYDGL06*) for high-mountain Asia between 2002 and 2018

“…Snow cover mapping is generally crucial for areas densely populated downstream and where snowmelt dominates the discharge (Smith et al, 2017). In the topographically complex high mountains of Asia, snow covers a vast spatial extent which is difficult to measure in the field (Immerzeel et al, 2009).…”