Numerical solution of the advectionedispersion equation, used to evaluate transport of solutes in porous media, requires discretization schemes for space and time stepping. We examine use of quadratic upstream interpolation schemes QUICK, QUICKEST, and the total variation diminution scheme ULTIMATE, and compare these with UPSTREAM and CENTRAL schemes in the HYDRUS-1D model. Results for purely c...