Preprints
https://doi.org/10.5194/tcd-8-5511-2014
https://doi.org/10.5194/tcd-8-5511-2014
28 Oct 2014
 | 28 Oct 2014
Status: this preprint was under review for the journal TC. A revision for further review has not been submitted.

Model calibration for ice sheets and glaciers dynamics: a general theory of inverse problems in glaciology

M. Giudici, F. Baratelli, A. Comunian, C. Vassena, and L. Cattaneo

Abstract. Numerical modelling of the dynamic evolution of ice sheets and glaciers requires the solution of discrete equations which are based on physical principles (e.g. conservation of mass, linear momentum and energy) and phenomenological constitutive laws (e.g. Glen's and Fourier's laws). These equations must be accompanied by information on the forcing term and by initial and boundary conditions (IBCs) on ice velocity, stress and temperature; on the other hand the constitutive laws involve many physical parameters, some of which depend on the ice thermodynamical state. The proper forecast of the dynamics of ice sheets and glaciers requires a precise knowledge of several quantities which appear in the IBCs, in the forcing terms and in the phenomenological laws. As these quantities cannot be easily measured at the study scale in the field, they are often obtained through model calibration by solving an inverse problem (IP). The objective of this paper is to provide a thorough and rigorous conceptual framework for IPs in cryospheric studies and in particular: to clarify the role of experimental and monitoring data to determine the calibration targets and the values of the parameters that can be considered to be fixed; to define and characterise identifiability, a property related to the solution to the forward problem; to study well-posedness in a correct way, without confusing instability with ill-conditioning or with the properties of the method applied to compute a solution; to cast sensitivity analysis in a general framework and to differentiate between the computation of local sensitivity indicators with a one-at-a-time approach and first-order sensitivity indicators that consider the whole possible variability of the model parameters. The conceptual framework and the relevant properties are illustrated by means of a simple numerical example of isothermal ice flow, based on the shallow-ice approximation.

M. Giudici, F. Baratelli, A. Comunian, C. Vassena, and L. Cattaneo
 
Status: closed (peer review stopped)
Status: closed (peer review stopped)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
Printer-friendly Version - Printer-friendly version Supplement - Supplement
 
Status: closed (peer review stopped)
Status: closed (peer review stopped)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
Printer-friendly Version - Printer-friendly version Supplement - Supplement
M. Giudici, F. Baratelli, A. Comunian, C. Vassena, and L. Cattaneo
M. Giudici, F. Baratelli, A. Comunian, C. Vassena, and L. Cattaneo

Viewed

Total article views: 1,864 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
1,179 560 125 1,864 97 124
  • HTML: 1,179
  • PDF: 560
  • XML: 125
  • Total: 1,864
  • BibTeX: 97
  • EndNote: 124
Views and downloads (calculated since 28 Oct 2014)
Cumulative views and downloads (calculated since 28 Oct 2014)

Cited

Saved

Latest update: 28 Mar 2024
Download
Short summary
Inverse problems (IPs) are fundamental to calibrate numerical models, but are not yet widely diffused in glaciology, even if some theroretical developments and excellent applications were published. This paper provides a conceptual framework for IP in cryospheric sciences, which permits to clearly define and discuss the properties of IP. The examined topics include: the role of data and parameters: the link of ill-posedness with ill-conditioning and identifiability; sensitivity analysis.