Next Article in Journal
An Elevation Ambiguity Resolution Method Based on Segmentation and Reorganization of TomoSAR Point Cloud in 3D Mountain Reconstruction
Next Article in Special Issue
Decadal Changes in Glacier Area, Surface Elevation and Mass Balance for 2000–2020 in the Eastern Tanggula Mountains Using Optical Images and TanDEM-X Radar Data
Previous Article in Journal
Integrating Satellite, UAV, and Ground-Based Remote Sensing in Archaeology: An Exploration of Pre-Modern Land Use in Northeastern Iraq
Previous Article in Special Issue
ACFNet: A Feature Fusion Network for Glacial Lake Extraction Based on Optical and Synthetic Aperture Radar Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glacier Area and Snow Cover Changes in the Range System Surrounding Tarim from 2000 to 2020 Using Google Earth Engine

1
State Key Laboratory of Remote Sensing Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Faculty of Civil Engineering and Earth Sciences, Delft University of Technology, 2628 CN Delft, The Netherlands
4
Key Laboratory for Geographical Process Analysis & Simulation of Hubei Province, College of Urban and Environmental Sciences, Central China Normal University, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(24), 5117; https://doi.org/10.3390/rs13245117
Submission received: 11 October 2021 / Revised: 8 December 2021 / Accepted: 13 December 2021 / Published: 16 December 2021
(This article belongs to the Special Issue Remote Sensing in Glaciology and Cryosphere Research)

Abstract

:
Glacier and snow are sensitive indicators of regional climate variability. In the early 21st century, glaciers in the West Kunlun and Pamir regions showed stable or even slightly positive mass budgets, and this is anomalous in a worldwide context of glacier recession. We studied the evolution of snow cover to understand whether it could explain the evolution of glacier area. In this study, we used the thresholding of the NDSI (Normalized Difference Snow Index) retrieved with MODIS data to extract annual glacier area and snow cover. We evaluated how the glacier trends related to snow cover area in five subregions in the Tarim Basin. The uncertainty in our retrievals was assessed by comparing MODIS results with the Landsat-5 TM in 2000 and Landsat-8 OLI in 2020 glacier delineation in five subregions. The glacier area in the Tarim Basin decreased by 1.32%/a during 2000–2020. The fastest reductions were in the East Tien Shan region, while the slowest relative reduction rate was observed in the West Tien Shan and Pamir, i.e., 0.69%/a and 1.08%/a, respectively, during 2000–2020. The relative glacier stability in Pamir may be related to the westerlies weather system, which dominates climate in this region. We studied the temporal variability of snow cover on different temporal scales. The analysis of the monthly snow cover showed that permanent snow can be reliably delineated in the months from July to September. During the summer months, the sequence of multiple snowfall and snowmelt events leads to intermittent snow cover, which was the key feature applied to discriminate snow and glacier.

1. Introduction

The Tarim Basin is the largest inland drainage basin in China, located in the north-west of the Tibetan Plateau and far away from the ocean. The basin is surrounded by the mountain ranges of the Tien Shan in the north, Pamir in the west, West Kunlun in the south and East Kunlun in the southeast [1]. These high mountains are the ideal place for glacial development. Glaciers and snow meltwater account for about 75% of the runoff [2], thus being the main contributor to water security in the Tarim Basin, particularly for irrigated agriculture along the fringes of the Taklamakan Desert [3,4,5]. Therefore, changes in upstream glaciers have an effect on the whole Tarim Basin runoff [6,7]. In the context of global warming, the annual average temperature has increased at a rate of 0.2 °C per decade in the Tarim Basin [8]. Extreme hydrological events have intensified, such as glacier lake outburst floods (GLOFs) and glacier surging [9,10,11]. Thus, it is essential to investigate glacier variability and the potential impacts on the long-term sustainability of water supply.
Snow is an essential component of the water cycle and accounts for a large part of total precipitation in high elevation mountain areas [12]. It plays a key role in nourishing glaciers and in reducing the absorption of solar irradiance by the underlaying glaciers, because of its very high reflectance. Mapping snow cover from satellite images provides useful information on mass accumulation in glaciers [13,14]. Snow cover varies significantly both spatially and temporally over the Tarim Basin, but the details of the temporal variability in snow cover have received less attention [12].
Glaciers in the Tarim Basin show a very diverse response to the variability of snow cover in different mountain ranges. The remote position and dispersion of glaciers restricted glacier monitoring in the Tarim Basin, and most of the previous glacier studies focused on several single glaciers [15,16]. Thus, glacier variability across the entire Tarim Basin during the past two decades is still unclear, as, for example, shown by two studies relating to largely overlapping periods. Shangguan et al. [17] and Liu et al. [18] investigated changes of about 7665 glaciers in the Tarim Basin by comparing the Chinese glacier inventory compiled between 1960 and 1970 with Landsat ETM+ images in 1999–2001, respectively, which indicated that the total glacier area decreased by 3.3% (456.2 km2) in this region. Zhang et al. [19] concluded that the glaciers in the Aksu River sub-basin (West Tien Shan Region) declined 965.7 km2 from 1975 to 2016. Most glaciers are retreating, as expected, but, conversely, West Kunlun (W Kunlun) glaciers had stable or even slightly positive mass budgets in the early 21st century, and this is anomalous in a global context of glacier reduction [16]. In addition, the Pamir Mountains are also characterized by relatively stable glaciers [20]. The westerly winds feed most of the total annual precipitation in the Pamir Mountains [21]. Previous studies suggested multiple causes of the Pamir anomaly, including summertime cooling [22,23,24], increasing winter snowfall [25,26,27] and summertime precipitation [28]. W Kunlun glaciers are more sensitive to changes in winter precipitation, rather than to temperature variability, contrary to other regions of the Tarim Basin [29]. From ICES at laser altimetry measurements, Kääb et al. [30] found that the W Kunlun glacier stability has persisted for at least 40 years.
The increasing annual precipitation may result from the strengthening westerlies [31]. Therefore, the increased winter precipitation probably led to snow accumulation, which might have offset summer warming and ice ablation [16]. It has also been suggested that intensive regional irrigation led to increased cloudiness during summer and reduced irradiance and ice ablation [28]. These studies have tried to analyze the Pamir and W Kunlun glaciers anomaly from a specific perspective, identifying possible factors that could bring about glacier growth. A more systematic analysis on the relationship between variability in snow cover and glacier extent in the Tarim Basin is lacking.
Both snow cover and glaciers can be monitored using multispectral satellite image data [32], which offers a practical approach to monitor glacier and snow variability in mountain ranges. Various band ratios have been applied to map snow cover and glacier extents, such as the near infrared (NIR)/shortwave infrared (SWIR) [33] and the red/SWIR [34] ratios or the Normalized Difference Snow Index (NDSI, see equation in Section 3.2), in combination with a threshold [35]. It is difficult, however, to distinguish snow and glaciers with these method under all conditions, as, e.g., glacier ice and compressed snow have similar spectral properties. The NDSI threshold algorithm is able to effectively distinguish snow and ice from other surface features [36,37], including the discrimination of snow and glaciers from clouds and shadows [38]. Mapping snow cover has a long history, starting in 1989 with early tests using Landsat multispectral scanner (MSS) satellite images [35], to the recent use of daily MODIS data to derive time series of snow cover [39]. NDSI is a widely used method to delineate glaciers and snow cover [32,40]; therefore, we applied thresholding NDSI to identify glaciers and snow in this study. The debris-covered glaciers are very difficult to recognize from satellite images [32,41], while the debris-covered glacier in High Mountain Asia was estimated to be about 10% of the total glacier area based on the manually delineated glacier boundary from images [42]. Most of the debris-covered portion is located at the glacier terminus, which is only a limited portion of glacier extent [43]. Taking into account that debris-covered ice cannot be detected by NDSI, we only distinguished snow and clear ice in this study. The usability of Moderate-Resolution Imaging Spectroradiometer (MODIS) images is improved by the sub-daily temporal resolution and global coverage [44]. The Terra and Aqua satellites orbit the earth and cross the equator about 3 h apart, increasing the likelihood of acquiring cloud-free images [45]. In order to diminish data volumes and to eliminate cloud-affected observations, composite products are produced from observations collected over the course of 8 days [46]. At present, the MODIS snow cover products are generated and made available by the National Aeronautics and Space Administration (NASA). The snow products generated with the MODIS data acquired from Terra (MOD10A2) and Aqua (MYD10A2) are an eight-day combined dataset available at a spatial resolution of 500 m, which contains maximum snow cover information over the compositing period [47].
Snow and glacier ice are characterized by high reflectance in the visible bands and low reflectance in the near-infrared [48], which explains why thresholding NDSI is widely used for snow cover mapping [36,49]. Mishra et al. [50] evaluated the range of NDSI values of snow in the Himalayan region and found that NDSI can vary between 0.04 and 0.92 as the amount of snow cover increases. Zhang et al. [51] suggested that the NDSI threshold of 0.35 is the optimal threshold for use in the Tibetan Plateau based on daily snow-depth measurements during 2003–2013. For global analyses, an NDSI threshold of 0.4 is still recommended [44]. In this study, we also use the standard NDSI threshold of 0.4.
Google Earth Engine (GEE) is a platform that provides a large collection of satellite data and products and allows users to visualize and analyze data through a web-based system [52]. The platform is an application programming interface, which uses a combination of JavaScripts and functions to process and analyze the data that are provided to the user. GEE is a powerful geospatial processing and analysis tool not only because it is computationally time-efficient but also due to the expanding availability of many open-source geospatial datasets. It has been widely used to map rock glaciers [53], water areas [54], wildfire detection [54], among other applications. To explore the snow cover relationship with changes in glacier area, we carried out this study using Google Earth Engine.
We used MODIS reflectance products and spatial analysis techniques to investigate the evolution of glacier area from 2000 to 2020 in mountain ranges in the Tarim Basin. According to the drainage basin delineation from the Hydrological data and maps based on Shuttle Elevation Derivatives at multiple Scales (HydroSHEDS), we divided the whole Tarim Basin into five regions: East Tien Shan, West Tien Shan, Pamir, West Kunlun, and East Kunlun. We examined the evolution of snow cover in relation to its effects on the glacier area in different regions in the Tarim Basin, and further explored whether the Pamir-W Kunlun glacier anomaly was related to the regional evolution of snow cover. We evaluated the MODIS glacier area using higher spatial resolution glacier area retrieved from Landsat 5—Thematic Mapper (TM) in 2000 and L8—Operational Land Imager (OLI) in 2020. Finally, we combined our results with digital elevation models and previous studies of glaciers to evaluate local glacial changes and their effects on water resources.

2. Study Area

The range system surrounding Tarim Basin (34–44°N, 75–92°E) covers an area of about 1,020,000 km2, with elevations ranging between 773 and 7500 m. The mountain streams originate from snow and glacier meltwater at high elevations. The basin is mainly fed by four major streams: Aksu River from West Tien Shan to the north, the Karakax River and Yurungkax originating in the north slope of Kunlun Shan to the south, and the Yarkant River from the Pamir Mountains [55] (Figure 1). The Taklamakan Desert is situated in the middle of the Tarim Basin, encircled by mountains and glaciers. It is characterized by an extremely arid desert climate, with little precipitation and high potential evaporation [56]. Due to the extremely arid desert climate in the lower reach of these rivers, most of the rivers disappear in the Taklamakan Desert.
The climate of the Tarim Basin is continental warm-temperate, with a mean annual precipitation of approximately 50 mm and potential evapotranspiration of up to 2200 mm annually [57]. Under the dominance of westerlies, the precipitation at high elevation is above 300 mm and almost always in the form of snow. Temperature ranges between −3 °C in the winter and 40 °C in summer [58].

3. Data and Method

3.1. Data

3.1.1. MODIS 8-Day Surface Reflectance Product (MOD09A1)

MODIS 8-day 500 m Surface Reflectance product (MOD09A1) was used in this study to extract glacier and seasonal snow cover. To map annual glacier area, six tiles of MODIS, h23v4, h23v5, h24v4, h24v5, h25v4, and h25v5, were employed to construct mosaic images. The MODIS product version V006 is based on improved snow, cloud and cloud shadow detection algorithms [59]. These MODIS products were imported in the Google Earth Engine Application Programming Interface (GEE API). Data acquired between July and September, when snow cover was minimum, were used to delineate glacier area. Temporal composite image was created by applying the GEE reducer function to obtain the median radiance over all valid samples from July to September for each pixel. This composite image was then used to calculate the NDSI and then to extract the annual glacier outline by applying the prescribed threshold. Table 1 describes the properties of the MODIS spectral bands [59].

3.1.2. Landsat Satellite Data

Surface Reflectance retrievals from L5—Thematic Mapper (TM) and L8—Operational Land Imager (OLI) data are provided by the United States Geological Survey (USGS) and are available in GEE (see Table 1). Here, we used the USGS L5-TM and L8-OLI Surface Reflectance Tier 1 dataset to extract glacier areas to evaluate our results of the glacier delineation from MODIS data in 2000 and 2020 accordingly. Cloudless L5-TM in 2000 and L8-OLI in 2020 were selected for validating MODIS glacier delineation. As the reflectance of snow and ice is high in the visible band and near zero in the short-wave infrared band, bands 2 and 5 of L5-TM and bands 3 and 6 of L8-OLI were used to delineate snow and glaciers. The method to delineate glaciers was similar to the one applied to MODIS data.

3.1.3. Boundary of the Tarim Basin

The boundary of the Tarim Basin was defined using the HydroSHEDS data to determine the boundaries of the sub-basins and of the entire Tarim Basin. HydroSHEDS data provide hydrographic information on global and local scales [60] available at https://hydrosheds.org (accessed on 1 October 2021). The HydroSHEDS database provides a suite of geo-referenced datasets, including watershed boundaries, drainage directions, channel networks, flow accumulations, and river topology details [61]. The data are a good source to depict the catchment areas or the watershed boundaries on different scales.

3.1.4. Glacier Inventory

The Randolph Glacier Inventory (RGI 6.0) were applied as a reference to evaluate the glacier delineation based on MODIS data. The RGI 6.0 provided debris-covered information in our study area (Table 2). It was estimated that 3.91% of glaciers are covered by debris in Tarim Basin [42]. The RGI 6.0 dataset was motivated by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC AR5) [62], which utilized most of the Landsat Thematic Mapper (TM) and Enhanced TM Plus (ETM+) images between 2003 and 2010 [63] to delineate glaciers in China.

3.1.5. SRTM DEM

The Shuttle Radar Topography Mission (SRTM) DEM is a 30 m resolution void-filled digital elevation model, produced by NASA based on the STRM Interferometric Synthetic Aperture Radar (InSAR) data. This dataset was created and made available by the USGS and the German Aerospace Center (DLR). The accuracy of SRTM DEM is regarded as better than 20 m in the horizontal and 16 m in the vertical direction expressed as linear error at a 90% confidence level [64]. The SRTM DEM dataset was generated using data in February 2000, while our glacier area was determined for each year between February 2000 and February 2020. In this study, we employed the SRTM DEM to acquire the elevation attributes of glaciers, as many earlier studies where the SRTM DEM was applied to support analyses on the distribution of glaciers with elevation [65,66].

3.1.6. Meteorological Data

In order to analyze regional climate forcing, the China Meteorological Forcing Dataset (CMFD) was downloaded from the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn/zh-hans/data/8028b944-daaa-4511-8769-965612652c49, accessed on 3 January 2021). The CMFD is a high spatial-temporal resolution gridded near-surface meteorological dataset developed specifically for studies on land surface processes [67]. The dataset was generated through the fusion of remote sensing products, reanalysis data and in situ weather station data. The dataset includes 2 m air temperature and precipitation rate with a temporal resolution of three hours and a spatial resolution of 0.1° [68]. The annual mean temperature and precipitation from 1979 to 2018 were calculated for each sub-basin (see Section 5.2) in our study.

3.2. Method

Google Earth Engine is suitable to develop and implement long-term and large area time series remote sensing applications; its high-performance computation capabilities offer a state-of-the-art platform to explore and process earth observation data without actually downloading such large datasets [52].
To explore the relationship of snow cover with glacier area in different regions of the Tarim Basin, we evaluated snow cover based on the MODIS 8-day composite images between 18 February 2000 and 25 January 2021. Summer is the best season to separate glaciers from snow cover through satellite images because snow not contributing to glacier mass balance melts in summer due to very few snow events occurring in summer in High Mountain Asia [49]. However, glacier mapping is challenging because clouds frequently appear during summer in the Tarim Basin. To mitigate this problem, we developed a summer composite image to delineate glaciers, as done in [69]. We use summer observations to separate the glacier and snow cover area based on the minimum snow cover during summer, i.e., from July to September, to delineate glaciers [49].
The method primarily includes two parallel workflows: (1) glacier extraction and (2) snow cover extraction. The glacier extraction workflow includes: selection and retrieval of datasets, generating a summertime composite image, annual glacier extraction, water mask generation to exclude water pixels from glacier, and glacier area delineation. Most of the water pixels were located close to the glacier terminal or to the glacier boundary and can be excluded using NDSI, but there were still a few misidentified pixels. To ensure the accuracy of estimated glacier changes, we masked off the water pixels. The snow cover extraction workflow includes: retrieval of 8-day MODIS composite datasets, identifying pixels with end of summer snow cover and delineation of intermittent snow cover by thresholding the median NDSI (see Section 3.2.1 for details), water mask generation excluding water pixels from snow cover, and snow cover area calculation. The glacier area thus delineated includes end of summer snow cover, which contributes to ice accumulation. Figure 2 illustrates the two parallel work-flows.

3.2.1. Extracting Glacier Area

The boundaries of the Tarim Basin and its sub-basins were defined based on HydroSHEDS data (see Section 3.1.3). We uploaded the Boundary shapefile to the GEE driver and defined it as Region of Interest. The Date Range was used to select datasets within a specific time period.
The Datasets applied to annual delineation of glaciers were MODIS (called ‘MODIS/006/MOD09A1’ in GEE). The coarse spatial resolution of MODIS can only capture glaciers larger than 0.25 km2 and it is likely to overestimate the rate of decrease in small glaciers. Glaciers were extracted using summer composites of MODIS surface reflectance data. For example, to extract the glacier area in 2000, using the GEE Image Collection function, we selected, accessed and processed the MOD09A1 images from 1 July 2000 to 30 September 2000, i.e., the summer period to minimize the impact of snow cover. We called this step “median surface in-band reflectance”, which was applied to generate a summer composite image for each year. All reflectance retrievals for each pixel contributed to determine the median reflectance during summer. This compositing procedure image yielded combinations of in-band reflectance such that NDSI was lower when observing rapidly varying surface features, such as clouds and seasonal snow cover, and higher for stable snow and glaciers. Next, we calculated the NDSI based on the composite image.
NDSI can be applied to discriminate glacier/snow from many other land cover features and decreases the influence of atmospheric effects and viewing angles. This is achieved by normalizing the spectral difference between MODIS band 4 and band 6 to their sum [70]:
N D S I = G r e e n ( M O D I S   b a n d 4 ) S W I R ( M O D I S   b a n d 6 )   G r e e n ( M O D I S   b a n d 4 ) + S W I R ( M O D I S   b a n d 6 )
Glaciers were extracted by applying a threshold to the NDSI calculated with median summertime values of the at-surface reflectance in MODIS bands 4 and 6 (see Equation (1)). Pixels with NDSI ≥ 0.4 are identified as being glacier or stable snow [42,44]. A threshold of 0.11 on an NIR band is widely used to prevent water pixels from being classified as glaciers [37,48,71]. Furthermore, the Primary Glacier Boundary (Figure 2) is generated by applying these threshold settings. This operation is repeated for each year.
The snow products are affected, however, by multiple limitations, such as the misclassification of water pixels [44]. For this reason, we removed water pixels when re-generating glacier and snow cover areas. The water pixels were identified by applying a threshold to the Normalized Difference Water Index (NDWI). In addition, since melting occurs at the glacier surface during summer, liquid water may also exist on glaciers and will give reflectance spectra similar to liquid water. In order to further diminish the impact of water pixels, we generated a water mask based on the summertime composite image. The water pixels were identified by applying the Normalized Difference Water Index [72], due to the strong water absorption in the near-infrared band and relatively higher reflectance in the green band. The NDWI [72] is defined as:
N D W I = G r e e n ( M O D I S   b a n d 4 ) N I R ( M O D I S   b a n d 2 )   G r e e n ( M O D I S   b a n d 4 ) + N I R ( M O D I S   b a n d 2 )
Pixels are identified as water when NDWI > 0.2 and the near-infrared reflectance in band 2 is smaller than 0.2 [37,48]. This step is to generate a water mask according to the summer composite images and further remove the water pixels. The annual glacier area was calculated after removing the water mask. The annual glacier boundary was generated from 2000 to 2020 using the composite images during summer.
As we mentioned in Section 3.1.2, the L5-TM and L8-OLI Surface Reflectance product (namely, ‘LANDSAT/LT05/C01/T1_SR’ and ‘LANDSAT/LC08/C01/T1_SR’ in GEE) data were used to evaluate the delineation of glaciers based on the lower resolution MODIS data. Due to the limited number of cloudless L5-TM scenes, we used L5-TM scenes from 1 January 2000 to 30 December in 2000 to ensure the coverage of study area. Cloudless L8-OLI scenes from 1 July 2020 to 30 September in 2020 were mosaiced in GEE.
Next, the glacier area is obtained after removing the water mask in 2000 and 2020, respectively. We compared the glacier area in 2000/2020 generated by MODIS with the glacier area in 2000/2020 generated by L5-TM/L8-OLI data to evaluate our glacier delineation based on MODIS data (Figure 3).

3.2.2. Extracting Snow Cover

In this study, we separated the glacier and intermittent snow cover mainly based on the stability of summer snow cover using frequent observations. We considered summer as the season with the smallest snow area, and thus considered the area delineated by thresholding the median NDSI during summer as the glacier and stable snow cover area. Intermittent snow cover was calculated using data acquired at 8-day intervals during winter, i.e., from 18 February 2000 to 25 January 2021, and using the MODIS data accessed via GEE. Snow cover was delineated by applying the same threshold, as done with glaciers, to the 8-day NDSI images calculated with each 8-day MODIS in-band at-surface reflectance in bands 4 and 6. The snow cover generation method is similar to the process of glacier generation. The primary snow cover was generated when NDSI ≥ 0.4 and band 2 > 0.11, then the intermittent snow cover was the total number of pixels identified as snow in all available images during winter. The water mask was also generated when NDWI > 0.2 and the near-infrared reflectance in band 2 < 0.2. Next, the water mask was excluded from the first estimate of intermittent snow cover. The calculated intermittent snow cover area was output in the GEE table format. Finally, the glacier area, determined with the summertime data, was excluded from the first estimate of the intermittent snow cover area within the same year.

4. Results

4.1. Accuracy Assessment

Glaciated areas may be misclassified as intermittent snow due to similar spectral characteristics [73]. In our proposed method, intermittent snow is identified by evaluating the temporal variability of snow cover for each pixel and relating it to the median surface spectral reflectance. Estimating the glacier area accuracy is challenging, as reference data are limited. Here, glacier area delineated with the MODIS images by applying the NDSI threshold (see Section 3) was evaluated against the area of five subregions with the boundaries of glaciers in 2000 delineated with L5-TM and as with L8-OLI image data in 2020 (Table 2). The delineation with MODIS data of glaciers in different subregions was also compared with previously published results (see Section 5.1).
Glacier delineation by a spectral band combination is supposed to be one of the most efficient methods applicable to debris-free glaciers, but it is not appropriate to debris-covered glaciers [74]. As mentioned in the Introduction, the debris-covered glacier portions were not taken into account in our study due to the limited incidence of changes in the extension of debris facies on the total change in glacier area of the study region [42,43,75,76]. According to the published RGI 6.0, debris-covered glaciers take up 3.91% of the total glacier area in the Tarim Basin, i.e., errors of estimate have limited effect on the glacier area result.
The uncertainty in the MODIS 2000 glacier area was higher in East Tien Shan and East Kunlun, when compared with the estimates based on the L5-TM image data (Table 2). The results for the Pamir and West Kunlun regions were the most accurate, with −6.76% and −7.94% relative differences when comparing the MODIS 2000 estimates with the L5-TM 2000 data. The comparison of L8-OLI delineated glacier area and MODIS glacier area in 2020 indicated a higher uncertainty in East Tien Shan and Pamir. The results for the West Tien Shan and East Kunlun regions were the most accurate, with 2.55% and 6.19%.
To estimate the rate of change of glacier area, we fitted a linear relationship to the annual values in Figure 4 and Figure 5, then we took the slope of this relationship as estimated rate of change. To evaluate the accuracy of these estimates, we calculated the relative Root Mean Square Deviation as:
R M S D = i = 1 N ( x i x i ¯ ) 2 N
where x i is the actual value, x i ¯ is the linear regression value, and N is the number of non-missing data points. The Relative Root Mean Square Deviation was calculated is then:
R R M S D = R M S D / A v e r a g e ( x i ¯ )
The RRMSD has been added in Figure 4 and Figure 5. The error of estimate on glacier area in a region was −6.48% (weighted average of all MODSI cases in Table 2).

4.2. Variability of Glacier Area in the Tarim Basin

The annual glacier area in the five subregions of the Tarim Basin was retrieved from MODIS data from 2000 to 2020 by applying the method described in Section 3 (Figure 4). In the whole Tarim Basin, the glacier area decreased by 7975.71 km2 from 2000 to 2020, i.e., 25.11% of the initial glacier area, at a rate of 0.94%/a, and it decreased fastest during 2016–2020 at a rate of approximately 1.82%/a. In other words, the decreasing rate accelerated after 2016.
The evolution of glacier area in the five subregions of the Tarim Basin was further analyzed separately (Figure 5). It was found that glacier area decreased in all the five subregions during the 2000–2020 period (Figure 5), but there were clear differences in the evolution across the five subregions.
The rate of change in clean glacier area in E Tien Shan, W Tien Shan, Pamir, W Kunlun, and E Kunlun was 2.98%/a, 1.07%/a, 0.5%/a, 1.04%/a, and 0.81%/a, respectively, during 2000–2020. These rates imply that the total reduction in 2000–2020 was largest, i.e., 50.86%, i.e., 154.13 km2, in the E Tien Shan area. The area of glaciers in W Tien Shan decreased by 973.76 km2 (13.02%). In the Pamir region, the area of glacier decreased by 20.52%. The reduction in the W Kunlun and E Kunlun regions were 30.98% and 36.49%, respectively.

4.3. Changes in Glacier Area at Different Elevations

Most glaciers in the Tarim Basin are located at elevations in the range 3000–7500 m (Figure 6 and Figure 7). Many small glaciers are scattered at these elevations in the Tien Shan region. In the Pamir and West Kunlun, glacier elevation ranges from 5000 to 7500 m. West Kunlun is the largest glacial region in the Tarim Basin. In East Kunlun, glaciers are distributed above 6000 m. The terminus position of most glaciers is above 3000 m, while the snow line is between 3600 and 5700 m [17]. In general terms, temperature decreases and precipitation increases with elevation, but spatial variability is significant and glacier response to climate variability needs to be evaluated at different elevations.
Glacier area in the 5500–6000 range is the largest fraction of total glacier area in the Tarim Basin, and the area was 9150.77 km2 in 2000, and decreased to 7434.62 km2 in 2020. Our analysis shows that the decrease in glacier area was largest in the 3500–5500 m range: the glacier area was 18,358.29 km2 in 2000, and decreased to 12,568.57 km2 in 2020. Glacier area in the 6000–6500 range was 3926.68 km2 in 2000, and decreased slightly to 3511.61 km2 in 2020. However, the area of glaciers above 6500 remained stable between 2000 and 2020, and even increased slightly, although the glacier area above 6500 m accounts for 1.08% of total glacier area. The stability is most likely due to the low temperature at higher elevations, while air temperature increased at lower elevations.

4.4. Snow Cover Area in Sub-Basins

Glacier area is determined by mass accumulation and ablation, with the dominant contributor to ice accumulation being snowfall [77]. Snowfall variability can be observed through intermittent snow cover. Snow cover can be detected by optical remote sensing because of its distinct spectral signature [48], but frequent cloud cover and mixed water pixels in high mountain regions need to be dealt with. As we mentioned in Section 3.2, we applied a water mask to minimize the water effect and used 8-day composite surface reflectance products to minimize the impact of cloud cover in the retrieval of intermittent snow cover. The temporal evolution of snow cover area showed a typical monthly cycle with lowest area observed in July and August (Figure 8). The fastest expansion of snow cover happened in October and the fastest reduction occurred in April (Figure 8). Overall, no significant trend could be observed in the time series at full temporal resolution in Tarim Basin from 2000 to 2020.
The seasonality of snow cover in each subregion was better captured by the mean monthly snow cover area (Figure 9). The snow cover area started increasing in August and grew until February in the following year. At a lower temperature, the water vapor carried by the westerlies precipitated in solid form in the high mountain regions, thus causing an increase in snow cover area. The snow cover area decreased from February to August, while snow melting already started at lower elevations in February. Snow cover area in the E Kunlun, however, already started decreasing in January, which might be related to the strengthening westerlies [31]. In the W Tien Shan region, the annual evolution of snow cover was slightly different, with the largest snow cover area being reached in February. Also worth noting was the lower rate of decrease, especially from April to July, in the W Kunlun when compared with other subregions. In addition, the snow cover rate of increase from August to December was higher than in other subregions.
To detect and evaluate a possible multi-annual trend in snow cover, we defined the Remaining Snow Index (RSI) as:
R S I = A m i n A m a x × 100 %
where ‘ A m a x ’ is the yearly maximum value of the monthly snow cover area, with the maximum value usually being reached in January or February (Figure 9). The ‘ A m i n ’ is the yearly non-zero minimum value of the monthly snow cover area, usually being reached in July or August (Figure 9). The Remaining Snow Index is an indicator of the snow remaining annually at the end of the yearly cycle; thus, it is related to mass accumulation in glaciers, but it is a relative, not an absolute, indicator of snow area. Here, we used the Remaining Snow Index to analyze the snow variation in 2000–2020 in each subregion (Figure 10), which could not be captured by the data plotted in Figure 8. The Remaining Snow Index in W Kunlun was the largest among the five regions, which means the W Kunlun accumulated more snow. RSI does not indicate the area of remaining snow; it reflects the snow cover change within one year. For a given minimum snow cover area, the RSI is larger when the snow area is smaller in winter. Accordingly, an increasing trend of RSI does not imply that the area of remaining snow increases, only that the fraction of minimum snow area in summer over maximum snow area in winter is increasing. Further detailed discussion can be found in Section 5.2.

5. Discussion

5.1. Comparison with Previous Studies

Other authors have studied changes in glacier area in the Tarim Basin, and we evaluated our findings against this literature evidence. To date, regional investigations on changes in glacier area were limited in spatial coverage, so it is impossible to verify our results in detail. The comparisons are performed using reported results on the annual rate of change in glacier area within slightly different regions and period of time. In these studies, the glacier area shows a diversified variability in the W Tien Shan, E Tien Shan, Pamir, W Kunlun, and E Kunlun (Table 3).
A decrease in glacier area was recorded in the W Tien Shan during 1970–2018. The rate of shrinkage in glacier area varied from 0.11%/a to 3.7%/a within different sub-basins, with the Sary-Jaz River Basin having the largest reduction rate during 1990–2010 [78]. The reduction rate of the Sary-Jaz River Basin was slightly larger than our study result of 1.07%/a. Our study area included the Central Tien Shan, while other studies paid more attention to the China part of the Tien Shan [6,19,79]. Some disagreement in the rate of decrease estimated by different studies may be due to differences in the area of interest [12]. In addition, the large number of small glaciers (<1 km2) [19] may be another reason for the discrepancy in estimated annual rate. Differences in the observation period may clearly lead to significant differences, especially between earlier and recent studies, since the glacier rate of decrease may have been accelerating due to increasing temperature and reduced precipitation in recent years [49]. In addition, most studies were focused on one or two river catchments of the Tarim Basin, while our study covers the entire Tarim Basin. The glacier area shrinkage rate was 0.8%/a between 1990 and 2018 for the E Tien Shan in the study of Huang et al. (2021) [49], while the glaciers shrank 2.98%/a for the same area during 2000–2020 according to our results. In the Pamir and East Kunlun region, the rate of glacier decrease was smaller than in the East Tien Shan region. We only consider the part of the Pamir region draining into Tarim Basin. The sub-basin within Pamir shrank by 0.05% to 0.59%/a during 1990–2018 [49,80,81,82], while our result of 1.01%/a is larger than previous studies. In Pamir and W Kunlun, the glacier snouts have been retreating, stable or advancing [49,80]. The smallest rate of glacier contraction was observed in the Pamir regions [31]. The deviation is possibly caused by the different observation period, i.e., the rate of decrease was 0.46 during 2000–2010, while it was 0.59 during 1990–2018 according to Huang et al. (2021). More comparisons in the five subregions and periods are shown in Table 3.

5.2. Glacier Area Variation in Response to Climate Forcing and Snow Cover

Since the end of the 1980s, a transition from warm-dry to warm-wet conditions with both increasing temperatures and precipitation was observed in northwestern China [87]. Many factors have caused a reduction in glacier area, but precipitation and temperature were the most significant, where the summer temperature mainly regulated glacier melting and winter precipitation mainly regulated the accumulation of ice. Here, annual precipitation and temperature between 1979 and 2018 in each sub-basin are evaluated.
In E Tien Shan, both annual precipitation and temperature increased during 1979–2004, while they decreased during 2004–2018 (Figure 11a). Zhang et al. (2019) [19] showed that the temperature and precipitation in the Aksu River Basin (W Tien Shan) increased in 1975–2016, and the warming temperatures were the main reason for the glacier recession in Tien Shan [19]. A positive trend in annual precipitation and temperature was noted in W Tien Shan (Figure 11b). Increases in temperature and precipitation have led to an elevation-dependent snow cover response, with decreased snow cover area at lower elevations and increased snow cover area at higher elevations [88]. The increase in precipitation has been induced by increasing water vapor transport by the westerly airflow [89]. Westerly-affected regions such as the Pamir and W Kunlun have extensive winter snow cover [90,91]. The smallest rate of glacier contraction was observed in the Pamir, with a reduction rate of 0.50%/a [31]. Positive mass balances observed from geodetic studies (using DEMs and ICESat data) in W Kunlun in 2006–2010 were reported by [12].
According to [55], the mean annual precipitation in the southern part of the Tarim Basin is lower than in the northern and western regions. Higher temperature and higher water vapor content in the atmosphere may lead to larger mass accumulation than in earlier years. In the E Kunlun, both temperature and precipitation increased from 1979 to 2004, and decreased from 2004 to 2018. We found that the glacier area decreased from 2000 to 2004 (Figure 5e), which is in consistent with the noted trends in climatic condition.
No trend was detected in the annual snow-covered area, with significant interannual variability highlighted by calculating the RSI (see Figure 10) and generally a small fraction of winter snow-covered area remaining after summer.
We also analyzed the normalized cross-correlation functions of glacier area and temperature, glacier area and precipitation, and glacier area and snow cover in the five Tarim subregions. The detailed results are presented in the Supplementary Materials (Figure S1). We calculated the normalized cross-correlation function by shifting step by step the response time series with respect to the forcing time series ten years forwards and ten years backwards, using the matlab function (xcorr(x,y, ‘coeff’)). The results show that the cross-correlation is symmetric around a zero time-lag where the correlation is highest. This indicates that the response of glacier area to temperature, precipitation and snow area is within less than a year, which is the temporal resolution of the data applied in this analysis. Precipitation and snow cover show positive correlation with glacier area, as expected, while the correlation of glacier area with temperature is positive in W Tien Shan and negative in other subregions.
The mean elevation of the glacier area in the Tien Shan region is lower than in the W Kunlun region, which increases the vulnerability to climate variability of the glaciers in the former region. Fan et al. (2014) [92] concluded that the spatial variability of temperature and precipitation in the high elevation mountain zone is significant, while the lower elevation zone (3500 m a.s.l.) is uniformly affected by less precipitation, higher temperature and large irradiance [89]. Most glaciers in the Tien Shan region are distributed at elevation between 3000 and 5500 m and experienced a continuous glacier recession. In our study, the absolute loss is highest in W Kunlun, and the highest relative rate of decrease is found in the E Tien Shan. Shangguan et al. (2009) [17] found that the altitude of the maximum glacier reduction was 4900–5400 m a.s.l., while another peak in the decrease of glacier area was found at an altitude of 4100–4500 m a.s.l. Our results also show that the glacier area varies greatly between 4000 and 5500 m a.s.l. Shangguan et al. (2009) [17] also found little change in the area of glaciers at an altitude of 6100 m, which is not consistent with our result.
Some studies used hydrologic modelling to describe the response of glaciers to temperature and precipitation [93,94,95]. However, the temperature and precipitation (liquid water) data used for hydrologic modelling mainly depend on weather stations scattered in the high elevation mountain area [93,94]. Additionally, elevation, monsoon, and terrain vary significantly in the Tarim region [96]. The widely noted changes in glacier area raised the awareness of the urgency of studying the mountain glaciers in the Tarim Basin [97,98]. The estimated response of glacier area by hydrologic modelling, however, is affected by the coarse spatial resolution and by uncertainty in the parameterization of terrain and elevation effects [96].

5.3. Limitations and Advantages of the Method

The MODIS data products take advantage of the temporal frequency. In this study, we used this data product to evaluate the annual change rate in glacier area and in the seasonal snow cover area. This is a systematic way to explore the relationship between glacier and snow cover areas. The MODIS coarse spatial resolution, however, limits the accuracy in the estimation of glacier area. The estimated total glacier area is affected by the minimal glacier size in relation to the spatial resolution of applied image data. The minimal size of detectable glaciers with MODIS data is 0.25 km2 (500 m × 500 m). Our estimates of the area of small glaciers and of their distribution may differ from results based on high- or medium-resolution data. In addition, the examples of Muztag and Kungur Tagh were used to show the limitations and advantages of the proposed NDSI based extraction method. We used the standard NDSI threshold of 0.4 to extract glacier area and intermittent snow cover in our study. In Figure 3, the extracted glacier outline does not correctly capture the debris-covered region, which may cause an underestimation of glacier area [74]. In W and E Kunlun, the fractional debris-covered area is less than 5% of the glacier area, so the effect of debris is limited [42]. The limitation of our study is that we did not map the debris-covered region. A glacier is regarded as “perennial mass of ice, and possibly firn and snow, originating on the land surface by the recrystallization of snow or other forms of solid precipitation and showing evidence of past or present flow” [52]. “Firn is snow that has survived at least one ablation season but has not been transformed to glacier ice” [52]. However, glacier mapping is challenging because clouds frequently appear during summer in the Tarim Basin. To mitigate this problem, we developed a summer composite image to delineate glaciers, as done in [53]. We choose the data acquired in July and September to delineate the glacier area, but in high altitude regions, such as W Kunlun (more than 5500 m a.s.l.), the perennial snow (“snow that has survived at least one ablation season but has not been transformed yet to glacier ice” [99]) can remain for years at high elevation. We have evaluated the impact of perennial snow on our results as follows.
The perennial snow appears to be adjacent to the glaciers, as delineated in the RGI 6.0 (Figure 12). The perennial snow was mapped annually in the entire Tarim Basin. The perennial snow that occurred in more than 5 years out of 20 was analyzed to determine its frequency distribution (Figure 13). The snow persistent for the entire observation period, i.e., 2000–2020, was already counted in our analysis of the evolution of glacier area and has no impact on our results. The snow with shorter persistence was not present continuously and cannot be counted as glacier area.

6. Conclusions

Glacier area extraction from satellite images is an important aspect of glacier mapping and studies on glacier changes. The selection of cloud-free images is difficult when using low revisit images such as the ones acquired by L8-OLI. Our method relies on the frequency of usable image pixels to minimize the impact of clouds and to detect intermittent snow. The major innovation in our method is the use of the median surface reflectance to composite several images pixel by pixel. This does not require that cloud-free observations of all pixels are available on the same day, just that a sufficient number of cloud-free samples is acquired for each pixel. The composite image improves the reliability of estimated change in glacier area, because the median is a robust statistic. The NDSI method combined with threshold setting was effective, leading to systematic estimates of the glacier area and monthly snow cover.
Changes in glacier and intermittent snow area between 2000 and 2020 were first evaluated for the entire Tarim Basin. The results show that the glacier area decreased by 7975.71 km2 (25.11%), with the annual rate of decrease being 0.94%/a. Most glaciers are distributed between 5000 and 6000 m a.s.l., where the glacier area fluctuates most. Above 6000 m a.s.l., the glacier area was rather stable. The consistent spatial distribution of perennial snow suggests that the RGI 6.0 boundaries may need to be revisited.
The Tarim Basin glaciers show considerable heterogeneity in the five different subregions due to the different circulation of the atmosphere and environmental conditions. Overall, the five subregions show a decreasing trend during the 2000–2020 period. In the Tien Shan region, the glacier area decreased rapidly. The intra-annual variability of snow cover was largest for W Kunlun and smallest in the E Tian Shan. The analysis of the cross-correlation functions showed that the glacier area responds to temperature, precipitation and snow cover within the same year. Further studies on in-depth analysis of driving factors of glacier area change are required.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs13245117/s1, Figure S1: Normalized Cross-correlation function of Glacier Area and Temperature, Glacier Area and Precipitation, Glacier Area and Snow Cover vs. time lag in five Tarim subregions.

Author Contributions

Conceptualization, J.Z. (Jing Zhang) and L.J.; methodology, J.Z. (Jing Zhang); software, J.Z. (Jing Zhang); validation, J.Z. (Jing Zhang); formal analysis, J.Z. (Jing Zhang), L.J. and M.M.; investigation, J.Z. (Jing Zhang); resources, L.J. and M.M.; data curation, J.Z. (Jing Zhang) and L.J.; writing—original draft preparation, J.Z. (Jing Zhang); writing—review and editing, L.J., M.M. and J.Z. (Jie Zhou), S.R.; visualization, J.Z. (Jing Zhang); supervision, L.J. and M.M.; project administration, L.J.; funding acquisition, L.J. and M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was jointly funded by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDA19070102), the National Natural Science Foundation of China project (grant no. 91737205), and the Chinese Academy of Sciences President’s International. Fellowship Initiative (grant no. 2020VTA0001); the MOST High Level Foreign Expert program (grant no. GL20200161002).

Acknowledgments

We thank the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn/zh-hans/data/8028b944-daaa-4511-8769-965612652c49, accessed on 3 January 2021) for offering the meteorological data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rumbaur, C.; Thevs, N.; Disse, M.; Ahlheim, M.; Brieden, A.; Cyffka, B.; Duethmann, D.; Feike, T.; Frör, O.; Gärtner, P.; et al. Sustainable management of river oases along the Tarim River (SuMaRiO) in Northwest China under conditions of climate change. Earth Syst. Dyn. 2015, 6, 83–107. [Google Scholar] [CrossRef] [Green Version]
  2. Shi, Y.; Liu, C.; Kang, E. The Glacier Inventory of China. Ann. Glaciol. 2009, 50, 1–4. [Google Scholar] [CrossRef] [Green Version]
  3. Duethmann, D.; Menz, C.; Jiang, T.; Vorogushyn, S. Projections for headwater catchments of the Tarim River reveal glacier retreat and decreasing surface water availability but uncertainties are large. Environ. Res. Lett. 2016, 11, 054024, Erratum in 2019, 11, 049501. [Google Scholar] [CrossRef] [Green Version]
  4. Hagg, W.; Mayer, C.; Lambrecht, A.; Kriegel, D.; Azizov, E. Glacier changes in the Big Naryn basin, Central Tian Shan. Global Planet. Change 2013, 110, 40–50. [Google Scholar] [CrossRef]
  5. Kaser, G.; Grosshauser, M.; Marzeion, B. Contribution potential of glaciers to water availability in different climate regimes. Proc. Natl. Acad. Sci. USA 2010, 107, 20223–20227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Pieczonka, T.; Bolch, T.; Junfeng, W.; Shiyin, L. 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]
  7. Chen, Y.; Li, W.; Fang, G.; Li, Z. Review article: Hydrological modeling in glacierized catchments of central Asia—Status and challenges. Hydrol. Earth Syst. Sci. 2017, 21, 669–684. [Google Scholar] [CrossRef] [Green Version]
  8. Yongping, S. Response of Glacier Flash Flood to Global Warming in Tarim River Basin. Clim. Chang. Res. 2006, 2, 32–35. [Google Scholar]
  9. Bazai, N.A.; Cui, P.; Carling, P.A.; Wang, H.; Hassan, J.; Liu, D.; Zhang, G.; Jin, W. Increasing glacial lake outburst flood hazard in response to surge glaciers in the Karakoram. Earth-Sci. Rev. 2021, 212, 103432. [Google Scholar] [CrossRef]
  10. Lv, M.; Guo, H.; Yan, J.; Wu, K.; Liu, G.; Lu, X.; Ruan, Z.; Yan, S. Distinguishing Glaciers between Surging and Advancing by Remote Sensing: A Case Study in the Eastern Karakoram. Remote Sens. 2020, 12, 2297. [Google Scholar] [CrossRef]
  11. Lv, M.; Guo, H.; Lu, X.; Liu, G.; Yan, S.; Ruan, Z.; Ding, Y.; Quincey, D.J. Characterizing the behaviour of surge- and non-surge-type glaciers in the Kingata Mountains, eastern Pamir, from 1999 to 2016. Cryosphere 2019, 13, 219–236. [Google Scholar] [CrossRef] [Green Version]
  12. Bolch, T.; Shea, J.M.; Liu, S.; Azam, M.F.; Gao, Y.; Gruber, S.; Immerzeel, W.W.; Kulkarni, A.; Li, H.; Tahir, A.A.; et al. Status and Change of the Cryosphere in the Extended Hindu Kush Himalaya Region. In The Hindu Kush Himalaya Assessment; Springer: Cham, Switzerland, 2019; pp. 209–255. [Google Scholar] [CrossRef] [Green Version]
  13. Rabatel, A.; Dedieu, J.-P.; Vincent, C. Using remote-sensing data to determine equilibrium-line altitude and mass-balance time series: Validation on three French glaciers, 1994–2002. J. Glaciol. 2005, 51, 539–546. [Google Scholar] [CrossRef] [Green Version]
  14. Rabatel, A.; Letréguilly, A.; Dedieu, J.-P.; Eckert, N. Changes in glacier equilibrium-line altitude in the western Alps from 1984 to 2010: Evaluation by remote sensing and modeling of the morpho-topographic and climate controls. Cryosphere 2013, 7, 1455–1471. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, X.; Chen, H.; Chen, Y. Topography-Related Glacier Area Changes in Central Tianshan from 1989 to 2015 Derived from Landsat Images and ASTER GDEM Data. Water 2018, 10, 555. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.; Hou, S.; Huai, B.; An, W.; Pang, H.; Liu, Y. Glacier anomaly over the western Kunlun Mountains, Northwestern Tibetan Plateau, since the 1970s. J. Glaciol. 2018, 64, 624–636. [Google Scholar] [CrossRef] [Green Version]
  17. Shangguan, D.; Liu, S.; Ding, Y.; Ding, L.; Xu, J.; Jing, L. Glacier changes during the last forty years in the Tarim Interior River basin, northwest China. Prog. Nat. Sci. 2009, 19, 727–732. [Google Scholar] [CrossRef]
  18. Liu, S.; Ding, Y.; Shangguan, D.; Zhang, Y.; Li, J.; Han, H.; Wang, J.; Xie, C. Glacier retreat as a result of climate warming and increased precipitation in the Tarim river basin, northwest China. Ann. Glaciol. 2006, 43, 91–96. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Q.; Chen, Y.; Li, Z.; Li, Y.; Xiang, Y.; Bian, W. Glacier changes from 1975 to 2016 in the Aksu River Basin, Central Tianshan Mountains. J. Geogr. Sci. 2019, 29, 984–1000. [Google Scholar] [CrossRef] [Green Version]
  20. Kääb, A.; Berthier, E.; Nuth, C.; Gardelle, J.; Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nat. Cell Biol. 2012, 488, 495–498. [Google Scholar] [CrossRef] [PubMed]
  21. Seong, Y.B.; Owen, L.A.; Caffee, M.W.; Kamp, U.; Bishop, M.P.; Bush, A.; Copland, L.; Shroder, J.F. Rates of basin-wide rockwall retreat in the K2 region of the Central Karakoram defined by terrestrial cosmogenic nuclide Be-10. Geomorphology 2009, 107, 254–262. [Google Scholar] [CrossRef]
  22. Bocchiola, D.; Diolaiuti, G. Recent (1980–2009) evidence of climate change in the upper Karakoram, Pakistan. Theor. Appl. Climatol. 2013, 113, 611–641. [Google Scholar] [CrossRef]
  23. Forsythe, N.; Fowler, H.; Li, X.-F.; Blenkinsop, S.; Pritchard, D. Karakoram temperature and glacial melt driven by regional atmospheric circulation variability. Nat. Clim. Chang. 2017, 7, 664–670. [Google Scholar] [CrossRef]
  24. Hasson, S.U.; Böhner, J.; Lucarini, V. Prevailing climatic trends and runoff response from Hindukush–Karakoram–Himalaya, upper Indus Basin. Earth Syst. Dyn. 2017, 8, 337–355. [Google Scholar] [CrossRef] [Green Version]
  25. Cannon, F.; Carvalho, L.M.V.; Jones, C.; Bookhagen, B. Multi-annual variations in winter westerly disturbance activity affecting the Himalaya. Clim. Dyn. 2015, 44, 441–455. [Google Scholar] [CrossRef] [Green Version]
  26. Norris, J.; Carvalho, L.M.V.; Jones, C.; Cannon, F. Deciphering the contrasting climatic trends between the central Himalaya and Karakoram with 36 years of WRF simulations. Clim. Dyn. 2019, 52, 159–180. [Google Scholar] [CrossRef]
  27. Norris, J.; Carvalho, L.M.V.; Jones, C.; Cannon, F. WRF simulations of two extreme snowfall events associated with contrasting extratropical cyclones over the western and central Himalaya. J. Geophys. Res. Atmos. 2015, 120, 3114–3138. [Google Scholar] [CrossRef] [Green Version]
  28. De Kok, R.J.; Tuinenburg, O.A.; Bonekamp, P.N.J.; Immerzeel, W.W. Irrigation as a Potential Driver for Anomalous Glacier Behavior in High Mountain Asia. Geophys. Res. Lett. 2018, 45, 2047–2054. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Kapnick, S.; Delworth, T.; Ashfaq, M.; Malyshev, S.; Milly, P.C.D. Snowfall less sensitive to warming in Karakoram than in Himalayas due to a unique seasonal cycle. Nat. Geosci. 2014, 7, 834–840. [Google Scholar] [CrossRef]
  30. 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] [Green Version]
  31. 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]
  32. Paul, F.; Bolch, T.; Kääb, A.; Nagler, T.; Nuth, C.; Scharrer, K.; Shepherd, A.; Strozzi, T.; Ticconi, F.; Bhambri, R.; et al. The glaciers climate change initiative: Methods for creating glacier area, elevation change and velocity products. Remote Sens. Environ. 2015, 162, 408–426. [Google Scholar] [CrossRef] [Green Version]
  33. Hall, D.K.; Ormsby, J.P.; Bindschadler, R.A.; Siddalingaiah, H. Characterization of snow and ice reflectance zones on glaciers using Landsat TM data. Ann. Glaciol. 1987, 9, 104–108. [Google Scholar] [CrossRef] [Green Version]
  34. Paul, F.; Kääb, A. Perspectives on the production of a glacier inventory from multispectral satellite data in Arctic Canada: Cumberland Peninsula, Baffin Island. Ann. Glaciol. 2005, 42, 59–66. [Google Scholar] [CrossRef] [Green Version]
  35. Dozier, J. Spectral signature of alpine snow cover from the landsat thematic mapper. Remote Sens. Environ. 1989, 28, 9–22. [Google Scholar] [CrossRef]
  36. Tang, B.-H.; Shrestha, B.; Li, Z.-L.; Liu, G.; Ouyang, H.; Gurung, D.R.; Giriraj, A.; Aung, K.S. Determination of snow cover from MODIS data for the Tibetan Plateau region. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 356–365. [Google Scholar] [CrossRef]
  37. Hall, D.K.; Riggs, G.A.; Salomonson, V.V. Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data. Remote Sens. Environ. 1995, 54, 127–140. [Google Scholar] [CrossRef]
  38. Rabatel, A.; Bermejo, A.; Loarte, E.; Soruco, A.; Gomez, J.; Leonardini, G.; Vincent, C.; Sicart, J.E. Can the snowline be used as an indicator of the equilibrium line and mass balance for glaciers in the outer tropics? J. Glaciol. 2012, 58, 1027–1036. [Google Scholar] [CrossRef] [Green Version]
  39. Brun, F.; Dumont, M.; Wagnon, P.; Berthier, E.; Azam, M.F.; Shea, J.M.; Sirguey, P.; Rabatel, A.; Ramanathan, A. Seasonal changes in surface albedo of Himalayan glaciers from MODIS data and links with the annual mass balance. Cryosphere 2015, 9, 341–355. [Google Scholar] [CrossRef] [Green Version]
  40. Raup, B.; Kääb, A.; Kargel, J.S.; Bishop, M.P.; Hamilton, G.; Lee, E.; Paul, F.; Rau, F.; Soltesz, D.; Khalsa, S.J.; et al. Remote sensing and GIS technology in the Global Land Ice Measurements from Space (GLIMS) Project. Comput. Geosci. 2007, 33, 104–125. [Google Scholar] [CrossRef]
  41. Huang, L.; Li, Z.; Tian, B.-S.; Zhou, J.-M.; Chen, Q. Recognition of supraglacial debris in the Tianshan Mountains on polarimetric SAR images. Remote Sens. Environ. 2014, 145, 47–54. [Google Scholar] [CrossRef]
  42. Scherler, D.; Wulf, H.; Gorelick, N. Global Assessment of Supraglacial Debris-Cover Extents. Geophys. Res. Lett. 2018, 45, 11798–11805. [Google Scholar] [CrossRef] [Green Version]
  43. Anderson, L.S.; Anderson, R.S. Debris thickness patterns on debris-covered glaciers. Geomorphology 2018, 311, 1–12. [Google Scholar] [CrossRef]
  44. Riggs, G.A.; Hall, D.K.; Román, M.O. Overview of NASA’s MODIS and Visible Infrared Imaging Radiometer Suite (VIIRS) snow-cover Earth System Data Records. Earth Syst. Sci. Data 2017, 9, 765–777. [Google Scholar] [CrossRef] [Green Version]
  45. National Snow & Ice Data Center. Available online: https://nsidc.org/ (accessed on 1 March 2021).
  46. Justice, C.; Townshend, J.; Vermote, E.; Masuoka, E.; Wolfe, R.; Saleous, N.; Roy, D.; Morisette, J. An overview of MODIS Land data processing and product status. Remote Sens. Environ. 2002, 83, 3–15. [Google Scholar] [CrossRef]
  47. Jin, S.; Sader, S.A. MODIS time-series imagery for forest disturbance detection and quantification of patch size effects. Remote Sens. Environ. 2005, 99, 462–470. [Google Scholar] [CrossRef]
  48. Hall, D.K.; Riggs, G.A.; Salomonson, V.V.; DiGirolamo, N.E.; Bayr, K.J. MODIS snow-cover products. Remote Sens. Environ. 2002, 83, 181–194. [Google Scholar] [CrossRef] [Green Version]
  49. Huang, L.; Li, Z.; Zhou, J.M.; Zhang, P. An automatic method for clean glacier and nonseasonal snow area change estimation in High Mountain Asia from 1990 to 2018. Remote Sens. Environ. 2021, 258, 112376. [Google Scholar] [CrossRef]
  50. Mishra, V.D.; Negi, H.S.; Rawat, A.K.; Chaturvedi, A.; Singh, R.P. Retrieval of sub-pixel snow cover information in the Himalayan region using medium and coarse resolution remote sensing data. Int. J. Remote. Sens. 2009, 30, 4707–4731. [Google Scholar] [CrossRef]
  51. Zhang, H.; Zhang, F.; Zhang, G.; Che, T.; Yan, W.; Ye, M.; Ma, N. Ground-based evaluation of MODIS snow cover product V6 across China: Implications for the selection of NDSI threshold. Sci. Total Environ. 2019, 651, 2712–2726. [Google Scholar] [CrossRef] [PubMed]
  52. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  53. Pandey, P. Inventory of rock glaciers in Himachal Himalaya, India using high-resolution Google Earth imagery. Geomorphology 2019, 340, 103–115. [Google Scholar] [CrossRef]
  54. Yang, X.; Chen, Y.; Wang, J. Combined use of Sentinel-2 and Landsat 8 to monitor water surface area dynamics using Google Earth Engine. Remote Sens. Lett. 2020, 11, 687–696. [Google Scholar] [CrossRef]
  55. Fang, G.; Yang, J.; Chen, Y.; Li, Z.; Ji, H.; De Maeyer, P. How Hydrologic Processes Differ Spatially in a Large Basin: Multisite and Multiobjective Modeling in the Tarim River Basin. J. Geophys. Res. Atmos. 2018, 123, 7098–7113. [Google Scholar] [CrossRef]
  56. Xu, Z.; Liu, Z.; Fu, G.; Chen, Y. Trends of major hydroclimatic variables in the Tarim River basin during the past 50 years. J. Arid. Environ. 2010, 74, 256–267. [Google Scholar] [CrossRef]
  57. Keilholz, P.; Disse, M.; Halik, Ü. Effects of Land Use and Climate Change on Groundwater and Ecosystems at the Middle Reaches of the Tarim River Using the MIKE SHE Integrated Hydrological Model. Water 2015, 7, 3040–3056. [Google Scholar] [CrossRef] [Green Version]
  58. Tao, H.; Gemmer, M.; Bai, Y.; Su, B.; Mao, W. Trends of streamflow in the Tarim River Basin during the past 50years: Human impact or climate change? J. Hydrol. 2011, 400, 1–9. [Google Scholar] [CrossRef]
  59. Vermote, E. MOD09A1 MODIS/Terra Surface Reflectance 8-Day L3 Global 500m SIN Grid V006. 2015, Distributed by NASA EOSDIS Land Processes DAAC. Available online: https://doi.org/10.5067/MODIS/MOD09A1.006 (accessed on 3 March 2021).
  60. Lehner, B.; Grill, G. Global river hydrography and network routing: Baseline data and new approaches to study the world’s large river systems. Hydrol. Process. 2013, 27, 2171–2186. [Google Scholar] [CrossRef]
  61. Phan, V.H.; Lindenbergh, R.C.; Menenti, M. Geometric dependency of Tibetan lakes on glacial runoff. Hydrol. Earth Syst. Sci. 2013, 17, 4061–4077. [Google Scholar] [CrossRef] [Green Version]
  62. RGI Consortium. Randolph Glacier Inventory—A Dataset of Global Glacier Outlines: Version 6.0: Technical Report; Digital Media; US National Snow and Ice Data Center: Boulder, CO, USA, 2017. [Google Scholar] [CrossRef]
  63. Guo, W.; Liu, S.; Xu, J.; Wu, L.; Shangguan, D.; Yao, X.; Wei, J.; Bao, W.; Yu, P.; Liu, Q.; et al. The second Chinese glacier inventory: Data, methods and results. J. Glaciol. 2015, 61, 357–372. [Google Scholar] [CrossRef] [Green Version]
  64. Falorni, G.; Teles, V.; Vivoni, E.; Bras, R.L.; Amaratunga, K.S. Analysis and characterization of the vertical accuracy of digital elevation models from the Shuttle Radar Topography Mission. J. Geophys. Res. Space Phys. 2005, 110, F2. [Google Scholar] [CrossRef]
  65. Zhang, J.; Jia, L.; Menenti, M.; Hu, G. Glacier Facies Mapping Using a Machine-Learning Algorithm: The Parlung Zangbo Basin Case Study. Remote Sens. 2019, 11, 452. [Google Scholar] [CrossRef] [Green Version]
  66. Liu, J.; Yao, X.; Liu, S.; Guo, W.; Xu, J. Glacial changes in the Gangdisê Mountains from 1970 to 2016. J. Geogr. Sci. 2020, 30, 131–144. [Google Scholar] [CrossRef] [Green Version]
  67. He, J.; Yang, K.; Tang, W.; Lu, H.; Qin, J.; Chen, Y.; Li, X. The first high-resolution meteorological forcing dataset for land process studies over China. Sci. Data 2020, 7, 25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Yang, K.; He, J.; Tang, W.; Qin, J.; Cheng, C.C. On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan Plateau. Agric. For. Meteorol. 2010, 150, 38–46. [Google Scholar] [CrossRef]
  69. Zhu, Z.; Woodcock, C.E. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sens. Environ. 2012, 118, 83–94. [Google Scholar] [CrossRef]
  70. Salomonson, V.; Appel, I. Development of the Aqua MODIS NDSI fractional snow cover algorithm and validation results. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1747–1756. [Google Scholar] [CrossRef]
  71. Klein, A.G.; Hall, D.K.; Riggs, G.A. Improving snow cover mapping in forests through the use of a canopy reflectance model. Hydrol. Process. 1998, 12, 1723–1744. [Google Scholar] [CrossRef]
  72. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  73. Sorg, A.; Bolch, T.; Stoffel, M.; Solomina, O.; Beniston, M. Climate change impacts on glaciers and runoff in Tien Shan (Central Asia). Nat. Clim. Chang. 2012, 2, 725–731. [Google Scholar] [CrossRef]
  74. Pan, B.T.; Zhang, G.L.; Wang, J.; Cao, B.; Geng, H.P.; Zhang, C.; Ji, Y.P. Glacier changes from 1966–2009 in the Gongga Mountains, on the south-eastern margin of the Qinghai-Tibetan Plateau and their climatic forcing. Cryosphere 2012, 6, 1087–1101. [Google Scholar] [CrossRef] [Green Version]
  75. Paul, F.; Barrand, N.; Baumann, S.; Berthier, E.; Bolch, T.; Casey, K.; Frey, H.; Joshi, S.; Konovalov, V.; Le Bris, R.; et al. On the accuracy of glacier outlines derived from remote-sensing data. Ann. Glaciol. 2013, 54, 171–182. [Google Scholar] [CrossRef] [Green Version]
  76. Hanshaw, M.N.; Bookhagen, B. Glacial areas, lake areas, and snow lines from 1975 to 2012: Status of the Cordillera Vilcanota, including the Quelccaya Ice Cap, northern central Andes, Peru. Cryosphere 2014, 8, 359–376. [Google Scholar] [CrossRef] [Green Version]
  77. Kumar, P.; Saharwardi, S.; Banerjee, A.; Azam, M.F.; Dubey, A.K.; Murtugudde, R. Snowfall Variability Dictates Glacier Mass Balance Variability in Himalaya-Karakoram. Sci. Rep. 2019, 9, 18192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Osmonov, A.; Bolch, T.; Xi, C.; Kurban, A.; Guo, W. Glacier characteristics and changes in the Sary-Jaz River Basin (Central Tien Shan, Kyrgyzstan)—1990–2010. Remote Sens. Lett. 2013, 4, 725–734. [Google Scholar] [CrossRef] [Green Version]
  79. Narama, C.; Kääb, A.; Duishonakunov, M.; Abdrakhmatov, K. Spatial variability of recent glacier area changes in the Tien Shan Mountains, Central Asia, using Corona (~1970), Landsat (~2000), and ALOS (~2007) satellite data. Glob. Planet. Chang. 2010, 71, 42–54. [Google Scholar] [CrossRef]
  80. Zhang, Z.; Xu, J.L.; Liu, S.Y.; Guo, W.Q.; Wei, J.F.; Feng, T. Glacier changes since the early 1960s, eastern Pamir, China. J. Mt. Sci. 2016, 13, 276–291. [Google Scholar] [CrossRef]
  81. Zeng, L.; Yang, T.; Tian, H. Response of glacier variations in the eastern Pamirs plateau to climate change, during the last 40 years. J. Arid. Land Resour. Environ. 2013, 27, 144–150. [Google Scholar]
  82. Feng, T.; Liu, S.Y.; Xu, J.L.; Guo, W.; Wei, J.; Zhang, Z. Glacier changes of the Yarkant River Basin from 1968 to 2009 derived from the First and Second Glacier Inventories of China. J. Glaciol. Geocryol. 2015, 37, 1–13. (In Chinese) [Google Scholar]
  83. Li, B.; Zhu, A.-X.; Zhang, Y.; Pei, T.; Qin, C.-Z.; Zhou, C. Glacier change over the past four decades in the middle Chinese Tien Shan. J. Glaciol. 2006, 52, 425–432. [Google Scholar] [CrossRef] [Green Version]
  84. Bao, W.-J.; Liu, S.-Y.; Wei, J.-F.; Guo, W.-Q. Glacier changes during the past 40 years in the West Kunlun Shan. J. Mt. Sci. 2015, 12, 344–357. [Google Scholar] [CrossRef]
  85. Li, C.X.; Yang, T.B.; Tian, H.Z. Variation of West Kunlun Mountains glacier during 1990–2011. Prog. Geogr. 2013, 32, 548–559. [Google Scholar] [CrossRef]
  86. Shangguan, D.; Liu, S.; Ding, Y.; Li, J.; Zhang, Y.; Ding, L.; Wang, X.; Xie, C.; Li, G. Glacier changes in the west Kunlun Shan from 1970 to 2001 derived from Landsat TM/ETM+ and Chinese glacier inventory data. Ann. Glaciol. 2007, 46, 204–208. [Google Scholar] [CrossRef] [Green Version]
  87. Shi, Y.; Shen, Y.; Kang, E.; Li, D.; Ding, Y.; Zhang, G.; Hu, R. Recent and Future Climate Change in Northwest China. Clim. Chang. 2006, 80, 379–393. [Google Scholar] [CrossRef]
  88. Gao, J.; Williams, M.W.; Fu, X.; Wang, G.; Gong, T. Spatiotemporal distribution of snow in eastern Tibet and the response to climate change. Remote Sens. Environ. 2012, 121, 1–9. [Google Scholar] [CrossRef]
  89. Krysanova, V.; Wortmann, M.; Bolch, T.; Merz, B.; Duethmann, D.; Walter, J.; Huang, S.; Tong, J.; Buda, S.; Kundzewicz, Z. Analysis of current trends in climate parameters, river discharge and glaciers in the Aksu River basin (Central Asia). Hydrol. Sci. J. 2015, 60, 566–590. [Google Scholar] [CrossRef] [Green Version]
  90. Shen, S.S.P.; Yao, R.; Ngo, J.; Basist, A.M.; Thomas, N.; Yao, T. Characteristics of the Tibetan Plateau snow cover variations based on daily data during 1997–2011. Theor. Appl. Clim. 2015, 120, 445–453. [Google Scholar] [CrossRef]
  91. Tahir, A.A.; Adamowski, J.F.; Chevallier, P.; Haq, A.U.; Terzago, S. Comparative assessment of spatiotemporal snow cover changes and hydrological behavior of the Gilgit, Astore and Hunza River basins (Hindukush–Karakoram–Himalaya region, Pakistan). Theor. Appl. Clim. 2016, 128, 793–811. [Google Scholar] [CrossRef]
  92. Fan, Y.; Chen, Y.; Li, W. Increasing precipitation and baseflow in Aksu River since the 1950s. Quat. Int. 2014, 336, 26–34. [Google Scholar] [CrossRef]
  93. Wang, X.; Yang, T.; Xu, C.-Y.; Yong, B.; Shi, P. Understanding the discharge regime of a glacierized alpine catchment in the Tianshan Mountains using an improved HBV-D hydrological model. Glob. Planet. Chang. 2019, 172, 211–222. [Google Scholar] [CrossRef]
  94. Xu, M.; Wu, H.; Kang, S. Impacts of climate change on the discharge and glacier mass balance of the different glacierized watersheds in the Tianshan Mountains, Central Asia. Hydrol. Process. 2018, 32, 126–145. [Google Scholar] [CrossRef]
  95. 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]
  96. Wang, P.; Li, Z.; Li, H.; Zhang, Z.; Xu, L.; Yue, X. Glaciers in Xinjiang, China: Past Changes and Current Status. Water 2020, 12, 2367. [Google Scholar] [CrossRef]
  97. Farinotti, D.; Longuevergne, L.; Moholdt, G.; Duethmann, D.; Mölg, T.; Bolch, T.; Vorogushyn, S.; Güntner, A. Substantial glacier mass loss in the Tien Shan over the past 50 years. Nat. Geosci. 2015, 8, 716–722. [Google Scholar] [CrossRef]
  98. Yi, S.; Sun, W. Evaluation of glacier changes in high-mountain Asia based on 10 year GRACE RL05 models. J. Geophys. Res. Solid Earth 2014, 119, 2504–2517. [Google Scholar] [CrossRef]
  99. Cogley, J.G.; Hock, R.; Rasmussen, L.A.; Arendt, A.A.; Bauder, A.; Braithwaite, R.J.; Jansson, P.; Kaser, G.; Möller, M.; Nicholson, L.; et al. Glossary of Glacier Mass Balance and Related Terms, IHP-VII Technical Documents in Hydrology No. 86; International Hydrological Programme: Paris, France, 2011. [Google Scholar]
Figure 1. Map of the Tarim Basin. Glacierized areas and the boundary of the Tarim Basin are indicated. Red dots are the Muztag and Kungur Tagh locations.
Figure 1. Map of the Tarim Basin. Glacierized areas and the boundary of the Tarim Basin are indicated. Red dots are the Muztag and Kungur Tagh locations.
Remotesensing 13 05117 g001
Figure 2. Work-flow implemented in the Google Earth Engine (GEE) to delineate glaciers and intermittent snow cover. NDSI, Normalized Difference Snow Index; NDWI, Normalized Difference Water Index.
Figure 2. Work-flow implemented in the Google Earth Engine (GEE) to delineate glaciers and intermittent snow cover. NDSI, Normalized Difference Snow Index; NDWI, Normalized Difference Water Index.
Remotesensing 13 05117 g002
Figure 3. Delineation of glaciers: (a) Muztag 2000; (b) Muztag 2020; (c) Kungur Tagh 2000; (d) Kungur Tagh 2020. The background is the MODIS 2020 true color composite images (Red = Band 1, Green = Band 4, Blue = Band 3). The size of one image pixel is 500 m.
Figure 3. Delineation of glaciers: (a) Muztag 2000; (b) Muztag 2020; (c) Kungur Tagh 2000; (d) Kungur Tagh 2020. The background is the MODIS 2020 true color composite images (Red = Band 1, Green = Band 4, Blue = Band 3). The size of one image pixel is 500 m.
Remotesensing 13 05117 g003
Figure 4. Annual glacier area in the Tarim Basin from 2000 to 2020 based on MODIS 500 m summertime composite images.
Figure 4. Annual glacier area in the Tarim Basin from 2000 to 2020 based on MODIS 500 m summertime composite images.
Remotesensing 13 05117 g004
Figure 5. Evolution of annual glacier area in the subregions of the Tarim Basin in 2000–2020: (a) E Tien Shan; (b) W Tien Shan; (c) Pamir; (d) W Kunlun; (e) E Kunlun. The central map shows the study area and the five subregions with the location of glaciers (light blue color areas).
Figure 5. Evolution of annual glacier area in the subregions of the Tarim Basin in 2000–2020: (a) E Tien Shan; (b) W Tien Shan; (c) Pamir; (d) W Kunlun; (e) E Kunlun. The central map shows the study area and the five subregions with the location of glaciers (light blue color areas).
Remotesensing 13 05117 g005
Figure 6. Location map and elevation of glacier area; the Digital Elevation Model is the Shuttle Radar Topography Mission (SRTM) DEM at 30 m resolution (Section 3.1.5).
Figure 6. Location map and elevation of glacier area; the Digital Elevation Model is the Shuttle Radar Topography Mission (SRTM) DEM at 30 m resolution (Section 3.1.5).
Remotesensing 13 05117 g006
Figure 7. Distribution of glacier area vs. elevation in Tarim Basin in the year of 2000, 2005, 2010, 2015 and 2020.
Figure 7. Distribution of glacier area vs. elevation in Tarim Basin in the year of 2000, 2005, 2010, 2015 and 2020.
Remotesensing 13 05117 g007
Figure 8. Eight-day interval snow cover area retrieved with MODIS data in the five subregions of the Tarim Basin: (a) E Tien Shan; (b) W Tien Shan; (c) Pamir; (d) W Kunlun; (e) E Kunlun.
Figure 8. Eight-day interval snow cover area retrieved with MODIS data in the five subregions of the Tarim Basin: (a) E Tien Shan; (b) W Tien Shan; (c) Pamir; (d) W Kunlun; (e) E Kunlun.
Remotesensing 13 05117 g008
Figure 9. Mean monthly snow cover area in the five subregions of the Tarim Basin in the period of 2000–2020.
Figure 9. Mean monthly snow cover area in the five subregions of the Tarim Basin in the period of 2000–2020.
Remotesensing 13 05117 g009
Figure 10. Remaining Snow Index in the five subregions of the Tarim Basin during 2000–2020.
Figure 10. Remaining Snow Index in the five subregions of the Tarim Basin during 2000–2020.
Remotesensing 13 05117 g010
Figure 11. Time series of mean annual surface air temperatures and mean annual precipitation (1979–2018) database. Data acquired from the China Meteorological Forcing Dataset [67].
Figure 11. Time series of mean annual surface air temperatures and mean annual precipitation (1979–2018) database. Data acquired from the China Meteorological Forcing Dataset [67].
Remotesensing 13 05117 g011
Figure 12. Example of the location of long-lasting snow in W Kunlun glaciers. The colored pixels refer to the years of persistence of long-lasting snow.
Figure 12. Example of the location of long-lasting snow in W Kunlun glaciers. The colored pixels refer to the years of persistence of long-lasting snow.
Remotesensing 13 05117 g012
Figure 13. Frequency distribution of the long-lasting snow in Tarim Basin in 2000–2020.
Figure 13. Frequency distribution of the long-lasting snow in Tarim Basin in 2000–2020.
Remotesensing 13 05117 g013
Table 1. Overview of applied satellite data on surface reflectance.
Table 1. Overview of applied satellite data on surface reflectance.
ProductsBand Information (Only Used Bands Listed)
NameSpatial Resolution (m)Wavelength (nm)Spectral Domain
MOD09A1.006 Terra Surface Reflectance 8-Day Global 500 msur_refl_b02500841–876near-IR (NIR)
sur_refl_b04500545–565Green
sur_refl_b065001628–1652SWIR2
USGS Landsat 5 TM Surface Reflectance Tier 1B230520–600Green
B430760–900near-IR (NIR)
B5301550–1750SWIR1
USGS Landsat 8 OLI Surface Reflectance Tier 1B330533–590Gree
B530851–879near-IR (NIR)
B6301566–1651SWIR1
Table 2. Estimated area of five subregions: L5-TM 2000 vs. MODIS 2000 and L8-OLI 2020 vs. MODIS 2020; relative error was calculated by applying L5-TM 2000 and L8-OLI 2020 as reference, respectively.
Table 2. Estimated area of five subregions: L5-TM 2000 vs. MODIS 2000 and L8-OLI 2020 vs. MODIS 2020; relative error was calculated by applying L5-TM 2000 and L8-OLI 2020 as reference, respectively.
NameL5-TM 2000 (km2)MODIS-2000 (km2)Relative Error of Estimate (%)L8-OLI-2020 (km2)MODIS-2020 (km2)Relative Error of Estimate (%)
East Tien Shan254.00303.0419.31179.00148.91−16.81
West Tien Shan6577.007477.1613.696673.546503.41−2.55
Pamir6846.216383.51−6.766403.005073.40−20.77
West Kunlun17,436.7616,052.64−7.9412,631.3811,079.11−12.29
East Kunlun1935.021546.00−20.101046.66981.83−6.19
Table 3. Comparison of glacier area from various studies surrounding the Tarim Basin. The subregions correspond to Figure 6.
Table 3. Comparison of glacier area from various studies surrounding the Tarim Basin. The subregions correspond to Figure 6.
SubregionGlacier Area Change
PeriodSpecific RegionRate(%/a)Reference
East Tien Shan1990–2000 −0.8[49]
2000–2020 −2.98This study
West Tien Shan1970–2000Pskem region−0.63[79]
2000–2007 −0.71[79]
1975–2008Aksu Catchment−0.11[6]
1956–2007Borkoldoy−0.57[4]
1956–2007Naryntau−1.00[4]
1975–2016Aksu River Basin−0.63[19]
1986–2000Kaidu River Basin−0.35[83]
1990–2010Sary-Jaz River Basin−3.7[78]
2000–2010 −0.51[79]
2000–2010 −0.22[49]
2000–2010 −0.73This study
1990–2018 −0.44[49]
2000–2020 −1.07This study
Pamir2000–2011sEast Pamirs−0.05[81]
1963–2009Pamirs in China−0.24[80]
1968–2009Yarkant River Basin−0.36[82]
2000–2010 −0.39[49]
2000–2010 −0.46This study
1990–2018 −0.59[49]
2000–2020 −0.50This study
West Kunlun1970–2010 −0.09[84]
1999–2011 −0.07[85]
2000–2010 −2.8This study
1990–2000 0.12[49]
1990–2000 0.01[86]
2000–2020 −1.04This study
East Kunlun1990–2018 −0.23[49]
2000–2020 −0.81This study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Jia, L.; Menenti, M.; Zhou, J.; Ren, S. Glacier Area and Snow Cover Changes in the Range System Surrounding Tarim from 2000 to 2020 Using Google Earth Engine. Remote Sens. 2021, 13, 5117. https://doi.org/10.3390/rs13245117

AMA Style

Zhang J, Jia L, Menenti M, Zhou J, Ren S. Glacier Area and Snow Cover Changes in the Range System Surrounding Tarim from 2000 to 2020 Using Google Earth Engine. Remote Sensing. 2021; 13(24):5117. https://doi.org/10.3390/rs13245117

Chicago/Turabian Style

Zhang, Jing, Li Jia, Massimo Menenti, Jie Zhou, and Shaoting Ren. 2021. "Glacier Area and Snow Cover Changes in the Range System Surrounding Tarim from 2000 to 2020 Using Google Earth Engine" Remote Sensing 13, no. 24: 5117. https://doi.org/10.3390/rs13245117

APA Style

Zhang, J., Jia, L., Menenti, M., Zhou, J., & Ren, S. (2021). Glacier Area and Snow Cover Changes in the Range System Surrounding Tarim from 2000 to 2020 Using Google Earth Engine. Remote Sensing, 13(24), 5117. https://doi.org/10.3390/rs13245117

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