Geostatistics combined with geographical information system(GIS) techniques were applied to a treatment and analysis of the spatial heterogeneity of soil nutrients in topsoil(0—20cm) at Hailun County of Heilongjiang Province. Sixty-five soil samples were collected from the central farmland of Hailun County. The results indicate that the semivariogram of total C and total P are best described by the Gaussian model; total N and total K are best described by the exponential model. Analysis shows that the spatial correlation of total C and total N at this large scale is moderately dependent, with Nugget/Sill were between 25%～75%; therefore changes in total C and total N are affected by both structural and random factors. The spatial correlation of total K and total P is highly dependent, with Nugget/Sill were less than 25%; therefore, changes in total K and total P are mainly affected by structural factors. There is high spatial correlation between total C and P. Spatial variability results from the geostatistics method were compared to soil fertility survey in 1978 results from Hailun farmland. The analysis is useful to inform strategic planning of farm production, to guide fertilizer application and to control soil erosion and degradation of Blacksoil, in northeastern China.