the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Feedback mechanisms controlling Antarctic glacial cycle dynamics simulated with a coupled ice sheet–solid Earth model
Abstract. The dynamics of the ice sheets on glacial-interglacial time scales are highly controlled by interactions with the solid Earth, i.e., glacial isostatic adjustment (GIA). Particularly at marine ice sheets, competing feedback mechanisms govern the migration of the ice sheet's grounding line and hence the ice sheet stability. In this study, we run coupled ice sheet–solid Earth simulations over the last two glacial cycles. For the ice sheet dynamics we apply the Parallel Ice Sheet Model PISM and for the load response of the solid Earth we use the three-dimensional viscoelastic Earth model VILMA, which, in addition, considers the gravitationally consistent redistribution of water (the sea level equation). We decided on an offline coupling between the two model components. By convergence of trajectories of the Antarctic Ice Sheet deglaciation we determine optimal coupling time step and spatial resolution and compare patterns of inferred relative sea level change since the Last Glacial Maximum with the results from previous studies. With our coupling setup we evaluate the relevance of feedback mechanisms for the glaciation and deglaciation phases in Antarctica considering different 3D Earth structures resulting in a range of load-response time scales. For rather long time scales, in a glacial climate associated with far-field sea level low stand, we find grounding line advance up to the edge of the continental shelf mainly in West Antarctica, dominated by a self-amplifying GIA feedback, which we call the `forebulge feedback'. For the much shorter time scale of deglaciation, dominated by the Marine Ice Sheet Instability, our simulations suggest that the stabilizing GIA feedback can significantly slow-down grounding line retreat in the Ross sector, which is dominated by a very weak Earth structure (i.e. low mantle viscosity and thin lithosphere). This delaying effect prevents a Holocene grounding line retreat beyond its present-day location, which is discussed in the scientific community, supported by observational evidence at the Siple Coast and by previous model simulations. The described coupled framework, PISM-VILMA, allows for defining restart states to run multiple sensitivity simulations from. It can be easily implemented in Earth System Models (ESMs) and provides the tools to gain a better understanding of ice sheet stability on glacial time scales as well as in a warmer future climate.
- Preprint
(3912 KB) - Metadata XML
-
Supplement
(3141 KB) - BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2023-2990', Matt King, 28 Feb 2024
The authors present an important development in coupling of a 3D GIA code with the PISM ice sheet model. They focus on the period LGM to present although the model is run for >200ka before present. The authors in particular investigate dynamics of ice sheet retreat and advance in the context of evolving sea level and bedrock, noting the existence of a 'forebulge feedback'.
The paper is well written and an important advance on previous work. I expect the paper to be influential in future developments and research that couple GIA to ice sheets or indeed to a full coupled climate model. THere is a wealth of well presented figures and the datasets are promised at final acceptance and the codes are available including the coupling codes.
I have two concerns with the paper, although I think adding some modest text and discussion would be sufficient for now
1. the authors explore mantle viscosities only down to about 10^19 Pa s whereas geodetic and seismic based studies have indicated viscosities for the upper mantle of as low as 10^18 Pa s in the Amundsen Sea and northern Antarctic Peninsula (Barletta et al. 2018; Nield et al 2014). Are the authors arguing implicitly these are transient viscosities and hence higher steady state? Or something else? Regardless, some discussion is required as 1-2 orders of magnitude change in viscosity could have a big impact on the results especially with the very low thicknesses of the elastic lithosphere in these regions. would you speculate what those differences would be? Much larger forebulges much closer to the ice load I guess - but what would that do?
2. building on this, presumably the time step and the horizontal resolution findings are a function of minimum viscosity.
3. the viscosity model adopted by the authors indicates low viscosity mantle in the Siple Coast region. Nield et al 2016 suggested that the upper-most mantle in this region was >10^20 Pa s. I think some discussion is required on this point.
minor remarks
L38 the phrase from the semicolon doesn't make sense to me and i think requires a full sentence
L59 I think it appropriate to reference Gomez here also
L71 this paragraph - I think the work of Coulon et al JGR 2021 is appropriate to mention here since they had lateral-varying version of ELRA
L127 ice->Ice
L190 lacks->lags
L215 the part in brackets is not sufficiently clear to reproduce this.
FIgure 6 it would be good to have 2 further panels focused on just the last 1ka or so.
Late Holocene readvance. THis is an interesting topic of present interest. I personally would like to see present-day RSL or VLM in a figure in supplementary material but I leave that to you. I am wondering if in East Antarctica coastline we have present-day uplift or subsidence in the model (see King et al. 2022 - but it is self-indulgent and not essential)?
L300 the text says that 3dmin is described in Fig 8a,d but those are ref-max.
L395 could not the interior subsidence is consistent with increase holocene accumulation as per other studies
L417 more->steeper
L419 *the* literature
L422 *The* largest...
FIgure 11. could you add the present-day GL location please?
Matt King Feb 28 2024
Fig S1a - it is not clear what the data are 'relative' to
'mio.' abbreviation is used throughout (and in Fig S2) but is not an English abbreviation
reference
Nield et al 2016 https://academic.oup.com/gji/article/205/1/1/2594800
King et al 2022 https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2021GL097232
Citation: https://doi.org/10.5194/egusphere-2023-2990-RC1 -
RC2: 'Comment on egusphere-2023-2990', Holly Han, 29 Feb 2024
Growth or retreat of grounded ice sheets barystatically contributes to global sea level by exchanging mass between cryosphere and hydrosphere. In addition, it changes local relative sea level by causing deformation of the solid Earth and perturbing the Earth’s gravitational field and rotation vectors, imprinting spatially varying sea-level changes.
In a little over the past decade, numerous modeling studies have shown the presence and significant impact of the feedback mechanisms between ice sheets, sea level and the solid Earth on ice-sheet modeling, highlighting the importance of coupled ice sheet – solid Earth/sea level modeling, particularly for the marine-based West Antarctic Ice Sheet (WAIS), which is situated upon reverse-sloped bedrock and under threat by MICI, and thus its stability is closely linked to the ocean depth at the grounding line and configuration of bed topography underneath. Moreover, the structure of the Earth across Antarctica is laterally heterogenous; it is characterized by thinner lithosphere and lower mantle viscosity in West Antarctica (WA), which causes faster viscous response of the solid Earth to ice loading/unloading, requiring a rather short coupling interval between the ice sheet model and solid Earth/sea-level model. Thus, capturing the interactions between ice sheets, solid Earth and sea level, as well as considering the 3D structure of the Earth and having a “short” coupling interval in the coupled modeling are warranted for accurate modeling of the Antarctic Ice Sheet.
In their paper, Albrecht and colleagues develop a newly coupled ice sheet-solid Earth/sea level model using the PISM and VILMA, which are a published ice-sheet model and viscoelastic solid Earth model, respectively. They explore Antarctic Ice Sheet dynamics over the last two glacial cycles in simulation experiments with varying Earth structure profiles, number of topographic corrections for the initial topography, coupling interval, and spatial resolution of the solid Earth model. They find a combination of these parameters according to the convergence of the AIS volume configuration and explore sensitivity of AIS to different 3D Earth Structure profiles. In addition to observing the negative sea-level feedback mechanism on grounding line retreat that has been seen by other modeling groups, the authors also observe what they call the “forebulge feedback”, which helps the grounding line to advance by reducing the local ocean depth during the glacial buildup phase. The authors also explore and confirm the effect of far-field sea-level forcing from melting of the Northern Hemispheric Ice Sheet during the last deglaciation on the AIS dynamics, which has been previously verified by other studies.
The model development, simulation experiments and analysis done in this manuscript represent lots of hard work by the authors, and undoubtingly an important contribution to the ice sheet and GIA community. It is great to see independent groups around the world building the capability of modeling coupled ice sheet - solid Earth/sea level dynamics; together which will push for better representation of ice-sheet physics and better prediction of future sea level. While I support this topic material is suitable for publication in the journal The Cryosphere, I anticipate major revision needs to be done regarding several aspects including more thorough literature review and correct citation, more exhaustive discussion on outcomes and justifications on experimental choices, which might incur additional simulations.
Video supplement
Simulation of the relative sea level change rate around Antarctica with PISM-VILMA over the last 25 kyr Torsten Albrecht https://doi.org/10.5446/65479
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
525 | 104 | 15 | 644 | 30 | 15 | 11 |
- HTML: 525
- PDF: 104
- XML: 15
- Total: 644
- Supplement: 30
- BibTeX: 15
- EndNote: 11
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1