 source and the resulting estimated positions as a measure of the corresponding systematic errors for each of the above cases. We have adopted the usual astrometric frequency of $8.4$ GHz for this analysis.

Phase model used in simulation {#model}
==============================

The phase delay and group delay in VLBI are described in @Sov98. The phase $ \phi = \nu \tau$ at frequency $\nu $ is related to the interferometer delay $$\tau = \tau_g + \tau_{trop}
+\tau_{iono}+ \tau_R + \tau_{struc} + \tau_{clk}.$$ Specifically, the geometric delay is: $$\tau_g  = { { { [P][N][EOP]} { \vec b. \vec k
\over c }}}$$ with the precession matrix $[P]$, the nutation matrix $[N]$, the Earth Orientation Parameters matrix $[EOP]$, the baseline coordinates $\vec b$ in the terrestrial frame, the source direction coordinates $\vec k$ computed with source right ascension and declination in the celestial frame. The “retarded baseline correction” to account for Earth rotation during elapsed time $\tau_g$ must also be modelled [@Sov98]. The differential tropospheric delay $\tau _{trop}$ between the two stations is computed with a static tropospheric model and the simple mapping function $ 1/\sin E $ (where $E$ is the source elevation at station) to transform the zenith delay into the line-of-sight delay at each station. The differential ionospheric phase delay $\tau_{iono} = -{{k {\rm TEC}}/ {\nu^2}}$ is related to the total electronic content TEC in the direction of the source at each station. The General Relativity delay $\tau_R$ takes into account light propagation travel time in the gravitational potential of the Sun. The source structure contribution $\tau_{struc}$ can be computed according to the model by [@Cha90] but was not included in our simulations which are for point sources. The clock delay $\tau
_{clk}$ cancels in differenced VLBI phases.

The model above is that implemented in the SPRINT software used for our simulations. It is thought to be complete for narrow-angle astrometry and additional refinements, such as ocean loading, atmospheric loading, etc…, would not make difference into our results. We have not studied
source and the resulting estimated positions as a standard of the correspond systematic errors for each of the above cases. We have adopt the usual astrometric frequency of $ 8.4 $   GHz for this psychoanalysis. 

 Phase model used in simulation { # exemplar } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 The phase delay and group delay in VLBI are described in @Sov98. The phase $ \phi = \nu \tau$ at frequency $ \nu $ is related to the interferometer delay $ $ \tau = \tau_g + \tau_{trop } 
 + \tau_{iono}+ \tau_R + \tau_{struc } + \tau_{clk}.$$ Specifically, the geometric delay is: $ $ \tau_g   = { { { [ P][N][EOP ] } { \vec b. \vec k 
 \over vitamin c } } } $ $ with the precession matrix $ [ P]$, the nutation matrix $ [ N]$, the Earth Orientation Parameters matrix $ [ EOP]$, the baseline coordinates $ \vec b$ in the terrestrial frame, the source direction organize $ \vec k$ computed with source right rise and declination in the celestial frame. The “ retarded baseline correction ” to report for Earth rotation during elapsed prison term $ \tau_g$ must also be modelled [ @Sov98 ]. The differential tropospheric delay $ \tau _ { trop}$ between the two place is computed with a static tropospheric model and the simple mapping affair $ 1/\sin E $ (where $ E$ is the source elevation at station) to transform the zenith delay into the line - of - sight delay at each station. The differential ionospheric phase delay $ \tau_{iono } = -{{k { \rm TEC}}/ { \nu^2}}$ is related to the entire electronic content TEC in the direction of the generator at each post. The General Relativity delay $ \tau_R$ takes into account clean generation travel time in the gravitational potential of the Sun. The reference structure contribution $ \tau_{struc}$ can be computed according to the model by [ @Cha90 ] but was not include in our simulations which are for point sources. The clock delay $ \tau 
 _ { clk}$ cancels in differenced VLBI phase. 

 The model above is that implement in the SPRINT software used for our simulations. It is thought to be complete for narrow - slant astrometry and additional refinements, such as ocean loading, atmospheric loading, etc …, would not make difference into our results. We have not studied
 sokrce and the resulting ertimated positions as e measude of thd corresponding systematic ecrorw for each of the above casds. We havv adopted the ysual astrometric nxequehgy of $8.4$ JHz for this anslysis.

Phasa model used it rilulation {#model}
==============================

The phase delay and gwoup dekaj in VLBI are qescgifed jn @Sov98. The phase $ \phi = \nu \tau$ at fdequencj $\nu $ is related yo the interferometer delaj $$\tak = \tau_g + \tau_{trop}
+\tak_{iono}+ \tau_R + \tau_{frruc} + \tau_{clk}.$$ Specifically, the geomgtric delay is: $$\tau_g  = { { { [P][N][EOP]} { \xec b. \vec k
\over c }}}$$ wljh the precewsion matrix $[P]$, thc nutathon mattix $[N]$, the Eartm Orixntarion Parameters matrie $[EOP]$, the baseline cjordinatev $\rec b$ in the terrestruao frake, tve suyrcd djrxctjon colrdmnates $\vec i$ computed qith source right axcqbsion and decminatijn in the celestial frame. The “retarded baveljne correction” to accoubt for Earth rotation during ejapsed time $\tau_g$ must also be modelled [@Sov98]. The dixfereitkal tvoporphfric delay $\tau _{trop}$ between the two stations ya vokputed with a ftatic troppsohrtic model and jhe sim'me mapping function $ 1/\sin E $ (wherw $E$ is tht soutce elevation at station) to transform tke zenith delay into che line-of-siyht dekay ay each station. The diffzrentizl ionosphegic phase aelay $\tau_{iono} = -{{k {\rm TAC}}/ {\nu^2}}$ is related to the toeal electconic conteng TEV in tre directiln of the source at each statnon. Tve General Relativity delay $\tau_R$ takes invp account libhd pgopagatiou travcl time in the dravitational kotential of tfe Sun. The source structure sontribution $\dwu_{struc}$ can ue computqd axcoreing to ghe model by [@Cna90] but was not inclyded in our simulabions shich are for piinu wources. The clpck dejaj $\teu
_{clk}$ wancels in dhffefenzrd VLCI phastf.

Tme oodek above is that implamenfed in the SPRINT xontware usgd for ouw simulations. It is thought to he cokplxte fot nwrrow-angle astrometry and addifional revincments, such af occan joading, atkospheric loading, etc…, would not make dihference into our resulrs. We have not studnec
source and the resulting estimated positions as of corresponding systematic for each of adopted usual astrometric frequency $8.4$ GHz for analysis. Phase model used in simulation ============================== The phase delay and group delay in VLBI are described in @Sov98. phase $ \phi = \nu \tau$ at frequency $\nu $ is related to interferometer $$\tau \tau_g \tau_{trop} +\tau_{iono}+ \tau_R + \tau_{struc} + \tau_{clk}.$$ Specifically, the geometric delay is: $$\tau_g = { { [P][N][EOP]} { \vec b. \vec k \over c with the precession matrix the nutation matrix $[N]$, the Orientation matrix $[EOP]$, baseline $\vec in the terrestrial the source direction coordinates $\vec k$ computed with source right ascension and declination in the celestial frame. “retarded baseline account for rotation elapsed $\tau_g$ must also [@Sov98]. The differential tropospheric delay $\tau two stations is computed with a static tropospheric and the mapping function $ 1/\sin E $ $E$ is the source elevation at station) to the zenith delay into the line-of-sight delay at each station. The differential ionospheric phase delay -{{k {\rm TEC}}/ {\nu^2}}$ related to the electronic TEC the of the at each station. The General Relativity delay $\tau_R$ takes into account propagation travel time in the gravitational potential of the Sun. structure $\tau_{struc}$ can be according to the model [@Cha90] was not included in which for clock $\tau cancels in differenced VLBI The model above is that in the SPRINT software is thought to be complete for narrow-angle astrometry additional refinements, such as ocean loading, atmospheric etc…, would not make difference into our results. We have not studied
 source and the resulting estiMated positIons aS a mEasUrE of tHe coRresponding sysTEmatIc errors for each of the abOve caSeS. we haVE aDopteD the usuAL aSTRomEtRiC frEqUEnCy of $8.4$ Ghz fOr this aNalysis.

PhaSe mOdEl used in simuLAtIon {#model}
==============================

ThE phAse delay and gRouP delay In vLBi Are deScrIbed iN @Sov98. ThE Phase $ \pHi = \nu \tau$ at FrEQuency $\NU $ is relaTED tO the Interferometer delAY $$\tAU = \tau_g + \tau_{trop}
+\taU_{iono}+ \tAu_r + \TaU_{STruC} + \taU_{clk}.$$ SpecifIcAlly, tHE geometRIc DELAy iS: $$\Tau_g  = { { { [P][N][EOP]} { \vec b. \Vec k
\over c }}}$$ wiTH thE preceSsIon MAtrix $[P]$, The nuTaTIon Matrix $[N]$, the EArth orientatiOn ParaMEters maTRix $[EOP]$, tHe baseLinE coOrdiNAtEs $\Vec B$ iN The TErResTRiaL frame, thE sOuRce diRectION COordInaTes $\vEc k$ coMputed with souRce RighT AscEnsioN and dEcliNaTion iN the ceLestiAl Frame. The “retardeD basEline corrEctIoN” to AcCount FOr EartH roTatIon duriNg elapsED tiMe $\TAU_G$ mUst also be modelled [@SOv98]. tHE dIfferentIal troPOsPhERic delay $\TaU _{trOp}$ beTWEen thE two STaTions is cOmputeD WiTh A static TrOposphErIc mOdeL and tHE simPle mapPing funcTion $ 1/\sIN E $ (where $E$ is the sOUrce elevation AT sTATiON) to tRanSform the zenIth dELay iNto tHE lIne-OF-sighT delaY aT EaCH station. The differenTiAl ionoSpherIc phase delay $\tAu_{iono} = -{{k {\rm Tec}}/ {\NU^2}}$ is relatEd to THe TOtal electronic ConteNt TEC in the DIrection Of the Source at Each statiON. the GenerAl RElaTivIty DELaY $\tau_R$ takes intO ACcouNt Light prOpaGation tRavEl tIme In tHe GravitatiOnal poteNtIaL oF tHe SUn. The SOurce strUcTurE cOntRibutIOn $\tau_{sTruc}$ cAn be CoMpUTed AccordiNG tO THe moDeL bY [@Cha90] But WaS not iNcluDEd iN our simUlations wHicH Are fOr PoInt sourCes. The clock deLaY $\tau
_{clk}$ canCeLs iN diffeRENced VLBI Phases.

The model above is thAT implemEntEd in tHe SPrINT softwAre Used foR ouR SimulaTions. IT is thOuGht TO Be comPLEtE foR nArrow-angle ASTroMetry AnD addItional Refinements, such as oCEan Loading, atmospHerIc loADInG, etC…, WoULd nOt MAke DIFference into our Results. We hAvE NoT studied
 source and the resultingestimatedposit ion s a sa me asur e of the corre s pond ing systematic errorsfor e ac h  oft he  abov e cases .  W e hav ead opt ed th e usu alastrome tric frequ enc yof $8.4$ GHz fo r this ana lys is.

Phase m ode l used  i n s i mulat ion  {#mo del}
= = ====== ========= == = ====== = ===

Th e ph asedelay and group d e la y  in VLBI are d escrib ed in @ Sov 98.  The phase  $  \phi = \nu \ t au $ a t f r equency $\nu$ is relate d  to  the i nt erf e romete r del ay $$\ tau = \tau_ g +\tau_{tro p}
+\t a u_{iono } + \tau_ R + \t au_ {st ruc} +\t au_ {c l k}. $ $Spe c ifi cally, t he  g eomet ricd e l a y is : $ $\ta u_g = { { { [P][N ][E OP]} { \ vec b . \ve c k\o ver c  }}}$$  with  t he precession m atri x $[P]$,the  n uta ti on ma t rix $[ N]$ , t he Eart h Orien t ati on P a ra meters matrix $[EO P] $ ,  t he basel ine co o rd in a tes $\ve cb$in t h e  terr estr i al  frame,the so u rc edirecti on  coord in ate s $ \veck $ co mputed  with so urcer ight ascension and declinati o ni n  t h e ce les tial frame.  The “ret arde d  b ase l ine c orrec ti o n” to account for Eart hrotati on du ring elapsedtime $\tau _ g $  must al so b e  m o delled [@Sov98 ]. Th e differen t ial trop osphe ric dela y $\tau _ { t rop}$ be twe enthe  tw o st ations is com p u tedwi th a st ati c tropo sph eri c m ode land the s imple ma pp in gfu nct ion $ 1/\sin E  $  (w he re$E$ i s  the s ource  ele va ti o n a t stati o n) t o tr an sf ormthe  z enith  del a y i nto the  line-of- sig h t de la yat each  station. The  d ifferentia lion ospher i c  phase d elay $\tau_{iono} = -{{ k  {\rm T EC} }/ {\ nu^2 }}$ is re lat ed tothe totalelectr onicco nte n t  TECi n  t hedi rection of t hesourc eat e ach sta tion. The GeneralR ela tivity delay$\t au_R $ ta kes in t o a cc o unt l ight propagatio n travel t im e  i n the grav i tat io nal pot entialof th e  Sun. T he source  structur econt r i but ion $\tau_ {struc}$  can be c o mpute d  a ccord ing  to th emod el by  [@Cha 9 0]but w as not  i nclude d inou r simula tions which are for poi nt sou rces.  Th e clock d ela y  $\ tau
_{clk }$ c ancels indif fer enced  VL B I pha ses. 
T hem odelabov e  is thati mp lem e n te d in the SP R I N T s oftwa reu sed fo r ou r simulations. It is thought tobe c o m ple tef or n ar row-angle astr ome tr y and addi ti onal refine ments, s uc h  as o cean l oading , atmos p h er i c load ing,  et c…, would  no tm ake dif fe re n ce int o ou rresult s. Weh aven o t studied
 source_and the_resulting estimated positions as_a measure_of_the corresponding_systematic_errors for each_of the above_cases. We have adopted_the usual astrometric_frequency_of $8.4$ GHz for this analysis.

Phase model used in simulation {#model}
==============================

The phase delay and group_delay_in VLBI_are_described_in @Sov98. The phase $_\phi = \nu \tau$ at_frequency $\nu_$ is related to the interferometer delay $$\tau_=_\tau_g + \tau_{trop}
+\tau_{iono}+_\tau_R + \tau_{struc} + \tau_{clk}.$$ Specifically, the geometric delay_is: $$\tau_g  = { {_{ [P][N][EOP]} {_\vec_b._\vec k
\over c }}}$$_with the precession matrix $[P]$, the_nutation matrix $[N]$, the Earth Orientation_Parameters matrix $[EOP]$, the baseline coordinates $\vec_b$ in the terrestrial frame, the_source direction coordinates $\vec k$_computed with_source right ascension and declination_in the celestial_frame. The_“retarded baseline correction”_to account for Earth rotation during_elapsed time $\tau_g$_must also be modelled [@Sov98]. The_differential_tropospheric delay $\tau__{trop}$_between_the two_stations is computed_with_a static_tropospheric_model and the simple mapping function_$_1/\sin E $ (where $E$ is the_source elevation at station)_to_transform the zenith delay_into the line-of-sight delay at_each station. The differential ionospheric phase_delay $\tau_{iono}_= -{{k_{\rm TEC}}/ {\nu^2}}$ is related to the total electronic content TEC_in the direction of the source_at each station. The_General Relativity_delay_$\tau_R$ takes into_account_light propagation_travel time in the gravitational potential of_the Sun._The source structure contribution $\tau_{struc}$ can_be computed according to_the_model by [@Cha90] but was not_included in our simulations which are_for point sources. The clock_delay_$\tau
_{clk}$_cancels in differenced VLBI phases.

The_model above is that implemented in_the SPRINT software_used for our simulations. It is thought_to_be complete for narrow-angle astrometry and_additional_refinements, such as ocean loading, atmospheric_loading,_etc…,_would not make difference into_our results. We have not studied