Pith. sign in

REVIEW 2 major objections 4 minor 37 references

Probing Variations in Earth's Ionosphere Using Pulsars

T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper shows that 10-minute rotation measure monitoring of one pulsar can reconstruct daily ionospheric electron content variations, including noontime bite-outs, in agreement with global maps.

desk verdict Useful high-cadence pulsar RM dataset, but the TEC reconstruction has a partially circular zero-point that weakens the quantitative agreement and the bite-out detections. read the letter →

arxiv 2608.06744 v1 pith:SWYPFYPL submitted 2026-08-07 astro-ph.HE astro-ph.EPastro-ph.IM

Yanqing Cai , Ziwei Wu , Weiwei Zhu , Joris P. W. Verbiest , Ningbo Wang , Moochickal Ambalappat Krishnakumar , Lars Künkel , Jörn Künsemöller
show 19 more authors
Yulan Liu Nataliya Porayko Golam M. Shaifullah Caterina Tiburzi Marcus Brüggen Ralf-Jürgen Dettmar Ziyao Fang Qiuyang Fu Jiawei Jin Caisong Liu Lingqi Meng Xueli Miao Jiarui Niu Dominik J. Schwarz Christian Vocks Mengyao Xue Mao Yuan Youling Yue Zhen Zhang
This is my paper · ORCID
classification astro-ph.HEastro-ph.EPastro-ph.IM
keywords pulsarsionosphererotationmeasureFaradaytotalelectroncontentnoontimebite-outlow-frequencyradioastronomyspaceweather
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to show that high-cadence rotation measure (RM) monitoring of a single pulsar can act as a probe of the Earth's ionosphere. Using 32 days of 10-minute observations of PSR J0814+7429 at low radio frequencies, it separates the ionosphere's contribution from the interstellar one and reconstructs the line-of-sight total electron content (TEC_LoS) with the World Magnetic Model. The reconstructed TEC_LoS follows the expected daily rise and fall and reproduces the noontime bite-out on two days, matching independent global ionosphere maps to within about 1–2 TECU. If the interpretation holds, pulsars would offer a complementary, all-sky way to monitor ionospheric variability, including over oceans and high latitudes where GNSS coverage is thin.

What carries the argument

The load-bearing machinery is the conversion from observed Faraday rotation to line-of-sight electron content. Equation 2, $\phi_{\rm ion} = 2.6\times10^{-17}\, \mathrm{TEC}_{\rm LoS}\, B_{\rm LoS}$, relates the ionospheric rotation measure to the product of TEC_LoS and the geomagnetic field $B_{\rm LoS}$ along the line of sight, and Equation 3 inverts it: $\mathrm{TEC}_{\rm LoS} = 1\times10^{17}\, \phi_{\rm obs-ISM} / (2.6\, B_{\rm LoS})$. This inversion assumes the ionosphere is a thin spherical shell at a fixed effective height, that $B_{\rm LoS}$ is evaluated at that height using the World Magnetic Model, and that the interstellar RM is constant over the campaign. The same shell model also underlies the independent GIM-based estimates, so the comparison tests the consistency of the reconstruction.

What would settle it

Compare the same pulsar's reconstructed TEC_LoS against an independent co-located GNSS receiver's vertical TEC converted to the same line of sight. If the two agree within 1–2 TECU over many days, the shell assumption holds; a persistent deviation that grows when the pulsar is at low elevation or during a geomagnetic storm would falsify the thin-shell conversion.

Watch

Extended reading notes

Core claim

The central claim is that the observed 0.5–1 rad $m^{{-2}}$ diurnal variation in the RM of PSR J0814+7429 is dominated by the Earth's ionosphere, not by the interstellar medium or any companion. After subtracting a constant interstellar RM of −13.60 ± 0.06 rad $m^{{-2}}$, the paper converts the residual rotation measure into line-of-sight ionospheric electron content using a thin-shell model and the World Magnetic Model. The resulting TEC_LoS curves agree with independent GNSS-based global ionosphere maps, reproduce the solar-driven shape of the diurnal cycle, and on two of the 32 days show a clear noontime bite-out with a valley near 14:00 local time and an average duration of about 8 hours. The paper concludes that pulsar RM monitoring is a feasible probe of temporal variations in ionospheric electron density.

Load-bearing premise

The reconstruction depends on treating the ionosphere as a thin spherical shell at a single fixed height, with the magnetic field evaluated at that height; if the real electron layer is thicker, lower, or shifted, the RM-to-TEC conversion is systematically biased.

Editorial extensions

If this is right

  • If correct, the method means a single pulsar's polarized signal can deliver ionospheric TEC along that line of sight at 10-minute cadence, without any satellite.
  • The 0.5–1 rad m^{-2} daily RM swing sets the scale of the ionospheric correction needed for precision pulsar timing and absolute polarimetry at low frequencies.
  • Observing more pulsars in different directions would extend ionospheric monitoring to regions where GNSS receivers are sparse, such as oceans and high latitudes.
  • The detection of two noontime bite-outs in 32 days suggests that pulsar RM time series can capture small-scale, solar-driven ionospheric structures that daily models may smooth over.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural extension the paper does not pursue is to check whether the reconstructed TEC_LoS is independent of the pulsar's elevation angle; if a fixed shell height misrepresents the true profile, a systematic elevation-dependent bias should appear.
  • Multi-station simultaneous observations, like the two shown in Figure 3, could be turned into a triangulation of small-scale ionospheric structures if more than two stations observe the same pulsar at once.
  • The same RM-to-TEC conversion could be applied to fast radio bursts or other polarized transients to obtain instantaneous ionospheric probes along arbitrary lines of sight, though burst RM includes the host galaxy and intergalactic contributions.
  • A practical testable extension would be to compare pulsar-derived TEC_LoS against co-located GNSS TEC during a known geomagnetic storm, where the thin-shell assumption is most likely to break down.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper presents 32 days of LOFAR/GLOW monitoring of PSR J0814+7429 with 10-minute time resolution, showing diurnal rotation-measure (RM) variations of roughly 0.5-1 rad/m2. It attributes these variations to the Earth's ionosphere, based on the amplitude being four to five orders of magnitude above expected interstellar turbulence and on the pulsar being isolated. The authors subtract a constant interstellar RM, calibrated using the CODE global ionosphere map, to isolate the ionospheric RM, and then use a thin-shell model with the World Magnetic Model to convert this RM into line-of-sight total electron content (TEC_LoS, Eq. 3). They report good agreement with TEC_LoS from several independent GIM products, identify two days with noontime bite-outs, and use two-station comparisons to suggest a solar-activity effect on July 9, 2017.

Significance. If the quantitative reconstruction issues can be resolved, the paper would demonstrate a new, complementary use of existing pulsar RM monitoring for ionospheric science. The strengths include the high-cadence, nearly full-day sampling, the multi-station dataset, the explicit amplitude argument separating ionospheric from interstellar contributions, and the use of public software and GIM products. The claim is falsifiable through the stated bite-out criterion. However, the quantitative TEC_LoS results are currently entangled with the CODE-based zero-point calibration and an unspecified shell height, so the significance is conditional on those points being addressed. The paper is timely given the growth of low-frequency arrays and the need for ionospheric monitoring over oceans and high latitudes.

major comments (2)
  1. [3.2, Eqs. (2)-(3), Fig. 2, Sec. 3.3]
  2. [3.2, Eqs. (2)-(3)]
minor comments (4)
  1. [3.1, Fig. 1 caption]
  2. [2, data processing]
  3. [3.4]
  4. [3.2, Fig. 2]

Circularity Check

1 steps flagged · score 6.0 of 10

Pulsar TEC zero-point is anchored to the CODE GIM used to set the ISM offset, so the later 'independent' GIM agreement and the noontime bite-out detections are partly a restatement of that same model.

  1. fitted input called prediction [Section 3.2, 'Probing the Ionospheric Electron Density', Eqs. (2)-(3)]
    "We first compute the daily mean residuals (RMobs − RMion) between RMobs and RMion to estimate the LoS ISM RM for each day. ... We finally take the mean of these values to obtain the ISM RM, yielding a value of −13.60±0.06 rad m−2. ... the ionospheric RM variations along the LoS, hereafter denoted as φobs−ISM, are derived by subtracting the final ISM RM from the RMobs. ... By equating φion with φobs−ISM in Equation 2, the value of TECLoS can be derived as follows: TECLoS = 1×10^17 φobs−ISM / (2.6 × BLoS) m−2."

    RM_ISM is defined as mean(RM_obs − RM_ion_CODE), so φobs−ISM = RM_obs − RM_ISM has the same daily mean as RM_ion_CODE by construction. Plugging into Eq. (3) gives TEC_PSR(t) = T_true(t) + [mean(T_CODE B) − mean(T_true B)]/B(t); the CODE-GIM offset is therefore imprinted on every TEC point and modulated by the diurnal B_LoS(t). The 'good agreement' with GIM TEC in Figure 2 is partly the reappearance of the same model used to set the zero-point, and the offset can produce or enhance the 1/B-shaped bite-out signal. Independent content is limited to the diurnal shape of RM_obs, not the absolute TEC level or the GIM comparison.

full rationale

The paper contains one genuine partial circularity: the zero-point of the pulsar-derived TEC is fitted to the CODE GIM via the ISM-RM estimate, and that same model family (GIM products) is then used as the independent comparison. Because BLoS varies over the day, the constant CODE-vs-true offset is converted into a diurnal 1/B term that can mimic or enhance a noontime bite-out at the reported 1-2 TECU level. This warrants a score of 6 under the 'partial circularity' criterion. The remaining derivations are not circular: the thin-shell/Heff approximation and WMM field are external inputs with stated uncertainties; the ISM-turbulence bound cited from Porayko et al. (2019) is an external peer-reviewed measurement, not a result of this paper; and the 10-minute RM time series themselves are independent observables. The paper's limitation to one LoS and the unstated Heff are correctness risks, not circularity.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central reconstruction rests on a small set of external models and fitted quantities: a constant interstellar RM estimated from the data, an implicit thin-shell height, and the CODE GIM/WMM used to calibrate the ionospheric contribution. No new physical entities are introduced. The most important caveat is that the zero-point of TEC_LoS is set by an ionosphere model, so absolute values are not fully independent, although the diurnal variations are data-driven.

free parameters (2)
  • Interstellar RM (RM_ISM) = -13.60 +/- 0.06 rad/m^2
    Estimated as the mean over 32 days of daily (RM_obs minus CODE-GIM RM_ion) residuals; used to isolate the ionospheric contribution in Eq. 3, so the zero-point of the derived TEC depends on this fit and on the ionosphere model.
  • Effective ionospheric shell height H_eff = Not stated (implicit RMEXTRACT default)
    Required to evaluate B_LoS from WMM and to map RM to TEC_LoS through the thin-shell approximation; never given in the paper.
assumptions (4)
  • domain assumption The ionosphere is a thin spherical shell at a fixed effective height H_eff.
    Invoked in Section 3.2, Eq. 2 and Eq. 3; allows B_LoS to be evaluated at a single point. If the electron distribution is not shell-like, the TEC_LoS values are biased.
  • domain assumption The interstellar RM is constant over the 32-day campaign.
    Section 3.2 takes the mean of daily residuals as a constant RM_ISM after a Spearman trend test. No intraday ISM variability is considered, justified only by the pulsar being isolated and by ISM turbulence estimates.
  • domain assumption CODE GIM and WMM correctly describe the ionosphere and geomagnetic field for calibration purposes.
    Section 2 uses RMEXTRACT with CODE and WMM to compute RM_ion, which sets the RM_ISM zero-point. Model errors become systematic offsets in the derived TEC.
  • domain assumption Faraday contributions from the ISM and ionosphere add linearly, with no other significant contributors.
    Used in Eq. 1 and Section 3.1 to separate ionospheric from interstellar RM. The isolated pulsar makes a companion contribution unlikely, but source-intrinsic RM variability is not discussed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Probing Variations in Earth's Ionosphere Using Pulsars." pith.science (2026). https://pith.science/paper/SWYPFYPL

@misc{pith2026260806744,
  author       = {Pith},
  title        = {Pith review of: Probing Variations in Earth's Ionosphere Using Pulsars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SWYPFYPL}},
  note         = {Machine review of arXiv:2608.06744}
}
read the original abstract

We present high cadence 10-minute rotation measure (RM) monitoring of PSR~J0814+7429 using the LOw-Frequency ARray, aiming to probe ionospheric variability along the pulsar line of sight (LoS). By separating the ionospheric contribution from the observed RM, we quantitatively reconstruct the diurnal variation of the ionospheric electron density along the pulsar LoS based on the World Magnetic Model. The derived variations exhibit clear solar-driven modulation, including the ionospheric noontime bite-out phenomenon, and show good agreement with the LoS total electron content reconstructed from independent global vertical total electron content maps. These results demonstrate the feasibility of using pulsars as probes of temporal variations in the electron density of the Earth's ionosphere.

Figures

Figures reproduced from arXiv: 2608.06744 by the authors.

Figure 1
Figure 1. The time series of RMs as a function of local time during a day for nine observations made in 2017. For the top panel in each plot, the red points show the RM of PSR J0814+7429 (left axis label, RMobs) and the black line shows the ionospheric RM (right axis label, RMion) along the LoS to this pulsar (as derived by the RMEXTRACT software package, based on the CODE+WMM). The bottom panel in each plot, the gray points … view at source ↗
Figure 2
Figure 2. The time series of electron density within the ionosphere along the LoS (TECLoS) as a function of local time during a day for nine representative observations, where the red points represent the TECLoS calculated by the pulsar signal (labelled as PSR J0814+7429), different colors and line styles (labeled on the top of the Figure) correspond to TECLoS obtained from different global ionospheric maps. The meanings of t… view at source ↗
Figure 3
Figure 3. The time series of RMs as a function of local time during a day, where the left and right plots correspond to the observation results from 01 July and 09 July, respectively. The red and black points in each plot were observed using the station DE609 and station DE605. The meanings of the gray area, yellow area, and the boundaries between them are consistent with [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

37 extracted references · 18 canonical work pages

  1. [1]

    Appleton, E. V. 1946, Nature, 157, 691, doi: 10.1038/157691a0

    doi:10.1038/157691a0 1946

  2. [2]

    A., & de Bruyn, A

    Brentjens, M. A., & de Bruyn, A. G. 2005, A&A, 441, 1217, doi: 10.1051/0004-6361:20052990

    doi:10.1051/0004-6361:20052990 2005

  3. [3]

    S., & Popov, M

    Burgin, M. S., & Popov, M. V. 2024, Astronomy Reports, 68, 257, doi: 10.1134/S1063772924700276

    doi:10.1134/s1063772924700276 2024

  4. [4]

    E., Hall, P

    Dewdney, P. E., Hall, P. J., Schilizzi, R. T., & Lazio, T. J. L. W. 2009, IEEE Proceedings, 97, 1482, doi: 10.1109/JPROC.2009.2021005

    arXiv 2009

  5. [5]

    Y., Verbiest, J

    Donner, J. Y., Verbiest, J. P. W., Tiburzi, C., et al. 2020, A&A, 644, A153, doi: 10.1051/0004-6361/202039517

    doi:10.1051/0004-6361/202039517 2020

  6. [6]

    2003, Advances in Space Research, 31, 635, doi: https://doi.org/10.1016/S0273-1177(03)00029-2 8

    Feltens, J. 2003, Advances in Space Research, 31, 635, doi: https://doi.org/10.1016/S0273-1177(03)00029-2 8

    work page doi:10.1016/s0273-1177(03)00029-2 2003

  7. [7]

    L., Zhou, D

    Han, J. L., Zhou, D. J., Wang, C., et al. 2025, Research in Astronomy and Astrophysics, 25, 014001, doi: 10.1088/1674-4527/ada3b7

    doi:10.1088/1674-4527/ada3b7 2025

  8. [8]

    2025, Satellite Navigation, 6, 11, doi: 10.1186/s43020-025-00164-x

    He, L., Zhu, Q., & Wang, C. 2025, Satellite Navigation, 6, 11, doi: 10.1186/s43020-025-00164-x

    work page doi:10.1186/s43020-025-00164-x 2025

Show all 37 references
  1. [9]

    2023, Advances in Space Research, 71, 1818, doi: 10.1016/j.asr.2022.09.042

    Jiang, C., Wang, W., Yang, G., & Zhao, Z. 2023, Advances in Space Research, 71, 1818, doi: 10.1016/j.asr.2022.09.042

    2023 doi

  2. [10]

    2024, Astronomical Techniques and Instruments, 1, 84, doi: 10.61977/ati2024012

    Jiang, P., Chen, R., Gan, H., et al. 2024, Astronomical Techniques and Instruments, 1, 84, doi: 10.61977/ati2024012

    2024 doi

  3. [11]

    2025, NewA, 121, 102460, doi: 10.1016/j.newast.2025.102460

    Kaur, D., Hobbs, G., Zic, A., et al. 2025, NewA, 121, 102460, doi: 10.1016/j.newast.2025.102460

    2025

  4. [12]

    2016, MNRAS, 458, 868, doi: 10.1093/mnras/stw189

    Lazarus, P., Karuppusamy, R., Graikou, E., et al. 2016, MNRAS, 458, 868, doi: 10.1093/mnras/stw189

    2016 doi

  5. [13]

    2023, Nature, 618, 484, doi: 10.1038/s41586-023-05983-z

    Li, D., Bilous, A., Ransom, S., Main, R., & Yang, Y.-P. 2023, Nature, 618, 484, doi: 10.1038/s41586-023-05983-z

    2023 doi

  6. [14]

    2025, Satellite Navigation, 6, 30, doi: 10.1186/s43020-025-00187-4

    Li, Z., Zhong, J., Wang, N., et al. 2025, Satellite Navigation, 6, 30, doi: 10.1186/s43020-025-00187-4

    2025 doi

  7. [15]

    Liang, P. H. 1947, Nature, 160, 642, doi: 10.1038/160642a0

    1947 doi

  8. [16]

    2026, Satellite Navigation, 7, 20, doi: 10.1186/s43020-026-00209-9

    Liu, A., Wang, N., Li, Z., et al. 2026, Satellite Navigation, 7, 20, doi: 10.1186/s43020-026-00209-9

    2026 doi

  9. [17]

    R., & Kramer, M

    Lorimer, D. R., & Kramer, M. 2012, Handbook of Pulsar Astronomy

    2012

  10. [18]

    Lynn, K. J. W., Gardiner-Garden, R. S., & Heitmann, A. 2014, Journal of Geophysical Research: Space Physics, 119, 10,294, doi: https://doi.org/10.1002/2014JA020617

    2014 doi

  11. [19]

    Manchester, R. N. 1971, ApJS, 23, 283, doi: 10.1086/190240

    1971 doi

  12. [20]

    N., Hobbs, G

    Manchester, R. N., Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993, doi: 10.1086/428488

    2005 doi

  13. [21]

    1998, Radio science, 33, 565

    Mannucci, A., Wilson, B., Yuan, D., et al. 1998, Radio science, 33, 565

    1998

  14. [22]

    2018, RMextract: Ionospheric Faraday Rotation calculator,, Astrophysics Source Code Library, record ascl:1806.024 http://ascl.net/1806.024

    Mevius, M. 2018, RMextract: Ionospheric Faraday Rotation calculator,, Astrophysics Source Code Library, record ascl:1806.024 http://ascl.net/1806.024

    2018

  15. [23]

    2025, Space Weather, 23, e2025SW004579, doi: 10.1029/2025SW004579 NCEI Geomagnetic Modeling Team, & British Geological Survey

    Nair, M., Fillion, M., Chulliat, A., & Califf, S. 2025, Space Weather, 23, e2025SW004579, doi: 10.1029/2025SW004579 NCEI Geomagnetic Modeling Team, & British Geological Survey. 2020, World Magnetic Model 2020, NOAA National Centers for Environmental Information, doi: 10.25921/...

    2025 doi

  16. [24]

    I., et al

    Noutsos, A., Sobey, C., Kondratiev, V. I., et al. 2015, A&A, 576, A62, doi: 10.1051/0004-6361/201425186

    2015 doi

  17. [25]

    1967, Nature, 216, 567, doi: 10.1038/216567a0

    Pacini, F. 1967, Nature, 216, 567, doi: 10.1038/216567a0

    1967 doi

  18. [26]

    A., Butler, B

    Perley, R. A., Butler, B. J., Greisen, E. W., et al. 2026, ApJS, 283, 82, doi: 10.3847/1538-4365/ae503c

    2026 doi

  19. [27]

    K., Noutsos, A., Tiburzi, C., et al

    Porayko, N. K., Noutsos, A., Tiburzi, C., et al. 2019, MNRAS, 483, 4100, doi: 10.1093/mnras/sty3324

    2019 doi

  20. [28]

    Rickett, B. J. 1977, ARA&A, 15, 479, doi: 10.1146/annurev.aa.15.090177.002403

    1977

  21. [29]

    Scheuer, P. A. G. 1968, Nature, 218, 920, doi: 10.1038/218920a0

    1968 doi

  22. [30]

    2024, Nature, 635, 365, doi: 10.1038/s41586-024-08072-x

    Smith, J., Kast, A., Geraschenko, A., et al. 2024, Nature, 635, 365, doi: 10.1038/s41586-024-08072-x

    2024 doi

  23. [31]

    Sotomayor-Beltran, C., Sobey, C., Hessels, J. W. T., et al. 2013, A&A, 552, A58, doi: 10.1051/0004-6361/201220728

    2013 doi

  24. [32]

    Tiburzi, C., Verbiest, J. P. W., Shaifullah, G. M., et al. 2019, MNRAS, 487, 394, doi: 10.1093/mnras/stz1278

    2019 doi

  25. [33]

    Seredkina, A. A. 2013, Baltic Astronomy, 22, 53, doi: 10.1515/astro-2017-0147 van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2, doi: 10.1051/0004-6361/201220873 van Straten, W., Demorest, P., & Oslowski, S. 2012, Astronomical Research and Technology, 9,...

    2013 doi arXiv

  26. [34]

    M., Manchester, R

    Yan, W. M., Manchester, R. N., Hobbs, G., et al. 2011, Ap&SS, 335, 485, doi: 10.1007/s10509-011-0756-0

    2011 doi

  27. [35]

    W., & Liu, Z

    Yu, S. W., & Liu, Z. Z. 2021, Satellite Navigation, 2, 20, doi: 10.1186/s43020-021-00051-1

    2021 doi

  28. [36]

    2021, Journal of Geodesy, 95, 35, doi: 10.1007/s00190-021-01487-8

    Zhao, J., Hern´ andez-Pajares, M., Li, Z., Wang, N., & Yuan, H. 2021, Journal of Geodesy, 95, 35, doi: 10.1007/s00190-021-01487-8

    2021 doi

  29. [37]

    I., Yermolaev, Y

    Zhuravlev, V. I., Yermolaev, Y. I., & Andrianov, A. S. 2020, MNRAS, 491, 5843, doi: 10.1093/mnras/stz3370

    2020 doi

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.