Abstract. We present the detailed construction of a manufactured analytical solution to time-dependent and steady-state isothermal full-Stokes ice sheet problems. The solutions are constructed for two-dimensional flowline and three-dimensional full-Stokes ice sheet models with variable viscosity. The construction is done by choosing for the specified ice surface and bed a velocity distribution that satisfies both mass conservation and the kinematic boundary conditions. Then a compensatory stress term in the conservation of momentum equations and their boundary conditions is calculated to make the chosen velocity distributions as well as the chosen pressure field into exact solutions. By substituting different ice surface and bed geometry formulas into the derived solution formulas, analytical solutions for different geometries can be constructed.

The boundary conditions can be specified as essential Dirichlet conditions or as periodic boundary conditions. By changing a parameter value, the analytical solutions allow investigation of algorithms for a different range of aspect ratios as well as for different, frozen or sliding, basal conditions. The analytical solutions can also be used to estimate the numerical error of the method in the case when the effects of the boundary conditions are eliminated, that is, when the exact solution values are specified as inflow and outflow boundary conditions.