We consider operators of the form L = n i=1 X 2 i + X0 in a bounded domain of R p where X0, X1, . . . , Xn are nonsmooth Hörmander's vector fields of step r such that the highest order commutators are only Hölder continuous. Applying Levi's parametrix method we construct a local fundamental solution γ for L and provide growth estimates for γ and its first derivatives with respect to the vector fields. Requiring the existence of one more derivative of the coefficients we prove that γ also possesses second derivatives, and we deduce the local solvability of L, constructing, by means of γ, a solution to Lu = f with Hölder continuous f . We also prove C 2,α X,loc estimates on this solution.