The latest implementation of the international global coordinate system ITRF2020 enables taking into account not only linear time parameters of points’ coordinates changes participating in its construction, but also periodic terms reflecting their seasonal (intra-annual) variation, as well as post-seismic effects in tectonically active areas. These shifts for control points are not only caused by the actual physical movement of the point itself relative to neighboring ones or together with them on the tectonic plate. It also includes the influence of joint processing: the mentioned changes may contain common periodic terms distributed even among points that do not have real intra-annual variations. In this study, a harmonic analysis of the velocity field’s horizontal components in the ITRF2020 coordinate system was performed, the expansion was made in terms of spherical functions from the degree of 0 to that of 360; smoothed global features of the velocity field are obtained by summing only up to the level of 5 and 10. Although the use of expansion in terms of spherical functions is a methodologically more correct tool for interpolating fields on the surface of bodies close to a sphere, the rates of change by their nature are a discontinuous function and their harmonic analysis is difficult