Consistent Point Data Assimilation in Firedrake and Icepack

When estimating quantities and fields that are difficult to measure directly, such as the fluidity of ice, from point data sources, such as satellite altimetry, it is important to solve a numerical inverse problem that is formulated with Bayesian consistency. Otherwise, the resultant probability density function for the difficult to measure quantity or field will not be appropriately clustered around the truth. In particular, the inverse problem should be formulated by evaluating the numerical solution at the true point locations for direct comparison with the point data source. If the data are first fitted to a gridded or meshed field on the computational grid or mesh, and the inverse problem formulated by comparing the numerical solution to the fitted field, the benefits of additional point data values below the grid density will be lost. We demonstrate, with examples in the fields of groundwater hydrology and glaciology, that a consistent formulation can increase the accuracy of results and aid discourse between modellers and observationalists. To do this, we bring point data into the finite element method ecosystem as discontinuous fields on meshes of disconnected vertices. Point evaluation can then be formulated as a finite element interpolation operation (dual-evaluation). This new abstraction is well-suited to automation, including automatic differentiation. We demonstrate this through implementation in Firedrake, which generates highly optimised code for solving Partial Differential Equations (PDEs) with the finite element method. Our solution integrates with dolfin-adjoint/pyadjoint, allowing PDE-constrained optimisation problems, such as data assimilation, to be solved through forward and adjoint mode automatic differentiation.

[1]  Barry F. Smith,et al.  PETSc/TAO Users Manual , 2022 .

[2]  David A. Ham,et al.  Code Generation for Productive, Portable, and Scalable Finite Element Simulation in Firedrake , 2021, Computing in Science & Engineering.

[3]  I. Joughin,et al.  icepack: a new glacier flow modeling package in Python, version 1.0 , 2021, Geoscientific Model Development.

[4]  C. Richardson,et al.  LEoPart: a particle library for FEniCS , 2019, Comput. Math. Appl..

[5]  Won Sang Lee,et al.  Deep glacial troughs and stabilizing ridges unveiled beneath the margins of the Antarctic ice sheet , 2019, Nature Geoscience.

[6]  B. Scheuchl,et al.  Continent‐Wide, Interferometric SAR Phase, Mapping of Antarctic Ice Velocity , 2019, Geophysical Research Letters.

[7]  Joel Nothman,et al.  SciPy 1.0-Fundamental Algorithms for Scientific Computing in Python , 2019, ArXiv.

[8]  Simon W. Funke,et al.  dolfin-adjoint 2018.1: automated adjoints for FEniCS and Firedrake , 2019, J. Open Source Softw..

[9]  Lawrence Mitchell,et al.  Thetis coastal ocean model: discontinuous Galerkin discretization for the three-dimensional hydrostatic equations , 2017, Geoscientific Model Development.

[10]  David A. Ham,et al.  Exposing and exploiting structure: optimal code generation for high-order finite element methods , 2017, ArXiv.

[11]  Lawrence Mitchell,et al.  TSFC: a structure-preserving form compiler , 2017, SIAM J. Sci. Comput..

[12]  I. Joughin,et al.  Basal resistance for three of the largest Greenland outlet glaciers , 2016 .

[13]  Anders Logg,et al.  The FEniCS Project Version 1.5 , 2015 .

[14]  Matthew G. Knepley,et al.  Efficient Mesh Management in Firedrake Using PETSc DMPlex , 2015, SIAM J. Sci. Comput..

[15]  Andrew T. T. McRae,et al.  Firedrake: automating the finite element method by composing abstractions , 2015, ACM Trans. Math. Softw..

[16]  Anders Logg,et al.  Unified form language: A domain-specific language for weak formulations of partial differential equations , 2012, TOMS.

[17]  Patrick E. Farrell,et al.  Directional integration on unstructured meshes via supermesh construction , 2012, J. Comput. Phys..

[18]  David A. Ham,et al.  Automated Derivation of the Adjoint of High-Level Transient Finite Element Programs , 2012, SIAM J. Sci. Comput..

[19]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[20]  Anders Logg,et al.  Efficient representation of computational meshes , 2009, Int. J. Comput. Sci. Eng..

[21]  H. Blatter,et al.  Dynamics of Ice Sheets and Glaciers , 2009 .

[22]  Matthew D. Piggott,et al.  Conservative interpolation between unstructured meshes via supermesh construction , 2009 .

[23]  T. Painter,et al.  MODIS-based Mosaic of Antarctica (MOA) data sets: Continent-wide surface morphology and snow grain size , 2007 .

[24]  A. Vieli,et al.  Numerical modelling and data assimilation of the Larsen B ice shelf, Antarctic Peninsula , 2006, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[25]  Robert C. Kirby,et al.  Algorithm 839: FIAT, a new paradigm for computing finite element basis functions , 2004, TOMS.

[26]  Ian Joughin,et al.  Basal shear stress of the Ross ice streams from control method inversions , 2004 .

[27]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[28]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[29]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[30]  D. Macayeal The basal stress distribution of Ice Stream E, Antarctica, inferred by control methods , 1992 .

[31]  R. Dennis Cook,et al.  Cross-Validation of Regression Models , 1984 .

[32]  Sebastian Kenji Mitusch,et al.  An Algorithmic Differentiation Tool for FEniCS , 2018 .

[33]  David A. Ham,et al.  Mesh dependence in PDE-constrained optimisation , 2017 .

[34]  S. Barnett,et al.  Philosophical Transactions of the Royal Society A : Mathematical , 2017 .

[35]  Martin Sandve Alnæs,et al.  UFL: a finite element form language , 2012 .

[36]  M. Truffer,et al.  Reconstruction of basal properties in ice sheets using iterative inverse methods , 2012, Journal of Glaciology.

[37]  D. Macayeal,et al.  Basal friction of Ice Stream E, West Antarctica , 1995, Journal of Glaciology.

[38]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[39]  Grace L. Yang,et al.  On Bayes Procedures , 1990 .