The objective of this work was to investigate glyphosate adsorption by soils and its relationship with unoccupied binding sites for phosphate adsorption. Soil samples of three Chilean soils series Valdivia (Andisol), Clarillo (Inceptisol) and Chicureo (Vertisol) were incubated with different herbicide concentrations. Glyphosate remaining in solution was determined by adjusting a HPLC method with a UV detector. Experimental maximum adsorption capacity were 15,000, 14,300 and 4,700 µg g -1 for Valdivia, Clarillo, and Chicureo soils, respectively. Linear, Freundlich, and Langmuir models were used to describe glyphosate adsorption. Isotherms describing glyphosate adsorption differed among soils. Maximum adjusted adsorption capacity with the Langmuir model was 231,884, 17,874 and 5,670 µg g -1 for Valdivia, Clarillo, and Chicureo soils, respectively. Glyphosate adsorption on the Valdivia soil showed a linear behavior at the range of concentrations used and none of the adjusted models became asymptotic. The high glyphosate adsorption capacity of the Valdivia soil was probably a result of its high exchangeable Al, extractable Fe, and alophan and imogolite clay type. Adsorption was very much related to phosphate dynamics in the Valdivia soil, which showed the larger unoccupied phosphate binding sites. However relationship between unoccupied phosphate binding sites and glyphosate adsorption in the other two soils (Clarillo and Chicureo) was not clear.