Subtraction of point sources from interferometric radio images through an algebraic forward modelling scheme

We present a method for subtracting point sources from interferometric radio images via forward modelling of the instrument response and involving an algebraic non-linear minimization. The method is applied to simulated maps of the Murchison Wide-field Array but is generally useful in cases where only image data are available. After source subtraction, the residual maps have no statistical difference to the expected thermal noise distribution at all angular scales, indicating high effectiveness in the subtraction. Simulations indicate that the errors in recovering the source parameters decrease with increasing signal-to-noise ratio, which is consistent with the theoretical measurement errors. In applying the technique to simulated snapshot observations with the Murchison Wide-field Array, we found that all 101 sources present in the simulation were recovered with an average position error of 10 arcsec and an average flux density error of 0.15 per cent. This led to a dynamic range increase of approximately 3 orders of magnitude. Since all the sources were deconvolved jointly, the subtraction was not limited by source sidelobes but by thermal noise. This technique is a promising deconvolution method for upcoming radio arrays with a huge number of elements and a candidate for the difficult task of subtracting foreground sources from observations of the 21-cm neutral hydrogen signal from the epoch of reionization.

[1]  Andrew J. Benson,et al.  Decontamination of cosmological 21-cm maps , 2007, 0712.0497.

[2]  Max Tegmark,et al.  Twenty-one centimeter tomography with foregrounds , 2005 .

[3]  W. D. Cotton,et al.  Pixelization and dynamic range in radio interferometry , 2008, 0809.1447.

[4]  Max Tegmark,et al.  21 cm Tomography with Foregrounds , 2006 .

[5]  L. Knox,et al.  Multifrequency Analysis of 21 Centimeter Fluctuations from the Era of Reionization , 2005 .

[6]  K. Gorski,et al.  HEALPix: A Framework for High-Resolution Discretization and Fast Analysis of Data Distributed on the Sphere , 2004, astro-ph/0409513.

[7]  Abhirup Datta,et al.  BRIGHT SOURCE SUBTRACTION REQUIREMENTS FOR REDSHIFTED 21 cm MEASUREMENTS , 2010 .

[8]  Matias Zaldarriaga,et al.  Will point sources spoil 21-cm tomography? , 2008, 0807.3952.

[9]  Coryn A. L. Bailer-Jones,et al.  The ILIUM forward modelling algorithm for multivariate parameter estimation and its application to derive stellar parameters from Gaia spectrophotometry , 2009, ArXiv.

[10]  Max Tegmark,et al.  A model of diffuse Galactic radio emission from 10 MHz to 100 GHz , 2008, 0802.1525.

[11]  S. Zaroubi,et al.  Foreground simulations for the LOFAR-epoch of reionization experiment , 2008, 0804.1130.

[12]  Judd D. Bowman,et al.  IMPROVING FOREGROUND SUBTRACTION IN STATISTICAL OBSERVATIONS OF 21 cm EMISSION FROM THE EPOCH OF REIONIZATION , 2006 .

[13]  Mervyn J. Lynch,et al.  THE PRECISION ARRAY FOR PROBING THE EPOCH OF RE-IONIZATION: EIGHT STATION RESULTS , 2009, 0904.2334.

[14]  Saleem Zaroubi,et al.  Non-parametric foreground subtraction for 21-cm epoch of reionization experiments , 2009 .

[15]  Ue-Li Pen,et al.  The GMRT EoR Experiment: Limits on Polarized Sky Brightness at 150 MHz , 2008, 0807.1056.

[16]  S. Zaroubi,et al.  Foregrounds for observations of the cosmological 21 cm line - I. First Westerbork measurements of Galactic emission at 150 MHz in a low latitude field , 2009, 0904.0404.

[17]  C. Carilli,et al.  Synthesis Imaging in Radio Astronomy II , 1999 .

[18]  W. A. Coles,et al.  Interferometric Imaging with the 32 Element Murchison Wide-Field Array , 2010, 1010.1733.

[19]  Miguel F. Morales,et al.  Toward Epoch of Reionization Measurements with Wide-Field Radio Observations , 2003 .

[20]  Matias Zaldarriaga,et al.  Cosmological Parameter Estimation Using 21 cm Radiation from the Epoch of Reionization , 2005, astro-ph/0512263.

[21]  Matias Zaldarriaga,et al.  An improved method for 21‐cm foreground removal , 2009, 0903.4890.

[22]  Garching,et al.  Detection and extraction of signals from the epoch of reionization using higher-order one-point statistics , 2008, 0809.2428.

[23]  A. Bruyn,et al.  High dynamic range mapping of strong radio sources, with application to 3C84 , 1982, Nature.

[24]  S. Bhatnagar,et al.  DETECTION OF SIGNALS FROM COSMIC REIONIZATION USING RADIO INTERFEROMETRIC SIGNAL PROCESSING , 2009, 0908.2639.

[25]  Judd D. Bowman,et al.  FOREGROUND CONTAMINATION IN INTERFEROMETRIC MEASUREMENTS OF THE REDSHIFTED 21 cm POWER SPECTRUM , 2008, 0807.3956.

[26]  Chris Biemesderfer,et al.  Astronomical Data Analysis Software and Systems X , 2001 .

[27]  J. Chengalur,et al.  Foregrounds for redshifted 21-cm studies of reionization: Giant Meter Wave Radio Telescope 153-MHz observations , 2008, 0801.2424.