In this paper, we analyze the nonlinear Stokes equations which is appeared in glaciology problems with least-squares finite element method. The modified Picard iteration is used to linearize the given ...