Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of ionospheric refraction on radio interferometric data
113
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The Epoch of Reionisation (EoR) is the period within which the neutral universe transitioned to an ionised one. This period remains unobserved using low-frequency radio interferometers which target the 21 cm signal of neutral hydrogen emitted in this era. The Murchison Widefield Array (MWA) radio telescope was built with the detection of this signal as one of its major science goals. One of the most significant challenges towards a successful detection is that of calibration, especially in the presence of the Earths ionosphere. By introducing refractive source shifts, distorting source shapes and scintillating flux densities, the ionosphere is a major nuisance in low-frequency radio astronomy. We introduce SIVIO, a software tool developed for simulating observations of the MWA through different ionospheric conditions estimated using thin screen approximation models and propagated into the visibilities. This enables us to directly assess the impact of the ionosphere on observed EoR data and the resulting power spectra. We show that the simulated data captures the dispersive behaviour of ionospheric effects. We show that the spatial structure of the simulated ionospheric media is accurately reconstructed either from the resultant source positional offsets or from parameters evaluated during the data calibration procedure. In turn, this will inform on the best strategies of identifying and efficiently eliminating ionospheric contamination in EoR data moving into the Square Kilometre Array era.
rate research
Read More
The short-spacing problem describes the inherent inability of radio-interferometric arrays to measure the integrated flux and structure of diffuse emission associated with extended sources. New interferometric arrays, such as SKA, require solutions to efficiently combine interferometer and single-dish data. We present a new and open source approach for merging single-dish and cleaned interferometric data sets requiring a minimum of data manipulation while offering a rigid flux determination and full high angular resolution. Our approach combines single-dish and cleaned interferometric data in the image domain. This approach is tested for both Galactic and extragalactic HI data sets. Furthermore, a quantitative comparison of our results to commonly used methods is provided. Additionally, for the interferometric data sets of NGC4214 and NGC5055, we study the impact of different imaging parameters as well as their influence on the combination for NGC4214. The approach does not require the raw data (visibilities) or any additional special information such as antenna patterns. This is advantageous especially in the light of upcoming radio surveys with heterogeneous antenna designs.
New and upcoming radio interferometers will produce unprecedented amounts of data that demand extremely powerful computers for processing. This is a limiting factor due to the large computational power and energy costs involved. Such limitations restrict several key data processing steps in radio interferometry. One such step is calibration where systematic errors in the data are determined and corrected. Accurate calibration is an essential component in reaching many scientific goals in radio astronomy and the use of consensus optimization that exploits the continuity of systematic errors across frequency significantly improves calibration accuracy. In order to reach full consensus, data at all frequencies need to be calibrated simultaneously. In the SKA regime, this can become intractable if the available compute agents do not have the resources to process data from all frequency channels simultaneously. In this paper, we propose a multiplexing scheme that is based on the alternating direction method of multipliers (ADMM) with cyclic updates. With this scheme, it is possible to simultaneously calibrate the full dataset using far fewer compute agents than the number of frequencies at which data are available. We give simulation results to show the feasibility of the proposed multiplexing scheme in simultaneously calibrating a full dataset when a limited number of compute agents are available.
A tomographic method is described to quantify the three-dimensional power-spectrum of the ionospheric electron-density fluctuations based on radio-interferometric observations by a two-dimensional planar array. The method is valid to first-order Born approximation and might be applicable to correct observed visibilities for phase variations due to the imprint of the full three-dimensional ionosphere. It is shown that not the ionospheric electron density distribution is the primary structure to model in interferometry, but its autocorrelation function or equivalent its power-spectrum. An exact mathematical expression is derived that provides the three dimensional power-spectrum of the ionospheric electron-density fluctuations directly from a rescaled scattered intensity field and an incident intensity field convolved with a complex unit phasor that depends on the w-term and is defined on the full sky pupil plane. In the limit of a small field of view, the method reduces to the single phase screen approximation. Tomographic self-calibration can become important in high-dynamic range observations at low radio frequencies with wide-field antenna interferometers, because a three-dimensional ionosphere causes a spatially varying convolution of the sky, whereas a single phase screen results in a spatially invariant convolution. A thick ionosphere can therefore not be approximated by a single phase screen without introducing errors in the calibration process. By applying a Radon projection and the Fourier projection-slice theorem, it is shown that the phase-screen approach in three dimensions is identical to the tomographic method. Finally we suggest that residual speckle can cause a diffuse intensity halo around sources, due to uncorrectable ionospheric phase fluctuations in the short integrations, which could pose a fundamental limit on the dynamic range in long-integration images.
Calibration of radio interferometric observations becomes increasingly difficult towards lower frequencies. Below ~300 MHz, spatially variant refractions and propagation delays of radio waves traveling through the ionosphere cause phase rotations that can vary significantly with time, viewing direction and antenna location. In this article we present a description and first results of SPAM (Source Peeling and Atmospheric Modeling), a new calibration method that attempts to iteratively solve and correct for ionospheric phase errors. To model the ionosphere, we construct a time-variant, 2-dimensional phase screen at fixed height above the Earths surface. Spatial variations are described by a truncated set of discrete Karhunen-Loeve base functions, optimized for an assumed power-law spectral density of free electrons density fluctuations, and a given configuration of calibrator sources and antenna locations. The model is constrained using antenna-based gain phases from individual self-calibrations on the available bright sources in the field-of-view. Application of SPAM on three test cases, a simulated visibility data set and two selected 74 MHz VLA data sets, yields significant improvements in image background noise (5-75 percent reduction) and source peak fluxes (up to 25 percent increase) as compared to the existing self-calibration and field-based calibration methods, which indicates a significant improvement in ionospheric phase calibration accuracy.
We introduce a new pipeline for analyzing and mitigating radio frequency interference (RFI), which we call Sky-Subtracted Incoherent Noise Spectra (SSINS). SSINS is designed to identify and remove faint RFI below the single baseline thermal noise by employing a frequency-matched detection algorithm on baseline-averaged amplitudes of time-differenced visibilities. We demonstrate the capabilities of SSINS using the Murchison Widefield Array (MWA) in Western Australia. We successfully image aircraft flying over the array via digital television (DTV) reflection detected using SSINS and summarize an RFI occupancy survey of MWA Epoch of Reionization data. We describe how to use SSINS with new data using a documented, publicly available implementation with comprehensive usage tutorials.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
