In this work, we investigate the existence and nonexistence of nonnegative solutions to a class of nonlocal elliptic systems set in a bounded open subset of R N , where the gradients of the unknowns act as source terms (see (S) below). Our approach can be also used to treat other nonlinear systems with different structures. This work extends previous results obtained in the local case by the fourth author and his coworkers, and points to significant differences between the local and the nonlocal cases.