intrinsic equilibrium constants of 17 representative brazilian oxisols were estimated from potentiometric titration measuring the adsorption of h + and oh − on amphoteric surfaces in suspensions of varying ionic strength. Equilibrium constants were fitted to two surface complexation models: diffuse layer and constant capacitance. The former was fitted by calculating total site concentration from curve fitting estimates and pH-extrapolation of the intrinsic equilibrium constants to the pznpc (hand calculation), considering one and two reactive sites, and by the FITEQL software. The latter was fitted only by FITEQL, with one reactive site. Soil chemical and physical properties were correlated to the intrinsic equilibrium constants. Both surface complexation models satisfactorily fit our experimental data, but for results at low ionic strength, optimization did not converge in fiteQl. data were incorporated in Visual minteQ and they provide a modeling system that can predict protonation-dissociation reactions in the soil surface under changing environmental conditions.