Degradation of soil due to unsuitable land management practices is a chief impairment of optimum land productivity. The spatial variability of soil properties is needed for agricultural productivity, food safety and environmental modeling. The present study was conducted in lateritic soils of West Bengal, India to understand the spatial variability of soil properties using a geostatistical model. Nitrogen (N), soil pH, electrical conductivity (EC), Phosphorus (P), Potassium (K) and organic carbon (OC) were measured. Surface maps of soil properties were prepared using the semivariogram model through Kriging techniques. A positive correlation was observed between OC and N. The Quantile-quantile plots showed a normal distribution of EC, K, pH, N, and OC. The value for nugget/sill of K, N, and EC were 0.25–0.75 indicating moderate spatial autocorrelation among the variables. Phosphorus (P) was highly concentrated in the eastern part, whereas the agglomeration of higher EC was found in the north east and south west corner of the study site. The cross validation results illustrated the smoothing effect of the spatial prediction. The present study suggests that the geostatistical model can directly reveal the spatial variability of lateritic soils and will help farmers and decision makers for improving soil-water management. Keywords: Soil properties, Kriging techniques, Semivariogram, Lateritic soil, India