<p>Soil contamination is a problem that negatively impacts the environment and, consequently, human health worldwide. To control and/or mitigate the environmental hazard generated by soil contamination, it is important to understand the spatial behavior of the contaminants and to model the level of contamination. In the present work, the geostatistical simulation technique is applied to model different pollutants present in the soils of the autonomous community of Cantabria, Spain. The results from the geostatistical simulation allow the determination of the probability that pollutant elements exceed allowable thresholds. The findings demonstrated that As concentrations exceed the permissible limits established for soils designated for mining/industrial and urban uses, as defined by Generic Reference Levels. These results highlight the need not only for preventive and corrective measures in the affected areas but also for strategies to mitigate the adverse effects of contamination. The results also suggest that geostatistical simulation algorithms can be effectively applied to support environmental management.</p>