Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
doi:10.5194/gmdd-6-5251-2013
© Author(s) 2013. This work is distributed
under the Creative Commons Attribution 3.0 License.
Model description paper
09 Oct 2013
Review status
A revision of this discussion paper 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. McGovern1,*, I. Rutt1, J. Utke2, and T. Murray1 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.


Citation: McGovern, J., Rutt, I., Utke, J., and Murray, T.: ADISM v.1.0: an adjoint of a thermomechanical ice-sheet model obtained using an algorithmic differentiation tool, Geosci. Model Dev. Discuss., 6, 5251-5288, doi:10.5194/gmdd-6-5251-2013, 2013.
J. McGovern et al.
J. McGovern et al.

Viewed

Total article views: 712 (including HTML, PDF, and XML)

HTML PDF XML Total BibTeX EndNote
368 325 19 712 33 29

Views and downloads (calculated since 09 Oct 2013)

Cumulative views and downloads (calculated since 09 Oct 2013)

Saved

Discussed

Latest update: 23 Mar 2017
Publications Copernicus
Download
Share