The concentrations of nine heavy metals (Fe, Mn, Zn, Cu, Ni, Cd, Cr, Co, and Pb) in soil samples of Arabian Gulf coast, Saudi Arabia, were investigated. Sediment Quality Guideline (SQGs), SQG-Quotient (SQG-Q), toxicity degree index (TDI), enrichment factor (EF), and multivariate analysis, including principle component analysis (PCA) and hierarchical cluster analysis (HCA), were used to measure heavy metals of concern in the study area and to identify their possible sources. The results showed that the concentrations of different heavy metals were 530–5,700 mg kg-1 for Fe; 9-150 mg kg−1 for Mn and 8-69 mg kg−1 for Zn; 1-21 mg kg−1 for Cu and not detectable-17 mg kg−1 for Ni; 6.9-130 mg kg−1 for Cr and not detectable-5 mg kg−1 for Co; and not detectable-24 mg kg−1 for Pb. Based on SQGs, only the maximum and mean Cr concentrations of study area were in heavy and moderate rates of pollution, respectively. The values of SQG-Q and TDI revealed that the investigated sampling points showed the lowest potential of adverse biological effects. The considerable number of collected soil samples has relatively higher EF values of 5–20 for Zn, Cu, Cr, and Pb, suggesting that these four metals may be derived from anthropogenic origin. Multivariate analysis also confirmed this finding that the sources of Zn, Cu, and Pb resulted primarily from anthropogenic sources, whereas Co, Ni, Fe, and Mn were mainly attributed to lithogenic sources. It could be generally concluded that it is possible to use multivariate analyses in combination with EF values as useful tools to identify the natural or anthropogenic sources of heavy metals in soils.