We consider the numerical simulation of contaminant transport in groundwater where the mathematical model includes a nonlinear adsorption term. The method we describe combines a higher order Godunov scheme for advection with a mixed finite element method for diffusion. The method is formulated in one space dimension, and numerical results for equilibrium and nonequilibrium adsorption are presented.