The Freundlich model and the Constant Capacitance Model (CCM) were used to describe silicon (Si) and phosphorus (P) sorption, both individually and for binary P-Si systems, on two Andisols with different chemical properties: Freire soil (FS) and Piedras Negras soil (PNS). Silicon sorption kinetics were examined through the Elovich equation, revealing that the initial sorption rate was 16 times greater in PNS. The Freundlich equation provides a good fit to the sorption data for both Andisols. When compared with FS, larger Si sorption capacity and lower Si affinity for the surface sites were observed in PNS; nevertheless, Si sorption decreased in both soils as P sorption increased. Slight reductions in P sorption capacity due to the presence of Si were found, whereas there was no apparent effect on P bonding intensity. The CCM was able to describe Si adsorption, and potentiometric titrations support that Si seems to be specifically sorbed mainly onto sites of negative charge. Comparable logK(Si)(int) values were obtained for both soils, indicating that Si was bound on similar sites. Phosphorus sorption was well described by the CCM, and logK(P)(int) denoted strong interactions of P with the surface sites. For binary systems, logK(P)(int) did not vary with increasing Si concentration; comparatively, logK(S1)(int) scarcely decreased with increasing P concentration in PNS, but a 28% reduction was found in FS at the highest initial P concentration. (C) 2015 Elsevier Ltd. All rights reserved.