Abstract: We present a calculation of next-to-leading order QCD corrections to QCD-induced ZZ production in association with two jets at hadron colliders. Both Z bosons decay leptonically with all off-shell effects, virtual photon contributions and spin-correlation effects fully taken into account. This process is an important background to weak boson scattering and to searches for signals of new physics beyond the Standard Model. As expected, the next-to-leading order corrections reduce significantly the scale uncertainty and show a non-trivial phase space dependence in kinematic distributions. Our code will be publicly available as part of the parton level Monte Carlo program VBFNLO.