this paper is based on a talk at the jenam2003 meeting in budapest , hungary in august 2003 , and will cover several topics . in  2 , reference systems are defined and a brief historical sketch is given . in  3 , the fundamental synthesis formula used for the determination of radio sources positions , and two basic observed quantities , the phase and the group delay , are discussed .
the international celestial reference frame ( icrf ) is described in some detail in  4 .
the measurement of relative positions to obtain the motion of radio sources , and calibrator choice concerns are then given in  5 , and in  6 , the description of four experiments highlight current techniques .
a brief conclusion is given in  7 , and an appendix concerning the controversy of the speed of gravity ends the paper .
in order to study the position and motion of objects , a reference frame is needed . in a normal three - dimensional space ( a flat minkowski space ) , three coordinate numbers specify the location of an object .
these coordinate numbers are defined on a reference frame which is given by an origin ( zero point ) and the direction of two of three orthogonal axes . in principle
, one can choose an arbitrary origin and axes directions , but some reference frames are _
better _ than others .
if a rotating or accelerating reference frame is used , then the objects in this frame will have peculiar motions which can not be understood by application of dynamical laws to their motions .
for example , within the earth - centered reference frame ( right ascension and declination ) , the stars show an annual parallactic motion which depends on their distance and direction in the sky .
these apparent motions are , of course , caused by the earth orbital motion ; hence our adopted reference frame was rotating in space .
in fact , the best determination of the quality of a reference frame is to look for un - dynamic and strange correlated motions of objects in the frame , and try to ascertain what frame motion is causing the aberrant behavior of the stars .
non - rotating and non - accelerating reference frames are called _ inertial _ frames , and all of the laws of mechanics are frame - independent as long as an inertial reference frame is used .
but , how can these frame properties be determined without some knowledge of an _ absolute space _ with no motion and rotation  a preferred reference frame ?
it is generally agreed that such a frame does not exist according to mach s principle and einstein s theory of general relativity ( gr ) ; nevertheless , we can choose a reference frame which is quasi - inertial in that its origin has negligible acceleration and the reference axes have little rotation  at least to the accuracy of the measurements .
for the studies of celestial bodies , the most convenient celestial reference frame has its origin at the barycenter position of the solar system .
the motion of the sun around the center of the galaxy , the gravitational bending from the bulge of our galaxy , and the motion of the galaxy in space produce both large static corrections ( which can be removed from the source positions ) and small differential corrections over the sky at the 0.01 mas level ( sovers , fanslow & jacobs @xcite ) .
the precise position of the barycenter is in error by several km because of the uncertainty in the masses of the major solar system members , predominantly the sun and jupiter .
the two axes of the barycenter frame ( usually called the pole and principle plane directions ) of the celestial reference frame are generally defined by the orientation of the earth in space at some specified time . since most observations are made from a slowly moving crust on the surface of the earth which also orbits the sun and rotates somewhat non - uniformly around a wobbling axis , the difference of the earth orientation at the observation time to that at the fiducial time is a major challenge in obtaining the accurate position of celestial objects on the adopted quasi - inertial reference frame .    after defining a suitable reference frame
, a method must be developed by which observations of celestial objects can be located on this frame .
this is commonly done in two ways : a _ kinematic frame _ is one in which the positions of a suitable set of celestial objects are known at any time ( some of the distant objects can be considered fixed in the sky ) .
thus , by comparing observations of target object with these fiducial objects , the target coordinates can be placed on the adopted reference frame . in principle , only two fiducial objects are needed to specify the frame completely , but observations of both objects and the target will often be impossible ; hence , the number of fiducial objects should be much larger so that several can be accessible near the time and position of the target observation .
a _ dynamic frame _ is defined , not by a fiducial set of objects in the sky , but by an accurate ephemeris of the solar system bodies and the earth rotation and orientation .
this information can be obtained , for example , by the motion of the stars which reflect the earth motion . in the 1970 s and 1980 s ,
the use of spacecraft tracking , laser - ranging to non - earth surfaces and planetary radar observations added significant accuracy in the determination of the orbital , rotational and orientation motions of the earth .
of course , the realization of a reference frame can consist of a combination of both kinematic and dynamic information .
newcomb s studies of the relative motion of stars in the @xmath0 century produced accurate values of the precessional motion of the earth s pole , the length of a day and year , by measuring the transit time of about 1000 bright cataloged stars . with these constants , the fk3 catalog ( fundamental katalog 3 ; kopff @xcite ) of the star positions and the assumed precession and nutation constants defined a celestial reference system which could be used to determine the coordinates of any star or solar system object to an accuracy of @xmath1 .
the extension of catalogs to more than 10000 bright stars and more accurate proper motion determinations led to the compilation of the fk4 catalog ( fricke @xcite ) along with better precessional constants .
the accuracy of this improved reference frame increase to about @xmath2 , with some degradation in the southern sky .
both frames are dynamical because they relied on the modeling of the earth s orbital and spin motion using the residual motion of many stars .
the fk4 system was not keeping up with the accuracy of the observations , and created pressure for continuing revision of the fundamental constants of the system .
thus , the fk5 system ( fricke  @xcite ) was developed around 1980 .
the procedures which improved the fk5 frame were better measurements of the proper motion of thousands of stars , and the specific assumption that distant galaxies were fixed in the sky .
this added a kinematical frame component to fk5 reference frame definition and improved the accuracy to about @xmath3 by 1990 .
other astronomical observations  the motion of near - earth asteroids and occultations of stars by planets  also increased the accuracy of this frame .
but , the seeds were planted for a kinematical approach to future reference frames .    beginning around 1980 , the jet propulsion laboratory ( jpl ) determined a celestial reference using many types of observations , mostly of the dynamic type .
the data were gathered from : optical observations of the planets , viking spacecraft range observations , radar observations of mercury , venus and mars , lunar - laser ranging , asteroid perturbations , and precision lunar modeling .
this de200 ephemeris was tied to the fk5 system by proper overlap of planetary objects ( standish @xcite ) .
the technique of very long baseline interferometry ( vlbi ) demonstrated in the 1970 s that many radio sources have significant emission within a component ( radio core ) less than 1 mas .
these sources are identified with quasars ( galaxies with intense point - like stellar nucleii ) at distance on the order of 1 gpc .
thus , it is an excellent assumption that these sources are virtually fixed in the sky , and they could define an inertial reference frame to much higher accuracy than the optical - based fk5 frame or the planetary ephemerides .    over the next 20 years , improvements in the stability and observing procedures of vlbi , with more precise modeling of the earth rotation ,
nutation , orientation and planetary motions , led to quasar position accuracy approaching 1 mas over the sky . in the early 1990 s
the astronomical community formed working groups in order to define an international celestial reference system ( icrs ) ( ferraz - mello  @xcite ; feissel & mignard @xcite ) which could best utilize the accurate astrometric theory and results .
the icrs was a set of rules and conventions , with the modeling required , to define at any time the orientation of the three coordinate axes ( only two are independent ) , located at the barycenter of the solar system .
the axis directions were fixed relative to a suitable number of distant extragalactic sources . for reasons of continuity with the fk5 system , the icrs - defined pole direction at epoch j2000.0
was set to the fk5 pole at that epoch , and the origin of right ascension was defined by a small radio component in the source 3c273 with an accurate position measurement determined by a lunar occultation ( hazard  @xcite ) .
the catalog of positions of the fiducial objects needed to realize the icrs is called the international celestial reference frame ( icrf ) ( ma  @xcite ) .
since many telescopes around the globe participate in vlbi observations , a complementary terrestrial reference frame was needed .
thus , the international terrestrial reference system ( itrs ) was formulated to describe the rules for determining specific locations on or near the earth . in order to realize the itrs ,
the international terrestrial reference frame ( itrf ) consists of a catalog of about 50 fiducial locations of the earth .
global position service ( gps ) observations over the last 10 years have also provided high accuracy in the determination of the itrf ( see http://www.iers.org/iers/products/itrf ) .
this paper will describe two aspects of astrometry : all - sky _ absolute _ astrometric techniques to determine the accurate positions used in defining the icrf , and _ relative _ astrometric techniques used to determine the motion of individual objects .
however , a brief introduction to radio interferometry is needed in order to understand the calibration methods and the different observing strategies used for absolute and relative astrometry .
a radio source signal , intercepted by two telescopes pointed in the same direction of sky , is essentially identical except for the different travel time ( delay ) of the radio wave from the quasar to each telescope .
the correlation ( multiplying and averaging ) of the two signals produces a response called the spatial coherence function .
its amplitude is related to the strength of the source and its angular size ; its phase is related to the delay of the signal between the telescope , with a minor contribution from the source structure . for an array of many telescopes , each telescope - pair ( baseline ) gives an independent measure of the spatial coherence function .    using a  priori calibrations of the telescope amplification properties
, the coherence amplitude can be converted into true energy units ( flux density ) , and is denoted as the visibility amplitude .
the coherence phase is a rapidly changing function of time because the earth rotation changes the delay between each baseline .
however , an accurate model of the parameters describing the observations ( the location of the telescopes on the earth surface , the orientation and rotation of the earth in space , the position of the radio source , propagation delays in the troposphere and ionosphere , etc ) can be calculated at any time for any baseline . when this model coherence phase is subtracted from the observed spatial coherence phase , a slowly changing residual phase , the observed visibility phase , is obtained .
the basic components of the observed visibility phase , @xmath4 observed at frequency @xmath5 for source @xmath6 at time @xmath7 , between telescopes @xmath8 and @xmath9 are      where the structure phase of the @xmath11 source is given by @xmath12 , and is zero for a point source .
@xmath13 is the model location of telescope @xmath8 , and @xmath14 is the unknown location offset ; @xmath15 is the model position of the @xmath6th source , and @xmath16 is the unknown position offset , @xmath17 is the residual instrumental and clock delay error for telescope @xmath8 , and @xmath18 is the residual tropospheric propagation delay in the direction to the @xmath11 source for telescope @xmath8 . the residual ionospheric refraction for telescope @xmath8 is @xmath19 .
since this delay is produced by plasma refraction , the phase varies as @xmath20 .
those quantities in bold face are vectors on the ground or directions in the plane of the sky . however , _
the measured visibility phase is only defined between @xmath21 and @xmath22 turn . _    the first line in eq .
( [ eq1 ] ) gives the terms which are associated with the source properties : its structure and position . when combined with the visibility amplitude , the ensemble of visibility amplitudes and phases for all baselines and times is the two - dimensional fourier - pair of the source brightness distribution in the sky ( thompson , moran & swenson @xcite ) .
hence , the fourier transform of the visibility amplitude and phase will produce an image from which positional information can be obtained .    the second line in eq .
( [ eq1 ] ) shows the major error terms which have significantly different spatial and temporal properties .
for example , @xmath23 represents the slowly changing location of the telescope ( with respect to the a  priori model ) which are caused by continental drift , or earth orientation , rotation and nutation uncertainties .
this error term ( by virtue of the dot product with the source position ) has a period of 24 hours .
the @xmath24-term represents all temporal delay errors associated with a telescope , including that from the independent maser clocks , with a typical one hour time variation time scale .
the propagation delay through the troposphere @xmath25 ( or the ionosphere @xmath26 at low observing frequencies ) is the most intractable of the error sources . even with detailed ground meteorological measurements and tropospheric / ionospheric path - length monitoring from gps satellites , the a  priori delay model can be significantly different than the actual path delay above each telescope .
see treuhaft and lanyi ( @xcite ) , niell ( @xcite ) and macmillan & ma ( @xcite ) for more information on modeling the troposphere .    in order to obtain _ absolute _ positions ,
all of the error contributions must be determined and parameterized as accurately as possible , and the necessary observing strategies are described in  4 . on the other hand _ relative _ positions are generally obtained by alternating observations of the target with a nearby fixed radio source ( calibrator ) with known properties , and assuming that the errors affecting the observations of the calibrator are nearly the same as that for the target source . this technique is described in  5 .      for small arrays ( @xmath27 km ) ,
the experimental a  priori model is sufficiently accurate so that cycle ambiguities of the phase are not a problem , and eq .
(  [ eq1 ] ) can be used directly to obtain accurate positions ( wade @xcite ) . however , for most vlbi experiments with baselines well in excess of 1000 km , even the most accurate a priori model now available produce residual delays at telescope baselines of 5000 km which are @xmath28 psec , equivalent to a 4 cm path - length , more than one cycle of phase at @xmath29ghz . at frequencies above 2 ghz ,
the dominant error is caused by the troposphere ; at lower frequencies the ionospheric delay becomes dominant , and delay changes over short times scales and in small patches of sky can exceed 5 cm .
thus , the measured visibility phases have cycle ambiguities and are difficult to use directly for astrometric analysis .
the phase ambiguity problem can be overcome by the measurement of the derivative of the visibility phase with frequency , which is called the group delay , @xmath30 .
it can be determined by observing at many frequencies simultaneously , or nearly simultaneously , to derive the phase slope . as long as some selected observing frequencies
are not too widely separated , the group delay is not affected by phase cycle ambiguities .
the use of the group delay for astrometric work and the need for special calibrations and frequency sampling were first described in detail by rogers ( @xcite ) , and the technique is often called bandwidth - synthesis . the derivative of eq .
( [ eq1 ] ) with respect to @xmath5 is trivial and simply removes the @xmath5 in front of the large bracket .
two additional terms are produced by source structure ( sovers  @xcite ) , which causes errors at the level of 5 psec , and the ionospheric refraction which can be removed by observing with a wide frequency range .
hence , the solution form using the phase or the group delay is identical upon differentiation with time is more complicated .
the phase rate is less accurate than the group delay because of the short - term contamination by the tropospheric phase changes , although at low elevations the rate term is a important contribution to the data .
the rate is also useful for narrow - band line emission which have poorly defined group delays . ] .
however , the group delay can only be adequately measured for relatively strong sources for which an accurate visibility phase can be determined in most of the simultaneously - observed channels .
astrometry using weak sources is not possible ( see  5 for more details ) .
the majority of the vlbi observations used for the icrf comes from observations by the nasa crustal dynamics project ( cdp ) , the united states naval observatory ( usno ) , the jet propulsion laboratory ( jpl ) , and other groups . although astrometric observations began as early as 1971 at jpl , the wide - band , dual - frequency observations which are the basis of the icrf solutions started in 1979 .
this set of data will be collectively designated in this paper as cdp observations .
the observations were made simultaneously at 2.3 ghz and 8.4 ghz with at least four frequency channels ( each with 16 subchannels or 16 delay lags ) at both frequencies .
the spanned bandwidth at 2.3 ghz was about 0.1 ghz , and at 8.4 ghz about 0.4 ghz , from which accurate group delays were determined . by a suitable combination of the two - frequency data ,
the ionospheric refraction was determined and removed , to produce 8.4 ghz data which is then ionospheric - free ( sovers , fanslow & jacobs @xcite ) .
the observing strategy was carefully planned in order to obtain robust solutions for all of the unknown parameters ( ma  @xcite ) . each observing period ( denoted as a session ) lasted 24-hours in order to fully sample the @xmath31 terms in eq .
( [ eq1 ] ) . by about 1990
each session had evolved into a schedule which consisted of about 150 observations in which about 50 good quality sources were observed for one to ten minutes with no set number of scans per source .
sources were scheduled over the entire sky in a relatively short time in order to determine the variable tropospheric refraction which produces the largest source of error .
determination of each telescope clock delay and mean atmosphere zenith - path delay and gradient were determined at hourly intervals .
many different arrays and telescopes have been used over the years , but the observational plan has remained essentially unchanged . the group delay and rates are analyzed by several available software packages ( modest  sovers , fanslow & jacobs @xcite , calc / solve  ryan  @xcite ) in order to obtain all of the parameters associated with the experiment .
some parameters are global , unchanged over the entire period covered by the observations ( source positions , some telescope properties ) , some parameters are slowly variable but considered as constant over each 24-hour session ( precession / nutation , ut1 offset and rate , telescope locations ) , and some parameters are extremely variable with parameter estimates made every hour ( zenith path delays and gradients , telescope clock drifts ) . the typical residual delay error for each 24-hour session after obtaining the best solution is about 25 psec , and is dominated by residual tropospheric delay residuals .
this corresponds to a phase error at 8 ghz frequency of @xmath32 .
cos@xmath33 and @xmath33 ( from ma  1998 , fig .
6a and 6b ) .
* ( right ) * the measured nutation terms , with respect to the iau 1980 model , observed for each session .
the plotted curve is fit to the data using known nutation periods with 18 yr , 9 yr , 1 yr , 0.5 yr , and 14 days period ( from ma  1998 , fig .
1a and 1b).,title="fig:",width=302 ] cos@xmath33 and @xmath33 ( from ma  1998 , fig .
6a and 6b ) . *
( right ) * the measured nutation terms , with respect to the iau 1980 model , observed for each session .
the plotted curve is fit to the data using known nutation periods with 18 yr , 9 yr , 1 yr , 0.5 yr , and 14 days period ( from ma  1998 , fig .
1a and 1b).,title="fig:",width=151 ]    in the fall of 1995 , a combined solution of all previous cdp observations determined thousands of parameters .
the global parameters of the source positions and some telescope properties were held fixed over the whole set of data . the details concerning this complex least - squares fit to determine the most consistent set of source positions
are given elsewhere ( ma  @xcite ) .    the resulting catalog contains three sets of sources .
the ` best ' 212 sources have frequent observations , relatively small position errors , and are dominated by a small - diameter radio component .
they were used to define the orientation of the axes of the icrf .
in addition , there are about 300 sources with a smaller number of observations , less accurate positions , or obviously unstable positions ( due to large source structure changes ) .
about 100 sources with bright optical counterparts were included for potential frame - tie connections .
the sky distribution of the defining 212 defining sources in the icrf is shown in fig .  1 ( left ) .
the distribution is relatively evenly spread north of the equator .
the distribution for the entire set of icrf sources is given in fig .  1 ( right ) .
observations are continuing and an additional 59 good quality sources have been added to the icrf - extension1 catalog .    the distribution of the root - mean - square ( rms ) radio source position error is given in fig .  2 ( left ) .
in order to obtain realistic , but conservative error estimates , the internal errors , generated from the scatter among the sessions , were multiplied by 1.5 and then added in quadrature to a 0.25 mas error .
these errors should then include structure effects over time and frequency .
the accuracy of the icrf realization of the icrs axes is about 0.02 mas .
an indication of the accuracy of the cdp observations can be seen in fig .  2 ( right ) which shows the derived correction to the pole orientation ( nutation ) obtained for each 24-hour cdp session .
it was fit to a realistic model known from geophysical and gravitational considerations with an rms of the fit of about 0.3 mas .
this is another indication of the accuracy of the frame realization for a 24-hour long session .      since the icrf frame is accurate at the level of a few 0.02 mas and is stable over time scales of decades , it is advantageous to link other frames to it .
the link to the fk5 frame was discussed in connection with continuity of the direction of the pole and origin of right ascension , with the icrf adopting the fk5 pole at a specific epoch of time .
the hipparcos catalog ( perryman  @xcite ) contains observations between 1989.85 and 1993.21 of over 100,000 stars
brighter than 12-mag .
these observations were placed on the hipparcos reference frame ( hg , e.  @xcite ) to an accuracy of 1.0 mas .
the link between the hipparcos optical frame and the icrf was accomplished in several ways .
first , the few bright stars with significant radio emission were observed with the vlbi , the vla and with merlin in order to determine the radio positions with respect to the icrf ( lestrade  @xcite ) .
accurate parallaxes and proper motions were also obtained for these radio emitting stars ( boboltz  @xcite ) . secondly ,
comparison of the hipparcos optical with hst observations linked these stars to extragalactic objects , and then to nearby icrf objects which could be detected by hst .
next , additional optical observations using ground based telescopes also helped link hipparcos positions to nearby galaxies . finally , the earth orientation parameters were obtained from the hipparcos observations and then compared with that obtained from that from nearly concurrent vlbi observations ( lindegren & kovalevsky @xcite ) . with these comparisons ,
the hipparcos optical catalog is now consistent with the icrf to about 0.6 mas in position and 0.25 mas / yr for axis rotation .
however , the position of an individual star and the reference frame tie decreases in accuracy with time because of the proper motion uncertainties of the stars of 1 mas per year .
the link of the icrf to the planetary dynamic ephemeris was accomplished using several techniques .
first , the time of arrival of a pulsar signal is sensitive to the orbital motion of the earth . by comparing the pulsar position derived from timing analysis with that obtained with vlbi observations , improvements in the dynamics of the earth motion
could be obtained ( dewey  @xcite ) .
secondly , occultations of radio sources with solar system objects also tie the icrf to that of the planetary ephemeris .
finally , radio interferometry of spacecraft with respect to nearby icrf source also improved the frame tie between the two systems ; for example , the observations the magellan and pioneer spacecraft with respect to nearby quasars ( folkner  @xcite ) .
comparison of the lunar ranging with vlbi earth orientation results were also important in tying the two frames together ( folkner  @xcite ) .
the resulting jpl de405 ephemeris is connected to the icrf to about 1.0 mas in position and about 0.05 mas / yr in rotation ( standish @xcite ) , and the tie has recently been improved by a factor of two or three ( border 2003 , private communication )      the icrf will be the best realization of a quasi - inertial system at least until 2015 when the space interferometry mission ( sim ) has been in operation for several years .
hence , improvements in the icrs and the accuracy of the icrf positions are needed , and may be obtained in several ways .
first , many more observations of the icrf sources have been made since the original solution was made in 1997 , and improved positions are designated as the icrf - ext1 catalog . a further update , icrf - ext2 ,
is expected ( fey  2004 , in preparation ) .
secondly , vlba observations , in its search for suitable calibrators needed for phase referencing ( beasley  @xcite ; fomalont  @xcite ) , have produced a catalog of over 2000 radio sources , some of which are candidates for additions to the icrf catalog .
also , the number of experiments for sources in the southern hemisphere ( roopesh  @xcite ) has enlarged over the last ten years .
a new defining list of icrf sources may be generated in the next five years , with an expected improvement in the grid accuracy of about a factor of two .
however , discontinuities between the old and new frame definitions must be avoided .
a major source of error in the parameters associated with the icrf is caused by the variable tropospheric delay ; it is the major contributor to the current limit of 0.25 mas accuracy for the good quality sources .
although the troposphere zenith path - delay and gradient are modeled every hour during most 24-hour sessions , the structure and kinematics of the troposphere are extremely complicated , and the typical rms delay scatter , after obtaining the best solution , is about 25 psec at a 5000-km baseline , which corresponds to a positional error of about 0.2 mas at 8 ghz .
more sophisticated models , better monitoring of the tropospheric delay using gps satellites , and ground measurements of water vapor emission in the line of sight to a source could decrease this uncertainty to one - half of the present level , but progress has been slow .
an example of the gain in precision by an improvement in the tropospheric modeling is shown in fig .  3 ( left ) .
the plot shows the difference in the derived declination of sources between two reductions of all of the 24-hour sessions in 1995 : one reduction used an azimuthally - symmetric tropospheric delay model over each telescope and the other included a tropospheric gradient term ( which is now routinely used ) .
the differences as a function of source declination are up to 0.4 mas and show that systematic errors would have been present in the icrf frame if the simpler tropospheric model had been used .
any remaining systematic error could be as large as 0.1 mas .
another source of error is associated with the evolution of quasars  their changing structure .
while these distant objects are fixed in the sky , the radio emission does not emanate from the galactic nucleus , but from the inner jet region which may be located 0.1 mas from the core at 8 ghz , and could be somewhat variable in position .
another source of error is caused by radio - emitting clouds of magnetic plasma which propagate down the jet and are often much brighter than the radio core nearer the galaxy nucleus .
3 ( right ) shows the apparent position versus time for several quasars which have undergone such evolution .
the typical time - scale of changes are several months to years , and apparent position movements of 0.1 mas are common , and occasional changes as large as 1 mas can occur . with our knowledge of the nature of quasar radio emission , and with frequent monitoring of the radio source emission structure , it should be possible to model these evolutionary changes to remove the effects of these apparent positional changes ( fey & charlot @xcite ; sovers  @xcite ) ) .    finally , the icrf positions are based on radio source properties at 8.4 ghz , but observations and studies are now underway to define a complementary icrf at a frequency of 32 ghz ( lanyi  @xcite ) . at this higher frequency , the ionospheric refraction is smaller and the radio source structures are more compact .
however , radio variability and structure changes will still be a problem .
nasa - jpl plans to use 32 ghz for spacecraft telemetry , beginning in 2005 . with an icrf defined at this frequency , the tie between the planetary ephemerides and the quasar reference reference could be strengthened .
however , such a tie is now perturbed by the limited modeling of the asteroids at the level of 0.2 mas per year which may be the dominant error ( standish & fienga @xcite ) .    the advances in absolute astrometric measurements over the next 20 years ( including space interferometry ) may reach accuracies at the level of 0.01 mas , when a true space - time reference frame must be considered .
this means that the astrometric system must become four - dimensionally based rather than the current system which is three dimensional , but with correction factors associated with the variable gravitational field of the solar system .
also , there is some discussion that the origin of the icrf should be moved to the barycenter of the earth - moon system , rather than that of the solar system , because of the lower gravitation potential and diminished gravitational space distortion near the earth compared with that near the sun ( soffel  @xcite ) .
many astrophysical phenomena are associated with the accurate space motion of celestial objects , as well as their accurate position . for a source in the milky way ,
the parallax is the only direct method of determining the distance , and its proper motion is often related to the evolution and formation of the object . from the orbital motion of a source , many properties of a binary or multiple system , including the detection of planet - like objects , can be made . for solar system objects
the accuracy of radio interferometry is similar to that of lunar - laser ranging and range / doppler measurements , and is complementary since interferometry determines positions on the celestial sphere and the other techniques determine the distance .
finally , the effect of the gravitational field in the solar system on the propagation of radio waves can be measured and compared with that predicted by gr or other theories of gravity .
the determination of the motion of a radio source requires less complicated observations and reductions than that associated with all - sky icrf observations .
the reason is that the motion of a radio source can be obtained by measuring its position with respect to any detectable radio source , called a calibrator , that is _ fixed _ in the sky .
if the calibrator - target separation in the sky is small and observations are switched rapidly between the calibrator and target , then the major error terms , shown in eq .
( [ eq1 ] ) , are similar for the two sources and nearly cancel when their observed visibility phases are differed .
this observation scheme is called _ phase referencing _ ( beasley & conway @xcite ) .    from eq .
( [ eq1 ] ) the phase difference , @xmath34 at any time ( t ) between the target source @xmath35 and the calibrator @xmath36 , for the baseline between telescopes @xmath8 and @xmath9 is @xmath37 + \frac{1}{\nu c}\delta i^{s - c}_{l , m}(t)\end{aligned}\ ] ] the first line shows the phases which dependent on the target and calibrator source structures and positions .
the next two lines contains the differential delay errors .
the small separation of the two sources , @xmath38 , significantly decreases the effect of the telescope - location errors ( in the most general sense ) , for example , by a factor of 50 for a one degree target - calibrator separation .
similarly , @xmath39 is the difference of the tropospheric delay error in the direction of the two sources , and @xmath40 is the ionospheric delay
. the quasi - random short - term delays from small clouds will not cancel particularly well , even for close calibrators , but larger angular - scale deviations from the a  priori model will cancel . this term will be discussed in more detail using multi - calibrator observations . finally , the purely temporal ( mostly clock ) delay variations would cancel precisely if the calibrator and target were observed simultaneously , simultaneously ( honma  @xcite ) . ] .
otherwise , time - interpolation between calibrator observations is needed , and small second order clock errors will remain in the difference .
the switching time between observations of the calibrator and target depend on temporal characteristics of the @xmath24 and @xmath25 terms in eq .
( [ eq2 ] ) .
the longest time for which interpolation between two observations will produce an accurate phase is called the _ coherence _ time , and varies from 30 sec at 23 ghz to 5 min at 1.4 ghz . during periods of inclement
whether or strong ionospheric activity , the coherence time can decrease to 10 sec , making phase referencing impossible unless the calibrator and target are observed simultaneously .
the differential delay error between the calibrator and target are often less than 5 psec . even at a high frequency of 23 ghz ,
the corresponding phase difference in @xmath41 is less than one cycle and ambiguities of the phase are unlikely to occur .
( changes more than one cycle can be used as long as the phase between subsequent calibrator observations are connected properly . )
the group delays can still be obtained by observing at several frequencies ; however , this is not recommended , unless necessary , for several reasons .
first , the precision of the visibility phase is greater than that of the group delay by the factor @xmath42 where @xmath5 is the observing frequency and @xmath43 is the spanned frequency range used to determine the group delay . for a reasonably strong sources , the typical delay error using the measured phase is @xmath44 psec , while that for the group delay is about 10 psec . since the all - sky observations , described for the cdp observations have typical residuals of 25 psec , which are dominated by the troposphere errors , the use of the much more accurate phase data ( if it were possible to sort out cycle ambiguities ) would not provide more accurate solutions than use of the group delay . on the other hand ,
the delay errors from a phase - referencing experiment are often @xmath45 psec , which is less than the inherent group delay accuracy , but not the phase accuracy .
another reason for using the phase , rather than the group delay , is to preserve the imaging capability afforded by the visibility phase .
if the error terms in the second line of eq .
( [ eq2 ] ) can be determined ( or assumed to be negligible ) , then fourier imaging of the residual phase will produce an image of the source with its position relative to that of the calibrator . even if the target source is extremely weak and
can not be detected during a single observation , such fourier imaging of the entire target data set will produce an image from which accurate astrometric information can be obtained .
for these weak sources , the group delay can not be measured .
the two most important properties of a calibrator source are its proximity in the sky to the target , and the strength of a compact radio component which must be detectable within a coherence time . in order for the calibrator phase to be accurately determined in a few minutes , the radio core must contain no less than  20 mjy ( using the vlba at 8 ghz at 64 mhz bandwidth which gives an rms noise of 4 mjy ) . the use of a much stronger calibrator source which is more distant from the target than a fainter , but otherwise suitable calibrator , is not recommended .
the decrease in the measurement phase error from a stronger source is more than balanced by the increase in differential errors between the calibrator and target .
many calibrators contain very extended emission and even most compact radio cores are not true point sources .
the phase effect of this structure is contained in the @xmath46 term in eq .
( [ eq2 ] ) . using the self - calibration techniques for which images of strong sources can be made ( cornwell & fomalont @xcite ) , the structure phase can be obtained and removed .
thus , imaging of all calibrators should be done routinely .
a more serious problem is the variable position offset of the calibrator @xmath47 between observing sessions , since this change is translated directly into an apparent position change of the target .
this problem is associated with extended sources where the maximum radio brightness or centroid brightness moves along the jet and with respect to the ( fixed ) galaxy nucleus ( see fig .  3 ( right ) for some examples ) . even for compact sources ,
it is possible that its position shifts with time since the radio emission probably emanates from the inner part of the radio jet which may not be a stationary position with respect to the galaxy nucleus . when trying to reach astrometric levels @xmath48 mas in a target over a period of time longer than a few months , it is crucial to image the calibrator in order to determine the probable angular change from the calibrator evolution .
only a limited number of suitable calibrators are contained in the icrf catalog .
over 2000 additional calibrators can be found in the vlba surveys of calibrations ( beasley  @xcite ; fomalont  @xcite ) , and most have derived positions which are within 1 mas of the icrf grid because they have been observed in one or more of the cdp 24-hour sessions .
independent searches for calibrator sources , which are near a desired target source , also can be undertaken using preliminary observations of km - sized arrays ( vla , wsrt , atca ) to search for faint , compact flat - spectrum , candidate sources . at 1.4 ghz ,
the field of view of many vlbi arrays are sufficiently large so that a calibrator can usually be found sufficiently close to the desired target , and both can then be observed simultaneously .     away .
the fit is shown by the curve .
* ( right ) * the motion of pulsar b0950 + 08 using a calibrator that is @xmath49 away , but using a multi - frequency technique for removing the differential ionospheric refraction . the best fit proper motion and parallax are shown by the curve.,title="fig:",width=226 ]   away .
the fit is shown by the curve . *
( right ) * the motion of pulsar b0950 + 08 using a calibrator that is @xmath49 away , but using a multi - frequency technique for removing the differential ionospheric refraction .
the best fit proper motion and parallax are shown by the curve.,title="fig:",width=113 ]
in this section several phase - referencing astrometric experiments illustrate the techniques which are currently being used .
the first example for b0656 + 14 illustrates a simple phase - referencing ( in - beam ) experiment .
the second example for b0950 + 08 describes a scheme for removing ionospheric refraction using multi - frequency observations .
the third example for hd8703 illustrates observations over many years to determine the orbital motion of a resolved radio source , and the problems associating with calibrator structure changes over time .
all three of these examples assume that the phase error of the calibrator is a perfect measure of the phase errors of the target .
the last example , the measurement of the deflection of a quasar by the jovian gravitational field , shows that by using more than one calibrator , relative positional accuracies less than 0.01 mas can be obtained from a six - hour observation .
the pulsar psr b0656 + 14 was observed with the vlba at 1.67 ghz for five sessions , separated by about six months each , between epochs 2000.9 to 2002.5 ( briskin  @xcite ) .
a calibrator , j0658 + 1410 with flux density 35 mjy , was found by preliminary vla observations near the pulsar , only @xmath50 away , within the reception area of the vlba antennas , so that both ` in - beam ' calibrator and target were observed simultaneously .
because of the proximity of the calibrator and target , the differential delay errors were assumed to be zero .
the calibrator image was sufficiently point - like ; hence , its structure phase was assumed to be zero and its apparent position was assumed fixed in the sky since it is almost certainly an extra - galactic radio source .
images of the pulsar were then made from @xmath51 using the standard fourier - techniques .
the pulsar position was determined from the location of the peak of the pulsar image , with error estimates .
its peak flux density was 3.6 mjy and was far too weak to be detected in a single observation and to derive a group delay .    using the images of the pulsar made from each session ,
the proper motion and parallax were accurately determined . the motion associated with the parallax motion ( after removal of the proper motion )
is shown in fig .  4 ( left ) .
the parallax @xmath52 mas , corresponding to a distance of @xmath53 pc .
the rms of the fit for each epoch was about 0.6 mas which is impressive for a relatively weak radio source .
another experiment , also a measurement of the motion of a pulsar , did not have the luxury of an available in - beam calibrator . using the vlba , alternating observations of the pulsar b0950 + 08 and calibrator j0946 + 1017 , about @xmath49 away , were made every minute over a period of 7 hours , for 4 sessions between 1998.33 and 1999.85 ( briskin  @xcite ) . with the relatively large calibrator - target separation , the ionospheric refraction at 1.6 ghz was expected to be the dominant source of error . since the pulsar is strong , about 64 mjy ,
the observations of the pulsar and calibrator were made at eight frequencies , simultaneously , between 1.31 ghz and 1.71 ghz .
the differential phase error terms in eq .
( [ eq2 ] ) can then be written more generally as @xmath54 where @xmath55 contains all of the non - dispersive delay terms ( independent of frequency ) , and @xmath56 is the expected ionospheric delay component  between the two sources . by fitting
the differential phase observed over the eight frequencies to a function of the form @xmath57 , both delay terms can be determine , and the ionospheric contribution removed from the observed differential phase . for telescopes more than 1000 km apart , the ionospheric - induced position shift can be as large as 10 mas per degree of source - calibrator separation , changing significantly within a time - scale of a few minutes ; however a shift of 1 mas per degree over five to ten minutes of time is more typical for the source position jitter .
pulsar images , made from the phase before the ionosphere correction , were severely distorted .
the pulsar position was then obtained from the ionospheric - free images , and the position of b0950 + 08 determined in the same manner as for b0656 + 14 . the structure phase of the calibrator j0946 +
1017 was determined at each observation epoch and removed from the differential phase .
the calibrator appearance did not change significantly , hence it was assumed to be stationary in the sky .
the derived motion of the pulsar is shown in fig .  4 ( right ) .
the proper motion is @xmath58 mas yr@xmath59 , @xmath60 mas yr@xmath59 , and @xmath61 mas .
the gravity probe b ( gp - b ) mission , developed by nasa and stanford university , will fly four precision gyroscopes in earth orbit in order to measure two general relativistic precessions , a geodetic effect and a frame - dragging effect ( turneaure  @xcite ) .
these precessions will be measured with respect to a suitably bright star whose position must be known to an accuracy of @xmath62 mas .
thus , the goal of the associated radio observations was to find a suitable bright star with sufficient radio emission , and to determine its position and motion with the necessary accuracy .
after exploratory observations , the rs can ven binary - star system , hr8703 , was chosen .
the radio emission comes from a main - sequence f - g star and the companions is a fainter k star with a known period of 24.65 days .
this project illustrates some of the problems associated with both the target and the calibrator with variable source structure .
the radio position of hr8703 has now been measured from more than 20 observation sessions since early 1997 , and is still continuing .
the vlba , the vla and the three nasa deep - space network telescopes , at 8.4 ghz , had sufficient sensitivity that accurate group delays could be obtained for the star when it was stronger than 5 mjy , which was most of the time .
the primary calibrator used was 3c454.3 ( 2251 + 158 ) , about @xmath63 from the star .
a second calibrator source , b2250 + 194 , about @xmath49 away , was also used , as a control object to determine the accuracy of the observations .
there are two groups performing the data reduction .
one group is using the differential phase between hr8703 and 3c454.3 , and the preliminary result reported here is from this group ( ransom ( @xcite ) .
another group is analyzing the experiment using the group delay .    using phase - referencing of hd8703 with 3c454.3 as the calibrator
, good quality images of the binary system were made for each session .
for some of the sessions , the radio emission showed two radio components separated by 1 mas , in some cases only one component was observed determined .
the best fit to the parallax and proper motion ( fig .  5 ( left ) ) , was obtained by using the peak position of the radio emission when one component was detected , and using the geometric center when two components were detected .
( this interpretation is consistent with the radio emission mechanism of fron rs can ven stars . )
also , the residual motion of the star , after removing the parallax and proper motion fit , showed a repeating pattern of 24.65 day period and , hence , this motion was associated with the binary orbit of the radio emission ( fig .  5 ( right ) . the estimated error for the proper motion , parallax and orbital elements is about 0.1 mas .
the positions of hd8703 are with respect to that of 3c454.3
. however comparison of the relative positions of 3c454.3 ( see fig .  4 ( right , bottom ) ) and b2250 + 194 over the experiment period show relative motion between them in the order of 0.3 mas with a time scale of one year . since both calibrators have extended , evolving structure , it is difficult to determine which source contributes to the apparent motion , and some of this shift may appear in the positions derived for hd8703 .
taken over the five years of observations , these position errors lead to an proper motion uncertainty of @xmath64 mas yr@xmath59 , which is still within the limited needed for gp - b mission . however , a more detailed analysis of the calibrator changes and other phase error terms in eq .
( [ eq2 ] ) is currently underway .
on september 8 , 2002 when jupiter passed within @xmath65 of a background quasar j0842 + 1835 , its gravitational field produced an apparent change ( deflection ) in the position of the quasar . at closest approach , according to gr , the deflection has two major components : a radial deflection of 1.19 mas , and a smaller deflection in the direction of motion of jupiter of 0.051 mas .
this smaller deflection component is caused by the gravitational aberrational produced by the relative motion of jupiter and the earth , and is related to the speed of propagation of gravity ( kopeikin @xcite ; frittelli @xcite ) .
thus , by measuring this aberrational component to an accuracy less than 0.01 mas , an estimate of the speed of gravity could be obtained . however ,
this precision was about a factor of three to five better than had been previously obtained ( see previous experiments as good examples ) .
such a gain in positional accuracy could only be obtained by determining the phase error terms associated with eq .
( [ eq2 ] ) more accurately .
although phase referencing between a calibrator and target removes much of the temporal dependence of the phase error , any coherent angular dependence in the sky will produce systematic phase errors between the target and calibrator , and produce distorted images and systematic errors in the position of the peak of the emission .
this error can be diminished by choosing a calibrator closer to the target ( if you can find one ) , or by using _ more than one calibrator _ to determine the angular dependence of the phase error .
preliminary vlba test observations in 2001 suggested that such a multi - calibrator scheme was effective .
additional details on multi - source calibrations are given by fomalont ( @xcite ) .
the design of the jupiter deflection experiment is shown in fig .  6 ( left ) .
the observations were made at 8.4 ghz with the vlba and effelsberg , germany telescope . in order to determine the phase errors toward the target source j0842 , two calibrators , j0839 and j0854 on opposite sides of the target , were observed .
each source was observed for 1.5 min in turn , and the trio of scans were repeated every 4.5 min for ten hours over the day .
although the deflection of j0842 was strongest on 2002 september 8 , identical observations on five days ( september 4 , 7 , 8 , 9 , 12 ) were made in order to obtain sufficient redundancy to estimate realistic positional errors of the deflection .
more details concerning the design and analysis of this experiment are given by fomalont & kopeikin ( @xcite ) .
6 ( right , top ) shows the measured phase for each of the three sources for a typical observing day and baseline .
the overall temporal behavior of the phases , dominated by the clock and gross troposphere delay errors , are similar for the three sources . however , after ut=17h the phases for the three sources show displacements which are consistent with a simple phase gradient error covering the three sources : the displacement for j0839 from j0842 is in the opposite sense and about 25% of that between j0854 and j0842 , as expected from their relative positions in the sky .
thus , the interpolation of the phase of j0839 and j0854 ( weighted by their inverse distance from j0842 ) gives a very good estimate of the phase associated with j0842 .
the corrected phase of j0842 , @xmath66 , after the two source calibration , is shown in fig .  6 ( right , bottom ) .
the phase scatter is small with a rms less than 0.05 cycle , or about 0.040 mas .
the slight offset from zero phase is associated with the sum of the position offsets from the a  priori values , as seen in the expression for @xmath66 , @xmath67 \nonumber \\
+ \hbox{terms of 2$^{nd}$ order in time and angle } \nonumber\end{aligned}\ ] ] where the 0.80 and 0.20 are the calibrator weighting factors for this particular configuration .
the position of the j0842 can be obtained by fourier inversion of the corrected phase and measured visibility amplitude , or by a least - square fit of the phases if the source is sufficiently strong .
the resultant displacement of the radio from the image center is a measure of the linear combination of offset positions ( from the a  priori positions ) , shown in eq .
( [ eq4 ] ) . assuming that the calibrator positions have not varied
, any changes in the target position with time reflect a change in the target position .
the goal of the september 2002 observation was to determine the speed of propagation of gravity , @xmath68 , from the determination of the aberrational part of the quasar deflection , as outlined in fig .
7 ( left ) . for @xmath69 ,
the expected deflection in the direction of jupiter s motion is 0.051 mas at closest approach .
more generally , the deflection is proportional to @xmath70 .
a graphical display of the experimental determination of this deflection term is given in fig .
7 ( right ) which shows the measured position of j0842 every hour on september 8 , relative to the other four observing days .
these positions were determined from images of the residual phases and the observed visibility amplitudes .
the measured aberrational deflection , averaged over september 8 , was @xmath71 mas , which gives a deflection that is @xmath72 times that predicted by gr .
the error estimates are consistent with the position of the target on the four off - jupiter days when no change of position is expected .
when converted into the speed of propagation of gravity , the result is @xmath73c .
a similar analysis on the data using only one calibrator , j0839 , gave a positional error of 0.027 mas compared with 0.009 mas achieved using a second calibrator .
after a discussion of inertial frames and reference systems , the methods for determining the position of radio sources were discussed .
first , all - sky observations used for the icrf were used to determine absolute positions of radio sources which formed the definition of the celestial inertial reference system now in use .
secondly , phase referencing was described using four different experiments to illustrate several aspects of measuring the motion of radio sources to accuracies approaching 0.010 mas per session .
although the experimental result is not in doubt , there is disagreement over the correct role of gravity in this experiment .
kopeikin ( @xcite ) , using the einstein equations directly , interprets the retarded deflection component as a measure of the speed of propagation of gravity ( see also kopeikin & fomalont ( @xcite ) ) . will ( @xcite ) concludes , however , that the retarded deflection component is associated with the aberration of light , not of gravity .
other interpretations are given by asada ( @xcite ) and samuel ( @xcite ) .    at the jenam2003 meeting in budapest ,
a lively give - and - take session took place during the late afternoon of august 28 .
there were about 20 - 30 participants .
the main points of view were presented by g. schfer ( institute for theoretical physics , friedrich schiller university , jena , germany ) , sergei kopeikin ( university of missouri - columbia , mo , usa ) , and gabor lanyi ( jet propulsion laboratories , pasadena , ca , usa )    schfer supported the point of view of asada and will that the experiment measured the speed of radio waves used in observations .
the main arguments were :    1 .
gravity and light move along null geodesics and their speed can not be disentangle in the experiment ; 2 .
the speed of gravity is associated with accelerated motion of light - ray deflecting bodies and thus , only with gravity waves ; 3 .
the propagation of gravity effect shows up only in terms of order of @xmath74 in the metric tensor but the experiment is sensitive to @xmath75 terms only .
kopeikin interpretation is that that the experiment measured the speed of propagation of gravity in the near zone of the solar system . in answer
to schfer s points , kopeikin pointed out that :    1 .
the relativistic deflection of light by jupiter depends on the euclidean product of two null vectors at the point of observation  the vector of light propagation from the quasar and that of gravity propagation from jupiter .
because the quasar and jupiter are located at different positions on the sky , the gravity propagation direction can be separated from the light propagation direction and , thus , the speed of gravity could be measured .
2 .   there are different types of gravitational fields according to algebraic classification b. y. petrov
. only gravitational fields of type n ( plane gravitational waves ) are generated by the accelerated motion of bodies , and they exist only in the radiative zone of the solar system .
other types of gravitational fields are more general and can be associated with the velocity of the bodies .
they dominate in the near zone of the solar system where the jovian deflection experiment has been done .
nevertheless , gravitational fields of the other petrov s types also propagate with finite speed which can be measured .
it is true that propagation of gravity effect shows up only in terms of order of @xmath74 in the metric tensor .
but the metric tensor is not directly observable quantity in the jovian deflection experiment because it depends on specific choice of coordinates .
therefore , any reference to the metric tensor in connection with measured gravitational effects is irrelevant .
what is observed is the phase of radio waves coming from the quasar which is perturbed by gravitational field of jupiter from its retarded position , in accordance with einstein s equation predicting that gravity travels with finite speed equal numerically to the speed of light in vacuum .
this relativistic effect of the retardation of gravity in the electromagnetic phase is gauge - independent and shows up already in terms of order of @xmath75 which we have observed .
the sections on reference frames and the icrf were significantly improved from the suggestions and corrections by chris jacobs .
i also thank chopo ma for his contributions , and ryan ranson for the preliminary hr8703 results .
ryan , j.w . ,
clark , t.a . , ma , c. , gordon , d. , capretee , d. & himwich , w.e . 1993 , in contributions of space geodesy to geodynamics : earth dynamics , ed d.e .
smith & d. l. turcotte ( washington : am.geophys .
unions ) , 37