Abstract. In this paper we give a sub-supersolution method for nonlinear elliptic singular systems with quadratic gradient whose model system is the followingwhere Ω is a smooth bounded domain of R N (N ≥ 3), β, µ ≥ 0, 0 < α, γ < 1 and regular f 1 , f 2 functions. Moreover, we apply it to prove existence of solution for some systems, including the classical Lotka-Volterra models with gradient terms. Specifically, we study the competition and the symbiotic LoktaVolterra systems.