This paper presents a locally conservative finite element method based on enriching the approximation space of the continuous Galerkin method with elementwise constant functions. The proposed method has a smaller number of degrees of freedom than the discontinuous Galerkin method. Numerical examples on coupled flow and transport in porous media are provided to illustrate the advantages of this method. We also present a theoretical analysis of the method and establish optimal convergence of numerical solutions.