skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, January 16 until 2:00 AM ET on Friday, January 17 due to maintenance. We apologize for the inconvenience.

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. ; ; ; ( , 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
  3. ; ; ; ; ; ( , 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
  4. ; ( , 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
  5. ; ; ; ( , The Cryosphere)
    Abstract. Antarctic ice shelves are vulnerable to warming ocean temperatures, and some have already begun thinning in response to increased basal melt rates.Sea level is therefore expected to rise due to Antarctic contributions, but uncertainties in its amount and timing remain largely unquantified. In particular, there is substantial uncertainty in future basal melt rates arising from multi-model differences in thermal forcing and how melt rates depend on that thermal forcing. To facilitate uncertainty quantification in sea level rise projections, we build, validate, and demonstrate projections from a computationally efficient statistical emulator of a high-resolution (4 km) Antarctic ice sheet model, the Community Ice Sheet Model version 2.1. The emulator is trained to a large (500-member) ensemble of 200-year-long 4 km resolution transient ice sheet simulations, whereby regional basal melt rates are perturbed by idealized (yet physically informed) trajectories. The main advantage of our emulation approach is that by sampling a wide range of possible basal melt trajectories, the emulator can be used to (1) produce probabilistic sea level rise projections over much larger Monte Carlo ensembles than are possible by direct numerical simulation alone, thereby providing better statistical characterization of uncertainties, and (2) predict the simulated ice sheet response under differing assumptions about basal melt characteristics as new oceanographic studies are published, without having to run additional numerical ice sheet simulations. As a proof of concept, we propagate uncertainties about future basal melt rate trajectories, derived from regional ocean models, to generate probabilistic sea level rise estimates for 100 and 200 years into the future. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email