The adsorption isotherms of acetonitrile, ethanol, 2-propanol, and THF were measured using frontal analysis on six columns packed with octadecyl RP stationary phase. The effect of the bonding density of the end-capped octadecyl bonded phase on the adsorption properties was measured. Adsorption isotherm data were collected from aqueous solutions of the four organic modifiers. The isotherm model for solvent adsorption was selected using two independent parameter estimation methods, the regression analysis and adsorption energy distribution. The fitted isotherm parameters were tested by modeling of overloaded elution bands with the aid of the equilibrium-dispersive model of chromatography. The surface heterogeneity estimations and the effect of the silanol groups on the adsorption of solvents were based on those data.