Multigrid methods for two first-order system least squares (FOSLS) using bilinear finite elements are developed for the pure traction problem of planar linear elasticity. They are two-stage algorithms that first solve for the gradients of displacement, then for the displacement itself. In this paper, concentration is given on solving for the gradients of displacement only. Numerical results show that the convergences are uniform even as the material becomes nearly incompressible. Computations for convergence rates are included.