In this paper, we investigate the decay properties of an axisymmetric D-solutions to stationary incompressible Navier-Stokes systems in $mathbb{R}^3$. We obtain the optimal decay rate $|{bf u}(x)|leq frac{C}{|x|+1}$ for axisymmetric flows without swirl. Furthermore, we find a dichotomy for the decay rates of the swirl component $u_{theta}$, that is, either $O(frac{1}{r+1})leq |u_{theta}(r,z)|leq frac{Clog(r+1)}{(r+1)^{1/2}}$ or $|u_{theta}(r,z)|leq frac{C r}{(rho+1)^3}$, where $rho=sqrt{r^2+z^2}$. In the latter case, we can further deduce that the other two components of the velocity field also attain the optimal decay rates: $|u_r(r,z)|+ |u_{z}(r,z)|leq frac{C}{rho+1}$. We do not require any small assumptions on the forcing term.