Adsorption of levulinic acid onto a basic polymeric adsorbent was experimentally measured and numerical analyzed by batch and column methods. The adsorption equilibrium was correlated by a highly favorable Toth isotherm. Experimental adsorption rates in batch adsorber were fitted by the pore diffusion model and parallel pore/surface diffusion model. The resultant isotherm equation and diffusion coefficients were incorporated into the general rate model to predict the experimental breakthrough data. The film resistance and the axial dispersion in the column were estimated by Wilson-Geankoplis equation and Chung-Wen equation, respectively. Simulation elucidated that the parallel diffusion model successfully described the breakthrough behaviors and gave a better prediction than the pore diffusion model.