scholarly journals Modelling of a marine glacier and ice-sheet-ice-shelf transition zone based on asymptotic analysis

1996 ◽  
Vol 23 ◽  
pp. 59-67 ◽  
Author(s):  
Vladimir A. Chugunov ◽  
Alexander V. Wilchinsky

All parts of a two-dimensional, isothermal, stationary marine glacier (grounded ice sheet, ice shelf and transition zone) with constant viscosity are analysed by perturbation methods. In so doing, all zones of different flow patterns can be considered separately. Correlations between spatial scales for all parts can be expressed in terms of the typical ice-surface slope distant from the ocean, which reflects exterior conditions of the glacier’s existence. In considering the ice-sheet–ice-shelf transition zone, a small parameter characterizing the difference between ice and water densities is used. Such an analysis allows us to find boundary conditions at the grounding line for the grounded ice mass. Glacier-surface profiles are determined by numerical methods. The grounding-line position found by using the boundary conditions derived in this paper differs from that obtained by using Thomas and Bentley’s (1978) boundary conditions by about 10% of the grounded ice-stream length.

1996 ◽  
Vol 23 ◽  
pp. 59-67 ◽  
Author(s):  
Vladimir A. Chugunov ◽  
Alexander V. Wilchinsky

All parts of a two-dimensional, isothermal, stationary marine glacier (grounded ice sheet, ice shelf and transition zone) with constant viscosity are analysed by perturbation methods. In so doing, all zones of different flow patterns can be considered separately. Correlations between spatial scales for all parts can be expressed in terms of the typical ice-surface slope distant from the ocean, which reflects exterior conditions of the glacier’s existence. In considering the ice-sheet–ice-shelf transition zone, a small parameter characterizing the difference between ice and water densities is used. Such an analysis allows us to find boundary conditions at the grounding line for the grounded ice mass. Glacier-surface profiles are determined by numerical methods. The grounding-line position found by using the boundary conditions derived in this paper differs from that obtained by using Thomas and Bentley’s (1978) boundary conditions by about 10% of the grounded ice-stream length.


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.


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.


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.


2014 ◽  
Vol 5 (2) ◽  
pp. 271-293 ◽  
Author(s):  
A. Levermann ◽  
R. Winkelmann ◽  
S. Nowicki ◽  
J. L. Fastook ◽  
K. Frieler ◽  
...  

Abstract. The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02–0.14 m; 90% range: 0.0–0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04–0.21 m; 90% range: 0.01–0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66% range: 0.04–0.17 m; 90% range: 0.02–0.25 m) for RCP-2.6 and 0.15 m (66% range: 0.07–0.28 m; 90% range: 0.04–0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.


2018 ◽  
Vol 857 ◽  
pp. 648-680 ◽  
Author(s):  
Samuel S. Pegler

A long-standing open question in glaciology concerns the propensity for ice sheets that lie predominantly submerged in the ocean (marine ice sheets) to destabilise under buoyancy. This paper addresses the processes by which a buoyancy-driven mechanism for the retreat and ultimate collapse of such ice sheets – the marine ice sheet instability – is suppressed by lateral stresses acting on its floating component (the ice shelf). The key results are to demonstrate the transition between a mode of stable (easily reversible) retreat along a stable steady-state branch created by ice-shelf buttressing to tipped (almost irreversible) retreat across a critical parametric threshold. The conditions for triggering tipped retreat can be controlled by the calving position and other properties of the ice-shelf profile and can be largely independent of basal stress, in contrast to principles established from studies of unbuttressed grounding-line dynamics. The stability and recovery conditions introduced by lateral stresses are analysed by developing a method of constructing grounding-line stability (bifurcation) diagrams, which provide a rapid assessment of the steady-state positions, their natures and the conditions for secondary grounding, giving clear visualisations of global stabilisation conditions. A further result is to reveal the possibility of a third structural component of a marine ice sheet that lies intermediate to the fully grounded and floating components. The region forms an extended grounding area in which the ice sheet lies very close to flotation, and there is no clearly distinguished grounding line. The formation of this region generates an upsurge in buttressing that provides the most feasible mechanism for reversal of a tipped grounding line. The results of this paper provide conceptual insight into the phenomena controlling the stability of the West Antarctic Ice Sheet, the collapse of which has the potential to dominate future contributions to global sea-level rise.


2021 ◽  
Author(s):  
◽  
Sanne M Maas

<p>Sediment Cores collected from the shallow sub-sea floor beneath the Ross Ice Shelf at Coulman High have been analysed using sedimentological techniques to constrain the retreat history of the Last Glacial Maximum (LGM) ice sheet in the Ross Embayment, and to determine when the modern-day calving line location of the Ross Ice Shelf was established. A characteristic vertical succession of facies was identified in these cores, that can be linked to ice sheet and ice shelf extent in the Ross Embayment. The base of this succession consists of unconsolidated, clast rich muddy diamicts, and is interpreted to be deposited subglacially or in a grounding line proximal environment on account of a distinct provenance in the clast content which can only be attributed to subglacial transport from the Byrd Glacier 400 km to the south of the drill site. This is overlain by a mud with abundant clasts, similar in character to a granulated facies that has been documented previously in the Ross Sea, and is interpreted as being a characteristic grounding line lift-o facies in a sub-ice shelf setting. These glacial proximal facies pass upward into a mud, which comprises three distinctive units. i) Muds with sub-mm scale laminae resulting from traction currents occurring near the grounding line in a sub-ice shelf environment overlain by, ii) muds with sub-mm scale laminae and elevated biogenic content (diatoms and foraminifera) and sand/gravel clasts, interpreted as being deposited in open water conditions, passing up into a iii) bioturbated mud, interpreted as being deposited in sub-ice shelf environment, proximal to the calving line. The uppermost facies consists of a 20 cm thick diatom ooze with abundant clasts and pervasive bioturbation, indicative of a condensed section deposited during periodically open marine conditions. During post-LGM retreat of the ice sheet margin in western Ross Sea, and prior to the first open marine conditions at Coulman High, it is hypothesized that the grounding and calving line were in relative close proximity to each other. As the calving line became "pinned" in the Ross Island region, the grounding line likely continued its retreat toward its present day location. New corrected radiocarbon ages on the foraminifera shells in the interval of laminated muds with clasts, provide some of the first inorganic ages from the Ross Sea, and strengthen inferences from previous studies, that the first open marine conditions in the vicinity of Ross Island were 7,600 14C yr BP. While retreat of the calving line south of its present day position is implied during this period of mid-Holocene warmth prior to its re-advance, at present it is not possible to constrain the magnitude of retreat or attribute this to climate change rather than normal calving dynamics.</p>


2021 ◽  
Author(s):  
Chen Zhao ◽  
Rupert Gladstone ◽  
Ben Galton-Fenzi ◽  
David Gwyther

&lt;p&gt;The ocean-driven basal melting has important implications for the stability of ice shelves in Antarctic, which largely affects the ice sheet mass balance, ocean circulation, and subsequently global sea level rise. Due to the limited observations in the ice shelf cavities, the couple ice sheet ocean models have been playing a critical role in examining the processes governing basal melting. In this study we use the Framework for Ice Sheet-Ocean Coupling (FISOC) to couple the Elmer/Ice full-stokes ice sheet model and the Regional Ocean Modeling System (ROMS) ocean model to model ice shelf/ocean interactions for an idealised three-dimensional domain. Experiments followed the coupled ice sheet&amp;#8211;ocean experiments under the first phase of the Marine Ice Sheet&amp;#8211;Ocean Model Intercomparison Project (MISOMIP1). A periodic pattern in the simulated mean basal melting rates is found to be highly consistent with the maximum barotropic stream function and also the grounding line retreat row by row, &amp;#160;which is likely to be related with the gyre break down near the grounding line caused by some non-physical instability events from the ocean bottom. Sensitivity tests are carried out, showing that this periodic pattern is not sensitive to the choice of couple time intervals and horizontal eddy viscosities but sensitive to vertical resolution in the ocean model, the chosen critical water column thickness in the wet-dry scheme, and the tracer properties for the nudging dry cells at the ice-ocean interface boundary. Further simulations are necessary to better explain the mechanism involved in the couple ice-ocean system, which is very significant for its application on the realistic ice-ocean systems in polar regions.&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document