Next Article in Journal
Disinfection Methods for Swimming Pool Water: Byproduct Formation and Control
Previous Article in Journal
Variability of Short-Term Diel Water Temperature Amplitudes in a Mountain Lake
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Recent Glacier Mass Balance and Area Changes from DEMs and Landsat Images in Upper Reach of Shule River Basin, Northeastern Edge of Tibetan Plateau during 2000 to 2015

1
Shaanxi Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Northwest University, Xi’an 710127, China
2
College of Urban and Environmental Science, Northwest University, Xi’an 710127, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(6), 796; https://doi.org/10.3390/w10060796
Submission received: 8 May 2018 / Revised: 14 June 2018 / Accepted: 14 June 2018 / Published: 16 June 2018
(This article belongs to the Section Hydrology)

Abstract

:
Glacier changes in the Upper Reach of the Shule River Basin (URSRB) serve as a good indicator of climate change in the western part of the Qilian Mountains, located on the northeastern edge of the Tibetan Plateau. However, information on recent glacier changes in the URSRB is limited. In this study, the changes in ice surface elevation were determined using geodetic methods based on digital elevation models (DEMs) derived from the Shuttle Radar Topography Mission (SRTM) (2000), and from pairs of Third Resources Satellite (ZY-3) of China (taken around 2013). In addition, glacier area changes from 2000–2015, were derived from Landsat TM/ETM+/OLI images. The results suggest that 478 glaciers with an area of 375.1 ± 2.68 km2 remained in the URSRB in 2015. Ice cover diminished by 57.5 ± 2.68 km2 (11.9 ± 0.60%), or 0.79 ± 0.04% a−1 and 35 small glaciers disappeared from 2000 to 2015 in the URSRB. The most pronounced glacier shrinkage occurred during 2004 to 2009. The average ice surface elevation of the URSRB from 1999 to 2013 reduced by about 4.98 ± 0.6 m, which is equal to a mass loss of 0.383 ± 0.046 m·a−1. This reduction indicates that the ice storage loss has accelerated since 1999, compared to a mass loss of 0.21 ± 0.04 m·a−1 around Tuanjiefeng from 1966 to 1999, as reported by Xu et al. (2013).

1. Introduction

Mountain glaciers are a critical freshwater resource for populations inhabiting the downstream region or valleys at the base of mountains [1,2]. Worldwide glacier changes and the associated change in future runoff raise major concerns for the sustainability of global water resources [3]. The meltwater from glaciers on the Tibetan Plateau (TP), which is called “Asia’s Water Tower”, is extremely important to the vitality of large rivers, such as the Indus and Brahmaputra Rivers [1]. Therefore, reduction in meltwater can be a major threat to the food security of an estimated 60 million people who live in the downstream regions [1,4]. The glacier mass balance is a very useful indicator of climate variability [5], as it provides information about glacier changes in high Asia as an important issue in cryospheric and climate change studies [6].
Due to recent climate warming [7], many mountain glaciers in Asia have shrunk in area and lost mass, such as glaciers of the Tianshan Mountains [8,9,10,11], Qiangtang Plateau [12], Aksu-Tarim catchment [13], central Himalaya [14,15], and Qilian Mountains [16]. However, some stabilization or positive mass balance has been reported in the central Karakoram Range [17], including the eastern Pamir Mountains [18]. Glacier mass balance response to climate change, its impact on meltwater and river runoff, and the risk of glacier-related disasters in different regions have received many concerns from both decision-makers and local communities that are affected.
The Upper Reach of the Shule River Basin (URSRB) is located in the western Qilian Mountains, which is at the northeastern edge of the Tibetan Plateau (TP). The Hexi Corridor region, which includes Dunhuang city, is fed by runoff from the Qilian Mountains, and is significant for trade and cultural exchanges in Eurasia [19]. The yearly average annual glacial meltwater from the Qilian Mountains contributes about 30% of river runoff in the URSRB [20]. Thus, glaciers have important roles in augmenting river runoff during the arid years, with less spring and summer river discharge [21], which is the main water source for populations that live downstream of the Shule River Basin, and for agriculture irrigation. Thus, understanding and predicting changes in glaciers of the URSRB is highly significant to manage water resources for the whole Shule River Basin.
Although glaciological, hydrological, and geodetic methods have been used to determine the mass balance of a glacier [22], the high altitude and harsh climatic conditions of the URSRB make fieldwork very difficult. Fortunately, satellite remote sensing has become a promising alternative, even in remote mountainous terrain like the URSRB, to assess glaciers [23]. In this study, changes in the ice surface elevation during different periods were obtained using digital elevation models (DEMs). Currently, the main sources of DEMs include topographic maps, Shuttle Radar Topography Mission (SRTM), and optical remote sensing stereo image pairs, which offer advantages of short acquisition time, wide coverage, low cost, and sufficient data.
In a similar study, Pieczonka et al. [13] investigated changes in glacier storage volume in the southern slope area of Tomor Peak in the Tianshan Mountains during the periods 1976–1999, 1999–2009, and 1976–2009, by analyzing Hexagon KH-9 data, SPOT 5 stereoscopic image pairs, and SRTM DEM. Other research has been conducted in the western Kunlun Mountains [24], central Tianshan Mountains [11,25], Bangong Co basin [26], and Altai region [27]. For instance, Zhang et al. [18] investigated glacier surface elevation changes and surging of Karayaylak glacier in the Pamir Mountains from 2013 to 2015 using ASTER image pairs.
Some studies have also investigated glacier changes in the Qilian Mountains, such as Du et al. [28], who reported that the terminus position of Laohugou glacier No. 12 retreated 100 m during 1957–1976, and 140.12 m during 1985–2005. It was further determined that the ice surface elevation of the latter glacier decreased by 18.6 ± 5.4 m between 1957 and 2007, by comparing topographic maps and measurements of real-time kinematic GPS (global positioning system) [29]. Typical glacier changes have also been analyzed at the Yanglong River Basin [30,31] and Ningchan glacier No. 1 [32] by kinematic GPS, topographic map, and stereoscopic image pairs. Xu et al. [33] studied the changes of glacier area and ice elevation of glaciers surrounding Tuanjiefeng from the 1960s to 2000 by comparing DEMs derived from topographic maps and SRTM DEMs. Sun et al. [19] pointed out that glacier area and volume decreased by 20.88% and 20.26%, respectively, in the Qilian Mountains in the past half-century, based on the revised First and Second Chinese Glacier Inventory.
Despite the numerous studies conducted on glacier change in China, there remains minimal understanding of the recent changes in glacier area and ice surface elevation in the URSRB. Therefore, in this study, recent changes in glacier ice surface elevation in the URSRB were investigated from DEMs, including SRTM DEM and ZiYuan-3 surveying satellite (ZY-3) stereo pairs. Glacier area changes were studied using Landsat satellite image series.

2. Data and Methods

2.1. Study Area

The Shule River originates in Mulligan Ridge between Shule Nanshan and Taole Nanshan in the Qilian Mountains. The URSRB sits above Changma hydrological station (96.6°–99.0° E, 38.2°–40.0° N) and has an area of about 1.14 × 104 km2. The URSRB is characterized by the continental climate of the plateau and an elevation between 2081–5764 m a.s.l. (Figure 1). The yearly annual average air temperature is about −4 °C, and the average annual precipitation is about 378.4 mm. Precipitation mainly occurs during the glacier melting season (May–September), accounting for more than 90% of the annual precipitation [34].

2.2. Data

Satellite images from Landsat 5, 7, and 8 were used to investigate the changes in the glacier distribution and area. Landsat 5 Thematic Mapper (TM) images contained seven bands, of which the spatial resolution of band 1–5 and band 7 is 30 m. Landsat 7 Enhanced Thematic Mapper (ETM+) images included 8 bands, and the spatial resolutions of bands 1–5, 7, and 8 were 30, 30, and 15 m, respectively. Landsat 8 Operational Land Imager (OLI) images consisted of nine bands, including 8 bands with a spatial resolution of 30 m and one pan band of 15 m. All satellite images were derived from the geospatial data cloud platform (http://www.gscloud.cn/) of the Computer Network Information Center (CNIC) of the Chinese Academy of Sciences. All products were Level L1T and were calibrated by radiation correction, geometry correction by ground control points, and terrain correction by the digital elevation model (DEM). All selected Landsat images had cloud cover less than 15%. Additional information about the images are listed in Table 1. All Landsat data on glacier area changes from 2000 to 2015 were extracted between July and September, mainly in August, to reduce the effect of seasonal snow cover in the images.
SRTM was launched at the Kennedy Space Center on 11 February 2000, which used microwaves to build a global DEM. The average elevation accuracy of the SRTM DEM is better than 16 m (90% confidence interval), but the accuracy is closely related to topographic factors, such as surface slope and slope direction [35]. The unfilled version SRTM DEM (containing data voids) was obtained from interferometry of the C-band of a microwave sensor and post-processing with 90 m resolution. Information from the SRTM DEM can be found on the website http://glcf.umd.edu/data/srtm. The X-band SAR system carried with the SRTM C-band SAR system had a swath width of 45 km, leaving large data gaps in the derived X-band DEM [36]. SRTM X-band DEMs only cover 21% of URSRB glaciers, which can be downloaded from DLR multi-mission web portal https://geoservice.dlr.de/egp/. The SRTM DEM and X-band DEM were resampled with 30 m resolution by a bilinear interpolation method.
The ZY-3, launched on 9 January 2012, is China’s first civilian high-resolution stereo mapping satellite [37]. ZY-3 contains three panchromatic cameras that are organized in a triple-linear array imaging geometry, and positioned front-facing at 22°, ground-facing, and rear-facing at 22°. The ground pixel resolutions of the front-facing and rear-facing cameras are about 3.5 m, and that of the ground-facing camera is about 2.1 m. We selected ZY-3 images that showed the least cloud cover and that were acquired in the same season. The information on the 5 selected ZY-3 images is list in Table 2. A small part of the URSRB was not imaged properly by the ZY-3 stereo pair cameras, due to the influence of cloud cover. ZY-3 DEMs were resampled to 30 m in comparison to SRTM DEMs.
There are 12 national meteorological stations in and around the URSRB, which have obtained observation data from 1970 to 2013. Some precipitation stations with T200B and total rain gauges, are set up in alpine regions, which have acquired observed precipitation data after 2008. All precipitation data observed by 40 total rain gauges at elevations ranging from 1139 to 4156 m a.s.l. were used to obtain the precipitation gradient. The monthly precipitation gradient of 12.0 mm/100 m and air temperature lapse rate were calculated using linear regression between meteorological and elevation data [38,39], and the regional averages of the annual average air temperature and annual precipitation were used to analyze the glacier response to long-term climate changes in the URSRB.

2.3. Method

2.3.1. Glacier Area Extraction

A semi-automated approach, using the band ratio between different Landsat bands and corrected by visual interpretation, was applied to outline the glaciers. The band ratio of the Landsat TM to Landsat ETM+ with TM3/TM5 band was used. The extraction threshold of the glacier was 2.1 [40]. Due to the Band 1 of Landsat 8 OLI image is new, all other bands were shifted compared to the older Landsat series, and the band ratio between Red and SWIR1 was shifted to TM4/TM6 for Landsat 8 OLI with a glacier extraction threshold of 2 [41]. To eliminate ice patches less than 0.01 km2 in area from the images, a 3 × 3 median filter was first applied to eliminate isolated pixels [42], and then “Eliminate” tool with the area 0.01 km2 in ARCGIS software was used.
The derived glacier polygons were checked manually against images from adjacent years with little or no snow and cloud coverage, to discriminate supraglacial lakes, seasonal snow, surface boulders, and debris-covered ice. The glacier outlines from all Landsat images were combined to produce the glacier polygons of five periods (2000/2001, 2004/2006, 2008/2009, 2012/2013, 2014/2015).
The distribution of mountain ridge was derived from the elevation of DEMs according the method of Guo et al. [43], and then, the boundary of each single glacier was determined by dividing the glacier polygon by ridge. Finally, the outlines of each glacier were manually compared using pseudo-color composed of bands 4/3/2 from Landsat images and glacier boundaries from the Second Chinese Glacier Inventory(SCGI), which contains glacier outlines from 2006 to 2010 [43].
The identification number (ID) of each glacier was obtained from the Chinese Glacier Inventory (CGI). The area changes of each single glacier were calculated by comparing the area in the five periods, and CGI and SCGI with the same ID.

2.3.2. Glacier Ice Surface Elevation Generation

The stereo pair images composed by rear-facing and ground-facing sensors were used to generate DEMs, and the stereo pair images composed by front-facing and ground-facing sensors were taken as assist images. In the process of extracting DEMs from the ZY-3 stereo pair images, the rational polynomial coefficients (RPCs) from the image files were automatically identified by ENVI software. The same homologous points on both images were manually identified and combined, which included about 80 to 100 homologous points on each pair of images, then the selection stopped when the largest y-parallax was less than 2°. Both the RPCs and selected homologous points were used for map orthorectification. The average maximum y-parallax of the five ZY-3 image pairs was about 1.2 degrees, suggesting that the orthorectification has high precision. Then, the absolute elevation was obtained by adding the measured control points in the URSRB.
The absolute orientation parameters of high-resolution satellite stereo image pairs, such as ZY-3, can be calculated from precise ephemeris data and sensor model’s parameters without ground control points, to obtain a high-precision DEM [44]. The dynamic exterior systematic errors (including camera installation errors) and static interior distortion of ZY-3 images were eliminated using 1:2000 digital orthophotomaps and DEMs of the Dengfeng (Henan) and Tianjin areas of China as control data. The plan and height accuracy of ZY-3 DEMs without ground control points were found to be better than 3 and 2 m over the Anping area in Hebei Province and Zhaodong area in Heilongjiang Province, respectively [45]. Our study also proves that the ZY-3 can obtain good quality DEMs.
For changes in glacier surface elevation, differences between the ZY-3 DEM and SRTM DEM were used to construct a map of the ice surface elevation change. Relative horizontal and vertical distortions between the two datasets were corrected with statistical approaches based on the relationship between elevation difference, slope, and aspect [46]. Elevation differences in off-glacier regions were used to analyze the consistency of the ZY-3 DEM and SRTM DEM.
We first compared the spatial patterns of elevation differences with that of a greyscale hill shadow map by visual interpretation, where consistency in the off-glacier regions indicates that there is no obvious horizontal spatial mismatch between DEMs [47]. We determined that the relationship between the elevation difference, which divides the tangent of the slope of terrain and the cosinusoidal of the terrain aspect [47], was far from ideal. In addition, a close relationship was found between the elevation difference and terrain aspect (Figure 2). Since the 5 ZY-3 satellites have different orbits, the relationships between the elevation difference and terrain aspect varied for different images, and the following corrections, Equations (1) to (5), were used for the elevation differences between the 5 ZY-3 DEM of a-e in Table 2 and SRTM DEM, respectively.
y = ( 14.74335 ) × s i n [ π × ( x 170.53084 ) 188.27206 ] 65.10709
y = 19.7106 × s i n [ π × ( x 148.70926 ) 175.55774 ] 30.15101
y = 8.05073 × s i n [ π × ( x 161.30038 ) 185.23823 ] 53.56976
y = 14.53019 × s i n [ π × ( x 215.71857 ) 186.57741 ] 22.05558
y = 8.62003 × s i n [ π × ( x 127.57363 ) 151.02383 ] 29.8892
where y is the elevation difference between SRTM DEM and ZY-3 DEM in off-glacier regions, and x is terrain aspect derived from SRTM DEM.
After co-registration, the histogram statistics showed that median elevation difference in off-glacier regions changed from −64.7, −33.6, −53.6, −22.9, and −28.4 m (left of Figure 3) before co-registration to 0.00, 0.00, 0.00, 0.00, and 0.13 m after co-registration (right of Figure 3) for a, b, c, d, and e images, respectively, and the elevation difference was concentrated between −5 and 5 m. It is concluded that elevation difference in off-glacier regions have been stabilized, and that the pre-processed DEMs were suitable for the estimation of changes in glaciers mass balance. Then, Equations (3) to (7) were applied to the glacierized regions in each image.
Considering the penetrating effect of C band of radar microwaves on snow, SRTM DEM is thought to reflect the glacier ice surface elevation at the end of melting season in October 1999 [33]. The penetration depth of the SRTM C-band radar beam into snow and ice must be considered when assessing glacier elevation changes [13,48,49]. The penetration depth can range from 0 to 10 m depending on snow temperature, density, and water content [48].
Kaab et al. [50] developed a first-order method to correct regionally-based elevation change measurements for the average penetration depth of the C-band signal into ice and snow. In this work, we used an average penetrate depth of 2.4 m in the eastern Karakoram region [50], which has a similar climate to URSRB due to the influence of Westerly winds.

2.3.3. Ice Mass Change

The change in glacier ice volume can be obtained according to the elevation change in glacier surface, which considers glacier area and elevation change in ice surface, and is calculated by
Δ V = S × Δ H
where Δ V represents the ice volume change, S is the glacier area, and Δ H is the change in the elevation of the ice.

2.4. Uncertainty Measurements

2.4.1. Glacier Area

The uncertainty in the glacier area extraction is composed of a system error and random error [51]. The system error, such as identifying snow and mountain shadow, can be reduced by selecting high quality remote sensing images and glacier expert interpretation, but the system error is difficult to quantitatively estimate. In the process of glacier outlier extraction, the random error is inevitably affected by the image spatial resolution. The uncertainty in the glacier outlier is considered to be equal to the perimeter of the grid number through the glacier boundary. The glacier outlier from raster to vector includes the echelon form, straight line, and tooth shape forms. The majority of the glacier’s boundary is in echelon form, while only a few are in straight line form. The tooth shape boundary is usually at the end of the glacier, and the influence of this type boundary on the number of total glacier grids can be ignored [33]. Therefore, when comparing the total glacier area, the area error σS caused by the remote sensing image resolution is expressed as
σ s = ( 1 / 2 × Perimeter + Perimeter 2 ) × resolution
where Perimeter is the perimeter of the stepped glacier outlier; Perimeter 2 represents the perimeter of the straight line of the glacier outlier; and resolution represents the spatial resolution of the remote sensing image.

2.4.2. Ice Surface Elevation

Generally, the best way to evaluate the elevation error is to use a more precise ground control point. Since ZY-3 has a very high resolution, we could not achieve a higher accuracy of the ground control point, so we used the residual error of the off-glacier region to estimate the error. The error of the sinusoidal model was calculated according to the standard error of each parameter [33], as follows:
σ 1 = k 2 σ 2 + ( 14.74335 ) 2 × s i n 2 ( X ) + σ s i n 2 × c o s 2 ( X )
where k represents the slope of the linear model, σ is the standard deviation of the regression model, and the last two terms on the right give the standard error of the sinusoidal part. The relative errors of the five regions were calculated using Equation (9).
The average error of the elevation difference σH in the basin was calculated as
σ H = i = 1 n σ i 2 / n
where n is the number of images, and σi is the standard error of the sinusoidal curve for the ith image.
The ZY-3 image displayed snow cover and mountain shadows in parts of the image, which can cause high elevation overestimation in areas with thick snow, and underestimation in areas in the mountain’s shadow. Therefore, these anomalies need to be removed before analyzing the changes. In this study, the threshold scheme was adopted [52,53], and the thread difference between DEMs was ±100 m.

2.4.3. Ice Storage Volume Change

Considering that the change in ice storage volume is affected by the differences between the elevation error and glacial area error, the error of ice volume change ε was estimated as
ε = ( S × σ H ) 2 + ( Δ H × σ S ) 2
where σH is the error of elevation change, and σ S the error in the area estimation.
The average density of 850 ± 60 kg·m−3 was used as the volume–mass conversion parameter [54], so that the ice storage change could be converted into water equivalent change.

3. Results

3.1. Glacier Area Change

There were 516 glaciers with areas larger than 0.01 km2, with a total area of 481.5 km2 according to the revised CGI. This is less than the 660 glaciers with a total area of 509.9 km2 in the entire Shule River Basin [19]. The difference is because parts of some glaciers, including Laohugou glacier No. 12, were not above Changma hydrological station, so they were not included. A previous study [19] reported that the total glacier area in the Shule River Basin decreased by about 24.5 km2 per decade, from the 1970s to 2010, according to the CGI (1970s) and SCGI (2006 to 2010). Comparatively, our results suggest that the glacier area decreased by about 21.9 km2 per decade, from the 1970s to 2008, which is consist with Sun’s [19] results.
Between 2000 and 2001, there were 513 glaciers with a total area of 432.6 km2 in the URSRB (Figure 4), and 478 glaciers with a total area of 375.6 ± 2.68 km2 during 2014/2015. This finding reveals that at least 35 glaciers disappeared between 2000 and 2015, with a decrease of 57.5 ± 2.68 km2 in total area (about 11.9% of 1970s), and an average decrease of 3.83 ± 0.18 km2·a−1 (0.79 ± 0.04% a−1). The glacier shrinkage was most pronounced from 2004/2006 to 2008/2009 by 9.21 km2·a−1 (1.90% a−1), which is more than twice the average area decrease from 2000/2001 to 2014/2015. In addition, the glacier area during 2014/2015 was nearly the same with that of 2012/2013.
It is interesting to find that although the 6 glaciers with areas greater than 8 km2 remained unchanged, the median glacier area obvious decreased (Figure 5). The number of medium-size glaciers with areas between 3 and 8 km2 decreased. The glacier areas between 1 and 3 km2 displayed similar patterns with those between 3 to 8 km2. The number of glaciers with area less than 1 km2 decreased from 417 to 393 from 2000/2001 to 2014/2015. These results indicate that the change in large glaciers is mainly displayed as area shrinkage, where 35 glaciers with an area less than 1 km2 disappeared. The average size of the 35 glaciers that disappeared was ~0.24 km2 in 2000/2001. The elevation of these glaciers ranged from 4422 to 5526 m a.s.l., indicating that the small glaciers are affected by multiple factors, excluding elevation.

3.2. Spatial Pattern of Glacier Area Change

According to the spatial patterns of glacier change, the whole URSRB can be divided into (I) Northwest, (II) Middle, and (III) Southeast subregions (Figure 6).
The total glacial area of region (I) was about 45.0 km2 in 2000/2001, accounting for 10.3% of the glacier area in the URSRB, and is characterized by a large number of small glaciers. The average elevation of glaciers in this region was 4828 m a.s.l. The glacier area in subregion (I) had a pronounced decrease by 21.5 ± 1.0 km2 (47.8%) between 2000/2001 and 2014/2015. Five glaciers disappeared, with an average area of 0.05 km2, and elevation ranging from 4422 to 5026 m a.s.l., in 2000/2001.
The total glacier area in region (II) was about 342.42 km2, accounting for 78.7% of the total glacier area in the URSRB. In this region, the glacier size was relatively large, among which, 6 glaciers displayed areas larger than 10 km2. The average elevation of glaciers in this region was 4983 m a.s.l. The glacier shrinkage mainly occurred at the terminus of the glacier in the elevation ranging from 4700 to 5100 m a.s.l., and the average area was reduced by about 5.4% from 2000/2001 to 2014/2015. Twenty-five glaciers disappeared, with an average area of 0.03 km2, and elevation ranging from 4599 to 5256 m a.s.l., in 2000/2001.In region (III), there were 59 glaciers with a total area of about 54.2 km2, which accounts for 11% of the total glacier area. The glaciers in this region were relatively small and concentrated. The average elevation of glaciers in this region was 4807 m a.s.l. The total glacier area decreased by about 8.95 ± 0.42 km2 (16.5%) between 2000 and 2015. Two glaciers disappeared, with an average area of 0.03 km2, and elevation range of 4592 to 4776 m a.s.l., in 2000/2001.

3.3. Glacier Surface Elevation Change

Due to the fact that 5 ZY-3 image pairs were used to calculate the ice surface elevation changes in the study, and there are still some areas which were not covered by the 5 ZY-3 image pairs in the URSRB, we therefore use the 5 subregions (a–e) which are in line with each ZY-3 image pair to investigate the characteristics of glacier surface elevation change in the URSRB.
The glacier surface elevation in the URSRB decreased from 1999 to around 2013, and large differences were found between the a, b, c, d, and e subregions that correspond with each ZY-3 image (Figure 7, Table 3). The average glacier surface elevation decreased by about 4.764 ± 0.506 m in subregion a in central UBSRB (97°36′–97°48′ E, 38°30′–38°45′ N) from 2000 to 2012, which is equal to about 0.397 ± 0.042 m·a−1. The average glacier surface elevation of subregion b (96°36′–96°50′ E, 39°21′–39°27′ N) reduced by 4.466 ± 0.906 m from 2000 to 2013, or 0.344 ± 0.070 m·a−1. The average glacier surface elevation of subregion c (97°10′–97°25′ E, 38°35′–38°47′ N) decreased by 3.620 ± 0.254 m from 2000 to 2013, or about 0.279 ± 0.020 m·a−1. The change in average glacier surface elevation of subregion d (97°48′–98°09′ E, 38°20′–38°35′ N) from 2000 to 2013 is different than the other subregions, as it only decreased by 2.386 ± 0.611 m, a slight decrease by 0.184 ± 0.047 m·a−1. The average glacier surface elevation of subregion e (98°33′–98°51′ E, 38°12′–38°18′ N), from 2000 to 2015, had the most pronounced decline among all subregions in the URSRB, decreasing by 8.113 ± 0.333 m, or 0.541 ± 0.022 m·a−1 (Figure 7e). The average glacier surface elevation in the URSRB from 2000 to 2013 reduced by about 4.98 ± 0.6 m, or 0.383 ± 0.046 m·a−1. The ice volume in the URSRB was reduced by about 2.459 ± 0.307 km3, which is equal to 0.189 ± 0.02 km3·a−1 and to a decrease of 2.089 ± 0.354 × 1012 kg·a−1.
It is pertinent to note that the change rates of glacier areas have not always been in line with those of ice surface elevation (Table 3). For instance, in the ice elevation change in d subregion was only −0.184 m·a−1, where the middle has changed in glacier area by −0.87% a−1. The large variations in ice surface elevation in each subregion reflect the complexity of glacier mass change.

4. Discussion

4.1. Uncertainty

Uncertainty in the delineation of glacier outlines is the result of positional and processing errors [55] in which seasonal snow, cloud, and debris cover complicate the precision of the mapping [56]. There are limited debris-covered glaciers in URSRB. The accuracy of the outlines in this study was assessed by comparisons with a previous study by Sun et al. [19], and the results in CGI and SCGI. The uncertainty due to the resolution of the remote sensing data was also evaluated.
The penetration depth of the SRTM C-band radar beam contributed to some uncertainties when the SRTM DEM was used for geodetic mass-balance calculations. A previous study [57] indicates that the penetration depth decreases as the temperature and water content of the surface snow increases. When the penetration depth was 1.4 m, the ice surface elevation decrease rate was 0.31 m·a−1, which increased to 0.46 m·a−1 when the penetration depth increased to 3.4 m.
Another challenge faced in glaciers studies is data voids in SRTM DEM, especially in the glacier accumulation area. In this study, information was available for elevation change from 4300 to 5500 m a.s.l. Therefore, we compared the unfilled version SRTM DEM and SRTM DEM4.1, found that the area of data voids was too small to affect the ice surface elevation significantly (1.0% of the area above 5000 m a.s.l.), so this issue was neglected.

4.2. Changes in Glacier Area and Ice Surface Elevation

The glacier area in the URSRB decreased by 57.5 km2 (11.91% relate to CGI) from 2000/2001 to 2014/2015, with an average annual shrinkage rate of 0.79% a−1, which is about 1.7 times the shrinkage rate from CGI to SCGI. This result indicates that the glacier shrinkage accelerated after 2000, which is consistent with studies on Tuanjiefeng region near to the URSRB [33], where the glacier area decreased 0.20% a−1 between 2000–2006, and 0.45% a−1 during 2006–2010. Thist also agrees with the shrinkage of Laohugou glacier No. 12, where the area decreased by 0.008 km2·a−1 from 1957 to 1994, 0.043 km2·a−1 from 1994 to 2000, and 0.041 km2·a−1 from 2000 to 2007 [29]. It is also similar to Yanglonghe glaciers No. 1 and No. 5 in the middle Qilian Mountain, where the glacier experienced accelerated retreat since 1999 [30]. A similar pattern was also found in Kangri Karpo region, where the glacier area decreased by 0.41% a−1 from 1980 to 2000, and decreased by 0.52% a−1 from 2000 to 2015 [58]. According to the spatial pattern of glacier change for the entire Qilian Mountain range [19], the glacier shrink rate in the URSRB is relatively lower than other parts of the Qilian Mountains, which indicates that accelerated glacier shrinkage is probably more pronounced in other regions of the Qilian Mountains.
In terms of ice surface elevation changes in the Qilian Mountains, an average decrease of 0.64 m·a−1 was observed for nine glaciers in the Ningchan and Shuiguan River regions in the eastern Qilian Mountains from 1972 to 2010 [41]. Average elevation decreases of 0.4 ± 0.22 m·a−1 and 0.33 ± 0.22 m·a−1 were measured in the ablation area of Yanglonghe glacier No. 1 and No. 5 in the eastern Qilian Mountains from 1956 to 2007 [30]. During 2010–2015, an elevation decrease of 0.9 ± 0.5 m a−1 (water equivalent) was observed for Ningchan glacier No. 1 [32]. These results demonstrate that the mass loss of both a single glacier and multiple glaciers in one region has accelerated in the 21st century.
Xu et al. [33] determined the average ice thinning rate of 0.21 ± 0.04 m·a−1 for the glacier near Tuanjiefeng based on SRTM DEM and topographic mapping during 1966–1999. Similarly, our study suggests that the average ice thinning on the URSRB was about 0.38 m ± 0.05 a−1 during 1999–2013, which also indicates that glacier thinning has accelerated. This trend was also observed in the Kangri Karpo Mountain region in southeastern TP, where glaciers thinned on average by 0.24 ± 0.16 m w.e.a−1 during 1980–2000 and 0.71 ± 0.10 m w.e.a−1 during 2000–2014 [59]. This pattern of accelerated loss matches the global record, as well [60].
The recent mass loss of glaciers in the URSRB is larger than the loss of 0.24 ± 0.22 m·a−1 in central Tianshan Mountain, based on SRTM DEM and TanDEM-X during 1999 to 2012 [11]. In addition, glacier mass loss in URSRB is less than the 0.71 ± 0.10 m·a−1 that occurred in the Kangri Karpo Mountain region in southeastern TP from 2000 to 2014 [58], which is consistent with the observed glacier area change pattern in western China [59].

4.3. Glacier Response to Climate Change

The annual average air temperature has significantly (p = 0.01) increased by 0.465 °C per decade between 1971 and 2013. The average air temperature during the snow and ice melting season (May to September) has significantly (p = 0.01) increased by 0.433 °C per decade. No significant (p = 0.05) trend was found for annual precipitation (Figure 8), suggesting that the recent glacier area loss and ice storage volume losses are mainly controlled by increased air temperature. The air temperature increase is more than twice than the global average of warming (0.12 °C per decade, 1951–2012) [7], which reveals that the URSRB has experienced pronounced warming, like many other mountain areas in the world [61].

5. Conclusions

The following conclusions can be made from the investigation of recent glacier changes in URSRB using multisource remote sensing data.
(1)
There were 478 glaciers with a total area of 375.1 ± 2.68 km2 in URSRB during 2014–2015. Thirty-five glaciers disappeared from 2000 to 2015, with a total glacier area shrinkage of 57.5 ± 2.68 km2 (about 11.9% of the area in the 1970s), which equals an annual decrease of 3.83 ± 0.18 km2·a−1 (0.79 ± 0.04% a−1). From 2004/2006 to 2008/2009, glacier shrinkage was most pronounced in which it decreased by 9.21 km2·a−1 (1.90% a−1). The glacier shrink rate in the URSRB is relatively lower than other parts of the Qilian Mountains, which indicates that glacier shrinkage is probably more pronounced in other regions of the Qilian Mountains.
(2)
The average glacier surface elevation of the URSRB from 1999 to 2013 was reduced by about 4.98 ± 0.6 m, which is equal to a mass loss of 0.383 ± 0.046 m·a−1. The ice storage loss has accelerated since 1999. This pattern of accelerated loss is in agreement with the global records reported. The mass loss of glaciers in the URSRB is larger than that in Tianshan, and less than that in southeastern TP, which is consistent with the observed glacier area change pattern in China.
The pronounced glacier shrinkage and mass loss probably had a significant impact on the annual volume and seasonal distribution of river runoff in the URSRB, increasing the difficulty of water supply management of the URSRB and downstream of Shule River. The potential changes in glacier and the subsequent impact on river runoff in URSRB require further study.

Author Contributions

Conceptualization, X.Z. and S.Z.; Methodology, H.L. and Z.Z.; Software, H.L. and Z.Z.; Validation, Q.W., and X.Z.; Formal Analysis, S.Z.; Investigation, X.Z. and Z.Z.; Resources, S.Z.; Data Curation, H.L.; Writing—Original Draft Preparation, X.Z. and H.L.; Writing—Review & Editing, S.Z.; Visualization, X.Z. and Z.Z.; Supervision, S.Z.; Project Administration, S.Z.; Funding Acquisition, S.Z.

Funding

This work was supported by the National Nature Science Foundation of China (41671056, 41730751, 41471060).

Acknowledgments

Thanks to Andrew Stumpf and the three anonymous referees for their helpful suggestions, which significantly improved the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Immerzeel, W.; van Beek, L.; Bierkens, M. Climate change will affect the Asian water towers. Science 2010, 328, 1382–1385. [Google Scholar] [CrossRef] [PubMed]
  2. Immerzeel, W.; Bierkens, M. Asia’s water balance. Nat. Geosci. 2012, 5, 841–842. [Google Scholar] [CrossRef]
  3. Huss, M.; Hock, R. Global-scale hydrological response to future glacier mass loss. Nat. Clim. Chang. 2018, 8, 135–140. [Google Scholar] [CrossRef] [Green Version]
  4. Pritchard, H.D. Asia’s glaciers are a regionally important buffer against drought. Nature 2017, 545, 169–174. [Google Scholar] [CrossRef] [PubMed]
  5. Yao, T.; Thompson, L.; Yang, W.; Yu, W.; Gao, Y.; Guo, X.; Yang, X.; Duan, K.; Zhao, H.; Xu, B.; et al. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings. Nat. Clim. Chang. 2012, 2, 663–667. [Google Scholar] [CrossRef]
  6. Qin, D.; Ding, Y.; Xiao, C.; Kang, S.; Ren, J.; Yang, J.; Zhang, S. Cryospheric Science: Research framework and disciplinary system. Natl. Sci. Rev. 2018, 5, 255–268. [Google Scholar] [CrossRef]
  7. IPCC (Intergovernmental Panel on Climate Change). Climate Change 2013: The Physical Science Basis; Cambridge University Press: Cambridge, MA, USA, 2013. [Google Scholar]
  8. Aizen, V.; Kuzmichenok, V.; Surazakov, A.; Aizen, E. Glacier changes in the central and northern Tien Shan during the last 140 years based on surface and remote-sensing data. Ann. Glaciol. 2006, 43, 202–213. [Google Scholar] [CrossRef]
  9. Bolch, T. Climate change and glacier retreat in northern Tien Shan (Kazakhstan/Kyrgyzstan) using remote sensing data. Glob. Planet. Chang. 2007, 56, 1–12. [Google Scholar] [CrossRef]
  10. Liu, Q.; Liu, S. Response of glacier mass balance to climate change in the Tianshan Mountains during the second half of the twentieth century. Clim. Dyn. 2016, 46, 303–316. [Google Scholar] [CrossRef]
  11. Li, J.; Li, Z.; Zhu, J.; Li, X.; Xu, B.; Wang, Q.; Huang, C.; Hu, J. Early 21st century glacier thickness changes in the Central Tien Shan. Remote Sens. Environ. 2017, 192, 12–29. [Google Scholar] [CrossRef]
  12. Wei, J.; Liu, S.; Guo, W.; Yao, X.; Xu, J.; Bao, W.; Jiang, Z. Surface-area changes of glaciers in the Tibetan Plateau interior area since the 1970s using recent Landsat images and historical maps. Ann. Glaciol. 2014, 55, 213–222. [Google Scholar] [Green Version]
  13. Pieczonka, T.; Bolch, T.; Wei, J.; Liu, S. Heterogeneous mass loss of glaciers in the Aksu-Tarim Catchment (Central Tien Shan) revealed by 1976 KH-9 Hexagon and 2009 SPOT-5 stereo imagery. Remote Sens. Environ. 2013, 130, 233–244. [Google Scholar] [CrossRef] [Green Version]
  14. Bolch, T. The state and fate of Himalayan Glaciers. Science 2012, 336, 310–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Shangguan, D.; Liu, S.; Ding, Y.; Wu, L.; Deng, W.; Guo, W.; Wang, Y.; Xu, J.; Yao, X.; Guo, Z. Glacier changes in the Koshi River basin, central Himalaya, from 1976 to 2009, derived from remote-sensing imagery. Ann. Glaciol. 2014, 55, 61–68. [Google Scholar] [CrossRef]
  16. Wang, P.; Li, Z.; Yu, G.; Li, H.; Wang, W.; Huai, B.; Zhou, P.; Jin, S.; Wang, L.; Zhang, H. Glacier shrinkage in the Daxue and Danghenan ranges of the western Qilian Mountains, China, from 1957 to 2010. Environ. Earth Sci. 2016, 75, 1–11. [Google Scholar] [CrossRef]
  17. Kääb, A.; Treichler, D.; Nuth, C.; Berthier, E. Brief Communication: Contending estimates of 2003–2008 glacier mass balance over the Pamir–Karakoram–Himalaya. Cryosphere 2015, 9, 557–564. [Google Scholar] [CrossRef]
  18. Zhang, Z.; Liu, S.; Wei, J.; Xu, J.; Guo, W.; Bao, W.; Jiang, Z. Mass Change of Glaciers in Muztag Ata-Kongur Tagh, Eastern Pamir, China from 1971/76 to 2013/14 as Derived from Remote Sensing Data. PLoS ONE 2016, 11, e0147327. [Google Scholar] [CrossRef] [PubMed]
  19. Sun, M.; Liu, S.; Yao, X.; Guo, W.; Xu, J. Glacier changes in the Qilian Mountains in the past half century: Based on the revised First and Second Chinese Glacier Inventory. Acta Geogr. Sin. 2018, 28, 206–220. [Google Scholar] [CrossRef]
  20. Gao, X.; Ye, B.; Zhang, S.; Qiao, C.; Zhang, X. Glacier runoff variation and its influence on river runoff during 1961–2006 in the Tarim River Basin, China. Sci. China Earth Sci. 2010, 53, 880–891. [Google Scholar] [CrossRef]
  21. Thayyen, R.J.; Gergan, J.T. Role of glaciers in watershed hydrology: A preliminary study of a “Himalayan catchment”. Cryosphere 2010, 4, 115–128. [Google Scholar] [CrossRef]
  22. Ye, Q.; Bolch, T.; Naruse, R.; Wang, Y.; Zong, J.; Wang, Z.; Zhao, R.; Yang, D.; Kang, S. Glacier mass changes in Rongbuk catchment on Mt. Qomolangma from 1974 to 2006 based on topographic maps and ALOS PRISM data. J. Hydrol. 2015, 530, 273–280. [Google Scholar] [CrossRef]
  23. Paul, F.; Haeberli, W. Spatial variability of glacier elevation changes in the Swiss Alps obtained from two digital elevation models. Geophys. Res. Lett. 2008, 35, L21502. [Google Scholar] [CrossRef]
  24. Bao, W.; Liu, S.; Wei, J.; Guo, W. Glacier changes during the past 40 years in the West Kunlun Shan. J. Mt. Sci. 2015, 12, 344–357. [Google Scholar] [CrossRef]
  25. Shangguan, D.; Bolch, T.; Ding, Y.; Kröhnert, M.; Pieczonka, T.; Wetzel, H.; Liu, S. Mass changes of Southern and Northern Inylchek Glacier, Central Tian Shan, Kyrgyzstan, during ∼1975 and 2007 derived from remote sensing data. Cryosphere 2015, 9, 703–717. [Google Scholar] [CrossRef]
  26. Wei, J.; Liu, S.; Guo, W.; Xu, J.; Bao, W.; Shangguan, D. Changes in glacier volume in the north bank of the Bangong Co Basin from 1968 to 2007 based on historical topographic maps, SRTM, and ASTER stereo images. Arct. Antarct. Alpine Res. 2015, 47, 301–311. [Google Scholar]
  27. Wei, J.; Liu, S.; Xu, J.; Guo, W.; Bao, W.; Shangguan, D.; Jiang, Z. Mass loss from glaciers in the Chinese Altai Mountains between 1959 and 2008 revealed based on historical maps, SRTM, and ASTER images. J. Mt. Sci. 2015, 12, 330–343. [Google Scholar] [CrossRef]
  28. Du, W.; Qin, X.; Liu, Y.; Wang, X. Variation of the Laohugou glacier No. 12 in the Qilian mountains (in Chinese with English Abstract). J. Glaciol. Geocryol. 2008, 30, 373–379. [Google Scholar]
  29. Zhang, Y.; Liu, S.; Shangguan, D.; Li, J.; Zhao, J. Thinning and shrinkage of Laohugou No. 12 glacier in the Western Qilian Mountains, China, from 1957 to 2007. J. Mt. Sci. 2012, 9, 343–350. [Google Scholar] [CrossRef]
  30. Shangguan, D.; Liu, S.; Ding, Y.; Zhang, Y.; Li, J.; Li, X.; Wu, Z. Changes in the elevation and extent of two glaciers along the Yanglonghe River, Qilian Shan, China. J. Glaciol. 2010, 56, 309–317. [Google Scholar] [CrossRef]
  31. Cao, B.; Pan, B.; Wang, J.; Shangguan, D.; Wen, Z.; Qi, W.; Cui, H.; Lu, Y. Changes in the glacier extent and surface elevation along the Ningchan and Shuiguan river source, eastern Qilian Mountains, China. Quat. Res. 2014, 81, 531–537. [Google Scholar] [CrossRef]
  32. Cao, B.; Pan, B.; Guan, W.; Wang, J.; Wen, Z. Changes in ice volume of the Ningchan No. 1 Glacier, China, from 1972 to 2014, as derived from in situ measurements. J. Glaciol. 2017, 63, 1025–1033. [Google Scholar] [CrossRef]
  33. Xu, J.; Liu, S.; Zhang, S.; Guo, W.; Wang, J. Recent Changes in Glacial Area and Volume on Tuanjiefeng Peak region of Qilian Mountains, China. PLoS ONE 2013, 8, e70574. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, H.; Ye, B.; Wu, J.; Li, M.; Qin, J.; Wang, X.; Wang, J. Observation and study on the CO2 flux in an alpine meadow ecosystem in the upper reach of Shule River Basin (in Chinese with English Abstract). Acta Prataculturae Sin. 2013, 22, 18–26. [Google Scholar]
  35. Rodriguez, E.; Morris, C.S.; Belz, J.E. A global assessment of the SRTM performance. Photogram. Eng. Remote Sens. 2006, 72, 249–260. [Google Scholar] [CrossRef]
  36. Rabus, B.; Eineder, M.; Roth, A.; Bamler, R. The shuttle radar topography mission—A new class of digital elevation models acquired by spaceborne radar. ISPRS J. Photogram. Remote Sens. 2003, 57, 241–262. [Google Scholar] [CrossRef]
  37. Li, D. China’s first civilian three-line-array stereo mapping satellite: ZY-3 (in Chinese with English Abstract). Acta Geod. Cartograph. Sin 2012, 41, 317–322. [Google Scholar]
  38. Zhao, Q.; Zhang, S.; Ding, Y.; Wang, J.; Han, H.; Xu, J.; Zhao, C.; Guo, W.; Shangguan, D. Modeling hydrologic response to climate change and shrinking glaciers in the highly Glacierized Kunma Like River catchment, Central Tian Shan. J. Hydrometeorol. 2015, 16, 2383–2402. [Google Scholar] [CrossRef]
  39. Zhao, Q.; Chen, R.; Zhang, S.; Qin, J. The comprehensive assessment of river runoff changes in cold regions catchments, China. In Impact of Cryospheric Changes on the River Runoff in the Cold Regions of Western China; Chen, R., Zhang, S., Yang, Y., Liu, J., Zhao, Q., Eds.; Science Press: Beijing, China, 2018. [Google Scholar]
  40. Xu, J.; Zhang, S.; Shangguan, D. Glacier change in headwaters of the Yangtze River in recent three decades (in Chinese with English Abstract). Arid Zone Res. 2013, 30, 919–926. [Google Scholar]
  41. Cao, B.; Pan, B.; Gao, H.; Jiang, S.; Wen, Y.; Shangguan, D. Glacier variation in the Lenglongling Range of eastern Qilian Mountains from 1972 to 2007 (in Chinese with English Abstract). J. Glaciol. Geocryol. 2010, 32, 242–248. [Google Scholar]
  42. Bolch, T.; Yao, T.; Kang, S.; Buchroithner, M.F.; Scherer, D.; Maussion, F.; Huintjes, E.; Schneider, C. A glacier inventory for the western Nyainqentanglha Range and Nam Co Basin, Tibet, and glacier changes 1976–2009. Cryosphere 2010, 4, 419–433. [Google Scholar] [CrossRef]
  43. Guo, W.; Liu, S.; Xu, J.; Wu, L.; Shangguan, D.; Yao, X.; Wei, J.; Bao, W.; Yu, P.; Liu, Q. The second Chinese glacier inventory: Data, methods and results. J. Glaciol. 2015, 61, 357–372. [Google Scholar] [CrossRef]
  44. Berthier, E.; Toutin, T. SPOT5-HRS digital elevation models and the monitoring of glacier elevation changes in North-West Canada and South-East Alaska. Remote Sens. Environ. 2008, 112, 2443–2454. [Google Scholar] [CrossRef]
  45. Zhang, G.; Jiang, Y.; Li, D.; Huang, W.; Pan, H.; Tang, X.; Zhu, X. In Orbit Geometric Calibration and Validation of Zy-3 Linear Array Sensors. Photogram. Rec. 2014, 29, 68–88. [Google Scholar] [CrossRef]
  46. Peduzzi, P.; Herold, C.; Silverio, W. Assessing high altitude glacier thickness, volume and area changes using field, GIS and remote sensing techniques: The case of Nevado Coropuna(Peru). Cryosphere 2010, 4, 313–323. [Google Scholar] [CrossRef]
  47. Nuth, C.; Kääb, A. Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change. Cryosphere 2011, 5, 271–290. [Google Scholar] [CrossRef] [Green Version]
  48. Berthier, E.; Arnaud, Y.; Vincent, C.; Re´my, F. Biases of SRTM in high-mountain areas. Implications for the monitoring of glacier volume changes. Geophys. Res. Lett. 2006, 33, L08502. [Google Scholar] [CrossRef]
  49. Gardelle, J.; Berthier, E.; Arnaud, Y. Impact of resolution and radar penetration on glacier elevation changes computed from DEM differencing. J. Glaciol. 2012, 58, 419–422. [Google Scholar] [CrossRef] [Green Version]
  50. Kaab, A.; Berthier, E.; Nuth, C.; Gardelle, J.; Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 2012, 488, 495–498. [Google Scholar] [CrossRef] [PubMed]
  51. Frenierre, J.L.; Mark, B.G. A review of methods for estimating the contribution of glacial meltwater to total watershed discharge. Prog. Phys. Geogr. 2014, 38, 173–200. [Google Scholar] [CrossRef]
  52. Berthier, E.; Schiefer, E.; Clarke, G.K.C.; Menounos, B.; Remy, F. Contribution of Alaskan glaciers to sea-level rise derived from satellite imagery. Nat. Geosci. 2010, 3, 92–95. [Google Scholar] [CrossRef] [Green Version]
  53. Bolch, T.; Pieczonka, T.; Benn, D. Multi-decadal mass loss of glaciers in the Everest area (Nepal Himalaya) derived from stereo imagery. Cryosphere 2011, 5, 349–358. [Google Scholar] [CrossRef] [Green Version]
  54. Huss, M. Density assumptions for converting geodetic glacier volume change to mass change. Cryosphere 2013, 7, 877–887. [Google Scholar] [CrossRef] [Green Version]
  55. Bolch, T.; Menounos, B.; Wheate, R. Landsat-based inventory of glaciers in western Canada, 1985–2005. Remote Sens. Environ. 2010, 114, 127–137. [Google Scholar] [CrossRef]
  56. Paul, F.; Barrand, N.E.; Baumann, S.; Berthier, E.; Bolch, T.; Casey, K.; Frey, H.; Joshi, S.; Konovalov, V.; Le Bris, R. On the accuracy of glacier outlines derived from remote-sensing data. Ann. Glaciol. 2013, 54, 171–182. [Google Scholar] [CrossRef] [Green Version]
  57. Surdyk, S. Using microwave brightness temperature to detect short-term surface air temperature changes in antarctica: An analytical approach. Remote Sens. Environ. 2002, 80, 256–271. [Google Scholar]
  58. Wu, K.; Liu, S.; Jiang, Z.; Xu, J.; Wei, J.; Guo, W. Recent glacier mass balance and area changes in the Kangri Karpo Mountains from DEMs and glacier inventories. Cryosphere 2018, 12, 103–121. [Google Scholar] [CrossRef] [Green Version]
  59. Ding, Y.; Liu, S.; Li, J.; Shangguan, D. The retreat of glaciers in response to recent climate warming in western China. Ann. Glaciol. 2006, 43, 97–105. [Google Scholar] [CrossRef] [Green Version]
  60. Zemp, M.; Frey, H.; Gärtner-Roer, I.; Nussbaumer, S.U.; Hoelzle, M.; Paul, F.; Haeberli, W.; Denzinger, F.; AhlstrØm, A.P.; Anderson, B.; et al. Historically unprecedented global glacier decline in the early 21st century. J. Glaciol. 2015, 61, 745–763. [Google Scholar] [CrossRef]
  61. Mountain Research Initiative EDW Working Group. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 2015, 5, 424–430. [Google Scholar] [Green Version]
Figure 1. Location of Upper Reach of Shule River Basin, the blue box in left panel is Tuanjiefeng region.
Figure 1. Location of Upper Reach of Shule River Basin, the blue box in left panel is Tuanjiefeng region.
Water 10 00796 g001
Figure 2. The sinusoidal relationship between the elevation difference between ZY-3 digital elevation model (DEM) and Shuttle Radar Topography Mission (SRTM) DEM and aspects in (a) sub-region a; (b) sub-region b; (c) sub-region c; (d) sub-region d; (e) sub-region e of Upper Reach of Shule River Basin
Figure 2. The sinusoidal relationship between the elevation difference between ZY-3 digital elevation model (DEM) and Shuttle Radar Topography Mission (SRTM) DEM and aspects in (a) sub-region a; (b) sub-region b; (c) sub-region c; (d) sub-region d; (e) sub-region e of Upper Reach of Shule River Basin
Water 10 00796 g002
Figure 3. Histogram of the elevation difference between ZY-3 DEM and SRTM DEM before co-registration correction of (a) sub-region a; (b) sub-region b; (c) sub-region c; (d) sub-region d; (e) sub-region e; and after co-registration correction of (f) sub-region a; (g) sub-region b; (h) sub-region c; (i) sub-region d; (j) sub-region e of Upper Reach of Shule River Basin
Figure 3. Histogram of the elevation difference between ZY-3 DEM and SRTM DEM before co-registration correction of (a) sub-region a; (b) sub-region b; (c) sub-region c; (d) sub-region d; (e) sub-region e; and after co-registration correction of (f) sub-region a; (g) sub-region b; (h) sub-region c; (i) sub-region d; (j) sub-region e of Upper Reach of Shule River Basin
Water 10 00796 g003
Figure 4. Variations in the glacier area in the Upper Reach of Shule River Basin from Chinese Glacier Inventory (CGI) (1970), 2000 to 2015.
Figure 4. Variations in the glacier area in the Upper Reach of Shule River Basin from Chinese Glacier Inventory (CGI) (1970), 2000 to 2015.
Water 10 00796 g004
Figure 5. Box-whisker plot of different glacier area classes in Upper Reach of Shule River Basin; the red total is the number of glaciers in this class; the periods in the columns from left to right are 2000/2001, 2004/2005, 2008/2009, 2012/2013, and 2014/2015, respectively. The top and bottom of the whisker represent the maximum and minimum glacier area in each class, respectively; the top, median, and bottom of the box represent the 75%, 50%, and 25% quantiles of the glacier area in the class, respectively.
Figure 5. Box-whisker plot of different glacier area classes in Upper Reach of Shule River Basin; the red total is the number of glaciers in this class; the periods in the columns from left to right are 2000/2001, 2004/2005, 2008/2009, 2012/2013, and 2014/2015, respectively. The top and bottom of the whisker represent the maximum and minimum glacier area in each class, respectively; the top, median, and bottom of the box represent the 75%, 50%, and 25% quantiles of the glacier area in the class, respectively.
Water 10 00796 g005
Figure 6. The glacier area changes in different subregions in the Upper Reach of Shule River Basin.
Figure 6. The glacier area changes in different subregions in the Upper Reach of Shule River Basin.
Water 10 00796 g006
Figure 7. Glacier ice surface elevation changes from 2000 to 2015 in sub-region a, b, c, d, and e in the Upper Reach of Shule River Basin.
Figure 7. Glacier ice surface elevation changes from 2000 to 2015 in sub-region a, b, c, d, and e in the Upper Reach of Shule River Basin.
Water 10 00796 g007
Figure 8. Annual average air temperature (°C), average air temperature in melt season (°C), and annual precipitation (mm) in the URSRB during 1971 to 2013, where k is the slope of the linear trend, and p is the significance level of the t test.
Figure 8. Annual average air temperature (°C), average air temperature in melt season (°C), and annual precipitation (mm) in the URSRB during 1971 to 2013, where k is the slope of the linear trend, and p is the significance level of the t test.
Water 10 00796 g008
Table 1. List of remote sensing images of multiperiods for deriving glacier area changes in the Upper Reach of Shule River Basin.
Table 1. List of remote sensing images of multiperiods for deriving glacier area changes in the Upper Reach of Shule River Basin.
YearSatellitePath/RowAcquisition Date (Year-Month-Day)
2000/2001Landsat 7134/332000-10-04
135/332000-05-20
136/322000-10-02
136/332000-07-14
136/332001-09-03
135/332001-06-24
134/332001-07-03
2004/2006Landsat 5134/332004-08-04
135/332004-08-27
136/322004-10-21
136/332004-08-02
134/332005-08-23
135/332005-05-10
136/322005-07-20
136/332005-07-20
134/332006-09-11
135/332006-09-18
136/322006-09-25
136/332006-10-11
2008/2009Landsat 5134/332008-08-31
135/332008-08-22
136/332008-08-29
134/332009-07-17
135/332009-08-09
136/322009-09-17
136/332009-09-01
2012/2013Landsat 7135/332012-08-25
136/322012-09-01
136/332012-09-01
Landsat 8134/332013-09-30
135/332013-10-07
136/322013-10-14
136/332013-10-14
2014/2015Landsat 8134/332014-08-16
135/332014-09-08
136/322014-08-30
136/332014-08-30
134/332015-09-04
135/332015-08-10
136/322015-08-17
136/332015-08-17
Table 2. List of remote sensing images for deriving glacier elevation changes in the Upper Reach of Shule River Basin.
Table 2. List of remote sensing images for deriving glacier elevation changes in the Upper Reach of Shule River Basin.
SatellitePath/RowAcquisition Date (Year-Month-Day)Remark
ZY-347/1282012-08-19a
50/1262013-08-03b
49/1282013-08-08c
46/1282013-10-26d
44/1292015-02-15e
Note: a, b, c, d, and e in the Remark column correspond to each correction equation in the Equations (3)–(7), which correspond with the elevation changes in the glacier sub-regions.
Table 3. Changes in glacier area and ice surface elevation of each subregion of the Upper Reach of Shule River Basin.
Table 3. Changes in glacier area and ice surface elevation of each subregion of the Upper Reach of Shule River Basin.
SubregionPeriodAverage Elevation (m)Area Change (km2)Area Change (% a−1)Elevation Change (m)Elevation Change (m·a−1)
a2000–20124989−6.27 −0.54 −4.764 ± 0.506−0.397 ± 0.042
b2000–20134859−12.44 −3.07 −4.466 ± 0.906 −0.344 ± 0.070
c2000–20135044−4.69 −0.52 −3.620 ± 0.254−0.279 ± 0.020
d2000–20134939−8.03 −0.87 −2.386 ± 0.611−0.184 ± 0.047
e2000–20154807−6.60 −0.90 −8.113 ± 0.333−0.541 ± 0.022

Share and Cite

MDPI and ACS Style

Zhang, X.; Li, H.; Zhang, Z.; Wu, Q.; Zhang, S. Recent Glacier Mass Balance and Area Changes from DEMs and Landsat Images in Upper Reach of Shule River Basin, Northeastern Edge of Tibetan Plateau during 2000 to 2015. Water 2018, 10, 796. https://doi.org/10.3390/w10060796

AMA Style

Zhang X, Li H, Zhang Z, Wu Q, Zhang S. Recent Glacier Mass Balance and Area Changes from DEMs and Landsat Images in Upper Reach of Shule River Basin, Northeastern Edge of Tibetan Plateau during 2000 to 2015. Water. 2018; 10(6):796. https://doi.org/10.3390/w10060796

Chicago/Turabian Style

Zhang, Xiaowen, Haojie Li, Zhihua Zhang, Qianxin Wu, and Shiqiang Zhang. 2018. "Recent Glacier Mass Balance and Area Changes from DEMs and Landsat Images in Upper Reach of Shule River Basin, Northeastern Edge of Tibetan Plateau during 2000 to 2015" Water 10, no. 6: 796. https://doi.org/10.3390/w10060796

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop