In this study, iron oxide (Fe3O4) was coated with ZrO2, and doped with three rare earth elements((Y/La/Ce), and a multi-staged rare earth doped zirconia adsorbent was prepared by using uniform design U14, Response Surface methodology, and orthogonal design, to remove As3+ and As5+ from the aqueous solution. Based on the results of TEM, EDS, XRD, FTIR, and N2-adsorption desorption test, the best molar ratio of Fe3O4:TMAOH:Zirconium butoxide:Y:La:Ce was selected as 1:12:11:1:0.02:0.08. The specific surface area and porosity was 263 m2/g, and 0.156 cm3/g, respectively. The isothermal curves and fitting equation parameters show that Langmuir model, and Redlich Peterson model fitted well. As per calculations of the Langmuir model, the highest adsorption capacities for As3+ and As5+ ions were recorded as 68.33 mg/g, 84.23 mg/g, respectively. The fitting curves and equations of the kinetic models favors the quasi second order kinetic model. Material regeneration was very effective, and even in the last cycle the regeneration capacities of both As3+ and As5+ were 75.15%, and 77.59%, respectively. Adsorption and regeneration results suggest that adsorbent has easy synthesis method, and reusable, so it can be used as a potential adsorbent for the removal of arsenic from aqueous solution.