Towards tackling ice-sheet ocean interaction with Finite Element Methods

2020 ◽  
Author(s):  
Stefano Ottolenghi ◽  
Josefin Ahlkrona
Keyword(s):  
Finite Element ◽  
Ocean Circulation ◽  
Finite Element Methods ◽  
Ice Sheet ◽  
Ocean Dynamics ◽  
Ice Shelf ◽  
Computing Platform ◽  
Fjord Circulation ◽  
Careful Handling

<p><span>Ice sheet-ocean interaction is important to properly understand phenomena such as ice sheet melting and ocean circulation. While the long term goal of this project is to fully couple the ice and ocean in one single numerical framework, we here start by modelling the ocean. We use the full non-hydrostatic equations in order to accurately model the complex ocean dynamics near the ice sheets. As numerical method, we employ finite element methods due to their capability of representing a complex fjord geometry and locally refining the mesh in the areas which require more careful handling, and its strong mathematical foundation. This will allow to overcome classical problems such as representing a moving ice shelf in a discretized setting. We here present an example of modeled fjord circulation obtained simulating the model with the FEniCS computing platform. <br></span></p>

scholarly journals Antarctic sub-shelf melt rates via PICO

2017 ◽  
Author(s):  
Ronja Reese ◽  
Torsten Albrecht ◽  
Matthias Mengel ◽  
Xylar Asay-Davis ◽  
Ricarda Winkelmann
Keyword(s):  
Ocean Circulation ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Amundsen Sea ◽  
Ice Shelves ◽  
Continental Scale ◽  
Wide Range ◽  
The Antarctic ◽  
Fully Coupled

Abstract. Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of sub-shelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice-ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice-shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. The two-dimensional melt rate fields provided by the model reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 m a−1 for cold sub-shelf cavities, for example underneath Ross or Ronne ice shelves, to 12 m a−1 for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally-feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.

scholarly journals Antarctic sub-shelf melt rates via PICO

2018 ◽  
Vol 12 (6) ◽  
pp. 1969-1985 ◽  
Author(s):  
Ronja Reese ◽  
Torsten Albrecht ◽  
Matthias Mengel ◽  
Xylar Asay-Davis ◽  
Ricarda Winkelmann
Keyword(s):  
Ocean Circulation ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Amundsen Sea ◽  
Ice Shelves ◽  
Continental Scale ◽  
Wide Range ◽  
The Antarctic ◽  
Fully Coupled

Abstract. Ocean-induced melting below ice shelves is one of the dominant drivers for mass loss from the Antarctic Ice Sheet at present. An appropriate representation of sub-shelf melt rates is therefore essential for model simulations of marine-based ice sheet evolution. Continental-scale ice sheet models often rely on simple melt-parameterizations, in particular for long-term simulations, when fully coupled ice–ocean interaction becomes computationally too expensive. Such parameterizations can account for the influence of the local depth of the ice-shelf draft or its slope on melting. However, they do not capture the effect of ocean circulation underneath the ice shelf. Here we present the Potsdam Ice-shelf Cavity mOdel (PICO), which simulates the vertical overturning circulation in ice-shelf cavities and thus enables the computation of sub-shelf melt rates consistent with this circulation. PICO is based on an ocean box model that coarsely resolves ice shelf cavities and uses a boundary layer melt formulation. We implement it as a module of the Parallel Ice Sheet Model (PISM) and evaluate its performance under present-day conditions of the Southern Ocean. We identify a set of parameters that yield two-dimensional melt rate fields that qualitatively reproduce the typical pattern of comparably high melting near the grounding line and lower melting or refreezing towards the calving front. PICO captures the wide range of melt rates observed for Antarctic ice shelves, with an average of about 0.1 m a−1 for cold sub-shelf cavities, for example, underneath Ross or Ronne ice shelves, to 16 m a−1 for warm cavities such as in the Amundsen Sea region. This makes PICO a computationally feasible and more physical alternative to melt parameterizations purely based on ice draft geometry.

scholarly journals Viscosity and elasticity: a model intercomparison of ice-shelf bending in an Antarctic grounding zone

2017 ◽  
Vol 63 (240) ◽  
pp. 573-580 ◽  
Author(s):  
CHRISTIAN T. WILD ◽  
OLIVER J. MARSH ◽  
WOLFGANG RACK
Keyword(s):  
Young’S Modulus ◽  
Ocean Circulation ◽  
Young's Modulus ◽  
Viscoelastic Model ◽  
Ice Sheet ◽  
Global Ocean ◽  
Elastic Model ◽  
Ice Shelf ◽  
Bending Stresses ◽  
Global Ocean Circulation

ABSTRACTGrounding zones are vital to ice-sheet mass balance and its coupling to the global ocean circulation. Processes here determine the mass discharge from the grounded ice sheet, to the floating ice shelves. The response of this transition zone to tidal forcing has been described by both elastic and viscoelastic models. Here we examine the validity of these models for grounding zone flexure over tidal timescales using field data from the Southern McMurdo Ice Shelf (78° 15′S, 167° 7′E). Observations of tidal movement were carried out by simultaneous tiltmeter and GPS measurements along a profile across the grounding zone. Finite-element simulations covering a 64 d period reveal that the viscoelastic model fits best the observations using a Young's modulus of 1.6 GPa and a viscosity of 1013.7 Pa s (≈ 50.1 TPa s). We conclude that the elastic model is only well-constrained for tidal displacements >35% of the spring-tidal amplitude using a Young's modulus of 1.62 ± 0.69 GPa, but that a viscoelastic model is necessary to adequately capture tidal bending at amplitudes below this threshold. In grounding zones where bending stresses are greater than at the Southern McMurdo Ice Shelf or ice viscosity is lower, the threshold would be even higher.

scholarly journals An improved bathymetry compilation for the Bellingshausen Sea, Antarctica, to inform ice-sheet and ocean models

2010 ◽  
Vol 4 (4) ◽  
pp. 2079-2101 ◽  
Author(s):  
A. G. C. Graham ◽  
F. O. Nitsche ◽  
R. D. Larter
Keyword(s):  
Ocean Circulation ◽  
Ice Sheet ◽  
West Antarctica ◽  
Ice Shelf ◽  
Ice Streams ◽  
Ice Shelves ◽  
Sea Floor ◽  
Bellingshausen Sea ◽  
Warm Deep Water ◽  
Basal Melting

Abstract. The southern Bellingshausen Sea (SBS) is a rapidly-changing part of West Antarctica, where oceanic and atmospheric warming has led to the recent basal melting and break-up of the Wilkins ice shelf, the dynamic thinning of fringing glaciers, and sea-ice reduction. Accurate sea-floor morphology is vital for understanding the continued effects of each process upon changes within Antarctica's ice sheets. Here we present a new bathymetric grid for the SBS compiled from shipborne echo-sounder, spot-sounding and sub-ice measurements. The 1-km grid is the most detailed compilation for the SBS to-date, revealing large cross-shelf troughs, shallow banks, and deep inner-shelf basins that continue inland of coastal ice shelves. The troughs now serve as pathways which allow warm deep water to access the ice fronts in the SBS. Our dataset highlights areas still lacking bathymetric constraint, as well as regions for further investigation, including the likely routes of palaeo-ice streams. The new compilation is a major improvement upon previous grids and will be a key dataset for incorporating into simulations of ocean circulation, ice-sheet change and history. It will also serve forecasts of ice stability and future sea-level contributions from ice loss in West Antarctica, required for the next IPCC assessment report in 2013.

Representation of basal melting in idealised coupled ice sheet ocean models

2021 ◽  
Author(s):  
Chen Zhao ◽  
Rupert Gladstone ◽  
Ben Galton-Fenzi ◽  
David Gwyther
Keyword(s):  
Ocean Circulation ◽  
Critical Role ◽  
Ice Sheet ◽  
Ocean Model ◽  
Periodic Pattern ◽  
Regional Ocean Modeling System ◽  
Ice Shelf ◽  
Grounding Line ◽  
Ocean Models ◽  
Basal Melting

<p>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–ocean experiments under the first phase of the Marine Ice Sheet–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,  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.</p>

scholarly journals Modelling the long-term response of the Antarctic ice sheet to global warming

1998 ◽  
Vol 27 ◽  
pp. 161-168 ◽  
Author(s):  
Roland C. Warner ◽  
W.Κ. Budd
Keyword(s):  
Global Warming ◽  
Mass Loss ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Antarctic Ice Sheet ◽  
Ice Shelves ◽  
Basal Melting ◽  
The Antarctic ◽  
Term Response

The primary effects of global warming on the Antarctic ice sheet can involve increases in surface melt for limited areas at lower elevations, increases in net accumulation, and increased basal melting under floating ice. For moderate global wanning, resulting in ocean temperature increases of a few °C, the large- increase in basal melting can become the dominant factor in the long-term response of the ice sheet. The results from ice-sheet modelling show that the increased basal melt rates lead to a reduction of the ice shelves, increased strain rates and flow at the grounding lines, then thinning and floating of the marine ice sheets, with consequential further basal melting. The mass loss from basal melting is counteracted to some extent by the increased accumulation, but in the long term the area of ice cover decreases, particularly in West Antarctica, and the mass loss can dominate. The ice-sheet ice-shelf model of Budd and others (1994) with 20 km resolution has been modified and used to carry out a number of sensitivity studies of the long-term response of the ice sheet to prescribed amounts of global warming. The changes in the ice sheet are computed out to near-equilibrium, but most of the changes take place with in the first lew thousand years. For a global mean temperature increase of 3°C with an ice-shelf basal melt rate of 5 m a−1 the ice shelves disappear with in the first few hundred years, and the marine-based parts of the ice sheet thin and retreat. By 2000 years the West Antarctic region is reduced to a number of small, isolated ice caps based on the bedrock regions which are near or above sea level. This allows the warmer surface ocean water to circulate through the archipelago in summer, causing a large change to the local climate of the region.

The role of ocean circulation in the propagation of rifts on ice shelves.

2020 ◽  
Author(s):  
Mattia Poinelli ◽  
Eric Larour ◽  
Riccardo Riva
Keyword(s):  
Ocean Circulation ◽  
Regional Climate ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Ice Shelves ◽  
Rift Propagation ◽  
Break Up ◽  
The Stability

<p>The break-up of large ice shelves and the associated loss of ice are thought to play a destabilizing role in the ice sheet dynamics. Although ice shelves are a substantial buttressing source in the stability of continental ice sheets, the propagation of large rifts eventually leads to the break-up of icebergs into the ocean. As consequence, this loss of ice would trigger further glacier acceleration and ice sheets retreat, destabilizing the ice cap. Retreat and collapse of ice sheets are also thought to be related to regional climate warming. Indeed, satellite observations suggest that a warming surrounding would induce the ice sheet to progressive thinning and weakening.</p><p>The prolongation of un-grounded ice into the ocean is often interrupted by the propagation of fractures that eventually separates large icebergs from the ice shelf. These fractures are called rifts and range from dimensions of 10 to 100 km. A recent example of such phenomena is the massive break-up of the Larsen C in July, 2017 which followed the disintegration of Larsen A in 1995 and the partial break-up of Larsen B in 2002. The tabular iceberg formed by Larsen C was limited by the propagation of a large rift that began in summer 2016, although the ice shelf had already been thinning since 1992.</p><p>Rift initiation and propagation are thought to be the result of glaciological and oceanographic sources that trigger ice to break. Nonetheless, exact mechanisms remain elusive. The on-going project focuses on ice-ocean interactions in ice shelves that accommodate rifts by using oceanographic models. The goal is to couple rift propagation and ocean circulation underneath ice cavities in order to infer how basal melting affects the development of rifts. The numerical framework is developed within the capabilities of the MITgcm. We aim to identify the sensitivity of propagation rate and opening rate of rifts to variations in the ocean circulation that have occurred during the separation of part of the ice shelf.</p><p>On a larger scale, we are interested in the role of rifting in the stability of Antarctic shelves. Therefore, we work toward a better understanding of which processes are involved in the triggering of rift propagation.</p>

Investigating the climate variability in the Totten area using NEMO-LIM regional model.

2020 ◽  
Author(s):  
Guillian Van Achter ◽  
Charles Pelletier ◽  
Thierry Fichefet
Keyword(s):  
Sea Ice ◽  
Ocean Circulation ◽  
Ice Sheet ◽  
Regional Model ◽  
East Antarctica ◽  
Shelf Area ◽  
Ice Shelf ◽  
Ice Shelves ◽  
Landfast Ice ◽  
The Impact

<p>The Totten ice shelf drains over 570 000 km² of East Antarctica. Most of the ice sheet that drains through the Totten ice-shelf is from Aurora Subglacial Basin and is marine based making the region potentially vulnerable to rapid ice sheet colapse.<br>Understanding how the changes in ocean circulation and properties are causing increased basal melt of Antarctic ice shelves is crucial for predicting future sea level rise.<br>In the context of the The PARAMOUR project (decadal predictability and variability of polar climate: the role of atmosphere-ocean-cryosphere multiscale interaction), we use a high resolution NEMO-LIM 3.6 regional model to investigate the variability and the predictability of the coupled climate system over the Totten area in East Antarctica.<br>In this poster, we will present our on-going work about the impact of landfast ice over the variability of the system. Landfast ice is sea ice that is fastened to the coastline, to the sea floor along shoals or to grouded icebergs. Current sea ice models are unable to represent very crudely the formation, maintenance and decay of coastal landfast ice. We applyed several parameterization for modeling landfast ice over the Totten ice shelf area.</p>

scholarly journals The Response of Ice Shelf Basal Melting to Variations in Ocean Temperature

2008 ◽  
Vol 21 (11) ◽  
pp. 2558-2572 ◽  
Author(s):  
Paul R. Holland ◽  
Adrian Jenkins ◽  
David M. Holland
Keyword(s):  
Ocean Circulation ◽  
General Circulation ◽  
Ice Sheet ◽  
Circulation Model ◽  
Ocean General Circulation Model ◽  
Flow Speed ◽  
Ocean Warming ◽  
Ocean Temperature ◽  
Ice Shelf ◽  
Ice Shelves

Abstract A three-dimensional ocean general circulation model is used to study the response of idealized ice shelves to a series of ocean-warming scenarios. The model predicts that the total ice shelf basal melt increases quadratically as the ocean offshore of the ice front warms. This occurs because the melt rate is proportional to the product of ocean flow speed and temperature in the mixed layer directly beneath the ice shelf, both of which are found to increase linearly with ocean warming. The behavior of this complex primitive equation model can be described surprisingly well with recourse to an idealized reduced system of equations, and it is shown that this system supports a melt rate response to warming that is generally quadratic in nature. This study confirms and unifies several previous examinations of the relation between melt rate and ocean temperature but disagrees with other results, particularly the claim that a single melt rate sensitivity to warming is universally valid. The hypothesized warming does not necessarily require a heat input to the ocean, as warmer waters (or larger volumes of “warm” water) may reach ice shelves purely through a shift in ocean circulation. Since ice shelves link the Antarctic Ice Sheet to the climate of the Southern Ocean, this finding of an above-linear rise in ice shelf mass loss as the ocean steadily warms is of significant importance to understanding ice sheet evolution and sea level rise.

scholarly journals A ground-based radar for measuring vertical strain rates and time-varying basal melt rates in ice sheets and shelves

2015 ◽  
Vol 61 (230) ◽  
pp. 1079-1087 ◽  
Author(s):  
Keith W. Nicholls ◽  
Hugh F.J. Corr ◽  
Craig L. Stewart ◽  
Lai Bun Lok ◽  
Paul V. Brennan ◽  
...  
Keyword(s):  
Ice Sheets ◽  
Ice Sheet ◽  
Present System ◽  
Strain Rates ◽  
Time Varying ◽  
Ice Shelf ◽  
Ross Ice Shelf ◽  
Vertical Strain ◽  
Phase Sensitive

AbstractThe ApRES (autonomous phase-sensitive radio-echo sounder) instrument is a robust, lightweight and relatively inexpensive radar that has been designed to allow long-term, unattended monitoring of ice-shelf and ice-sheet thinning. We describe the instrument and demonstrate its capabilities and limitations by presenting results from three trial campaigns conducted in different Antarctic settings. Two campaigns were ice sheet-based – Pine Island Glacier and Dome C – and one was conducted on the Ross Ice Shelf. The ice-shelf site demonstrates the ability of the instrument to collect a time series of basal melt rates; the two grounded ice applications show the potential to recover profiles of vertical strain rate and also demonstrate some of the limitations of the present system.

Sign in / Sign up

Export Citation Format

Share Document