Papers › Surface elevation errors in finite element Stokes models for glacier evolution

Surface elevation errors in finite element Stokes models for glacier evolution

12 Aug 2024arXiv:2408.06470links table onlyarchive 2025-07-28

Ed Bueler

The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.

The primary data which determine the evolution of glaciation are the bedrock elevation and the surface mass balance. From this data, which we assume is defined over a fixed land region, the glacier's geometry solves a free boundary problem which balances the time derivative of the surface elevation, the surface velocity from the Stokes flow of the ice, and the surface balance rate. This problem can be posed in weak form as a variational inequality over a cone of admissible surface elevation functions, those which are above the bedrock topography. After some preparatory theory for the Stokes problem, we conjecture that the corresponding continuous-space, implicit time-step variational inequality problem is well-posed if the surface kinematical equation is appropriately regularized. This conjecture is supported by physical arguments and numerical evidence. We then prove a general theorem which bounds the numerical error made by finite element approximations of nonlinear-operator variational inequalities in Banach spaces. This bound is a sum of error terms of different types, special to variational inequalities. When it is applied to the implicit time-step glacier problem there are three terms in the bound: an error from discretizing the bed elevation, an error from numerically solving for the Stokes velocity, and finally an expected error which is quasi-optimal in the finite element space representation of the surface elevation. The design of glacier models is then reconsidered based on this a priori error analysis.

PaperPDFCode

Code

bueler/glacier-fe-estimate officialmentioned on GitHub report
bueler/stokes-implicit mentioned on GitHub report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

Not run by Syntology. Nothing on this page verifies that the listed code works.

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections