A topographically controlled tipping point for complete Greenland ice sheet melt

The SMB calculation in CESM2 has been evaluated with contemporary observations and high-resolution modelling and includes an advanced representation of surface melt and snow–firn processes. We find a positive SMB threshold for complete GrIS melt of 230 ± 84 Gt yr−1, corresponding to a 60 % decrease in SMB and to a global mean warming of +3.4 K compared to pre-industrial CESM2–CISM2 simulated values. In our simulations, a small change in the initial SMB forcing (from 255 to 230 Gt yr−1) and global mean warming above pre-industrial levels (from +3.2 to +3.4 K) causes an abrupt change in the GrIS final volume (from 50 % mass to nearly complete deglaciation).

This nonlinear behaviour is caused by the SMB–elevation feedback, which responds to changes in surface topography due to surface melt and GIA. The GrIS tips from ∼ 50 % mass towards nearly complete melt when the impact of melt-induced surface lowering outweighs that of GIA-induced bedrock uplift and the (initially positive) SMB becomes and remains negative for at least a few thousand years. We also find that the GrIS tips towards nearly complete melt when the ice margin in the central west unpins from a coastal region with high topography and SMB.

We show that if we keep the SMB fixed (i.e. no SMB–elevation feedback) in this relatively confined region, the ice sheet retreat is halted and nearly complete GrIS melt is prevented even though the initial SMB forcing is past the threshold. Based on the minimum GrIS configuration in previous paleo-ice-sheet modelling studies, we suggest that the surface topography in the central west might have played a role in preventing larger GrIS loss during the last interglacial period ∼ 130–115 kyr BP.