Abstract: To analyze the spatial variability of soil nutrients in Huaxian, based on GIS and geostatistics, descriptive statistics analysis, semivariogram analysis and Kriging analysis were carried out on 924 soil samples in 2012. The results showed that the average contents of soil organic matter, total nitrogen, available phosphorus and available potassium in Huaxian were 14.81 g/kg, 0.87 g/kg, 16.02 mg/kg and 96.72 mg/kg, respectively, and the coefficient of variation was between 21.8% and 57%, which belonged to intermediate variability. The nugget/sill ratios were 50%, 54.6%, 89.2% and 90.1%, respectively. Soil organic matter and total nitrogen had a moderate spatial correlation, and available phosphorus and available potassium had a weak spatial correlation. The optimal interpolation model of four soil nutrients was the exponential model. The optimal search directions of soil organic matter and available phosphorus were 8 sectors, and the total nitrogen and available potassium were 4 sectors. The contents of organic matter and total nitrogen in Huaxian were below the moderate level, and the contents of available phosphorus and available potassium were at a moderate level. All the soil nutrients had low content regions respectively. Reasonable regulation was required in soil management and fertilization.