<div dir="ltr">Hi Adam and David,<div><br></div><div>Based on your mails and references, and from this reference (<a href="https://ipnpr.jpl.nasa.gov/progress_report2/VII/VIIG.PDF">https://ipnpr.jpl.nasa.gov/progress_report2/VII/VIIG.PDF</a>), I tried to construct what we require to estimate the geometric delay only (I'm not bringing in the other delays for the sake of simplicity).</div><div><br></div><div>So, first, we have the raw data, which has to be correlated. Before correlation, we have to estimate the delays that the data streams have to be given. </div><div><br></div><div>For this I need two vectors: the source direction vector & the baseline vector. The baseline vector for both the observing antennas when referenced to the earth's center will be [x1 , y1, z1]  & [x2 , y2 , z2 ] of ITRF system respectively. Let the position of the source in the celestial coordinate be some RA and dec. In the earth centered coordinate system like ITRF, the positions of the antennas are not time variable, but the source vector \hat(s) is time variant (?). </div><div><br></div><div>So for calculating the hour angle of the source, if I'm observing from India and my local sidereal time is say 20 hrs, since my observations are referenced to earth center, the hour angle wrt earth center is 300 degrees (Is this correct ?). </div><div><br></div><div>Now for the obtained <b>b.s </b>vector, I can multiply the earth orientation parameter matrix to factor in the earth rotation effects & add the excess relativistic delays to this. Have I understood the process correctly ? </div><div><br></div><div>Now, what is the ideal rate at which these delays have to be updated ? Is this determined by the fringe rate of the baseline ? Lets say if the fringe rate of my particular baseline is 1 second, then is it sufficient that my source position is updated for every 0.5 s ? Or calculation and correction of these parameters every sampling instant is necessary? </div><div><br></div><div>Thank you,</div><div>Regards,</div><div>Mugundhan</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 25, 2017 at 10:12 AM, Mugundhan vijayaraghavan <span dir="ltr"><<a href="mailto:v.vaishnav151190@gmail.com" target="_blank">v.vaishnav151190@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear David,<div><br></div><div>Thank you so much for the reference. I'll go through it first and come back with more questions (if any) :)</div><div><br></div><div>Thanks,</div><div>Best regards,</div><div>Mugundhan</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 24, 2017 at 8:56 PM, David Gordon <span dir="ltr"><<a href="mailto:geovlbi@gmail.com" target="_blank">geovlbi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div>Hi Mugundhan,<br><br></div>Maybe I can add a little to what Adam has said. I recommend you download the IERS Conventions (2010) technical note<br></div>(<a href="https://www.iers.org/TN36" target="_blank">https://www.iers.org/TN36</a>) for more details. But briefly, the calc programs compute a terrestrial-to-celestial rotation matrix<br></div>at the desired epoch using the value of UT1 at that epoch and smaller rotations for polar motion and precession/nutation.<br></div>It also computes the first and second time derivatives of that rotation matrix. The X,Y, Z site positions, velocities and accelerations<br>are then rotated into the J2000 celestial reference frame. Then some elastic and other corrections are computed and applied at<br>each station - such as the solid earth tide, the ocean loading, the pole tide, short period UT1 and polar motion corrections, etc. <br>An atmosphere delay is also computed. We use the 'concensus' delay model, described in chapter 11 of the conventions.<br></div><div>The source vector used is the un-abberrated direction to the source as seen from the solar system barycenter. Corrections for <br></div><div>aberration and the retarded baseline are contained in the concensus model. The 'geocenter'  is just the (0,0,0) point in the<br></div><div>VLBI reference frame (which we try to align with the ITRF), but does not strictly represent the earth's center of mass. It only <br>needs to be the same point for all the antennas in a correlation. <br><br></div><div>Cheers,<br></div><div>David (<a href="mailto:David.Gordon-1@nasa.gov" target="_blank">David.Gordon-1@nasa.gov</a>)<br></div><div><br></div></div><div class="m_8045491150957767375HOEnZb"><div class="m_8045491150957767375h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 24, 2017 at 3:19 AM, Mugundhan vijayaraghavan <span dir="ltr"><<a href="mailto:v.vaishnav151190@gmail.com" target="_blank">v.vaishnav151190@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Adam,<div><br></div><div>Thank you so much for the clarifications.</div><div><br></div><div>Regards,</div><div><br></div><div>Mugundhan</div></div><div class="m_8045491150957767375m_2834593170064819091HOEnZb"><div class="m_8045491150957767375m_2834593170064819091h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 24, 2017 at 12:44 PM, Adam Deller <span dir="ltr"><<a href="mailto:adeller@astro.swin.edu.au" target="_blank">adeller@astro.swin.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Mugundhan,<div><br></div><div>Since the x axis of the ITRF is nominally aligned with longitude 0 on the Earth, if you want to think of things in terms of hour angles then yes the hour angle will be GMST.  You don't then make any further adjustment for the latitude and longitude of the observing stations.  To be properly precise, though, you will have to take into account precession and nutation as well as polar offset (orientation of the Earth in space c.f. the best model including precession and nutation) and UT1-UTC (Earth's rotational phase c.f. "average").  They are changing the tilt and rotation of the xyz ITRF frame with respect to the J2000 positions, which you can "undo" by offsetting the hour angle and declination of the source very slightly.</div><div><br></div><div>Cheers,</div><div>Adam</div></div><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081HOEnZb"><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081h5"><div class="gmail_extra"><br><div class="gmail_quote">On 24 April 2017 at 14:41, Mugundhan vijayaraghavan <span dir="ltr"><<a href="mailto:v.vaishnav151190@gmail.com" target="_blank">v.vaishnav151190@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Adam,<div><br></div><div>Thank you for your clarifications.</div><div><br></div><div>One more question:</div><div>The unit vector \hat(s) originates from the center of the earth, in this case, the hour angles of the sources will have to be calculated as GMST-RA ? In this case, do we have to offset the obtained hour angles and the declination by the lat and long of the observing stations ? <br></div><div><br></div><div>Thank you,</div><div><br></div><div>Mugundhan</div></div><div class="gmail_extra"><div><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464h5"><br><div class="gmail_quote">On Mon, Apr 24, 2017 at 6:10 AM, Adam Deller <span dir="ltr"><<a href="mailto:adeller@astro.swin.edu.au" target="_blank">adeller@astro.swin.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Mugundhan,<div class="gmail_extra"><br><div class="gmail_quote"><span>On 21 April 2017 at 15:30, Mugundhan vijayaraghavan <span dir="ltr"><<a href="mailto:v.vaishnav151190@gmail.com" target="_blank">v.vaishnav151190@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear All,<div><br></div><div>I have a few queries about how delay modelling is carried out in VLBI for compensating the same.</div><div><br></div><div>1.) The geometric delay is calculated as tg=<b>b.s</b>, where <b>b </b>is the baseline vector and <b>s </b>is the source vector. Lets say I have two antennas, both located about 100 kms apart. How do standard VLBI delay modelling software calculate this delay ? Based on some preliminary reading I understood that the baseline distance are first calculated referenced to the earth center, if this is done, are delays estimated assuming the earth center to be the phase reference ? How is this earth centered reference then transformed to the celestial frame ? because both <b>b </b>and <b>s </b>must be in the same coordinate system for carrying out a dot product operation, right ?</div></div></blockquote><div><br></div></span><div>VLBI delay modeling is very complicated, involving considerably more than just a <b>b.s</b> operation.  Other propagation effects are taken into account too, and the length of the baseline <b>b</b> is changing with time due to tidal forces and what-not, plus the whole system is wobbling around due to the changing earth orientation.</div><div><br></div><div>But stripping it back to the minimum: yes, the Earth centre is usually used as the reference.  Look up the International Terrestrial Reference Frame (ITRF) to see the definition of the axes.  Then you obviously need to know the <b>time</b> (and the Earth orientation parameters) to figure out where the unit vector \hat(s) that points at the direction of the source is pointing in this reference frame.  For each telescope, we then compute the station-based delay from the telescope back to the geocentre at the desired instant of time, and each telescope's data stream is delayed by the computed amount (rather than shifting only one data stream by the difference between \tau_a and \tau_b).  That's what it means to use the geocentre as the reference.<br></div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>2.) In some books/articles i find a reference to a RA and Dec of Baseline ? What does this physically mean ? I'm not able to visualize this clearly. any help will be greatly appreciated ! </div></div></blockquote><div><br></div></span><div>Like I said above, it makes more sense to figure out where the source unit vector is pointing relative to a telescope coordinate system.  You can equivalently rotate the telescope coordinates and keep the source unit vector fixed, but that is (I think) less intuitive.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div>3.) In the complete delay model, tm, which is the sum of geometric delay+clock delay+ionospheric/atmospheric delay+fixed delays due to analog component, the fastest varying component will be geometric delay only, once this is compensated, if the other quantities are contributing to some excess time varying delay, this will be seen as a residual fringe. Now, for clock delay, is this estimated using the allan deviation of the clock being used? Lets say my clock loses 10^-9 seconds in 30 minutes, and if I sample my signal at 16 MHz which is ~ 62.5 ns, will I be able to integrate the data without any degradation due to clock upto 30 minutes ?<br clear="all"><div><br></div></div></div></blockquote><div><br></div></span><div>The sampling time is irrelevant.  It's the sky frequency that determines the visibility phase.  Your signal might only be 16 MHz wide, but if you were observing at 100 GHz then a change of 1 nanosecond translates to 100 turns of phase.  So in that example you could only integrate for maximally a fraction of a second.  Normally VLBI clock drifts are monitored to a level of at worst a few ns/day or so. If they are unknown then a test correlation is performed to determine the clock offset and drift, and then the observation is recorrelated having applied the best available clock model.</div><span><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div></div><div>4.) There is also an associated baseline velocity component which will lead to a time difference between the wavefronts received at both the antennas. Is this baseline velocity the same as the orbital velocity of the earth ? Or is this modelled differently ?</div><div><br></div></div></div></blockquote><div><br></div></span><div>By delay tracking to the geocentre, this problem is naturally taken into account. When you use the geocentre, you are automatically forced to account for the rotation of the reference frame between the time the signal is received at the antenna and the time that it would pass through the geocentre.  So you've corrected for the velocity of both of the stations, rather than their difference.  The process is known as retarded baseline correction.<br><br></div><div>Unfortunately the documentation for VLBI delay packages is not extensive.  You can look up CALC (<a href="https://lupus.gsfc.nasa.gov/software_calc_solve.htm" target="_blank">https://lupus.gsfc.nasa.gov/s<wbr>oftware_calc_solve.htm</a>) or VTD (<a href="http://astrogeo.org/vtd/" target="_blank">http://astrogeo.org/vtd/</a>) but neither have an excellent explanation of the theory.<br></div><div><br></div><div>Cheers,<br></div><div>Adam<br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span><div dir="ltr"><div><div></div><div>I would greatly appreciate if the experts here clarify my doubts. Kindly do point me to references that may lead to clarification of these doubts too !</div><div><br></div><div>Thanking you,</div><div>With best regards,</div><div><br></div><div>Mugundhan V.</div><br>
</div></div>
<br></span>______________________________<wbr>_________________<br>
Difx-users mailing list<br>
<a href="mailto:Difx-users@listmgr.nrao.edu" target="_blank">Difx-users@listmgr.nrao.edu</a><br>
<a href="https://listmgr.nrao.edu/mailman/listinfo/difx-users" rel="noreferrer" target="_blank">https://listmgr.nrao.edu/mailm<wbr>an/listinfo/difx-users</a><br>
<br></blockquote></div><span class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464m_-7958279512305540176HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464m_-7958279512305540176m_4722421876091630126gmail-m_-1367036963978668279gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px">!=============================<wbr>==============================<wbr>==!<br>Dr. Adam Deller         </div><div dir="ltr" style="font-size:12.8px">ARC Future Fellow, Senior Lecturer</div><div style="font-size:12.8px">Centre for Astrophysics & Supercomputing </div><div dir="ltr" style="font-size:12.8px">Swinburne University of Technology    <br>John St, Hawthorn VIC 3122 Australia</div><div style="font-size:12.8px">phone: <a href="tel:+61%203%209214%205307" value="+61392145307" target="_blank">+61 3 9214 5307</a></div><div style="font-size:12.8px">fax: <a href="tel:+61%203%209214%208797" value="+61392148797" target="_blank">+61 3 9214 8797</a></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">office days (usually): Mon-Thu<br>!=============================<wbr>==============================<wbr>==!</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div></div></div><span class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464HOEnZb"><font color="#888888">-- <br><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464m_-7958279512305540176gmail_signature" data-smartmail="gmail_signature">the giver of moksha</div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081m_-3578998529350541464gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr" style="font-size:12.8000001907349px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px"><div dir="ltr" style="font-size:12.8px">!=============================<wbr>==============================<wbr>==!<br>Dr. Adam Deller         </div><div dir="ltr" style="font-size:12.8px">ARC Future Fellow, Senior Lecturer</div><div style="font-size:12.8px">Centre for Astrophysics & Supercomputing </div><div dir="ltr" style="font-size:12.8px">Swinburne University of Technology    <br>John St, Hawthorn VIC 3122 Australia</div><div style="font-size:12.8px">phone: <a href="tel:+61%203%209214%205307" value="+61392145307" target="_blank">+61 3 9214 5307</a></div><div style="font-size:12.8px">fax: <a href="tel:+61%203%209214%208797" value="+61392148797" target="_blank">+61 3 9214 8797</a></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">office days (usually): Mon-Thu<br>!=============================<wbr>==============================<wbr>==!</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_8045491150957767375m_2834593170064819091m_-2927150448516471081gmail_signature" data-smartmail="gmail_signature">the giver of moksha</div>
</div>
</div></div><br>______________________________<wbr>_________________<br>
Difx-users mailing list<br>
<a href="mailto:Difx-users@listmgr.nrao.edu" target="_blank">Difx-users@listmgr.nrao.edu</a><br>
<a href="https://listmgr.nrao.edu/mailman/listinfo/difx-users" rel="noreferrer" target="_blank">https://listmgr.nrao.edu/mailm<wbr>an/listinfo/difx-users</a><br>
<br></blockquote></div><br></div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="m_8045491150957767375gmail_signature" data-smartmail="gmail_signature">the giver of moksha</div>
</div>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">the giver of moksha</div>
</div>