Francesca Pellicciotti

The importance of observed gradients of air temperature and precipitation for modeling runoff from a glacierized watershed in the Nepalese Himalayas

The performance of glaciohydrological models which simulate catchment response to climate variability depends to a large degree on the data used to force the models. The forcing data become increasingly important in high-elevation, glacierized catchments where the interplay between extreme topography, climate, and the cryosphere is complex. It is challenging to generate a reliable forcing data set that captures this spatial heterogeneity. In this paper, we analyze the results of a 1 year field campaign focusing on air temperature and precipitation observations in the Langtang valley in the Nepalese Himalayas. We use the observed time series to characterize both temperature lapse rates (LRs) and precipitation gradients (PGs). We study their spatial and temporal variability, and we attempt to identify possible controlling factors. We show that very clear LRs exist in the valley and that there are strong seasonal differences related to the water vapor content in the atmosphere. Results also show that the LRs are generally shallower than the commonly used environmental lapse rates. The analysis of the precipitation observations reveals that there is great variability in precipitation over short horizontal distances. A uniform valley wide PG cannot be established, and several scale-dependent mechanisms may explain our observations. We complete our analysis by showing the impact of the observed LRs and PGs on the outputs of the TOPKAPI-ETH glaciohydrological model. We conclude that LRs and PGs have a very large impact on the water balance composition and that short-term monitoring campaigns have the potential to improve model quality considerably.

Glaciers as a Proxy to Quantify the Spatial Distribution of Precipitation in the Hunza Basin

Abstract Accurate quantification of the spatial distribution of precipitation in mountain regions is crucial for assessments of water resources and for the understanding of high-altitude hydrology, yet it is one of the largest unknowns due to the lack of high-altitude observations. The Hunza basin in Pakistan contains very large glacier systems, which, given the melt, cannot persist unless precipitation (snow input) is much higher than what is observed at the meteorological stations, mostly located in mountain valleys. Several studies, therefore, suggest strong positive vertical precipitation lapse rates; in the present study, we quantify this lapse rate by using glaciers as a proxy. We assume a neutral mass balance for the glaciers for the period from 2001 to 2003, and we inversely model the precipitation lapse by balancing the total accumulation in the catchment area and the ablation over the glacier area for the 50 largest glacier systems in the Hunza basin in the Karakoram. Our results reveal a vertical precipitation lapse rate that equals 0.21 ± 0.12% m−1, with a maximum precipitation at an elevation of 5500 masl. We showed that the total annual basin precipitation (828 mm) is 260% higher than what is estimated based on interpolated observations (319 mm); this has major consequences for hydrological modeling and water resource assessments in general. Our results were validated by using previously published studies on individual glaciers as well as the water balance of the Hunza basin. The approach is more widely applicable in mountain ranges where precipitation measurements at high altitude are lacking.

Transmission of solar radiation through clouds on melting glaciers: a comparison of parameterizations and their impact on melt modelling

We explore the robustness and transferability of parameterizations of cloud radiative forcing used in glacier melt models at two sites in the Swiss Alps. We also look at the rationale behind some of the most commonly used approaches, and explore the relationship between cloud transmittance and several standard meteorological variables. The 2 m air-temperature diurnal range is the best predictor of variations in cloud transmittance. However, linear and exponential parameterizations can only explain 30-50% of the observed variance in computed cloud transmittance factors. We examine the impact of modelled cloud transmittance factors on both solar radiation and ablation rates computed with an enhanced temperature-index model. The melt model performance decreases when modelled radiation is used, the reduction being due to an underestimation of incoming solar radiation on clear-sky days. The model works well under overcast conditions. We also seek alternatives to the use of in situ ground data. However, outputs from an atmospheric model (2.2 km horizontal resolution) do not seem to provide an alternative to the parameterizations of cloud radiative forcing based on observations of air temperature at glacier automatic weather stations. Conversely, the correct definition of overcast conditions is important.

Refined energy-balance modelling of a supraglacial pond, Langtang Khola, Nepal

Abstract Supraglacial ponds on debris-covered glaciers present a mechanism of atmosphere/glacier energy transfer that is poorly studied, and only conceptually included in mass-balance studies of debris-covered glaciers. This research advances previous efforts to develop a model of mass and energy balance for supraglacial ponds by applying a free-convection approach to account for energy exchanges at the subaqueous bare-ice surfaces. We develop the model using field data from a pond on Lirung Glacier, Nepal, that was monitored during the 2013 and 2014 monsoon periods. Sensitivity testing is performed for several key parameters, and alternative melt algorithms are compared with the model. The pond acts as a significant recipient of energy for the glacier system, and actively participates in the glacier’s hydrologic system during the monsoon. Melt rates are 2-4 cm d-1 (total of 98.5 m3 over the study period) for bare ice in contact with the pond, and <1 mmd-1 (total of 10.6m3) for the saturated debris zone. The majority of absorbed atmospheric energy leaves the pond system through englacial conduits, delivering sufficient energy to melt 2612 m3 additional ice over the study period (38.4 m3 d-1). Such melting might be expected to lead to subsidence of the glacier surface. Supraglacial ponds efficiently convey atmospheric energy to the glacier’s interior and rapidly promote the downwasting process.

Changes of glaciers in the Andes of Chile and priorities for future work

Glaciers in the Andes of Chile seem to be shrinking and possibly loosing mass, but the number and types of studies conducted, constrained mainly by data availability, are not sufficient to provide a synopsis of glacier changes for the past or future or explain in an explicit way causes of the observed changes. In this paper, we provide a systematic review of changes in glaciers for the entire country, followed by a discussion of the studies that have provided evidence of such changes. We identify a missing type of work in distributed, physically-oriented modelling studies that are needed to bridge the gap between the numerous remote sensing studies and the specific, point scale works focused on process understanding. We use an advanced mass balance model applied to one of the best monitored glaciers in the region to investigate four main research issues that should be addressed in modelling studies for a sound assessment of glacier changes: 1) the use of physically-based models of glacier ablation (energy balance models) versus more empirical models (enhanced temperature index approaches); 2) the importance of the correct extrapolation of air temperature forcing on glaciers and in high elevation areas and the large uncertainty in model outputs associated with it; 3) the role played by snow gravitational redistribution; and 4) the uncertainty associated with future climate scenarios. We quantify differences in model outputs associated with each of these choices, and conclude with suggestions for future work directions.

A grid-based model of backwasting of supraglacial ice cliffs on debris-covered glaciers

Abstract Ice cliffs might be partly responsible for the high mass losses of debris-covered glaciers in the Hindu Kush-Karakoram-Himalaya region. The few existing models of cliff backwasting are point-scale models applied at few locations or assume cliffs to be planes with constant slope and aspect, a major simplification given the complex surfaces of most cliffs. We develop the first grid-based model of cliff backwasting for two cliffs on debris-covered Lirung Glacier, Nepal. The model includes an improved representation of shortwave and longwave radiation, and their interplay with the glacier topography. Shortwave radiation varies considerably across the two cliffs, mostly due to direct radiation. Diffuse radiation is the major shortwave component, as the direct component is strongly reduced by the cliffs’ aspect and slope through self-shading. Incoming longwave radiation is higher than the total incoming shortwave flux, due to radiation emitted by the surrounding terrain, which is 25% of the incoming flux. Melt is highly variable in space, suggesting that simple models provide inaccurate estimates of total melt volumes. Although only representing 0.09% of the glacier tongue area, the total melt at the two cliffs over the measurement period is 2313 and 8282 m3, 1.23% of the total melt simulated by a glacio-hydrological model for the glacier’s tongue.

On the role of subsurface heat conduction in glacier energy-balance modelling

Abstract We discuss the inclusion of the subsurface heat-conduction flux into the calculation of the energy balance and ablation at the glacier–atmosphere interface. Data from automatic weather stations are used to force an energy-balance model at several locations on alpine glaciers and at one site in the dry Andes of central Chile. The heat-conduction flux is computed using a two-layer scheme, assuming that 36% of the net shortwave radiation is absorbed by the surface layer and that the rest penetrates into the snowpack. We compare simulations conducted with and without subsurface heat flux. Results show that assuming a surface temperature of zero degrees leads to a larger overestimation of melt at the sites in the accumulation area (10.4–13.3%) than in the ablation area (0.5–2.8%), due to lower air temperatures and the presence of snow. The difference between simulations with and without heat conduction is also high at the beginning and end of the ablation season (up to 29% for the first 15 days of the season), when air temperatures are lower and snow covers the glacier surface, while they are of little importance during periods of sustained melt at all the locations investigated.

Seasonal surface velocities of a Himalayan glacier derived by automated correlation of unmanned aerial vehicle imagery

Abstract. Debris-covered glaciers play an important role in the high-altitude water cycle in the Himalaya, yet their dynamics are poorly understood, partly because of the difficult fieldwork conditions. In this study we therefore deploy an unmanned aerial vehicle (UAV) three times (May 2013, October 2013 and May 2014) over the debris-covered Lirung Glacier in Nepal. The acquired data are processed into orthomosaics and elevation models by a Structure from Motion workflow, and seasonal surface velocity is derived using frequency cross-correlation. In order to obtain optimal surface velocity products, the effects of different input data and correlator configurations are evaluated, which reveals that the orthomosaic as input paired with moderate correlator settings provides the best results. The glacier has considerable spatial and seasonal differences in surface velocity, with maximum summer and winter velocities 6 and 2.5 m a-1, respectively, in the upper part of the tongue, while the lower part is nearly stagnant. It is hypothesized that the higher velocities during summer are caused by basal sliding due to increased lubrication of the bed. We conclude that UAVs have great potential to quantify seasonal and annual variations in flow and can help to further our understanding of debris-covered glaciers.

Contrasting climate change impact on river flows from high-altitude catchments in the Himalayan and Andes Mountains

Significance Changes in the hydrology of high-altitude catchments may have major consequences for downstream water supply. Based on model projections with a higher spatiotemporal resolution and degree of process complexity than any previous intercontinental comparative study, we show that the impacts of climate change cannot be generalized. These impacts range from a high climatic sensitivity, decreasing runoff, and significant seasonal changes in the Central Andes of Chile to increasing future runoff, limited seasonal shifts, but increases in peak flows in the Nepalese Himalaya. This study constrains uncertainty about response times and mechanisms controlling glacier and runoff response to climate and sets a benchmark for process-based modeling of the climate change impact on the hydrology of high-altitude catchments. Mountain ranges are the world’s natural water towers and provide water resources for millions of people. However, their hydrological balance and possible future changes in river flow remain poorly understood because of high meteorological variability, physical inaccessibility, and the complex interplay between climate, cryosphere, and hydrological processes. Here, we use a state-of-the art glacio-hydrological model informed by data from high-altitude observations and the latest climate change scenarios to quantify the climate change impact on water resources of two contrasting catchments vulnerable to changes in the cryosphere. The two study catchments are located in the Central Andes of Chile and in the Nepalese Himalaya in close vicinity of densely populated areas. Although both sites reveal a strong decrease in glacier area, they show a remarkably different hydrological response to projected climate change. In the Juncal catchment in Chile, runoff is likely to sharply decrease in the future and the runoff seasonality is sensitive to projected climatic changes. In the Langtang catchment in Nepal, future water availability is on the rise for decades to come with limited shifts between seasons. Owing to the high spatiotemporal resolution of the simulations and process complexity included in the modeling, the response times and the mechanisms underlying the variations in glacier area and river flow can be well constrained. The projections indicate that climate change adaptation in Central Chile should focus on dealing with a reduction in water availability, whereas in Nepal preparedness for flood extremes should be the policy priority.

Suitability of a constant air temperature lapse rate over an Alpine glacier: testing the Greuell and Böhm model as an alternative

Abstract Near-surface air temperature, typically measured at a height of 2 m, is the most important control on the energy exchange and the melt rate at a snow or ice surface. It is distributed in a simplistic manner in most glacier melt models by using constant linear lapse rates, which poorly represent the actual spatial and temporal variability of air temperature. In this paper, we test a simple thermodynamic model proposed by Greuell and Böhm in 1998 as an alternative, using a new dataset of air temperature measurements from along the flowline of Haut Glacier d’Arolla, Switzerland. The unmodified model performs little better than assuming a constant linear lapse rate. When modified to allow the ratio of the boundary layer height to the bulk heat transfer coefficient to vary along the flowline, the model matches measured air temperatures better, and a further reduction of the root-mean-square error is obtained, although there is still considerable scope for improvement. The modified model is shown to perform best under conditions favourable to the development of katabatic winds – few clouds, positive ambient air temperature, limited influence of synoptic or valley winds and a long fetch – but its performance is poor under cloudy conditions.