We report a molecular simulation study of aqueous solutions at montmorillonite clay surfaces. Unlike most previous studies, ours does not focus on the interlayer nanopores, but looks at both kinds of external surfaces of clay particles: basal surfaces along the clay layers, and lateral surfaces through which interlayer and larger interparticle pores are linked. We present results on structural, dynamic and thermodynamic properties and phenomena, including hydration complexes of ions, H bonding networks, modification of the water dynamics with respect to the bulk, and the role of water in the cation exchange between interlayer and interparticle pores.