The contribution of the Antarctic Ice Sheet to barystatic sea-level rise could be as high as eight metres around 2300 but remains deeply uncertain. Ice sheet retreat causes bedrock uplift, which can exert a stabilising effect on the grounding line. Yet, sea-level projections exclude bedrock adjustment, use simplified Earth structures or omit the uncertainty in climate response and Earth structure. We show that the grounding line retreat is delayed by 50 to 130 years and the barystatic sea-level contribution reduced by 9–23% when the heterogeneity of the solid Earth is included in a coupled ice – bedrock model under different emission scenarios till 2500. The effect of the solid Earth feedback in ice sheet projections can be twice as large as the uncertainty due to differences between climate models. We emphasise that realistic Earth structures should be considered when projecting the Antarctic contribution to barystatic sea-level rise on centennial time scales.

Observations and models agree that the Greenland Ice Sheet (GrIS) surface mass balance (SMB) has decreased since the end of the 1990s due to an increase in meltwater runoff and that this trend will accelerate in the future. However, large uncertainties remain, partly due to different approaches for modelling the GrIS SMB, which have to weigh physical complexity or low computing time, different spatial and temporal resolutions, different forcing fields, and different ice sheet topographies and extents, which collectively make an inter-comparison difficult. Our GrIS SMB model intercomparison project (GrSMBMIP) aims to refine these uncertainties by intercomparing 13 models of four types which were forced with the same ERA-Interim reanalysis forcing fields, except for two global models. We interpolate all modelled SMB fields onto a common ice sheet mask at 1 km horizontal resolution for the period 1980-2012 and score the outputs against (1) SMB estimates from a combination of gravimetric remote sensing data from GRACE and measured ice discharge; (2) ice cores, snow pits and in situ SMB observations; and (3) remotely sensed bare ice extent from MODerate-resolution Imaging Spectroradiometer (MODIS). Spatially, the largest spread among models can be found around the margins of the ice sheet, highlighting model deficiencies in an accurate representation of the GrIS ablation zone extent and processes related to surface melt and runoff. Overall, polar regional climate models (RCMs) perform the best compared to observations, in particular for simulating precipitation patterns. However, other simpler and faster models have biases of the same order as RCMs compared with observations and therefore remain useful tools for long-term simulations or coupling with ice sheet models. Finally, it is interesting to note that the ensemble mean of the 13 models produces the best estimate of the present-day SMB relative to observations, suggesting that biases are not systematic among models and that this ensemble estimate can be used as a reference for current climate when carrying out future model developments. However, a higher density of in situ SMB observations is required, especially in the south-east accumulation zone, where the model spread can reach 2 m w.e. yr-1 due to large discrepancies in modelled snowfall accumulation.