Efficient GLS map-making and assessment of 1/f noise impact for future CMB satellite missions
Abstract
Time-correlated 1/f noise is a major source of large-scale contamination in measurements of the polarized cosmic microwave background (CMB) that can limit sensitivity to primordial B modes. Generalized least-squares (GLS) map-making is optimal for stationary Gaussian noise when its covariance is known, but practical implementations rely on approximations to make the computation tractable. In particular, approximating the noise covariance as circulant enables efficient map-making using fast Fourier transform, but introduces artificial edge effects when applied to finite data chunks. We introduce STITCH (STatistical Inpainting for Timeline CHunks), a preprocessing method that mitigates these effects by extending time-ordered data chunks with synthetic samples that make the extended timelines statistically consistent with the assumed circulant noise model. They can then be processed with standard circulant GLS solvers, preserving their computational efficiency while reducing residuals associated with chunk boundaries. Simulations show that STITCH substantially suppresses these edge effects, bringing the resulting noise power spectra close to those obtained when the circulant approximation holds. We demonstrate negligible impact of our procedure on the recovered signal, such as E→ B leakage, even at large scales. Finally, we investigate how temporal 1/f noise projects onto angular scales for representative LiteBIRD-, \emph{Planck}-, and PICO-like scanning strategies. Using noise simulations spanning a range of temporal noise parameters, we fit an effective model to the resulting angular noise spectra. These results provide an empirical characterization of the temporal-to-angular noise mapping for the configurations studied and a basis for further investigations of correlated noise in future CMB polarization missions.
Related
- A measurement-based approach to EPFD calculation, to quantify the impact of satellite mega-constellations on low-frequency radio astronomy
- Exploring non-Poisson satellite occupation in HOD models and its impact on 2- and 3-point galaxy clustering
- Satellite Streaks in HETDEX: Spectra, Orbital Identification, and Survey Impact
- A measurement-based approach to equivalent power flux density calculation, to quantify the impact of satellite mega-constellations on low-frequency radio astronomy
- Testing black holes with interstellar missions: I. Orbiting probes
BibTeX
@misc{monelli2026efficient,
title = {{Efficient GLS map-making and assessment of $1/f$ noise impact for future CMB satellite missions}},
author = {Marta Monelli and Guillaume Patanchon and Avinash Anand},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.39581},
eprint = {2610.39581},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.39581},
}
RIS
TY - UNPB
TI - Efficient GLS map-making and assessment of 1/f noise impact for future CMB satellite missions
AU - Marta Monelli
AU - Guillaume Patanchon
AU - Avinash Anand
PY - 2026
DA - 2026/10/29
PB - arXiv
DO - 10.48550/arxiv.2610.39581
UR - https://arxiv.org/abs/2610.39581
AB - Time-correlated 1/f noise is a major source of large-scale contamination in measurements of the polarized cosmic microwave background (CMB) that can limit sensitivity to primordial B modes. Generalized least-squares (GLS) map-making is optimal for stationary Gaussian noise when its covariance is known, but practical implementations rely on approximations to make the computation tractable. In particular, approximating the noise covariance as circulant enables efficient map-making using fast Fourier transform, but introduces artificial edge effects when applied to finite data chunks. We introduce STITCH (STatistical Inpainting for Timeline CHunks), a preprocessing method that mitigates these effects by extending time-ordered data chunks with synthetic samples that make the extended timelines statistically consistent with the assumed circulant noise model. They can then be processed with standard circulant GLS solvers, preserving their computational efficiency while reducing residuals associated with chunk boundaries. Simulations show that STITCH substantially suppresses these edge effects, bringing the resulting noise power spectra close to those obtained when the circulant approximation holds. We demonstrate negligible impact of our procedure on the recovered signal, such as E→ B leakage, even at large scales. Finally, we investigate how temporal 1/f noise projects onto angular scales for representative LiteBIRD-, \emph{Planck}-, and PICO-like scanning strategies. Using noise simulations spanning a range of temporal noise parameters, we fit an effective model to the resulting angular noise spectra. These results provide an empirical characterization of the temporal-to-angular noise mapping for the configurations studied and a basis for further investigations of correlated noise in future CMB polarization missions.
ER -