A Eulerian-Lagrangian scheme is used to solve the two-dimensional advection-dispersion equation. Concentration and its partial differential operator are decomposed into advection and dispersion terms. Thus, advection is formally decoupled from dispersion and solved by continuous forward particle tracking. Dispersion is handled by implicit finite elements on a fixed Eulerian grid. Translation of steep gradients of concentration in advection-dominated flow regimes, is done without numerical distortion. Continuous spatial distribution of velocities are evaluated by using Galerkin's approach in conjunction with Darcy's law based on hydraulic input data from each element. The method was implemented on coarse FE grid with linear shape functions, demonstrating no over/under shooting and practically no numerical dispersion. Simulations, covering a wide range of Peclet numbers, yield high agreement with analytic and practical results.