The article deals with assessment of the Earth`s gravity field variations (GFV) affecting the repeated leveling results using modern geospatial analysis tools in the study of natural and engineering objects’ deformations. The basis of the calculation is put in the well-known formula that connects reference level of the point displacement with potential of gravity variation on it, caused by the volume of the redistributed surrounding masses. The desired change in the potential is determined by the method of numerical integration through the sum of elementary masses redistributed volume’s contributions to it. At estimating we used the possibility of obtaining the most reliable digital model of various landforms in the ArcGIS software environment and determining the dimensions of a complex body with a fairly high accuracy. In this case, the advantage of the rectangular partition into elementary volumes of the quasi-cone layer base, which is close in geometry to the shapes of igneous or man-made rocks model over the radial one is shown. The developed method for calculating the impact of GFV on leveling using the capabilities of ArcGIS is fully automated through compiling the corresponding script in python. Specific experimental material shows that it is possible to assess the above mentioned impact with any detailing: from the reading on the rod and the elevation at the station to the assessment of the corrections of the marker frame relative to the reference start outside the influence zone of GFV redistributed from the masses. Hence, estimating the GFV affecting the results of high-precision leveling should be coordinated with the method of geometric leveling, since the real corrections for gravitational changes in the reference mark may be less than the displacement of the level surface caused by the GFV on the estimated reference point.