Soil contamination is a major issue that endangers the ecology in most countries. Total concentrations of As, Cd, Co, Cr, Cu, Mn, Ni, Pb, VFe, and Zn were determined by analyzing soil samples from 32 surface soil samples in southwest Saudi Arabia, including certain areas of Al-Baha. Kriging techniques were used to create maps of the distribution of metal. To assess the levels of soil contamination in the research area, principal component analysis (PCA), contamination factors (CF), and pollution load index were used. The results show the stable model gave the best fit to the As and Zn semivariograms. The circular model fits the Cd, Co, and Ni semivariograms the best, while the exponential model fits the Cr, V, and Fe semivariograms the best. For Ni and Pb, respectively, spherical and Gaussian models are fitted. The findings demonstrated two clusters containing different soil heavy metal concentrations. According to the data, there were two different pollution levels in the research region: 36.58% of it is strongly contaminated, while 63.41% of it has a moderate level of contamination (with average levels of these metals 5.28 ± 5.83, 0.81 ± 0.19, 18.65 ± 6.22, 45.15 ± 23.25, 60.55 ± 23.74, 972.30 ± 223.50, 33.45 ± 14.11, 10.05 ± 5.13, 84.15 ± 30.72, 97.40 ± 30.05, and 43,245.00 ± 8942.95 mg kg-1 for As, Cd, Co, Cr, Cu, Mn, Ni, Pb, V, Fe, and Zn, respectively). The research area's poor management practices are reflected in the current results, which raised the concentration of harmful elements in the soil's surface layers. Ultimately, the outcomes of pollution concentration and spatial distribution maps could aid in informing decision-makers when creating suitable heavy metal mitigation strategies.