Parameter and state estimation with a time-dependent adjoint marine ice sheet model
To date, assimilation of observations into large-scale ice models has consisted predominantly of time-independent inversions of surface velocities for basal traction, bed elevation, or ice stiffness, and has relied primarily on analytically derived adjoints of glaciological stress balance models. To...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2013-11-01
|
| Series: | The Cryosphere |
| Online Access: | http://www.the-cryosphere.net/7/1659/2013/tc-7-1659-2013.pdf |
| _version_ | 1818540557607108608 |
|---|---|
| author | D. N. Goldberg P. Heimbach |
| author_facet | D. N. Goldberg P. Heimbach |
| author_sort | D. N. Goldberg |
| collection | DOAJ |
| description | To date, assimilation of observations into large-scale ice models has
consisted predominantly of time-independent inversions of surface velocities
for basal traction, bed elevation, or ice stiffness, and has relied primarily
on analytically derived adjoints of glaciological stress balance models. To
overcome limitations of such "snapshot" inversions – i.e., their inability to
assimilate time-dependent data for the purpose of constraining transient flow
states, or to produce initial states with minimum artificial drift and
suitable for time-dependent simulations – we have developed an adjoint of a
time-dependent parallel glaciological flow model. The model implements a
hybrid shallow shelf–shallow ice stress balance, solves the continuity
equation for ice thickness evolution, and can represent the floating,
fast-sliding, and frozen bed regimes of a marine ice sheet. The adjoint is
generated by a combination of analytic methods and the use of algorithmic
differentiation (AD) software. Several experiments are carried out with
idealized geometries and synthetic observations, including inversion of
time-dependent surface elevations for past thicknesses, and simultaneous
retrieval of basal traction and topography from surface data. Flexible
generation of the adjoint for a range of independent uncertain variables is
exemplified through sensitivity calculations of grounded ice volume to
changes in basal melting of floating and basal sliding of grounded ice. The
results are encouraging and suggest the feasibility, using real observations,
of improved ice sheet state estimation and comprehensive transient
sensitivity assessments. |
| first_indexed | 2024-12-11T21:56:57Z |
| format | Article |
| id | doaj.art-bddcb24761e94074a5b03aae1380441c |
| institution | Directory Open Access Journal |
| issn | 1994-0416 1994-0424 |
| language | English |
| last_indexed | 2024-12-11T21:56:57Z |
| publishDate | 2013-11-01 |
| publisher | Copernicus Publications |
| record_format | Article |
| series | The Cryosphere |
| spelling | doaj.art-bddcb24761e94074a5b03aae1380441c2022-12-22T00:49:16ZengCopernicus PublicationsThe Cryosphere1994-04161994-04242013-11-01761659167810.5194/tc-7-1659-2013Parameter and state estimation with a time-dependent adjoint marine ice sheet modelD. N. Goldberg0P. Heimbach1Earth, Atmospheric and Planetary Sciences, M.I.T., Cambridge, USA.Earth, Atmospheric and Planetary Sciences, M.I.T., Cambridge, USA.To date, assimilation of observations into large-scale ice models has consisted predominantly of time-independent inversions of surface velocities for basal traction, bed elevation, or ice stiffness, and has relied primarily on analytically derived adjoints of glaciological stress balance models. To overcome limitations of such "snapshot" inversions – i.e., their inability to assimilate time-dependent data for the purpose of constraining transient flow states, or to produce initial states with minimum artificial drift and suitable for time-dependent simulations – we have developed an adjoint of a time-dependent parallel glaciological flow model. The model implements a hybrid shallow shelf–shallow ice stress balance, solves the continuity equation for ice thickness evolution, and can represent the floating, fast-sliding, and frozen bed regimes of a marine ice sheet. The adjoint is generated by a combination of analytic methods and the use of algorithmic differentiation (AD) software. Several experiments are carried out with idealized geometries and synthetic observations, including inversion of time-dependent surface elevations for past thicknesses, and simultaneous retrieval of basal traction and topography from surface data. Flexible generation of the adjoint for a range of independent uncertain variables is exemplified through sensitivity calculations of grounded ice volume to changes in basal melting of floating and basal sliding of grounded ice. The results are encouraging and suggest the feasibility, using real observations, of improved ice sheet state estimation and comprehensive transient sensitivity assessments.http://www.the-cryosphere.net/7/1659/2013/tc-7-1659-2013.pdf |
| spellingShingle | D. N. Goldberg P. Heimbach Parameter and state estimation with a time-dependent adjoint marine ice sheet model The Cryosphere |
| title | Parameter and state estimation with a time-dependent adjoint marine ice sheet model |
| title_full | Parameter and state estimation with a time-dependent adjoint marine ice sheet model |
| title_fullStr | Parameter and state estimation with a time-dependent adjoint marine ice sheet model |
| title_full_unstemmed | Parameter and state estimation with a time-dependent adjoint marine ice sheet model |
| title_short | Parameter and state estimation with a time-dependent adjoint marine ice sheet model |
| title_sort | parameter and state estimation with a time dependent adjoint marine ice sheet model |
| url | http://www.the-cryosphere.net/7/1659/2013/tc-7-1659-2013.pdf |
| work_keys_str_mv | AT dngoldberg parameterandstateestimationwithatimedependentadjointmarineicesheetmodel AT pheimbach parameterandstateestimationwithatimedependentadjointmarineicesheetmodel |