In this paper we consider an undisturbed soil column that has been
exposed to polluted water through both infiltration and precipitation.
By identifying the dominant reactions that are taking place in the
soil a reaction diffusion model for the ionic species present in the
soil is developed. Parameters in the model are tuned by fitting
time dependent concentration data for a single reaction but the
presence of other ions was found to significantly affect these
parameters. The full model does seem capable of modeling the
experimentally observed break-through curves but a much greater
understanding of the underlying chemistry is required to develop
a more refined model.