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

DOI10.1016/J.APM.2022.03.001MaRDI QIDQ2109570zbMATHOpenAlexWikidataFDO

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)

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)

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.

Cites work

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

Describes a project that uses

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)