With rapid economic and social development, soil contamination arising from heavy metals has become a serious problem in many parts of China. We collected a total of 445 samples (0-20 cm) at the nodes of a 2 km×2 km grid in surface soils of Rizhao city, and analyzed sources and risk pattern of 10 heavy metals (As, Cd, Co, Cr, Cu, Hg, Mn, Ni, Pb and Zn). The combination of Multivariate statistics analysis and Geostatistical methods was applied to identify the sources and hazardous risk of heavy metals in soils. The result indicated that Cr, Ni, Co, Mn, Cu, and As were mainly controlled by parent materials and came from natural sources. Cd and Hg originated from anthropogenic sources. Pb and Zn, belonging to different groups in multivariate analysis, were associated with joint effect of parent materials and human inputs. Ordinary Kriging and Indicator Kriging suggested that single element and elements association from the same principal components had similar spatial distribution. Through comprehensive assessment on all elements, we also found the high risk areas were located in the populated urban areas and western study area, which could be attributed to the higher geological background in the western part and strong human interference in the eastern part.