SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD
We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD. The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective fu...
Published in: | Geoscientific Model Development |
---|---|
Main Authors: | , , , |
Format: | Article in Journal/Newspaper |
Language: | English |
Published: |
Copernicus Publications
2020
|
Subjects: | |
Online Access: | https://doi.org/10.5194/gmd-13-1845-2020 https://doaj.org/article/11cf4086df184848a310540e865ccc23 |
id |
ftdoajarticles:oai:doaj.org/article:11cf4086df184848a310540e865ccc23 |
---|---|
record_format |
openpolar |
spelling |
ftdoajarticles:oai:doaj.org/article:11cf4086df184848a310540e865ccc23 2023-05-15T13:52:54+02:00 SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD L. C. Logan S. H. K. Narayanan R. Greve P. Heimbach 2020-04-01T00:00:00Z https://doi.org/10.5194/gmd-13-1845-2020 https://doaj.org/article/11cf4086df184848a310540e865ccc23 EN eng Copernicus Publications https://www.geosci-model-dev.net/13/1845/2020/gmd-13-1845-2020.pdf https://doaj.org/toc/1991-959X https://doaj.org/toc/1991-9603 doi:10.5194/gmd-13-1845-2020 1991-959X 1991-9603 https://doaj.org/article/11cf4086df184848a310540e865ccc23 Geoscientific Model Development, Vol 13, Pp 1845-1864 (2020) Geology QE1-996.5 article 2020 ftdoajarticles https://doi.org/10.5194/gmd-13-1845-2020 2022-12-31T04:00:41Z We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD. The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters. Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances: (i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community; (ii) the AD tool used, OpenAD, is an open-source tool; (iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations. A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large. As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet. Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that ... Article in Journal/Newspaper Antarc* Antarctic Antarctica Greenland Ice Sheet Ice Shelves Queen Mary land Directory of Open Access Journals: DOAJ Articles Antarctic Austral Greenland Queen Mary Land ENVELOPE(96.000,96.000,-68.000,-68.000) Geoscientific Model Development 13 4 1845 1864 |
institution |
Open Polar |
collection |
Directory of Open Access Journals: DOAJ Articles |
op_collection_id |
ftdoajarticles |
language |
English |
topic |
Geology QE1-996.5 |
spellingShingle |
Geology QE1-996.5 L. C. Logan S. H. K. Narayanan R. Greve P. Heimbach SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
topic_facet |
Geology QE1-996.5 |
description |
We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD. The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters. Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances: (i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community; (ii) the AD tool used, OpenAD, is an open-source tool; (iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations. A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large. As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet. Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that ... |
format |
Article in Journal/Newspaper |
author |
L. C. Logan S. H. K. Narayanan R. Greve P. Heimbach |
author_facet |
L. C. Logan S. H. K. Narayanan R. Greve P. Heimbach |
author_sort |
L. C. Logan |
title |
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
title_short |
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
title_full |
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
title_fullStr |
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
title_full_unstemmed |
SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD |
title_sort |
sicopolis-ad v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool openad |
publisher |
Copernicus Publications |
publishDate |
2020 |
url |
https://doi.org/10.5194/gmd-13-1845-2020 https://doaj.org/article/11cf4086df184848a310540e865ccc23 |
long_lat |
ENVELOPE(96.000,96.000,-68.000,-68.000) |
geographic |
Antarctic Austral Greenland Queen Mary Land |
geographic_facet |
Antarctic Austral Greenland Queen Mary Land |
genre |
Antarc* Antarctic Antarctica Greenland Ice Sheet Ice Shelves Queen Mary land |
genre_facet |
Antarc* Antarctic Antarctica Greenland Ice Sheet Ice Shelves Queen Mary land |
op_source |
Geoscientific Model Development, Vol 13, Pp 1845-1864 (2020) |
op_relation |
https://www.geosci-model-dev.net/13/1845/2020/gmd-13-1845-2020.pdf https://doaj.org/toc/1991-959X https://doaj.org/toc/1991-9603 doi:10.5194/gmd-13-1845-2020 1991-959X 1991-9603 https://doaj.org/article/11cf4086df184848a310540e865ccc23 |
op_doi |
https://doi.org/10.5194/gmd-13-1845-2020 |
container_title |
Geoscientific Model Development |
container_volume |
13 |
container_issue |
4 |
container_start_page |
1845 |
op_container_end_page |
1864 |
_version_ |
1766257751550328832 |