A note on hysteresis in glaciology Original Research Article

Recent studies on the mechanism governing the Laurentide ice sheet oscillations of the Last Ice Age focus on the most critical effect of the basal hydraulic processes enhanced when the ice is sliding along soft deformable beds. To understand the import of this, we consider Fowler and Johnsonʹs 0-D hydrological flow model describing the sudden and rapid movements forward (surges) of a till-based 1-D ice sheet sliding on a flat soft bed. The basic idea is that the interplay between the ice sheets dynamics, the basal drainage system, and the sliding law can generate a surging behaviour. Mathematically this means that a multiple valued relationship between the ice flux and the ice thickness arises and the mass conservation equation turns out to be of multivalued type for some special values of the dimensionless parameters involved in the model. Assuming that a multiple valued ice flux law of the Fowler and Johnson type holds, we prove the existence of a weak bounded discontinuous solution to the system which becomes periodic after a suitable time.