Albedo, which impacts the surface heat budget of a glacier, is a dominant factor governing snow and ice melt variability. We analyzed variations in albedo on the Dongkemadi Glacier (92°4.94′E, 33°05.88′N) using MODIS daily snow albedo products. The results show that the summer albedo over the accumulation zone and the ablation zone all presented a decreasing trend from 2002 to 2012, and both precipitation and summer temperature were critical factors that impacted variations in the albedo of the glacier. We linked the albedo variation with the measured mass balance of the glacier and found that 2006 and 2010 were high-melt years at the Dongkemadi Glacier, with lower albedo and a more negative mass balance. The area around the equilibrium line altitude (ELA) on the glacier was very sensitive to albedo feedback. The albedo change rate reached approximately -34% when the summer average temperature increased about 1 °C and the summer precipitation declined 94.5 mm. We built a linear model between the annual lowest albedo and the annual mass balance of the Xiao Dongkemadi Glacier, one branch of the Dongkemadi Glacier. This model was used to evaluate the annual mass balance of another branch of the Dongkemadi Glacier, the Da Dongkemadi Glacier. The mass balance of the Da Dongkemadi Glacier was significantly negative in 2006 (-768.7 mm w.e.) and 2010 (-696.9 mm w.e.). This work indicates that the glacier surface annual lowest albedo may be a reasonable proxy for evaluating the inter-annual variability of the mass balance of the glacier.
Surface albedo is referred to as bihemispherical broadband albedo and is the ratio of the radiant flux reflected from a unit surface area into the whole hemisphere in relation to the incident radiant flux of hemispherical angular extent (Schaepman-Strub et al., 2006). The albedo of glaciers varies both temporally and spatially. It can range from more than 90% for fresh snow down to approximately 20% for dirty ice (Oerlemans and Knap, 1998). The variation in albedo on a glacier primarily depends on two aspects, the physical properties of snow and ice, and the solar radiation as affected by clouds (Kylling et al., 2000). Correctly measuring and predicting albedo is critical not only because solar radiation is a significant aspect of the surface heat budget, but also because the ice albedo feedback can amplify ongoing changes in the cover (Perovich and Polashenski, 2012). Thus, it is necessary to monitor variations of surface albedo on glaciers to understand cryospheric changes and climate feedback. Furthermore, albedo is an important parameter in climate models and in the energy balance models of glaciers or ice sheets. With the development of these models, spatial and temporal estimates of the albedo of ice and snow are required (Liang et al., 2005; Pedersen and Winther, 2005; Munneke et al., 2011).
Field measurements can provide reliable data for glacier albedo variation analysis. Usually, automatic weather stations (AWSs) installed at the glacier surface are used to monitor the temporal variation of albedo. Field measurements are invaluable but are difficult to acquire because glaciers or ice sheets are usually distributed in high-altitude or high-latitude regions. Another limitation of point field measurements is that they only capture the data in a point on the glacier surface, but the glacier surface is heterogeneous. In recent years, remote sensing has been shown to be a desirable and efficient way to detect the characteristics of glaciers, such as for monitoring snow cover (Immerzeel et al., 2009; Fassnacht et al., 2012), ice thickness (Kwok and Rothrock, 2009), and mass balance (Ramillien et al., 2006). Stroeve et al. (2001) provided temporally and spatially consistent estimates of the Greenland Ice Sheet albedo using the AVHRR Polar Pathfinder (APP) albedo data set; they analyzed albedo variation spanning 1981–1998. Jiang et al. (2010) investigated the spatiotemporal variation in albedo on the Qiyi Glacier in the Qilian Mountains of China, based on spectral pyranometers and data sets from an AWS. Tedesco et al. (2011) analyzed the role of low albedo in the 2010 melting record in Greenland using remote-sensing data, surface observations, and output from a regional atmosphere model. Flanner et al. (2011) concluded that the albedo feedback from the Northern Hemisphere cryosphere falls between 0.3 and 1.1 W m-2 K-1, using MODIS MCD43C3 data in conjunction with CMIP3 models.
Some researchers have deliberated linking glacier albedo and glacier annual mass balance. Greuell et al. (2007) estimated annual anomalies of the surface mass balance of glaciers on Svalbard, a Norwegian archipelago in the Arctic Ocean, for the period 2000–2005, by calculating the so-called “satellite-derived mass balance” from a time series of satellite-derived surface albedos. Historical mass balance gradients of glaciers were combined with MODIS-derived equilibrium line altitutdes (ELAs) to estimate annual mass change at the Columbia, Lillooet, and the Sittakanay ice fields in British Columbia, Canada (Shea et al., 2012). Dumont et al. (2012) processed one decade (2000–2009) of MODIS data to create a time series of albedo maps of the Saint Sorlin Glacier in the western Alps of France during the ablation season, and concluded that monitoring of the glacier surface albedo may provide a useful means to evaluate the inter-annual variability of the glacier mass balance.
In China, previous research primarily focused on monitoring variation in glacier albedo based on point field measurements. For example, based on AWS data, Yang et al. (2010) analyzed the temporal variation in albedo on the Rongbuk Glacier in the Himalaya Mountains over one year. Sun et al. (2012) analyzed the correlation of albedo variation and glacier melting in the accumulation zone of the Laohugou Glacier No. 12 in the western Qilian Mountains, China. However, to our knowledge, less research has been done concerning the spatial distribution and temporal variation of albedo on glaciers.
The goals of this study were, first, to assess the spatiotemporal variations in albedo on the Dongkemadi Glacier on the Tibetan Plateau based on remote sensing data and to analyze the impact factors and feedback mechanism of the albedo on that glacier. Second, we aimed to build a valid model between albedo and glacier mass balance based on the albedo feedback mechanism on the Xiao Dongkemadi Glacier. Finally, we wanted to produce a mass balance estimation of another branch of the Dongkemadi Glacier (Da Dongkemadi) based on the built model.
The Dongkemadi River basin is a typical river basin of the headwaters of the Yangtze River on the Tibetan Plateau, where the elevation is over 5000 m a.s.l. The basin is a tributary at the upper reaches of the Buqu River near Tanggula Mountain Pass. There are 117 glaciers with a total area of 195.73 km2 and volume of 14.4 km3 in the Buqu River region (Pu et al., 1995). The Dongkemadi Glacier (92°4.94′E, 33°05.88′N) is located at the headwaters of the Dongkemadi River (Fig. 1) and is a subcontinental type of glacier. It is composed of the Da (Greater) Dongkemadi Glacier and the Xiao Dongkemadi Glacier, converging at the glacier terminus. The Da Dongkemadi Glacier's area is 14.63 km2 with a length of 5.4 km. In 1993 the elevation of the terminus of the glacier was 5275 m and the ELA of the glacier was 5600 m (Pu and Yao, 1993). The Xiao Dongkemadi Glacier's area is 1.767 km2 with a length of 2.8 km and an average width of 0.5–0.6 km. The summit and the terminus of the glacier are at elevations of 5926 m and 5380 m a.s.l., respectively, with a relative height difference of 546 m (Yao et al., 2012).
The Dongkemadi Glacier accumulates during the summer monsoon season. The melting period of the glaciers in the Tanggula Mountains is short, being only 21 days, with a daily average air temperature above 0 °C and a maximum air temperature of 3.5 °C (Pu et al., 1995). From AWS data observed near the equilibrium line (5600 m a.s.l.) on the glacier, the average annual air temperature and precipitation are -9.8 °C and 302 mm, respectively (Pu et al., 2008). With a gentle surface and without any avalanche or surface moraines, the glacier is an ideal place for mass balance study. The mass balance of the Xiao Dongkemadi Glacier has been monitored since 1989, with recording still ongoing.
Data and Methods
MODIS ALBEDO PRODUCTS AND PROCESSING
MOD10A1 daily products were retrieved from the Terra platform MODIS sensor from the National Aeronautics and Space Administration (NASA) taken from 5 March 2000, which is available from the National Snow and Ice Data Center (NSIDC). The product contains snow extent, snow albedo, fractional snow cover, and quality assessment data at 500 m resolution, gridded in a sinusoidal map projection. More frequent and higher spatial resolution albedo maps are provided by the product, which is suitable for cryospheric research. Validation efforts have shown that MOD10A1 albedo values are within 10% of in situ values in the Karasu Basin, Turkey (Tekeli et al., 2006). Comparisons with in situ observations at five AWSs in Greenland indicated an overall RMS (root mean square) error of 0.067 for the Terra instrument and an RMS error of 0.075 on Aqua. Stroeve et al. (2006) have suggested that MOD10A1 albedo products capture the natural seasonal cycle in albedo but exhibit significantly more temporal and spatial variability than recorded by ground observations. The products on the continental-type glaciers are also reasonably accurate, with an average error of less than 0.026 for albedo, and an RMS error less than 0.068 for albedo (Wang et al., 2011). In this study area, in situ observations had been taken in field work on 11 May 2011. Our open observation site was located at the end of the Dongkemadi Glacier, at 33°3.8′N, 92°4.7′E with an elevation of 5260 m. Albedo was measured by an albedometer (CMA11, Kipp & Zonen, Delft, Netherlands) with a wavelength range from 310 nm to 2800 nm. The absolute difference between MOD10A1 albedo (0.592) and in situ observation (0.638) was small (0.046) on 11 May 2011.
In this research, MOD10A1 daily albedo products were compiled over the Dongkemadi Glacier from March 2002 to October 2012. Only images in which the selected pixels were not cloudy were considered. About 215 images per year were selected over the period from April to October, covering the whole ablation period. In order to detect the variation of the albedo averaged over the whole glacier surface, only the data of pixels entirely within the glacier boundary were used, and potential mixed pixels on the side of the ablation zone were discarded to reduce biases. A total of 41 pixels (Fig. 2) were chosen for computing the average glacier albedo spatially.
The summer average albedo (from June to August) for the 11 years was calculated to detect spatiotemporal variation in albedo of the glacier during the ablation season. During the melt period, the glacier surface was generally comprised of bare ice at the ablation zone while the surface on the accumulation zone was still covered by snow. Therefore, albedo pixels on the glacier were divided into two parts, the ablation zone and the accumulation zone, based on ELA location and 30 m Aster DEM. We calculated the average albedo for the ablation zone and the accumulation zone, respectively.
The temperature and precipitation in the Dongkemadi River basin were derived from an empirical model because of the lack of long-term meteorological data in the study glacier region. The model was built by polynomial interpolation based on AWS data from 2005 to 2008 in the Dongkemadi River basin and meteorological data of national meteorological stations near the basin. This method improved the accuracy of interpolation, and the fitted linear equation for temperature is (Gao et al., 2011):
where T0 represents the temperature of the Dongkemadi River basin, and T1, T2, T3 and T4 are temperatures recorded at the Anduo, Naqu, Tuotuohe, and Wudaoliang meteorological stations, respectively. The fitted equation for precipitation is:
where P0 represents the precipitation of the Dongkemadi River basin, and P1 and P2 are precipitation recorded at the Anduo and Naqu meteorological stations, respectively. In accordance with the above models, we reconstructed the summer and annual average temperature and precipitation of the Dongkemadi River basin from 2002 to 2012.
The summer (from June to August) mean albedo values of the whole glacier for every year were calculated using the arithmetic mean method. Linear regression analysis was used to assess the albedo trends for each of the selected pixels on the glacier during 2002–2012:
where T(t) is the calculated albedo series, t is the year (t = 1, 2, …, n), and a is a trend term. A positive value of a indicates that there is an increased tendency of albedo, and a negative value of a indicates that there is a decreased trend of albedo. The magnitude of a suggests albedo changes the amount of degree. The values of a and b can be obtained by the least-squares method and the significance can be tested by t-tests.
The average values of the albedo for a 10-day moving window were calculated every fifth day for the daily albedo series. Based on this, the inter-annual variability and seasonal variations of albedos at the Dongkemadi Glacier could be obtained. Albedos in anomaly years were compared with multi-year means in summer, and change rates were calculated to quantify albedo variation.
INTER-ANNUAL AND SEASONAL VARIABILITY IN ALBEDO ON DONGKEMADI GLACIER
Figure 2 shows the spatial variation tendency in the summer average albedo at the glacier. The slope of the albedo variation was gradually larger from the accumulation to the ablation zone. Decreasing tendencies that exceeded the 90% confidence level occurred for 21.9% of the glacier, which was primarily distributed on the terminus and the top of the glacier.
The inter-annual variation in average albedo on the whole Dongkemadi Glacier is shown in Figure 3. Over the ablation zone, albedo declined from 0.508 in 2002 to 0.411 in 2012 (time correlation coefficient = -0.328. p = 0.325) (Fig. 3, part a). The absolute change value (-0.097) of the albedo was more than one standard deviation (σ, ±0.089) and regression residuals (0.071). Over the accumulation zone, the albedo (time correlation coefficient = -0.353, p = 0.287) declined from 0.636 to 0.561 (absolute change value: -0.075) over the same period, and was twice as large as the regression residuals (0.034) and also exceeded one standard deviation (±0.061). The annual decline of albedo on the glacier is clearly shown in Fig. 3, part b, and the variation in average albedo in the ablation zone was more pronounced than that at the accumulation zone because of differences in meteorological conditions and in the physical properties of snow and ice.
Seasonal variations in albedo for all 11 years followed a similar pattern (Fig. 4). The mean albedo on the glacier during that time period declined from 0.73 in April to the lowest value of 0.5 by early August. Then the albedo began to increase until October. 2006, 2010, and 2012 had extremely low albedo values on the glacier. The annual lowest albedos occurred in late July, with a mean value of 0.5 from 2002 to 2012. We defined albedo anomalous years as the years in which the annual lowest albedo values were more than one standard deviation (±1σ) from the mean value. On this basis, 2006, 2010, and 2012 were anomalous years with low albedo, and the annual lowest albedo values were all less than 0.4 in July. Moreover, lower albedo values (0.54 in 2006, 0.35 in 2010, and 0.45 in 2012) appeared in June in the anomalous years.
RECONSTRUCTED TEMPERATURE AND PRECIPITATION AND LINKAGE WITH THE GLACIER ALBEDO
Table 1 shows the relationship between the mean summer albedo in the ablation and accumulation zones with precipitation and temperature. The annual and summer precipitations were significantly positively correlated with the albedo in the ablation and accumulation zones (p < 0.01). There was also a significant correlation between mean summer temperature and the albedo in the ablation zone (r = -0.88, p < 0.005) and the accumulation zone (r = -0.93, p < 0.005), respectively, whereas there was no significant correlation between mean annual temperature and the albedo in either the ablation zone or the accumulation zone. These relationships suggest that both the precipitation amount and the mean summer temperature affected the summer albedo on the glacier. Precipitation affected the total mass accumulation of a glacier in a year, and summer temperature affected glacier melt.
There were apparent high summer temperatures in 2006 (4.5 °C) and 2010 (4.7 °C), which were more than 1σ from the mean (Fig. 5, part a). Also, less summer precipitation occurred in 2006 (242 mm) and 2010 (263 mm) (Fig. 5, part b), which were lower than 1σ from the mean. Concurrently, lower summer albedos on the glacier also occurred in 2006 and 2010 as compared to other years (Fig. 5).
Figure 6 shows the spatial distribution of the albedo change rates on the glacier in 2003, 2006, 2008, and 2010. The albedo change rates in 2006 and 2010 were negative, ranging between -34% and -4% (Fig. 6, parts a and b). This suggests that the glacier surface albedo was anomalously low in anomalous warm and dry years. However, when the summer temperature was relatively low and the annual precipitation was relatively rich (i.e., in 2003 and 2008), the albedo on the glacier was high. The change rates of albedo on the glacier ranged from 0.8% to 32% in 2003 and from 2% to 16% in 2008 (Fig. 6, parts c and d).
We also found that the area close to the equilibrium line was sensitive for albedo change. The largest albedo variation was distributed in the middle of the glacier (5500–5600 m), which is around the ELA (Fig. 6). Specifically, the albedo change rate around the ELA reached about -34% in 2006 when the summer temperature was about 1 °C higher and summer precipitation was 94.5 mm lower than the mean values of the 11 years. In 2010, the albedo change rate reached approximately -32% around the ELA, when the summer temperature was about 1.2 °C higher and summer precipitation was 73.4 mm lower than the mean values of the 11 years.
LINKING ALBEDO AND GLACIER MASS BALANCE
Mass balance is a direct and reliable indicator of glacier status. Annual variations of the mass balance observed in the Xiao Dongkemadi Glacier showed that negative mass balance appeared in the past decade. There is a decreasing tendency of mass balance on the Xiao Dongkemadi Glacier, although the trend is not statistically significant (Fig. 7, part a). The annual mass balances were close to -1000 mm w.e. in 2006 and 2010, which were much lower than the average level (-406 mm w.e.) of 2002–2010.
Relationships between summer albedo and temperature and precipitation in the accumulation zone and the ablation zone in summer and the whole year.
Associated with the albedo variations previous described, the annual lowest values of albedo on the glacier reached their minimum values of the 11-years closely corresponding to the year when major negative mass balance appeared. There was a statistically significant positive correlation between annual mass balance and the annual lowest albedos on the Xiao Dongkemadi Glacier (r = 0.94, p < 0.01). Therefore, we established a linear model using the annual lowest albedo as the dependent variable and the annual mass balance of the Xiao Dongkemadi Glacier as the independent variable (Fig.7, part b). This model can explain the relationship between the annual lowest albedo of the glacier and mass balance with a good fit (R2 = 0.88, p < 0.001).
Based on the linear model on the Xiao Dongkemadi Glacier, we reconstructed the mass balance of the Da Dongkemadi Glacier (Table 2). The trend of mass balance of Da Dongkemadi Glacier was similar to that of Xiao Dongkemadi Glacier, but the amount of variation was not as large as that of Xiao Dongkemadi Glacier. The reason may be that the responses of glaciers with different sizes to climate change vary: the smaller of the glacier area, generally the more sensitive, is the climate response of the glacier (Zhang et al., 2011a, 2011b).
ALBEDO DECLINE ON GLACIERS AND POSSIBLE IMPACT FACTORS
The declining linear trend of albedo on the Dongkemadi Glacier from 2002 to 2012 was not statistically significant, but it cannot be neglected. Over the ablation zone of the glacier, albedo declined from 0.508 in 2002 to 0.411 in 2012; over the accumulation zone, the albedo declined from 0.636 to 0.561 over the same period. The declining trend of albedo on other glaciers and impacts of albedo variation on glaciers have also been reported. A persistent drop in Greenland Ice Sheet surface albedo has also been observed in satellite data and climate models from 2000 to 2011, and the ablation zone albedo declined from 0.715 in 2000 to 0.632 in 2011; over the accumulation area, there was a significant linear decline from 0.817 to 0.766 over the same period (As et al., 2012). Oerlemans et al. (2009) reported that surface albedo on the Morteratsch Glacier (Swiss Alps) in summer dropped from 0.32 to 0.15 over the period 2003 to 2006, which resulted in additional removal of about 3.5 m of ice thickness in that period. Decreased albedos on glacier surfaces evidently enhance the absorption of radiation, which may also be responsible for the acceleration of glacier melting in the past 20 years (Li et al., 2007). Climate models suggest plausible estimates for the impact of albedo variation on glaciers. Albedo declines (1.5% in the Arctic and 3% in northern hemisphere land areas) yield a climate forcing of +0.3 W m-2 in the northern hemisphere (Hansen and Nazarenko, 2004), which enhances glacier melting.
Seasonal variation in albedo on the study glacier is clear. Albedo values in April were the highest, when cold, thick snow covered the glacier. As temperatures increased, the snow cover began to warm and melt, thereby decreasing the snow albedo. Until the end of July, when air temperature reached the annual maximum value, meltwater production peaked and ice was exposed to the greatest extent in the ablation zone, and the albedo declined to its lowest value. Then, at the end of the melt season (in late August), snow accumulation brightened the surface of the glacier, and the albedo of the glacier began to increase until October. In the Dongkemadi Glacier, extremely low albedo values were detected at the end of July (in 2006 and 2010), which means that the glacier could have experienced high melt at those years. The lower albedo values that appeared in early June in anomalous years also indicated that the snow melt was stronger than in other years during the same period.
Both precipitation and summer temperature are the dominant factors for albedo variation. The summer temperature was significantly negatively correlated with the summer albedo on the Dongkemadi Glacier, in both the accumulation zone and the ablation zone, during the period of 2002–2012. Air temperature physically affects the reflectance properties of the surface material on the glacier, such as snow metamorphism and increases in the concentration of light-absorbing impurities in snow or on ice (Brock et al., 2000a). Temperature rise in the summer leads to snow grain metamorphosis change in the accumulation zone. The thickness of the firn layer becomes thin, the amount of fine firn with high reflectivity is reduced, and the amount of coarse firn with low reflectivity is increased. Thus, high temperature would cause a decline in albedo in the accumulation zone as a whole. In the ablation zone, higher summer temperature would accelerate ice melting and increase the water content of the glacier surface, lowering the surface albedo.
Estimated mass balance on the Da Dongkemadi Glacier using minimum 10-day running mean albedo.
Precipitation primarily affects the annual mass accumulation of glaciers. Decreasing precipitation reduces the glacial accumulation in the summer and exposes ice over time, which accelerates the melting process of glaciers. At the same time, meltwater on the surface of glaciers contains many dust particles during transport, and microorganisms attached to the surface of glaciers proliferate with the rise of temperature. These deepen the ice color, thereby significantly reducing the surface reflectivity. This phenomenon is more likely to happen over the ablation zone. Figure 8 compares two groups of photos taken on the Dongkemadi Glacier in May 2004 and May 2010. The surface of the glacier was dirtier in May 2010 compared to May 2004. From research in this study, summer average albedos and lowest albedos all were very low in 2010 compared to 2004; and the mass balances also presented a major negative balance. So, lower albedo in 2010 was more likely to be affected by dust on the surface of the glacier. A similar result was also found in a previous study. Fujita et al. (2000) observed a drastic melt on the Dongkemadi Glacier from May to June in 1995, which they considered to have been caused by low surface albedo due to dust deposition, even under low-air-temperature conditions. Since 2000, declining trends in the albedo of Himalayan glaciers have been detected; their upper areas have become darker, and light-absorbing constituents (e.g., black carbon and dust) could be partly responsible for this phenomenon (Ming et al., 2012).
MECHANISM OF ALBEDO CHANGE AROUND THE ELA
The ELA is an important parameter in describing a glacier because it helps to define the relation between the local climate and the distribution of glacier mass (Fountain et al., 1999). Compared to other areas on the study glacier, the change rates in albedo around the ELA were larger. The area near the ELA on the glacier is sensitive and may fall or rise due to temperature changes during the ablation season. A glacier advances when the ELA falls and retreats when the ELA rises (Wang et al., 2010). In a positive mass balance year when the glacier ELA falls, the original bare ice surface is replaced by snow or firn and the surface albedo rises steeply (Yao, 1987), and the energy of the glacier thus decreases. Conversely, in a negative mass balance year when the glacier ELA rises, the snow or firn area is uplifted again and the surface albedo decreases to the original area. This mechanism can explain why the albedo change in the area near the ELA is larger than in other areas on the glacier.
POTENTIAL TO ESTIMATE ANNUAL MASS BALANCE OF THE GLACIER BY THE ANNUAL LOWEST ALBEDO
The change in glacier mass balance on the Tibetan Plateau can serve as a signal of local climatic change. Herein a strong correlation between the measured annual mass balance and the annual lowest albedo of the glacier was established based on the physical mechanism. Glacier surface melt is highly sensitive to albedo variations, since surface energy balances during the ablation season are commonly dominated by receipt of solar radiation (Brock et al., 2000b). Similar results also reported that the modeled surface mass balance is highly sensitive to surface albedo (van Angelen et al., 2012). Positive melt-albedo feedback results in lower albedo on the glacier in high-melt years. For the Dongkemadi Glacier, the years of 2006, 2010, and 2012 were high-melt years, which can be proved by the field measurement of ELA. The higher observed ELA (5840 m; unpublished data) occurred on the Xiao Dongkemadi Glacier in the lower mass balance year 2005/2006, which was about 200 m higher than the average value of 5650 m ( 1989–2002 ) (Pu et al., 2008). Similar results were reported on the other glaciers on the Tibetan Plateau, and a higher observed ELA also occurred on the Qiyi Glacier located in the Qilian Mountains in the north of the Tibetan Plateau in 2005/2006, reaching the highest level (5131 m) close to the glacier summit in recent years (Wang et al., 2010). The observed ELA on the Xiao Dongkemadi Glacier in 2009/2010 reached 5780 m (unpublished data), which was about 130 m higher than the average value. It can be concluded that the year when the annual lowest albedo reached to minimum was consistent with the year when the major negative mass balance appeared and the ELA reached the highest position during the period of 2002–2012.
Even though annual mass balance on the Da Dongkemadi Glacier was estimated using the linear model based on annual lowest albedo of the glacier, this method needs to be tested on several glaciers in China. Empirical models for different glaciers in the world may vary, and the study methods to connect albedo and mass balance are slightly different; this is a valid subject for future glacier researches.
Monitoring the spatiotemporal variation in albedo on glaciers is a significant task, and it is valuable for simulating glacier mass balance and reflecting regional climate change. Our 11-year series of albedo values derived from satellite images and 8 years of annual measured mass balance on the Xiao Dongkemadi Glacier were analyzed for the albedo tendency of the glacier surface on multi-annual timescales and to link the two variables. The summer mean albedo over the Dongkemadi Glacier declined from 2002 to 2012, and albedo variation in the ablation zone during the period was more pronounced than that in the accumulation area because of differences in meteorological conditions and the physical properties of snow and ice. Precipitation, summer temperature, and dust on the glacier surface all play critical roles in summer mean albedo variation on the Dongkemadi Glacier. In this study, several anomalous years were found on the Dongkemadi Glacier, which were high-melt years with lower albedo and major negative mass balance. In addition, we found that the region around the ELA on the glacier was very sensitive to albedo change. There was a strong positive correlation between annual mass balance and the annual lowest albedo on the Xiao Dongkemadi Glacier during the ablation season. Consequently, it can be considered that the annual lowest albedo of a glacier could be a reasonable proxy for the annual mass balance of a glacier. If global warming continues in the future, glacier albedos in summer would become lower and mass balance would become more negative, resulting in continuous glacier retreat.
Admittedly, much uncertainty is introduced in the study process. For instance, the MODIS albedo algorithm does not provide a surface albedo estimate under cloudy skies. Our inspection of the raw MOD10A1 images revealed an abundance of residual cloud artifacts (shadows, contrails, and thin clouds) in the albedo product. The combination of other remote sensing data and a proper suggested model is necessary for further cryospheric research.
The authors thank NASA for providing image data, and the China Meteorological Data Sharing Service System for providing meteorological data. This work was supported by the National Natural Science Foundation of China (grant no. 41190084), Strategic Priority Research Program (B) of Chinese Academy of Sciences (grant no. XDB03030200), and the Young Talents funding of the Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences (grant no. Y35111001). We also thank the anonymous reviewers and editors, whose comments and suggestions helped us to improve the quality of this paper.