Inference of the bottom topography in anisothermal mildly-sheared shallow ice flows

• Main Authors: Monnier, Jerome, Zhu, Jiamin
• Other Authors: Institut de Mathématiques de Toulouse UMR5219 (IMT), Université Toulouse Capitole (UT Capitole), Université de Toulouse (UT)-Université de Toulouse (UT)-Institut National des Sciences Appliquées - Toulouse (INSA Toulouse), Institut National des Sciences Appliquées (INSA)-Université de Toulouse (UT)-Institut National des Sciences Appliquées (INSA)-Université Toulouse - Jean Jaurès (UT2J), Université de Toulouse (UT)-Université Toulouse III - Paul Sabatier (UT3), Université de Toulouse (UT)-Centre National de la Recherche Scientifique (CNRS)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2019
• Subjects: [MATH.MATH-OC]Mathematics [math]/Optimization and Control [math.OC]; [MATH.MATH-NA]Mathematics [math]/Numerical Analysis [math.NA]; [SDU.STU.GL]Sciences of the Universe [physics]/Earth Sciences/Glaciology; [SPI.MECA.MEFL]Engineering Sciences [physics]/Mechanics [physics.med-ph]/Fluids mechanics [physics.class-ph]; Antarc*; Antarctica
• Online Access: https://hal.science/hal-01827991; https://hal.science/hal-01827991v3/document; https://hal.science/hal-01827991v3/file/Monnier-Zhu-CMAME-2019.pdf

International audience This study proposes an inverse method to infer the bed topography beneath ice flows from the surface observations (e.g. altimetry elevations and InSAR velocities) and sparse depth measurements (e.g. from airborne campaigns). The flow model is valid for highly to mildly-sheared regimes (hence mildly-rapid) and takes into account varying vertical thermal profiles; it is depth-integrated (long-wave assumption). The inverse problem is particularly challenging since the assimilated surface signatures integrate the bottom features (bed elevation and friction-slip amount) and the internal deformation due to non constant rate factor vertical profile. The first key step of this multi-physics flow inversion is a re-derivation of the anisothermal xSIA model (lubrication type model for generalized Newtonian fluids) leading to a Reduced Uncertainty (RU) version presenting a single uncertain multi-physic parameter γ; that is the so-called RU-SIA equation. The next key steps are advanced Variational Data Assimilation (VDA) formulations combined with a stochastic extension of γ based on the trend observed in the in-situ measurements (e.g. along the flight tracks). The resulting method provides the first physical-based depth (ice thickness) inversions in mildly-sheared mildly-slippery shallow flows. Numerical results are presented in a poorly monitored inland Antarctica area. The uncertainty of the estimated bedrock elevation is noticeably reduced compared to the current estimations uncertainties. The robustness of the inversion process is demonstrated through numerous numerical experiments and empirical sensitivity analyses. The new RU-SIA model may provide a-posteriori estimations of the thermal basal boundary layer too.