Summary: The issue of a posteriori error estimation is discussed. In particular, we derive a theorem showing the dependence of the effectivity index for the Zienkiewicz-Zhu error estimator on the convergence rate of the recovered solution. This shows that with superconvergent recovery the effectivity index tends asymptotically to unity. The superconvergent recovery technique developed in the first part of the paper [see the foregoing entry] is used in the computation of the Zienkiewicz-Zhu error estimator to demonstrate accurate estimation of the exact error attainable. Numerical tests are shown for various element types illustrating the excellent effectivity of the error estimator in the energy norm and pointwise gradient (stress) error estimation. Several examples of the performance of the error estimator in adaptive mesh refinement are also presented.