Simultaneous adsorption modeling of four ions was predicted with a strict net charge-neutral ion-exchange theory and its corresponding equilibrium and mass balance equations. An important key to the success of this approach was the proper collection of all the data, particularly the proton adsorption data, and the inclusion of variable concentrations of conjugate ions from the experimental pH adjustments. Using IExFit software, the ion-exchange model used here predicted the competitive retention of several ions on goethite by assuming that the co-adsorption or desorption of all ions occurred in the correct stoichiometries needed to maintain electroneutrality. This approach also revealed that the retention strength of Cl− ions on goethite increases in the presence of phthalate ions. That is, an anion-anion enhancement effect was observed. The retention of Cl− ions was much weaker than phthalate ions, and this also resulted in a higher sensitivity of the Cl− ions toward minor variations in the surface reactivity. The proposed model uses four goethite surface sites. The drop in retention of phthalate ions at low pH was fully described here as resulting from competitive Cl− reactions, which were introduced in increasing concentrations into the matrix as the conjugate base to the acid added to lower the pH.