Interactions between stationary waves and ice sheets: Linear versus nonlinear atmospheric response

This study examines the mutual interaction between topographically-forced atmospheric stationary waves and continental-scale ice sheets using a thermomechanical ice-sheet model coupled to a linear as well as a fully-nonlinear dry atmospheric primitive equation model. The focus is on how the stationary-wave induced ablation feeds back on the ice sheet. Simulations are conducted in which an embryonal ice mass, on an idealised "North American" continent, evolves to an equilibrium ice sheet. Under the coupling to the linear atmospheric model, the equilibrium ice sheet is primarily controlled by the ratio between the wavelength of the stationary waves and the zonal continental extent. When this ratio is near two, the ice sheet has its center of mass shifted far eastward and its shape is broadly reminiscent of the Laurentide ice sheet at LGM. For wavelengths comparable to the continental extent, however, the ice margin extends far equatorward on the central continent but is displaced poleward near the eastern coast. Remarkably, the coupling to the nonlinear atmospheric model yields equilibrium ice sheets that are virtually identical to the ones obtained in uncoupled simulations, i. e. a symmetric ice sheet with a zonal southern margin. Thus, the degree of linearity of the atmospheric response should control to what extent topographically-forced stationary waves can reorganise the structure of ice sheets. If the stationary-wave response is linear, the present results suggest that spatial reconstructions of past ice sheets can provide some information on the zonal-mean atmospheric circulation that prevailed.