Let SS be a surface with pg(S)=q(S)=0pg(S)=q(S)=0 and endowed with a very ample line bundle 𝒪S(h) such that h1(S,𝒪S(h))=0. We show that S supports special (often stable) Ulrich bundles of rank 2, extending a recent result by A. Beauville. Moreover, we show that such an S supports families of dimension p of pairwise non-isomorphic, indecomposable, Ulrich bundles for arbitrary large p except for very few cases. We also show that the same is true for each linearly normal non-special surface with pg(S)=q(S)=0 in ℙ4 of degree at least 4, Enriques surface and anticanonical rational surface.