Preprints
https://doi.org/10.5194/gmdd-6-5251-2013
https://doi.org/10.5194/gmdd-6-5251-2013

Submitted as: model description paper 09 Oct 2013

Submitted as: model description paper | 09 Oct 2013

Review status: this preprint was under review for the journal GMD. A revision for further review has not been submitted.

ADISM v.1.0: an adjoint of a thermomechanical ice-sheet model obtained using an algorithmic differentiation tool

J. McGovern*,1, I. Rutt1, J. Utke2, and T. Murray1 J. McGovern et al.
  • 1Department of Geography, College of Science, Swansea University, Singleton Park, Swansea SA2 8PP, UK
  • 2Argonne National Laboratory, Argonne, IL, USA
  • *now at: Laboratoire de Glaciologie et Géophysique de l'Environnement, UJF-Grenoble, CNRS – UMR5183, Saint-Martin-d'Hères, France

Abstract. A number of problems in contemporary glaciology could benefit from the application of adjoint models. On a simple level, adjoint models can be used to calculate ice-sheet sensitivities with respect to spatially varying parameters such as the basal sliding coefficient. At a more sophisticated level, adjoint models may be used as components of variational data assimilation schemes, allowing problems of model initialization and data-constrained evolution to be tackled.

Fundamentally, adjoint models calculate the sensitivity of a cost function to a suite of control parameters. Such model sensitivities can alternatively be obtained by running the model many times, perturbing each control parameter separately in turn, and calculating the resulting sensitivity in each case. For large numbers of control parameters, however, such as the case where a control parameter corresponds to each point in the model domain, the computational cost becomes prohibitive. The use of adjoint models allows sensitivities to be obtained more efficiently – adjoint model sensitivities are obtained in a single run – and more accurately, since the differentiation of the model is done with machine precision.

We present a finite-difference shallow ice approximation (SIA), thermomechanical ice-sheet model (the forward model), and its adjoint, as generated by using the OpenAD algorithmic differentiation tool. We verify the ice-sheet model using standard SIA benchmark tests and check the consistency between derivatives computed by OpenAD and certain numerically approximated derivatives. Typical adjoint calculations are demonstrated by application to the Greenland ice sheet.

J. McGovern et al.

 
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

J. McGovern et al.

Viewed

Total article views: 1,235 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
710 470 55 1,235 75 74
  • HTML: 710
  • PDF: 470
  • XML: 55
  • Total: 1,235
  • BibTeX: 75
  • EndNote: 74
Views and downloads (calculated since 09 Oct 2013)
Cumulative views and downloads (calculated since 09 Oct 2013)

Cited

Saved

Latest update: 31 Jul 2021
Download