We develop implicit a posteriori error estimators for elliptic boundary value problems. Local problems are formulated for the error and the corresponding Neumann type boundary conditions are approximated using a new family of gradient averaging procedures. Convergence properties of the implicit error estimator are discussed independently of residual type error estimators, and this gives a freedom in the choice of boundary conditions. General assumptions are elaborated for the gradient averaging which define a family of implicit a posteriori error estimators. We will demonstrate the performance and the favor of the method through numerical experiments.