www.the-cryosphere-discuss.net/3/181/2009/ doi:10.5194/tcd-3-181-2009 © Author(s) 2009. This work is distributed under the Creative Commons Attribution 3.0 License. Estimating basal properties of glaciers from surface measurements: a non-linear Bayesian inversion approach 1Versuchsanstalt für Wasserbau, Hydrologie und Glaziologie (VAW), ETH Zürich, 8092 Zürich, Switzerland 2British Antarctic Survey, High Cross, Madingley Rd., Cambridge CB3 0ET, UK Abstract. We propose a new approach to indirectly estimate basal properties of a glacier, i.e. bedrock topography and basal slipperiness, from observations of surface topography and surface velocities. We demonstrate how a maximum a posteriori estimate of basal conditions can be determined using a Bayesian inference approach in a combination with an analytical linearisation of the forward model. Using synthetic data we show that for non-linear media and non-linear sliding law only a few forward-step model evaluations are needed for convergence. The forward step is solved with a numerical finite-element model using the full Stokes equations. Forward Fréchet derivatives are approximated through analytical small-perturbation solutions. This approximation is a key feature of the method and the effects of this approximation on model performance are analyzed. The number of iterations needed for convergence increases with the amplitude of the basal perturbations, but generally less than ten iterations are needed. Discussion Paper (PDF, 790 KB) Interactive Discussion (Closed, 3 Comments) Final Revised Paper (TC) Citation: Raymond, M. J. and Gudmundsson, G. H.: Estimating basal properties of glaciers from surface measurements: a non-linear Bayesian inversion approach, The Cryosphere Discuss., 3, 181-222, doi:10.5194/tcd-3-181-2009, 2009. Bibtex EndNote Reference Manager XML |
Search TCDRecent Papers |