[en] Prediction of radionuclides migration through a geosphere is of fundamental interest in order to evaluate the safety of an underground radwastes repository. The authors have studied uranyl and europium (III) ion sorption mechanisms onto three phosphate compounds: ZrP₂O₇, Zr₂O(PO₄)₂, and Th₄(PO₄)₄P₂O₇. A spectroscopic investigation which has been previously reported has allowed one to determine all species involved in the sorption processes under study. This paper presents the modeling of experimental retention data using FITEQLv3.2 code (constant capacitance model). All sorption isotherms were successfully simulated with respect to the structural constraints. The authors have shown that uranyl sorption onto the solids under study is less influenced by electrostatic interactions between aqueous species and charged on the surface than is that for europium(III) ion. Because many experimental constraints, obtained from independent spectroscopic techniques, have been taken into account for the fitting procedure, the sorption constant values can be determined accurately