Initialization of ice-sheet forecasts viewed as an inverse Robin problem
As simulations of 21st-century climate start to include components with longer timescales, such as ice sheets, the initial conditions for those components will become critical to the forecast. This paper describes an algorithm for specifying the initial state of an ice-sheet model, given spatially c...
Main Authors: | , |
---|---|
Format: | Article in Journal/Newspaper |
Language: | unknown |
Published: |
International Glaciological Society
2010
|
Subjects: | |
Online Access: | http://nora.nerc.ac.uk/id/eprint/10838/ http://www.igsoc.org/journal/56/197/t09J101.pdf |
id |
ftnerc:oai:nora.nerc.ac.uk:10838 |
---|---|
record_format |
openpolar |
spelling |
ftnerc:oai:nora.nerc.ac.uk:10838 2023-05-15T16:40:30+02:00 Initialization of ice-sheet forecasts viewed as an inverse Robin problem Arthern, Robert J. Gudmundsson, Hilmar 2010 http://nora.nerc.ac.uk/id/eprint/10838/ http://www.igsoc.org/journal/56/197/t09J101.pdf unknown International Glaciological Society Arthern, Robert J. orcid:0000-0002-3762-8219 Gudmundsson, Hilmar orcid:0000-0003-4236-5369 . 2010 Initialization of ice-sheet forecasts viewed as an inverse Robin problem. Journal of Glaciology, 56 (197). 527-533. Meteorology and Climatology Glaciology Publication - Article PeerReviewed 2010 ftnerc 2023-02-04T19:26:58Z As simulations of 21st-century climate start to include components with longer timescales, such as ice sheets, the initial conditions for those components will become critical to the forecast. This paper describes an algorithm for specifying the initial state of an ice-sheet model, given spatially continuous observations of the surface elevation, the velocity at the surface and the thickness of the ice. The algorithm can be viewed as an inverse procedure to solve for the viscosity or the basal drag coefficient. It applies to incompressible Stokes flow over an impenetrable boundary, and is based upon techniques used in electric impedance tomography; in particular, the minimization of a type of cost function proposed by Kohn and Vogelius. The algorithm can be implemented numerically using only the forward solution of the Stokes equations, with no need to develop a separate adjoint model. The only requirement placed upon the numerical Stokes solver is that boundary conditions of Dirichlet, Neumann and Robin types can be implemented. As an illustrative example, the algorithm is applied to shear flow down an impenetrable inclined plane. A fully three-dimensional test case using a commercially available solver for the Stokes equations is also presented. Article in Journal/Newspaper Ice Sheet Journal of Glaciology Natural Environment Research Council: NERC Open Research Archive |
institution |
Open Polar |
collection |
Natural Environment Research Council: NERC Open Research Archive |
op_collection_id |
ftnerc |
language |
unknown |
topic |
Meteorology and Climatology Glaciology |
spellingShingle |
Meteorology and Climatology Glaciology Arthern, Robert J. Gudmundsson, Hilmar Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
topic_facet |
Meteorology and Climatology Glaciology |
description |
As simulations of 21st-century climate start to include components with longer timescales, such as ice sheets, the initial conditions for those components will become critical to the forecast. This paper describes an algorithm for specifying the initial state of an ice-sheet model, given spatially continuous observations of the surface elevation, the velocity at the surface and the thickness of the ice. The algorithm can be viewed as an inverse procedure to solve for the viscosity or the basal drag coefficient. It applies to incompressible Stokes flow over an impenetrable boundary, and is based upon techniques used in electric impedance tomography; in particular, the minimization of a type of cost function proposed by Kohn and Vogelius. The algorithm can be implemented numerically using only the forward solution of the Stokes equations, with no need to develop a separate adjoint model. The only requirement placed upon the numerical Stokes solver is that boundary conditions of Dirichlet, Neumann and Robin types can be implemented. As an illustrative example, the algorithm is applied to shear flow down an impenetrable inclined plane. A fully three-dimensional test case using a commercially available solver for the Stokes equations is also presented. |
format |
Article in Journal/Newspaper |
author |
Arthern, Robert J. Gudmundsson, Hilmar |
author_facet |
Arthern, Robert J. Gudmundsson, Hilmar |
author_sort |
Arthern, Robert J. |
title |
Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
title_short |
Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
title_full |
Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
title_fullStr |
Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
title_full_unstemmed |
Initialization of ice-sheet forecasts viewed as an inverse Robin problem |
title_sort |
initialization of ice-sheet forecasts viewed as an inverse robin problem |
publisher |
International Glaciological Society |
publishDate |
2010 |
url |
http://nora.nerc.ac.uk/id/eprint/10838/ http://www.igsoc.org/journal/56/197/t09J101.pdf |
genre |
Ice Sheet Journal of Glaciology |
genre_facet |
Ice Sheet Journal of Glaciology |
op_relation |
Arthern, Robert J. orcid:0000-0002-3762-8219 Gudmundsson, Hilmar orcid:0000-0003-4236-5369 . 2010 Initialization of ice-sheet forecasts viewed as an inverse Robin problem. Journal of Glaciology, 56 (197). 527-533. |
_version_ |
1766030909745659904 |