The spatial variability of organic matter (OM), available nitrogen (AN), available phosphorus (AP), and available potassium (AK) in 0-20 cm horizon in cultivated soils of Sujiatun District, Shenyang City were investigated using geostatistics and geographic information system (GIS). Data of each variable fitted normal distribution or normal distribution after log transformation, and a spherical model for semivariogram fitted all the variables. Semivariograms indicated that OM content strongly spatially dependent with the ratio of nugget to sill of 24.11%. The contents of AN, AP and AK were moderately spatially dependent with the ratio of nugget to sill of 29.53%, 60.77% and 58.82%, respectively. The ranges of spatially dependent for OM, AN, AP and AK were 26.051, 28.293, 15.132 and 23.813 km, respectively. The spatial variability for OM and AN was fundamentally affected by structural factors such as soil formation factors, and that for AP and AK was by random factors such as fertilization. Kriging maps showed the spatial distributions of the four variables, which were useful in monitoring nutrient changes and providing quantitative support to decision and policy making for agriculture and environment management.