Soil represents a complex structure of specific chemical and physical properties. Spatial variability in the soil depends on many factors: soil age, climate, parent material, topography and human activities. These factors, individually or in combination, affect soil parameters differently which affects his heterogeneity on micro, regional and national scale. Therefore, it is necessary to make decisions that will take into account the variability and the only appropriate way is precise mapping of its properties. In this thesis, the investigated area includes 1000 samples which were collected from the area Orahovica in a period from 2010.-2012. All the data were used to describe spatial variability based on four parameters: pH, organic matter, plant-available phosphorus and plant-available potassium. The aim of this investigation is to determine spatial variability as well as potential correlations between the investigated parameters in order to estimate the quality of geostatistical visualisation of investigated parameters. In order to accomplish this aim, descriptive statistic and geostatistical analytics will be applied. The results of this investigation will provide a better overview on the spatial variability of chemical properties on the investigated area of Orahovica. Spatial correlation of basic soil parameters will be described in order to create interpolation maps with significant accuracy, as well to reduce required number of samples to describe and make good decisions in soil management of investigated area.