“…The morphology of the illite crystal and the complexity of its surface suggest a number of surface reactive sites, including weakly acidic basal planes and amphoteric silanol and aluminol sites on the edge surface [16][17][18][19]. For the most part, however, because of the clay structural complexity and the limited ability of current surface complexation models (SCMs), and optimization routines such those included, for example, in the FITEQL Ò code [20], to predict the number of reactive sites, two binding sites have been assumed on clay mineral surfaces [16,17,21,22]. In order to better understand the nature of metal and organic matter complexation to clays, not only the number of available sites, but their most accurate apparent pK a values and concentrations must be first defined.…”