Resolving the atomic structure of the surface of ice particles within clouds, over the temperature range encountered in the atmosphere and relevant to understanding heterogeneous catalysis on ice, remains an experimental challenge. By using first-principles calculations, we show that the surface of crystalline ice exhibits a remarkable variance in vacancy formation energies, akin to an amorphous material. We find vacancy formation energies as low as similar to 0.1-0.2 eV, which leads to a higher than expected vacancy concentration. Because a vacancy's reactivity correlates with its formation energy, ice particles may be more reactive than previously thought. We also show that vacancies significantly reduce the formation energy of neighbouring vacancies, thus facilitating pitting and contributing to pre-melting and quasi-liquid layer formation. These surface properties arise from proton disorder and the relaxation of geometric constraints, which suggests that other frustrated materials may possess unusual surface characteristics. 2 Despite ice being a ubiquitous and well-studied substance, it is surprising that some basic questions about its properties and structure are still debated. For example, Faraday contentiously proposed that the surface of hexagonal ice (Ih) was liquid-like to explain a