An augmented Lagrangian algorithm for recovery of ice thickness in unidirectional flow using the shallow ice approximation

DOI10.1016/J.APM.2022.03.001zbMATH Open1503.86003arXiv2108.00854OpenAlexW3188078541WikidataQ114208677 ScholiaQ114208677MaRDI QIDQ2109570FDOQ2109570

Authors: Elizabeth K. McGeorge, Phillip L. Wilson, M. Sellier, M. Moyers-Gonzalez 

Publication date: 21 December 2022

Published in: Applied Mathematical Modelling (Search for Journal in Brave)

Abstract: A key parameter in ice flow modelling is basal slipping at the ice-bed interface as it can have a large effect on the resultant ice thickness. Unfortunately, its contribution to surface observations can be hard to distinguish from that of bed undulations. Therefore, inferring the ice thickness from surface measurements is an interesting and non-trivial inverse problem. This paper presents a method for recovering dually the ice thickness and the basal slip using only surface elevation and speed measurements. The unidirectional shallow ice approximation is first implemented to model steady state ice flow for given bedrock and basal slip profiles. This surface is then taken as synthetic observed data. An augmented Lagrangian algorithm is then used to find the diffusion coefficient which gives the best fit to observations. Combining this recovered diffusion with observed surface velocity, a simple Newton's method is used to recover both the ice thickness and basal slip. The method was successful in each test case and this implies that it should be possible to recover both of these parameters in two-dimensional cases also.

Full work available at URL: https://arxiv.org/abs/2108.00854

zbMATH Keywords

inverse problemsshallow ice approximationbasal slipice flowsice thickness

Mathematics Subject Classification ID

Inverse problems in geophysics (86A22) Inverse problems in optimal control (49N45) Hydrology, hydrography, oceanography (86A05) Glaciology (86A40) Inverse problems in fluid mechanics (76M21)

Cites Work

• Solving PDEs in Python. The FEniCS tutorial I
• Title not available (Why is that?)
• Automated solution of differential equations by the finite element method. The FEniCS book
• An efficient linear solver for nonlinear parameter identification problems
• An Augmented Lagrangian Method for Identifying Discontinuous Parameters in Elliptic Systems
• Inverse acoustic and electromagnetic scattering theory
• Inverse rheometry and basal properties inference for pseudoplastic geophysical flows
• Inference of the bottom topography in anisothermal mildly-sheared shallow ice flows
• Bedrock reconstruction from free surface data for unidirectional glacier flow with basal slip
• Inference of the bottom properties in shallow ice approximation models
• Accelerating the Uzawa algorithm
• On the cost of solving augmented Lagrangian subproblems
• Title not available (Why is that?)

Uses Software

• FEniCS
• UFC

This page was built for publication: An augmented Lagrangian algorithm for recovery of ice thickness in unidirectional flow using the shallow ice approximation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2109570)