A one-dimensional model problem for 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 space–time element in the sense of a finite-volume 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 third-order (QUICK) upstreaming as well as modified TVD Lax–Friedrichs schemes (TVDLFs). the performance of these schemes in terms of convergence and accuracy is discussed. Second-order upstreaming, themodifiedTVDLF scheme with Minmod slope limiter and, with limitations of the accuracy directly at the base, first-order upstreaming prove to be the most suitable for numerical age computations in ice-sheet models.