This paper introduces a case study of 6 townships in Lingyuan City,Liaoning Province.1042 surface soil samples were collect ed with the aid of GPS positioning.Traditional statistics combined with geostatistics was applied to analyze the spat ial variation of soil nutrients.Results show ed that the variation coefficients of four kinds of soil nutrients ranged from 29.45%to 67.25%.Soil nutrients,interms of the spatial heterogeneity,was in the descendant order of AP,AN,OM,and AK.AN,A K,and OM had m oderate spatial dependence.However,the C0/(C0+ C)for AP was 100%,so the relevance of space was weak.Soil nutrients,in terms of the fractal dimension,was in the descendant order of AP,AN,AK,and OM.Analysis by Krig ing method indicated that spatial pattern of soil nutrients in the study area was mainly affected by landuse type.