Let f be an automorphism of a complex Enriques surface Y and let p f denote the characteristic polynomial of the isometry f * of the numerical Néron-Severi lattice of Y induced by f . We apply a modification of McMullen's method to prove that the modulo-2 reduction (p f (x) mod 2) is a product of modulo-2 reductions of (some of) the five cyclotomic polynomials Φm, where m ≤ 9 and m is odd. We study Enriques surfaces that realize modulo-2 reductions of Φ7, Φ9 and show that each of the five polynomials (Φm(x) mod 2) is a factor of the modulo-2 reduction (p f (x) mod 2) for a complex Enriques surface.