skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: glaciome1D: python code for modeling ice mélange flow, 2024
glaciome1D is a quasi-one-dimensional continuum model for modeling the flow of dense packs of icebergs (ice mélange) found in some fjords. In many respects the model is similar to one-dimensional models of ice streams and ice shelves, except that it uses the nonlocal granular fluidity rheology of Henann and Kamrin (2013). The model was created with the intention of developing coupled glacier-ocean-melange models. This is reflected in the modeling framework, which mimics that used for ice streams and ice shelves (Schoof, 2007). Using the model involves creating an instance of the glaciome class, which contains information on the glacier velocity, viscosity (granular fluidity), and geometry as well as model parameters and various external forcings. The glaciome class includes several basic and easy to use functions, such as: self.diagnostic(), self.prognostic(), self.steadystate(), self.save(). The model physics and numerics are described in detail in Amundson et al. (in press). This data set includes a single python module that includes all of the functions for setting up and running the model, an example script that runs the model, and a conda environment list that contains python modules that the code has been tested on.  more » « less
Award ID(s):
2025692
PAR ID:
10625262
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
NSF Arctic Data Center
Date Published:
Subject(s) / Keyword(s):
iceberg ice mélange dynamics granular materials
Format(s):
Medium: X Other: text/xml
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, The Cryosphere)
    Field and remote sensing studies suggest that ice mélange influences glacier–fjord systems by exerting stresses on glacier termini and releasing large amounts of freshwater into fjords. The broader impacts of ice mélange over long timescales are unknown, in part due to a lack of suitable ice mélange flow models. Previous efforts have included modifying existing viscous ice shelf models, despite the fact that ice mélange is fundamentally a granular material, and running computationally expensive discrete element simulations. Here, we draw on laboratory studies of granular materials, which exhibit viscous flow when stresses greatly exceed the yield point, plug flow when the stresses approach the yield point, and exhibit stress transfer via force chains. By implementing the nonlocal granular fluidity rheology into a depth- and width-integrated stress balance equation, we produce a numerical model of ice mélange flow that is consistent with our understanding of well-packed granular materials and that is suitable for long-timescale simulations. For parallel-sided fjords, the model exhibits two possible steady-state solutions. When there is no calving of icebergs or melting of previously calved icebergs, the ice mélange is pushed down-fjord by the advancing glacier terminus, the velocity is constant along the length of the fjord, and the thickness profile is exponential. When calving and melting are included and treated as constants, the ice mélange evolves into another steady state in which its location is fixed relative to the fjord walls, the thickness profile is relatively steep, and the flow is extensional. For the latter case, the model predicts that the steady-state ice mélange buttressing force depends on the surface and basal melt rates through an inverse power-law relationship, decays roughly exponentially with both fjord width and gradient in fjord width, and increases with the iceberg calving flux. The buttressing force appears to increase with calving flux (i.e., glacier thickness) more rapidly than the force required to prevent the capsizing of full-glacier-thickness icebergs, suggesting that glaciers with high calving fluxes may be more strongly influenced by ice mélange than those with small fluxes. 
    more » « less
  2. ; ; ; (, EGU)
    Abstract. Field and remote sensing studies suggest that ice mélange influences glacier-fjord systems by exerting stresses on glacier termini and releasing large amounts of freshwater into fjords. The broader impacts of ice mélange over long time scales are unknown, in part due to a lack of suitable ice mélange flow models. Previous efforts have included modifying existing viscous ice shelf models, despite the fact that ice mélange is fundamentally a granular material, and running computationally expensive discrete element simulations. Here, we draw on laboratory studies of granular materials, which exhibit viscous flow when stresses greatly exceed the yield point, plug flow when the stresses approach the yield point, and stress transfer via force chains. By implementing the nonlocal granular fluidity rheology into a depth- and width-integrated stress balance equation, we produce a numerical model of ice mélange flow that is consistent with our understanding of well-packed granular materials and that is suitable for long time-scale simulations. For parallel-sided fjords, the model exhibits two possible steady state solutions. When there is no calving of new icebergs or melting of previously calved icebergs, the ice mélange is pushed down fjord by the advancing glacier terminus, the velocity is constant along the length of the fjord, and the thickness profile is exponential. When calving and melting are included, the ice mélange evolves to another steady state in which its location is fixed relative to the fjord walls, the thickness profile is relatively steep, and the flow is extensional. For the latter case, the model predicts that the steady-state ice mélange buttressing force depends on the surface and basal melt rates through an inverse power law relationship, decays roughly exponentially with both fjord width and gradient in fjord width, and increases with the iceberg calving flux. The increase in buttressing force with the calving flux, which depends on glacier thickness, appears to occur more rapidly than the force required to prevent the capsize of full-glacier-thickness icebergs, suggesting that glaciers with high calving fluxes may be more strongly influenced by ice mélange than those with small fluxes. 
    more » « less
  3. ; ; ; (, Journal of Glaciology)
    Abstract The Antarctic Ice Sheet loses mass via its ice shelves predominantly through two processes: basal melting and iceberg calving. Iceberg calving is episodic and infrequent, and not well parameterized in ice-sheet models. Here, we investigate the impact of hydrostatic forces on calving. We develop two-dimensional elastic and viscous numerical frameworks to model the ‘footloose’ calving mechanism. This mechanism is triggered by submerged ice protrusions at the ice front, which induce unbalanced buoyancy forces that can lead to fracturing. We compare the results to identify the different roles that viscous and elastic deformations play in setting the rate and magnitude of calving events. Our results show that, although the bending stresses in both frameworks share some characteristics, their differences have important implications for modeling the calving process. In particular, the elastic model predicts that maximum stresses arise farther from the ice front than in the viscous model, leading to larger calving events. We also find that the elastic model would likely lead to more frequent events than the viscous one. Our work provides a theoretical framework for the development of a better understanding of the physical processes that govern glacier and ice-shelf calving cycles. 
    more » « less
  4. ; ; ; ; ; (, Geoscientific Model Development)
    Abstract. Geoscientific models are facing increasing challenges to exploit growing datasets coming from remote sensing. Universal differential equations (UDEs), aided by differentiable programming, provide a new scientific modelling paradigm enabling both complex functional inversions to potentially discover new physical laws and data assimilation from heterogeneous and sparse observations. We demonstrate an application of UDEs as a proof of concept to learn the creep component of ice flow, i.e. a nonlinear diffusivity differential equation, of a glacier evolution model. By combining a mechanistic model based on a two-dimensional shallow-ice approximation partial differential equation with an embedded neural network, i.e. a UDE, we can learn parts of an equation as nonlinear functions that then can be translated into mathematical expressions. We implemented this modelling framework as ODINN.jl, a package in the Julia programming language, providing high performance, source-to-source automatic differentiation (AD) and seamless integration with tools and global datasets from the Open Global Glacier Model in Python. We demonstrate this concept for 17 different glaciers around the world, for which we successfully recover a prescribed artificial law describing ice creep variability by solving ∼ 500 000 ordinary differential equations in parallel. Furthermore, we investigate which are the best tools in the scientific machine learning ecosystem in Julia to differentiate and optimize large nonlinear diffusivity UDEs. This study represents a proof of concept for a new modelling framework aiming at discovering empirical laws for large-scale glacier processes, such as the variability in ice creep and basal sliding for ice flow, and new hybrid surface mass balance models. 
    more » « less
  5. (, William & Mary. Virginia Institute of Marine Science)
    A three-dimensional numerical model of the Amundsen Sea (Antarctica) was used to simulate the period Jan.2006-Mar.2022 under consistent atmospheric/oceanic forcings, bathymetry/ice shelf topography, and model equations/parameters. The model is an implementation of the Regional Ocean Modeling System (ROMS, https://www.myroms.org/) with extensions for sea ice (Budgell 2005) and ice shelves (Dinniman et al. 2011). It simulates the ocean hydrography and circulation, sea ice thermodynamics and dynamics, and the basal melt of the ice shelves, with a uniform horizontal mesh of 1.5km and 20 topography-following vertical levels. Forcings include the ERA5 reanalysis (3-hourly), 10 tidal constituents from CATS 2008, and ocean/sea ice conditions at the edges of the model domain taken from the 5km-resolution circumpolar model of Dinniman et al. 2020 and from daily SSM/I satellite images. The model outputs are divided into nine directories each containing two years worth of model results (run661-669) in the NetCDF format. Each directory contains: daily-averaged model fields (roms_avg_xxxx.nc), instantaneous snapshots every 3 hours for select fields (roms_qck_xxxx.nc), and instantaneous snapshots every 30 days (roms_his_xxxx.nc). All the metadata information necessary for the interpretation of the model outputs (dimensions, units, etc) is included inside the NetCDF files. The NetCDF files follow the CF conventions and can be opened with various software that are open source and freely available over the Internet. In addition to the model outputs, this archive includes the computer code as well as the input files necessary for reproducing the model outputs of this archive. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email