Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Impact of Air Pollution (PM2.5) on Child Mortality: Evidence from Sixteen Asian Countries
Previous Article in Journal
Multi-Dimensional Accessibility Barriers in Care Services for the Rural Elderly with Disabilities: A Qualitative Study in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monitoring the Spatiotemporal Difference in Glacier Elevation on Bogda Mountain from 2000 to 2017

1
School of Surveying and Land Information Engineering, Henan Polytechnic University, Jiaozuo 454003, China
2
State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China
3
Yellow River Institute of Hydraulic Research, Yellow River Conservancy Commission of the Ministry of Water Resources, Zhengzhou 450003, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(12), 6374; https://doi.org/10.3390/ijerph18126374
Submission received: 22 April 2021 / Revised: 4 June 2021 / Accepted: 10 June 2021 / Published: 12 June 2021
(This article belongs to the Section Environmental Earth Science and Medical Geology)

Abstract

:
The difference in glacier surface elevation is a sensitive indicator of climate change and is also important for disaster warning and water supply. In this paper, 25 glaciers on Bogda Mountain, in the eastern Tianshan Mountains, are selected as the study object as they are typical of glaciers in arid or semi-arid areas with importance for water supply. The Repeat Orbit Interferometry (ROI) method is used to survey the surface elevation of these glaciers using Sentinel-1A Radar data from 2017. Using data from the Shuttle Radar Topography Mission (SRTM) and a Digital Elevation Model (DEM), the difference in the glacier surface elevation between 2000 and 2017 is obtained. A scheme to evaluate the accuracy of estimated variations in glacier surface elevation is proposed in this article. By considering the surfaces of lakes in the study region as ideal horizontal planes, the average standard deviation (SD) value of the lake elevation is taken as the error caused by the radar sensor and observing conditions. The SD of the lake elevation is used as an index to evaluate the error in the estimated variation of the glacier surface elevation, and the obtained SD values indicate that the result obtained using the ROI method is reliable. Additionally, the glacier surface elevation variation pattern and a Logarithmic Fitting Model (LFM) are used to reduce the error in high-altitude glacial accumulation areas to improve the estimation of the difference in the glacier surface elevation obtained using ROI. The average SD of the elevation of the 12 lakes is ±2.87 m, which shows that the obtained glacier surface elevations are reliable. This article concludes that, between 2000 and 2017, the surface elevation of glaciers on Bogda Mountain decreased by an average of 11.6 ± 1.3 m, corresponding to an average decrease rate of 0.68 m/a, and glaciers volume decreased by an average of 0.504 km3. Meanwhile, the surface elevations of the lakes increased by an average of 8.16 m. The decrease of glacier surface elevation leads to the expansion of glacial lakes. From the north slope clockwise to the south slope, the glacier elevation variation showed a decreasing trend, and the elevation variation gradually increased from the south slope to the north slope. With the increase of glacier altitude, the variation of glacier surface elevation gradually changed from negative to positive. The findings of this article suggest that the rate of glacier retreat on Bogda Mountain increased from 2000 to 2017.

1. Introduction

Glacier variation is a sensitive indicator of climate change [1] and is also important for disaster warning and water supply in arid and semi-arid areas [2]. Bogda glacier is a typical glacier with important water supply significance in arid and semi-arid areas [3]. The melting of glaciers in the Tianshan area has an important influence on the freshwater resources and economic development in Xinjiang and the whole of China. The glaciers in the east Tianshan mountains are very sensitive to climate change in Western China and their retreat and response to climate play an important indicating role. The monitoring of glaciers can identify their change characteristics, including their response to climate change. Moreover, the change of glacier melting and accumulation have an important impact on rivers and therefore also on ecosystems and economic production [4]. Thus, the advance and retreat of glaciers is of great significance to human activities and economic construction [5,6]. Additionally, the rapid melting of glaciers makes some lakes unstable, which leads to a series of ecological and environmental problems such as ice lake breakwaters, glacial floods, and changes in the coastal ecological environment, thereby affecting people’s health. Therefore, glacier monitoring is of great significance [7]. Glacier monitoring includes glacier area monitoring, height monitoring, volume monitoring, and mass balance monitoring. Among these, the monitoring of glacier elevation can reflect the volume change and elevation change of glaciers [8] and also provides support for the study of glacier mass balance (GMB).
The routine monitoring of changes in glacier elevation mainly involves field stereophotogrammetry, the Differential Global Positioning System technique, and ground-penetrating radar analysis [9,10]. However, although these methods have high precision, they can only be applied over a limited area [11]. Furthermore, the harsh environmental conditions in high-altitude areas can make the onsite measurement of glaciers dangerous, which represents another limitation of such measurements [12].
In recent years, with the rapid development of remote sensing technology, an increasing number of satellites that can be used for glacier elevation monitoring have become available and new elevation monitoring methods have been developed. Stereo mapping using visible-wavelength images, Interferometric Synthetic-aperture Radar (InSAR), satellite gravimetry, and laser altimetry are the main methods for glacier elevation monitoring. Moholdt et al. obtained recent and long-term elevation changes of Svalbard glacier based on laser altimetry measurements from ICESat [13,14]. Farintotti et al. used three independent methods based on satellite gravity measurements, laser altimetry, and glaciological modeling to estimate the change of the area of all glaciers in the Tianshan Mountains [15]. Furthermore, ICESat laser altimetry data and a Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) were used to study the GMB in the Qilian Mountains, the Qinghai-Tibet Plateau, the Himalayans, and the Tianshan Mountains [16,17]. Some studies using the Repeat Orbit Interferometry (ROI) method estimated changes in GMB in the Tianshan Mountains using TanDEM-X SAR data to obtain a high-accuracy DEM [18]. These applications have shown that the ROI method using SAR data can provide elevation variation information with a high spatiotemporal resolution, thus providing comprehensive information about glacier morphological changes [19]. However, remote sensing methods have large uncertainties, and much work is needed to remove errors associated with remote sensing observations. In order to obtain reliable results for glacier elevation change, Farinotti [20] et al. used a simple dynamic model to obtain the spatially distributed thickness of individual glaciers by inverting their surface topography. The results were validated against a comprehensive set of thickness observations for 300 glaciers from most of the glaciated regions of the world. Huang et al. [21] proposed a new method based on the polynomial fitting of the GLAS footprint of the Sparse ICESat—GLAS to estimate the elevation change of the surface of two glaciers in Naimona’nyi glacier and Yanong glacier; the results showed that the polynomial fitting method gives a reliable estimation of elevation change. These methods offer new methods to monitor the elevation of glaciers on Bogda Mountain.
Researchers have invested large amounts of financial and human resources in monitoring the elevations of glaciers on Bogda Mountain [22]. For example, in the early 1980s, a joint Chinese–Japanese team estimated the accumulation, ablation characteristics, and material balance of the Sigonghe glacier No. 4. Additionally, researchers at the Tianshan Glacier Station conducted an extensive scientific investigation of a glacier on Bogda Mountain [23,24]. In this paper, Sentinel-1A radar data and Landsat images are used to obtain the surface elevations of glaciers on Bogda Mountain using the ROI method; the terrain law analysis of glacial central line was used to improve the accuracy of InSAR observations; the standard deviation (SD) of glacial lake elevation was used to evaluate the InSAR observational errors; the difference in the surface elevations of glaciers and glacial lakes was determined; and the response of glaciers to climate change was discussed.

2. Data and Methods

2.1. Study Area

Bogda Mountain, the highest peak in the eastern Tianshan Mountains, is located at 43.8° N, 88.3° E, around 50 km from Urumqi City, the capital of Xinjiang Province. The mountain has an average altitude of 5445 m. The location of the study area is shown in Figure 1. Prior to 1983, there were 113 glaciers on Bogda Mountain. However, this number has now decreased to 110 [25]. Many of these glaciers are of the valley type, such as the Heigou glacier No. 8 and the Sigonghe glacier No. 4, which are the sources of many rivers [26]. The Bogda Mountain area is cold and dry in winter and is influenced by westerly circulation in summer, which gives the area a continental climate typical of the mid-latitude westerlies. The annual and daily temperature varies greatly. In summer, precipitation is low, evaporation is high, and the variation of river runoff depends on glacial melting to a great extent [27].
As shown in Figure 1c, 25 typical glaciers on Bogda Mountain were investigated in this study. These include Heigou glacier No. 8 (which in this article is denoted Glacier No. 17); Sigonghe glacier No. 4, (here, Glacier No. 11); and fan-shaped distributary glaciers (here, glacier nos. 13 and 14). The elevations of these 25 glaciers were extracted using 25 glacial vector boundaries by using the vector file of CGI2 (the Second China Glacier Inventory) as a mask file. Additionally, as shown in Figure 1c, 12 lakes were also selected as study objects. Except for Lake 7, which is far from the glacier area, the rest of the lakes are within 100 m of a glacier and are closely related to the variation of glacier elevation. More information about the glaciers is shown in Table 1.

2.2. Data

The Sentinel-1A data were downloaded from the website of the European Space Agency (https://scihub.copernicus.eu/dhus, last access: 10 June 2021). To ensure the satisfactory coherence of the complex image pair, two complex SAR image pairs with intervals of 12 days were selected, as well as precise orbit data corresponding to SAR images [28]. Sentinel-1A is in a sun-synchronous orbit with an orbit height of 693 km and an inclination angle of 98.15°, and the introduction of precise orbit data can remove the systematic errors caused by orbit errors [29]. Table 2 shows the relevant information of the SAR image pairs that were used in this study to generate the DEM using the ROI method.
The elevation datum version 4 of the SRTM 3 DEM data is the geoid of EGM96 (Earth Gravitational Model 1996), and the plane datum is WGS84 (World Geodetic System 1984). The SRTM is a global DEM obtained using InSAR and has a vertical accuracy of ±16 m [30]. Additionally, the locations, names, sizes, and other general information of the glaciers on Bogda Mountain were obtained from CGI2. The CGI2 glacier information is based on remote sensing images and is compiled using a glacial delineation method and GIS technology [31]. A Landsat-8 image from 29 September 2017 and a Landsat-5 image from 2 September 2000 were selected to delineate lake boundaries and facilitate the analysis of lake elevation obtained from ROI.
This paper used temperature and precipitation data from 2000 to 2017 recorded at three meteorological stations near the study area, namely Dabancheng (Station No. 51477; 43.21° N, 88.19° E; altitude 1103.5 m), Kumish (Station No. 51526; 42.14° N, 88.13° E; altitude 922.4 m), and Turpan (Station No. 51573; 42.56° N, 89.12° E; altitude 34.5 m); these data were obtained from the China Meteorological Data Network (http://data.cma.cn/user/info.html, last access: 10 June 2021).
Synthetic-aperture radar (SAR) data acquired in summer in vertical–vertical polarization mode were selected. Summer images were chosen to avoid snow cover in non-glaciated areas, which would lead to serious decoherence in these areas. A SAR image pair for 2018 (see Table 2) was compared with a SAR image pair for 2017 (see Table 2). It was thus found that the accuracy of the 2017 pair was better than that of the 2018 pair; therefore, the SAR image pair for 2017 was chosen as a reference to study the difference in glacier elevation between 2000 and 2017.
The central flow of data processing and results analysis is shown in Figure 2.

2.3. Methods

2.3.1. Using the ROI Method to Generate a DEM

The ROI method was employed to compute the variation in glacier surface elevation using the SARscape software (Sarmap, Caslano, Switzerland). This method includes baseline estimation, interferogram generation, filtering, coherence calculation, phase unwrapping, track refinement, repeated removal of the phase of the flat earth, phase-to-elevation, and geocoding [32]. Following the above processes, using the DEM generated with the ROI method, the elevation of each pixel is calculated. Thus, the elevation of the glaciers in the study area is obtained. The glacier vector boundaries of CGI2 and glacier boundaries delineated from Landsat OLI images were used to clip the DEM image so that only glaciers and lakes were obtained [33]. Finally, the elevations of the 25 glaciers obtained from Sentinel-1A data were compared with the glacier elevations obtained from the SRTM DEM to obtain the glacier surface variation between 2000 and 2017.
This study used an accuracy index based on the principle of interference theory. The ROI process is affected by many factors, such as snow scattering characteristics, decoherence effects, and radar imaging characteristics. It is difficult to quantify the influence of each factor separately. Therefore, in this study, a formula for the calculation of elevation accuracy [34] based on a large amount of practical modeling experience was used.
p r e = 1 γ 2 2 γ 2 λ R sin θ 4 π B
where γ is the coherence coefficient, λ is the wavelength, R is the slant distance, θ is the local incidence angle, and B is the baseline length of the two radar images. However, further analysis is needed to assess the accuracy of the calculated elevation variation; this is because the backscatter of the radar signal of water is weak, which means that, in glacial and lake areas, the DEM accuracy is relatively low. In this study, in Section 2.3.2, a Logarithmic Fitting Model (LFM) was introduced to improve the accuracy of the DEM in glacial and lake areas.
To evaluate the calculated variation of glacier elevation, the SD of lake elevation was used as an index to measure the error of the elevation change of nearby glaciers. Each lake was extracted from Landsat 8 OLI images using the Normalized Difference Water Index (NDWI) [35]. Because the lake surface is a plane, the theoretical lake surface elevation is a unified value. However, due to the measurement error and the scattering characteristics of the lake itself, the lake surface elevation varies slightly. Assuming that the lake surface is an ideal horizontal plane [36], the SD value shown in Formula (2) was used as the error due to the radar sensor and observing environment.
S D = ± 1 N i = 1 N ( x i u )
where N is the number of pixels occupied by the lake, x i is the elevation value of pixel i , and u is the average of the elevation values of all the lake pixels. It was found that the average SD of the elevations of the 12 lakes was ±2.87 m, which is far smaller than the calculated variation of glacial elevation, suggesting that the obtained difference in glacier surface elevation is reliable.

2.3.2. Establishment of Logarithmic Fitting Model (LFM)

Here, the pixels whose glacier elevation variation values are within 1% of the actual variation are selected as the original sampling pixels. In Figure 3, Glacier No. 14 is selected as an example to illustrate the selection of the original sampling pixels. The selection rules were as follows: (1) avoid speckle areas (as can be seen in Figure 3); (2) select pixels located on a curve perpendicular to the isoline (as shown by the black sampling line in the enlarged view of Figure 3a); (3) select pixels where the elevation variation value decreases as the elevation increases. This was performed to make the original sampling pixels meet the prerequisite that the glacial accumulation area be located at a higher altitude and the glacial melting area be located at a lower altitude; and (4) select as many pixels as possible. The more pixels, the better the LFM analysis can reflect the real glacial elevation variation. The improved Glacier No. 14 elevation variation map after applying the LFM is shown in Figure 3b.
As shown in Figure 4, using the glacier elevation in 2000 as a reference point, the glacier elevation variation rule from 2000 to 2017 was determined by producing a scatter plot between the SRTM and the mean elevation variation, and subsequently, the LFM curve, which can reflect the characteristics of elevation variation, was fitted using these elevation variation points. Taking the value of the SRTM DEM as the independent variable and the elevation variation value of the corresponding pixel as the dependent variable, the LFM function reflects the change relationship between these two values. Finally, the SRTM DEM and the LFM function were used to calculate the glacier elevation change.
The glacier surface elevation variation [37] pattern is described as follows. For a glacier, normally, the GMB line can be defined as a contour line whose altitude is the average altitude of the whole glacier. The GMB line divides the glacier into two parts. Part 1 is the upper part, which represents the accumulation area with higher altitude; in this part, the glacier surface elevation increases because of the increased snowpack depth. If the glacier experiences accumulation, the maximum glacier surface elevation variation pixel is located in the accumulation part, and, as elevation decreases, the absolute elevation variation value reduces. Part 2 is the lower part, which represents the ablation area with lower altitude; in this part, the glacier surface elevation decreases because of the higher temperature. If the glacier retreats, the maximum glacier surface elevation variation pixel is located in the ablation part, and, as elevation increases, the absolute elevation variation value reduces. As shown in Figure 5, in the study area, between 2000 and 2017, the glacier variation trend was dominated by retreat. In the LFM, the dependent variable decreases as the independent variable increases, and the LFM was used to fit the elevation value and the elevation variation value. Here, the selected pixel takes the SRTM DEM value as the abscissa axis, and the elevation variation value of the corresponding position is the ordinate axis. Using the MATLAB software (MathWorks, Natick, MA, USA), the coordinates of the pixels were read and a scatter diagram was produced, as shown in Figure 4. From the scatter plot, it can be seen that the higher the elevation value, the higher the elevation variation value is and the larger the error trend. The LFM can reduce the elevation variation error of the glacier accumulation area, and the LFM result follows the glacier surface elevation variation rule. The MATLAB Curve Fitting Toolbox was used for fitting analysis, and the LFM was obtained using Formula (3).
y = a * ln x + c
where y is the elevation variation value from 2000 to 2017, x is the elevation value for 2000 obtained from the SRTM DEM, and a and c are the fitting coefficients, with values of 117.3 and −982.5, respectively. Using the LFM, the coefficient of determination and the SD were found to be 0.907 and ±3.713 m, respectively. Because the coefficient of determination is close to 1 and the value of SD is not large, the LFM result is acceptable. The curve fitting is also shown in Figure 4.

3. Results

3.1. Threshold Analysis of Glacier Elevation Variation

Due to the weak backscattering characteristics of water, the accuracy of DEMs produced using the ROI method is lower for glacial and lake areas. The error can reach hundreds of meters due to the complex radar scattering mechanism and environmental factors [38]. Therefore, in this study, a threshold of glacier elevation variation was determined and was taken as the maximum of the glacier elevation variation. This threshold was used as a reference to reduce the error. The results of previous studies of the elevation variation of glaciers on Bogda Mountain—which are summarized later in this paragraph—were used as a reference for estimating the maximum threshold for glacier variation. For example, in [39], researchers showed that the elevation of Heigou glacier No. 8 (Glacier No. 17 in Figure 6) decreased by 13 m between 1986 and 2009 and that the elevation of Sigonghe glacier No. 4 (Glacier No. 11 in Figure 6) decreased by 15 m between 1962 and 2009. The geographic locations of two fan-shaped distributary glaciers (glaciers No. 13 and No. 14 in Figure 6) and Sigonghe glacier No. 4 were similar, and their water and heat conditions were also similar; therefore, the elevations of the two glaciers decreased at similar rates from 1962 to 2006, reducing by nearly 15 m.
Considering that the melting of glaciers may have accelerated in recent years, the threshold value of elevation variation was set to 1.5 times the experientially derived glacier elevation change introduced above, that is, the threshold is defined as ±25 m. Therefore, the elevation pixels whose absolute value of elevation variation is more than 25 m were removed, and the Kriging interpolation method was used to interpolate for the null-value pixels after this initial noise removal. The results of the estimated glacier surface elevation variation are shown in Figure 5. As shown in the figure, most of the glacial areas show elevation reduction in the ablation area with lower altitude under the GMB line and a slight elevation increase in the accumulation area with higher altitude above the GMB line. These observations are consistent with glacier dynamics theory [40]. However, the glacier surface elevation variation obtained in this way did not conform well to the variation obtained from CGI2. Therefore, the following improvements were made to improve the accuracy of the estimated glacier surface elevation variation.

3.2. LFM Results

Using the LFM and the elevation map for 2000 as the benchmark, the improved elevation variation results for the 25 glaciers and 12 lakes were obtained, as shown in Figure 6. From the above research, the following results can be obtained. Between 2000 and 2017, the glaciers on Bogda Mountain exhibited a retreating trend, and the glacier elevation decreased. During this period, the glacier volume decreased by an average of 0.504 km3, and the glacier elevation decreased by an average of 11.6 ± 1.3 m. In the following, the results obtained in this article are compared with previous results in the literature to prove the correctness of the proposed method.
(1) Glacier No. 17. The results of this study suggest that, from 2000 to 2017, the glacier’s volume decreased by 0.024 km3, and its average elevation decreased by 13 ± 0.9 m, representing an average decrease of 0.76 ± 0.05 m/a. This conclusion is consistent with the conclusions of a previous study [41], which found that, between 1986 and 2009, the volume of Glacier No. 17 decreased by 0.026 ± 0.0118 km3, the average thickness of the ice tongue decreased by 13 ± 6 m, and the thickness decreased by an average of 0.57 ± 0.26 m/a. Taking these results together with those obtained in the present study indicates that the thinning rate of Glacier No. 17 increased from 2000 to 2017, which indicates that the glacier retreat is becoming increasingly serious. (2) Glacier No. 11. The results of the present study suggest that, from 2000 to 2017, its volume decreased by 0.076 km3, and its average elevation decreased by 18 ± 1.9 m, representing an average decrease of 1.06 ± 0.11 m/a. This conclusion is consistent with the finding [42] that, between 1962 and 2009, the ice reserves of Glacier No. 11 decreased by 0.014 ± 0.0095 km3, the average thickness of the ice tongue decreased by 15 ± 10 m, and the thickness of the ice tongue decreased by an average of 0.32 ± 0.2 m/a. By comparing these results with those obtained in the present study, it can be seen that the thinning rate of Glacier No. 11 increased from 2000 to 2017 which indicates that the glacier is retreating faster. (3) Glacier No. 14. The results of this study indicate that, from 2000 to 2017, its volume decreased by 0.18 km3 and its average elevation decreased by 14 ± 1.1 m. Because the location of this glacier is similar to that of Glacier No. 11, and the water and heat conditions are similar, it has similar characteristics of elevation change. Therefore, from the above analysis, it can be seen that the estimated elevation variation of Glacier No. 14 is reliable. By comparing these results with the results of the present study, it can be found that the thinning rate of Glacier No. 14 increased from 2000 to 2017, which indicates that the glacier is retreating faster. The above comparison with previous literature results suggests that the results obtained in this paper are credible. The results of the elevation variations of other glaciers determined in the present study are shown in Figure 7. Figure 7a shows the maximum and minimum elevation variation of the 25 glaciers during the period between 2000 and 2017, and Figure 7b shows average elevation variation and SD.

3.3. Meteorological Data Processing Results

Using the annual precipitation anomaly and the annual average temperature change measured at the Dabancheng, Kumish, and Turpan meteorological stations in the area surrounding Bogda Mountain, a map showing the annual precipitation anomaly and annual average temperature change for 2000 to 2018 in the Bogda Mountain area was produced, as shown in Figure 8. From this figure, the following information can be obtained. From 2000 to 2018, the average temperature of the three meteorological stations was 11.07 °C and the average temperature increase rate was 0.42 °C·(10a)−1. The average precipitation was 50.6 mm and the decrease rate of precipitation was 8.75 mm·(10a)−1. Zhang Mingjun et al. found that for every 1 °C temperature rise, corresponding precipitation needs to increase by 47% or 56% to compensate for the melting of glaciers [43]. Therefore, the increase of precipitation cannot make up for the ablation caused by temperature rise. Using wavelet analysis, He et al. [44] showed that the Bogda Mountain area was in a period of high temperature and low precipitation from 1972 to 2013. Therefore, it can be inferred that temperature has a significant effect on glacier variation in the Bogda Mountain area.

4. Discussion

The variation of glacier surface elevation is related to many factors, such as glacier surface area, aspect, altitude, temperature, and precipitation. Therefore, based on the average elevation variation, this paper analyzes the relationship between the surface elevation change of glaciers in the Bogda Mountain area and the aforementioned five quantities.

4.1. Effects of Glacier Area on Glacier Elevation

By combining the average glacier elevation variation with the information of glacier area and altitude in CGI2, a relationship diagram was established between glacier area, glacier altitude, and average glacier elevation variation from 2000 to 2017, as shown in Figure 9. From the figure, it can be seen that most of the glaciers with large elevation variation are located at lower altitudes, while the glaciers with small elevation variation are located at higher altitudes. For glaciers of a similar size, glaciers with a lower altitude experienced more variation than glaciers with a higher altitude; that is, there is a relationship between glacier area variation and elevation. The results of this analysis suggest that, in general, the size of glaciers affects the overall elevation difference.

4.2. Effects of Glacier Aspect on Glacier Elevation

By combining the average glacier elevation variation with the information of glacier aspect in CGI2, the relationship between the average glacier elevation variation and aspect from 2000 to 2017 was established. The results are shown in Figure 10. The aspect is divided into eight areas, of which (0–22.5°, 337.5–360°) is north and (22.5–67.5°) is northeast. According to this classification rule, the clockwise direction is east, southeast, south, southwest, west, and northwest. According to the information shown in Figure 10, overall, the variation of the north slope aspect is relatively fast, and the variation of the south slope aspect is relatively slow. From the north slope clockwise to the south slope, the range of variation shows a decreasing trend and gradually expands from the south slope to the north slope. The rate of change is highest in the north direction and lowest in the south direction. This trend is consistent with the conclusions of Zhou Yuangang [45]. This phenomenon may be caused by the temperature difference caused by the influence of the aspect on the direction of water vapor movement and sunlight irradiance [46]. Additionally, because the studied glaciers are located at different altitudes, elevation is also an important factor causing the glacier surface elevation variation. In this paper, the relationship between glacier aspect and altitude from 2000 to 2017 was established, as shown in Figure 10 and Figure 11. It can be appreciated that the glaciers with fast ablation are mostly located at low altitudes, while glaciers with slow ablation or even accumulation are mostly in high-altitude areas. This proves that the influence of altitude on glacier elevation is greater than that of aspect on glacier variation.

4.3. Effects of Glacier Altitude on Glacier Elevation

By combining the average glacier elevation variation with the information of glacier altitude in CGI2, we established the relationship between the average elevation variation of glaciers and the altitude between 2000 and 2017, as shown in Figure 12. As shown in the figure, the elevation variations of the glaciers range from −25 to 20 m. Because the temperature rises faster in low-altitude areas than in high-altitude areas, glaciers at low altitudes experience more melting, so there is a decrease in the ablation area. With the increase of altitude, affected by snowfall, glaciers accumulate, thus the glacier elevation tends to increase in parts of the glaciers at higher altitude. With the increase of altitude, the glacier elevation variation gradually changes from negative to positive. The trend of glacier elevation difference is consistent with the glacier surface elevation variation pattern.

4.4. Effect of Lakes on Glacier Elevation

After a long period of melting, a lake will appear in the area before the glacier terminus. As shown in Figure 6, 11 lakes within 100 m of glaciers and one independent lake were selected. The basic information of these lakes is shown in Table 3. From Figure 6 and Table 3, it can be seen that Glacier No. 13 and Glacier No. 14 are closely related to Lake No. 1. Lake No. 7 is an independent natural lake, and the other lakes have a glacier meltwater supply. As shown in Table 3, the SD of the lake elevation variation was used as an index to evaluate the reliability of the glacier elevation variation. As shown in the table, the absolute SD values are all smaller than the lake elevation variation values and the glacier elevation variation values; this shows that the deformation information extracted using the ROI method is reliable. It was found that, for the lakes, there is a relationship between the lake area and the average glacier elevation variation. Between 2000 and 2017, the lake area expanded, the lake elevation increased, and correspondingly the glacier surface elevation decreased. This is due to the influence of glacial melting, which caused the glacial lake to expand and the lake’s average elevation to increase.

4.5. Effects of Meteorological Factors on Glacier Elevation

Glacier change is caused by climate change. Temperature and precipitation affect the development of glaciers and are the main climatic factors affecting glacier change; temperature determines the amount of melting and precipitation determines the amount of recharge [47,48,49]. Zhang Qibing et al. studied the relationship between glacier change and meteorological factors in the Qilian Mountains, an arid area in Western China. They found that glacier shrinkage is related to temperature rise [50]. In the present paper, through the processing of data from three meteorological stations around the Bogda Mountain area, we reached the conclusion that temperature has a significant effect on glacier variation in this area. The reliability of the results is indicated by the fact that Zhou et al. obtained the same results in their analysis of the influencing factors of glaciers in the Bogda peak area; that is, they found that the temperature has a significant effect on the change of the glacier area and snow cover area in this area [45]. Therefore, similar topography and environment may cause the same changes in glaciers.

5. Conclusions

In this paper, the surface elevation of Bogda Mountain in 2017 was obtained using orbit interferometry. Then, by comparison with the SRTM DEM of the mountain for 2000, the difference in the surface elevations of 25 glaciers on the mountain between 2000 and 2017 was obtained. It was found that the elevations of the glaciers decreased, and the glacier reserves decreased. From 2000 to 2017, the glacier volume decreased by an average of 0.504 km3, the elevation decreased by an average of 11.6 ± 1.3 m, and the average glacier elevation decrease rate was 0.68 m/a. As the elevations of the GMB lines of the glaciers rose, the glaciers tended to retreat faster. From 2000 to 2017, for the typical glaciers numbered Glacier No. 17, No. 11, and No. 14, the glacier volumes decreased by 0.024 km3, 0.076 m3, and 0.18 km3, the average glacier elevation decreased by 13 ± 0.9 m, 18 ± 1.9 m, and 14 ± 1.1 m, respectively, and the glaciers’ annual average elevation decreases were 0.76 ± 0.05 m/a, 1.06 ± 0.11 m/a, and 0.82 ± 0.064 m/a, respectively. These observations for these three glaciers are consistent with the conclusions of a previous study [41]. The comparison between the results presented in this study and previous results indicates that the glacial thinning rate increased from 2000 to 2017, which indicates that glacial retreat is becoming increasingly serious.
Finally, we analyzed the causes of the surface elevation variations of glaciers on Bogda Mountain from the aspects of glacier surface area, aspect, altitude, pre-glacial lake, temperature, and precipitation. The variations of glacier area, lake, and glacier elevation were small. From the north slope clockwise to the south slope, the glacier elevation variation showed a decreasing trend, and the elevation variation gradually increased from the south slope to the north slope. With the increase of glacier altitude, the variation of glacier surface elevation gradually changed from negative to positive. However, temperature has a significant effect on glaciers in the Bogda Mountain area.
The variation of glacier elevation is caused by many factors. Therefore, further studies are necessary to obtain more accurate estimates of glacier elevation changes.

Author Contributions

W.D. designed the research. N.S. processed the data, analysed the results, and wrote the manuscript. L.X., S.Z., D.M. and S.W. provided analysis tools and technical assistance. All authors contributed to the final version of the manuscript by proofreading and providing constructive ideas. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly funded by the Special Funds of the National Natural Science Foundation of China, grant numbers 42041004 and 42041006, the National Natural Science Foundation of China, grant numbers 41601364, 51709123, 41601392 and 41975036, the Scientific and Technological Projects of Henan Province, grant number 172102210280, the Doctoral Foundation of Henan Polytechnic University, grant number B2016-11, and the Special Basic Research Fund for Central Public Research Institutes, grant number HKY-JBYW-2018-03. The APC was funded by Henan Polytechnic University, grant number 41975036.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express sincere gratitude to the European Space Agency for the provision of Sentinel-1A data from their website at https://scihub.copernicus.eu/dhus, last access: 10 June 2021. We acknowledge the U.S. Geological Survey for sharing the SRTM DEM and Landsat image data from their website at https://earthexplorer.usgs.gov, last access: 10 June 2021. We thank the China Cold and Arid Regions Scientific Data Center for sharing the Second Glacier Inventory Dataset of China Version 1.0 (CGI2). We also thank Tongtong Zhu from the China University of Geosciences for the data collection and preprocessing work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, W.B.; Li, J.L.; Bao, A.M.; Wang, B.S. Mapping changes in the glaciers of the eastern Tienshan Mountains during 1977–2013 using multitemporal remote sensing. J. Appl. Remote Sens. 2014, 8, 084683. [Google Scholar] [CrossRef] [Green Version]
  2. Rossman, N.R.; Zlotnik, V.A.; Rowe, C.M.; Szilagyi, J. Vadose zone lag time and potential 21st century climate change effects on spatially distributed groundwater recharge in the semi-arid Nebraska Sand Hills. J. Hydrol. 2014, 519, 656–669. [Google Scholar] [CrossRef] [Green Version]
  3. Yang, X.H.; Zhao, J.D.; Han, H. Mass balance of glaciers in the Harlik mountain area of Eastern Tianshan Mountains from 1972 to 2016. J. Glaciol. Geocryol. 2019, 41, 1–11. [Google Scholar]
  4. Hoelzle, M.; Azisov, E.; Barandun, M.; Huss, M.; Zemp, M. Re-establishing glacier monitoring in Kyrgyzstan and Uzbekistan, Central Asia. Geosci. Instrum. Methods Data Syst. Discuss. 2017, 6, 397–418. [Google Scholar] [CrossRef] [Green Version]
  5. Qin, Y.; Yang, T.B.; Ji, Q.; Wang, J. Analysis of modern glacier change in Harlik mountain from 1973 to 2016. J. LanZhou Univ. 2018, 54, 569–576. [Google Scholar]
  6. Wang, N.L.; Yao, T.D.; Xu, B.Q.; Chen, A.A.; Wang, W.C. Temporal and spatial patterns, trends and impacts of glacier changes in the Tibetan Plateau and its surrounding areas under global warming. Chin. Acad. Sci. 2019, 64, 2865–2875. [Google Scholar]
  7. Kulshrestha, A.; Bahuguna, I.M.; Rathore, B.P.; Iyer, K.V. Resolving Biases in DEM Differencing for Estimation of Change in Elevation of Glacier Surfaces Using Cartosat-I Stereo Data. J. Indian Soc. Remote Sens. 2020, 48, 1443–1453. [Google Scholar] [CrossRef]
  8. Patel, A.; Prajapati, R.; Dharpure, J.K.; Mani, S.; Chauhan, D. Mapping and monitoring of glacier areal changes using multispectral and elevation data: A case study over Chhota-Shigri glacier. Earth. Sci. Inf. 2019, 12, 489–499. [Google Scholar] [CrossRef]
  9. Wang, Y.; Hou, S.; Liu, Y. Glacier changes in the Karlik Shan, eastern Tien Shan, during 1971/72–2001/02. Ann. Glaciol. 2009, 50, 39–45. [Google Scholar] [CrossRef] [Green Version]
  10. Du, W.; Ji, W.; Xu, L.; Wang, S. Deformation time series and driving-force analysis of glaciers in the Eastern Tienshan Mountains using the SBAS InSAR method. Int. J. Environ. Res. Public Health 2020, 17, 2836. [Google Scholar] [CrossRef] [Green Version]
  11. Minallah, N.; Khan, W. Comparison of neural networks and support vector machines for the mass balance ablation observation of glaciers in Baltoro region. J. Inf. Commun. Technol. Robot. Appl. 2019, 9, 37–45. [Google Scholar]
  12. Prudnikova, T.N. The Tarim Basin and the transformation of its landscapes. Arid. Ecosyst. 2019, 9, 157–165. [Google Scholar] [CrossRef]
  13. Moholdt, G.; Nuth, C.; Hagen, J.O.; Kohler, J. Recent elevation changes of Svalbard glaciers derived from ICESat laser altimetry. Remote Sens. Environ. 2010, 114, 2756–2767. [Google Scholar] [CrossRef]
  14. Nuth, C.; Moholdt, G.; Kohler, J.; Hagen, J.O.; Kääb, A. Svalbard glacier elevation changes and contribution to sea level rise. J. Geophys. Res. Earth Surf. 2010, 115. [Google Scholar] [CrossRef] [Green Version]
  15. 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]
  16. Chao, N.F.; Chen, G.; Wang, Z.T. Impact of Eastern Tibetan Plateau Glacier melt on land water storage change across the Yangtze river basin. J. Hydrol. Eng. 2020, 25, 05020002. [Google Scholar] [CrossRef]
  17. 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, 419, 495–498. [Google Scholar] [CrossRef]
  18. Parrella, G.; Hajnsek, I. TanDEM-X dual-polarization time series for the estimation of snow accumulation on an alpine glacier. In Proceedings of the TerraSAR-X/TanDEM-X Science Team Meeting 2019, Oberpfaffenhofen, Germany, 21–24 October 2019. [Google Scholar]
  19. Strozzi, T.; Caduff, R.; Jones, N.; Barboux, C.; Delaloye, R.; Bodin, X.; Kaab, A.; Matzler, E.; Schrott, L. Monitoring rock glacier kinematics with satellite synthetic aperture radar. Remote Sens. 2020, 12, 559. [Google Scholar] [CrossRef] [Green Version]
  20. Huss, M.; Farinotti, D. Distributed ice thickness and volume of all glaciers around the globe. J. Geophys. Res. Earth Surf. 2012, 117. [Google Scholar] [CrossRef]
  21. Huang, T.; Li, J.; Menenti, M.; Jing, L.; Hu, G. A new method to estimate changes in glacier surface elevation based on polynomial fitting of sparse ICESat—GLAS footprints. Sensors 2017, 17, 1803. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, X.Y.; Edwards, R. Anthropogenic sulfate and nitrate signals in snow from Bogda Glacier, Eastern Tianshan. J. Earth Sci. 2011, 22, 490–502. [Google Scholar] [CrossRef]
  23. Li, B.L.; Zhu, A.X.; Zhang, Y.C.; Pei, T.; Qin, C.Z.; Zhou, C.H. Glacier change over the past four decades in the middle Chinese Tien Shan. J. Glaciol. 2006, 52, 425–432. [Google Scholar] [CrossRef] [Green Version]
  24. Wang, P.Y.; Li, Z.Q.; Li, H.L.; Cao, M.; Wang, W.B.; Wang, F.T. Glacier No. 4 of Sigong River over Mt. Bogda of eastern Tianshan, central Asia: Thinning and retreat during the period 1962–2009. Environ. Earth Sci. 2012, 66, 265–273. [Google Scholar] [CrossRef]
  25. Guo, W.Q.; Liu, S.Y.; Xu, L.; Wu, L.Z.; Shangguan, D.H.; Yao, X.J.; Wei, J.F.; Bao, W.J.; Yu, P.C.; Liu, Q.; et al. The second Chinese glacier inventory: Data, methods and results. J. Glaciol. 2015, 61, 357–372. [Google Scholar] [CrossRef] [Green Version]
  26. Shen, Y.P.; Wang, G.Y.; Ding, Y.J.; Su, H.C.; Mao, W.F.; Wang, S.D.; Mamatkanov, D.M. Changes in Merzbacher Lake of Inylchek Glacier and glacial flash floods in Aksu River Basin, Tianshan during the period of 1903–2009. J. Glaciol. Geocryol. 2009, 31, 993–1002. [Google Scholar]
  27. Li, K.M.; Li, H.L.; Wang, L.; Gao, W.Y. On the Relationship between Local Topography and Small Glacier Change under Climatic Warming on Mt. Bogda, Eastern Tian Shan, China. J. Earth Sci. 2011, 22, 515–527. [Google Scholar]
  28. Scheuchl, B.; Mouginot, J.; Rignot, E.; Morlighem, M.; Khazendar, A. Grounding line retreat of Pope, Smith, and Kohler Glaciers, West Antarctica, measured with Sentinel-1a radar interferometry data. Geophys. Res. Lett. 2016, 43, 8572–8579. [Google Scholar] [CrossRef] [Green Version]
  29. Sanchez-Gamez, P.; Navarro, F.J. Glacier surface velocity retrieval using D-InSAR and offset tracking techniques applied to ascending and descending passes of Sentinel-1 data for Southern Ellesmere Ice Caps, Canadian Arctic. Remote Sens. 2017, 9, 442. [Google Scholar] [CrossRef] [Green Version]
  30. Surazakov, A.B.; Aizen, V.B. Estimating volume change of mountain glaciers using SRTM and map-based topographic data. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2991–2995. [Google Scholar] [CrossRef]
  31. Ye, Q.H.; Zong, J.B.; Tian, L.D.; Cogley, J.G.; Song, C.Q.; Guo, W.Q. Glacier changes on the Tibetan Plateau derived from Landsat imagery: Mid-1970s-2000-13. J. Glaciol. 2017, 63, 273–287. [Google Scholar] [CrossRef] [Green Version]
  32. Nicolas, J.-M.; Vasile, G.; Gay, M.; Tupin, F.; Trouve, E. SAR processing in the temporal domain: Application to direct interferogram generation and mountain glacier monitoring. Can. J. Remote. Sens. 2007, 33, 52–59. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, X.W.; Liu, L.; Zhao, L.; Wu, T.H.; Li, Z.Q.; Liu, G.X. Mapping and inventorying active rock glaciers in the northern Tien Shan of China using satellite SAR interferometry. Cryosphere 2017, 11, 997–1014. [Google Scholar] [CrossRef] [Green Version]
  34. Martone, M.; Braeutigam, B.; Rizzoli, P.; Gonzalez, C.; Bachmann, M.; Krieger, G. Coherence evaluation of TanDEM-X interferometric data. ISPRS J. Photogramm. 2012, 73, 21–29. [Google Scholar] [CrossRef] [Green Version]
  35. Li, J.L.; Sheng, Y.W. An automated scheme for glacial lake dynamics mapping using Landsat imagery and digital elevation models: A case study in the Himalayas. Int. J. Remote Sens. 2012, 33, 5194–5213. [Google Scholar] [CrossRef]
  36. Zhang, G.Q.; Li, J.L.; Zheng, G.X. Lake-area mapping in the Tibetan Plateau: An evaluation of data and methods. Int. J. Remote Sens. 2017, 38, 742–772. [Google Scholar] [CrossRef]
  37. Wang, D.; Kaab, A. Modeling Glacier Elevation Change from DEM Time Series. Remote Sens. 2015, 7, 10117–10142. [Google Scholar] [CrossRef] [Green Version]
  38. Fattahi, H.; Amelung, F. DEM Error Correction in InSAR Time Series. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4249–4259. [Google Scholar] [CrossRef]
  39. Wang, P.Y.; Li, Z.Q.; Wang, W.B.; Li, H.L.; Wang, F.T. Glacier volume calculation from ice-thickness data for mountain glaciers—A case study of Glacier No. 4 of Sigong River over Mt. Bogda, Eastern Tianshan, Central Asia. J. Earth Sci. 2014, 25, 371–378. [Google Scholar] [CrossRef]
  40. Fistook, J.L.; Head, J.W.; Marchant, D.R.; Forget, F. Tropical mountain glaciers on Mars: Altitude-dependence of ice accumulation, accumulation conditions, formation times, glacier dynamics, and implications for planetary spin-axis/orbital history. Icarus 2008, 198, 305–317. [Google Scholar] [CrossRef]
  41. Wang, P.Y.; Li, Z.Q.; Wang, W.B.; Li, H.L.; Wu, L.Z.; Huai, B.J.; Zhou, P.; Jin, S.; Wang, L. Comparison of changes in glacier area and thickness on the northern and southern slopes of Mt. Bogda, Eastern Tianshan Mountains. J. Appl. Geophys. 2016, 132, 164–173. [Google Scholar] [CrossRef]
  42. Wang, P.Y.; Li, Z.Q.; Huai, B.J.; Wang, W.B.; Li, H.L.; Wang, L. Spatial variability of glacial changes and their effects on water resources in the Chinese Tianshan Mountains during the last five decades. J. Arid Land 2015, 7, 717–727. [Google Scholar] [CrossRef] [Green Version]
  43. Zhang, M.J.; Wang, S.J.; Li, Z.Q.; Wang, F.T. Analysis of glacier area in China under the background of climate change in the past 50 years. Acta Geogr. Sin. 2011, 66, 1155–1165. [Google Scholar]
  44. He, Y.; Yang, T.B.; Chen, J.; Huang, Q. Remote sensing monitoring of glacier change in Bogda peak area of the East Tianshan Mountains from 1972 to 2013. Geogr. Sci. 2015, 35, 925–932. [Google Scholar]
  45. Zhou, Y.G.; Zhao, R.F.; Zhang, L.H.; Zhao, M. Remote sensing monitoring of glacier and snow change and analysis of influencing factors in Bogda peak area. Arid Land Geogr. 2019, 42, 1395–1403. [Google Scholar]
  46. Wang, H.W.; Huang, C.L.; Hao, X.H.; Zhang, P.; Hou, J.L. Analysis of the influencing factors of the temporal and spatial changes of snow cover in the northern Xinjiang. J. Glaciol. Geocryol. 2014, 36, 508–516. [Google Scholar]
  47. Yao, Y.B.; Xiao, G.J.; Wang, R.Y.; Zhang, X.Y. Characteristics of climate change in Northwest semi arid region in recent 50 years. Arid Land Geogr. 2009, 32, 159–165. [Google Scholar]
  48. Lu, A.G.; Kang, S.C.; Pang, D.Q.; Wang, T.M. Study on the asynchrony of temperature change in China under global warming. Arid Land Geogr. 2009, 32, 506–511. [Google Scholar]
  49. Ye, Z.X.; Chen, Y.N.; Zhang, X.; Mao, X.H.; Luo, J.Y. Analysis of Tarim River cutoff and future hydrological situation. Arid Land Geogr. 2009, 32, 841–849. [Google Scholar]
  50. Zhang, Q.; Kang, S. Glacier elevation change in the Western Qilian mountains as observed by TerraSAR-X/TanDEM-X images. Geocarto Int. 2019, 36, 1365–1377. [Google Scholar] [CrossRef]
Figure 1. The research area of Bogda Mountain. (a) The administrative zoning map of China. (b) The administrative zoning map of Xinjiang Province. (c) A satellite image of Bogda Mountain in 2000. The red and black numbers show the reference numbers of the glaciers and lakes investigated in this study, respectively.
Figure 1. The research area of Bogda Mountain. (a) The administrative zoning map of China. (b) The administrative zoning map of Xinjiang Province. (c) A satellite image of Bogda Mountain in 2000. The red and black numbers show the reference numbers of the glaciers and lakes investigated in this study, respectively.
Ijerph 18 06374 g001
Figure 2. The central flow of data processing and results analysis.
Figure 2. The central flow of data processing and results analysis.
Ijerph 18 06374 g002
Figure 3. The elevation variation map of Glacier No. 14 before and after the application of the Logarithmic Fitting Model (LFM). (a) the elevation variation results obtained using the Repeat Orbit Interferometry (ROI) method. (b) the improved results by LFM.
Figure 3. The elevation variation map of Glacier No. 14 before and after the application of the Logarithmic Fitting Model (LFM). (a) the elevation variation results obtained using the Repeat Orbit Interferometry (ROI) method. (b) the improved results by LFM.
Ijerph 18 06374 g003aIjerph 18 06374 g003b
Figure 4. The LFM curve that was used to improve the ROI result. The green points are the sampling points and the red line is the fitting line obtained using the LFM.
Figure 4. The LFM curve that was used to improve the ROI result. The green points are the sampling points and the red line is the fitting line obtained using the LFM.
Ijerph 18 06374 g004
Figure 5. A map showing the change in elevation of glaciers on Bogda Mountain from 2000 to 2017. The red boundaries are the lake boundaries and the black lines are the isoline at an altitude of 500 m.
Figure 5. A map showing the change in elevation of glaciers on Bogda Mountain from 2000 to 2017. The red boundaries are the lake boundaries and the black lines are the isoline at an altitude of 500 m.
Ijerph 18 06374 g005
Figure 6. The elevation variation of the target glaciers from 2000 to 2017. The red boundaries are the lake boundaries.
Figure 6. The elevation variation of the target glaciers from 2000 to 2017. The red boundaries are the lake boundaries.
Ijerph 18 06374 g006
Figure 7. The elevation variation of the 25 studied glaciers in 2000–2017. (a) The maximum and minimum elevation variation. (b) The average elevation variation and SD.
Figure 7. The elevation variation of the 25 studied glaciers in 2000–2017. (a) The maximum and minimum elevation variation. (b) The average elevation variation and SD.
Ijerph 18 06374 g007
Figure 8. The annual average temperature and precipitation anomaly change in the Bogda Mountain area from 2000 to 2017. (a) The annual average precipitation anomaly change in the Dabancheng Station from 2000 to 2017. (b) The annual average temperature change in the Dabancheng Station from 2000 to 2017. (c) The annual average precipitation anomaly change in the Kumish Station from 2000 to 2017. (d) The annual average temperature change in the Kumish Station from 2000 to 2017. (e) The annual average precipitation anomaly change in the Turpan Station from 2000 to 2017. (f) The annual average temperature change in the Turpan Station from 2000 to 2017.
Figure 8. The annual average temperature and precipitation anomaly change in the Bogda Mountain area from 2000 to 2017. (a) The annual average precipitation anomaly change in the Dabancheng Station from 2000 to 2017. (b) The annual average temperature change in the Dabancheng Station from 2000 to 2017. (c) The annual average precipitation anomaly change in the Kumish Station from 2000 to 2017. (d) The annual average temperature change in the Kumish Station from 2000 to 2017. (e) The annual average precipitation anomaly change in the Turpan Station from 2000 to 2017. (f) The annual average temperature change in the Turpan Station from 2000 to 2017.
Ijerph 18 06374 g008
Figure 9. A graph showing the relationship between glacier area, glacier altitude, and average glacier elevation variation between 2000 and 2017.
Figure 9. A graph showing the relationship between glacier area, glacier altitude, and average glacier elevation variation between 2000 and 2017.
Ijerph 18 06374 g009
Figure 10. The relationship between glacier aspect and glacier altitude from 2000 to 2017.
Figure 10. The relationship between glacier aspect and glacier altitude from 2000 to 2017.
Ijerph 18 06374 g010
Figure 11. The relationship between glacier aspect and glacier elevation from 2000 to 2017.
Figure 11. The relationship between glacier aspect and glacier elevation from 2000 to 2017.
Ijerph 18 06374 g011
Figure 12. The relationship between the average elevation variation and altitude between 2000 and 2017. The green points are the average values for each glacier, and the red line is the fitting line that was used to study the relationships between the altitude and average elevation variation for each glacier.
Figure 12. The relationship between the average elevation variation and altitude between 2000 and 2017. The green points are the average values for each glacier, and the red line is the fitting line that was used to study the relationships between the altitude and average elevation variation for each glacier.
Ijerph 18 06374 g012
Table 1. Information about the glaciers used in the study.
Table 1. Information about the glaciers used in the study.
Glacier No.Glacier IDMaximum Elevation (m)Minimum Elevation (m)Average Elevation (m)Average Aspect (°)
15Y725B00104218.4 3846.9 4032.7 97.3
25Y725B00105393.7 3440.6 4049.7 46.7
35Y725B00105031.3 3405.2 3844.6 28.4
45Y725B00095096.1 3636.5 4196.7 29.5
55Y725B00113939.0 3670.6 3790.2 11.3
65Y725B00075304.5 4638.2 5055.2 29.8
75Y725B00084638.3 3344.7 3630.1 16.3
85Y725B00064319.3 3570.7 3908.5 17.2
95Y725C00074020.7 3592.9 3809.0 350.3
105Y725C00034067.2 3663.5 3840.9 29.2
115Y725D00044323.3 3657.4 3956.4 340.8
125Y812B00104381.9 3553.9 3875.8 16.2
125Y812B00114704.8 3548.1 3906.8 325.7
145Y725D00055428.6 3532.6 4150.7 322.0
155Y812B00224674.5 3604.1 4199.3 330.6
165Y812B00154634.1 3667.8 4027.3 156.7
175Y813B00085208.1 3401.4 4200.6 187.3
185Y813B00114994.0 3719.0 4164.0 197.2
195Y813B00154290.1 3955.1 4185.8 314.2
205Y813B00124421.7 3905.1 4136.6 279.3
215Y813B00175293.8 3684.4 4229.1 194.6
225Y813C00064469.0 3600.8 3933.5 47.6
235Y813C00105321.8 3400.4 3917.7 57.5
245Y813C00014140.4 3603.2 3909.4 36.5
255Y813C00124173.2 3730.9 3916.1 47.6
Table 2. Information of the Sentinel-1A and Landsat satellite data that were used in this article.
Table 2. Information of the Sentinel-1A and Landsat satellite data that were used in this article.
SAR DateTemporal BaselineSpatial BaselineLandsat Image DateSensorPath/Row
29 August 201712 days47.679 m2 September 2000TM142/030
10 September 2017
24 August 201812 days100.872 m29 September 2017OLI142/030
05 September 2018
Note: SAR: Synthetic-aperture radar.
Table 3. A comparison between the elevation variations of glaciers and lakes and the SD of the lake elevation variation.
Table 3. A comparison between the elevation variations of glaciers and lakes and the SD of the lake elevation variation.
Glacier No.Lake No.Lake Area (m2)Lake Elevation Variation (m)Glacier Elevation Variation (m)SD of Lake Elevation (m)
1, 2, 3641,6007.09−14.12±0.5
712136,0004.15−23.13±0.08
111142,0003.26−12.80±0.004
128248,80015.71−14.96±0.9
17350,4007.77−7.59±0.04
21436,40016.32−5.95±0.00015
23529,2008.08−14.56±0.44
23918,80014.50−14.56±2.30
14245,60018−9.40±1.40
141049,2002.66−9.40±0.07
131659,6007.56−10.81±3.75
141659,6007.56−17.43±3.75
-7380,000−6.52-±0.00042
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Du, W.; Shi, N.; Xu, L.; Zhang, S.; Ma, D.; Wang, S. Monitoring the Spatiotemporal Difference in Glacier Elevation on Bogda Mountain from 2000 to 2017. Int. J. Environ. Res. Public Health 2021, 18, 6374. https://doi.org/10.3390/ijerph18126374

AMA Style

Du W, Shi N, Xu L, Zhang S, Ma D, Wang S. Monitoring the Spatiotemporal Difference in Glacier Elevation on Bogda Mountain from 2000 to 2017. International Journal of Environmental Research and Public Health. 2021; 18(12):6374. https://doi.org/10.3390/ijerph18126374

Chicago/Turabian Style

Du, Weibing, Ningke Shi, Linjuan Xu, Shiqiong Zhang, Dandan Ma, and Shuangting Wang. 2021. "Monitoring the Spatiotemporal Difference in Glacier Elevation on Bogda Mountain from 2000 to 2017" International Journal of Environmental Research and Public Health 18, no. 12: 6374. https://doi.org/10.3390/ijerph18126374

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