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: The transferability of adjoint inversion products between different ice flow models
Abstract. Among the most important challenges faced by ice flow models is how to represent basal and rheological conditions, which are challenging to obtain from direct observations. A common practice is to use numerical inversions to calculate estimates for the unknown properties, but there are many possible methods and not one standardised approach. As such, every ice flow model has a unique initialisation procedure. Here we compare the outputs of inversions from three different ice flow models, each employing a variant of adjoint-based optimisation to calculate basal sliding coefficients and flow rate factors using the same observed surface velocities and ice thickness distribution. The region we focus on is the Amundsen Sea Embayment in West Antarctica, the subject of much investigation due to rapid changes in the area over recent decades. We find that our inversions produce similar distributions of basal sliding across all models, despite using different techniques, implying that the methods used are highly robust and represent the physical equations without much influence by individual model behaviours. Transferring the products of inversions between models results in time-dependent simulations displaying variability on the order of or lower than existing model intercomparisons. Focusing on contributions to sea level, the highest variability we find in simulations run in the same model with different inversion products is 32 %, over a 40-year period, a difference of 3.67 mm. There is potential for this to be improved with further standardisation of modelling processes, and the lowest variability within a single model is 13 %, or 1.82 mm over 40 years. While the successful transfer of inversion outputs from one model to another requires some extra effort and technical knowledge of the particular models involved, it is certainly possible and could indeed be useful for future intercomparison projects.  more » « less
Award ID(s):
2152622
PAR ID:
10319977
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
The Cryosphere
Volume:
15
Issue:
4
ISSN:
1994-0424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Glaciology)
    Abstract The time-evolution of glacier basal motion remains poorly constrained, despite its importance in understanding the response of glaciers to climate warming. Athabasca Glacier provides an ideal site for observing changes in basal motion over long timescales. Studies from the 1960s provide an in situ baseline dataset constraining ice deformation and basal motion. We use two complementary numerical flow models to investigate changes along a well-studied transverse profile and throughout a larger study area. A cross-sectional flow model allows us to calculate transverse englacial velocity fields to simulate modern and historical conditions. We subsequently use a 3-D numerical ice flow model, Icepack, to estimate changes in basal friction by inverting known surface velocities. Our results reproduce observed velocities well using standard values for flow parameters. They show that basal motion declined significantly (30–40%) and this constitutes the majority (50–80%) of the observed decrease in surface velocities. At the same time, basal resistive stress has remained nearly constant and now balances a much larger fraction of the driving stress. The decline in basal motion over multiple decades of climate warming could serve as a stabilizing feedback mechanism, slowing ice transport to lower elevations, and therefore moderating future mass loss rates. 
    more » « less
  2. ; ; ; ; ; (, 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
  3. ; ; ; (, Journal of Glaciology)
    Abstract Recent observations of increased discharge through fast-flowing outlet glaciers and ice streams motivate questions concerning the inland migration of regions of fast flow, which could increase drawdown of the ice-sheet interior. To investigate one process that could lead to inland migration we conduct experiments with a two-dimensional, full-stress, transient ice-flow model. An initial steady state is perturbed by initiating a jump in sliding speed over a fraction of the model domain. As a result, longitudinal-stress gradients increase frictional melting upstream from the slow-to-fast sliding transition, and a positive feedback between longitudinal-stress gradients, basal meltwater production and basal sliding causes the sliding transition to migrate upstream over time. The distance and speed of migration depend on the magnitude of the perturbation and on the degree of non-linearity assumed in the link between basal stress and basal sliding: larger perturbations and/or higher degrees of non-linearity lead to farther and faster upstream migration. Migration of the sliding transition causes the ice sheet to thin over time and this change in geometry limits the effects of the positive feedback, ultimately serving to impede continued upstream migration. 
    more » « less
  4. ; ; (, Journal of Geophysical Research: Machine Learning and Computation)
    Abstract Predicting the future contribution of the ice sheets to sea level rise over the next decades presents several challenges due to a poor understanding of critical boundary conditions, such as basal sliding. Traditional numerical models often rely on data assimilation methods to infer spatially variable friction coefficients by solving an inverse problem, given an empirical friction law. However, these approaches are not versatile, as they sometimes demand extensive code development efforts when integrating new physics into the model. Furthermore, this approach makes it difficult to handle sparse data effectively. To tackle these challenges, we use the Physics‐Informed Neural Networks (PINNs) to seamlessly integrate observational data and governing equations of ice flow into a unified loss function, facilitating the solution of both forward and inverse problems within the same framework. We illustrate the versatility of this approach by applying the framework to two‐dimensional problems on the Helheim Glacier in southeast Greenland. By systematically concealing one variable (e.g., ice speed, ice thickness, etc.), we demonstrate the ability of PINNs to accurately reconstruct hidden information. Furthermore, we extend this application to address a challenging mixed inversion problem. We show how PINNs are capable of inferring the basal friction coefficient while simultaneously filling gaps in the sparsely observed ice thickness. This unified framework offers a promising avenue to enhance the predictive capabilities of ice sheet models, reducing uncertainties, and advancing our understanding of poorly constrained physical processes. 
    more » « less
  5. ; (, Journal of Geophysical Research: Earth Surface)
    Abstract Ice shelves regulate the ice‐ocean boundary by buttressing the flux of grounded ice into the ocean and are vulnerable to basal melt, which can lead to ice‐shelf thinning and loss of buttressing. Localized, enhanced basal melt can form basal channels, which may impact ice‐shelf stability. Here we investigate the evolution of the Getz Ice Shelf Basal Channel (GISBC) in West Antarctica using a novel suite of geophysical data, including Reference Elevation Model of Antarctica (REMA) digital elevation models, ICESat‐1 and ‐2 altimetry, Operation IceBridge altimetry and radar, and InSAR‐derived ice flow velocities. We describe basal‐channel and ice‐shelf change in both Eulerian and Lagrangian frameworks and document changes in the channel's shape and its lateral motion and estimate basal melting. We find a high degree of spatial and temporal variability in GISBC evolution, with several locations of active basal incision. Incision occurs at rates of up to 22 m a−1at the head of the channel, which is extending toward the grounding line at a rate of ~1 km a−1. Freeboard heights over areas of rapid basal incision are out of hydrostatic equilibrium. The GISBC is also migrating to the northwest, perpendicular to the northeasterly ice flow direction, at an average rate of 70–80 m a−1. The spatiotemporal variability of evolution of the GISBC motivates further characterization of basal channels and their impact on ice‐shelf stability, so that these effects may more readily be incorporated in ice‐ocean models predicting ice flow and sea‐level rise. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email