Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 4.713
IF4.713
IF 5-year value: 4.927
IF 5-year
4.927
CiteScore value: 8.0
CiteScore
8.0
SNIP value: 1.425
SNIP1.425
IPP value: 4.65
IPP4.65
SJR value: 2.353
SJR2.353
Scimago H <br class='widget-line-break'>index value: 71
Scimago H
index
71
h5-index value: 53
h5-index53
Volume 8, issue 2
The Cryosphere, 8, 721–741, 2014
https://doi.org/10.5194/tc-8-721-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
The Cryosphere, 8, 721–741, 2014
https://doi.org/10.5194/tc-8-721-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.

Research article 24 Apr 2014

Research article | 24 Apr 2014

Adjoint accuracy for the full Stokes ice flow model: limits to the transmission of basal friction variability to the surface

N. Martin and J. Monnier N. Martin and J. Monnier
  • INSA and Institut de Mathématiques de Toulouse, 135 Avenue de Rangueil, 31400, Toulouse, France

Abstract. This work focuses on the numerical assessment of the accuracy of an adjoint-based gradient in the perspective of variational data assimilation and parameter identification in glaciology. Using noisy synthetic data, we quantify the ability to identify the friction coefficient for such methods with a non-linear friction law. The exact adjoint problem is solved, based on second-order numerical schemes, and a comparison with the so-called "self-adjoint" approximation, neglecting the viscosity dependence on the velocity (leading to an incorrect gradient), common in glaciology, is carried out. For data with a noise of 1%, a lower bound of identifiable wavelengths of 10 ice thicknesses in the friction coefficient is established, when using the exact adjoint method, while the "self-adjoint" method is limited, even for lower noise, to a minimum of 20 ice thickness wavelengths. The second-order exact gradient method therefore provides robustness and reliability for the parameter identification process. In another respect, the derivation of the adjoint model using algorithmic differentiation leads to the formulation of a generalization of the "self-adjoint" approximation towards an incomplete adjoint method, adjustable in precision and computational burden.

Publications Copernicus
Download
Citation