The aim of the study was to develop a methodology for determining homogeneous territorial zones for precision farming. In this study we took into account the national land use system which provides for the absence of private ownership of agricultural land. The algorithm for determining management-zones provides for: establishing zones of spatial heterogeneity; determining the presence of clusters and emissions; modelling the spatial distribution of soil quality indicators. It is recommended to use data from agrochemical soil studies which are conducted centrally every 4 years for each agricultural enterprise as input parameters. These data include: the humus content in the soil, the content of available phosphorus and potassium and soil pH. The data should be carefully examined using spatial statistics tools to provide a more accurate delineation of the management-zones boundaries. The developed technique makes it possible to determine fertile and marginal areas within each individual field and differentiate the use of fertilizers, taking into account the presence of intra-field heterogeneity. This will save from 2.5 to 21.8 kg P haE−1 and from 0.9 to 26.7 kg K haE−1 due to the redistribution of the fertilizer dose calculated for the planned yield, taking into account the identified site-specific management zones. The differentiated use of mineral fertilizers will increase the profitability of growing winter cereals by 2.2%, sugar beets by 1.3%, rapeseed by 1.1%, and malting barley by 0.8%.

The main objective of this study is to review and evaluate three common interpolation methods namely: Inverse Distance Weighting (IDW), Radial Basis Function (RBF) and Ordinary Kriging (OK), and generate maps of soil pH using these methods. The accuracy and efficiency of the generated maps have been examined as well as the most fitting technique for estimating spatial distribution of soil pH in the study area is identified. Studies were conducted within the limits of land use of RUP “Uchkhoz BGSHA” (Republic of Belarus, Mogilev region, Goretsky district). The total area of the surveyed territory is 3197.89 hectares. For the analysis data is used about pHKCl of soil solution obtained from materials of an agrochemical survey executed in 2014. Forecasting and visualization of the spatial distribution of pH sub(KCl) was carried out using the Geostatistical Analyst module of the ArcGIS software. The experimental anisotropic variograms were calculated to determine the possible spatial structure of soil pH. Based on cross-validation results, a polynomial function was identified as the best variogram model. The model created by the method of radial basis functions turned out to be the most suitable for forecasting purposes (the value of the root-mean-square error was 0.763). In terms of interpolation accuracy, the investigated deterministic and geostatistical methods are located in the next descending row: RBF greater than IDW greater than OK.