The study on spatial variability of soil properties performed through geostatistical techniques allow us to identify the spatial distribution of phenomena by means of a spatial model that considers degree of dependence among observed data, depending on distance and also the direction that separate them, if there is geometric anisotropy, in other words, a directional trend in spatial continuity. However, the main difficulty in decision making regarding the use of anisotropic spatial model focuses on its relevance to the parameters that express the geometric anisotropy in a spatial model exercise in relation to the estimation space. This study aims at identifying the degree of influence of geometric anisotropy on the accuracy of spatial estimation using simulated data sets with different sample sizes and soil chemical properties such as: Fe, potential acidity (H + Al), organic matter and Mn. Comparing the isotropic and anisotropic models, especially for smaller sample sizes (100 and 169) showed an increased sum of squares of differences between predictions anisotropy factor (Fα) equals 2. Furthermore, from Fα equals 2.5, over 50% of the simulations showed values of overall accuracy (OA) of less than 0.80 and values for the concordance index Kappa (K) and Tau (T) from 0.67 to 0.80, indicating differences between thematic maps. Similar conclusions were obtained for chemical properties of the soil, from Fα equals 2, showing that there are relevant differences regarding the inclusion or not of geometric anisotropy.