We study Hardy spaces on C 1 and Lipschitz domains in Riemannian manifolds. Hardy spaces, originally introduced in 1920 in complex analysis setting, are invaluable tool in harmonic analysis. For this reason these spaces have been studied extensively by many authors.Our main result is an equivalence theorem proving that the definition of Hardy spaces by conjugate harmonic functions is equivalent to the atomic definition of these spaces. We establish this theorem in any dimension if the domain is C 1 , in case of a Lipschitz domain the result holds if dim M ≤ 3. The remaining cases for Lipschitz domains remain open. This result is a nontrivial generalization of flat (R n ) equivalence theorems due to Fefferman, Stein, Dahlberg and others.The material presented here required to develop potential theory approach for C 1 domains on Riemannian manifolds in the spirit of earlier works by Fabes, Jodeit and Rivière and recent results by Mitrea and Taylor. In particular, the first part of this work is of interest in itself, since we consider the boundary value problems for the Laplace-Beltrami operator. We prove that both Dirichlet and Neumann problem for Laplace-Beltrami equation are solvable for any given boundary data in L p (∂Ω), where 1 < p < ∞. Same remains true in Hardy spaces p (∂Ω) for (n − 1)/n < p ≤ 1.In the whole work we work with Riemannian metric g with smallest possible regularity. In particular, mentioned results for the Laplace-Beltrami equation require Hölder class regularity of the metric tensor; the equivalence theorem requires g in C 1,1 .