2019 marked the 20th anniversary of the International VLBI Service for Geodesy and Astrometry (IVS). This service is the largest and most authoritative organization that coordinates international activities in radio astrometry and VLBI sub-system of space geodesy. Currently, about 60 antennas located in many countries on all continents participate in the IVS observing programs. The IVS Data Centers have accumulated more than 18 million observations obtained during more than 17'000 sessions, including more than 10,000 Intensive sessions for rapid determination of Universal Time. The paper traces the dynamics of IVS development based on statistical processing of the array of observations collected in the IVS Data Centers for the period of 1979-2018. Various statistics by the years, stations, baselines, and radio sources are provided. The evolution of the IVS observational data and the accuracy of results obtained from processing VLBI observations is considered.

Lidia Dmitrievna Kostina and Natalia Romanovna Persiyaninova left a bright mark in the history of the Pulkovo Observatory, as well as in the history of the domestic and international latitude services. In the first place, they were absolute leaders in the latitude observations with the famous zenith telescope ZTF-135. In 1954-2001, they obtained together more than 66'000 highly accurate latitudes, which make about 2/3 of all the observations made by 23 observers with the ZTF-135 after the WW2. They also provided a large contribution to investigation of the instrumental errors, methods of the data analysis, developing of the observing programs. Their results in studies of the latitude variations and polar motion were also highly recognized by the community.

A new method Spherical Rectangular Equal-Area Grid (SREAG) was proposed in Malkin (2019) for splitting spherical surface into equal-area rectangular cells. In this work, some more detailed features of SREAG are presented. The maximum number of rings that can be achieved with SREAG for coding with 32-bit integer is $N_{ring}$=41068, which corresponds to the finest resolution of $\sim$16$''$. Computational precision of the SREAG is tested. The worst level of precision is $7\cdot10^{-12}$ for large $N_{ring}$. Simple expressions were derived to calculate the number of rings for the desired number of cells and for the required resolution.

A new method SREAG (spherical rectangular equal-area grid) is proposed to divide a spherical surface into equal-area cells. The method is based on dividing a sphere into latitudinal rings of near-constant width with further splitting each ring into equal-area cells. It is simple in construction and use, and provides more uniform width of the latitudinal rings than other methods of equal-area pixelization of a spherical surface. The new method provides a rectangular grid cells with the latitude- and longitude-oriented boundaries, near-square cells in the equatorial rings, and the closest to uniform width of the latitudinal rings as compared with other equal-area isolatitudinal grids. The binned data is easy to visualize and interpret in terms of the longitude-latitude rectangular coordinate system, natural for astronomy and geodesy. Grids with arbitrary number of rings and, consequently, wide and theoretically unlimited range of cell size can be built by the proposed method. Comparison with other methods used in astronomical research showed the advantages of the new approach in sense of uniformity of the ring width, a wider range of grid resolution, and simplicity of use.

A new version of the Optical Characteristics of Astrometric Radio Sources (OCARS) catalog is presented. This compiled catalog includes radio sources observed in different VLBI programs and experiments that result in source position determination, their redshift, and photometry in the visible and near-infrared bands. A cross-identification table between the OCARS and other catalogs is also provided. The status of the catalog as of 2018 September 7 is described in this paper. The OCARS catalog currently contains 6432 sources, of which 3895 have redshift data and 5479 have photometric data. Compared with the previous version, the current version has been enriched with extended redshift and photometry information, and cross-identification with several catalogs in radio, optical, infrared, ultraviolet, X-ray, and gamma-ray bands. The OCARS catalog is updated every few weeks on average to incorporate new data that appear in the NASA/IPAC Extragalactic Database (NED), SIMBAD database, and in the literature.

The link problem between radio (VLBI/ICRF) and optical (Gaia/GCRF) celestial reference frames is analyzed. Both systems should be a realization of the ICRS (International Celestial Reference System) at microarcsecond level of accuracy. Therefore, the link between the ICRF and GCRF should be obtained with similar accuracy, which is not a trivial task due to relatively large systematic and random errors in source positions at different frequency bands. In this presentation, additional possibilities to improve the GCRF-ICRF link accuracy are discussed. In particular, a possibility to increase the number of ICRF and GCRF common objects is considered using advanced scheduling of the regular IVS sessions such as R1 and R4. It is shown that inclusion of supplement prospective southern sources in these sessions allows enriching southern ICRF zone without noticeable loss of accuracy of geodetic results. Another topic discussed in this presentation is using the correlations between radio source coordinates, which can impact the orientation angles between two frames at a level of a few tens of microarcseconds.

A new method is proposed to divide a spherical surface into equal-area cells. The method is based on dividing a sphere into several latitudinal bands of near-constant span with further division of each band into equal-area cells. It is simple in construction and provides more uniform latitude step be-tween latitudinal bands than other methods of isolatitudinal equal-area tessellation of a spherical surface.

The current state of the link problem between radio and optical celestial reference frames is considered. The main objectives of the investigations in this direction during the next few years are the preparation of a comparison and the mutual orientation and rotation between the optical {\it Gaia} Celestial Reference Frame (GCRF) and the 3rd generation radio International Celestial Reference Frame (ICRF3), obtained from VLBI observations. Both systems, ideally, should be a realization of the ICRS (International Celestial Reference System) at micro-arcsecond level accuracy. Therefore, the link accuracy between the ICRF and GCRF should be obtained with similar error level, which is not a trivial task due to relatively large systematic and random errors in source positions at different frequency bands. In this paper, a brief overview of recent work on the GCRF--ICRF link is presented. Additional possibilities to improve the GCRF--ICRF link accuracy are discussed. The suggestion is made to use astrometric radio sources with optical magnitude to 20$^m$ rather than to 18$^m$ as currently planned for the GCRF--ICRF link. In addition, the use of radio stars is also a prospective method to obtain independent and accurate orientation between the Gaia frame and the ICRF.

The Allan variance (AVAR) was introduced 50 years ago as a statistical tool for assessing of the frequency standards deviations. For the past decades, AVAR has increasingly being used in geodesy and astrometry to assess the noise characteristics in geodetic and astrometric time series. A specific feature of astrometric and geodetic measurements, as compared with the clock measurements, is that they are generally associated with uncertainties; thus, an appropriate weighting should be applied during data analysis. Besides, some physically connected scalar time series naturally form series of multi-dimensional vectors. For example, three station coordinates time series $X$, $Y$, and $Z$ can be combined to analyze 3D station position variations. The classical AVAR is not intended for processing unevenly weighted and/or multi-dimensional data. Therefore, AVAR modifications, namely weighted AVAR (WAVAR), multi-dimensional AVAR (MAVAR), and weighted multi-dimensional AVAR (WMAVAR), were introduced to overcome these deficiencies. In this paper, a brief review is given of the experience of using AVAR and its modifications in processing astro-geodetic time series.

A possible method for linking the optical Gaia Celestial Reference Frame (GCRF) to the VLBI-based International Celestial Reference Frame (ICRF) is to use radio stars in a manner similar to that in the linking of the Hipparcos Celestial Reference Frame (HCRF) to ICRF. In this work, an obtainable accuracy of the orientation angles between GCRF and ICRF frames was estimated by Monte Carlo simulation. If the uncertainties in the radio star positions obtained by VLBI are in the range of 0.1-4 mas and those obtained by Gaia are in the range of 0.005-0.4 mas, the orientation angle uncertainties are 0.018-0.72 mas if 46 radio stars are used, 0.013-0.51 mas if 92 radio stars are used, and 0.010-0.41~mas if 138 radio stars are used. The general conclusion from this study is that a properly organized VLBI programme for radio star observation with a reasonable load on the VLBI network can allow for the realization of GCRF-ICRF link with an error of about 0.1 mas.

The goal of this presentation is to report the latest progress in creation of the next generation of VLBI-based International Celestial Reference Frame, ICRF3. Two main directions of ICRF3 development are improvement of the S/X-band frame and extension of the ICRF to higher frequencies. Another important task of this work is the preparation for comparison of ICRF3 with the new generation optical frame GCRF expected by the end of the decade as a result of the Gaia mission.

The highly-accurate optical reference frame GCRF (Gaia Celestial Reference Frame) is expected to be available in several years. By the same time, a new version of radio reference frame ICRF (International Celestial Reference Frame) will be also published. The link of GCRF to ICRF will be defined by means of computation of the orientation angles between the two frames using the common extragalactic objects observed in both radio (VLBI) and optics (Gaia). Taking into account the expected accuracy of ICRF and GCRF of the first tens microarcseconds, the link between them should be defined at a microarcsecond level, which requires using the most accurate algorithms and models. One of such models is the Galactic aberration in proper motions, which is not included in the data processing yet. In this paper its impact on the ICRF-Gaia link is estimated. Preliminary results showed that this impact is at a level of about 1 microarcsecond.

Source position time series produced by International VLBI Service for Geodesy and astrometry (IVS) Analysis Centers were analyzed. These series was computed using different software and analysis strategy. Comparison of this series showed that they have considerably different scatter and systematic behavior. Based on the inspection of all the series, new sources were identified as sources with irregular (non-random) position variations. Two statistics used to estimate the noise level in the time series, namely RMS and ADEV were compared.

During the last years, much attention has been paid to the astrometric implications of the galactic aberration in proper motions (GA). This effect causes systematic errors in astrometric measurements at a microarcsecond level. Some authors consider it so serious that it requires redefinition of the celestial reference system (CRF). We argue that such attention to the GA is too much exaggerated. It is just a small astrometric correction that must be taken into account during highly accurate astrometric and geodetic data processing. The accuracy of this correction depends on accuracy of the Galactic rotation parameters and, for most application, on the accuracy of the rotation matrix between Galactic and equatorial systems. Our analysis has shown that our today knowledge of these two factors is sufficient to compute the GA correction with accuracy of better than 10%. The remaining effect at a level of few tenths microarcsecond/yr is negligible nowadays. Another consequence of introducing the GA correction is necessity to return to classical astrometric modeling of the VLBI-derived extragalactic radio source position by the linear trend model. Changing the current paradigm of VLBI-derived CRF based on assumption of zero motion of radio sources to classical one leads to bias in the radio source positions up to several tens microarcseconds for catalog at epoch J2000.0.

In this paper we outline several problems related to the realization of the international celestial and terrestrial reference frames ICRF and ITRF at the millimeter level of accuracy, with emphasis on ICRF issues. The main topics considered are: analysis of the current status of the ICRF, mutual impact of ICRF and ITRF, and some considerations for future ICRF realizations.

Laboratory of Space Geodesy and Earth Rotation (LSGER) of the Institute of Applied Astronomy (IAA) of the Russian Academy of Sciences has been carrying on, since its creation, the computation of geodynamical products: Earth Orientation Parameters (EOP) and station coordinates (TRF) based on observations of space geodesy techniques: Very Long Baseline Interferometry (VLBI), Satellite Laser Ranging (SLR), Global Positioning System (GPS). Principal software components, used for these investigations, include: package GROSS for processing of SLR observations, package Bernese for processing of GPS observations, package OCCAM for processing of VLBI observations, software for data exchange, and software for combination of space geodesy products.

In this paper an attempt is made to extract a systematic part from the source motions obtained from the position time series provided by several IVS Analysis Centers in the framework of the ICRF-2 project. Our preliminary results show that the radio source velocities and the parameters of the systematic part of the velocity field differ substantially between the source position time series, which does not allow us to get a reliable solution for the coefficients of spherical harmonics.

In this study, we compared results of determination of the orientation angles between celestial reference frames realized by radio source position catalogues using three methods of accounting for correlation information: using the position errors only, using additionally the correlations be-tween the right ascension and declination (RA/DE correlations) reported in radio source position catalogues published in the IERS format, and using the full covariance matrix. The computations were performed with nine catalogues computed at eight analysis centres. Our analysis has shown that using the RA/DE correlations only slightly influences the computed rotational angles, whereas using the full correlation matrices leads to substantial change in the orientation parameters be-tween the compared catalogues.

In this presentation, several publicly available empiric models of the celestial pole offset (CPO) and free core nutation (FCN), included those developed by the author, are investigated and compared each other from different points of view, such as representation of the observation data, FCN parameters variation, prediction accuracy. Based on this study, some practical recommendations are proposed.

In this paper we outline several problems related to the realization of the international celestial and terrestrial reference frames --- the ICRF and ITRF --- at the millimeter level of accuracy, with emphasis on ICRF issues. We consider here the current status of the ICRF, the interrelationship between the ICRF and ITRF, and considerations for future ICRF realizations.

In this paper, a new method of investigation of the external radio source position catalogs RSPCs stochastic errors is presented. Using this method the stochastic errors of nine recently published RSPCs were evaluated. It has been shown that the result can be affected by the systematic differences between catalogs if the latter are not accounted for. It was also found that the formal uncertainties of the source position in the RSPCs correlate with the external errors. We also investigated several topics related to the formal uncertainties and systematic errors of RSPC.

Assessing the external stochastic errors of radio source position catalogues derived from VLBI observations is important for tasks such as estimating the quality of the catalogues and their weighting during combination. One of the widely used methods to estimate these errors is the three-cornered-hat technique, which can be extended to the N-cornered-hat technique. A critical point of this method is how to properly account for the correlations between the compared catalogues. We present a new approach to solving this problem that is suitable for simultaneous investigations of several catalogues. To compute the correlation between two catalogues $A$ and $B$, the differences between these catalogues and a third arbitrary catalogue $C$ are computed. Then the correlation between these differences is considered as an estimate of the correlation between catalogues $A$ and $B$. The average value of these estimates over all catalogues $C$ is taken as a final estimate of the target correlation. In this way, an exhaustive search of all possible combinations allows one to compute the paired correlations between all catalogues. As an additional refinement of the method, we introduce the concept of weighted correlation coefficient. This technique was applied to nine recently published radio source position catalogues. We found large systematic differences between catalogues, that significantly impact determination of their stochastic errors. Finally, we estimated the stochastic errors of the nine catalogues.

Several attempts to discover the FICN signal in VLBI nutation series made during last years failed. In this paper, we present some results of our further steps in this direction, unfortunately not successful either. We investigated several VLBI CPO series by means of spectral and wavelet analysis. It has been shown that there are several periodic signals with close amplitude around the expected FICN period without any prevailing one that can be reliably associated with the FICN. The most interesting for further analysis is a relatively stable oscillation with period of 800 days, which is, however, beyond the intervals predicted in other studies cited above. It seems to be necessary to improve theoretical estimates of the FICN period to make its search in the observational data more promising.

In this work, we analyzed all available determinations of the Galactic rotation parameters $R_0$ and $\Omega_0$ made during last 10 years to derive the most probable value of the Galactic aberration constant $A = R_0 \Omega_0^2/c$. We used several statistical methods to obtain reliable estimates of $R_0$ and $\Omega_0$ and their realistic errors. In result, we have obtained the value of $A = 5.0 \pm 0.3$ $\mu$as/yr as the current best estimate of the GA constant. We suggest that the proposed value of the GA constant can be safely used in practice during coming years.

Based on several tens of Ro measurements made during the past decades, several studies have been performed to derive the best estimate of Ro. Some used just simple averaging to derive a result, whereas others provided comprehensive analyses of possible errors in published results. In either case, detailed statistical analyses of data used were not performed. However, a computation of the best estimates of the Galactic rotation constants is not only an astronomical but also a metrological task. Here we perform an analysis of 53 Ro measurements (published in the past 20 years) to assess the consistency of the data. Our analysis shows that they are internally consistent. It is also shown that any trend in the Ro estimates from the last 20 years is statistically negligible, which renders the presence of a bandwagon effect doubtful. On the other hand, the formal errors in the published Ro estimates improve significantly with time.

Radio astronomy observations of close approaches of the Solar system planets to compact radio sources as well as radio source occultations by planets may be of large interest for planetary sciences, dynamical astronomy, and testing gravity theories. In this paper, we present extended lists of occultations of astrometric radio sources observed in the framework of various astrometric and geodetic VLBI programs by planets, and close approaches of planets to radio sources expected in the nearest years. Computations are made making use of the EPOS software package.

This work is devoted to investigation of low frequency variations in Polar motion (PM). It has been shown that the main PM features can be effectively investigated using not only time series of the Pole coordinates, but also using series of latitude variations obtained from observations at one observatory. Such an approach allows us to increase the length of observation series available for analysis. In our study, we extended the IERS PM series back to 1840. We investigate trends and (quasi) harmonic oscilla-tions with periods from one year to decades. The main results were obtained making use of the Singular Spectrum Analysis. Other methods were also used for specific analysis and independent check. The most interesting results are detection of two new large phase jumps in the 1840s and 2000s, and revealing of 80-year period in the PM variations.

The paper investigates the question of whether or not determinations of the distance between the Sun and the Galactic center R0 are affected by the so-called "bandwagon effect", leading to selection effects in published data that tend to be close to expected values, as was suggested by some authors. It is difficult to estimate numerically a systematic uncertainty in R0 due to the bandwagon effect; however, it is highly probable that, even if widely accepted values differ appreciably from the true value, the published results should eventually approach the true value despite the bandwagon effect. This should be manifest as a trend in the published R0 data: if this trend is statistically significant, the presence of the bandwagon effect can be suspected in the data. Fifty two determinations of R0 published over the last 20 years were analyzed. These data reveal no statistically significant trend, suggesting they are unlikely to involve any systematic uncertainty due to the bandwagon effect. At the same time, the published data show a gradual and statistically significant decrease in the uncertainties in the R0 determinations with time.

In this paper, the current best estimate of a fundamental Galactic parameter Galactocentric distance of the Sun $R_0$ has been evaluated using all the available estimates published during the last 20 years. Unlike some other studies, our analysis of these results showed no statistically significant trend in $R_0$ during this period. However we revealed a statistically valuable improvement in the $R_0$ uncertainties with time of about 0.2 kpc for 20 years. Several statistical techniques have been used and compared to obtain the most reliable common mean of 52 determinations made in 1992--2011 and its realistic uncertainty. The statistical methods used include unweighted mean, seven variants of the weighted mean, and two variants of median technique. The $R_0$ estimates obtained with these methods range from 7.91 to 8.06 kpc with uncertainties varying from 0.05 to 0.08 kpc. The final value derived in this analysis is $R_0 = 7.98 \pm 0.15\,|_{stat} \pm 0.20\,|_{syst}$ kpc, which can be recommended as the current best estimate of the Galactocentric distance of the Sun. For most of the practical applications the value of $R_0 = 8.0 \pm 0.25$ kpc can be used.

The influence of proper motions of sources due to Galactic aberration on precession models based on VLBI data is determined. Comparisons of the linear trends in the coordinates of the celestial pole obtained with and without taking into account Galactic aberration indicate that this effect can reach 20 $\mu$as per century, which is important for modern precession models. It is also shown that correcting for Galactic aberration influences the derived parameters of low-frequency nutation terms. It is therefore necessary to correct for Galactic aberration in the reduction of modern astrometric observations.

Allan variance (AVAR) was first introduced more than 40 years ago as a estimator of the stability of frequency standards, and now it is actively used for investigations of time series in astronomy, geodesy and geodynamics. This method allows one to effectively explore the noise characteristics for various data, such as variations of station and source coordinates, etc. Moreover, this technique can be used to investigate the spectral and fractal structure of the noise in measured data. To process unevenly weighted and multidimensional data, which are usual for many astronomy and geodesy applications, AVAR modifications are proposed by the author. In this paper, a brief overview is given of using of classical and modified AVAR method in astronomy and geodynamics.

Astrometric observations of the radio source occultations by solar system bodies may be of large interest for testing gravity theories, dynamical astronomy, and planetary physics. In this paper, we present an updated list of the occultations of astrometric radio sources by planets expected in the nearest years. Such events, like the solar eclipses, generally speaking, can be only observed in a limited region. The map of the shadow path is provided for the events occurred in regions with several VLBI stations and hence the most interesting for radio astronomy experiments.

From observations with the 6-m BTA telescope at SAO RAS, we have determined spectroscopic redshifts of seven optical objects whose coordinates coincide with those of radio sources from the list of IVS (International VLBI Service for Geodesy and Astrometry). When compared to radio data, the obtained spectra and redshifts provide evidence for reliable identification of four observed objects; the other three require further study. The distances to the sources derived from our measurements will make it possible to refine the current estimates for parameters of cosmological models based on proper motions of these objects, which are determined from geodetic VLBI observations.

A new list of physical characteristics of 4261 astrometric radio sources, including all 717 ICRF-Ext.2 sources has been compiled. Comparison of our data of optical characteristics with the official International Earth Rotation and Reference Systems Service (IERS) list showed significant discrepancies for about half of 667 common sources. We also found that asymmetry in the radio sources distribution between hemispheres could cause significant correlation between the vector spherical harmonics, especially if the case of sparse distribution of the sources with high redshift. We identified radio sources having many-year observation history and lack redshift. This sources should be urgently observed at large optical telescopes. The list of optical characteristics created in this paper is recommended for use as a supplement material for the next International Celestial Reference Frame (ICRF) realization. It can be also effectively used for cosmological studies and planning of observing programs both in radio and optics.

A new list of physical characteristics of 3914 astrometric radio sources, including all 717 ICRF-Ext.2 sources, observed during IVS and NRAO VCS sessions have been compiled. The list includes source type, redshift and visual magnitude (if available). In case of doubt detailed comment is provided. The list of sources with their positions was taken from the Goddard VLBI astrometric catalog with addition of two ICRF-Ext.2 sources. At this stage the source characteristics were mainly taken from the NASA/IPAC Extragalactic Database (NED). 667 sources from our list are included into the IERS list. Comparison has shown a significant difference in characteristics for about half of these 667 common sources. We compiled a list of frequently observed sources without known physical characteristics for urgent optical identification and spectrophotometric observations with large optical telescopes. This presented list of physical characteristics can be used as a supplement material for the ICRF-2, as well as a database for kinematic studies of the Universe and other related works, including scheduling of dedicated IVS programs.

In 2007, a joint IERS/IVS Working Group has been established to consider practical issues of creating the next ICRF generation, ICRF-2. The goal of the WG is to seek after ways to improve the existing ICRF. In this study we investigate a possibility of ICRF improvement by means of using combined ICRF catalogue instead of a catalogue computed in a single analysis centre, even though using most advanced models and software. In this work, we present a new version of Pulkovo combined catalogue of radio source positions computed using the method proposed in \cite{SokMal07}. Radio source catalogues that were submitted in 2007 in the framework of the WG activity were used as input for mutual comparison and combination. Four combined catalogues have been calculated: Two first of them provide a stochastic improvement of the ICRF, and last two of them allow us to account also for systematic errors in the current ICRF version.

In this paper, several issues are considered, related to the construction of the next ICRF generation, ICRF-2. Between them, the following points are touched: ICRF-2 structure, ICRF Core sources selection, and some expected user's requirements.

Free core nutation (FCN) is the main factor that limits the accuracy of the modeling of the motion of Earth's rotational axis in the celestial coordinate system. Several FCN models have been proposed. A comparative analysis is made of the known models including the model proposed by the author. The use of the FCN model is shown to substantially increase the accuracy of the modeling of Earth's rotation. Furthermore, the FCN component extracted from the observed motion of Earth's rotational axis is an important source for the study of the shape and rotation of the Earth's core. A comparison of different FCN models has shown that the proposed model is better than other models if used to extract the geophysical signal (the amplitude and phase of FCN) from observational data.

In this paper, a new geometry index of Very Long Baseline Interferometry (VLBI) observing networks, the volume of network $V$, is examined as an indicator of the errors in the Earth orientation parameters (EOP) obtained from VLBI observations. It has been shown that both EOP precision and accuracy can be well described by the power law $\sigma=aV^c$ in a wide range of the network size from domestic to global VLBI networks. In other words, as the network volume grows, the EOP errors become smaller following a power law. This should be taken into account for a proper comparison of EOP estimates obtained from different VLBI networks. Thus performing correct EOP comparison allows us to accurately investigate finer factors affecting the EOP errors. In particular, it was found that the dependence of the EOP precision and accuracy on the recording data rate can also be described by a power law. One important conclusion is that the EOP accuracy depends primarily on the network geometry and to lesser extent on other factors, such as recording mode and data rate and scheduling parameters, whereas these factors have stronger impact on the EOP precision.

During close angular approaches of solar system planets to astrometric radio sources, the apparent positions of these sources shift due to relativistic effects and, thus, these events may be used for testing the theory of general relativity; this fact was successfully demonstrated in the experiments on the measurements of radio source position shifts during the approaches of Jupiter carried out in 1988 and 2002. An analysis, performed within the frames of the present work, showed that when a source is observed near a planet's disk edge, i.e., practically in the case of occultation, the current experimental accuracy makes it possible to measure the relativistic effects for all planets. However, radio occultations are fairly rare events. At the same time, only Jupiter and Saturn provide noticeable relativistic effects approaching the radio sources at angular distances of about a few planet radii. Our analysis resulted in the creation of a catalog of forthcoming occultations and approaches of planets to astrometric radio sources for the time period of 2008-2050, which can be used for planning experiments on testing gravity theories and other purposes. For all events included in the catalog, the main relativistic effects are calculated both for ground-based and space (Earth-Moon) interferometer baselines.

In this study, three empirical Free Core Nutation (FCN) models developed to the present time, MHB2000, Malkin's and Lambert's ones, are compared on the basis of representation of variations of the FCN amplitude and phase predicted by these models. It is possible to conclude, that the model of the author provides the most realistic representation of the FCN variations. However, the specified models are based on representation about single FCN rotational mode. At the same time, some results of processing of the VLBI observations made during last years, specify possible presence of two close FCN periods. A theoretical explanation to presence of a second FCN frequency FCN has been given by G. Krasinsky in his theory of rotation of the Earth with two-layer liquid core, ERA2005. In the present work, for more detailed studying this phenomenon, the IVS time series of celestial pole offset, and also those predicted by the ERA2005 theory, have been investigated by means of several statistical methods which confidently show presence of two fluctuations in nutational movement of an Earth's rotation axis with the periods about -452 and -410 solar days.

A new international Pilot Project for the re-determination of the ICRF was initiated by the International VLBI Service for Geodesy and Astrometry (IVS) in January 2005. The purpose of this project is to compare the individual CRF solutions and to analyze their systematic and random errors with focus on the selection of the optimal strategy for the next ICRF realization. Eight radio source catalogues provided by the IVS Analysis Centers GA, SHAO, DGFI, GIUB-BKG, JPL, MAO NANU, GSFC, USNO were analyzed. In present study, four analytical models were used to investigate the systematic differences between solutions: solid rotation, rotation and deformation (IERS method), and expansion in orthogonal functions: Legendre-Fourier polynomials and spherical functions. It was found that expansions by orthogonal function describe the differences between individual catalogues better than the two former models. Finally, the combined CRF was generated. Using the radio source positions from this combined catalogue for estimation of EOP has shown improvement of the uncertainty of the celestial pole offset time series.

In this paper, we propose to use the scatter of celestial pole offset (CPO) series obtained from VLBI observations as a measure of the accuracy of the celestial reference frame (CRF) realizations. Several scatter indices (SIs) including those proposed for the first time are investigated. The first SI is based on analysis of residuals of CPO series w.r.t. Free Core Nutation (FCN) model. The second group of SIs includes Allan deviation and its extensions, which allow one to treat unequal and multidimensional observations. Application of these criteria to several radio source catalogues showed their ability to perform a preliminary assessment of the quality of radio source catalogues, and 2D Allan deviation estimate seems to be a most sensitive measure. However, in common case, the sensitivity of tested criteria is yet not sufficient to distinguish between radio source catalogues of the highest quality. Proposed extensions of Allan deviation, weighted and multidimensional, can be effectively used also for statistical analysis of astro-geodetic and other time series.