We implemented a numerical model to simulate transport of multiple species and geochemical reactions occurring during electrokinetic remediation of metal-contaminated porous media. The main phenomena described by the model were: (1) species transport by diffusion, electromigration and electroosmosis, (2) pH-dependent buffering of H+, (3) adsorption of metals onto particle surfaces, (4) aqueous speciation, (5) formation and dissolution of solid precipitates. The model was applied to simulate the electrokinetic extraction of heavy metals (Pb, Zn and Ni) from marine harbour sediments, characterized by a heterogeneous solid matrix, high buffering capacity and aged pollution. A good agreement was found between simulations of pH, electroosmotic flow and experimental results. The predicted residual metal concentrations in the sediment were also close to experimental profiles for all of the investigated metals. Some removal overestimation was observed in the regions close to the anode, possibly due to the significant metal content bound to residual fraction. (C) 2016 Elsevier B.V. All rights reserved.