scholarly journals Snowfall versus sub-shelf melt: response of an idealized 3D ice-sheet-shelf system to mass redistribution

2018 ◽  
Author(s):  
Johannes Feldmann ◽  
Ronja Reese ◽  
Ricarda Winkelmann ◽  
Anders Levermann

Abstract. Surface accumulation and sub-ice-shelf melting are key drivers for the flow dynamics of the Antarctic Ice Sheet and are most likely to change under future warming which leads to 1) higher snowfall and 2) stronger melting below ice shelves. Here we carry out conceptual simulations in which an equilibrium ice-sheet-shelf system is perturbed such that the increased sub-shelf melting is compensated by enhanced snowfall. Although the net surface mass balance of the whole system remains unchanged, the redistribution of mass leads to a dynamic response of the ice sheet due to changes in ice thickness, surface slope, ice-shelf backstress and ice discharge. In particular, we show that such forcing can lead to the counter-intuitive situation of a retreating ice sheet which gains mass, thus having a negative sea-level contribution but smaller ice-sheet extent. The ice-sheet evolution and the corresponding steady states are investigated varying relevant parameters that affect ice properties and bed geometry as well as for different magnitudes of mass redistribution. Furthermore, the ice-sheet response is analyzed with respect to the pattern of applied melting, i.e., the area over which melting is distributed and the location where it is applied. We find throughout the ensemble of simulations that after two decades, melting at the lateral ice-shelf margins induces more ice-shelf thinning, resulting in stronger grounding line retreat and transient ice discharge compared to melting adjacent to the central grounding-line section. Analyzing changes in ice-shelf backstress with respect to changes in the ice-shelf length and mean thickness, respectively, we show that a thickness change has up to four times more influence on the backstress of the ice shelf than a length change.

2018 ◽  
Vol 12 (1) ◽  
pp. 49-70 ◽  
Author(s):  
Werner M. J. Lazeroms ◽  
Adrian Jenkins ◽  
G. Hilmar Gudmundsson ◽  
Roderik S. W. van de Wal

Abstract. Basal melting below ice shelves is a major factor in mass loss from the Antarctic Ice Sheet, which can contribute significantly to possible future sea-level rise. Therefore, it is important to have an adequate description of the basal melt rates for use in ice-dynamical models. Most current ice models use rather simple parametrizations based on the local balance of heat between ice and ocean. In this work, however, we use a recently derived parametrization of the melt rates based on a buoyant meltwater plume travelling upward beneath an ice shelf. This plume parametrization combines a non-linear ocean temperature sensitivity with an inherent geometry dependence, which is mainly described by the grounding-line depth and the local slope of the ice-shelf base. For the first time, this type of parametrization is evaluated on a two-dimensional grid covering the entire Antarctic continent. In order to apply the essentially one-dimensional parametrization to realistic ice-shelf geometries, we present an algorithm that determines effective values for the grounding-line depth and basal slope in any point beneath an ice shelf. Furthermore, since detailed knowledge of temperatures and circulation patterns in the ice-shelf cavities is sparse or absent, we construct an effective ocean temperature field from observational data with the purpose of matching (area-averaged) melt rates from the model with observed present-day melt rates. Our results qualitatively replicate large-scale observed features in basal melt rates around Antarctica, not only in terms of average values, but also in terms of the spatial pattern, with high melt rates typically occurring near the grounding line. The plume parametrization and the effective temperature field presented here are therefore promising tools for future simulations of the Antarctic Ice Sheet requiring a more realistic oceanic forcing.


1990 ◽  
Vol 14 ◽  
pp. 17-19 ◽  
Author(s):  
W.J. Böhmer ◽  
K. Herterich

We present a simplified numerical three-dimensional ice-sheet/ice-shelf model with a coarse horizontal resolution (100 km), designed for simulations of ice-volume changes on ice-age time scales (100 000 years and longer). The ice-sheet part uses the shallow-ice approximation to determine the flow, and includes a three-dimensional temperature calculation. The ice shelf is described in a quasi-stationary way. Ice-shelf thickness depends only on the thicknesses at the grounding line and the distances to the grounding line. The effect of the transition zone between ice sheet and ice shelf (assuming a width ≪100 km) is parameterized in terms of the ice thicknesses defined on the coarse grid. The characteristics at the base of the transition zone formally enter through a friction coefficient μ. We performed a series of sensitivity experiments with the coupled system, by integrating over 10 000 model years, starting from the present (modelled) state of the Antarctic and forcing the model by currently-observed accumulation rates. The position of the grounding line of the ice-sheet/ice-shelf model is quite sensitive to the choice of the friction parameter μ (in the range 0.025 > μ > 0.01). With μ = 0.05, the grounding line was maintained at the currently-observed position in the model.


2020 ◽  
Author(s):  
Christoph Kittel ◽  
Charles Amory ◽  
Cécile Agosta ◽  
Nicolas Jourdain ◽  
Stefan Hofer ◽  
...  

<p><span>The surface mass balance (SMB) of the Antarctic ice sheet is often considered as a negative contributor to the sea level rise as present snowfall accumulation largely compensate</span><span>s</span><span> for ablation through wind erosion, sublimation and runoff. The latter is even almost negligible since current Antarctic surface melting is limited to relatively scarce events over generally peripheral areas and refreezes almost entirely into the snowpack. However, melting can significantly affect the stability of ice shelves through hydrofracturing, potentially leading to their disintegration, acceleration of grounded ice and increased sea level rise. Although a large increase in snowfall is expected in a warmer climate, more numerous and stronger melting events could conversely lead to a larger risk of ice shelf collapse. In this study, we provide an estimation of the SMB of the Antarctic ice sheet for the end of the 21st century by forcing the state-of-the-art regional climate model MAR with three different global climate models. We chose the models (from both the Coupled Model Intercomparison Project Phase 5 and 6 - CMIP5 and CMIP6) providing the best metrics for representing the current Antarctic climate. While the increase in snowfall largely compensates snow ablation through runoff in CMIP5-forced projections, CMIP6-forced simulations reveal that runoff cannot be neglected in the future as it accounts for a maximum of 50% of snowfall and becomes the main ablation component over the ice sheet. Furthermore, we identify a tipping point (ie., a warming of 4°C) at which the Antarctic SMB starts to decrease as a result of enhanced runoff particularly over ice shelves. Our results highlight the importance of taking into account meltwater production and runoff and indicate that previous model studies neglecting these processes yield overestimated SMB estimates, ultimately leading to underestimated Antarctic contribution to sea level rise. Finally, melt rates over each ice shelf are higher than those that led to the collapse of the Larsen A and B ice shelves, suggesting a high probability of ice shelf collapses all over peripheral Antarctica by 2100.</span></p>


2017 ◽  
Author(s):  
Werner M. J. Lazeroms ◽  
Adrian Jenkins ◽  
G. Hilmar Gudmundsson ◽  
Roderik S. W. van de Wal

Abstract. Basal melting below ice shelves is a major factor in the decline of the Antarctic ice sheet, which can contribute significantly to possible future sea-level rise. Therefore, it is important to have an adequate description of the basal melt rates for use in ice-dynamical models. Most current ice models use rather simple parametrizations based on the local balance of heat between ice and ocean. In this work, however, we use a recently derived parametrization of the melt rates based on a buoyant meltwater plume travelling upward beneath an ice shelf. This plume parametrization combines a nonlinear ocean temperature sensitivity with an inherent geometry dependence, which is mainly described by the grounding-line depth zgl and the local slope α of the ice-shelf base. For the first time, this type of parametrization is evaluated on a two-dimensional grid covering the entire Antarctic continent. In order to apply the essentially one-dimensional parametrization to realistic ice-shelf geometries, we present an algorithm that determines effective values for zgl and α for any point beneath an ice shelf. Furthermore, since detailed knowledge of temperatures and flow patterns in the ice-shelf cavities is sparse or absent, we construct an effective ocean temperature field from observational data with the purpose of matching (area-averaged) melt rates from the model with observed present-day melt rates. The result is a realistic map of basal melt rates around Antarctica, not only in terms of average values, but also in terms of the spatial pattern, with high melt rates typically occurring near the grounding line. The plume parametrization and the effective temperature field are therefore promising tools for future simulations of the Antarctic ice sheet.


1990 ◽  
Vol 14 ◽  
pp. 17-19 ◽  
Author(s):  
W.J. Böhmer ◽  
K. Herterich

We present a simplified numerical three-dimensional ice-sheet/ice-shelf model with a coarse horizontal resolution (100 km), designed for simulations of ice-volume changes on ice-age time scales (100 000 years and longer). The ice-sheet part uses the shallow-ice approximation to determine the flow, and includes a three-dimensional temperature calculation. The ice shelf is described in a quasi-stationary way. Ice-shelf thickness depends only on the thicknesses at the grounding line and the distances to the grounding line. The effect of the transition zone between ice sheet and ice shelf (assuming a width ≪100 km) is parameterized in terms of the ice thicknesses defined on the coarse grid. The characteristics at the base of the transition zone formally enter through a friction coefficient μ. We performed a series of sensitivity experiments with the coupled system, by integrating over 10 000 model years, starting from the present (modelled) state of the Antarctic and forcing the model by currently-observed accumulation rates. The position of the grounding line of the ice-sheet/ice-shelf model is quite sensitive to the choice of the friction parameter μ (in the range 0.025 > μ > 0.01). With μ = 0.05, the grounding line was maintained at the currently-observed position in the model.


2021 ◽  
Author(s):  
Sainan Sun ◽  
Frank Pattyn

<p>Mass loss of the Antarctic ice sheet contributes the largest uncertainty of future sea-level rise projections. Ice-sheet model predictions are limited by uncertainties in climate forcing and poor understanding of processes such as ice viscosity. The Antarctic BUttressing Model Intercomparison Project (ABUMIP) has investigated the 'end-member' scenario, i.e., a total and sustained removal of buttressing from all Antarctic ice shelves, which can be regarded as the upper-bound physical possible, but implausible contribution of sea-level rise due to ice-shelf loss. In this study, we add successive layers of ‘realism’ to the ABUMIP scenario by considering sustained regional ice-shelf collapse and by introducing ice-shelf regrowth after collapse with the inclusion of ice-sheet and ice-shelf damage (Sun et al., 2017). Ice shelf regrowth has the ability to stabilize grounding lines, while ice shelf damage may reinforce ice loss. In combination with uncertainties from basal sliding and ice rheology, a more realistic physical upperbound to ice loss is sought. Results are compared in the light of other proposed mechanisms, such as MICI due to ice cliff collapse.</p>


2017 ◽  
Vol 63 (240) ◽  
pp. 731-744 ◽  
Author(s):  
JORGE BERNALES ◽  
IRINA ROGOZHINA ◽  
MAIK THOMAS

ABSTRACTIce-shelf basal melting is the largest contributor to the negative mass balance of the Antarctic ice sheet. However, current implementations of ice/ocean interactions in ice-sheet models disagree with the distribution of sub-shelf melt and freezing rates revealed by recent observational studies. Here we present a novel combination of a continental-scale ice flow model and a calibration technique to derive the spatial distribution of basal melting and freezing rates for the whole Antarctic ice-shelf system. The modelled ice-sheet equilibrium state is evaluated against topographic and velocity observations. Our high-resolution (10-km spacing) simulation predicts an equilibrium ice-shelf basal mass balance of −1648.7 Gt a−1 that increases to −1917.0 Gt a−1 when the observed ice-shelf thinning rates are taken into account. Our estimates reproduce the complexity of the basal mass balance of Antarctic ice shelves, providing a reference for parameterisations of sub-shelf ocean/ice interactions in continental ice-sheet models. We perform a sensitivity analysis to assess the effects of variations in the model set-up, showing that the retrieved estimates of basal melting and freezing rates are largely insensitive to changes in the internal model parameters, but respond strongly to a reduction of model resolution and the uncertainty in the input datasets.


2018 ◽  
Author(s):  
Niall Gandy ◽  
Lauren J. Gregoire ◽  
Jeremy C. Ely ◽  
Christopher D. Clark ◽  
David M. Hodgson ◽  
...  

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. A valuable case to examine these processes is the last deglaciation of the British-Irish Ice Sheet. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British-Irish Ice Sheet during the last deglaciation, is well constrained, with abundant empirical data which could be used to inform, validate and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We simulate retreat from the shelf edge under constant "warm" surface mass balance and subshelf melt, to isolate the role of internal ice dynamics from external forcings. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at a "pinning-point" between mainland Scotland and the Isle of Lewis. At this stage, the presence of ice shelves became a major control on deglaciation, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch Strait produces an acceleration in retreat, leading to a "collapsed" state, even when the climate returns to the initial "cold" conditions. Our simulations demonstrate the importance of the Marine Ice Sheet Instability and ice shelf buttressing during the deglaciation of parts of the British-Irish Ice Sheet. Thus, geological data could be used to constrain these processes in ice sheet models used for projecting the future of our contemporary ice sheets.


2021 ◽  
Author(s):  
Erwin Lambert ◽  
André Jüling ◽  
Paul Holland ◽  
Roderik van de Wal

<p>The contact between ice shelves and relatively warm ocean waters causes basal melt, ice shelf thinning, and ultimately ice sheet mass loss. This basal melt, and its dependence on ocean properties, is poorly understood due to an overall lack of direct observations and a difficulty in explicit simulation of the circulation in sub-shelf cavities. In this study, we compare a number of parameterisations and models of increasing complexity, up to a 2D ‘Layer’ model. Each model is aimed at quantifying basal melt rates as a function of offshore temperature and salinity. We test these models in an idealised setting (ISOMIP+) and in a realistic setting for the Amundsen Sea Embayment. All models show a comparable non-linear sensitivity of ice-shelf average basal melt to ocean warming, indicating a positive feedback between melt and circulation. However, the Layer model is the only one which explicitly resolves the flow direction of the buoyant melt plumes, which is primarily governed by rotation and by the basal topography of the ice shelves. At 500m resolution, this model simulates locally enhanced basal melt near the grounding line, in topographical channels, and near the western boundary. The simulated melt patterns for the Amundsen Sea ice shelves are compared to satellite observations of ice shelf thinning and to 3D numerical simulations of the sub-shelf cavity circulation. As detailed melt rates near the grounding line are essential for the stability of ice sheets, spatially realistic melt rates are crucial for future projections of ice sheet dynamics. We conclude that the Layer model can function as a relatively cheap yet realistic model to downscale 3D ocean simulations of ocean properties to sub-kilometer scale basal melt fields to provide detailed forcing fields to ice sheet models.</p>


2019 ◽  
Vol 13 (5) ◽  
pp. 1441-1471 ◽  
Author(s):  
Hélène Seroussi ◽  
Sophie Nowicki ◽  
Erika Simon ◽  
Ayako Abe-Ouchi ◽  
Torsten Albrecht ◽  
...  

Abstract. Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.


Sign in / Sign up

Export Citation Format

Share Document