Effects of first-order stress gradients in an ice sheet evaluated by a three-dimensional thermomechanical coupled model

A three-dimensional thermomechanically coupled ice-sheet model including calculation of the first-order stress gradients is developed to perform numerical studies on the effects of stress components neglected in the zeroth-order shallow-ice approximation. Steady-state solutions are obtained with both a shallow-ice and a first-order ice-sheet model for several idealistic, radially symmetric ice-sheet configurations. The results show that the effects of the normal deviatoric-stress gradients on ice thickness are generally small, but the influence on basal temperatures and age profiles at an ice divide is significant.