Simulation models for the fate and transport of groundwater contaminants are important tools for testing our understanding of transport phenomena at long-term contaminated sites and for designing remedial action plans. A finite difference formulation for contaminant transport including a distribution of contaminant mass-transfer rates between the water and soil is developed. Optimal model simulations using both log-normal and gamma distributions of mass transfer rates are compared to the two-site equilibrium/kinetic model. In all cases, optimal sorption parameters were determined by best fit to laboratory data. For desorption of trichloroethene from long-term contaminated soils, the distributed mass-transfer rate model provided significantly improved simulations of aqueous concentrations, as compared to the two-site model, for both batch and soil column experiments. However, use of an apparent partition coefficient demonstrated that the performance of the two-site model was very sensitive to the value of the partition coefficient, while the performances of the distributed models were robust over a wide range of partition coefficients. Desorption studies in continuous-flow stirred tank reactors with laboratory-contaminated soils demonstrated that as the length of the contamination period increases, the simulation capability of the two-site model decreases.