The first aim is to develop a model equation that describes the spread of heavy metals in soil vertically, for any depth and time then solving this model equation by using ADM. The second aim is suggesting a numerical model which describes the spread of heavy metals in soil horizontally. The suggested model can be used to estimates the concentration of heavy metals in soil for any time. So, this model would be a method to determine the contamination levels of some heavy metals in soil such Copper, Lead, Zinc, Cadmium, Cobalt, and Nickel with the potential for this contamination sources with their impact. The third aim is to develop a conceptual theory of interpolation in two dimensions then illustrating the accuracy, efficiency, rapidly and easy implementation of suggested methods. The fourth aim is to compare the performance of the suggested model with traditional methods, by estimating the concentration of heavy metals in soil.