Source-to-source adjoint algorithmic differentiation of an ice sheet model written in C. (English) Zbl 1453.86048
Summary: Algorithmic Differentiation (AD) has become a powerful tool to improve our understanding of the Earth System, because it can generate adjoint code which permits efficient calculation of gradients that are essential to sensitivity studies, inverse problems, parameter estimation and data assimilation. Most source-to-source transformation tools, however, have been designed for FORTRAN and support for C remains limited. Here we use the Adjoinable Land Ice Flow model (ALIF), a C clone of the C++ Ice Sheet System Model (ISSM) and employ source-to-source AD to produce its adjoint code. We present the first running source-to-source adjoint of ALIF, and its application to basal drag inversion under Pine Island Glacier, West Antarctica. ALIF brought several challenges to AD tool development, such as the correct treatment of the context code, which does not compute the differentiable function, but controls this computation through the setup of data structures, including possible aliasing, as well as data-flow reversal in the presence of pointers and dynamic memory, which are ubiquitous in codes such as ISSM and ALIF. We present the strategies we have developed to overcome these challenges.
MSC:
| 86A40 | Glaciology |
| 68W30 | Symbolic computation and algebraic computation |
| 86-08 | Computational methods for problems pertaining to geophysics |
Keywords:
ice sheet model; ISSM; algorithmic differentiation; adjoint methods; dynamic memory; tapenadeReferences:
| [1] | data representation alternatives in semantically augmented numerical models6th IEEE International Workshop on Source Code Analysis and Manipulation, SCAM 2006, Philadelphia, PA, USA, 2006 |
| [2] | Generating recomputations in reverse mode ADAutomatic Differentiation of Algorithms: From Simulation to OptimizationChapter 33SpringerHeidelberg2002283291Available at |
| [3] | Giles, M. B.; Bischof, C. H.; Bücker, H. M.; Hovland, P. D.; Naumann, U.; Utke, J., Collected matrix derivative results for forward and reverse mode algorithmic differentiation, Advances in Automatic Differentiation, 35-44, (2008), Springer · Zbl 1154.65308 |
| [4] | Goldberg, D. N.; Heimbach, P., parameter and state estimation with a time-dependent adjoint marine ice sheet model, Cryosphere, 7, 1659-1678, (2013) · doi:10.5194/tc-7-1659-2013 |
| [5] | Goldberg, D. N.; Heimbach, P.; Joughin, I.; Smith, B., committed retreat of Smith, pope, and kohler glaciers over the next 30 years inferred by transient model calibration, Cryosphere, 9, 2429-2446, (2015) · doi:10.5194/tc-9-2429-2015 |
| [6] | Griewank, A.; Walther, A., Evaluating Derivatives: Principles and Techniques of Algorithmic Differentiation, (2008), SIAM, Philadelphia, PA · Zbl 1159.65026 |
| [7] | Hascoët, L.; Pascual, V., the tapenade automatic differentiation tool: principles, model, and specification, ACM Trans. Math. Softw., 39, (2013) · Zbl 1295.65026 · doi:10.1145/2450153.2450158 |
| [8] | Hascoët, L.; Utke, J., programming language features, usage patterns, and the efficiency of generated adjoint code, Optim. Methods Softw., 31, 885-903, (2016) · Zbl 1365.65054 · doi:10.1080/10556788.2016.1146269 |
| [9] | Heimbach, P.; Bugnion, V., greenland ice-sheet volume sensitivity to basal, surface and initial conditions derived from an adjoint model, Ann. Glaciol., 50, 67-80, (2009) · doi:10.3189/172756409789624256 |
| [10] | Larour, E.; Utke, J.; Csatho, B.; Schenk, A.; Seroussi, H.; Morlighem, M.; Rignot, E.; Schlegel, N.; Khazendar, A., inferred basal friction and surface mass balance of the northeast greenland ice stream using data assimilation of icesat (ice cloud and land elevation satellite) surface altimetry and ISSM (ice sheet system model), Cryosphere, 8, 2335-2351, (2014) · doi:10.5194/tc-8-2335-2014 |
| [11] | MacAyeal, D. R., large-scale ice flow over a viscous basal sediment: theory and application to ice stream B, antarctica, J. Geophys. Res., 94, 4071-4087, (1989) · doi:10.1029/JB094iB04p04071 |
| [12] | MacAyeal, D. R., A tutorial on the use of control methods in ice-sheet modeling, J. Glaciol., 39, 91-98, (1993) · doi:10.3189/S0022143000015744 |
| [13] | Morlighem, M.; Rignot, E.; Seroussi, H.; Larour, E.; Ben Dhia, H.; Aubry, D., spatial patterns of basal drag inferred using control methods from a full-Stokes and simpler models for pine island glacier, west antarctica, Geophys. Res. Lett., 37, 1-6, (2010) · doi:10.1029/2010GL043853 |
| [14] | Morlighem, M.; Seroussi, H.; Larour, E.; Rignot, E., inversion of basal friction in antarctica using exact and incomplete adjoints of a higher-order model, J. Geophys. Res., 118, 1746-1753, (2013) · doi:10.1002/jgrf.20125 |
| [15] | TAPENADE for C, in Advances in Automatic DifferentiationLecture Notes in Computational Science and Engineering, selected papers from AD2008, August 2008, Bonn, Springer2008199210 |
| [16] | Utke, J.; Naumann, U.; Fagan, M.; Tallent, N.; Strout, M.; Heimbach, P.; Hill, C.; Wunsch, C., openad/F: A modular open-source tool for automatic differentiation of Fortran codes, ACM Trans. Math. Softw., 34, 18:1-18:36, (2008) · Zbl 1291.65140 · doi:10.1145/1377596.1377598 |
| [17] | Utke, J.; Hascoët, L.; Heimbach, P.; Hill, C.; Hovland, P.; Naumann, U., toward adjoinable MPI, (2009) |
| [18] | Getting started with ADOL-CCombinatorial Scientific ComputingChapter 7, Chapman-Hall CRC Computational Science2012181202 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
