First-principles calculations were utilized to elucidate the complete defect equilibria of surfaces of proton conducting BaZrO3, encompassing charged species adsorbed to the surface, defects in the surface layer as well as in the subsurface space-charge region and bulk.