Fuzzy classification combined with spatial prediction was used to assess the state of soil pollution in the peri-urban Beijing area. Total concentrations of As, Cr, Cd, Hg, and Pb were determined in 220 topsoil samples (0-20 cm) collected using a grid design in a study area of 2 600 km2. Heavy metal concentrations were grouped into three classes according to the optimum number of classes and fuzziness exponent using the fuzzy c-mean (FCM) algorithm. Membership values were interpolated using ordinary kriging. The polluted soils of the study area induced by the measured heavy metals were concentrated in the northwest corner and eastern part, especially the southeastern part close to the urban zone, whereas the soils free of pollution were mainly distributed in the southwestern part. The soils with potential risk of heavy metal pollution were located in isolated spots mainly in the northern part and southeastern corner of the study region. The FCM algorithm combined with geostatistical techniques, as compared to conventional single geostatistical kriging methods,could produce a prediction with a quantitative uncertainty evaluation and higher reliability. Successful prediction of soil pollution achieved with FCM algorithm in this study indicated that fuzzy set theory had great potential for use in other areas of soil science.