The influence of Fermi level position and annealing ambient on the zinc vacancy V Zn generation and Al diffusion is studied in monocrystalline zinc oxide (ZnO). From secondary-ion mass spectrometry and positron annihilation spectroscopy results, a quadratic dependence between the concentrations of V Zn and Al is established, demonstrating the Fermi level dependence of the formation of the electrically compensating −2 charge state of V Zn in conductive n-type ZnO crystals. In contrast, thermal treatment in the zinc-rich ambient is shown to efficiently reduce the V Zn concentration and related complexes. Using a reaction-diffusion model, the diffusion characteristics of Al at different donor background concentrations are fully accounted for by mobile (Al Zn V Zn) − pairs. These pairs form via the migration and reaction of isolated V 2− Zn with the essentially immobile Al + Zn. We obtain a migration barrier for the (Al Zn V Zn) − pair of 2.4 ± 0.2 eV, in good agreement with theoretical predictions. In addition to strongly alter the shape of the Al diffusion profiles, increasing the donor background concentration also results in an enhanced effective Al diffusivity, attributed to a reduction in the V 2− Zn formation energy as the Fermi level position increases.