[1] A distributed energy balance model was applied to Zongo Glacier, Bolivia (16°S, 6000-4900 m above sea level, 2.4 km 2 ), to investigate atmospheric forcing that controls seasonal variations in the mass balance and in meltwater discharge of glaciers in the outer tropics. Surface energy fluxes and melt rates were simulated for each 20 × 20 m 2 grid cell at an hourly resolution, for the hydrological year 1999-2000, using meteorological measurements in the ablation area. Model outputs were compared to measurements of meltwater discharge, snow cover extent, and albedo at two weather stations set up on the glacier. Changes in melt rate in three distinct seasons were related to snowfall and cloud radiative properties. During the dry season (May to August), the low melt rate was mainly caused by low long-wave emission of the cloudless thin atmosphere found at these high altitudes. From September to December, meltwater discharge increased to its annual maximum caused by an increase in solar radiation, which was close to its summer peak, as well as a decrease in glacier albedo. From January on, melt was reduced by snowfalls in the core wet season via the albedo effect but was maintained thanks to high long-wave emission from convective clouds. The frequent changes in snow cover throughout the long ablation season lead to large vertical mass balance gradients. Annual mass balance depends on the timing and length of the wet season, which interrupts the period of highest melt rates caused by solar radiation.Citation: Sicart, J. E., R. Hock, P. Ribstein, M. Litt, and E. Ramirez (2011), Analysis of seasonal variations in mass balance and meltwater discharge of the tropical Zongo Glacier by application of a distributed energy balance model,

Snow sublimation is a loss of water from the snowpack to the atmosphere. So far, snow sublimation has remained unquantified in the Himalaya, prohibiting a full understanding of the water balance and glacier mass balance. Hence, we measured surface latent heat fluxes with an eddy covariance system on Yala Glacier (5,350 m a.s.l) in the Nepalese Himalaya to quantify the role snow sublimation plays in the water and glacier mass budget. Observations reveal that cumulative sublimation is 32 mm for a 32-day period from October to November 2016, which is high compared to observations in other regions in the world. Multiple turbulent flux parameterizations were subsequently tested against this observed sublimation. The bulk-aerodynamic method offered the best performance, and we subsequently used this method to estimate cumulative sublimation and evaporation at the location of the eddy covariance system for the 2016-2017 winter season, which is 125 and 9 mm respectively. This is equivalent to 21% of the annual snowfall. In addition, the spatial variation of total daily sublimation over Yala Glacier was simulated with the bulk-aerodynamic method for a humid and non-humid day. Required spatial fields of meteorological variables were obtained from high-resolution WRF simulations of the region in combination with field observations. The cumulative daily sublimation at the location of the eddy covariance system equals the simulated sublimation averaged over the entire glacier. Therefore, this location appears to be representative for Yala Glacier sublimation. The spatial distribution of sublimation is primarily controlled by wind speed. Close to the ridge of Yala Glacier cumulative daily sublimation is a factor 1.7 higher than at the location of the eddy covariance system, whereas it is a factor 0.8 lower at the snout of the glacier. This illustrates that the fraction of snowfall returned to the atmosphere may be much higher than 21% at wind-exposed locations. This is a considerable loss of water and illustrates the importance and need to account for sublimation in future hydrological and mass balance studies in the Himalaya.

Temperature index (TI) models are convenient for modelling glacier ablation since they require only a few input variables and rely on simple empirical relations. The approach is generally assumed to be reliable at lower elevations (below 3500 m above sea level, a.s.l) where air temperature ( T a ) relates well to the energy inputs driving melt. We question this approach in High Mountain Asia (HMA). We study in-situ meteorological drivers of glacial ablation at two sites in central Nepal, between 2013 and 2017, using data from six automatic weather stations (AWS). During the monsoon, surface melt dominates ablation processes at lower elevations (between 4950 and 5380 m a.s.l.). As net shortwave radiation ( SW net ) is the main energy input at the glacier surface, albedo ( α ) and cloudiness play key roles while being highly variable in space and time. For these cases only, ablation can be calculated with a TI model, or with an Enhanced TI (ETI) model that includes a shortwave radiation ( SW ) scheme and site specific ablation factors. In the ablation zone during other seasons and during all seasons in the accumulation zone, sublimation and other wind-driven ablation processes also contribute to mass loss, and remain unresolved with TI or ETI methods.