In this paper, we give a sharp estimate on the cardinality of the set generating the convex hull containing the image of harmonic maps with polynomial growth rate on a certain class of manifolds into a Cartan-Hadamard manifold with sectional curvature bounded by two negative constants. We also describe the asymptotic behavior of harmonic maps on a complete Riemannian manifold into a regular ball in terms of massive subsets, in the case when the space of bounded harmonic functions on the manifold is finite dimensional.