It is shown that cubic spline method is effective for restraining numerical dispersion in numerical simulation of hydrodynamicdispersion equations with radioactive decay.
The variability of the pressure head, velocity and hydrodynamicdispersion coefficient obtained from conditional simulation are smaller than those obtained from unconditional simulation.