Mountains are important suppliers of freshwater to downstream areas, affecting large populations in particular in High Mountain Asia (HMA). Yet, the propagation of water from HMA headwaters to downstream areas is not fully understood, as interactions in the mountain water cycle between the cryo-, hydro- and biosphere remain elusive. We review the definition of blue and green water fluxes as liquid water that contributes to runoff at the outlet of the selected domain (blue) and water lost to the atmosphere through vapor fluxes, that is evaporation from water, ground, and interception plus transpiration (green) and propose to add the term white water to account for the (often neglected) evaporation and sublimation from snow and ice. We provide an assessment of models that can simulate the cryo-hydro-biosphere continuum and the interactions between spheres in high mountain catchments, going beyond disciplinary separations. Land surface models are uniquely able to account for such complexity, since they solve the coupled fluxes of water, energy, and carbon between the land surface and atmosphere. Due to the mechanistic nature of such models, specific variables can be compared systematically to independent remote sensing observations–providing vital insights into model accuracy and enabling the understanding of the complex watersheds of HMA. We discuss recent developments in spaceborne earth observation products that have the potential to support catchment modeling in high mountain regions. We then present a pilot study application of the mechanistic land surface model Tethys & Chloris to a glacierized watershed in the Nepalese Himalayas and discuss the use of high-resolution earth observation data to constrain the meteorological forcing uncertainty and validate model results. We use these insights to highlight the remaining challenges and future opportunities that remote sensing data presents for land surface modeling in HMA.

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.