Comparison of numerical schemes for the solution of the advective age equation in ice sheetsR. GREVE, Y. WANG and B. MÜGGE AbstractA onedimensional model problem for the computation of the age field in ice sheets, which is of great importance for dating deep ice cores, is considered. The corresponding partial differential equation (PDE) is of purely advective (hyperbolic) type, which is notoriously difficult to solve numerically. By integrating the PDE over a spacetime element in the sense of a finitevolume approach, a general difference equation is constructed from which a hierarchy of solution schemes can be derived. Iteration rules are given explicitly for central differences, first, second and thirdorder (QUICK) upstreaming as well as modified TVD LaxFriedrichs schemes. The performance of these schemes in terms of convergence and accuracy is discussed. It turns out that secondorder upstreaming and the modified TVD LaxFriedrichs scheme with Minmod slope limiter are most suitable for numerical age computations in icesheet models. Annals of Glaciology, 35, 487494 (2002). 
