Surface complexation modeling was performed using the Generalized Two-Layer Model for a series of low molecular weight organic acids. Sorption of these organic acids to goethite was investigated in a previous study to assess the influence of particular structural features on sorption. Here, the ability to describe the observed sorption behavior for compounds with similar structural features using surface complexation modeling was investigated. A set of surface reactions and equilibrium constants yielding optimal data fits was obtained for each organic acid over a range of total sorbate concentrations. Surface complexation modeling successfully described sorption of a number of the simple organic acids, but an additional hydrophobic component was needed to describe sorption behavior of some compounds with significant hydrophobic character. These compounds exhibited sorption behavior that was inconsistent with ligand exchange mechanisms since sorption did not decrease with increasing total sorbate concentration and/or exceeded surface site saturation. Hydrophobic interactions appeared to be most significant for the compound containing a 5-carbon aliphatic chain, Comparison of optimized equilibrium constants for similar surface species showed that model results were consistent with observed sorption behavior: equilibrium constants were highest for compounds having adjacent carboxylic groups, lower for compounds with adjacent phenolic groups, and lowest for compounds with phenolic groups in the ortho position relative to a carboxylic group. Surface complexation modeling was also performed to fit sorption data for Suwannee River fulvic acid. The data could be described well using reactions and constants similar to those for pyromellitic acid. This four-carboxyl group compound may be useful as a model for fulvic acid with respect to sorption. Other simple organic acids having multiple carboxylic and phenolic functional groups were identified as potential models for humic substances.