Reliability of equilibrium state evaluation about settlement slopes in the context of natural and human-made hazards is a complex issue. The geological structure of the vicinity of the upland slope in the urban environment of Warsaw is characterised by a significant spatial diversification of the layers. This is especially due to the glacitectonics in the Mio-Pliocene clays, which are located shallowly under the sandy tills’ formations. With substantial variability in the clay roof surface, point recognition by drilling is often insufficient. The use of electrical resistivity imaging (ERI) in the quasi-3D variant provides accurate images of the real ground conditions, which is crucial in optimal geotechnical design. In forecasting the behaviour of the slope, it is necessary to quantify the impact of spatially differentiated systems of disturbed layers on changes in the safety factor (SF), which corresponds to the observed landslide activity of the Warsaw Slope. This study concerns numerous calculation model analyses of the optional clay position in the context of slope stability conditions. A wide range of soil properties variability was taken into account, resulting from both lithogenesis and subsequent processes disintegrating the original soil structure. Regarding the geological conditions of the slip surface, the use of classical computational methods and numerical modelling (FEM) was considered for comparative purposes. The results indicated that local changes in equilibrium conditions were affected by the different morphology of the clay roof surface of the slope and the alternation in strength characteristics on the slip surfaces. The findings of the study contribute to sustainable spatial planning of near-slope regions.