TY - JOUR

T1 - Estimating glacier basal properties from surface measurements

T2 - A non-linear inversion approach applied to the Rutford Ice Stream, West Antarctica

AU - Raymond, Mélanie

AU - Gudmundsson, G. H.

N1 - This chapter of Melanie Raymond's thesis was co-authored by Hilmar Gudmundsson. The entire thesis was published in Mitteilungen. The chapter later appeared in revised form in Cryosphere (not JGR as suggested in the original) PB 15-2-19.

PY - 2007/12/1

Y1 - 2007/12/1

N2 - We introduce a non-linear Bayesian inversion approach (e.g., Tarantola, 2005; Rodgers, 2000) to estimate the basal properties, that is bedrock topography and basal slipperiness, along a flow line of Rutford ice stream from observations of surface topography and surface velocities. The forward function describing the relationship between the basal properties and the observations is solved numerically with a plane-strain finite-element model. As there exist no explicit solution for the inversion of this non-linear forward function, the solution of the inverse problem is sought numerically and iteratively by using a nonlinear Gauss-Newton procedure. The first order forward model derivatives needed for inversion are approximated by analytical linear transfer functions (Gudmundsson, 2003). This approximation is attractive as it greatly enhances the numerical efficiency of the method by sparing the time-consuming evaluation of the numerical derivatives (Chapter 5 in this thesis; Raymond and Gudmundsson (2007)). Inversions performed with synthetic data showed that the inversion procedure behaves correctly and converges to the correct solution. The basal properties we compute for Rutford ice stream are consistent both with the surface observations and the radar measurements of bedrock topography.

AB - We introduce a non-linear Bayesian inversion approach (e.g., Tarantola, 2005; Rodgers, 2000) to estimate the basal properties, that is bedrock topography and basal slipperiness, along a flow line of Rutford ice stream from observations of surface topography and surface velocities. The forward function describing the relationship between the basal properties and the observations is solved numerically with a plane-strain finite-element model. As there exist no explicit solution for the inversion of this non-linear forward function, the solution of the inverse problem is sought numerically and iteratively by using a nonlinear Gauss-Newton procedure. The first order forward model derivatives needed for inversion are approximated by analytical linear transfer functions (Gudmundsson, 2003). This approximation is attractive as it greatly enhances the numerical efficiency of the method by sparing the time-consuming evaluation of the numerical derivatives (Chapter 5 in this thesis; Raymond and Gudmundsson (2007)). Inversions performed with synthetic data showed that the inversion procedure behaves correctly and converges to the correct solution. The basal properties we compute for Rutford ice stream are consistent both with the surface observations and the radar measurements of bedrock topography.

M3 - Article

AN - SCOPUS:77957267106

SP - 107

EP - 139

JO - Mitteilungen der Versuchsanstalt fur Wasserbau, Hydrologie und Glaziologie an der Eidgenossischen Technischen Hochschule Zurich

JF - Mitteilungen der Versuchsanstalt fur Wasserbau, Hydrologie und Glaziologie an der Eidgenossischen Technischen Hochschule Zurich

SN - 0374-0056

IS - 202

ER -