Novel catalytic adsorbent (ruthenium on carbon) was employed for the treatment of pulp mill effluent in the presence of hydrogen peroxide. Mathematical model and optimization of the process regarding the most favorable COD (%), TOC (%) and color (%) removal rates was developed and performed with experimental design taking into account catalytic adsorption process kinetics. As the initial experimental design, 3(3-1) half-fractional factorial design (H-FFD) was accomplished at two levels to study the significance of the main effects, such as catalytic adsorbent (g l(-1)) and hydrogen peroxide (ppm) concentrations using the response surface methodology (RSM). Finally, a four factor-three coded level central composite design (CCD) with 28 runs was performed in order to fit a second-order polynomial model. Validation of the model was accomplished by different criteria including coefficient of determination and the corresponding analysis of variance. The achieved removal rates for TOC (up to 75%), COD (up to 73%) and color (up to 68%) were observed for the defined optimal conditions: 1g l(-1) of ruthenium on carbon, 7 ppm of hydrogen peroxide, pH = 4 and ambient temperature. The proposed method benefited significantly improved TOC, COD and color removal efficiency, regenerability and reusability of the catalytic adsorbent and unaltered initial pH of an effluent in comparison to traditional adsorption or oxidation processes.