[Objective]The spatial variabilities and spatial distribution patterns of topsoil organic matter(SOM), total nitrogen(TN), available phosphorus(P) and available potassium(K) in Zengchen District of Guangzhou City were analyzed in order to provide some hints on how to improve fertilizer use efficiency and to prevent water and soil pollution in analogous area with lateritic red soil.[Method] Geostatistic and GIS methods are used in this study.[Result] The results showed that P content is high overall in the district. The maximum coefficient of spatial variation comes with K, up to 107.78%; whereas the minimum nutrient is SOM, about 32.51%. Coefficients of variation have a rank of K >P >TN >SOM. SOM and total N have distinct spatially isotropy auto-correlations; On the contrary, available P and K behave obvious anisotropy auto-correlations. The contents of SOM and TN are higher in the southeast and northwest, but lower in the central; P content of the district increases from north to south gradually. A general increase from north to south displays for the spatial distribution of K content, but it seems to be more complicated.[Conclusion] If observed within a certain space, topsoil SOM, TN, P and K have their different spatial correlations, upon which, a suitable directional semi-variance function model need to be find out to describe their spatial variabilities better.