The goal of this paper is to investigate the ability of Amberlite XAD-7 (AXAD-7) resin impregnated with CYANEX-272 (di-2,4,4-trimethylpentyl phosphonic acid) to remove cerium (Ce(ΙΙΙ)) and lanthanum (La(ΙΙΙ)) ions from aqueous solutions in the batch scheme. The prepared adsorbent material was determined utilizing FTIR, SEM–EDX, and BET methods. The impact of three individual process variable factors involving feed solution pH (2–6), adsorbent dose (0.05–0.65), and process temperature (15–55 °C) on the simultaneous removal of Ce(ΙΙΙ) and La(ΙΙΙ) ions was evaluated via response surface methodology (RSM) according to the central composite design (CCD). The modeling of Ce(ΙΙΙ) and La(ΙΙΙ) ions adsorption was performed using the quadratic model and was evaluated using a coefficient of determination for both ions. The optimization data revealed that the adsorption amount of Ce(ΙΙΙ) and La(ΙΙΙ) ions removal under optimal conditions were 99.75% and 69.98%, respectively. Equilibrium and kinetic investigations were also conducted to define the removal performance of the calculated adsorbent for Ce(ΙΙΙ) and La(ΙΙΙ) ions removal. Various isotherms models such as Langmuir, Freundlich, Temkin, and Sips were examined at 25 °C to analyze the equilibrium isotherm data. The data revealed that the Sips approach is compatible with the experimental data. The highest adsorption capacity of the resin for Ce(ΙΙΙ) and La(ΙΙΙ) ions were 11.873 mg g−1 and 7.324 mg g−1, correspondingly. The kinetic study of the Ce(ΙΙΙ) and La(ΙΙΙ) adsorption process was conducted via pseudo-first-order, pseudo-second-order, and intraparticle diffusion models(IDMs). Based on the data obtained, kinetic data were fitted well to a pseudo-second-order rate correlation. According to the obtained results, the (AXAD-7) resin impregnated with CYANEX-272 performed well in removing both Ce(ΙΙΙ) and La(ΙΙΙ) ions from aqueous solutions with well stability during several adsorption–desorption cycles and well regeneration and excellent metallic ions recovery.