In this paper the author will be describe phenomena of advection dispersion in subsurface flow by using Eulerian Lagrangian Finite Element Method. Where Finite Element Method with Galerkin formulation and weigthed residual method is used to solve seepage and advection dispersion equation. The problem of one dimensional and two dimensional rectangular wave are analyzed in this paper. And the result of numerical analyses will be compared with analytical solutions. The numerical results showed the very good agreement with the analytical solutions.