The dynamic removal of copper by Purolite C100-MB cation exchange resin was studied in packed bed columns. The values of column parameters are predicted as a function of flow rate and bed height. Batch experiments were performed using the Na-form resin to determine equilibrium and kinetics of copper removal. The uptake of Cu(II) by this resin follows first-order kinetics. The effect of stirring speed and temperature on the removal kinetics was studied. The activation energy for the exchange reaction is 13.58 kJ mol(-1). The equilibrium data obtained in this study have been found to fit both the Langmuir and Freundlich isotherm equations. A series of column tests were performed to determine the breakthrough curves with varying bed heights and flow rates. To predict the breakthrough curves and to determine the characteristic parameters of the column useful for process design, four kinetic models; Bohart-Adams, Bed Depth Service Time (BDST). Clark and Wolborska models are applied to experimental data. All models are found suitable for describing the whole or a definite part of the dynamic behavior of the column with respect to flow rate and bed height. The simulation of the whole breakthrough curve is effective with the Bohart-Adams and the Clark models, but the Bohart-Adams model is better. The breakthrough is best predicted by the Wolborska model. The breakthrough data gave a good fit to the BDST model, resulting in a bed exchange capacity very close to the value determined in the batch process. (C) 2008 Elsevier B.V. All rights reserved.