Search

Recent progress has been made in quantifying snowmelt in the Himalaya. Although the conditions are favorable for refreezing, little is known about the spatial variability of meltwater refreezing, hindering a complete understanding of seasonal snowmelt dynamics. This study aims to improve our understanding about how refreezing varies in space and time. We simulated refreezing with the seNorge (v2.0) snow model for the Langtang catchment, Nepalese Himalaya, covering a 5-year period. Meteorological forcing data were derived from a unique elaborate network of meteorological stations and high-resolution meteorological simulations. The results show that the annual catchment average refreezing amounts to 122 mm w.e. (21% of the melt), and varies strongly in space depending on elevation and aspect. In addition, there is a seasonal altitudinal variability related to air temperature and snow depth, with most refreezing during the early melt season. Substantial intra-annual variability resulted from fluctuations in snowfall. Daily refreezing simulations decreased by 84% (annual catchment average of 19 mm w.e.) compared to hourly simulations, emphasizing the importance of using sub-daily time steps to capture melt–refreeze cycles. Climate sensitivity experiments revealed that refreezing is highly sensitive to changes in air temperature as a 2°C increase leads to a refreezing decrease of 35%.

Terrain challenges the prediction of near-surface atmospheric conditions, even in kilometre-scale numerical weather prediction (NWP) models. In this study, the ALADIN-HIRLAM NWP system with 0.5 km horizontal grid spacing and an increased number of vertical levels is compared to the 2.5-km model system similar to the currently operational NWP system at the Norwegian Meteorological Institute. The impact of the increased resolution on the forecasts’ ability to represent boundary-layer processes is investigated for the period from 12 to 16 February 2018 in an Arctic fjord-valley system in the Svalbard archipelago. Model simulations are compared to a wide range of observations conducted during a field campaign. The model configuration with sub-kilometre grid spacing improves both the spatial structure and overall verification scores for the near-surface temperature and wind forecasts compared to the 2.5-km experiment. The sub-kilometre experiment successfully captures the wind channelling through the valley and the temperature field associated with it. In a situation of a cold-air pool development, the sub-kilometre experiment has a particularly high near-surface temperature bias at low elevations. The use of measurement campaign data, however, reveals some encouraging results, e.g. the sub-kilometre system has a more realistic vertical profile of temperature and wind speed, and the surface temperature sensitivity to the net surface energy is closer to the observations. This work demonstrates the potential of sub-kilometre NWP systems for forecasting weather in complex Arctic terrain, and also suggests that the increase in resolution needs to be accompanied with further development of other parts of the model system.

Seasonal snow is an important component of the Himalayan hydrological system, but a lack of observations at high altitude hampers understanding and forecasting of water availability in this region. Here, we use a passive gamma ray sensor that measures snow water equivalent (SWE) and complementary meteorological instruments installed at 4962 m a.s.l. in the Nepal Himalayas to quantify the evolution of SWE and snow depth over a 2-year period. We assess the accuracy, spatial representativeness and the applicability of the SWE and snow depth measurements using time-lapse camera imagery and field observations. The instrument setup performs well for snowpacks >50 mm SWE, but caution must be applied when interpreting measurements from discontinuous, patchy snow cover or those that contain lenses of refrozen meltwater. Over their typical ∼6-month lifetime, snowpacks in this setting can attain up to 200 mm SWE, of which 10–15% consists of mixed precipitation and rain-on-snow events. Precipitation gauges significantly underrepresent the solid fraction of precipitation received at this elevation by almost 40% compared to the gamma ray sensor. The application of sub-daily time-lapse camera imagery can help to correctly interpret and increase the reliability and representativeness of snowfall measurements. Our monitoring approach provides high quality, continuous, near-real time information that is essential to develop snow models in this data scarce region. We recommend that a similar instrument setup be extended into remote Himalayan environments to facilitate widespread snowpack monitoring and further our understanding of the high-altitude water cycle.

Surface energy balance models are common tools to estimate melt rates of debris-covered glaciers. In the Himalayas, radiative fluxes are occasionally measured, but very limited observations of turbulent fluxes on debris-covered tongues exist to date. We present measurements collected between 26 September and 12 October 2016 from an eddy correlation system installed on the debris-covered Lirung Glacier in Nepal during the transition between monsoon and post-monsoon. Our observations suggest that surface energy losses through turbulent fluxes reduce the positive net radiative fluxes during daylight hours between 10 and 100%, and even lead to a net negative surface energy balance after noon. During clear days, turbulent flux losses increase to over 250 W m−2 mainly due to high sensible heat fluxes. During overcast days the latent heat flux dominates the turbulent losses and together they reach just above 100 W m−2. Subsequently, we validate the performance of three bulk approaches in reproducing the observations from the eddy correlation system. Large differences exist between the approaches, and accurate estimates of surface temperature, wind speed, and surface roughness are necessary for their performance to be reasonable. Moreover, the tested bulk approaches generally overestimate turbulent latent heat fluxes by a factor 3 on clear days, because the debris-covered surface dries out rapidly, while the bulk equations assume surface saturation. Improvements to bulk surface energy models should therefore include the drying process of the surface. A sensitivity analysis suggests that, in order to be useful in distributed melt models, an accurate extrapolation of wind speed, surface temperature and surface roughness in space is a prerequisite. By applying the best performing bulk model over a complete melt period, we show that turbulent fluxes reduce the available energy for melt at the debris surface by 17% even at very low wind speeds. Overall, we conclude that turbulent fluxes play an essential role in the surface energy balance of debris-covered glaciers and that it is essential to include them in melt models.

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.