Acid-base equilibria of a humic acid were characterized by considering ionic group heterogeneity and electrostatic effects. The continuous pK distribution model, in which acidic functional groups having two centers of distribution were assumed in the humic molecule, was adapted to the experimental titration data by non-linear least square regression analysis.
Excellent curve-fitting was obtained, and two centers of pK distribution ( p1 and pz) and the total amounts of functional groups (C,, and CA,) were evaluated. It was suggested that there are strong and weak acidic groups in the humic molecule and that they continuously deprotonate according to their pK distribution. The /.L values were evaluated in media of various ionic strength (I = 0.004-l); p, (strong groups): 3.4-5.1, pL2 (weak groups): 6.9-7.9, and the ratio of C,, to C,, was ca. 5:2. The electrostatic model including the Debye-Hiickel theory was adapted to evaluate the intrinsic dissociation constants, pKi,int. From the relationship between the p values and I'/', the intrinsic p Ki values of two centers of distributions, pin,, could be evaluated; /_Q~: 3.03, /_L~,~,,~: 5.12. The intrinsic pK distribution and the species distribution curves were calculated under electrostatic-free conditions. Investigation of Fourier transform infrared spectra and the species distribution curves suggested that the two different types of carboxylic groups would contribute to the p K distribution of each functional group.