In this paper, we present a multivariate spectral gradient projection method for nonlinear monotone equations with convex constraints, which can be viewed as an extension of multivariate spectral gradient method for solving unconstrained optimization problems. The proposed method does not need the computation of the derivative as well as the solution of some linear equations. Under some suitable conditions, we can establish its global convergence results. Preliminary numerical results show that the proposed method is efficient and promising.