Acid–base titration data for two goethites samples in sodium nitrate and sodium chloride media are discussed. The data are modeled based on various surface complexation models in the framework of the multi site complexation (MUSIC) model. Various assumptions with respect to the goethite morphology are considered in determining the site density of the surface functional groups. The results from the various model applications are not statistically significant in terms of goodness of fit. More importantly, various published assumptions with respect to the goethite morphology (i.e., the contributions of different crystal planes and their repercussions on the “overall” site densities of the various surface functional groups) do not significantly affect the final model parameters within simple 1-pK approximations. The simultaneous fit of the chloride and nitrate data results in electrolyte binding constants, which are applicable over a wide range of electrolyte concentrations including mixtures of chloride and nitrate. Model parameters for the goethite sample with 90 m2/g specific surface area are in excellent agreement with parameters that were independently obtained by another group on different goethite titration data sets.