Surface modelling has been a widely used methodology for interdisciplinary facilities in all kinds of earth-related studies. There are many interpolation methods applied for model generation using the measured points, samples, on the ground. The quality of the outcomes of an interpolation method is highly related to the accuracy, quantity, and distribution of the selected samples reflecting the topography within the study area. This study aims to examine the quality of four interpolation methods, namely the methods Kriging, Modified Shepard's, Inverse distance weighting, and Radial Basis Function, considering height differences between the neighbour stations. To check the quality of height components within the study area derived applying different interpolation models, four artificial surfaces with sudden height changes were created. The standard deviations used for comparison of the quality of interpolation models were determined using differences between the height values of control points. These values were set as true and interpolated values for the same points.