Abstract
In this manuscript we present a superconvergent discontinuous Galerkin method equipped with an element residual error estimator applied to scalar first-order hyperbolic problems using tetrahedral meshes. We present a local error analysis to derive a discontinuous Galerkin orthogonality condition for the leading term of the discretization error and establish new superconvergence points, lines and surfaces. We also derive new basis functions spanning the error and propose an implicit error estimation procedure by solving a local problem on each tetrahedron. The DG method combined with the a posteriori error estimation procedure yields both accurate error estimates and superconvergent solutions. Computations validate our theory.