Faezeh M. Nick

Dynamics of tidewater glaciers: comparison of three models

A minimal model of a tidewater glacier based solely on mass conservation is compared with two one-dimensional numerical flowline models, one with the calving rate proportional to water depth, and the other with the flotation criterion as a boundary condition at the glacier terminus. The models were run with two simplified bed geometries and two mass-balance formulations. The models simulate the full cycle of length variations and the equilibrium states for a tidewater glacier. This study shows that the branching of the equilibrium states depends significantly on the bed geometry. The similarity between the results of the three models indicates that if there is a submarine undulation at the terminus of a tidewater glacier, any model in which the frontal ice loss is related to the water depth yields qualitatively the same non-linear behaviour. For large glaciers extending into deep water, the flotation model causes unrealistic behaviour.

A minimal model of a tidewater glacier

Abstract We propose a simple, highly parameterized model of a tidewater glacier. The mean ice thickness and the ice thickness at the glacier front are parameterized in terms of glacier length and, when the glacier is calving, water depth. We use a linear relation between calving rate and water depth. The change in glacier length is determined by the total change in the mass budget (surface balance and calving flux), but not by the details of the glacier profile and the related velocity field. We show that this may still yield relatively rapid rates of retreat for an idealized bed geometry with a smooth overdeepening. The model is able to simulate the full cycle of ice-free conditions, glacier terminus on land, tidewater glaciers terminus, and backwards. We study two cases: (i) a glacier with a specific balance (accumulation) that is spatially uniform, and (ii) a glacier in a warmer climate with the specific balance being a linear function of altitude. Equilibrium states exhibit a double branching with respect to the climatic forcing (equilibrium-line altitude). One bifurcation is related to the dependence of the calving process on the bed profile; the other bifurcation is due to the height–mass-balance feedback. We discuss the structure of the solution diagram for different values of the calving-rate parameter. The model results are similar to those of Vieli and others (2001), who combined a fairly sophisticated two-dimensional (vertical plane) numerical ice-flow model with the modified flotation criterion suggested by Van der Veen (1996). With regard to the global dynamics of a tidewater glacier, we conclude that the details of the glacier profile or velocity field are less significant than the bed profile and the relation between the water depth and the calving rate.