Abstract. In this paper we consider the application of least squares principles to the approximate solution of the Stokes equations cast into a first-order velocity-vorticity-pressure system. Among the most attractive features of the resulting methods are that the choice of approximating spaces is not subject to the LBB condition and a single continuous piecewise polynomial space can be used for the approximation of all unknowns, that the resulting discretized problems involve only symmetric, positive definite systems of algebraic equations, that no artificial boundary conditions for the vorticity need be devised, and that accurate approximations are obtained for all variables, including the vorticity. Here we study two classes of least squares methods for the velocity-vorticitypressure equations. The first one uses norms prescribed by the a priori estimates of Agmon, Doughs, and Nirenberg and can be analyzed in a completely standard manner. However, conforming discretizations of these methods require C1 continuity of the finite element spaces, thus negating the advantages of the velocity-vorticity-pressure formulation. The second class uses weighted L2-norms of the residuals to circumvent this flaw. For properly chosen meshdependent weights, it is shown that the approximations to the solutions of the Stokes equations are of optimal order. The results of some computational experiments are also provided; these illustrate, among other things, the necessity of introducing the weights.