The first-principles density functional theory is used to study the interaction of atomic hydrogen with ZnO surfaces. We find that atomic-hydrogen environments significantly reduce ZnO surface formation energies with hydrogen adsorption on the surface. The negative surface energy of the O-terminated ZnO(0001) surface is demonstrated in the O- and H-rich limits. The roughening and damage of ZnO surfaces are discussed in the context of the fluctuation of the surface formation energy.