EGU21-3506
https://doi.org/10.5194/egusphere-egu21-3506
EGU General Assembly 2021
© Author(s) 2023. This work is distributed under
the Creative Commons Attribution 4.0 License.

Assimilating sparse data in glaciological inverse problems

Daniel Shapero1 and Reuben Nixon-Hill2
Daniel Shapero and Reuben Nixon-Hill
  • 1Applied Physics Lab, University of Washington, Seattle, United States of America (shapero@uw.edu)
  • 2Department of Mathematics, Imperial College, London, United Kingdom (reuben.nixon-hill10@imperial.ac.uk)

Most of the existing work on solving inverse problems in glaciology has assumed that the observational data used to constrain the model are spatially dense. This assumption is very convenient because it means that the model-data misfit term in the objective functional can be written as an integral. In many scenarios, however, the computational mesh can locally be much finer than the observational grid, or the observations can have large patches of missing data. Moreover, pretending as if the observations are a globally-defined continuous field obscures valuable information about the number of independent measurements we have. It is then impossible to apply a posteriori sanity checks on the expected model-data misfit from regression theory. Here we'll describe some recent work we've done on assimilating sparse point data into ice flow models and how this allows us to be more rigorous about the statistical interpretation of our results. For now we are focusing on the kinds of inverse problems that have been solved in the glaciology literature for a long time -- inferring rheology and basal friction from surface velocities. But these developments open up the possibility of assimilating new sources of data, such as measurements from strain gauges or ice cores.

How to cite: Shapero, D. and Nixon-Hill, R.: Assimilating sparse data in glaciological inverse problems, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-3506, https://doi.org/10.5194/egusphere-egu21-3506, 2021.