Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.gao.spb.ru/english/personal/malkin/publ/zm08j.pdf Äàòà èçìåíåíèÿ: Tue May 11 13:11:02 2010 Äàòà èíäåêñèðîâàíèÿ: Tue Oct 2 03:09:28 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: http astrokuban.info astrokuban

J Geod (2008) 82:325­329 DOI 10.1007/s00190-007-0181-x

ORIGIN AL AR TICLE

On the accuracy assessment of celestial reference frame realizations
Z. Malkin

Received: 4 March 2007 / Accepted: 19 July 2007 / Published online: 24 August 2007 © Springer-Verlag 2007

Abstract The scatter of the celestial pole offset (CPO) time-series obtained from very long baseline interferometry (VLBI) observations is used as a measure of the accuracy of celestial reference frame (CRF) realizations. Several scatter indices (SI), including some proposed for the first time, are investigated. The first SI is based on residual analysis of CPO series with respect to a free core nutation (FCN) model. The second group of SIs includes Allan deviation and its extensions, which allow the treatment of unequally weighted and multidimensional observations. Application of these criteria to several radio source catalogues (RSCs) showed their ability to perform a preliminary assessment of the quality of each RSC. The 2D Allan deviation estimate seems to be the most sensitive measure. The proposed extensions of Allan deviation, weighted and multidimensional, can also be used for the statistical analysis of other time-series. Keywords Reference systems · Celestial reference frame · VLBI · Statistical data analysis · Accuracy assessment

1 Introduction Very long baseline interferometry (VLBI) is the base technique for the realization of the International Celestial Reference Frame (ICRF) given as a set of radio source coordinates (Ma et al. 1998). Improvement of the ICRF's accuracy is one of the primary tasks of the International VLBI Service for Geodesy and Astrometry (IVS), and the assessment of the
Z. Malkin (B) Central Astronomical Observatory at Pulkovo, Pulkovskoe Ch. 65/1, St. Petersburg 196140, Russia e-mail: malkin@gao.spb.ru

ICRF's accuracy is of primary importance in its improvement (SchlÝter and Behrend 2007). Many studies of the accuracy of the ICRF have been performed. Some authors have investigated the time behavior of radio source positions (e.g., Gontier et al. 2001; FeisselVernier 2003; MacMillan and Ma 2007). Others have studied the accuracy of different CRF realizations by means of investigation of stability of the coordinate axis (e.g., Arias et al. 1988; Arias and Bouquillon 2004; Feissel-Vernier et al. 2006). However, existing methods allow us to investigate only differences between celestial reference frame (CRF) realizations or, in other words, radio source position catalogues (RSCs). There is no evident and generally accepted methods to assess the absolute RSC, including ICRF, accuracy. For nearly the past three decades, VLBI uniquely provides highly accurate measurements of the movement of the Earth's rotation axis in space, precession and nutation. Currently, the most accurate precession-nutation model recommended by the International Astronomical Union (IAU) comprises the IAU2000A nutation model and the IAU2006 precession model (see McCarthy and Petit 2004, Chap. 5 and cited papers for details), hereafter referred to as IAU2000A. Observed corrections to the precession-nutation model are called celestial pole offsets (CPOs). One of the main factors, which affect the quality of the CPO obtained from VLBI observations, is the accuracy of the RSC used during data processing. In this paper, possible criteria to assess the RSC accuracy through its impact on the CPO time-series are considered. The criteria tested in this study are based on investigation of the scatter of CPO time-series obtained from VLBI observations as a measure of the accuracy of the CRF realizations. Several scatter indices (SIs), including some proposed for the first time, were applied to the RSC obtained in the framework

123


326

Z. Malkin

of the joint project of the International Earth Rotation and Reference Systems Service (IERS) and the IVS on the next ICRF realization (Ma 2005), as well as the RSC recently submitted by IVS Analysis Centers (ACs) to the IVS Data Centers. The latest ICRF realization, ICRF-Ext.2 (Fey et al. 2004), was also tested. Based on the obtained results, one can make a conclusion on the ability of various SIs to estimate the accuracy (quality) of the CRF realizations provided by VLBI ACs.

depends on results of VLBI data processing, and thus, such an estimate implicitly depends on underlying CRF realizations. 2.2 Allan deviation The Allan variance (AVAR) was originally developed for investigation of noise parameters of frequency standards (Allan 1966; Allan et al. 1991). In the last years, AVAR-based methods are often applied to the analysis of geodetic and geodynamic time-series. For example, the IERS EOP Product Center at the Paris Observatory employs this method for the assessment of statistical characteristics of the EOP series (Gambis 2002, 2004). Furthermore, this method is widely used to investigate stochastic properties of station coordinates and baseline time-series (e.g., Malkin and Voinov 2001; Roberts et al. 2002; Le Bail and Feissel-Vernier 2003; Le Bail 2006). Gontier et al. (2001) and Feissel-Vernier (2003) have used the AVAR to analyze the noise in the time-series of radio source positions. An advantage of using the AVAR to investigate the noise of time-series is its independence, in most of the cases, from the systematic errors, trends and low-frequency variations. Usually, the square root of the AVAR, the Allan deviation (ADEV), is used in geodetic time-series analysis. The ADEV estimate for time series of measurements (observations) y1 , y2 ,..., yn is given by (Allan 1966) ADEV = 1 2(n - 1)
n -1

2 Scatter indices Investigation of the scatter of a geodetic time-series is a powerful tool for the assessment of its quality and statistical characteristics (e.g., Gontier et al. 2001; Gambis 2004; Trenkamp et al. 2004; Feissel-Vernier and Le Bail 2005; MacMillan and Ma 2007). The scatter of the time-series can be computed in different ways. The term "scatter index" (SI) will be used to distinguish different methods. 2.1 Residuals of CPO series Both precision and accuracy of the CPO estimates directly depend on the accuracy of radio source positions used during VLBI data processing (Sovers et al. 1998). There is no external absolutely accurate reference series that may be used for comparison, since only the VLBI technique provides highly accurate CPO measurements. For many years, the scatter of Earth orientation parameters (EOPs) time-series is widely used to compare EOP series computed by various ACs. Here, this will be applied to the comparison of CRF realizations. A similar test was used in Feissel-Vernier et al. (2006), where they compared several RSCs obtained in the course of their work using the weighted root-mean-square (WRMS) of the differences with respect to the IAU2000A model. However, it seems that such a test can be improved using IAU2000A supplemented with a FCN model as a reference. In this case, the residuals between VLBI nutation series and the model (CPO) become smaller by a factor of 2­3, and the test should be more sensitive. Thus, the IAU2000A model was used with the addition of the observed corrections computed by means of smoothing the differences between the IVS combined EOP series and the IAU2000A model.1 Since the FCN is the largest part of the observed corrections, for conciseness, hereafter, the latter will be referred to as the FCN contribution. However, this SI cannot provide an independent estimate of the CPO scatter, since any FCN model obtained from a comparison of observations and the IAU2000A model
1

( yi - yi
i =1

+1

)2 ,

(1)

where n is the number of measurements. Generally, Eq. (1) is valid for equally weighted measurements. However, in geodetic and astrometric practice, we mostly deal with measurements of different precision. In this case, we have the measurements y1 , y2 ,..., yn with the associated uncertainties s1 , s2 ,..., sn . To account for the difference in precision of the measured values, they should be weighted. In our case, to treat unequally weighted data, an extension to the ADEV, weighted Allan deviation (WADEV), is proposed and given by WADEV =
n -1

1 2p

n -1

pi ( yi - yi
i =1

+1

)2 , (2)

p=
i =1

pi ,
-1

pi = (si2 + si2+1 )

,

http://www.gao.spb.ru/english/as/persac/fcn2.dat.

where pi are the weights. Figure 1 displays an example of actual weekly series of a station coordinates, which illustrates that the difference between ADEV and WADEV estimates may be significant.

123


Accuracy assessment of CRF realizations
20 Station height, mm 15 10 5 0 -5 -10 1 3 5 7 9 11 Week number 13 15 17 19 ADEV = 4.18 WADEV = 2.66

327 Table 1 IVS analysis centers providing the RSCs used in this study Abbreviation Full name AUS BKG DGF GSF JPL MAO Geoscience Australia Bundesamt fÝr Kartographie und GeodÄsie, Germany Deutsches GeodÄtisches Forschungsinstitut, Germany NASA Goddard Space Flight Center, USA Caltech/NASA Jet Propulsion Laboratory, USA Main Astronomical Observatory of the National Academy of Sciences of Ukraine Shanghai Astronomical Observatory, China US Naval Observatory, USA

Fig. 1 An example of estimation of classical (ADEV) and weighted (WADEV) Allan deviation for unequally weighted measurements: a case of weekly series of a station height

SHA USN

When the statistical parameters of astrogeodetic timeseries are investigated, we often deal with multidimensional values, e.g., terrestrial coordinates and/or velocities (3D or 6D), celestial coordinates and/or proper motions (2D or 4D). Further extension of the WADEV estimator, weighted multidimensional Allan deviation (WMADEV) is useful for analysis of such types of time-series. In this case, we can consider multidimensional measured values yi = ( yi1 , yi2 ,..., yik ) with the associated k -dimensional uncertainties si = (si1 , si2 ,..., sik ). The main idea of the WMADEV extension of the classical Allan deviation is the consideration ADEV as a statistics of the differences between the adjacent measurements, which can be taken as the Euclidean length between the adjacent measured points considered as the coordinates in some space, one-dimensional in the classical case, and k -dimensional in the extended one. Thus, we can define the k -dimensional Allan deviation as WMADEV = di = | yi - yi
n -1 i1 = +1

No significant difference between results of WMADEV estimates computed with weights given by Eqs. (3) and (4) was found during processing of real data.

3 Comparison of catalogues For this study, several CPO time-series were computed with the OCCAM software in the Kalman filter mode. For all computed CPO series, the processing options were kept the same, except different RSCs were used. During the processing, station and source positions were kept fixed, and only EOPs (universal time and terrestrial and celestial pole coordinates), clocks and troposphere parameters, including troposphere gradients (one value per session), were estimated. 504 IVS R1 and R4 VLBI sessions observed during the period 2002.0­2007.0 were processed. Then, the SI described above was computed for all the series. Three groups of VLBI CRF realizations were included. First, eight RSCs obtained in the framework of the joint IERS-IVS project mentioned above.2 The RSCs were provided by eight IVS ACs listed in Table 1. The organizers of the IERS-IVS project requested the participating ACs to submit two versions of RSCs obtained using all available data since 1979 (version 1), and the data from 1900 only (version 2). The intention was to investigate the impact of data selection on CRF realization. Six RSCs submitted by IVS ACs in 2006­2007 in the framework of their routine operations3 were then used. The latest ICRF realization, ICRF-Ext.2 (Fey et al. 2004), was also included in this comparison. Results of this test are presented in Table 2. From Table 2, the accuracy of all the compared RSCs as estimated by this method is generally about the same level. Nevertheless, some discrepancies can be seen. In particular, the ICRF-Ext.2 shows the worst result, i.e largest
2 3

1 2p |,

n -1

pi di2 ,
i =1

p= pi =

pi ,
k

(3) -1 ( yi - yi
j j +1

)/di

2

(si )2 + (si

j

j =1

j 2 +1 )



,

where | ... | denotes the Euclidean norm. The expression for the weights pi is derived from the error propagation law. One can see, it has a singular point in the case of di = 0, i.e. of equal adjacent measurements. Hence, this case require special treatment, e.g. assigning unit weight for (near-)zero di . To avoid the singularity, a simplified formula for computing pi is pi =
k

-1 (si )2 + (si
j j 2 +1 )



.

(4)

ftp://cddis.gsfc.nasa.gov/vlbi/ivs-special/icrf-next/. ftp://cddis.gsfc.nasa.gov/vlbi/ivsproducts/.

j =1

123


328 Table 2 Scatter indices for the CPO time-series obtained with different RSCs Catalogue FCN ADEV WADEV WMADEV

Z. Malkin

IERS/IVS project catalogues, 1979­2005 AUS1 BKG1 DGF1 GSF1 JPL1 MAO1 SHA1 USN1 AUS2 BKG2 DGF2 GSF2 JPL2 MAO2 SHA2 USN2 Latest catalogues AUS2006B BKG2006C CGS2006A GSF2007A IAA2006A USN2006C ICRF-Ext.2 97 96 95 96 97 95 99 112 112 110 111 114 110 117 103 103 102 103 106 102 108 147 146 144 146 150 144 154 99 96 96 96 98 96 96 96 97 95 96 95 96 95 96 95 113 111 113 111 113 113 114 114 113 112 111 111 110 111 112 111 104 104 104 103 105 104 104 104 104 103 103 103 102 103 104 102 150 147 148 146 149 147 148 147 148 146 146 146 145 145 147 145

IERS/IVS project catalogues, 1990­2005

software and the similar processing options, and the main difference in the data processing strategy was handling of the ICRF defining sources, whose positions were adjusted under the No-Net-Rotation condition in the USN catalog and were kept to the first ICRF in the ICRF-Ext.2 (Alan Fey, 2006, private communication). From our test, one can see that the catalogues of version 2 constructed using only the observations made from 1990.0 on show small but steady improvement of the accuracy. This effect can also be explained by a larger weight of the IVS R1/R4 sessions in the RSC of version 2 with respect to the RSC of version 1. On the other hand, previous studies (Malkin 2004a,b) have shown that both EOP and precession parameters derived from VLBI observations show smaller uncertainty when only data from 1990 on are used for analysis. Comparing the SIs considered in this test, the 2D WMADEV estimate seems to be the most sensitive criterion, which can be seen from the ratio of the smallest and largest SI values in each column of Table 2, i.e. for different RSCs. In other words, the WMADEV shows the most scatter among the different CRF realizations. The same ratio comparison shows that 'FCN' SI is the least sensitive.

4 Conclusion In this paper, the scatter of CPO series obtained from VLBI observations was used as a measure of the accuracy of the CRF realizations. Several scatter indices (SIs) have been investigated. The first SI is based on analysis of residuals of CPO series with respect to the IAU2000A precession-nutation model supplemented with the FCN contribution. It turned out to be less sensitive than the other SIs tested. Besides this, the SI cannot be considered as fully independent, since both the IAU2000A model and FCN contribution determined by previous VLBI data analysis rely on CRF realizations. Conversely, the SIs based on Allan deviation and its extensions, weighted Allan deviation (WADEV) and weighted multidimensional Allan deviation (WMADEV), proposed in this paper, provide an independent estimate of the quality of CRF realizations. The application of these criteria to several RSCs has shown their ability to perform a preliminary assessment of the RSC quality. The WMADEV seems to be most sensitive to differences in different CRF realizations (see Table 2). Proposed extensions of the classical Allan deviation estimator, WADEV and WMADEV, can also be used for statistical analysis of other geodetic time-series, such as EOP, station positions and troposphere parameter time-series. Both WADEV and WMADEV estimates can be also applied, like the classical ADEV estimate, for different sampling intervals

FCN column shows the WRMS of the differences between CPO series computed with given RSC and IAU2000A+FCN reference series, ADEV and WADEV columns contain classical and weighted Allan deviation estimates. The values presented in the FCN, ADEV and WADEV columns are computed as the mean of the corresponding independent estimates for the x p and y p celestial pole coordinates. The WMADEV column contains 2D estimates computed with the Eq. (3), k = 2, for both x p and y p coordinates. Unit: µas

scatter, most probably because in ICRF-Ext.2, the 212 defining source positions were kept the same as in the first ICRF realization. In the analyzed data, 186 defining sources were observed of total number of 507 (36.6%), and there were 307155 observations of defining sources of total number of 953637 (32.2%). Thus, positions of the defining sources used for EOP computations have a large impact on the result. Since using ICRF yields the largest scatter of the CPO series as compared to the other catalogues tested, one can conclude that it might be inadvisable to keep coordinates of the defining ICRF sources in successive ICRF realizations. In particular, such a conclusion follows from comparison of the SI values obtained for the ICRF-Ext.2 and USN RSC. Both the RSCs were computed at the USNO using the same

123


Accuracy assessment of CRF realizations

329 Le Bail K, Feissel-Vernier M (2003) Time series statistics of the DORIS and GPS colocated observations. EGS-AGU-EUG Joint Assembly, Nice, April Le Bail K (2006) Estimating the noise in space-geodetic positioning: the case of DORIS. J Geod 80:541­565 Ma C, Arias EF, Eubanks TM, Fey AL, Gontier AM, Jacobs CS, Sovers OJ, Archinal BA, Charlot P (1998) The international celestial reference frame as realized by very long baseline interferometry. Astron J 116:516­546 Ma C (2005) Steps towards the next radio realization of the ICRS. In: Capitaine N (ed) Proceedings of the JournÈes systÕmes de rÈfÈrence spatio-temporels 2004, Paris, September, pp 3­7 MacMillan DS, Ma C (2007) Radio source instability in VLBI analysis. J Geod 81:443­453 Malkin ZM, Voinov AV (2001) Preliminary results of processing EUREF network observations using a non-fiducial strategy. Phys Chem Earth (A) 26:579­583 Malkin Z (2004a) Comparison of VLBI nutation series with the IAU2000A model. In: Finkelstein A, Capitaine N (eds) Proceedings of the JournÈes systÕmes de rÈfÈrence spatio-temporels 2005, St. Petersburg, September, pp 24­31 Malkin Z (2004b) 5,000,000 Delays­some statistics. In: Vandenberg NR, Baver KD (eds) Proceedings of the IVS 2004 general meeting, Ottawa, February, pp 47­51 McCarthy DD, Petit G (eds) (2004) IERS Conventions (2003). IERS Technical Note No. 32. Verlag des Bundesamts fÝr Kartographie und GeodÄsie, Frankfurt am Main Roberts CA, Morgan P, Rizos C (2002) Allan variance applied to time series baseline results for GPS-based deformation monitoring applications. In: Kahmen H, Niemeier W, Retscher G (eds) Proceedings of the 2nd symposium on geodesy for geotechnical and structural applications, Berlin, May, pp 299­311 Sovers OJ, Fanselow JL, Jacobs CS (1998) Astrometry and geodesy with radio interferometry: experiments, models, results. Rev Mod Phys 70:1393­1454 Trenkamp R, Mora H, Salcedo E, Kellogg J (2004) Possible rapid strain accumulation rates near Cali, Colombia, Determined from GPS measurements (1996­2003). Earth Sci Res J 8:25­33 SchlÝter W, Behrend D (2007) The International VLBI Service for Geodesy and Astrometry (IVS): current capabilities and future prospects. J Geod 81:379­387 Williams SDP (2003) The effect of coloured noise on the uncertainties of rates estimated from geodetic time-series. J Geod 76:486­494

as well as overlapping intervals for more detailed investigation of the various time-series noise (see, e.g., Gambis 2002, 2004; Feissel-Vernier 2003; Williams 2003; Le Bail 2006).
Acknowledgments The author greatly acknowledge all the IVS Analysis Centers providing their RSCs used in this study, and all the IVS components whose hard work made this possible to have got the VLBI data for scientific analysis. I am also very grateful to Volker Tesmer and anonymous reviewers for their valuable comments, which helped to improve the manuscript.

References
Allan DW (1966) Statistics of atomic frequency standards. Proc IEEE 54:221­230 Allan DW, Weiss MA, Jespersen JL (1991) A frequency-domain view of time-domain characterization of clocks and time and frequency distribution systems. In: Proceedings of the 45th annual symposium on frequency control, Boulder, May, pp 667­678 Arias EF, Feissel M, Letsrade J-F (1988) Comparison of VLBI celestial reference frames. Astron Astrophys 199:357­363 Arias EF, Bouquillon S (2004) Representation of the International Celestial Reference System (ICRS) by different sets of selected radio sources. Astron Astrophys 422:1105­1108 Feissel-Vernier M (2003) Selecting stable extragalactic compact radio sources from the permanent astrogeodetic VLBI program. Astron Astrophys 403:105­110 Feissel-Vernier M, Le Bail K (2005) Datum definition within and across techniques: station weight factors reflecting their long-term stability. a talk at the IERS workshop on combination, Vienna, April, http://www.iers.org/documents/wg3/Feissel_LeBail.pdf Feissel-Vernier M, Ma C, Gontier AM, Barache C (2006) Analysis strategy issues for the maintenance of the ICRF axes. Astron Astrophys 452:1107­1112 Fey AL, Ma C, Arias EF, Charlot P, Feissel-Vernier M, Gontier AM, Jacobs CS, Li J, MacMillan DS (2004) The second extension of the International Celestial Reference Frame: ICRF-Ext.2. Astron J 127:3587­3608 Gambis D (2002) Allan variance in Earth rotation time series analysis. Adv Space Res 30:207­212 Gambis D (2004) Monitoring Earth orientation using space-geodetic techniques: state-of-the-art and prospective. J Geod 78:295­303 Gontier A-M, Le Bail K, Feissel M, Eubanks TM (2001) Stability of the extragalactic VLBI reference frame. Astron Astrophys 375:661­ 669

123