scholarly journals A Simplified Three-Dimensional Ice-Sheet Model Including Ice Shelves

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.

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.


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.


2019 ◽  
Vol 13 (3) ◽  
pp. 1043-1049 ◽  
Author(s):  
Tyler Pelle ◽  
Mathieu Morlighem ◽  
Johannes H. Bondzio

Abstract. Basal melting at the bottom of Antarctic ice shelves is a major control on glacier dynamics, as it modulates the amount of buttressing that floating ice shelves exert onto the ice streams feeding them. Three-dimensional ocean circulation numerical models provide reliable estimates of basal melt rates but remain too computationally expensive for century-scale projections. Ice sheet modelers therefore routinely rely on simplified parameterizations based on either ice shelf depth or more sophisticated box models. However, existing parameterizations do not accurately resolve the complex spatial patterns of sub-shelf melt rates that have been observed over Antarctica's ice shelves, especially in the vicinity of the grounding line, where basal melting is one of the primary drivers of grounding line migration. In this study, we couple the Potsdam Ice-shelf Cavity mOdel (PICO, Reese et al., 2018) to a buoyant plume melt rate parameterization (Lazeroms et al., 2018) to create PICOP, a novel basal melt rate parameterization that is easy to implement in transient ice sheet numerical models and produces a melt rate field that is in excellent agreement with the spatial distribution and magnitude of observations for several ocean basins. We test PICOP on the Amundsen Sea sector of West Antarctica, Totten, and Moscow University ice shelves in East Antarctica and the Filchner-Ronne Ice Shelf and compare the results to PICO. We find that PICOP is able to reproduce inferred high melt rates beneath Pine Island, Thwaites, and Totten glaciers (on the order of 100 m yr−1) and removes the “banding” pattern observed in melt rates produced by PICO over the Filchner-Ronne Ice Shelf. PICOP resolves many of the issues contemporary basal melt rate parameterizations face and is therefore a valuable tool for those looking to make future projections of Antarctic glaciers.


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.


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.


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

Abstract. Basal ice-shelf melting is the key driver of Antarctica's increasing sea-level contribution. In diminishing the buttressing force of the ice shelves that fringe the ice sheet the melting increases the solid-ice discharge into the ocean. Here we contrast the influence of basal melting in two different ice-shelf regions on the time-dependent response of an idealized, inherently buttressed ice-sheet-shelf system. Carrying out three-dimensional numerical simulations, the basal-melt perturbations are applied close to the grounding line in the ice-shelf's 1) ice-stream region, where the ice shelf is fed by the fastest ice masses that stream through the upstream bed trough and 2) shear margins, where the ice flow is slower. The results show that melting below one or both of the shear margins can cause a decadal to centennial increase in ice discharge that is more than twice as large compared to a similar perturbation in the ice-stream region. We attribute this to the fact that melt-induced ice-shelf thinning in the central grounding-line region is attenuated very effectively by the fast flow of the central ice stream. In contrast, the much slower ice dynamics in the lateral shear margins of the ice shelf facilitate sustained ice-shelf thinning and thereby foster buttressing reduction. Regardless of the melt location, a higher melt concentration toward the grounding line generally goes along with a stronger response. Our results highlight the vulnerability of outlet glaciers to basal melting in stagnant, buttressing-relevant ice-shelf regions, a mechanism that may gain importance under future global warming.


2018 ◽  
Author(s):  
Tyler Pelle ◽  
Mathieu Morlighem ◽  
Johannes H. Bondzio

Abstract. Basal melt at the bottom of Antarctic ice shelves is a major control on glacier dynamics, as it modulates the amount of buttressing that floating ice shelves exert onto the ice streams feeding them. Three-dimensional ocean circulation numerical models provide reliable estimates of basal melt rates but remain too computationally expensive for century scale projections. Ice sheet modelers therefore routinely rely on simplified parameterizations either based on ice shelf depth or on more sophisticated box models. However, existing parameterizations do not accurately resolve the complex spatial patterns of sub-shelf melt rates that have been observed over Antarctica's ice shelves, especially in the vicinity of the grounding line, where basal melt is one of the primary drivers of grounding line migration. In this study, we couple the Potsdam Ice-shelf Cavity mOdel (PICO) to a buoyant Plume melt rate parameterization to create PICOP, a novel basal melt rate parameterization that is easy to implement in transient ice sheet numerical models and produces a melt rate field that is in excellent agreement with the spatial distribution and magnitude of observations for a wide variety of ocean basins. We test PICOP on the Amundsen Sea sector of West Antarctica, Totten and Moscow University ice shelves in Eastern Antarctica, and the Ronne-Filchner ice shelf and compare the results to PICO. We find that PICOP is able to reproduce the high melt rates near the grounding lines of Pine Island, Thwaites, and Totten glaciers (on the order of 100 m/yr) and removes the “banding” pattern observed in melt rates produced by PICO over the Ronne-Filchner ice shelf. PICOP resolves many of the issues contemporary basal melt rate parameterizations face and is therefore a valuable tool for those looking to make future projections of Antarctic glaciers.


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.


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>


Sign in / Sign up

Export Citation Format

Share Document