REVIEW 1 major objections 34 references
Station coordinate errors substantially reduce pulsar timescale construction accuracy when zenith angles vary long-term.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.3
2026-06-27 23:39 UTC pith:XWKCSQ3Z
load-bearing objection The paper runs straightforward TEMPO2 simulations of station coordinate errors but misinterprets its own Kendall correlation (r=1.67%, p=100%) as support for Roemer dominance when the numbers show no relationship. the 1 major comments →
The Analysis of the Influence of Coordinate Error of Observation Station On the Construction Accuracy of Pulsar Time
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Errors in observatory coordinates directly impact the precision of pulsar time-scale construction. Using TEMPO2 simulations of various station position errors for three millisecond pulsars over 13 days and 5 years, the analysis shows that station coordinate errors substantially reduce the accuracy of pulsar timescale construction when the zenith angle exhibits long-term variations. This holds independent of pulsar type and daily observable time. A linear relationship exists between station coordinate errors and the RMS of pulsar timing residuals, with the Roemer delay error caused by coordinate inaccuracies notably larger than other terms.
What carries the argument
Roemer delay error induced by inaccuracies in the three-dimensional terrestrial reference frame coordinates of the observation station.
Load-bearing premise
Long-term variations in zenith angle are present during the simulated observations.
What would settle it
Real observations with varying zenith angles that show no corresponding linear increase in RMS timing residuals with station coordinate errors would falsify the central relationship.
If this is right
- Station coordinate errors produce a linear effect on RMS of pulsar timing residuals with fitted coefficients from 1.36×10^{-11} to 1.61×10^{-9}.
- Errors along x- and y-axes have comparable influence on timing precision while z-axis errors have smaller effect.
- The degradation is independent of pulsar type and the daily observable time of the station antenna.
- At current timing precision, coordinate errors affect pulse arrival times mainly through the Roemer delay term.
- The reported effects may not apply under constant zenith angle or limited elevation angles.
Where Pith is reading between the lines
- Observatories located where zenith angles stay nearly constant could experience reduced sensitivity to coordinate errors.
- Routine high-precision station positioning or real-time corrections might improve long-term stability of pulsar-based timescales.
- The same coordinate-error mechanism could be tested in other radio timing applications that rely on geometric delays.
- Extending the simulations to include actual multi-year data sets with measured elevation variations would provide a direct check.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper simulates the effects of observatory station coordinate errors (in the ITRF) on pulsar timing residuals using TEMPO2 for three millisecond pulsars over 13-day and 5-year spans. It claims that such errors substantially degrade pulsar timescale construction accuracy whenever zenith angle exhibits long-term variations (independent of pulsar type), that a linear relationship exists between coordinate error magnitude and RMS residual with fitted slopes 1.36×10^{-11} to 1.61×10^{-9}, that the induced Roemer delay error dominates other delay/correction terms, and that x/y-axis errors affect timing more than z-axis errors. The sole quantitative support cited for Roemer dominance is a Kendall rank correlation between Roemer delay error and RMS yielding r=1.67% and p=100% in all cases, interpreted as confirming that coordinate errors act primarily through the Roemer term and are consistent with theory. The authors note the findings may not apply under constant zenith angle or limited-elevation conditions such as FAST.
Significance. If the simulation results and their interpretation were robust, the work would provide a concrete error-budget contribution for pulsar timing arrays and pulsar-based timescales, quantifying how station-position uncertainty propagates via Roemer delay and supplying linear coefficients that could be used in observation planning. The axis-dependent and zenith-variation dependence would also be useful for site selection and scheduling.
major comments (1)
- [Kendall correlation analysis (abstract and results)] Kendall correlation analysis (abstract and results section): the reported coefficient r = 1.67 % with p = 100 % is statistically indistinguishable from zero correlation and is interpreted in the text as evidence that “coordinate errors primarily affect the Roemer delay term … which is highly consistent with theoretical models.” A near-zero r with p-value = 1 indicates absence of monotonic association, directly contradicting the claim that Roemer delay is “notably larger than other delay and correction terms” and that the correlation supports the dominance conclusion. This internal tension is load-bearing for the central quantitative claim.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review. The single major comment identifies a clear inconsistency in our use of the Kendall correlation statistic. We address it directly below and agree that revision is required.
read point-by-point responses
-
Referee: Kendall correlation analysis (abstract and results section): the reported coefficient r = 1.67 % with p = 100 % is statistically indistinguishable from zero correlation and is interpreted in the text as evidence that “coordinate errors primarily affect the Roemer delay term … which is highly consistent with theoretical models.” A near-zero r with p-value = 1 indicates absence of monotonic association, directly contradicting the claim that Roemer delay is “notably larger than other delay and correction terms” and that the correlation supports the dominance conclusion. This internal tension is load-bearing for the central quantitative claim.
Authors: We agree with the referee that the reported Kendall tau of 1.67 % (p = 100 %) indicates no monotonic association and cannot support the stated interpretation. This is an error in our statistical analysis and its textual framing. We will remove all references to the Kendall correlation from the abstract and results. The claim that Roemer delay error is notably larger will be retained only where it is directly supported by the magnitude comparisons performed in the TEMPO2 simulations; the correlation statistic will no longer be invoked. The linear RMS–coordinate-error relations and axis-dependent findings are unaffected by this change. revision: yes
Circularity Check
No significant circularity; results derive from forward simulation
full rationale
The paper conducts forward simulations in the established public TEMPO2 package, computes RMS residuals and Roemer delays from those runs, then reports fitted linear coefficients and a Kendall correlation directly from the simulation outputs. No claimed prediction reduces to its own inputs by construction, no self-citation is used as a load-bearing premise, and the central claims (linear error-RMS relation, Roemer dominance) are statistical summaries of the simulated data rather than self-referential derivations. The work is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (1)
- linear coefficients relating coordinate error to RMS residual
axioms (1)
- domain assumption TEMPO2 accurately models all relevant delay terms including the Roemer delay for perturbed station coordinates
read the original abstract
\abstract{Errors in observatory coordinates directly impact the precision of pulsar time-scale construction. Using the pulsar timing software TEMPO2, this study simulates various station position errors within the three-dimensional terrestrial reference frame for three different types of millisecond pulsars, over periods of 13 days and 5 years, and analyzes their effects on pulsar timing results.The findings demonstrate that,for both 13-day and 5-year observation spans, station coordinate errors substantially reduce the accuracy of pulsar timescale construction when the zenith angle exhibits long-term variations. This effect is independent of pulsar type and the daily observable time of the station antenna for the pulsar. A linear relationship is found between station coordinate errors and the Root-Mean-Square (RMS) of pulsar timing residuals, with fitted linear coefficients ranging from $1.36 \times 10^{-11}$ to $1.61 \times 10^{-9}$ for the three pulsars. The Roemer delay error caused by coordinate inaccuracies is notably larger than other delay and correction terms. Errors along the x- and y-axes have comparable influences on timing precision, whereas errors along the z-axis have a relatively smaller effect. Kendall correlation analysis between station error-induced Roemer delay and RMS yields a correlation coefficient $r = 1.67\%$ and $p = 100\%$ in all cases, indicating that, at current timing precision levels, coordinate errors primarily affect the Roemer delay term and thus the pulse arrival times, which is highly consistent with theoretical models.While these findings offer valuable insights into the key factors influencing pulsar timescale accuracy and related applications, they may not hold under conditions of a constant zenith angle or limited elevation angles, such as those at FAST.}
Figures
Reference graph
Works this paper leans on
-
[1]
Taylor, J. H. 1991. Millisecond pulsars: nature’s most stable clocks. IEEE Proceedings, 79
1991
-
[2]
N., et al
Hobbs, G., Guo, L., Caballero, R. N., et al. 2020. A pulsar-based time-scale from the International Pulsar Timing Array. Monthly Notices of the Royal Astronomical Society, 491(4), 5951–5965
2020
-
[3]
Zhao, C., Chen, D., Cai, H., et al. 2011. Research on Available Target Sources for X-ray Pulsar Navigation. Progress in Astronomy , 29(3), 334–
2011
-
[4]
Shuai, P., Liu, Q., Huang, L., et al. 2019. The first pulsar navigation test satellite and its ob- servational results. Journal of Chinese Inertial Technology, 27(3), 281–287. (in Chinese)
2019
-
[5]
J., et al
Hobbs, G., Jenet, F., Lee, K. J., et al. 2009. TEMPO2: a new pulsar timing package - III. Gravitational wave simulation. Monthly Notices of the Royal Astronomical Society, 394(4), 1945–1955
2009
-
[6]
Yang, T. 2008. Measurement of X-ray Pulsar Pulse Time of Arrival at Spacecraft. Chinese Journal of Space Science, (4), 330–334. (in Chinese)
2008
-
[7]
Zhou, Q., Wei, Z., Yan, L., et al. 2021. Research on Space-Ground Integrated Pulsar Time for Comprehensive PNT System. Acta Physica Sinica, 70(13), 471–483. (in Chinese)
2021
-
[8]
Han, M. 2023. Research on Frequency Deviation Estimation Method of Atomic Clock Based on Pulsar. Ph.D. thesis, University of Chinese Academy of Sciences (National Time Service Center, Chinese Academy of Sciences). DOI: 10.27547/d.cnki.gkgsc.2023.000005. (in Chinese)
-
[9]
Ding, Y., Tong, M., Zhao, C., et al. 2017. Analysis of Factors Affecting the Stability of Pulsar Time. Journal of Time and Frequency, 40(4), 260–267. (in Chinese)
2017
-
[10]
Zhu, X., Fu, Y., Cai, F., et al. 2020. Research on the Influence Mechanism of Nonlinear Variation of GNSS Station Coordinates. Progress in Geophysics, 35(1), 79–85. (in Chinese)
2020
-
[11]
Liu, Y., Wang, H., He, L., et al. 2025. The impact of yaw attitude models on precise orbit determina- tion: The latest blocks of GNSS satellites and their yaw models. Advances in Space Research , 75(1), 1310–1329
2025
-
[12]
Guillory, J., Truong, D., Wallerand, J.-P., L¨ osler, M., Eschelbach, C., M¨ ahler, S., Kl¨ ugel, T. 2023. Determination of the reference point of a radio telescope using a multilateration-based coordinate measurement prototype. Precision Engineering, 83, 69–81
2023
-
[13]
Ma, X., Zhang, Z., Sun, Z., et al. 2023. Determination of the Reference Point of Radio Telescope by Unattended GNSS Method. Acta Astronomica Sinica, 64(2), 1–10. (in Chinese)
2023
-
[14]
Tian, F. 2011. Autonomous Positioning Error Estimation Method for Spacecraft Based on Pulsar Timing. Progress in Astronomy, 29(1), 97–104. (in Chinese)
2011
-
[15]
I., Ray, P
Sheikh, S. I., Ray, P. S., Wolff, M. T., et al
-
[16]
In: Proceedings of the 63rd Annual Meeting of The Institute of Navigation
Relative Navigation of Spacecraft Utilizing Bright, Aperiodic Celestial Sources. In: Proceedings of the 63rd Annual Meeting of The Institute of Navigation
-
[17]
A., Speyer, J
Emadzadeh, A. A., Speyer, J. L., et al. 2011. Relative Navigation Between Two Spacecraft Using X-ray Pulsars. IEEE Transactions on Control Systems Technology, 19(5), 1021–1035
2011
-
[18]
Li, L., Wang, G., Guo, L., et al. 2018. Astrometric Considerations for Pulsar Navigation. Journal of Deep Space Exploration , 5(3), 235–240. (in Chinese)
2018
-
[19]
Kaur, D., Hobbs, G., Zic, A., et al. 2025. Unlocking the hidden potential of pulsar astronomy. New Astronomy, 121, 102460
2025
-
[20]
Tong, M., Yang, T., Zhao, C., et al. 2017. Analysis and Estimation of Measurement Accuracy of Pulsar Timing Model Parameters. Scientia Sinica: Physica, Mechanica & Astronomica , 47(9), 103–
2017
-
[21]
T., Hobbs, G
Edwards, R. T., Hobbs, G. B., Manchester, R. N., et al. 2006. TEMPO2, a new pulsar timing package - II. The timing model and precision estimates. Monthly Notices of the Royal Astronomical Society , 372(4), 1549–1574
2006
-
[22]
Han, W., Wang, N., Wang, J., Yuan, J., He, D
-
[23]
Astrophysics and Space Science, 364(3), 48
Using single millisecond pulsar for terrestrial position determination. Astrophysics and Space Science, 364(3), 48
-
[24]
Luzum, B., Capitaine, N., Fienga, A., et al. 2011. The IAU 2009 system of astronomical constants: the report of the IAU working group on numerical standards for Fundamental Astronomy. Celestial Mechanics and Dynamical Astronomy, 110(4), 293– 304
2011
-
[25]
J., Kapur, A., et al
Zic, A., Reardon, D. J., Kapur, A., et al. 2023. The Parkes Pulsar Timing Array third data re- lease. Publications of the Astronomical Society of 572 www.ati.ac.cn Australia, 40, e049
2023
-
[26]
Perera, B. B. P., DeCesar, M. E., Demorest, P. B., et al. 2019. The International Pulsar Timing Array: second data release. Monthly Notices of the Royal Astronomical Society, 490(4), 4666–4687
2019
-
[27]
Nice, D., Demorest, P., Stairs, I., et al. 2015. Tempo: Pulsar timing data analysis. Astrophysics Source Code Library, ascl:1509.002
2015
-
[28]
V., Kopejkin, S
Doroshenko, O. V., Kopejkin, S. M., et al. 1990. High-precision timing measurement of single pul- sars. Astronomicheskii Zhurnal, 67, 986
1990
-
[29]
B., Edwards, R
Hobbs, G. B., Edwards, R. T., et al. 2006. TEMPO2, a new pulsar-timing package - I. An overview. Monthly Notices of the Royal Astronomical Society, 369(2), 655–672
2006
-
[30]
Hobbs, G., Archibald, A., Arzoumanian, Z., et al. 2010. The International Pulsar Timing Array project: using pulsars as a gravitational wave detec- tor. Classical and Quantum Gravity , 27(8), 084013
2010
-
[31]
EPTA Collaboration, InPTA Collaboration, Antoniadis, J., et al. 2023. The second data release from the European Pulsar Timing Array. III. Search for gravitational wave signals. Astronomy & Astrophysics, 678, A50
2023
-
[32]
J., Zic, A., Shannon, R
Reardon, D. J., Zic, A., Shannon, R. M., et al
-
[33]
The Astrophysical Journal Letters , 951(1), L6
Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array. The Astrophysical Journal Letters , 951(1), L6
-
[34]
Xu, H., Chen, S., Guo, Y., et al. 2023. Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I. Research in Astronomy and Astrophysics, 23(7), 075024. Astronomical Techniques and Instruments, 1(5), 560–573, 2025 573
2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.