The GNSS (Global Navigation Satellite System) coordinates time series are still used as a source for determining the velocities of GNSS permanent stations. These coordinates, apart from the geodynamical signals, also contain an interference signal. This paper shows the results of the comparative analysis of the GNSS coordinates time series with a deformation of the Earth's crust obtained from loading models. In the analysis, coordinates time series are used (CODE Repro2013) without loading models (Atmospheric Pressure Loading, Hydrology, Non-Tidal Ocean Loading) at the stage of the reprocessing of GNSS archival data. The analyses showed that in the case of the Up component there is a high correlation between the GNSS coordinates changes and deformations of the Earth's crust from the loading models (coefficient 0.5-0.8). Additionally, we noticed that for horizontal components (North, East) changes occur in the phase shift between coordinates, and the Earth’s crust deformations signals are accelerated or delayed each other (-150 to 200 days). This article shows new methods of iLSE (iteration Least Square Estimation) to determine periodic signals in the time series. Additionally, we compared the values of estimated amplitudes for GNSS and deformation time series. and Kaczmarek Adrian, Kontny Bernard.