We study the mobility and interaction under competing conditions observed for copper (Cu2+) and zinc (Zn2+) ions in the context of laboratory-scale experiments performed in natural soil columns. The experiments focus on the analysis of solute breakthrough curves (BTCs) obtained after injection of an aqueous solution containing similar concentrations of the two metal ions into a soil column fully saturated with double deionized water. Transport of the competing ions is tested for the same soil under aerobic and anaerobic conditions. Measurements show that the species with lower affinity for the soil, Zn2+, migrates occupying all available adsorption sites, and is then progressively replaced by the ion with higher affinity, Cu2+. The two ions are displaced in the system with different effective retardation. The slowest species replaces the sorbed ions, resulting in observed Zn2+ concentrations that display a non-monotonic behavior in time and which, for a certain period, are larger than the concentration supplied continuously at the inlet. In the absence of a complete geochemical characterization of the system, we show that the measured concentrations of both metals can be interpreted through simple models based on a set of coupled partial differential and algebraic equations, involving a small subset of aqueous and adsorbed species that are present in the system. Depending on the model considered, the relationship between aqueous and adsorbed ion concentrations is described at equilibrium by a Gaines–Thomas (GT) formulation, a competitive Sheindorf–Rebhun–Sheintuch (SRS) isotherm, or an Extended Langmuir (EL) isotherm, respectively. The GT formulation provides the best interpretation of the observed behavior among the models tested. We find that employing these simple models, which account only for the main governing reactive processes, allows reasonable estimation of the observed BTCs in experiments where only partial geochemical datasets are available.