The effect of the present-day imbalance on schematic and climate forced simulations of the West Antarctic Ice Sheet collapse
Abstract. Recent observations reveal that the West Antarctic Ice Sheet is rapidly thinning, particularly at its two largest outlet glaciers, Pine Island Glacier and Thwaites Glacier, while East Antarctica remains relatively stable. Projections give a mixed picture, some model project mass gain by increased surface mass balance, most models project some or severe mass loss by increasing ice discharge. In this study, we explore the effect of present-day ice thickness change rates on forced future simulations of the Antarctic Ice Sheet using the Community Ice Sheet Model (CISM). We start with a series of schematic, uniform ocean temperature perturbations to probe the sensitivity of the modelled present-day imbalance to ocean warming. We then apply ocean and atmospheric forcing from seven ESMs from the CMIP5 and CMIP6 ensemble to simulate the Antarctic Ice Sheet from 2015 to 2500. The schematic experiments suggest the presence of an ice-dynamical limit, TG cannot collapse before ~2100 without more than 2 degrees of schematic, suddcen and uniform ocean warming. Meanwhile, the maximum GMSL rise rate during the collapse increases linearly with ocean temperature, indicating that while earlier collapse timing shows diminishing returns, the rate of sea-level rise keeps on intensifying with stronger forcing. In the simulations driven with ESM forcing, including or excluding the present-day imbalance contributes for the West Antarctic Ice Sheet as much to the uncertainty in the mass loss rates in the coming 5 centuries as the choice of ESM forcing. For the East Antarctic Ice Sheet on shorter timescales (until 2100), adding the present-day observed mass change rates doubles its global mean sea level rise contribution. On longer timescales (2100–2500), the effect of the present-day observed mass change rates is smaller. Thinning of the West Antarctic Ice Sheet induced by the present-day imbalance is to a small degree partly compensated by present-day ice sheet thickening of the East Antarctic Ice Sheet over the coming centuries, which persists in our simulations. Moreover, these deviations are overshadowed by the mass losses induced by the projected ocean warming. The relative importance of including the observed present-day mass loss rates decreases for larger (ocean) warming under climate forcing, and decreases over time.
The paper presents a series of model experiments on the evolution of Thwaites Glacier over the next centuries using an ice sheet model forced by several climate models. The paper demonstrates that the historical imbalance of the glacier matters a lot for its future stability and its potential for 'collapse'. Thee authors also put forward a limit in global temperature increase for which the glacier could 'collapse'.
The paper is rather lengthy and could benefit from some trimming, which would make the message clearer. Especially the experiments of steady state versus transient initialization are of interest, followed by the forcing experiments. The introduction on the different ways of initializing models could be shortened, as the importance for the paper is to make the distinction between steady state and including imbalance.
Overall, I find this an interesting study that with some polishing and a few clarifications (see below) I would recommend for publication.
Line 14: model -> models
Line 14: what models are meant here. I guess climate models and not ice sheet models. Please specify.
Line 19: Collapse occurs 58 times in the text, but it is never defined what is exactly meant by collapse of the ice sheet. Later on (onset of collapse' is used, which also requires a clarification.
The introduction is quite long giving a complete overview of different methods of initialization. It is quite interesting in itself, but is not necessarily guiding the reader towards the core of the paper, i.e., that starting a historical simulation from an observed imbalance results in different response of TG compared to starting from steady-state conditions. It is not so much the way an initialization is done, but what the imbalance is that counts for understanding the remaining of the manuscript.
Line 91: Our null hypothesis is that the GMSL rise from the present-day mass loss rates is independent of the GMSL rise caused by an increase in ocean thermal forcing, i.e. that the present-day mass loss rates do not influence future forced projections.
Quite confusing. I would suggest to remove the mention to GMSL. It is about mass loss either caused by a given imbalance due to a grounding line retreat some time ago, or due to the current applied ocean forcing. I don't see how GMSL rise can be caused from thermal forcing (except thermal expansion, but that is not what you are talking about I presume).
Line 100: Is a spatial resolution of 4km sufficient to guarantee grounding line migration (see for instance discussion in Pattyn et al., 2013). Maybe briefly state what is done to facilitate grounding line migration at such spatial resolution.
Line 109: The regularized Coulomb friction law was already used in Joughin et al (2019), and is based on the work of Schoof and Gagliardini. (Joughin, I., Smith, B. E., and Schoof, C. G.: Regularized Coulomb Friction Laws for Ice Sheet Sliding: Application to Pine Island Glacier, Antarctica, Geophys. Res. Lett., 46, 4764–4771, https://doi.org/10.1029/2019gl082526, 2019.)
Line 121: The reference that marine sediments are likely more prevalent in submarine basins may be a bit outdated. There are more recent studies that have investigated the probability of sediment versus hard bed of Antarctica. See for instance: https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2021RG000767 and https://www.nature.com/articles/s41561-022-00992-5
It shows a larger diversity of possible outcomes for regions lying below sea level.
Line 127: What are these parameters (Ho, tau, r L)? How do they influence your optimization? It is not defined what these parameters are about.
Line 158: See my remark of Line 100: is this the way grounding line migration is dealt with? Interpolation of friction within partially floating cells AND subshelf melt as well? It has been shown in Seroussi and Morlighem (https://tc.copernicus.org/articles/12/3085/2018/) that it increases the sensitivity of grounding line retreat big time. Some discussion is needed.
Line 165: isn't this not too much different than keeping the calving front fixed, as you probably need quite high melt rates to have the front retreating through melting alone.
Line 181: Is an initialization of 10 ka enough for the temperature field to reach an equilibrium?
Table 1: Overall, I found the figures in the supplementary material more of interest than the first few figures shown in the manuscript. Therefore, some information of figures S4 and S5 could be transferred to the main manuscript and replace table 1. One way of representing this is as in Martin et al, (2011) Figure 15 (https://tc.copernicus.org/articles/5/727/2011/tc-5-727-2011.pdf), so that different regions of the ice sheet/ice shelf system are represented.
Figure 3: Is not showing integrated mass loss, but mass contribution to SL in terms of GSLR and % of VAF. Mass loss also comprises that mass that is lying underneath floatation level.
Line 307: I don't think that delta T can be considered an inverted parameter, as there is not inversion method used. Maybe use 'optimized'.
Figure 4 and Line 349: instead of pointing the readers to a supplementary figure S6 just to find out where a little line is drawn, it would be more informative to mention in the caption where this line is situated as a function of present-day GL position (i.e., XX km inland from the current GL position). This line is also defined as bedrock ridge and important as onset of collapse. What is meant by onset of collapse (see also remark on collapse in general)?
Line 400: I wouldn't call these simulations outliers. They are valid solutions for that given forcing. Just that these forcings are relatively low in melt and high in accumulation and therefore result in less mass loss than other forcings. This is not the definition of an outlier.