A major impact of anthropogenic climate change is the crossing of tipping points, which may have severe consequences such as the complete mass loss of the Greenland ice sheet (GrIS). At present, the GrIS is losing mass at an accelerated rate, largely due to a steep decrease in its surface mass balance (SMB; the balance between snow accumulation and surface ablation from melt and associated runoff). Previous work on the magnitude and nature of a threshold for GrIS complete melt remains controversial. Here, we explore a potential SMB threshold for complete melt of the GrIS; the impact and interplay of surface melt and glacial isostatic adjustment (GIA) in determining this threshold; and whether the GrIS exhibits characteristics commonly associated with tipping points, such as sensitivity to external forcing. To this end, we force the Community Ice Sheet Model v.2 (CISM2) by cycling different SMB climatologies previously calculated at multiple elevation classes with the Community Earth System Model v.2 (CESM2) in a two-way coupled CESM2-CISM2 transient simulation of the global climate and GrIS under high CO₂ forcing. 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 Gtyr^-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 Gtyr^-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 kyrBP.

The Greenland ice sheet (GrIS) is now losing mass at a rate of 0.7 mm of sea level rise (SLR) per year. Here we explore future GrIS evolution and interactions with global and regional climate under high greenhouse gas forcing with the Community Earth System Model version 2.1 (CESM2.1), which includes an interactive ice sheet component (the Community Ice Sheet Model v2.1 [CISM2.1]) and an advanced energy balance-based calculation of surface melt. We run an idealized 350-year scenario in which atmospheric CO₂ concentration increases by 1% annually until reaching four times pre-industrial values at year 140, after which it is held fixed. The global mean temperature increases by 5.2 and 8.5 K by years 131–150 and 331–350, respectively. The projected GrIS contribution to global mean SLR is 107 mm by year 150 and 1,140 mm by year 350. The rate of SLR increases from 2 mm yr⁻¹ at year 150 to almost 7 mm yr⁻¹ by year 350. The accelerated mass loss is caused by rapidly increasing surface melt as the ablation area expands, with associated albedo feedback and increased sensible and latent heat fluxes. This acceleration occurs for a global warming of approximately 4.2 K with respect to pre-industrial and is in part explained by the quasi-parabolic shape of the ice sheet, which favors rapid expansion of the ablation area as it approaches the interior “plateau.”.

The Greenland Ice Sheet (GrIS) mass balance is examined with an Earth system/ice sheet model that interactively couples the GrIS to the broader Earth system. The simulation runs from 1850 to 2100, with historical and SSP5-8.5 forcing. By the mid-21st century, the cumulative GrIS contribution to global mean sea level rise (SLR) is 23 mm. During the second half of the 21st century, the surface mass balance becomes negative in all drainage basins, with an additional SLR contribution of 86 mm. The annual mean GrIS mass loss in the last two decades is 2.7-mm sea level equivalent (SLE) year⁻¹. The increased SLR contribution from the surface mass balance (3.1 mm SLE year⁻¹) is partly offset by reduced ice discharge from thinning and retreat of outlet glaciers. The southern GrIS drainage basins contribute 73% of the mass loss in mid-century but 55% by 2100, as surface runoff increases in the northern basins.