Linear alkylbenzene sulfonate (LAS) anionic surfactants are widely used, present in wastewater, and can be discharged causing environmental damage. Without considering biodegradation, sorption and desorption reactions play an important role and depend on the characteristics of the aquifer medium such as the content of organic matter and clay. These processes are modelled with PHREEQC using laboratory-obtained data on columns with an input LAS pulse. LAS adsorption and desorption step experiments carried out in sand columns were also modelled with this software and distribution constants were similar to those obtained for pulse. The values of the starting distribution constant provided by batch tests with sand and agricultural soil mixtures establish the maximum value that must be obtained for column modelling. Distribution constants from column experiments obtained by Freundlich and linear isotherms provided very similar results, the latter with more consistent achievement. The modelling results showed longer elution times in the case of the heavier homologues and with a higher percentage of agricultural soil, which would translate into the great persistence of these compounds in the environment. In the case of soil contamination with wastewater-containing surfactants, these results would allow the development of mitigation strategies for these environmental impacts.