Competitive interaction of carbonate and phosphate on goethite has been studied quantitatively. Both anions are omnipresent in soils, sediments, and other natural systems. The PO4-CO3 interaction has been studied in binary goethite systems containing 0-0.5 M (bi)carbonate, showing the change in the phosphate concentration as a function of pH, goethite concentration, and carbonate loading. In addition, single ion systems have been used to study carbonate adsorption as a function of pH and initial (H)CO3 concentration. The experimental data have been described with the charge distribution (CD) model. The charge distributions of the inner-sphere surface complexes of phosphate and carbonate have been calculated separately using the equilibrium geometries of the surface complexes, which have been optimized with molecular orbital calculations applying density functional theory (MO/DFT). In the CD modeling, we rely for phosphate on recent parameters from the literature. For carbonate, the surface speciation and affinity constants have been found by modeling the competitive effect of CO3 on the phosphate concentration in CO3-PO4 systems. The CO3 constants obtained can also predict the carbonate adsorption in the absence of phosphate very well. A combination of inner- and outer-sphere CO3 complexation is found. The carbonate adsorption is dominated by a bidentate inner-sphere complex, (FeO)2CO. This binuclear bidentate complex can be present in two different geometries that may have a different IR behavior. At a high PO(4) and CO3 loading and a high Na+ concentration, the inner-sphere carbonate complex interacts with a Na+ ion, probably in an outer-sphere fashion. The Na+ binding constant obtained is representative of Na-carbonate complexation in solution. Outer-sphere complex formation is found to be unimportant. The binding constant is comparable with the outer-sphere complexation constants of, e.g., SO(2-)4 and SeO(2-)4.