-
Unsupervised generation of high dynamic range solar images: A novel algorithm for self-calibration of interferometry data
Authors:
Surajit Mondal,
Atul Mohan,
Divya Oberoi,
John S. Morgan,
Leonid Benkevitch,
Colin J. Lonsdale,
Meagan Crowley,
Iver H. Cairns
Abstract:
Solar radio emission, especially at metre-wavelengths, is well known to vary over small spectral ($\lesssim$100\,kHz) and temporal ($<1$\,s) spans. It is comparatively recently, with the advent of a new generation of instruments, that it has become possible to capture data with sufficient resolution (temporal, spectral and angular) that one can begin to characterize the solar morphology simultaneo…
▽ More
Solar radio emission, especially at metre-wavelengths, is well known to vary over small spectral ($\lesssim$100\,kHz) and temporal ($<1$\,s) spans. It is comparatively recently, with the advent of a new generation of instruments, that it has become possible to capture data with sufficient resolution (temporal, spectral and angular) that one can begin to characterize the solar morphology simultaneously along the axes of time and frequency. This ability is naturally accompanied by an enormous increase in data volumes and computational burden, a problem which will only become more acute with the next generation of instruments such as the Square Kilometre Array (SKA). The usual approach, which requires manual guidance of the calibration process, is impractical. Here we present the "Automated Imaging Routine for Compact Arrays for the Radio Sun (AIRCARS)", an end-to-end imaging pipeline optimized for solar imaging with arrays with a compact core. We have used AIRCARS so far on data from the Murchison Widefield Array (MWA) Phase-I. The dynamic range of the images is routinely from a few hundred to a few thousand. In the few cases, where we have pushed AIRCARS to its limits, the dynamic range can go as high as $\sim$75000. The images made represent a substantial improvement in the state-of-the-art in terms of imaging fidelity and dynamic range. This has the potential to transform the multi-petabyte MWA solar archive from raw visibilities into science-ready images. AIRCARS can also be tuned to upcoming telescopes like the SKA, making it a very useful tool for the heliophysics community.
△ Less
Submitted 28 February, 2019; v1 submitted 23 February, 2019;
originally announced February 2019.
-
Solar Imaging using Low Frequency Arrays
Authors:
C. Lonsdale,
L. Benkevitch,
I. Cairns,
M. Crowley,
P. Erickson,
M. Knapp,
K. Kozarev,
F. Lind,
P. McCauley,
J. Morgan,
D. Oberoi
Abstract:
Low frequency imaging radio arrays such as MWA, LWA and LOFAR have been recently commissioned, and significantly more advanced and flexible arrays are planned for the near term. These powerful instruments offer new opportunities for direct solar imaging at high time and frequency resolution. They can also probe large volumes of the heliosphere simultaneously, by virtue of very large fields of view…
▽ More
Low frequency imaging radio arrays such as MWA, LWA and LOFAR have been recently commissioned, and significantly more advanced and flexible arrays are planned for the near term. These powerful instruments offer new opportunities for direct solar imaging at high time and frequency resolution. They can also probe large volumes of the heliosphere simultaneously, by virtue of very large fields of view. They allow highly detailed, spatially resolved study of solar and heliospheric radio bursts, which are complemented by heliospheric propagation studies using both background astronomical radio emissions as well as the bursts themselves.
In this paper, the state of the art in such wide field solar and heliospheric radio studies is summarized, including recent results from the Murchison Widefield Array (MWA). The prospects for major advances in observational capabilities in the near future are reviewed, with particular emphasis on the RAPID system developed at Haystack Observatory
△ Less
Submitted 5 February, 2018;
originally announced February 2018.
-
Reconstruction of Static Black Hole Images Using Simple Geometric Forms
Authors:
Leonid Benkevitch,
Kazunori Akiyama,
Rusen Lu,
Shepherd Doeleman,
Vincent Fish
Abstract:
General Relativity predicts that the emission close to a black hole must be lensed by its strong gravitational field, illuminating the last photon orbit. This results in a dark circular area known as the black hole 'shadow'. The Event Horizon Telescope (EHT) is a (sub)mm VLBI network capable of Schwarzschild-radius resolution on Sagittarius A* (or Sgr A*), the 4 million solar mass black hole at th…
▽ More
General Relativity predicts that the emission close to a black hole must be lensed by its strong gravitational field, illuminating the last photon orbit. This results in a dark circular area known as the black hole 'shadow'. The Event Horizon Telescope (EHT) is a (sub)mm VLBI network capable of Schwarzschild-radius resolution on Sagittarius A* (or Sgr A*), the 4 million solar mass black hole at the Galactic Center. The goals of the Sgr A* observations include resolving and measuring the details of its morphology. However, EHT data are sparse in the visibility domain, complicating reliable detailed image reconstruction. Therefore, direct pixel imaging should be complemented by other approaches. Using simulated EHT data from a black hole emission model we consider an approach to Sgr A* image reconstruction based on a simple and computationally efficient analytical model that produces images similar to the synthetic ones. The model consists of an eccentric ring with a brightness gradient and a two-dimensional Gaussian. These elemental forms have closed functional representations in the visibility domain, which lowers the computational overhead of fitting the model to the EHT observations. For model fitting we use a version of the Markov chain Monte-Carlo (MCMC) algorithm based on the Metropolis-Hastings sampler with replica exchange. Over a series of simulations we demonstrate that our model can be used for determining geometric measures of a black hole, thus providing information on the shadow size, linking General Relativity with accretion theory.
△ Less
Submitted 31 August, 2016;
originally announced September 2016.
-
Observing the Sun with the Murchison Widefield Array
Authors:
D. Oberoi,
R. Sharma,
S. Bhatnagar,
C. J. Lonsdale,
L. D. Matthews,
I. H. Cairns,
S. J. Tingay,
L. Benkevitch,
A. Donea,
S. M. White,
G. Bernardi,
J. D. Bowman,
F. Briggs,
R. J. Cappallo,
B. E. Corey,
A. Deshpande,
D. Emrich,
B. M. Gaensler,
R. Goeke,
L. J. Greenhill,
B. J. Hazelton,
M. Johnston-Hollitt,
D. L. Kaplan,
J. C. Kasper,
E. Kratzenberg
, et al. (20 additional authors not shown)
Abstract:
The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged i…
▽ More
The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths.
Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations.
△ Less
Submitted 25 March, 2014;
originally announced March 2014.
-
A Method for Smooth Merging of Electron Density Distributions at the Chromosphere-Corona Boundary
Authors:
Leonid Benkevitch,
Divya Oberoi
Abstract:
The electron number density N_e distributions in solar chromosphere and corona are usually described with models of different nature: exponential for the former and inverse power law for the latter. Moreover, the model functions often have different dimensionality, e.g. the chromospheric distribution may depend solely on solar altitude, while the coronal number density may be a function of both al…
▽ More
The electron number density N_e distributions in solar chromosphere and corona are usually described with models of different nature: exponential for the former and inverse power law for the latter. Moreover, the model functions often have different dimensionality, e.g. the chromospheric distribution may depend solely on solar altitude, while the coronal number density may be a function of both altitude and latitude. For applications which need to consider both chromospheric and coronal models, the chromosphere-corona boundary, where these functions have different values as well as gradients, can lead to numerical problems. We encountered this problem in context of ray tracing through the corona at low radio frequencies, as a part of effort to prepare for the analysis of solar images from new generation radio arrays like the Murchison Widefield Array (MWA), Low Frequency Array (LOFAR) and Long Wavelength Array (LWA). We have developed a solution to this problem by using a {\em patch} function, a thin layer between the chromosphere and the corona which matches the values and gradients of the two regions at their respective interfaces. We describe the method we have developed for defining this patch function to seamlessly "stitch" chromospheric and coronal electron density distributions, and generalize the approach to work for any arbitrary distributions of different dimensionality. We show that the complexity of the patch function is independent of the stitched functions dimensionalities. It always has eight parameters (even four for univariate functions) and they may be found without linear system solution for every point. The developed method can potentially be useful for other applications.
△ Less
Submitted 18 October, 2011; v1 submitted 11 October, 2011;
originally announced October 2011.
-
First Spectroscopic Imaging Observations of the Sun at Low Radio Frequencies with the Murchison Widefield Array Prototype
Authors:
Divya Oberoi,
Lynn D. Matthews,
Iver H. Cairns,
David Emrich,
Vasili Lobzin,
Colin J. Lonsdale,
Edward H. Morgan,
T. Prabu,
Harish Vedantham,
Randall B. Wayth,
Andrew Williams,
Christopher Williams,
Stephen M. White,
G. Allen,
Wayne Arcus,
David Barnes,
Leonid Benkevitch,
Gianni Bernardi,
Judd D. Bowman,
Frank H. Briggs,
John D. Bunton,
Steve Burns,
Roger C. Cappallo,
M. A. Clark,
Brian E. Corey
, et al. (45 additional authors not shown)
Abstract:
We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emis…
▽ More
We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years.
△ Less
Submitted 3 January, 2011;
originally announced January 2011.
-
Interferometric imaging with the 32 element Murchison Wide-field Array
Authors:
S. M. Ord,
D. A. Mitchell,
R. B. Wayth,
L. J. Greenhill,
G. Bernardi,
S. Gleadow,
R. G. Edgar,
M. A. Clark,
G. Allen,
W. Arcus,
L. Benkevitch,
J. D. Bowman,
F. H. Briggs,
J. D. Bunton,
S. Burns,
R. J. Cappallo,
W. A. Coles,
B. E. Corey,
L. deSouza,
S. S. Doeleman,
M. Derome,
A. Deshpande,
D. Emrich,
R. Goeke,
M. R. Gopalakrishna
, et al. (34 additional authors not shown)
Abstract:
The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 M…
▽ More
The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 MHz to 300 MHz with an instantaneous field of view that is tens of degrees wide and a resolution of a few arcminutes. A 32-station prototype of the MWA has been recently commissioned and a set of observations taken that exercise the whole acquisition and processing pipeline. We present Stokes I, Q, and U images from two ~4 hour integrations of a field 20 degrees wide centered on Pictoris A. These images demonstrate the capacity and stability of a real-time calibration and imaging technique employing the weighted addition of warped snapshots to counter extreme wide field imaging distortions.
△ Less
Submitted 8 October, 2010;
originally announced October 2010.
-
The Murchison Widefield Array
Authors:
Daniel A. Mitchell,
Lincoln J. Greenhill,
Stephen M. Ord,
Gianni Bernardi,
Randall B. Wayth,
Richard G. Edgar,
Michael A. Clark,
Kevin Dal,
Hanspeter Pfister,
Stewart J. Gleadow,
W. Arcus,
F. H. Briggs,
L. Benkevitch,
J. D. Bowman,
J. D. Bunton,
S. Burns,
R. J. Cappallo,
B. E. Corey,
A. de Oliveira-Costa,
L. Desouza,
S. S. Doeleman,
M. F. Derome,
D. Emrich,
M. Glossop,
R. Goeke
, et al. (35 additional authors not shown)
Abstract:
It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.
It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.
△ Less
Submitted 15 August, 2010;
originally announced August 2010.
-
Algorithm for Tracing Radio Rays in Solar Corona and Chromosphere
Authors:
Leonid Benkevitch,
Igor Sokolov,
Divya Oberoi,
Thomas Zurbuchen
Abstract:
In this paper a new efficient algorithm for computation of radio wave ray trajectories is described. The algorithm is based on an original second-order difference scheme with a specific "length-conservation" property, which allows to resolve the ray shape even in the regions where its curvature is high. Besides the scheme, the algorithm includes a number of mechanisms ensuring its correctness and…
▽ More
In this paper a new efficient algorithm for computation of radio wave ray trajectories is described. The algorithm is based on an original second-order difference scheme with a specific "length-conservation" property, which allows to resolve the ray shape even in the regions where its curvature is high. Besides the scheme, the algorithm includes a number of mechanisms ensuring its correctness and stability. The algorithm is intended for obtaining multi-pixel solar images in wide range of radio frequencies, and it is designed to be used in studies of the solar environment with modern high-resolution radio interferometers and radio telescopes such as the Murchison Widefield Array
△ Less
Submitted 18 October, 2010; v1 submitted 29 June, 2010;
originally announced June 2010.
-
The Murchison Widefield Array: Design Overview
Authors:
Colin J. Lonsdale,
Roger J. Cappallo,
Miguel F. Morales,
Frank H. Briggs,
Leonid Benkevitch,
Judd D. Bowman,
John D. Bunton,
Steven Burns,
Brian E. Corey,
Ludi deSouza,
Sheperd S. Doeleman,
Mark Derome,
Avinash Deshpande,
M. R. Gopalakrishna,
Lincoln J. Greenhill,
David Herne,
Jacqueline N. Hewitt,
P. A. Kamini,
Justin C. Kasper,
Barton B. Kincaid,
Jonathan Kocz,
Errol Kowald,
Eric Kratzenberg,
Deepak Kumar,
Mervyn J. Lynch
, et al. (23 additional authors not shown)
Abstract:
The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the…
▽ More
The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.
△ Less
Submitted 10 March, 2009;
originally announced March 2009.