Differences
This shows you the differences between two versions of the page.
projects:mkat_sband:pub:fringefit [2020/11/05 16:32] wucknitz |
projects:mkat_sband:pub:fringefit [2020/11/10 14:56] (current) wucknitz |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== Fringe-fitting ====== | ====== Fringe-fitting ====== | ||
+ | |||
+ | ===== Procedure ===== | ||
Frequency ranges strongly affected by RFI are flagged (not included) before anything else. These ranges are chosen manually. | Frequency ranges strongly affected by RFI are flagged (not included) before anything else. These ranges are chosen manually. | ||
Line 28: | Line 30: | ||
Besides explicitly fitting the curvature to distinguish between dispersive and non-dispersive delays, we can thus combine the group delay with the phase to achieve the same result, actually with higher accuracy. Note that we can have arbitrary constants per baseline in each of the parameters (with the exception of phase), because they can be partly absorbed into the bandpass. Before applying the equations one should thus generally subtract the mean. | Besides explicitly fitting the curvature to distinguish between dispersive and non-dispersive delays, we can thus combine the group delay with the phase to achieve the same result, actually with higher accuracy. Note that we can have arbitrary constants per baseline in each of the parameters (with the exception of phase), because they can be partly absorbed into the bandpass. Before applying the equations one should thus generally subtract the mean. | ||
+ | |||
+ | |||
+ | ===== Test for dataset 1603186576 ===== | ||
+ | |||
As a test for this concept we use the dataset (see [[projects:mkat_sband:pub:imobs_list|the list]]) 1603186576, which is one hour on PKS 1934-638 in 4k mode in band S3. We fringe-fit over 60-sec blocks with curvature and delay rate. As example we use the baseline m045-m060, because it has the strongest dispersive delay variations. | As a test for this concept we use the dataset (see [[projects:mkat_sband:pub:imobs_list|the list]]) 1603186576, which is one hour on PKS 1934-638 in 4k mode in band S3. We fringe-fit over 60-sec blocks with curvature and delay rate. As example we use the baseline m045-m060, because it has the strongest dispersive delay variations. | ||
Line 51: | Line 57: | ||
- | We absorb time-averaged closure errors into the baseline-based bandpasses before trying to separate them into station-based functions. See [[projects:mkat_sband:pub:bandpass|Bandpass part.]] It turns out that remaining closure errors (as small as they are) can be minimised even more by excluding short baselines. For PKS 1934-63 a baseline limit (measured in antenna XY coordinates in order to be time-independent, not in UV) of 800 metres works well and keeps a sufficient number of baselines for the fits. | + | We absorb time-averaged closure errors into the baseline-based bandpasses before trying to separate them into station-based functions. See [[projects:mkat_sband:pub:bandpass|Bandpass part.]] It turns out that remaining closure errors (as small as they are) can be minimised even more by excluding short baselines. For PKS 1934-63 a baseline limit (measured in antenna XY coordinates or UV) of a few hundred metres works well and keeps a sufficient number of baselines for the fits. |
+ | |||
+ | |||
+ | |||
+ | ===== Interpretation of solutions, phase stability ===== | ||
+ | |||
+ | As basis for this we use the station-based solutions, because of their better SNR. Whenever needed, we can re-derive the baseline-based values from them. | ||
+ | |||
+ | All effects are largest on the longest baselines and thus on the outer stations. Here we show solutions for non-dispersive and dispersive delays for four stations representing the corners of the array. Both polarisations are combined, because they are consistent with each other. | ||
+ | |||
+ | non-dispersive delays (troposphere+instrumental): \\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_corners_nondisp.png?direct&500|non-dispersive delays}} | ||
+ | |||
+ | dispersive delays (ionosphere):\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_corners_disp.png?direct&500|dispersive delays}} | ||
+ | |||
+ | In the dispersive delays we notice that the large fluctuations in the first half are very similar but with opposite signs in the left and right panel. The signature of the second half has opposite signs in the top and bottom plots and is weaker in the bottom. We interpret this as follows: The ionosphere causes mostly large-scale gradients over the entire array. In the first half the direction of this gradient is east-west, in the second half it is more north-south. This is consistent with a large scale of ionospheric variations and with a high speed at which they cross the array. This means we do not see individual effects at each station, but are really dominated by one gradient at each time. | ||
+ | |||
+ | We do not see a similar behaviour for the non-dispersive delays. Their fluctuations are much faster and we know that they cannot cross the array faster than the wind speed. This means that we really see individual signatures and not just a large-scale gradient. | ||
+ | |||
+ | We can test the gradient-hypothesis by fitting such a gradient (independent for each time) to the data and check how well it can reproduce the observed signal. We can do this per station or per baseline. The former has the advantage that we have to check less plots. | ||
+ | |||
+ | Here are results for the antenna-based dispersive delay:\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_ant_disp.png?direct&600|gradient fit for antenna-based dispersive delay}} | ||
+ | |||
+ | We see that the gradient fit (green) fits the observations (red) almost exactly for all stations. For the innermost stations (left) the noise dominates the signal, but the systematic variations still agree very well. This means the ionosphere really acts mostly as a large-scale time-varying gradient. | ||
+ | |||
+ | Here is the result for the non-dispersive delay:\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_ant_nondisp.png?direct&600|gradient fit for antenna-based non-dispersive delay}} | ||
+ | |||
+ | For the outer stations (right) the agreement is still pretty good, but this is because we have only a few outer stations, which thus have a strong lever arm and dominate the solution. In contrast to the dispersive delay, this solutions does not fit the inner stations (left) well. Fluctuations there are dominated by small-scale and more local variations. | ||
+ | |||
+ | |||
+ | Here are PDF files of all these solutions for completeness: | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_ant_disp.pdf|gradient fit for antenna-based dispersive delay}}\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_ant_nondisp.pdf|gradient fit for antenna-based non-dispersive delay}}\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_bl_disp.pdf|gradient fit for baseline-based dispersive delay}}\\ | ||
+ | {{:projects:mkat_sband:pub:1603186576_linxy_bl_nondisp.pdf|gradient fit for baseline-based non-dispersive delay}} |