scholarly journals Predicting ocean-induced ice-shelf melt rates using a machine learning image segmentation approach

2022 ◽  
Author(s):  
Sebastian Harry Reid Rosier ◽  
Christopher Bull ◽  
G. Hilmar Gudmundsson

Abstract. Through their role in buttressing upstream ice flow, Antarctic ice shelves play an important part in regulating future sea level change. Reduction in ice-shelf buttressing caused by increased ocean-induced melt along their undersides is now understood to be one of the key drivers of ice loss from the Antarctic Ice Sheet. However, despite the importance of this forcing mechanism most ice-sheet simulations currently rely on simple melt-parametrisations of this ocean-driven process, since a fully coupled ice-ocean modelling framework is prohibitively computationally expensive. Here, we provide an alternative approach that is able to capture the greatly improved physical description of this process provided by large-scale ocean-circulation models over currently employed melt-parameterisations but with trivial computational expense. We introduce a new approach that brings together deep learning and physical modelling to develop a deep neural network framework, MELTNET, that can emulate ocean model predictions of sub-ice shelf melt rates. We train MELTNET on synthetic geometries, using the NEMO ocean model as a ground-truth in lieu of observations to provide melt rates both for training and to evaluate the performance of the trained network. We show that MELTNET can accurately predict melt rates for a wide range of complex synthetic geometries and outperforms more traditional parameterisations for > 95 % of geometries tested. Furthermore, we find MELTNET's melt rate estimates show sensitivity to established physical relationships such as a changes in thermal forcing and ice shelf slope. This study demonstrates the potential for a deep learning framework to calculate melt rates with almost no computational expense, that could in the future be used in conjunction with an ice sheet model to provide predictions for large-scale ice sheet models.

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

<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>


2015 ◽  
Vol 8 (11) ◽  
pp. 9859-9924 ◽  
Author(s):  
X. S. Asay-Davis ◽  
S. L. Cornford ◽  
G. Durand ◽  
B. K. Galton-Fenzi ◽  
R. M. Gladstone ◽  
...  

Abstract. Coupled ice sheet-ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects (MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions in the West Antarctic Ice Sheet (WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase (MISMIP+), ocean experiments under the ice shelf-ocean MIP second phase (ISOMIP+) and coupled ice sheet-ocean experiments under the MISOMIP first phase (MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations (MISOMIP1) to be compared directly to the individual component simulations (MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for evaluation of the participating models.


1998 ◽  
Vol 27 ◽  
pp. 99-104 ◽  
Author(s):  
K. Grosfeld ◽  
R. Gerdes

We investigate the sensitivity of the ocean circulation in the Filchner Trough to changes in the large-scale oceanic environment and its impact on the mass balance of the Filchner Ice Shelf, Antarctica. Three experiments with a three-dimensional ocean model describe (i) the current situation, (ii) a scenario with increased ocean temperatures, and (iii) a scenario with reduced sea-ice formation rates on the adjacent continental shelf. in the final discussion brief results of a combined scenario with increased ocean temperatures and reduced sea-ice formation are presented. The changes from the current situation affect the circulation in the Filchner Trough, and melting and freezing processes beneath the ice shelf. The latter affect the amount and properties of Ice Shelf Water (ISW), a component of Antarctic Bottom Water. Net basal melt rates provide an overall measure for the changes: while the control run yields 0.35 m a−1 net melting averaged over the Filchner Ice Shelf area, the warming scenario results in a more than twofold increase in ice-shelf mass loss. Reduced production of High Salinity ShelfWater due to smaller sea-ice formation rates in the second scenario leads, on the other hand, to a decrease in basal mass loss, because the deep cavity is less well ventilated than in the control run. ISW is cooled and the ice shelf is stabilized under this scenario, which is arguably the more likely development in the southern Weddell Sea.


2021 ◽  
Author(s):  
Darin Comeau ◽  
Xylar Asay-Davis ◽  
Carolyn Begeman ◽  
Matthew Hoffman ◽  
Wuyin Lin ◽  
...  

<p>The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS) -- ice-shelf basal melting and iceberg calving -- are generally poorly represented in current Earth System Models (ESMs). Here, we document the first effort to date at simulating the ocean circulation and exchanges of heat and freshwater within ice-shelf cavities in a coupled ESM, the Department of Energy's Energy Exascale Earth System Model (E3SM). As a step towards full ice-sheet coupling, we implemented static Antarctic ice-shelf cavities and the ability to calculate ice-shelf basal melt rates from the heat and freshwater fluxes computed by the ocean model. In addition, we added the capability to prescribe forcing from iceberg melt, allowing us to realistically represent the other dominant mass loss process from the AIS. In global, low resolution (i.e., non-eddying ocean) simulations, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a tipping point to high melt regimes under certain ice shelves, presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy-induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic and stable ice-shelf basal melt rates across the continent under pre-industrial climate forcing. We also show preliminary results using an ocean/sea-ice grid that makes use of E3SM’s regional-refinement capability, where increased resolution (down to 12km) is placed in the Southern Ocean around Antarctica, bypassing the need for parameterization of eddy-induced transport in this region. The accurate representation of these processes within a coupled ESM is an important step towards reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea-level rise.</p>


2021 ◽  
Author(s):  
Lars Zipf ◽  
Charles Pelletier ◽  
Konstanze Haubner ◽  
Sainan Sun ◽  
Frank Pattyn

<p>Sub-shelf melting is the main driver of Antarctica's ice sheet mass loss. However, sub-shelf melt rate parameterizations for standalone ice models lack the capability to capture complex ocean circulation within ice shelf cavities. To overcome drawbacks of standalone models and to improve melt parameterizations, high resolution coupling of ice sheet and ocean models are capable of hindcasting past decennia and be compared to observations.</p><p>Here, we present first results of a hindcast (1985-2018) of the new circumpolar coupled Southern Ocean – Antarctic ice sheet configuration, developed within the framework of the PARAMOUR project. The configuration is based on the ocean and sea ice model NEMO3.6-LIM3 and the ice sheet model f.ETISh v1.7. The coupling routine facilitates exchange of monthly sub-shelf melt rates (from ocean to ice model) and evolving ice shelf cavity geometry (from ice to ocean model).</p><p>We investigate the impact of the coupling frequency (more precisely, the frequency of updating the ice shelf cavity geometry within the ocean model) on the sub-shelf melt rates and its feedback on the ice dynamics. We further compare the sub-shelf melt rates of the coupled setup to those of the standalone ice sheet model with different sub-shelf melt rate parametrizations (ISMIP6, plume, PICO, PICOP) and investigate the sensitivity of the response of the ice sheet for the different basal melt rate patterns on decadal time scales.</p>


2021 ◽  
Vol 14 (6) ◽  
pp. 3697-3714
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Nicholas Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  

Abstract. The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.


2016 ◽  
Vol 9 (7) ◽  
pp. 2471-2497 ◽  
Author(s):  
Xylar S. Asay-Davis ◽  
Stephen L. Cornford ◽  
Gaël Durand ◽  
Benjamin K. Galton-Fenzi ◽  
Rupert M. Gladstone ◽  
...  

Abstract. Coupled ice sheet–ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects (MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet (WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase (MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase (ISOMIP+) and coupled ice sheet–ocean experiments under the MISOMIP first phase (MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations (MISOMIP1) to be compared directly to the individual component simulations (MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.


2017 ◽  
Author(s):  
Ronja Reese ◽  
Torsten Albrecht ◽  
Matthias Mengel ◽  
Xylar Asay-Davis ◽  
Ricarda Winkelmann

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.


2020 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Nicholas Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  

Abstract. The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degree, respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the entire Antarctic Ice Sheet and the Earth system over time spans on the order of centuries to millennia. In this study we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet model is calculated by PICO from modeled ocean temperatures and salinities at the depth of the continental shelf and, vice versa, the resulting mass and energy fluxes from the melting at the ice-ocean interface are transferred to the ocean model. Mass and energy fluxes are shown to be conserved to machine precision across the considered model domains. The implementation is computationally efficient as it introduces only minimal overhead. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions and time scales between the ice and ocean model in a generic way, and can thus be adopted to a wide range of model setups.


2018 ◽  
Vol 12 (6) ◽  
pp. 1969-1985 ◽  
Author(s):  
Ronja Reese ◽  
Torsten Albrecht ◽  
Matthias Mengel ◽  
Xylar Asay-Davis ◽  
Ricarda Winkelmann

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.


Sign in / Sign up

Export Citation Format

Share Document