Organoclays are sorbent materials prepared from clays by exchanging inorganic with organic cations. Their properties depend on the loading and conformational structure of the organic cations, but little information is available about the surface structures of organoclays. In this work, X-ray photoelectron spectroscopy (XPS) and classical molecular dynamics (MD) simulations are combined to characterize the external interface of an organoclay prepared from hexadecylpyridinium (HDPy+) and bentonite. The XPS survey spectra show well the varying elemental composition of the surface with increasing amount of surfactant, showing a decreasing contribution of clay-derived elements with increasing organic coverage. The high-resolution C 1s XPS spectra depict sensitively the surface arrangement of the surfactant. In combination with MD simulations, the results implied a monolayer coating for low surfactant coverage and a disordered bilayer arrangement at high surfactant uptakes. Molecular dynamics simulations showed that for very high cation uptake a quasi-paraffin-like configuration is also possible. The combination of experimental and modelling methods yielded congruent information on the molecular-scale arrangement of organic cations at the organoclay surfaces and the controlling mechanisms.