Comparison of Numerical Schemes for the Solution of the Ice-Thickness Equation in a Dynamic/Thermodynamic Ice-Sheet Model
✍ Scribed by Ralf Greve; Reinhard Calov
- Publisher
- Elsevier Science
- Year
- 2002
- Tongue
- English
- Weight
- 247 KB
- Volume
- 179
- Category
- Article
- ISSN
- 0021-9991
No coin nor oath required. For personal study only.
✦ Synopsis
A general finite-difference marching scheme for the numerical solution of the ice-thickness equation in ice sheets is considered. From this scheme, a variety of explicit, ADI, implicit and over-implicit methods can be derived. These methods are compared for stability and accuracy within the dynamic/thermodynamic ice-sheet model SICOPOLIS for two different problems: (i) a simple axi-symmetric steadystate ice sheet which rests on a flat bedrock, and (ii) the time-dependent paleoglaciation of the northern hemisphere. As expected, over-implicit methods turn out to be most stable. For the simple problem, all schemes provide a good accuracy, whereas for the northern hemisphere simulations, the accuracy of the over-implicit scheme is not satisfactory, so that the implicit technique without over-weighing appears favorable for this application.