Mathematical Foundation of Ice Sheet and Ice Shelf Dynamics. A Physicist’s View

2019 ◽  
pp. 192-203 ◽  
Author(s):  
K. Hutter
2012 ◽  
Vol 58 (208) ◽  
pp. 216-222 ◽  
Author(s):  
Wenxiao Pan ◽  
Alexandre M. Tartakovsky ◽  
Joe J. Monaghan

AbstractMathematical modeling of ice sheets is complicated by the nonlinearity of the governing equations and boundary conditions. Standard grid-based methods require complex front-tracking techniques and have a limited capability to handle large material deformations and abrupt changes in bottom topography. Consequently, numerical methods are usually restricted to shallow ice-sheet and ice-shelf approximations. We propose a new smoothed-particle hydrodynamics (SPH) model for coupled ice-sheet and ice-shelf dynamics. SPH, a fully Lagrangian particle method, is highly scalable and its Lagrangian nature and meshless discretization are well suited to the simulation of free surface flows, large material deformation and material fragmentation. In this paper, we use the SPH model to study ice-sheet/ice-shelf behavior, and the dynamics of the grounding line. The steady-state position of the grounding line obtained from SPH simulations is in good agreement with laboratory observations for a wide range of simulated bedrock slopes and density ratios, similar to those of ice and sea water. The numerical accuracy of the SPH algorithm is verified by simulating the plane-shear flow of two immiscible fluids and the propagation of a highly viscous blob of fluid along a horizontal surface. In the experiment, the ice was represented with a viscous Newtonian fluid. For consistency, in the described SPH model the ice is also modeled as a viscous Newtonian fluid. Typically, ice sheets are modeled as a non-Newtonian fluid, accounting for the changes in the mechanical properties of the ice. Implementation of a non-Newtonian rheology in the SPH model is the subject of our ongoing research.


2001 ◽  
Vol 47 (157) ◽  
pp. 271-282 ◽  
Author(s):  
Richard C.A. Hindmarsh ◽  
E. Le Meur

AbstractMarine ice sheets with mechanics described by the shallow-ice approximation by definition do not couple mechanically with the shelf. Such ice sheets are known to have neutral equilibria. We consider the implications of this for their dynamics and in particular for mechanisms which promote marine ice-sheet retreat. The removal of ice-shelf buttressing leading to enhanced flow in grounded ice is discounted as a significant influence on mechanical grounds. Sea-level rise leading to reduced effective pressures under ice streams is shown to be a feasible mechanism for producing postglacial West Antarctic ice-sheet retreat but is inconsistent with borehole evidence. Warming thins the ice sheet by reducing the average viscosity but does not lead to grounding-line retreat. Internal oscillations either specified or generated via a MacAyeal–Payne thermal mechanism promote migration. This is a noise-induced drift phenomenon stemming from the neutral equilibrium property of marine ice sheets. This migration occurs at quite slow rates, but these are sufficiently large to have possibly played a role in the dynamics of the West Antarctic ice sheet after the glacial maximum. Numerical experiments suggest that it is generally true that while significant changes in thickness can be caused by spatially uniform changes, spatial variability coupled with dynamical variability is needed to cause margin movement.


1998 ◽  
Vol 27 ◽  
pp. 275-280 ◽  
Author(s):  
Akira Nishimura ◽  
Toru Nakasone ◽  
Chikara Hiramatsu ◽  
Manabu Tanahashi

Based on sedimenlological and micropaleontological work on three sediment cores collected at about 167° Ε in the western Ross Sea, Antarctica, and accelerator mass spectrometer l4C ages of organic carbon, we have reconstructed environmental changes in the area during the late Quaternary. Since 38 ka BP at latest, this area was a marine environment with low productivity. A grounded ice sheet advanced and loaded the sediments before about 30-25 ka BP. After 25 ka BP, the southernmost site (76°46'S) was covered by floating ice (shelf ice), preventing deposition of coarse terrigenous materials and maintaining a supply of diatom tests and organic carbon until 20 ka BP. The northernmost site (74°33'S) was in a marine environment with a moderate productivity influenced by shelf ice/ice sheet after about 20 ka BP. Since about 10 ka BP, a sedimentary environment similar to the present-day one has prevailed over this area.


2021 ◽  
Author(s):  
Joanna Davies ◽  
Anders Møller Mathiasen ◽  
Kristiane Kristensen ◽  
Christof Pearce ◽  
Marit-Solveig Seidenkrantz

<p>The polar regions exhibit some of the most visible signs of climate change globally; annual mass loss from the Greenland Ice Sheet (GrIS) has quadrupled in recent decades, from 51 ± 65 Gt yr<sup>−1</sup> (1992-2001) to 211 ± 37 Gt yr<sup>−1</sup> (2002-2011). This can partly be attributed to the widespread retreat and speed-up of marine-terminating glaciers. The Zachariae Isstrøm (ZI) is an outlet glacier of the Northeast Greenland Ice Steam (NEGIS), one of the largest ice streams of the GrIS (700km), draining approximately 12% of the ice sheet interior. Observations show that the ZI began accelerating in 2000, resulting in the collapse of the floating ice shelf between 2002 and 2003. By 2014, the ice shelf extended over an area of 52km<sup>2</sup>, a 95% decrease in area since 2002, where it extended over 1040km<sup>2</sup>. Paleo-reconstructions provide an opportunity to extend observational records in order to understand the oceanic and climatic processes governing the position of the grounding zone of marine terminating glaciers and the extent of floating ice shelves. Such datasets are thus necessary if we are to constrain the impact of future climate change projections on the Arctic cryosphere.</p><p>A multi-proxy approach, involving grain size, geochemical, foraminiferal and sedimentary analysis was applied to marine sediment core DA17-NG-ST8-92G, collected offshore of the ZI, on  the Northeast Greenland Shelf. The aim was to reconstruct changes in the extent of the ZI and the palaeoceanographic conditions throughout the Early to Mid Holocene (c.a. 12,500-5,000 cal. yrs. BP). Evidence from the analysis of these datasets indicates that whilst there has been no grounded ice at the site over the last 12,500 years, the ice shelf of the ZI extended as a floating ice shelf over the site between 12,500 and 9,200 cal. yrs. BP, with the grounding line further inland from our study site. This was followed by a retreat in the ice shelf extent during the Holocene Thermal Maximum; this was likely to have been governed, in part, by basal melting driven by Atlantic Water (AW) recirculated from Svalbard or from the Arctic Ocean. Evidence from benthic foraminifera suggest that there was a shift from the dominance of AW to Polar Water at around 7,500 cal. yrs. BP, although the ice shelf did not expand again despite of this cooling of subsurface waters.</p>


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.


2017 ◽  
Author(s):  
Ralph Timmermann ◽  
Sebastian Goeller

Abstract. A Regional Antarctic and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner-Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in-situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase of ice-shelf basal melt rates.


2014 ◽  
Vol 8 (3) ◽  
pp. 1057-1068 ◽  
Author(s):  
Y. Gong ◽  
S. L. Cornford ◽  
A. J. Payne

Abstract. The interaction between the climate system and the large polar ice sheet regions is a key process in global environmental change. We carried out dynamic ice simulations of one of the largest drainage systems in East Antarctica: the Lambert Glacier–Amery Ice Shelf system, with an adaptive mesh ice sheet model. The ice sheet model is driven by surface accumulation and basal melt rates computed by the FESOM (Finite-Element Sea-Ice Ocean Model) ocean model and the RACMO2 (Regional Atmospheric Climate Model) and LMDZ4 (Laboratoire de Météorologie Dynamique Zoom) atmosphere models. The change of ice thickness and velocity in the ice shelf is mainly influenced by the basal melt distribution, but, although the ice shelf thins in most of the simulations, there is little grounding line retreat. We find that the Lambert Glacier grounding line can retreat as much as 40 km if there is sufficient thinning of the ice shelf south of Clemence Massif, but the ocean model does not provide sufficiently high melt rates in that region. Overall, the increased accumulation computed by the atmosphere models outweighs ice stream acceleration so that the net contribution to sea level rise is negative.


Sign in / Sign up

Export Citation Format

Share Document