The Non-Ideal Competitive Adsorption model (NICA) developed from the Langmuir adsorption theory plus the Freundlich empirical equation has some advantages in describing adsorption of ions on solid surface in solutions. The reason is that the NICA model separates the non-ideality of adsorption into the intrinsic chemical heterogeneity of the surface and the ion-specific non-ideality. We believe that the NICA model is more suitable to describe ion adsorption in complex systems, like soils, which have multiple ions, highly heterogeneous surfaces and a variety of sites on the particle surfaces to adsorb ions. Although it has been successfully applied in the binding of protons and metal ions to humic substances, its application for the adsorption of anions on the positively charged surfaces has not been documented. The purpose of this study was to find a suitable way to adapt the NICA model for describing hydroxyl and phosphate adsorption, which involves major assumptions and redefination of surface reaction. Results show that by considering 2 types of surface, the NICA model can excellently fit (R²>0.99) the hydroxyl data obtained from -11 to -4 of log [OH], which corresponds to the pH from 3 to 10. By using the parameters generated from hydroxyl adsorption, including adsorption maxima (Q_max), binding strength (K_OH) and non-ideality (m), the NICA model gave a remarkable goodness of fit (R²>0.98) for the P adsorption data obtained at different pH levels. The successful application of P adsorption combined with hydroxyl adsorption provides a promising tool for analyzing competitive adsorption among anions in soils.