Recent spectroscopic work has suggested that only one surface species of sulfate is dominant on hematite. Sulfate is therefore a very suitable anion to test and develop adsorption models for variable charge minerals. We have studied sulfate adsorption on goethite covering a large range of sulfate concentrations, surface coverages, pH values, and electrolyte concentrations. Four different techniques were used to cover the entire range of conditions. For characterization at low sulfate concentrations, below the detection limit of sulfate with ICP-AES, we used proton-sulfate titrations at constant pH. Adsorption isotherms were studied for the intermediate sulfate concentration range. Acid-base titrations in sodium sulfate and electromobility were used for high sulfate concentrations. All the data can be modeled with one adsorbed species if it is assumed that the charge of adsorbed sulfate is spatially distributed in the interface. The charge distribution of sulfate follows directly from modeling the proton-sulfate adsorption stoichiometry since this stoichiometry is independent of the intrinsic affinity constant of sulfate. The charge distribution can be related to the structure of the surface complex by use of the Pauling bond valence concept and is in accordance with the microscopic structure found by spectroscopy. The intrinsic affinity constant follows from the other measurements. Modeling of the proton-ion stoichiometry with the commonly used 2-pK models, where adsorbed ions are treated as point charges, is possible only if at least two surface species for sulfate are used.