In the present study, modeling of 19 adsorption isotherms, 8 kinetic models and thermodynamics of methylene blue (MB) adsorption on chitosan flakes synthesized using Archachatina marginata shell wastes was investigated in a batch mode. The operational parameters' effects on the MB adsorption were studied. The model parameters were statistically analyzed using 10 error functions. The choices of the best fitted adsorption and kinetic models were based on the comparison of the sum of normalized error (SNE) and two statistical tools of information-based criteria. The 5-p Fritz-Schlüender isotherm best fitted the experimental adsorption data of MB on chitosan flakes based on SNE whereby maximum adsorption capacity, ax m q , of 143.6660 mg/g was obtained. The adsorption rate of MB on chitosan flakes was kinetically described by pseudo second-order model at all initial concentrations of MB investigated, with film diffusion being the rate-controlling step and the adsorption process chemisorption-influenced. The calculated thermodynamic parameters, 0 H = 4.23 kJ/mol, 0 S = 0.4563 kJ/(mol K), negative 0 G values revealed that the adsorption of MB onto chitosan flakes was physical, endothermic, spontaneous, energetically favorable and exergonic. The reaction mechanism of the adsorption of MB onto chitosan flakes was proposed taking cognizance of the electrostatic force of attraction between the negatively charged surface of the chitosan (biosorbent) and the positively charged MB.