Glaciers are crucial water resources in the Third Pole (the Tibetan Plateau and its surroundings) and are shrinking in response to climate change. Glacier albedo is an expression of glacier interactions with climate and dust/black carbon, and albedo reduction enhances glacier mass loss, but its changes and potential drivers remain poorly quantified. We leverage satellite observations to explore the variability of glacier albedo and understand its sensitivity to potential drivers and its future evolution. We find that glacier albedo has declined during 2001–2020, but high interannual variability is also an important signal. These variations are highly sensitive to air temperature and snow conditions and to nearby dust/black carbon emission sources. Future changes to these drivers will lead to further decreases of 2.9%–12.5% in glacier albedo by 2100 under different warming scenarios. These findings highlight the importance of albedo in glacier future evolution and the urgency of action to mitigate climate warming.

Currently, about 12ĝ€¯%-13ĝ€¯% of High Mountain Asia's glacier area is debris-covered, which alters its surface mass balance. However, in regional-scale modelling approaches, debris-covered glaciers are typically treated as clean-ice glaciers, leading to a bias when modelling their future evolution. Here, we present a new approach for modelling debris area and thickness evolution, applicable from single glaciers to the global scale. We derive a parameterization and implement it as a module into the Global Glacier Evolution Model (GloGEMflow), a combined mass-balance ice-flow model. The module is initialized with both glacier-specific observations of the debris' spatial distribution and estimates of debris thickness. These data sets account for the fact that debris can either enhance or reduce surface melt depending on thickness. Our model approach also enables representing the spatiotemporal evolution of debris extent and thickness. We calibrate and evaluate the module on a selected subset of glaciers and apply GloGEMflow using different climate scenarios to project the future evolution of all glaciers in High Mountain Asia until 2100. Explicitly accounting for debris cover has only a minor effect on the projected mass loss, which is in line with previous projections. Despite this small effect, we argue that the improved process representation is of added value when aiming at capturing intra-glacier scales, i.e. spatial mass-balance distribution. Depending on the climate scenario, the mean debris-cover fraction is expected to increase, while mean debris thickness is projected to show only minor changes, although large local thickening is expected. To isolate the influence of explicitly accounting for supraglacial debris cover, we re-compute glacier evolution without the debris-cover module. We show that glacier geometry, area, volume, and flow velocity evolve differently, especially at the level of individual glaciers. This highlights the importance of accounting for debris cover and its spatiotemporal evolution when projecting future glacier changes.

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 obser-vations. 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 analyz-ing glacier albedo at large spatial scales.

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.