Volumes and Issues  Contents of Issue 1  
Geosci. Model Dev. Discuss., 2, 81-158, 2009
www.geosci-model-dev-discuss.net/2/81/2009/
doi:10.5194/gmdd-2-81-2009
© Author(s) 2009. This work is distributed
under the Creative Commons Attribution 3.0 License.


Derivation of a numerical solution of the 3D coupled velocity field for an ice sheet – ice shelf system, incorporating both full and approximate stress solutions

T. J. Reerink, R. S. W. van de Wal, and P.-P. Borsboom
Institute for Marine and Atmospheric Research Utrecht, Utrecht University, The Netherlands

Abstract. To overcome the mechanical coupling of an ice sheet with an ice shelf, one single set of velocity equations is presented covering both the sheet and the shelf. This set is obtained by applying shared sheet-shelf approximations. The hydrostatic approximation and a constant density are the only approximations that are applied to the full-Stokes momentum equations. The remaining stress terms from the momentum equations and the stress terms from the usual ice-flow law are multiplied by coefficients which can be put to zero or one, facilitating several stress approximations per domain within one model. In addition we derived a matrix format for the discretized set of the fully coupled velocity equations on a three-dimensional vertically scaled grid, in which all linear derivative terms are treated implicitly. The compact vector format of this sparse matrix equation is developed, including the boundary conditions.

Discussion Paper (PDF, 1025 KB)   Interactive Discussion (Closed, 5 Comments)   Publication in GMD not foreseen   

Citation: Reerink, T. J., van de Wal, R. S. W., and Borsboom, P.-P.: Derivation of a numerical solution of the 3D coupled velocity field for an ice sheet – ice shelf system, incorporating both full and approximate stress solutions, Geosci. Model Dev. Discuss., 2, 81-158, doi:10.5194/gmdd-2-81-2009, 2009.   Bibtex   EndNote   Reference Manager    XML