“…We set up a one-dimensional (1D) advective–dispersive reactive transport model to quantify the mobility and retention of glyphosate in the Tufa–sand packed columns. Sorption in the model was described via a combination of kinetic reversible sorption, accounting for Freundlich-type nonlinearity in the adsorption term and an irreversible sorption term after the multireaction model (MRM) for glyphosate sorption of Padilla and Selim. , We considered Freundlich-type sorption based on results from preliminary equilibrium sorption batch experiments, detailed in the SI, section 2. Equation is the advection–dispersion equation for a sorbing compound, in which C [mg L –1 ] is the aqueous phase concentration of glyphosate and S [mg g –1 ] is the concentration of solid-phase glyphosate, the sum of reversibly ( S r ) and irreversibly ( S ir ) bound fractions, n e [−] is the flow-effective porosity, ρ s [g L –1 ] is the mass density of the solids, v [m s –1 ] is the flow velocity, D [m 2 s –1 ] is the dispersion coefficient ( D = v α + D m , where α [m] is the dispersivity and D m [m 2 s –1 ] is the molecular diffusion coefficient in water), and t and x are the temporal and spatial coordinates. The rate of change of solid-phase glyphosate is the sum of the temporal derivatives of S r and S ir , which are given by the following differential equations where k att [mg (1–n) L –(1–n) s –1 ], k det [s –1 ], and k irr [s –1 ] are the rate constants for reversible sorption and desorption, and irreversible sorption, respectively.…”