The aim of the research is to diagnose the metrological characteristics of high-precision GNSS-observations by methods of non-classical error theory of measurements (NETM) based on Ukrainian reference stations. Methodology. We selected 72 GNSS reference stations, downloaded daily observation files from the LPI analysis center server, and created time series in the topocentric coordinate system. The duration of the time series is almost two years (March 24, 2019 - January 2, 2021). Using a specialized software package, the time series have been cleaned of offsets and breaks, seasonal effects, and the trend component has been removed. Verification of empirical distributions of errors was provided by the procedure of NETM on the recommendations offered by G. Jeffries and on the principles of hypothesis tests the theory according to Pearson's criterion. The main result of the research. It is established that the obtained time series of coordinates of reference GNSS stations do not confirm the hypothesis of their conformity to the normal Gaussian distribution law. NETM diagnostics of the accuracy of high-precision GNSS measurements, which is based on the use of confidence intervals for assessing the asymmetry and kurtosis of a significant sample, followed by the Pearson test, confirms the presence of weak, not removed from GNSS-processing, sources of systematic errors. Scientific novelty. The authors use the possibility of NETM to improve the processing of high-precision GNSS measurements and the need to take into account the sources of systematic errors. Failure to take into account certain factors creates the effect of shifting the time coordinate series, which, in turn, leads to subjective estimates of station velocity, i.e. their geodynamic interpretation. Practical significance. Research of the reasons for deviations of errors distribution from the established norms provides metrological literacy of carrying out high-precision GNSS measurements of large samples.