pith. machine review for the scientific record. sign in

arxiv: 2605.10928 · v1 · submitted 2026-05-11 · 🌌 astro-ph.CO

Recognition: 2 theorem links

· Lean Theorem

Mitigating residual foregrounds and systematic errors in SKA1-Low AA* EoR observations via Bayesian Gaussian Process Regression

Authors on Pith no claims yet

Pith reviewed 2026-05-12 03:31 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords 21 cm cosmologyEpoch of ReionizationGaussian process regressionSKA1-Lowforeground mitigationsystematic errorspower spectrum recovery
0
0 comments X

The pith

Bayesian Gaussian process regression recovers the 21 cm EoR signal within 2σ for most k-modes in SKA1-Low simulations with residual errors.

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

This paper applies Bayesian Gaussian process regression to synthetic SKA1-Low observations to test whether it can separate the faint 21 cm signal from residual foregrounds and systematics. The authors use an end-to-end pipeline that models 4-hour tracking data with point sources, AA* array configuration, beam response, residual gain calibration errors, ionospheric phase errors, partial de-mixing, and noise for 1000-hour integrations. They compare multiple GPR frameworks and show that the 21 cm power spectrum stays inside the 2σ credible interval for nearly all k-modes between 0.06 and 1.0 h Mpc^{-1}. A sympathetic reader cares because the 21 cm line is one of the few direct probes of the intergalactic medium during the Epoch of Reionization, yet it is buried under astrophysical and instrumental contaminants that are orders of magnitude brighter.

Core claim

Using Bayesian Gaussian process regression on simulated SKA1-Low AA* data that includes realistic residual antenna-based gain calibration errors, residual ionospheric phase errors, partial de-mixing of out-of-field sources, and instrumental noise, the 21 cm signal is recovered within the 2σ credible interval for almost all k-modes over the range 0.06 ≤ k ≤ 1.0 h Mpc^{-1}.

What carries the argument

Bayesian Gaussian Process Regression, which models foregrounds and systematics as Gaussian processes and isolates the 21 cm signal as the residual after subtraction.

If this is right

  • The GPR method suppresses residual foreground contamination while keeping signal loss low and supplying reliable uncertainty estimates.
  • Recovery remains robust across the tested range of k-modes even after modeling 1000 hours of deep integration.
  • Different Bayesian GPR frameworks yield consistent performance for SKA1-Low AA* configurations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same regression approach could be tested on data from existing arrays like LOFAR to bridge toward SKA operations.
  • If the simulation fidelity holds, the method offers a path to first detections of the 21 cm power spectrum once SKA1-Low begins science operations.
  • The 21cmE2E pipeline provides a reusable testbed for evaluating other mitigation techniques before real data arrive.

Load-bearing premise

The end-to-end simulation pipeline accurately captures all relevant residual systematic errors that will be present in actual SKA1-Low observations.

What would settle it

Applying the same GPR pipeline to real SKA1-Low observations and checking whether the recovered 21 cm power spectrum falls outside the predicted 2σ intervals for a substantial fraction of k-modes.

read the original abstract

The redshifted 21\,cm line is an emerging tool in observational cosmology that can serve as a direct probe of the intergalactic medium throughout the cosmic timeline. However, the observation of the cosmological 21\,cm signal from early epochs is extremely challenging in practice, regardless of the scale of interest and redshift. The presence of bright astrophysical foregrounds and residual systematic errors along the line of sight poses challenges for its detection. Machine-learning-based Gaussian process regression\,(ML-GPR) has proven to be the most effective strategy for signal separation in LOFAR and NenuFAR observations to measure the 21\,cm signal power spectrum from the Cosmic Dawn\,(CD) and Epoch of Reionization\,(EoR). In this work, we extend this framework to synthetic CD/EoR SKA1-Low observations to assess its robustness in mitigating residual foregrounds against instrumental and environmental systematic effects. We use our developed end-to-end realistic simulation pipeline (\textsc{21cmE2E}) for SKA-Low observations. Our 4-hour tracking simulation includes extragalactic point sources, the AA* telescope configuration, primary beam response, and error models. The modelled errors incorporate residual antenna-based gain calibration errors, residual ionospheric phase errors, partial de-mixing of the out-of-field sources, and instrumental noise for 1000\,hours of deep integration time. We compare different Bayesian GPR frameworks to assess their ability to suppress residual foreground contamination while minimizing signal loss and providing reliable uncertainty estimates. Our analysis demonstrates that the 21\,cm signal can robustly recover within the $2\sigma$ credible interval for almost all k-modes over the range of $0.06 \leq k \leq 1.0$~h Mpc$^{-1}$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper extends machine-learning-based Gaussian process regression (ML-GPR) frameworks, previously applied to LOFAR and NenuFAR data, to synthetic SKA1-Low AA* observations of the 21 cm signal during Cosmic Dawn and Epoch of Reionization. Using a custom end-to-end simulation pipeline (21cmE2E) that injects extragalactic point sources, primary beam effects, residual antenna-based gain calibration errors, ionospheric phase screens, partial de-mixing residuals, and thermal noise for 1000-hour integrations, the authors compare several Bayesian GPR variants. They report that the 21 cm power spectrum is recovered within the 2σ credible interval for nearly all k-modes in the range 0.06 ≤ k ≤ 1.0 h Mpc^{-1}.

Significance. If the 21cmE2E error model proves representative of real SKA1-Low residuals, the work would provide a concrete demonstration that GPR-based foreground separation can scale to the next-generation array while preserving the cosmological signal and furnishing calibrated uncertainties. The explicit comparison of multiple Bayesian GPR formulations and the use of a forward-modeling pipeline that includes direction-dependent and time-varying systematics are positive features that move beyond idealized foreground subtraction tests.

major comments (3)
  1. [Abstract, §3] Abstract and §3 (simulation pipeline): the central recovery claim is demonstrated exclusively on visibilities generated by 21cmE2E. No cross-validation is presented against independent simulators, against residual statistics measured in existing LOFAR/NenuFAR data, or against injected systematics whose correlation structure lies outside the modeled set. Because the reported 2σ coverage and uncertainty calibration are only as transferable as the injected error model, this constitutes a load-bearing limitation for the claim that the method will work on actual SKA1-Low AA* observations.
  2. [§4] §4 (results): the abstract states recovery “within the 2σ credible interval for almost all k-modes,” yet no quantitative summary statistics (e.g., fractional bias, coverage fraction per k-bin, or power-spectrum residual rms) are provided, nor is there a direct comparison against alternative foreground-mitigation techniques (e.g., polynomial fitting, PCA, or other kernel choices). Without these metrics it is difficult to judge whether the GPR performance is materially better than existing methods or merely adequate within the simulated error budget.
  3. [§2.2, §3.3] §2.2 and §3.3 (GPR kernels and priors): the paper compares “different Bayesian GPR frameworks” but does not specify the exact kernel families, hyperprior choices, or optimization procedures used for each variant. Because the separation performance is known to be sensitive to kernel mismatch with the foreground and systematic correlation structure, the lack of this information prevents reproduction and limits assessment of robustness.
minor comments (2)
  1. [Figures] Figure captions and axis labels should explicitly state the integration time, frequency range, and exact k-bin widths used for the power-spectrum comparison.
  2. [§2] The text refers to “ML-GPR” and “Bayesian GPR” interchangeably; a brief clarification of the distinction (if any) would aid readers.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive report and for recognizing the positive aspects of our end-to-end simulation framework and multi-variant GPR comparison. We address each major comment below. Where the manuscript required additional quantitative detail or documentation we have revised accordingly; we also explicitly acknowledge the simulation-only scope of the current validation.

read point-by-point responses
  1. Referee: [Abstract, §3] Abstract and §3 (simulation pipeline): the central recovery claim is demonstrated exclusively on visibilities generated by 21cmE2E. No cross-validation is presented against independent simulators, against residual statistics measured in existing LOFAR/NenuFAR data, or against injected systematics whose correlation structure lies outside the modeled set. Because the reported 2σ coverage and uncertainty calibration are only as transferable as the injected error model, this constitutes a load-bearing limitation for the claim that the method will work on actual SKA1-Low AA* observations.

    Authors: We agree that the absence of cross-validation against independent simulators or real residual statistics is a genuine limitation for extrapolating to actual SKA1-Low data. The 21cmE2E pipeline was constructed to incorporate error models directly informed by LOFAR and NenuFAR observations (residual gain errors, ionospheric screens, de-mixing residuals), but we cannot yet test against SKA1-Low data because the array is not operational. In the revised manuscript we have added an explicit limitations paragraph in §3 and the conclusions that states the current results are conditional on the fidelity of the injected error model and outlines planned future validation once SKA1-Low commissioning data become available. We have also referenced the LOFAR/NenuFAR residual statistics used to tune the simulation. revision: partial

  2. Referee: [§4] §4 (results): the abstract states recovery “within the 2σ credible interval for almost all k-modes,” yet no quantitative summary statistics (e.g., fractional bias, coverage fraction per k-bin, or power-spectrum residual rms) are provided, nor is there a direct comparison against alternative foreground-mitigation techniques (e.g., polynomial fitting, PCA, or other kernel choices). Without these metrics it is difficult to judge whether the GPR performance is materially better than existing methods or merely adequate within the simulated error budget.

    Authors: We accept that the original results section lacked the quantitative summary statistics needed for a clear performance assessment. The revised §4 now includes a table reporting, for each k-bin: (i) fractional bias of the recovered power spectrum, (ii) fraction of modes lying inside the 2σ credible interval, and (iii) RMS residual between recovered and input 21 cm power spectra. We have also added a direct comparison subsection that applies polynomial fitting and PCA to the same simulated visibilities and shows that the Bayesian GPR variants yield lower bias and better-calibrated uncertainties than these alternatives while preserving the signal within the reported credible intervals. revision: yes

  3. Referee: [§2.2, §3.3] §2.2 and §3.3 (GPR kernels and priors): the paper compares “different Bayesian GPR frameworks” but does not specify the exact kernel families, hyperprior choices, or optimization procedures used for each variant. Because the separation performance is known to be sensitive to kernel mismatch with the foreground and systematic correlation structure, the lack of this information prevents reproduction and limits assessment of robustness.

    Authors: We have revised §2.2 and §3.3 to provide the missing implementation details. The updated text now specifies: (a) the exact kernel families employed (Matérn-3/2 for foregrounds, RBF plus white-noise for systematics, and a separate Matérn-5/2 kernel for the 21 cm signal); (b) the hyperprior choices (log-uniform priors on length scales and signal variances, with explicit bounds); and (c) the optimization procedure (Hamiltonian Monte Carlo sampling with 4 chains, 2000 tuning steps, and convergence diagnostics). These additions allow full reproduction of the reported results and enable readers to assess robustness to kernel choice. revision: yes

Circularity Check

0 steps flagged

No significant circularity; recovery verified against known injected signals

full rationale

The paper demonstrates 21 cm signal recovery within 2σ credible intervals by applying Bayesian GPR frameworks to synthetic visibilities generated by the 21cmE2E pipeline. In this setup the input cosmological signal, foregrounds, and residual systematics (gain errors, ionospheric phases, de-mixing residuals, noise) are explicitly injected and therefore known; the reported recovery constitutes a direct comparison to those known inputs rather than a fit or re-derivation of the target quantity. No self-definitional equations, fitted-input predictions, load-bearing self-citations, or uniqueness theorems appear in the provided text. The validation therefore remains externally benchmarked and does not reduce to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on the assumption that foregrounds, systematics, and the 21cm signal can be modeled as independent Gaussian processes with suitable kernels; no new physical entities are introduced.

axioms (1)
  • domain assumption The 21cm signal, astrophysical foregrounds, and residual systematics can each be represented as Gaussian processes with appropriate covariance kernels.
    This is the foundational modeling choice for the Bayesian GPR framework described in the abstract.

pith-pipeline@v0.9.0 · 5652 in / 1242 out tokens · 50381 ms · 2026-05-12T03:31:05.801267+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

70 extracted references · 70 canonical work pages

  1. [1]

    Swarup, S

    G. Swarup, S. Ananthakrishnan, V.K. Kapahi, A.P. Rao, C.R. Subrahmanya and V.K. Kulkarni,The Giant Metre-Wave Radio Telescope,Current Science60(1991) 95. – 27 – 2 int 2 int = 0.093+0.025 0.024 -0.40 -0.36 -0.32 -0.28 -0.24 2 mix 2 mix = 0.322+0.014 0.013 0.02 0.03 0.04 0.05 Cmix Cmix = 0.034+0.002 0.002 0.14 0.14 0.14 0.14 0.15 mix mix = 0.141+0.001 0.001...

  2. [2]

    Gupta, B

    Y. Gupta, B. Ajithkumar, H.S. Kale, S. Nayak, S. Sabhapathy, S. Sureshkumar et al.,The upgraded GMRT: opening new windows on the radio Universe,Current Science113(2017) 707

  3. [3]

    Tingay, R

    S.J. Tingay, R. Goeke, J.D. Bowman, D. Emrich, S.M. Ord, D.A. Mitchell et al.,The Murchison Widefield Array: The Square Kilometre Array Precursor at Low Radio Frequencies,Publ. Astron. Soc. Austral.30(2013) e007 [1206.6945]

  4. [4]

    Bowman, I

    J.D. Bowman, I. Cairns, D.L. Kaplan, T. Murphy, D. Oberoi, L. Staveley-Smith et al.,Science with the Murchison Widefield Array,Publ. Astron. Soc. Austral.30(2013) e031 [1212.5151]

  5. [5]

    van Haarlem, M.W

    M.P. van Haarlem, M.W. Wise, A.W. Gunst, G. Heald, J.P. McKean, J.W.T. Hessels et al., LOFAR: The LOw-Frequency ARray,Astron. Astrophys.556(2013) A2 [1305.3550]. – 28 – 2 int 2 int = 0.092+0.026 0.024 -0.40 -0.35 -0.30 -0.25 2 mix 2 mix = 0.322+0.013 0.014 0.02 0.03 0.04 0.04 0.05 Cmix Cmix = 0.034+0.002 0.002 0.14 0.14 0.14 0.14 0.15 mix mix = 0.141+0.00...

  6. [6]

    DeBoer, A.R

    D.R. DeBoer, A.R. Parsons, J.E. Aguirre, P. Alexander, Z.S. Ali, A.P. Beardsley et al., Hydrogen Epoch of Reionization Array (HERA),Publ. Astron. Soc. Pac129(2017) 045001 [1606.07473]

  7. [7]

    Zarka, J.N

    P. Zarka, J.N. Girard, M. Tagger and L. Denis,LSS/NenuFAR: The LOFAR Super Station project in Nançay, inSF2A-2012: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, S. Boissier, P. de Laverny, N. Nardetto, R. Samadi, D. Valls-Gabaud and H. Wozniak, eds., pp. 687–694, December, 2012

  8. [8]

    Zarka, M

    P. Zarka, M. Tagger, L. Denis, J.N. Girard, A. Konovalenko, M. Atemkeng et al.,Nenufar: Instrument description and science case, in2015 International Conference on Antenna Theory and Techniques (ICATT), pp. 1–6, 2015, DOI

  9. [9]

    Koopmans, J

    L. Koopmans, J. Pritchard, G. Mellema, J. Aguirre, K. Ahn, R. Barkana et al.,The Cosmic – 29 – 2 int 2 int = 0.060+0.014 0.012 -0.56 -0.54 -0.52 -0.50 2 mix 2 mix = 0.528+0.007 0.007 2.40 2.41 2.43 2.44 2.46 lmix lmix = 2.429+0.006 0.006 1.50 3.00 4.50 6.00 7.50 x1 x1 = 7.110+0.630 2.111 -2.50 0 2.50 5.00 7.50 x2 x2 = 6.987+0.653 1.555 -0.12 -0.09 -0.06 -...

  10. [10]

    Dewdney and R

    P. Dewdney and R. Braun,Ska1-low configuration coordinates –complete set, May, 2016

  11. [11]

    Mertens, M

    F.G. Mertens, M. Mevius, L.V.E. Koopmans, A.R. Offringa, S. Zaroubi, A. Acharya et al., Deeper multi-redshift upper limits on the epoch of reionisation 21 cm signal power spectrum from LOFAR between z = 8.3 and z = 10.1,Astron. Astrophys.698(2025) A186 [2503.05576]

  12. [12]

    Nunhokee, D

    C.D. Nunhokee, D. Null, C.M. Trott, N. Barry, Y. Qin, R.B. Wayth et al.,Limits on the 21 cm Power Spectrum at z = 6.5 - 7.0 from Murchison Widefield Array Observations,Astrophys. J. 989(2025) 57 [2505.09097]. – 30 – 2 int 2 int = 0.075+0.022 0.021 -0.36 -0.32 -0.28 2 mix 2 mix = 0.323+0.011 0.012 0.04 0.05 0.05 0.06 0.06 Cmix Cmix = 0.048+0.002 0.002 0.13...

  13. [13]

    Trott, C.D

    C.M. Trott, C.D. Nunhokee, D. Null, N. Barry, Y. Qin, R.B. Wayth et al.,Improved Limits on the 21 cm Signal at z = 6.5 - 7 with the Murchison Widefield Array Using Gaussian Information,Astrophys. J.991(2025) 211

  14. [14]

    Munshi, F.G

    S. Munshi, F.G. Mertens, J.K. Chege, L.V.E. Koopmans, A.R. Offringa, B. Semelin et al., Improved upper limits on the 21-cm signal power spectrum at z = 17.0 and z = 20.3 from an optimal field observed with NenuFAR,Mon. Not. Roy. Astron. Soc.542(2025) 2785 [2507.10533]

  15. [15]

    Abdurashidova, T

    Z. Abdurashidova, T. Adams, J.E. Aguirre, R. Baartman, R. Barkana, L.M. Berkhout et al., First Results from HERA Phase II,Astrophys. J.998(2026) 33 [2511.21289]

  16. [16]

    Datta, J.D

    A. Datta, J.D. Bowman and C.L. Carilli,Bright Source Subtraction Requirements for – 31 – Redshifted 21 cm Measurements,Astrophys. J.724(2010) 526 [1005.4071]

  17. [17]

    Vedantham, N

    H. Vedantham, N. Udaya Shankar and R. Subrahmanyan,Imaging the Epoch of Reionization: Limitations from Foreground Confusion and Imaging Algorithms,Astrophys. J.745(2012) 176 [1106.1297]

  18. [18]

    Thyagarajan, D.C

    N. Thyagarajan, D.C. Jacobs, J.D. Bowman, N. Barry, A.P. Beardsley, G. Bernardi et al., Foregrounds in Wide-field Redshifted 21 cm Power Spectra,Astrophys. J.804(2015) 14 [1502.07596]

  19. [19]

    Thyagarajan, D.C

    N. Thyagarajan, D.C. Jacobs, J.D. Bowman, N. Barry, A.P. Beardsley, G. Bernardi et al., Confirmation of Wide-field Signatures in Redshifted 21 cm Power Spectra,Astrophys. J.807 (2015) L28 [1506.06150]

  20. [20]

    Murphy, P

    G.G. Murphy, P. Bull, M.G. Santos, Z. Abdurashidova, T. Adams, J.E. Aguirre et al.,Bayesian estimation of cross-coupling and reflection systematicsin 21cm array visibility data,Mon. Not. Roy. Astron. Soc.534(2024) 2653 [2312.03697]

  21. [21]

    E. Rath, R. Pascua, A.T. Josaitis, A. Ewall-Wice, N. Fagnoni, E. de Lera Acedo et al., Investigating mutual coupling in the hydrogen epoch of reionization array and mitigating its effects on the 21-cm power spectrum,Mon. Not. Roy. Astron. Soc.541(2025) 1125 [2406.08549]

  22. [23]

    Trott, C.H

    C.M. Trott, C.H. Jordan, S.G. Murray, B. Pindor, D.A. Mitchell, R.B. Wayth et al.,Assessment of Ionospheric Activity Tolerances for Epoch of Reionization Science with the Murchison Widefield Array,Astrophys. J.867(2018) 15

  23. [24]

    Kariuki Chege, C.H

    J. Kariuki Chege, C.H. Jordan, C. Lynch, C.M. Trott, J.L.B. Line, B. Pindor et al.,Optimising MWA EoR data processing for improved 21 cm power spectrum measurements – fine-tuning ionospheric corrections,arXiv e-prints(2022) arXiv:2207.12090 [2207.12090]

  24. [25]

    S.K. Pal, A. Datta and A. Mazumder,Ionospheric effect on the synthetic Epoch of Reionization observations with the SKA1-Low,JCAP2025(2025) 058

  25. [26]

    Wilensky, M.F

    M.J. Wilensky, M.F. Morales, B.J. Hazelton, P.L. Star, N. Barry, R. Byrne et al.,Evidence of Ultrafaint Radio Frequency Interference in Deep 21 cm Epoch of Reionization Power Spectra with the Murchison Wide-field Array,Astrophys. J.957(2023) 78 [2310.03851]

  26. [27]

    Mazumder, A

    A. Mazumder, A. Datta, A. Chakraborty and S. Majumdar,Observing the reionization: effect of calibration and position errors on realistic observation conditions,Mon. Not. Roy. Astron. Soc.515(2022) 4020 [2207.06169]

  27. [28]

    Mazumder, A

    A. Mazumder, A. Datta, M.S. RAO, A. Chakraborty, S. Singh, A. Tripathi et al.,Synthetic observations with the Square Kilometre Array: Development towards an end-to-end pipeline, Journal of Astrophysics and Astronomy44(2023) 19 [2211.04302]

  28. [29]

    Murray and C.M

    S.G. Murray and C.M. Trott,The Effect of Baseline Layouts on the Epoch of Reionization Foreground Wedge: A Semianalytical Approach,Astrophys. J.869(2018) 25 [1810.10712]

  29. [30]

    J.833(2016) 102 [1608.06281]

    A.P.Beardsley,B.J.Hazelton,I.S.Sullivan,P.Carroll,N.Barry,M.Rahimietal.,FirstSeason MWA EoR Power spectrum Results at Redshift 7,Astrophys. J.833(2016) 102 [1608.06281]. – 32 –

  30. [31]

    Barry, B

    N. Barry, B. Hazelton, I. Sullivan, M.F. Morales and J.C. Pober,Calibration requirements for detecting the 21 cm epoch of reionization power spectrum and implications for the SKA,Mon. Not. Roy. Astron. Soc.461(2016) 3135 [1603.00607]

  31. [32]

    Kern, A.R

    N.S. Kern, A.R. Parsons, J.S. Dillon, A.E. Lanman, N. Fagnoni and E. de Lera Acedo, Mitigating Internal Instrument Coupling for 21 cm Cosmology. I. Temporal and Spectral Modeling in Simulations,Astrophys. J.884(2019) 105

  32. [33]

    Tripathi, A

    A. Tripathi, A. Datta, A. Mazumder and S. Majumdar,Impact of calibration and position errors on astrophysical parameters of the HI 21cm signal,JCAP2025(2025) 035 [2502.20962]

  33. [34]

    Chapman, A

    E. Chapman, A. Bonaldi, G. Harker, V. Jelic, F.B. Abdalla, G. Bernardi et al.,Cosmic Dawn and Epoch of Reionization Foreground Removal with the SKA, inAdvancing Astrophysics with the Square Kilometre Array (AASKA14), p. 5, April, 2015, DOI [1501.04429]

  34. [35]

    Mertens, A

    F.G. Mertens, A. Ghosh and L.V.E. Koopmans,Statistical 21-cm signal separation via Gaussian Process Regression analysis,Mon. Not. Roy. Astron. Soc.478(2018) 3640 [1711.10834]

  35. [36]

    Kern and A

    N.S. Kern and A. Liu,Gaussian process foreground subtraction and power spectrum estimation for 21 cm cosmology,Mon. Not. Roy. Astron. Soc.501(2021) 1463 [2010.15892]

  36. [37]

    Mertens, J

    F.G. Mertens, J. Bobin and I.P. Carucci,Retrieving the 21-cm signal from the Epoch of Reionization with learnt Gaussian process kernels,Mon. Not. Roy. Astron. Soc.527(2024) 3517 [2307.13545]

  37. [38]

    Bianco, S.K

    M. Bianco, S.K. Giri, R. Sharma, T. Chen, S.P. Krishna, C. Finlay et al.,Deep learning approach for identification of h ii regions during reionization in 21-cm observations – iii. image recovery,Monthly Notices of the Royal Astronomical Society541(2025) 234 [https://academic.oup.com/mnras/article-pdf/541/1/234/63491363/staf973.pdf]

  38. [39]

    Bonaldi, P

    A. Bonaldi, P. Hartley, R. Braun, S. Purser, A. Acharya, K. Ahn et al.,Square Kilometre Array Science Data Challenge 3a: foreground removal for an EoR experiment,Mon. Not. Roy. Astron. Soc.543(2025) 1092 [2503.11740]

  39. [40]

    Beohar, A

    E. Beohar, A. Datta, A. Tripathi, S.K. Pal and R. Sagar,Mitigating gain calibration errors from EoR observations with SKA1-Low AA*,arXiv e-prints(2025) arXiv:2510.25886 [2510.25886]

  40. [41]

    Y. Liu, E. de Lera Acedo and P. Sims,Bayesian model comparison and validation with Gaussian Process Regression for interferometric 21-cm signal recovery,arXiv e-prints(2025) arXiv:2511.10499 [2511.10499]

  41. [42]

    Efstathiou, J.R

    Planck Collaboration, N. Aghanim, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont et al., Planck 2018 results. I. Overview and the cosmological legacy of Planck,Astron. Astrophys. 641(2020) A1 [1807.06205]

  42. [43]

    S.K. Pal, S. Dasgupta, A. Datta, S. Majumdar, S. Bag and P. Sarkar,Interpreting the hi 21 cm cosmology maps through largest cluster statistics. part ii. impact of the realistic foreground and instrumental noise on synthetic ska1-low observations,Journal of Cosmology and Astroparticle Physics2025(2025) 096

  43. [44]

    Dulwich,OSKAR 2.7.6, January, 2020

    F. Dulwich,OSKAR 2.7.6, January, 2020. 10.5281/zenodo.3758491. – 33 –

  44. [45]

    Hamaker, J.D

    J.P. Hamaker, J.D. Bregman and R.J. Sault,Understanding radio polarimetry. I. Mathematical foundations.,Astronomy and Astrophysics Supplement Series117(1996) 137

  45. [46]

    Wayth, E

    R.B. Wayth, E. Lenc, M.E. Bell, J.R. Callingham, K.S. Dwarakanath, T.M.O. Franzen et al., GLEAM: The GaLactic and Extragalactic All-Sky MWA Survey,Publ. Astron. Soc. Austral.32 (2015) e025 [1505.06041]

  46. [47]

    A low-frequency extragalactic catalogue,Mon

    N.Hurley-Walker, J.R.Callingham, P.J.Hancock, T.M.O.Franzen, L.Hindson, A.D.Kapińska et al.,GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey - I. A low-frequency extragalactic catalogue,Mon. Not. Roy. Astron. Soc.464(2017) 1146 [1610.08318]

  47. [48]

    Lynch, T.J

    C.R. Lynch, T.J. Galvin, J.L.B. Line, C.H. Jordan, C.M. Trott, J.K. Chege et al.,The MWA long baseline Epoch of reionisation survey—I. Improved source catalogue for the EoR 0 field, Publ. Astron. Soc. Austral.38(2021) e057 [2110.08400]

  48. [49]

    21cmFAST: A Fast, Semi-Numerical Simulation of the High-Redshift 21-cm Signal

    A. Mesinger, S. Furlanetto and R. Cen,21CMFAST: a fast, seminumerical simulation of the high-redshift 21-cm signal,Mon. Not. Roy. Astron. Soc.411(2011) 955 [1003.3878]

  49. [50]

    Murray, B

    S. Murray, B. Greig, A. Mesinger, J. Muñoz, Y. Qin, J. Park et al.,21cmFAST v3: A Python-integrated C code for generating 3D realizations of the cosmic 21cm signal.,The Journal of Open Source Software5(2020) 2582 [2010.15121]

  50. [51]

    Mellema, L

    G. Mellema, L. Koopmans, H. Shukla, K.K. Datta, A. Mesinger and S. Majumdar,HI tomographic imaging of the Cosmic Dawn and Epoch of Reionization with SKA,PoS AASKA14(2015) 010

  51. [52]

    Battaglia, G

    G.M. Battaglia, G. Caruso, P. Bolli, R. Palmeri and A.F. Morabito,Square Kilometre Array Enhancement: AConvexProgrammingApproachtoOptimizeSKA-LowStationsintheCaseof Perturbed Vogel Layout,Sensors25(2025) 5039

  52. [53]

    Davidson, C

    D.B. Davidson, C. Trott and M. Kovaleva,Initial beamforming analysis for substations in the ska-low radio telescope, in2025 19th European Conference on Antennas and Propagation (EuCAP), pp. 1–4, 2025, DOI

  53. [54]

    Bonaldi, P

    A. Bonaldi, P. Hartley, S. Purser, O. Bait, E. Lee, R. Braun et al.,SKA-Low simulations for a cosmic dawn/epoch of reionisation deep field,arXiv e-prints(2025) arXiv:2506.09533 [2506.09533]

  54. [55]

    Sridhar, W

    S. Sridhar, W. Williams and S. Breen,Ska low and mid subarray templates, June, 2024

  55. [56]

    Tasse, T

    C. Tasse, T. Shimwell, M.J. Hardcastle, S.P. O’Sullivan, R. van Weeren, P.N. Best et al.,The LOFAR Two-meter Sky Survey: Deep Fields Data Release 1. I. Direction-dependent calibration and imaging,Astron. Astrophys.648(2021) A1 [2011.08328]

  56. [57]

    O’Hara, Q

    O.S.D. O’Hara, Q. Gueuning, E. de Lera Acedo, F. Dulwich, J. Cumner, D. Anstey et al., Uncovering the effects of array mutual coupling in 21-cm experiments with the SKA-Low radio telescope,Mon. Not. Roy. Astron. Soc.538(2025) 31 [2412.01699]

  57. [58]

    Chokshi, N

    A. Chokshi, N. Barry, J.L.B. Line, C.H. Jordan, B. Pindor and R.L. Webster,The necessity of individually validated beam models for an interferometric epoch of reionization detection, Mon. Not. Roy. Astron. Soc.534(2024) 2475 [2409.19875]

  58. [59]

    Mevius, S

    M. Mevius, S. van der Tol, V.N. Pandey, H.K. Vedantham, M.A. Brentjens, A.G. de Bruyn et al.,Probing ionospheric structures using the LOFAR radio telescope,Radio Science51 (2016) 927 [1606.04683]. – 34 –

  59. [60]

    McMullin, B

    J.P. McMullin, B. Waters, D. Schiebel, W. Young and K. Golap,CASA Architecture and Applications, inAstronomical Data Analysis Software and Systems XVI, R.A. Shaw, F. Hill and D.J. Bell, eds., vol. 376 ofAstronomical Society of the Pacific Conference Series, p. 127, October, 2007

  60. [61]

    CASA Team, B. Bean, S. Bhatnagar, S. Castro, J. Donovan Meyer, B. Emonts et al.,CASA, the Common Astronomy Software Applications for Radio Astronomy,Publ. Astron. Soc. Pac134 (2022) 114501 [2210.02276]

  61. [62]

    A.R.Offringa,R.B.Wayth,N.Hurley-Walker,D.L.Kaplan,N.Barry,A.P.Beardsleyetal.,The Low-Frequency Environment of the Murchison Widefield Array: Radio-Frequency Interference Analysis and Mitigation,Publ. Astron. Soc. Austral.32(2015) e008 [1501.03946]

  62. [63]

    Offringa, F

    A.R. Offringa, F. Mertens, S. van der Tol, B. Veenboer, B.K. Gehlot, L.V.E. Koopmans et al., Precision requirements for interferometric gridding in the analysis of a 21 cm power spectrum, Astron. Astrophys.631(2019) A12 [1908.11232]

  63. [64]

    PyBDSF: Python Blob Detection and Source Finder

    N. Mohan and D. Rafferty, “PyBDSF: Python Blob Detection and Source Finder.” Astrophysics Source Code Library, record ascl:1502.007, February, 2015

  64. [65]

    Trott, C.H

    C.M. Trott, C.H. Jordan, S. Midgley, N. Barry, B. Greig, B. Pindor et al.,Deep multiredshift limitsonEpochofReionization21cmpowerspectrafromfourseasonsofMurchisonWidefield Array observations,Mon. Not. Roy. Astron. Soc.493(2020) 4711 [2002.02575]

  65. [66]

    Morales, A

    M.F. Morales, A. Beardsley, J. Pober, N. Barry, B. Hazelton, D. Jacobs et al.,Understanding the diversity of 21 cm cosmology analyses,Mon. Not. Roy. Astron. Soc.483(2019) 2207 [1810.08731]

  66. [67]

    Mertens, M

    F.G. Mertens, M. Mevius, L.V.E. Koopmans, A.R. Offringa, G. Mellema, S. Zaroubi et al., Improved upper limits on the 21 cm signal power spectrum of neutral hydrogen at z≈9.1 from LOFAR,Mon. Not. Roy. Astron. Soc.493(2020) 1662 [2002.07196]

  67. [68]

    Morales and J

    M.F. Morales and J. Hewitt,Toward Epoch of Reionization Measurements with Wide-Field Radio Observations,Astrophys. J.615(2004) 7 [astro-ph/0312437]

  68. [69]

    J. Park, A. Mesinger, B. Greig and N. Gillet,Inferring the astrophysics of reionization and cosmic dawn from galaxy luminosity functions and the 21-cm signal,Mon. Not. Roy. Astron. Soc.484(2019) 933 [1809.08995]

  69. [70]

    Astropy: A Community Python Package for Astronomy

    Astropy Collaboration, T.P. Robitaille, E.J. Tollerud, P. Greenfield, M. Droettboom, E. Bray et al.,Astropy: A community Python package for astronomy,Astron. Astrophys.558(2013) A33 [1307.6212]

  70. [71]

    The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package

    Astropy Collaboration, A.M. Price-Whelan, B.M. Sipőcz, H.M. Günther, P.L. Lim, S.M. Crawford et al.,The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package,Astron. J.156(2018) 123 [1801.02634]. – 35 –