The surface properties of the diatomite were investigated using nitrogen adsorption/deadsorption isotherms, TG-DSC, FTIR, and XPS, and surface protonation-deprotonation behavior was determined by continuous acid-base potentiometric titration technique. The diatomite sample with porous honeycomb structure has a BET specific surface area of 10.21 m(2)/g and large numbers of surface hydroxyl functional groups (i.e. Si-OH, Fe-OH, and Al-OH). These surface hydroxyls can be protonated or deprotonated depending on the pH of the suspension. The experimental potentiometric data in two different ionic strength solutions (0.1 and 0.05 mol/L NaCl) were fitted using ProtoFit GUI V2.1 program by applying diffuse double layer model (DLM) with three amphoteric sites and minimizing the sum of squares between a dataset derivative function and a model derivative function. The optimized surface parameters (i.e. surface dissociation constants (log K-1, log K-2) and surface site concentrations (log C)) of the sample were obtained. Based on the optimized surface parameters, the surface species distribution was calculated using Program-free PHREEQC 3.1.2. Thus, this work reveals considerable new information about surface protonation-deprotonation processes and surface adsorptive behaviors of the diatomite, which helps us to effectively use the cheap and cheerful diatomite clay adsorbent. (C) 2014 Elsevier B.V. All rights reserved.