  In this paper, we investigate the decay properties of axially symmetric
solutions to the steady incompressible magnetohydrodynamic equations in
$\mathbb{R}^3$ with finite Dirichlet integral. We first derive the decay rates
of general D-solutions to the axisymmetric MHD equations. In the special case
where the magnetic field only has the swirl component ${\textbf {h}}(r,z)=
h_{\th}(r,z) {\bf e_{\th}}$, we obtain better decay rates. The last result
examines the decay rates along the axis $Oz$ also within the special class of
D-solutions with only swirl magnetic field. The main tool in this paper is the
combination of the scaling argument, the \textit{Brezis-Gallouet inequality}
and the weighted energy estimate.
