Evaluation of a sub-kilometre NWP system in an Arctic fjord-valley system in winter

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 subkilometre 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

Glacier ablation and temperature indexed melt models in the Nepalese 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 (Ta) 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 (SWnet) 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

The importance of snow sublimation on a Himalayan glacier

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

Glacier ablation and temperature indexed melt models in the Nepalese 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 (Ta) 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 (SWnet) 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

The importance of snow sublimation on a Himalayan glacier

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

A Model Setup for Mapping Snow Conditions in High-Mountain Himalaya

Seasonal snow cover is an important source of melt water for irrigation and hydropower production in many regions of the world, but can also be a cause of disasters, such as avalanches and floods. In the remote Himalayan environment there is a great demand for up-to-date information on the snow conditions for the purposes of planned hydropower development and disaster risk reduction initiatives. We describe and evaluate a snow mapping setup for the remote Langtang Valley in the Nepal Himalayas, which can deliver data for snow and water availability mapping all year round. The setup utilizes (1) robust and almost maintenance-free in-situ instrumentation with satellite transmission, (2) a freely available numerical snow model, and (3) estimation of model key parameters from local meteorological and snow observations as well as from freely available climatological data. Novel features in the model include the estimation of melt parameters and solid precipitation from passive gamma-radiation based snow sensor data, as well as improved parameterization and estimation of melt water refreezing (36% of total melt) within, and sublimation/evaporation (57 mm yr−1) from the snow pack. Evaluation of the model results show a reasonable fit with snow cover data from satellite images. As many of the high-mountain regions in central and eastern Nepal show high correlation (>0.8) with the estimated snow line elevation in the Langtang catchment, the results may provide a first-order approximation of the snow conditions for these areas too

Evaluation of a sub-kilometre NWP system in an Arctic fjord-valley system in winter

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 subkilometre 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