The ultraconvergence property of a gradient recovery technique proposed by Zienkiewicz and Zhu is analyzed for the Laplace equation in the two dimensional setting. Under the assumption that the pollution effect is not present or is properly controlled, it is shown that the convergence rate of the recovered gradient at an interior node is two orders higher than the optimal global convergence rate when even-order finite element spaces and local uniform rectangular meshes are used.