In this talk, I will present two different aspects of the ice flow modelling, including a theoretical part and an applied part. In the theoretical part, I will derive some "mechanical error estimators'', i.e. estimators that can measure the mechanical error between the most accurate ice flow model (Glen-Stokes) and some approximations based on shallowness assumption. To do so, I will follow residual techniques used to obtain a posteriori estimators of the numerical error in finite element methods for non-linear elliptic problems. In the applied part, I will present some simulations of the ice flow generated by the Rhone Glacier, Switzerland, during the last glacial maximum (~ 22 000 years ago), analyse the trajectories taken by erratic boulders of different origins, and compare these results to geomorphological observations. In particular, I will show that erratic boulders, whose origin is known, constitute valuable data to infer information about paleo-climate, which is the most uncertain input of any paleo ice sheet model.
- Mathematical Geoscience Seminar