Author Archives: Philip Kraaijenbrink

Spatial precipitation patterns resolved with atmospheric modelling

Our new open access study, led by Pleun, shows the importance of subkilometer atmospheric modelling and correct surface boundary conditions in areas with complex topography to accurately estimate catchment-scale meteorological variability.

Frequently used gridded meteorological datasets poorly represent precipitation in the Himalayas because of their relatively low spatial resolution and the associated representation of the complex topography. Dynamical downscaling using high-resolution atmospheric models may improve the accuracy and quality of the precipitation fields. Therefore, we have used the Weather Research and Forecasting (WRF) Model to determine the resolution that is required to most accurately simulate monsoon and winter precipitation, 2-m temperature, and wind fields in the Nepalese Himalayas.

Results show that a high resolution of 500 m is computationally still feasible and provides the best match with the observations, gives the most plausible spatial distribution of precipitation, and improves the quality of the wind and temperature fields. Our findings suggest that, in combination with future improvements to atmospheric models for applications in complex terrain, subkilometer grid spacing may resolve catchment-scale meteorological variability more accurately. This will improve our capabilities to study glacio-hydrological changes at catchment and larger scale. Future modeling studies of High Mountain Asia should consider subkilometer grids to accurately estimate local meteorological variability.

P.N.J. Bonekamp, E. Collier, W.W. Immerzeel (2018), The Impact of Spatial Resolution, Land Use, and Spinup Time on Resolving Spatial Precipitation Patterns in the Himalayas, Journal of Hydrometeorology, 19, 1565-1581.


Model output compared to the observations for the summer (left panels) and winter period (right panels).


The effect of spatial resolution for the simulated winter and summer period.

PhD degree for Philip Kraaijenbrink

Wednesday September 19th Philip successfully defended his PhD during the formal Utrecht University defense ceremony and he may now call himself Dr. Kraaijenbrink.

Philip has received his PhD with distinction for his thesis titled High-resolution insights into the dynamics of Himalayan debris-covered glaciers in which he improved our understanding of debris-covered glaciers by studying them in detail with unmanned aerial vehicles and modelling. The thesis is available for open access download via the Utrecht University library.

Philip was supervised by Prof. Steven de Jong, Dr. Walter Immerzeel and Dr. Joseph Shea. The external defense committee was formed by Prof. Andreas Kääb, Prof. Etienne Berthier, Prof. Koji Fujita, Prof. Michiel van den Broeke and Dr. Francesca Pellicciotti.


Prof. dr. Steven de Jong provides Philip with his degree in the formal Utrecht University tradition.

Snow sublimation on a Himalayan glacier

The eddy covariance tower on Yala Glacier after installation in October 2016 (Photo: Walter Immerzeel).

Our new study, which was led by Emmy and published in Frontiers in Earth Science, shows that snow sublimation should no longer be ignored in future hydrological and mass balance studies in the Himalaya. We assessed the importance of snow sublimation to the water and mass budget of Yala Glacier in the Langtang Valley, Nepalese Himalaya.

From a hydrological and glaciological perspective, 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 (5350 m a.s.l) to quantify the role snow sublimation plays in the water and glacier mass budget.

The observed 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. The bulk-aerodynamic method was used 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.

A combination of meteorological observations and WRF simulations were used to estimate the spatial variability in sublimation. These simulations reveal that sublimation is primarily controlled by wind speed. The daily cumulative sublimation is a factor 1.7 higher at the ridge of Yala Glacier, which is wind-exposed, compared the location of the eddy covariance system. This is a considerable loss of water and illustrates the importance and need to account for sublimation in future studies in the Himalaya.

This work quantifies surface sublimation only. However, sublimation may be enhanced under conditions with wind-induced snow transport. Therefore, future research will focus on including this component to fully assess the importance of snow sublimation in the high-altitude water cycle.


Open access article:
Stigter, E. E., Litt, M., Steiner, J. F., Bonekamp, P. N. J., Shea, J. M., Bierkens, M. F., & Immerzeel, W. W. (2018). The importance of snow sublimation on a Himalayan glacier. Frontiers in Earth Science6 (108), 1-16.


Scatter plots of meteorological variables against sublimation rate, observed at AWS Yala Glacier. The color of the data points refers to the observed wind speed. Results show that vapor pressure deficit and wind speed are the best sublimation predictors.

Mapping surface temperatures using a UAV

A new study led by Philip presents a method to map surface temperatures of a debris-covered glacier with an unmanned aerial vehicle (UAV), which has potential to study melt processes of such glaciers. It was published open access today in Frontiers of Earth Sciences.

In the paper we map surface temperatures of Lirung Glacier in three flights on a morning in May 2016. We present a methodology to georeference and process the acquired thermal imagery, and correct for emissivity and sensor bias. Derived UAV surface temperatures are compared with distributed simultaneous ground-based temperature measurements and with Landsat 8 thermal satellite imagery.

Surface temperatures vary greatly both spatially and temporally and have a large range of 50 °C over course of the morning. Statistical analysis shows that the variability is largely independent of incoming radiation and topography, and that much of the signal in surface temperature originates from variation in properties of the debris. Future research of surface melt processes can utilize this data to further unravel heterogeneous melt patterns on debris-covered glaciers.


Kraaijenbrink, P. D. A., J. M. Shea, M. Litt, J. F. Steiner, D. Treichler, I. Koch, and W. W. Immerzeel (2018)
Mapping surface temperatures on a debris-covered glacier with an unmanned aerial vehicle
Frontiers in Earth Science, 6(64), 1–19.


Surface temperature orthomosaics of the three UAV flights on 1 May 2016 (A–C; 06:45, 09:20, and 10:35) and the brightness temperature of the Landsat 8 band 10 on 2 May 2016 at 10:32.


Comparison of the average warming rate (f) over the surveyed glacier surface area with five different DEM derivatives: aspect (a,g), slope (b,h), upstream area (c,i), relative local elevation (d,j), and mean incoming shortwave radiation (e,k). The relative importance of each variable as a predictor in a random forest regression is shown in (l)


Animation of ground-based thermal imaging of an ice cliff and its surroundings, performed synchonous to the UAV flights.


Mountain Hydrology Team at AGU 2017

The AGU Fall Meeting 2017 is coming up at a fast pace and we are excited to present some of our scientific work we have been busy with during the last year.

We are also exited and proud to be chairing a session on the High Mountain Hydrosphere with posters Wednesday afternoon (session) and an interesting oral session on Thursday morning (session). The session features talks from all over the globe aimed at improving our understanding of high mountain hydrology using a wide range of techniques.


We have something to present to you on every single day of the week!


  • In the morning Pleun will show results from her WRF modelling of high-altitude precipitation (abstract), looking at the impact of different resolutions on model outcomes.
  • In the afternoon Emmy will present her estimation of the fraction of high-altitude solid precipitation that returns back to the atmosphere through sublimation (abstract).

Clouds forming over Ganja-La in Langtang Valley.

Gamma ray sensor near Yala glacier used to measure snow water equivalent and derive sublimation. Captured by one of our time-lapse cameras.


  • In the morning Jakob will present his poster on the recent Khurdopin Glacier surge, which we documented using Planet satellite imagery (abstract). To learn more about this study, have a look at our brief communication in The Cryosphere.


  • Walter will give a presentation early morning about the drivers of the future water gap in the Ganges-Indus-Brahmaputra basins (abstract).
  • Maxime will present his recent findings on the differences in turbulent fluxes over debris-covered and debris-free glaciers (abstract) during our poster session in the afternoon.
  • In the same session, Pascal will present his efforts in modelling supraglacial ice cliffs on debris-covered glaciers and estimating their contribution to the glacier mass balance (abstract).

Walter crossing Langtang Khola.

Max admiring the beauty of Langtang Glacier.


  • As part of our own session, Remco will present his findings on the influences of cropland irrigation on high-altitude precipitation (abstract).
  • Later in the morning Philip will present his recent paper on the consequences of climate warming on the glaciers in High Mountain Asia (abstract).


  • Evan will talk about seasonal deformation of Khumbu Glacier as measured from remote sensing and field observations in the morning (abstract).
  • Jakob will talk about turbulent fluxes on a debris-covered glacier just after lunch (abstract).
  • Jakob will also present some recent investigations into the terrestrial ice margin in Greenland in the afternoon (abstract).
  • On Friday afternoon, Joseph and Philip will present their work on the multi-temporal UAV monitoring of a mountain snow pack in the Canadian Rockies using optical and thermal sensors (abstract).

A curious visitor at our time-lapse camera monitoring the debris-covered Shalbachum Glacier.

Joseph pounding in markers for the UAV surveys on Fortress Ridge.

Follow us on a trip into the world of mountain hydrology and visit one of our talks or posters!

1 2 3 8