Let $A=-(\nabla-i\vec{a})^2+V$ be a magnetic Schrödinger operator acting on$L$^{2}($R$^{$n$}),$n$≥1, where $\vec{a}=(a_1,\cdots,a_n)\in L^2_{\rm loc}$ and 0≤$V$∈$L$^{1}_{loc}. Following [1], we define, by means of the area integral function, a Hardy space$H$^{1}_{$A$}associated with$A$. We show that Riesz transforms (∂/∂$x$_{$k$}-$i$$a$_{$k$})$A$^{-1/2}associated with$A$, $k=1,\cdots,n$ , are bounded from the Hardy space$H$^{1}_{$A$}into$L$^{1}. By interpolation, the Riesz transforms are bounded on$L$^{$p$}for all 1<$p$≤2.