A mathematical model for galvanostatic soil electroremediation of cadmium-spiked kaolinite was built considering: chemical and electrical propel-ties from both phases, assumption of local equilibrium conditions, electroneutrality front phase interactions, a dynamic surface charge, and surface adsorbed concentrations defined by Boltzmann's equation. Model prediction trends agree with the expected phenomena to occur during an electroremediation process, as far as they show: a right electroosmotic fluid velocity with a suction effect at the anode position, a pH drop at the cathode position, soluble cadmium transport from anode to cathode, as well as cadmium accumulation before reaching the cathode, and a drop in electric potential gradient. Also, model strength allows the prediction of chemical species transient concentration; and provides an insight on species distribution through the spatial domain, while electroremediation takes place. (c) 2005 American Institute of Chemical Engineers.