The article presents the experience of using geostatistical methods for reconstruction of agro-zones near the borders of ancient Kerkinitis (VI century BC –II century BC) according to the geochemical soil analysis. In GIS application ArcGIS 10.1 cluster analysis of soil objects (n = 122) was made for six geochemical indicators and organic carbon content and integral indicator of soil properties was calculated. Received two spatially clustered arrays of points, which differ significantly on set of soil properties that reflect the duration of the agricultural load. As a result of geostatistical analysis a pronounced spatial data autocorrelation trends and field uniformity were detected. By probability kriging border of the transition zone between the different age arable land was substantiated, as well as «anomalies» areas were installed, that may be of interest for future archaeological research.