Dynamic regimes of the Greenland Ice Sheet emerging from interacting melt-elevation and glacial isostatic adjustment feedbacks
- 1Earth System Analysis, Potsdam Institute for Climate Impact Research, Member of the Leibniz Association, Telegrafenberg A31, 14473 Potsdam, Germany
- 2Institute for Physics and Astronomy, University of Potsdam, Potsdam, Germany
- 3Delft University of Technology, Faculty of Civil Engineering and Geosciences, Department of Geoscience and Remote Sensing, Delft, The Netherlands
- 4Stockholm Resilience Centre, Stockholm University, Kräftriket 2B, 11419 Stockholm, Sweden
- 1Earth System Analysis, Potsdam Institute for Climate Impact Research, Member of the Leibniz Association, Telegrafenberg A31, 14473 Potsdam, Germany
- 2Institute for Physics and Astronomy, University of Potsdam, Potsdam, Germany
- 3Delft University of Technology, Faculty of Civil Engineering and Geosciences, Department of Geoscience and Remote Sensing, Delft, The Netherlands
- 4Stockholm Resilience Centre, Stockholm University, Kräftriket 2B, 11419 Stockholm, Sweden
Abstract. The stability of the Greenland Ice Sheet under global warming is governed by a number of dynamic processes and interacting feedback mechanisms in the ice sheet, atmosphere and solid Earth. Here we study the long-term effects due to the interplay of the competing melt-elevation and glacial isostatic adjustment (GIA) feedbacks for different temperature step forcing experiments with a coupled ice-sheet and solid-Earth model. Our model results show that for warming levels above 2 °C, Greenland could become essentially ice-free on the long-term, mainly as a result of surface melting and acceleration of ice flow. These ice losses can be mitigated, however, in some cases with strong GIA feedback even promoting the partial recovery of the Greenland ice volume. We further explore the full-factorial parameter space determining the relative strengths of the two feedbacks: Our findings suggest distinct dynamic regimes of the Greenland Ice Sheets on the route to destabilization under global warming – from recovery, via quasi-periodic oscillations in ice volume to ice-sheet collapse. In the recovery regime, the initial ice loss due to warming is essentially reversed within 50,000 years and the ice volume stabilizes at 61–93 % of the present-day volume. For certain combinations of temperature increase, atmospheric lapse rate and mantle viscosity, the interaction of the GIA feedback and the melt-elevation feedback leads to self-sustained, long-term oscillations in ice-sheet volume with oscillation periods of tens to hundreds of thousands of years and oscillation amplitudes between 15–70 % of present-day ice volume. This oscillatory regime reveals a possible mode of internal climatic variability in the Earth system on time scales on the order of 100,000 years that may be excited by or synchronized with orbital forcing or interact with glacial cycles and other slow modes of variability. Our findings are not meant as scenario-based near-term projections of ice losses but rather providing insight into of the feedback loops governing the "deep future" and, thus, long-term resilience of the Greenland Ice Sheet.
- Preprint
(893 KB) -
Supplement
(395 KB) - BibTeX
- EndNote
Maria Zeitz et al.
Status: final response (author comments only)
-
RC1: 'Comment on esd-2021-100', Kristin Poinar, 23 Feb 2022
Summary and general comments:
This manuscript presents a discovery of unforced, long-term fluctuations in the size of the Greenland Ice Sheet. The fluctuations (which are not really oscillations, as they are not strictly regular or repeating) have periods ~80 - 300 kyr and originate from the interactions between the melt-elevation feedback (a positive feedback) and glacial isostatic adjustment (a negative feedback). This has not been previously studied on long (ice age) timescales in the absence of external triggers (e.g., Heinrich events initiated by ocean heat pulses) for a land-terminating ice sheet. The finding of these emergent cycles could be relevant for "deep future" states of the Greenland Ice Sheet, although it is a challenge to imagine a future without climate forcings that would presumably overshadow the internal variability. Regardless, it is an interesting discovery that merits reporting, and this paper is largely successful. I have only minor suggestions, and although they are somewhat numerous, they are all quite attainable.
Specific comments:The authors used a "power spectrum analysis" to identify periods in the ice volume time series. These methods should be explained, if only briefly, and some test for significance should be carried out. The authors state that "The oscillation times do not seem to show a clear dependence on the values for warming, lapse rate or mantle viscosity" (P11 L12). This seems troubling -- wouldn't we expect a clear pattern to emerge within the parameter space? If so, the authors should do additional thinking and put forth possible explanations for the scatter. If not, that is interesting too, and the authors should elaborate on the reasons why this system is not governed regularly.
Relatedly, Figure 2 shows that some of the parameter combinations do, apparently, have quite regular periods (especially in Figure 2b), while others do not (such as the higher lapse rates in Figure 2a). A short presentation of the values of the periods (and which are significant) should be done. The significant period values (kyr) could even simply be written inside the cyan blocks of Figure 4.
As alluded to in my summary, I suggest replacing "oscillation" throughout the manuscript with a similar word that does not imply regularity, such as "fluctuation" or even "variation". This is because the sequence of states does not always have a regular repeat interval.
The first six lines of the Discussion restate the results, as do lines 11-17 on this page. These are redundant to the rest of the manuscript and should be removed. The last three lines of the first paragraph describe one possible extended importance of this study, which is not actually studied or discussed, and therefore would be more appropriate in the Conclusion or Introduction.
Finally, I would suggest a different name than "recovery" for the state where the ice sheet reaches a new equilibrium size significantly smaller than its start. "Recovery" implies, to me, that the ice sheet returns to its initial state. More precise names could be "re-equilibration" or "new steady state".
Technical corrections:- P1 L5 - "Greenland could become essentially ice-free on the long-term" - I suggest stating the rough number of years found for this, instead of the vague "long-term".
- P1 L13 - "oscillation periods of tens to hundreds of thousand of years" - similarly, I suggest stating the rough number of years here. This is because your minimum period (80 kyr) is not that well described by "tens of thousand of years", so it is unintentionally misleading.
- P4 L4 - add Laurentide Ice Sheet, which is what Bassis et al. (2017) studied.
- P5 L8 - Please include a brief explanation, and/or citation, for why the enhancement factors (1 and 1.5) are different depending on which stress balance is used across the domain.
- P5 Sect 2.1 - The level of description of the ice sheet model (2.1.1) is much more general than the earth deformation model (2.1.2). The classic bending-beam PDE (Eq. 1) is included with all parameters described and values supplied, for instance, but the sliding law and till stress model used in PISM are only described in words, with no parameter values given. These should be enhanced to match the level of 2.1.2.
- P8 L12 - Missing reference (?).
- P11 L3 - specify meters global sea level rise; write 1 \times 10^{19} instead of 1e+19
- P11 L21 - typo "2astern"
- P12 L4 - I have never seen a zero-indexed "o/i/ii/iii" list before. I suggest standardizing to "i/ii/iii/iv".
- Table 1 - Specific values used for \Delta T are listed, which is helpful. Values for \Gamma and \nu, instead of their ranges, should be listed similarly.
- Figure 2 - Title of panel a is missing the "times" sign. X axis labels in kyr would make it more legible.
- Figure 5 - I suggest you outline or stipple the boxes that you classify as oscillating. As it is, the figure relies on the reader to interpret on their own which boxes show "significant difference" in color.
- Figure 6 - What is the mantle viscosity & climate change forcing (\Delta T) used here? It looks like it might be the same runs shown on Figure 2a, but that is only my guess.
- AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
-
RC2: 'Comment on esd-2021-100', Anonymous Referee #2, 06 Apr 2022
General comments
The paper documents intriguing dynamic behaviour of the Greenland ice sheet resulting from the interplay between the melt-elevation feedback and the GIA feedback. The material is generally well presented and easy to follow. By itself the results are very interesting and potentially provide a very novel insight into the longer-term internal dynamics of the coupled climate-ice sheet-bedrock system. At the same time, I am also very puzzled by the results, in particular by the self-sustained quasi-periodic oscillations the authors find for (a rather narrow range) of parameter combinations. Many Greenland ice sheet modelers have performed similar experiments already since the early 1990s by imposing a stepwise warming in very similar model setups involving quasi the same degree-day type of climate forcing and taking into account isostasy with state-of-the-art models, but none of these studies have ever found even a trace of the kind of oscillations described in the paper. This makes me conclude that indeed their oscillatory behaviour may well be an ‘artifact’ (to cite their own words) of their particular experimental design and parameter choice. In other words, their model behaviour is probably not a very robust type of behaviour, to say the least, and might be very difficult to replicate in other models. My suspicion is that their model behaviour is a result from the particular choice of the Lingle-Clarke isostatic model and will not show up for any other isostatic model, be it of the ‘ELRA’ type, or of the more sophisticated ‘self-gravitating visco-elastic earth model’ type.
I think the paper would be a very valuable addition to the literature of Greenland ice sheet dynamics, but first I would like to find out more on the robustness of the results and the specific role played by isostasy. A particular feature of the Lingle-Clarke model, and its implementation by Bueler et al. (2007) is that the relaxation time increases for wavelengths up to a few thousand km (a wavelength corresponding to the Greenland situation), which I believe is unrealistic. Full visco-elastic models show the contrary, the relaxation time decreases for a larger load. My guess is that it is exactly this specific behaviour of the LC model that is causing the oscillations. I suggest the authors make an effort to respond to this criticism by including material (figures and/or discussion) to prove or disprove this point.
Specific comments
page 2, line 5: a reference is needed to substantiate the 65/35% attribution of current ice losses of the Greenland ice sheet. As far as I am aware from comprehensive studies, the ratio is more like 50/50 for both SMB changes and ice calving changes (e.g . IMBIE team, 2020)
page 2, line 22: here, and elsewhere (page 6, line 5) it is stated that ‘to our knowledge’ their have been no previous studies coupling Greenland ice-sheet dynamics to bedrock dynamics. That is not entirely true. Le Meur and Huybrechts (1998, also in GJI in 2001) have done this for the glacial cycles, also in Zweck and Huybrechts (2005) Greenland ice sheet dynamics was included and was part of the sensitivity study.
page 4, line 11: explain what the ‘small ice cap instability’ is.
page 4, line 12: To what does ‘This’ refer?
page 4, line 17: explain why the factor 1/3 is expected.
page 6, section 2.1.2: a critical appraisal of the specific features of the LC model is in order here. A more thorough discussion of the dependence of the relaxation time on the wavelength of the load change and how this compares to other models is required here, as this may well be a crucial issue in this paper.
page 6, section 2.1.3: apparently the precipitation pattern from RACMO does not interact with climate change or ice sheet geometry as it seems to be fixed. Mention this explicitly and mention the shortcomings of such an approach.
page 7, section 2.1.3: is the rain fraction a function of the monthly mean temperature? If so, the transition temperature range between 0 and 2°C seems much too small. One would still expect rainfall during a month with a mean temperature below 0°C and snowfall for a mean temperature above 2°C. Please discuss the limitations of this approach.
page 7, line 16: it is mentioned that ice-ocean interaction is included via PICO. More information is needed here. Where is the ocean forcing coming from? At what resolution? What about water circulation in the fjords? How is oceanic forcing transferred to calving fronts? Does the model have a grounding line and attached ice shelves, and how are they treated? Does it matter to include ocean forcing for the type of experiments described here at all?
page 7, section 2.2.1, and associated figures in the supplement: it is puzzling to me that while the climatic mass balance from the model differs substantially from RACMO (Fig. S2), the simulated ice sheet domain almost exactly matches the observations (Figs. S1 and S3). Almost on view it can be seen that the ice-sheet wide average surface mass balance must be positive over the domains shown, yet there is hardly any advance of the margin for the initial state. How was the initial ice sheet constrained? What is the meaning of the row of black points (low or zero velocity) at the margin in Fig. S3? To me it is hard to believe that the initial state corresponds to a self-sustained steady-state ice sheet with a freely evolving margin, the latter of which is crucial in the experiments.
page 12 and further, section 3.2: A crucial issue is how realistic the bedrock model is. In the model only viscosities are changed to control the relaxation time scale. What about the effect of variations in flexural rigidity of the lithosphere?
page 14, figure 6: The figure is very difficult to read and understand, and should be improved. The colour saturation seems to represent time (but the caption does not say), however the pale parts of the lines are difficult to see. What is the meaning of both crosses? Lower axis: accumumlation-> accumulation. Left axis: are you sure the average level of topography has negative values? Please adapt the figure and the caption to increase readability.
Page 15, lines 18-20: it is impossible to discern on Figure 6 the clockwise or counterclockwise sense of the trajectories. Perhaps an arrow would help.
Page 16, line 31: Petrini et al. (2021) is a crucial reference to prove that the results are not an artifact of the specific experimental design. However, that is an EGU abstract, and cannot be checked. Remove the reference to Petrini et al. (2021).
Technical corrections
Page 3, line 14: solte -> solve
Page 3, line 32: sophisticates -> sophisticated
Page 3, line 33: year of Fettweis et al. publication missing
Page 4, line16: add ‘itself’ between ‘manifest’ and ‘on’
Page 8, table 1: the mantle viscosity value of 1x1**-19 cannot be right.
Page 8, line 12: there is a ‘?’ in the reference list
Page 8, line 22: remove the comma between ‘both’ and ‘the’.
Page 14: line 11: do not start a sentence with a capital after a semi-colon
Page 16, line 31: Hoever -> However
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
General comments
The paper documents intriguing dynamic behaviour of the Greenland ice sheet resulting from the interplay between the melt-elevation feedback and the GIA feedback. The material is generally well presented and easy to follow. By itself the results are very interesting and potentially provide a very novel insight into the longer-term internal dynamics of the coupled climate-ice sheet-bedrock system. At the same time, I am also very puzzled by the results, in particular by the self-sustained quasi-periodic oscillations the authors find for (a rather narrow range) of parameter combinations. Many Greenland ice sheet modelers have performed similar experiments already since the early 1990s by imposing a stepwise warming in very similar model setups involving quasi the same degree-day type of climate forcing and taking into account isostasy with state-of-the-art models, but none of these studies have ever found even a trace of the kind of oscillations described in the paper. This makes me conclude that indeed their oscillatory behaviour may well be an ‘artifact’ (to cite their own words) of their particular experimental design and parameter choice. In other words, their model behaviour is probably not a very robust type of behaviour, to say the least, and might be very difficult to replicate in other models. My suspicion is that their model behaviour is a result from the particular choice of the Lingle-Clarke isostatic model and will not show up for any other isostatic model, be it of the ‘ELRA’ type, or of the more sophisticated ‘self-gravitating visco-elastic earth model’ type.
I think the paper would be a very valuable addition to the literature of Greenland ice sheet dynamics, but first I would like to find out more on the robustness of the results and the specific role played by isostasy. A particular feature of the Lingle-Clarke model, and its implementation by Bueler et al. (2007) is that the relaxation time increases for wavelengths up to a few thousand km (a wavelength corresponding to the Greenland situation), which I believe is unrealistic. Full visco-elastic models show the contrary, the relaxation time decreases for a larger load. My guess is that it is exactly this specific behaviour of the LC model that is causing the oscillations. I suggest the authors make an effort to respond to this criticism by including material (figures and/or discussion) to prove or disprove this point.
Specific comments
page 2, line 5: a reference is needed to substantiate the 65/35% attribution of current ice losses of the Greenland ice sheet. As far as I am aware from comprehensive studies, the ratio is more like 50/50 for both SMB changes and ice calving changes (e.g . IMBIE team, 2020)
page 2, line 22: here, and elsewhere (page 6, line 5) it is stated that ‘to our knowledge’ their have been no previous studies coupling Greenland ice-sheet dynamics to bedrock dynamics. That is not entirely true. Le Meur and Huybrechts (1998, also in GJI in 2001) have done this for the glacial cycles, also in Zweck and Huybrechts (2005) Greenland ice sheet dynamics was included and was part of the sensitivity study.
page 4, line 11: explain what the ‘small ice cap instability’ is.
page 4, line 12: To what does ‘This’ refer?
page 4, line 17: explain why the factor 1/3 is expected.
page 6, section 2.1.2: a critical appraisal of the specific features of the LC model is in order here. A more thorough discussion of the dependence of the relaxation time on the wavelength of the load change and how this compares to other models is required here, as this may well be a crucial issue in this paper.
page 6, section 2.1.3: apparently the precipitation pattern from RACMO does not interact with climate change or ice sheet geometry as it seems to be fixed. Mention this explicitly and mention the shortcomings of such an approach.
page 7, section 2.1.3: is the rain fraction a function of the monthly mean temperature? If so, the transition temperature range between 0 and 2°C seems much too small. One would still expect rainfall during a month with a mean temperature below 0°C and snowfall for a mean temperature above 2°C. Please discuss the limitations of this approach.
page 7, line 16: it is mentioned that ice-ocean interaction is included via PICO. More information is needed here. Where is the ocean forcing coming from? At what resolution? What about water circulation in the fjords? How is oceanic forcing transferred to calving fronts? Does the model have a grounding line and attached ice shelves, and how are they treated? Does it matter to include ocean forcing for the type of experiments described here at all?
page 7, section 2.2.1, and associated figures in the supplement: it is puzzling to me that while the climatic mass balance from the model differs substantially from RACMO (Fig. S2), the simulated ice sheet domain almost exactly matches the observations (Figs. S1 and S3). Almost on view it can be seen that the ice-sheet wide average surface mass balance must be positive over the domains shown, yet there is hardly any advance of the margin for the initial state. How was the initial ice sheet constrained? What is the meaning of the row of black points (low or zero velocity) at the margin in Fig. S3? To me it is hard to believe that the initial state corresponds to a self-sustained steady-state ice sheet with a freely evolving margin, the latter of which is crucial in the experiments.
page 12 and further, section 3.2: A crucial issue is how realistic the bedrock model is. In the model only viscosities are changed to control the relaxation time scale. What about the effect of variations in flexural rigidity of the lithosphere?
page 14, figure 6: The figure is very difficult to read and understand, and should be improved. The colour saturation seems to represent time (but the caption does not say), however the pale parts of the lines are difficult to see. What is the meaning of both crosses? Lower axis: accumumlation-> accumulation. Left axis: are you sure the average level of topography has negative values? Please adapt the figure and the caption to increase readability.
Page 15, lines 18-20: it is impossible to discern on Figure 6 the clockwise or counterclockwise sense of the trajectories. Perhaps an arrow would help.
Page 16, line 31: Petrini et al. (2021) is a crucial reference to prove that the results are not an artifact of the specific experimental design. However, that is an EGU abstract, and cannot be checked. Remove the reference to Petrini et al. (2021).
Technical corrections
Page 3, line 14: solte -> solve
Page 3, line 32: sophisticates -> sophisticated
Page 3, line 33: year of Fettweis et al. publication missing
Page 4, line16: add ‘itself’ between ‘manifest’ and ‘on’
Page 8, table 1: the mantle viscosity value of 1x1**-19 cannot be right.
Page 8, line 12: there is a ‘?’ in the reference list
Page 8, line 22: remove the comma between ‘both’ and ‘the’.
Page 14: line 11: do not start a sentence with a capital after a semi-colon
Page 16, line 31: Hoever -> However
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
-
EC1: 'Comment on esd-2021-100', Michel Crucifix, 21 Apr 2022
Dear authors,
After a (somewhat) long wait, the two reviews are in. What you show in your paper can, in my view, be called 'oscilations' (even if not perfectly periodic), and oscillations in glacial/isostatic systems are rare, though not quite unprecedented (Oerlemans, J. Glacial cycles and ice-sheet modelling. Climate Change, 4, 353–374 (1982). https://doi.org/10.1007/BF02423468. The context was quite different, as well at the overall setup, but this old example suggests, as pointed out by reviewer 2, that assumptions invoved in the lithosphere/asthenosphere model are crucial. I would therefore invite you to consider the possibility of sensitivity experiments that would adress the question, though I would not be willing to substantially delay the publication of your study.
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
Dear authors,
After a (somewhat) long wait, the two reviews are in. What you show in your paper can, in my view, be called 'oscilations' (even if not perfectly periodic), and oscillations in glacial/isostatic systems are rare, though not quite unprecedented (Oerlemans, J. Glacial cycles and ice-sheet modelling. Climate Change, 4, 353–374 (1982). https://doi.org/10.1007/BF02423468. The context was quite different, as well at the overall setup, but this old example suggests, as pointed out by reviewer 2, that assumptions invoved in the lithosphere/asthenosphere model are crucial. I would therefore invite you to consider the possibility of sensitivity experiments that would adress the question, though I would not be willing to substantially delay the publication of your study.
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
-
EC2: 'Comment on esd-2021-100', Michel Crucifix, 21 Apr 2022
Dear authors,
After a (somewhat) long wait, the two reviews are in. What you show in your paper can, in my view, be called 'oscilations' (even if not perfectly periodic), and oscillations in glacial/isostatic systems are rare, though not quite unprecedented (Oerlemans, J. Glacial cycles and ice-sheet modelling. Climate Change, 4, 353–374 (1982). https://doi.org/10.1007/BF02423468. The context was quite different, as well at the overall setup, but this old example suggests, as pointed out by reviewer 2, that assumptions invoved in the lithosphere/asthenosphere model are crucial. I would therefore invite you to consider the possibility of sensitivity experiments that would adress the question, though I would not be willing to substantially delay the publication of your study.
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
Dear authors,
After a (somewhat) long wait, the two reviews are in. What you show in your paper can, in my view, be called 'oscilations' (even if not perfectly periodic), and oscillations in glacial/isostatic systems are rare, though not quite unprecedented (Oerlemans, J. Glacial cycles and ice-sheet modelling. Climate Change, 4, 353–374 (1982). https://doi.org/10.1007/BF02423468. The context was quite different, as well at the overall setup, but this old example suggests, as pointed out by reviewer 2, that assumptions invoved in the lithosphere/asthenosphere model are crucial. I would therefore invite you to consider the possibility of sensitivity experiments that would adress the question, though I would not be willing to substantially delay the publication of your study.
-
AC2: 'Reply on RC1, RC2 and EC1&EC2', Maria Zeitz, 04 May 2022
Maria Zeitz et al.
Model code and software
Parallel Ice Sheet Model (PISM) the PISM authors https://github.com/pism/pism
Maria Zeitz et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
779 | 146 | 20 | 945 | 44 | 10 | 4 |
- HTML: 779
- PDF: 146
- XML: 20
- Total: 945
- Supplement: 44
- BibTeX: 10
- EndNote: 4
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1