Search

Debris-covered glaciers have contrasting melting mechanisms and climate response patterns if compared with debris-free glaciers and thus show a unique influence on the hydrological process. Based on high-resolution satellite images and unpiloted aerial vehicle surveys, this study investigated the dynamic changes of Zhuxi Glacier, a thick debris-covered glacier in the southeastern Tibetan Plateau. Our result shows that the whole glacier can be divided into the active regime and stagnant regime along the elevation of 3400 m a.s.l. The mean surface velocity of the active regime was 13.1 m yr−1, which was five times higher than that of the stagnant regime. The surface-lowing rate of this debris-covered glacier reaches more than 1 m yr−1 and displays an accelerating trend. The majority of ice loss concentrates around ice cliffs and supraglacial ponds, the ablation hotspots. These hotspots can be roughly classified into three types, including persistent, expanding, and shrinking patterns, at different dynamic regimes on the Zhuxi Glacier. With the evolution of these hotpots and glacier dynamic changes, the supraglacial ponds showed significant change, with the total number fluctuating from 15 to 38 and the total area increasing from 1128 m2 to 95790 m2 during the past decade. The recent exponential expansion of the proglacial lake and the significant downwasting of stagnant ice inside the dammed terminus moraine possibly trigger the glacial lake outburst flood and thus threaten the security of livelihoods and infrastructure downstream.

Glaciers are sensitive indicators of climate change and have a significant influence on regional water cycle, human survival and social development. Global warming has led to great changes in glaciers over the High Mountain Asia (HMA) region. Glacier elevation change is a measure of glacier mass balance driven by the processes of energy and mass exchange between the glacier surface and the atmosphere which are influenced by climatic factors and glacier surface properties. In this study, we estimated the inter-annual and intra-annual elevation changes of glaciers in the HMA region in 2003–2020 using Ice, Cloud and land Elevation Satellite (ICESat) data and Shuttle Radar Terrain Mission (SRTM) digital elevation model (DEM) data by developing an “elevation-aspect bin analysis method” that considered the difference of glacier elevation changes in different elevations and aspects of glacier topography. The results showed that: (1) The inter-annual change of glacier elevation in 2003–2020 had large spatial heterogeneity. Glacier elevation reduction mainly occurred in the marginal region of the HMA with the maximum decline in the Nyainqentanglha region, while glacier elevation showed increase in the West Kunlun of inner HMA regions in 2003–2020. The glacier elevation change rate showed an accelerating reduction trend in most of the HMA regions, except in the west HMA where the glacier elevation reduction rate showed slowdown tendency. Specifically, the glacier elevation change rate in the entire HMA was −0.21 ± 0.12 m/year in 2003–2008 and −0.26 ± 0.11 m/year in 2003–2020, respectively. (2) The intra-annual change of HMA glacier elevation in 2019 and 2020 showed obvious spatiotemporal heterogeneity, and the glacier thickening period was gradually delayed from the marginal area to the inner area of the HMA. The glaciers in the western marginal part of the HMA (the Tienshan Mountains, Pamir and Hindu Kush and Spiti Lahaul) and Karakoram thickened in winter or spring, the glaciers in the Nyainqentanglha Mountains exhibited spring accumulation. The glaciers in West Kunlun accumulated in two time periods, i.e., from March to June and from July to September. The glaciers in the Inner Tibetan Plateau and Bhutan and Nepal areas experienced spring or summer accumulation, especially in June or July. Moreover, we found that the inter-annual and intra-annual change of glacier elevation could be explained by the changes in temperature and precipitation. A similar analysis can be extended to mountain glaciers in other regions of the world, and glacier change trends could be further explored over a longer time span with the continuous operation of ICESat-2.

This project explored the integrated use of satellite, ground observations and hydrological distributed models to support water resources assessment and monitoring in High Mountain Asia (HMA). Hydrological data products were generated taking advantage of the synergies of European and Chinese data assets and space-borne observation systems. Energy-budget-based glacier mass balance and hydrological models driven by satellite observations were developed. These models can be applied to describe glacier-melt contribution to river flow. Satellite hydrological data products were used for forcing, calibration, validation and data assimilation in distributed river basin models. A pilot study was carried out on the Red River basin. Multiple hydrological data products were generated using the data collected by Chinese satellites. A new Evapo-Transpiration (ET) dataset from 2000 to 2018 was generated, including plant transpiration, soil evaporation, rainfall interception loss, snow/ice sublimation and open water evaporation. Higher resolution data were used to characterize glaciers and their response to environmental forcing. These studies focused on the Parlung Zangbo Basin, where glacier facies were mapped with GaoFeng (GF), Sentinal-2/Multi-Spectral Imager (S2/MSI) and Landsat8/Operational Land Imager (L8/OLI) data. The geodetic mass balance was estimated between 2000 and 2017 with Zi-Yuan (ZY)-3 Stereo Images and the SRTM DEM. Surface velocity was studied with Landsat5/Thematic Mapper (L5/TM), L8/OLI and S2/MSI data over the period 2013–2019. An updated method was developed to improve the retrieval of glacier albedo by correcting glacier reflectance for anisotropy, and a new dataset on glacier albedo was generated for the period 2001–2020. A detailed glacier energy and mass balance model was developed with the support of field experiments at the Parlung No. 4 Glacier and the 24 K Glacier, both in the Tibetan Plateau. Besides meteorological measurements, the field experiments included glaciological and hydrological measurements. The energy balance model was formulated in terms of enthalpy for easier treatment of water phase transitions. The model was applied to assess the spatial variability in glacier melt. In the Parlung No. 4 Glacier, the accumulated glacier melt was between 1.5 and 2.5 m w.e. in the accumulation zone and between 4.5 and 6.0 m w.e. in the ablation zone, reaching 6.5 m w.e. at the terminus. The seasonality in the glacier mass balance was observed by combining intensive field campaigns with continuous automatic observations. The linkage of the glacier and snowpack mass balance with water resources in a river basin was analyzed in the Chiese (Italy) and Heihe (China) basins by developing and applying integrated hydrological models using satellite retrievals in multiple ways. The model FEST-WEB was calibrated using retrievals of Land Surface Temperature (LST) to map soil hydrological properties. A watershed model was developed by coupling ecohydrological and socioeconomic systems. Integrated modeling is supported by an updated and parallelized data assimilation system. The latter exploits retrievals of brightness temperature (Advanced Microwave Scanning Radiometer, AMSR), LST (Moderate Resolution Imaging Spectroradiometer, MODIS), precipitation (Tropical Rainfall Measuring Mission (TRMM) and FengYun (FY)-2D) and in-situ measurements. In the case study on the Red River Basin, a new algorithm has been applied to disaggregate the SMOS (Soil Moisture and Ocean Salinity) soil moisture retrievals by making use of the correlation between evaporative fraction and soil moisture.

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.

Glacier albedo determines the net shortwave radiation absorbed at the glacier surface and plays a crucial role in glacier energy and mass balance. Remote sensing techniques are efficient means to retrieve glacier surface albedo over large and inaccessible areas and to study its variability. However, corrections of anisotropic reflectance of glacier surface have been established for specific shortwave bands only, such as Landsat 5 Thematic Mapper (L5/TM) band 2 and band 4, which is a major limitation of current retrievals of glacier broadband albedo. In this study, we calibrated and evaluated four anisotropy correction models for glacier snow and ice, applicable to visible, near-infrared and shortwave-infrared wavelengths using airborne datasets of Bidirectional Reflectance Distribution Function (BRDF). We then tested the ability of the best-performing anisotropy correction model, referred to from here on as the ‘updated model’, to retrieve albedo from L5/TM, Landsat 8 Operational Land Imager (L8/OLI) and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and evaluated these results with field measurements collected on eight glaciers around the world. Our results show that the updated model: (1) can accurately estimate anisotropic factors of reflectance for snow and ice surfaces; (2) generally performs better than prior approaches for L8/OLI albedo retrieval but is not appropriate for L5/TM; (3) generally retrieves MODIS albedo better than the MODIS standard albedo product (MCD43A3) in both absolute values and glacier albedo temporal evolution, i.e., exhibiting both fewer gaps and better agreement with field observations. As the updated model enables anisotropy correction of a maximum of 10 multispectral bands and is implemented in Google Earth Engine (GEE), it is promising for observing and analyzing glacier albedo at large spatial scales.

Monitoring glacier flow is vital to understand the response of mountain glaciers to environmental forcing in the context of global climate change. Seasonal and interannual variability of surface velocity in the temperate glaciers of the Parlung Zangbo Basin (PZB) has attracted significant attention. Detailed patterns in glacier surface velocity and its seasonal variability in the PZB are still uncertain, however. We utilized Landsat-8 (L8) OLI data to investigate in detail the variability of glacier velocity in the PZB by applying the normalized image cross-correlation method. On the basis of satellite images acquired from 2013 to 2020, we present a map of time-averaged glacier surface velocity and examined four typical glaciers (Yanong, Parlung No.4, Xueyougu, and Azha) in the PZB. Next, we explored the driving factors of surface velocity and of its variability. The results show that the glacier centerline velocity increased slightly in 2017–2020. The analysis of meteorological data at two weather stations on the outskirts of the glacier area provided some indications of increased precipitation during winter-spring. Such increase likely had an impact on ice mass accumulation in the up-stream portion of the glacier. The accumulated ice mass could have caused seasonal velocity changes in response to mass imbalance during 2017–2020. Besides, there was a clear winter-spring speedup of 40% in the upper glacier region, while a summer speedup occurred at the glacier tongue. The seasonal and interannual velocity variability was captured by the transverse velocity profiles in the four selected glaciers. The observed spatial pattern and seasonal variability in glacier surface velocity suggests that the winter-spring snow might be a driver of glacier flow in the central and upper portions of glaciers. Furthermore, the variations in glacier surface velocity are likely related to topographic setting and basal slip caused by the percolation of rainfall. The findings on glacier velocity suggest that the transfer of winter-spring accumulated ice triggered by mass conservation seems to be the main driver of changes in glacier velocity. The reasons that influence the seasonal surface velocity change need further investigation.

Mountain glaciers are excellent indicators of climate change and have an important role in the terrestrial water cycle and food security in many parts of the world. Glaciers are the major water source of rivers and lakes in the Nyainqentanglha Mountains (NM) region, where the glacier area has the second largest extent on the Tibetan Plateau. The potential of the high spatial resolution ZiYuan-3 (ZY-3) Three-Line-Array (TLA) stereo images to retrieve glacier mass balance has not been sufficiently explored. In this study, we optimized the procedure to extract a Digital Elevation Model (DEM) from ZY-3 TLA stereo images and estimated the geodetic mass balance of representative glaciers in the two typical areas of the NM using ZY-3 DEMs and the C-band Shuttle Radar Topography Mission (SRTM) DEM in three periods, i.e., 2000−2013, 2013−2017 and 2000−2017. The results provide detailed information towards better understanding of glacier change and specifically show that: (1) with our new stereo procedure, ZY-3 TLA data can significantly increase point cloud density and decrease invalid data on the glacier surface map to generate a high resolution (5 m) glacier mass balance map; (2) the glacier mass balance in both the Western Nyainqentanglha Mountains (WNM) and Eastern Nyainqentanglha Mountains (ENM) was negative in 2000−2017, and experienced faster mass loss in recent years (2013−2017) in the WNM. Overall, the glaciers in the western and eastern NM show different change patterns since they are influenced by different climate regimes; the glacier mass balances in WNM was −0.22 ± 0.23 m w.e. a−1 and −0.43 ± 0.06 m w.e. a−1 in 2000−2013 and 2013−2017, respectively, while in 2000−2017, it was −0.30 ± 0.19 m w.e. a−1 in the WNM and −0.56 ± 0.20 m w.e. a−1 in the ENM; (3) in the WNM, the glaciers experienced mass loss in 2000−2013 and 2013−2017 in the ablation zone, while in the accumulation zone mass increased in 2000−2013 and a large mass loss occurred in 2013−2017; as regards the ENM, the glacier mass balance was negative in 2000−2017 in both zones; (4) glacier mass balance can be affected by the fractional abundance of debris and glacier slope; (5) the glacier mass balances retrieved by ZY-3 and TanDEM-X data agreed well in the ablation zone, while a large difference occurred in the accumulation zone because of the snow/firn penetration of the X-band SAR signal.