Here we present a thorough density functional study combining experiments on ZnO nanostructures aimed at the identification, by means of Infrared (IR) spectroscopies, of hydroxyl and hydride species formed on most stable low-index Miller surfaces of würtzite ZnO; namely, Zn-and O-terminated (0001) and (0001) polar surfaces, and nonpolar (1010) and (1120) surfaces. The Perdew-Burke-Ernzerhof functional was employed on the periodic slab calculations, and all possible H and OH adsorption modes were studied at medium and full coverage, while IR spectra were simulated for most favourable situations. This information was used to model the most likely surface arrangements when exposed to either H 2 or H 2 O. IR experiments on ZnO surfaces and nanoparticles are discussed based on the calculated adsorption energy values and simulated IR spectra. The study emphasizes the detailed assignment of OH moieties with the help of IR and their interpretation as fingerprints of surface morphology, allowing for a consistent interpretation of water adlayers stability and their corresponding vibrational fingerprints as a function of coverage, low-index Miller surface, and hydrogen source.