Taking the total contents of iodine in surface soils (0 –20 cm) of high-quality fruit producing area in Botou, Cangxian, Xianxian three counties of Cangzhou city as studied object, the spatial variations of iodine were analyzed by geostatistics and GIS techniques. The results showed that iodine in Cangxian and Xianxian were of normal distribution and fitted with exponential model, while iodine in Botou fitted with Gaussian model after logarithm conversion. Iodine in Cangxian and Xianxian showed higher spatial correlation, which suggested the greater influences from natural structural factors, such as climate, soil parent material, topography and soil type; iodine in Botou showed medium spatial correlation, which stated the greater influences from random factors, such as human activities. By adopting the most excellent Kriging method, the spatial distribution maps of iodine in the three countries were edited, the results showed that the spatial distributions of iodine in three countries were asymmetrical, it was mostly affected by hydrology, geology and topography. Iodine contents in both shores of the rivers and low-lands were higher than other areas, whereas iodine contents in some areas showed partial abnormality owing to the influence of human activities.