Quaternary geology and Paleoclimate

Reconstructing the Englacial Stratigraphy of the Greenland Ice Sheet

PhD candidate: Philipp Immanuel Voigt

snowy landscape
Philipp Voigt

Main content

I am a PhD student at the Department of Earth Science and the Bjerknes Centre for Climate Research, in the group of Andreas Born. I have a Bachelor and Master degree in Physics from Heidelberg University, Germany, mainly focused on environmental physics, image processing and numerical simulations. Both my theses were conducted in the field of small-scale air-sea interaction. In my PhD I will work with inverse modelling of the Greenland ice sheet, assimilating radiostratigraphy data, and ultimately reconstruct accumulation rates and patterns for the last glacial cycle. The contribution of the Greenland ice sheet to sea level change can also be constrained if the reconstruction is successful. The project is funded by the Research Council of Norway.


Briefly explained the stratigraphy of an ice sheet depends not only on accumulation rates and patterns, but also on the ice dynamics. This coupling between layer thickness and ice flow vastly complicates the reconstruction of the precipitation rates. These rates are therefore relatively poorly constrained, even for the last glacial cycle. A large amount of dated stratigraphy data obtained by air-borne ice penetrating radar during the Nata Operation IceBridge exists, however, holding the information on the accumulation fields. This project aims at deploying and advancing inverse modelling techniques for evaluating the fields from the radiostratigraphy. In essence this means that the model parameters, in this case the surface mass balance, are adjusted until the (measured) ice stratigraphy data is reproduced by the model as required by the data uncertainties. The problem will be approached by using the isochronal tracing scheme ELSA (Englacial Layer Simulation Architecture) in combination the non-steady 3D thermomechanical ice model Yelmo, together explicitly simulating the ice layers. This combined model will be forced with different realistic precomputed surface mass balance patterns, simulated using BESSI (the BErgen Snow Simulator). Thereby the sensitivity of the stratigraphy on the different patterns and accumulation rates can be determined. In turn this allows for the reconstruction of past mass balance fields from the measured stratigraphy. At the current stage, a linearized model between mass balance and layer thickness will be used for the inversion process, building upon the methodology of Alexios Theofilopoulos.