Bayesian power spectrum estimation with modelling of systematic effects in delay-fringe rate space
Pith reviewed 2026-05-07 12:10 UTC · model grok-4.3
The pith
A Bayesian model treats cable reflections as multiplicative factors in delay-fringe rate space inside a Gibbs sampler, recovering the 21cm delay power spectrum from simulated visibilities without extra filtering losses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We demonstrate the method on simulated visibility data for a single baseline, showing that the 21cm delay power spectrum can be recovered well regardless of the location of the systematics in delay-fringe rate space.
Load-bearing premise
The systematics can be accurately represented as a multiplicative effect in delay-fringe rate space whose parameters are jointly inferred with the 21cm signal and foregrounds; this is stated as an extension of the hydra-pspec model but its validity for real data is not demonstrated.
read the original abstract
Observing the Epoch of Reionisation using 21cm radio interferometry has proven to be a challenging task. Extraction of the extremely faint redshifted signal is complicated by the presence of bright foregrounds, radio frequency interference (RFI), and systematic artefacts. We discuss the challenge of accounting for systematic effects, particularly cable reflections, that appear in the visibility data obtained from 21cm interferometers. Cable reflections cause attenuated copies of the foreground signal to appear outside the 'foreground wedge' region in which foreground contamination is supposed to be localised. We build on the hydra-pspec Gibbs sampler to implement a model of the systematics as a multiplicative effect in delay-fringe rate space. We include this model in the inference of the joint posterior distribution, in addition to the 21cm signal, its power spectrum, and foregrounds. This allows the systematics contribution to be marginalised, rather than filtering it out and causing additional signal loss. We demonstrate the method on simulated visibility data for a single baseline, showing that the 21cm delay power spectrum can be recovered well regardless of the location of the systematics in delay-fringe rate space. Our implementation is suitable for modelling other multiplicative factors on the visibilities, e.g. residual gain errors.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript extends the hydra-pspec Gibbs sampler to jointly infer 21 cm signal, power spectrum, foregrounds, and systematic effects (e.g., cable reflections) modeled as multiplicative factors in delay-fringe rate space. This allows marginalization over the systematics rather than filtering. The method is demonstrated on simulated single-baseline visibility data, with the claim that the 21 cm delay power spectrum is recovered well independent of the location of the systematics in delay-fringe rate space. The approach is also positioned as suitable for other multiplicative visibility effects such as residual gain errors.
Significance. The joint Bayesian inference framework is a clear strength, as it propagates uncertainties from the systematics parameters into the power spectrum posterior and avoids the signal loss associated with aggressive filtering. If the multiplicative model proves adequate for real data, the method could improve robustness in Epoch of Reionization analyses. The paper provides a clear demonstration on matched simulations and notes extensibility to other effects, but the current validation scope limits the immediate impact.
major comments (1)
- [Results / Simulation demonstration] The results demonstration uses simulated visibilities generated from the exact multiplicative model in delay-fringe rate space that the sampler is constructed to fit. Successful recovery is therefore expected by construction upon marginalization, but this does not test whether the multiplicative representation adequately describes real cable reflections (physically closer to additive delayed echoes with frequency-dependent attenuation). This assumption is load-bearing for the central claim of unbiased 21 cm power spectrum recovery in the presence of systematics.
minor comments (2)
- [Abstract and Results] The abstract and results section would benefit from quantitative metrics (bias, error bars, or direct comparison to filtering approaches) rather than the qualitative statement that the power spectrum is 'recovered well'.
- [Discussion / Future work] The demonstration is restricted to a single baseline; discussion of scaling to multi-baseline arrays and baseline-dependent systematics would strengthen applicability.
Simulated Author's Rebuttal
We thank the referee for their constructive review of our manuscript. The major comment raises an important point about the scope of our simulation tests, which we address in detail below. We have revised the manuscript to better contextualize the results and acknowledge the limitations of the current validation.
read point-by-point responses
-
Referee: The results demonstration uses simulated visibilities generated from the exact multiplicative model in delay-fringe rate space that the sampler is constructed to fit. Successful recovery is therefore expected by construction upon marginalization, but this does not test whether the multiplicative representation adequately describes real cable reflections (physically closer to additive delayed echoes with frequency-dependent attenuation). This assumption is load-bearing for the central claim of unbiased 21 cm power spectrum recovery in the presence of systematics.
Authors: We agree with the referee that our simulation demonstration uses data generated from the exact model assumed in the sampler, making successful recovery expected when the model is correct. This serves to validate the joint inference framework and the marginalization over systematics parameters. However, we note that the manuscript positions the multiplicative model as an approximation suitable for certain systematic effects, such as residual gain errors, rather than claiming it exactly captures all aspects of cable reflections. In the revised version, we have added a new paragraph in the discussion section explicitly addressing the physical nature of cable reflections as potentially additive, the limitations of the multiplicative approximation in delay-fringe rate space, and the need for future studies with more realistic simulations including frequency-dependent effects. We have also tempered the language around the central claim to specify that the recovery is unbiased under the assumed model. This revision clarifies the scope without altering the core contribution of the Bayesian approach. revision: yes
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.