  In this paper, we investigate subelliptic harmonic maps with potential from
noncompact complete sub-Riemannian manifolds corresponding to totally geodesic
Riemannian foliations. Under some suitable conditions, we give the gradient
estimates of these maps and establish a Liouville type result.
