Research

Glaciers play a crucial role in the Earth System: they are important water suppliers to lower-lying areas during hot and dry periods, and they are major contributors to the observed present-day sea-level rise. Glaciers can also act as a source of natural hazards and have a major touristic value. Given their societal importance, there is large scientific interest in better understanding and accurately simulating the temporal evolution of glaciers, both in the past and in the future. An important part of my research focuses on modelling the temporal evolution of glaciers. Using insights from ice flow modelling, I also try to enhance our understanding of where and why meteorites tend to concentrate at very specific locations on the Antarctic ice sheet.

Modelling glacier evolution across various spatial scales: from case studies to global glacier modelling

I started my academic journey by modelling the evolution of individual glaciers and ice caps with detailed 3-D ice flow models. More specifically, my PhD (2013-2017) entitled ‘modelling the evolution of glaciers and ice caps in a changing climate’ consisted of modelling the dynamics of a Swiss mountain glacier (Vadret da Morteratsch) and a Greenlandic ice cap (Hans Tausen Iskappe). For this, I used 3-D higher-order ice flow models, which were coupled to surface mass balance models. Additionally, I also developed a visualization toolbox (TopoZeko), and performed a statistical regression analyses on the surface mass balance of the Morteratsch glacier complex.

As a postdoctoral scientist (since October 2017), the focus of my research has more shifted toward modelling larger ensembles of glaciers, i.e. a focus on the regional- to global-scale evolution of glaciers. For this, we use techniques and models that are simpler than those typically used when modelling the evolution of single glaciers at high spatial resolution, but which include ice-flow processes nevertheless. The aim of my first two postdoctoral projects (2017-2019 @ ETH Zürich-WSL Birmensdorf and 2019-2021 @ TU Delft) was to extend the GLObal Glacier Evolution Model (GloGEM) of Huss & Hock (2015, Frontiers in Earth Science) by incorporating a dynamic ice flow component. The newly developed flowline model, GloGEMflow, was applied to the European Alps in order to project the future evolution of all glaciers under the EURO-CORDEX RCM ensemble. GloGEMflow was also used to produce results for the second phase of GlacierMIP and formed the starting point to create a response time inventory of all glaciers in the European Alps. Other applications include studying the evolution of glaciers in the European Alps under 1.0, 1.5 and 2 degree warming, modelling the future evolution of all glaciers in Iceland and Scandinavia, modelling debris cover evolution in High-Mountain Asia and the location of possible future ice-dammed lakes. GloGEMflow is now being used by various groups of scientists through various collaborations.

More recently, I have started exploring the possibility to use 3-D models to simulate the evolution of glaciers at regional- to global scales. The goal of this work is to use techniques that were previously typically used when modelling individual glaciers in detail and now use these for large-scale glacier modelling. Using 3-D techniques is particularly relevant for large glaciers that cover the underlying landscape (often referred to as ice caps), which are difficult to represent with flowline approaches such as OGGM and GloGEMflow.

Glacier hydrology and hydropower

Several collaborations have been established with hydrologists, who use our simulated future glacier evolutions as an input for modelling the future runoff in Alpine regions. This includes partnerships with scientists from WSL Birmensdorf, University of Bern, the École Polytechnique Fédérale de Lausanne (Switzerland), the TU Delft (Netherlands) and the University of Bourgogne (France). I have also participated in a study that estimated the global hydropower and water storage potential that will arise from deglaciating basins during the 21st century. Publications on these topics can be found here.

Antarctic meteorites

From December 2012 to February 2013, I had the unique chance to join a Belgian-Japanese meteorite search expedition on the Nansen ice field on the Antarctic plateau (ca. 2900 m elevation), 120 km south of the Belgian Antarctic Station (Princes Elisabeth). During this expedition hundreds of blue ice samples surface samples were collected and these revealed clear climatic signals at the surface. To constrain the age of the ice, the terrestrial age of selected meteorites (i.e. the time since they entered the atmosphere) was determined. By combining this with satellite derived surface velocities and other sources of information, we tried to better constrain the mechanism behind the Nansen blue ice trap. As of October 2020, we are now combining many different datasets in a machine learning algorithm in order to predict where to find meteorites on the Antarctic ice sheet. This work is led by Veronica Tollenaar, funded through a PhD fellowship of the Fonds de Recherche Scientifique – FNRS.

Interests

  • Alpine glaciology
  • Arctic ice caps
  • Large-scale glacier evolution modelling
  • Machine learning
  • Numerical ice flow modeling
  • Palaeoglaciology
  • Stable isotope geochemistry of ice
  • Surface mass balance modelling and statistical techniques
  • Visualisation techniques

Research Projects

  • Alpine glacier modelling

    Modelling the time evolution of the European ice mass with a novel glacier model - GloGEMflow

    Glaciers in the European Alps play an important role in the hydrological cycle, act as a source for hydroelectricity and have a large touristic importance. The future evolution of these glaciers is driven by surface mass balance and ice flow processes, of which the latter is to date not included explicitly in regional glacier projections for the Alps. In my first postdoctoral project (@ETH Zürich & WSL Birmensdorf), we modelled the future evolution of glaciers in the European Alps with GloGEMflow, an extended version of the Global Glacier Evolution Model (GloGEM), in which both surface mass balance and ice flow are explicitly accounted for. The mass balance model was calibrated with glacier-specific geodetic mass balances, and forced with high-resolution regional climate model (RCM) simulations from the EURO-CORDEX ensemble. The evolution of the total glacier volume in the coming decades was found to be relatively similar under the various representative concentrations pathways (RCP2.6, 4.5 and 8.5), with volume losses of about 47-52% in 2050 with respect to 2017. We found that under RCP2.6, the ice loss in the second part of the 21st century is relatively limited and that about one-third (36.8% ± 11.1%, multi-model mean ± 1σ) of the present-day (2017) ice volume will still be present in 2100. Under a strong warming (RCP8.5) the future evolution of the glaciers is dictated by a substantial increase in surface melt, and glaciers are projected to largely disappear by 2100 (94.4±4.4% volume loss vs. 2017). For a given RCP, differences in future changes are mainly determined by the driving global climate model (GCM), rather than by the RCM, and these differences are larger than those arising from various model parameters (e.g. flow parameters and cross-section parameterization). We find that under a limited warming, the inclusion of ice dynamics reduces the projected mass loss and that this effect increases with the glacier elevation range, implying that the inclusion of ice dynamics is likely to be important for global glacier evolution projections.

  • Antarctica

    Blue Ice Area's and Meteorites

    From December 2012 to February 2013 I had the unique chance to join a Belgian-Japanese meteorite search expedition on the Nansen ice field on the Antarctic plateau (ca. 2900 m elevation), 120 km south of the Belgian Antarctic Station (Princes Elisabeth) (Imae et al., 2015). During this expedition hundreds of blue ice samples surface samples were collected and these revealed clear climatic signals at the surface. To constrain the age of the ice, the terrestrial age of selected meteorites (i.e. the time since they entered the atmosphere) was determined. By combining this with satellite derived surface velocities and other sources of information, we aimed at better understanding the mechanism behind the Nansen blue ice trap.

    Expedition report

    Scientific paper summarising our findings

  • TopoZeko

    A MATLAB function for 3-D and 4-D topographical visualization in geosciences

    TopoZeko is a MATLAB function for plotting a variety of natural environments with a pronounced topography, such as glaciers, volcanoes and lakes in mountainous regions. This function extends existing MATLAB plotting routines and allows for high-quality 3-D landscape visualization, with a single color defining a featured surface type or with a color scale defining the magnitude of a variable. More Info

  • Hans Tausen Iskappe (Greenland)

    Holocene and future evolution of the world's northernmost ice cap

    The second part of my PhD focussd on the dynamics of world’s northernmost ice cap: the Hans Tausen Iskappe in Northern Greenland. Several mass balance and ice thickness measurements occurred in the 1990’s and an ice-core was also drilled. From this we know that the ice cap (partly) disappeared during the Holocene Thermal Maximum and that it subsequently started building up again. By modelling the ice dynamics and the mass balance we will tried to understand how this happened, how stable the present-day ice cap is and how the ice cap may evolve in a warming climate.

    Study on present-day dynamics, sensitivity to climatic forcing and future evolution: Zekollari et al. (2017, The Cryosphere)

    Holocene evolution of the ice cap and implications for the palaeoclimatic evolution of the high Arctic: Zekollari et al. (2017, Quaternary Science Reviews)

  • Vadret da Morteratsch (Switzerland)

    Ice thickness modelling and glacier dynamics

    The first part of my PhD mainly focused on the modelling the dynamics of the Morteratsch (Engadine, Switzerland). This is a medium-sized alpine glacier on which we have field program since 2001 (coordinated by P. Huybrechts, first at the Alfred Wegener Institut – Bremerhaven, since 2005 at the Vrije Universiteit Brussel). Each year we visit the glacier at the end of the ablation season (end of September –  beginning of October) to determine the annual mass balance of the glacier. Besides this we also measure the surface velocities, the ice thickness and perform Unmanned Aerial Vehicle (UAV) flights to reconstruct the surface topography (since 2014). This field data is used to calibrate and validate our ice flow and surface mass balance model.

    Ice thickness reconstruction and ice flow model calibration: Zekollari et al. (2013, Annals of Glaciology)

    Glacier evolution between 1864 and 2100: Zekollari et al. (2014, Journal of Glaciology)

    Glacier climate-geometry imbalance, response time and volume-area scaling: Zekollari and Huybrechts (2015, Annals of Glaciology)

    Statistical modelling of the surface mass balance variability of the Morteratsch glacier, Switzerland: strong control of early melting season meteorological conditions: Zekollari and Huybrechts (accepted, Journal of Glaciology)