Major and trace elements in soils originate from natural processes and different anthropogenic activities which are difficult to discriminate. On a 17-ha impacted site in northern France, two industrial sources of soil contamination were xidentified: a former iron foundry and a current secondary lead smelter. To discriminate and map natural and anthropogenic sources of major and trace elements on this site, the rarely applied MULTISPATI-principal component analysis (PCA) method was used. Using a 20-m x 20-m grid, 247 topsoil horizons were sampled and analysed with a field-portable X-ray fluorescence analyser for screening soil contamination. The study site was heavily contaminated with Pb and, to a lesser degree, with Sn. Summary statistics and enrichment factors allowed the differentiation of the main lithogenic or anthropogenic origin of the elements. The MULTISPATI-PCA method, which explained 73.9% of the variability with the three first factors, evidenced strong spatial structures. Those spatial structures were attributed to different natural and artificial processes in the study area. The first axis can be interpreted as a lithogenic effect. Axes 2 and 3 reflect the two different contamination sources. Pb, Sn and S originated from the secondary lead smelter while Fe and Ca were mainly derived from the old iron foundry activity and the old railway built with foundry sand. This study demonstrated that the MULTISPATI-PCA method can be successfully used to investigate multicontaminated sites to discriminate the various sources of contamination.