Viscoelastic modelling of grounding line migration

In this thesis we investigate the influence of tides to the dynamics of ice sheet – ice shelf systems with numerical modelling. Tides play an important role in these systems by moving ice shelves and modulating the flow velocities of ice streams even far upstream. The grounding line as the boundary between the shelf and the ice sheet plays a crucial role in the mass balance and general stability of an ice sheet. It has been observed to migrate in response to tidal forcing, but the exact mechanisms and consequences are not yet understood in detail. On short timescales, as present in tidal forcing, we need to account for the viscoelastic character of glacier ice and choose a Maxwell model as an appropriate rheological representation. We develop and implement a viscoelastic full stokes ice flow model and implement it in the finite element software COMSOL Multiphysics. Two different test setups are used to verify our flow model and show good agreement. In our model we are able to identify two processes, which control ice flow variations with tides. Uplifting of the ice shelf leads to retreat of the grounding line and therefore less area of the ice base is in contact with the bedrock. This leads to smaller basal shear stress, resulting in an increase in flow velocity. Additionally high tide causes increased normal stress at the ice – water boundary, which slows the ice flow. When forced with the S2 (12 h) and M2 (12.42 h) tidal constituents, we observe a non-linear interaction, which leads to a perturbation of the horizontal flow velocity close to the Msf (14.76 d) constituent. By not including tides and viscoelasticity into ice models we commit significant errors for the estimation of the flux across the grounding line and the resulting mass balance. For our experimental setup this error depends on the elastic parameter and we obtain a maximal error of 3.75 %. We also observe a general retreat of the grounding line due to tidal forcing. This implies that tides possibly lead to a different equilibrium of the grounding line position