Murchison Widefield Array Research Articles

Image-based searches for pulsar candidates using MWA VCS data

Abstract Pulsars have been studied extensively over the last few decades and have proven instrumental in exploring a wide variety of physics. Discovering more pulsars emitting at low radio frequencies is crucial to further our understanding of spectral properties and emission mechanisms. The Murchison Widefield Array Voltage Capture System (MWA VCS) has been routinely used to study pulsars at low frequencies and discover new pulsars. The MWA VCS offers the unique opportunity of recording complex voltages from all individual antennas (tiles), which can be off-line beamformed or correlated/imaged at millisecond time resolution. Devising imaged-based methods for finding pulsar candidates, which can be verified in beamformed data, can accelerate the complete process and lead to more pulsar detections. Image-based searches for pulsar candidates can reduce the number of tied-array beams required, increasing compute resource efficiency. Despite a factor of $\sim$ 4 loss in sensitivity, searching for pulsar candidates in images from the MWA VCS, we can explore a larger parameter space, potentially leading to discoveries of pulsars missed by high-frequency surveys such as steep spectrum pulsars, exotic binary systems, or pulsars obscured in high-time resolution time series data by propagation effects. Image-based searches are also essential to probing parts of parameter space inaccessible to traditional beamformed searches with the MWA (e.g. at high dispersion measures). In this paper we describe the innovative approach and capability of dual-processing MWA VCS data, that is forming 1-s visibilities and sky images, finding pulsar candidates in these images, and verifying by forming tied-array beam. We developed and tested image-based methods of finding pulsar candidates, which are based on pulsar properties such as steep spectral index, polarisation and variability. The efficiency of these methodologies has been verified on known pulsars, and the main limitations explained in terms of sensitivity and low-frequency spectral turnover of some pulsars. No candidates were confirmed to be a new pulsar, but this new capability will now be applied to a larger subset of observations to accelerate pulsar discoveries with the MWA and potentially speed up future searches with the SKA-Low.

Measuring photometric redshifts for high-redshift radio source surveys

Abstract With the advent of deep, all-sky radio surveys, the need for ancillary data to make the most of the new, high-quality radio data from surveys like the Evolutionary Map of the Universe (EMU), GaLactic and Extragalactic All-sky Murchison Widefield Array survey eXtended, Very Large Array Sky Survey, and LOFAR Two-metre Sky Survey is growing rapidly. Radio surveys produce significant numbers of Active Galactic Nuclei (AGNs) and have a significantly higher average redshift when compared with optical and infrared all-sky surveys. Thus, traditional methods of estimating redshift are challenged, with spectroscopic surveys not reaching the redshift depth of radio surveys, and AGNs making it difficult for template fitting methods to accurately model the source. Machine Learning (ML) methods have been used, but efforts have typically been directed towards optically selected samples, or samples at significantly lower redshift than expected from upcoming radio surveys. This work compiles and homogenises a radio-selected dataset from both the northern hemisphere (making use of Sloan Digital Sky Survey optical photometry) and southern hemisphere (making use of Dark Energy Survey optical photometry). We then test commonly used ML algorithms such as k-Nearest Neighbours (kNN), Random Forest, ANNz, and GPz on this monolithic radio-selected sample. We show that kNN has the lowest percentage of catastrophic outliers, providing the best match for the majority of science cases in the EMU survey. We note that the wider redshift range of the combined dataset used allows for estimation of sources up to $z = 3$ before random scatter begins to dominate. When binning the data into redshift bins and treating the problem as a classification problem, we are able to correctly identify $\approx$ 76% of the highest redshift sources—sources at redshift $z > 2.51$ —as being in either the highest bin ( $z > 2.51$ ) or second highest ( $z = 2.25$ ).

The tracking tapered gridded estimator for the power spectrum from drift scan observations

ABSTRACT Intensity mapping with the redshifted 21-cm line is an emerging tool in cosmology. Drift scan observations, where the antennas are fixed to the ground and the telescope’s pointing centre (PC) changes continuously on the sky due to earth’s rotation, provide broad sky coverage and sustained instrumental stability needed for 21-cm intensity mapping. Here, we present the Tracking Tapered Grided Estimator (TTGE) to quantify the power spectrum of the sky signal estimated directly from the visibilities measured in drift scan radio interferometric observations. The TTGE uses the data from the different PC to estimate the power spectrum of the signal from a small angular region located around a fixed tracking centre (TC). The size of this angular region is decided by a suitably chosen tapering window function that serves to reduce the foreground contamination from bright sources located at large angles from the TC. It is possible to cover the angular footprint of the drift scan observations using multiple TC, and combine the estimated power spectra to increase the signal-to-noise ratio. Here, we have validated the TTGE using simulations of $154 \, {\rm MHz}$ Murchison Wide-field Array drift scan observations. We show that the TTGE can recover the input model angular power spectrum Cℓ within $20 {{\ \rm per\ cent}}$ accuracy over the ℓ range 40 < ℓ < 700.

MeerKAT follow-up of enigmatic GLEAM 4-Jy (G4Jy) sources

ABSTRACT We present the results from studying 140 radio sources in the GLEAM (GaLactic and Extragalactic All-sky MWA [Murchison Widefield Array]) 4-Jy (G4Jy) Sample. These sources were followed-up with MeerKAT to assess their radio morphology and enable host-galaxy identification, as existing radio images of 25 to 45-arcsec resolution do not provide sufficient information. We refer to these sources as the MeerKAT-2019 subset. The aim is to identify the host galaxy of these sources by visually inspecting the overlays comprising radio data from four surveys (at 150, 200, 843/1400, and 1300 MHz). Our morphological classification and host-galaxy identification relies upon the ∼7-arcsec resolution images from MeerKAT (1300 MHz). Through the visual inspection of the overlays, 14 radio sources in the MeerKAT-2019 subset have wide-angle tail morphology, 10 are head-tail, and 5 have X-, S-/Z-shaped morphology. Most of the remaining sources have the radio morphology of typical symmetric lobes. Of 140 sources, we find host galaxies for 98 sources, leaving 42 with no identified host galaxy. These 42 sources still have ambiguous identification even with higher resolution images from MeerKAT.

A targeted search for repeating fast radio bursts with the MWA

ABSTRACT We present a targeted search for low-frequency (144–215 MHz) fast radio burst (FRB) emission from five repeating FRBs using 23.3 h of archival data taken with the Murchison Widefield Array (MWA) voltage capture system (VCS) between 2014 September and 2020 May. This is the first time that the MWA VCS has been used to search for FRB signals from known repeaters, which enables much more sensitive FRB searches than previously performed with the standard MWA correlator mode. We performed a standard single-pulse search with a temporal and spectral resolution of $400\, \mu$s and 10 kHz, respectively, over a $100\, \text{pc}\, \text{cm}^{-3}$ dispersion measure (DM) range centred at the known DM of each studied repeating FRB. No FRBs exceeding a 6σ threshold were detected. The fluence upper limits in the range of 32–1175 and 36–488 Jy ms derived from 10 observations of FRB 20190711A and four observations of FRB 20201124A, respectively, allow us to constrain the spectral indices of their bursts to ≳−1 if these two repeaters were active during the MWA observations. If free–free absorption is responsible for our non-detection, we can constrain the size of the absorbing medium in terms of the electron temperature T to ${\lt} 1.00\times (T/10^4\,\text{K})^{-1.35}\, \text{pc}$, ${\lt} 0.92\times (T/10^4\,\text{K})^{-1.35}\, \text{pc}$, and ${\lt} [0.22\!-\!2.50]\times (T/10^4\,\text{K})^{-1.35}\, \text{pc}$ for FRB 20190117A, FRB 20190711A, and FRB 20201124A, respectively. However, given that the activities of these repeaters are not well characterized, our non-detections could also suggest they were inactive during the MWA observations.

Constraining the Nature of the 18 min Periodic Radio Transient GLEAM-X J162759.5-523504.3 via Multiwavelength Observations and Magneto-thermal Simulations

We observed the periodic radio transient GLEAM-X J162759.5-523504.3 (GLEAM-X J1627) using the Chandra X-ray Observatory for about 30 ks on 2022 January 22–23, simultaneously with radio observations from the Murchison Widefield Array, MeerKAT, and the Australia Telescope Compact Array. Its radio emission and 18 min periodicity led the source to be tentatively interpreted as an extreme magnetar or a peculiar highly magnetic white dwarf. The source was not detected in the 0.3–8 keV energy range with a 3σ upper limit on the count rate of 3 × 10−4 counts s−1. No radio emission was detected during our X-ray observations either. Furthermore, we studied the field around GLEAM-X J1627 using archival European Southern Observatory and DECam Plane Survey data, as well as recent Southern African Large Telescope observations. Many sources are present close to the position of GLEAM-X J1627, but only two within the 2″ radio position uncertainty. Depending on the assumed spectral distribution, the upper limits converted to an X-ray luminosity of L X < 6.5 × 1029 erg s−1 for a blackbody with temperature kT = 0.3 keV, or L X < 9 × 1029 erg s−1 for a power law with photon index Γ = 2 (assuming a 1.3 kpc distance). Furthermore, we performed magneto-thermal simulations for neutron stars considering crust- and core-dominated field configurations. Based on our multiband limits, we conclude that (i) in the magnetar scenario, the X-ray upper limits suggest that GLEAM-X J1627 should be older than ∼1 Myr, unless it has a core-dominated magnetic field or has experienced fast cooling; (ii) in the white dwarf scenario, we can rule out most binary systems, a hot sub-dwarf, and a hot magnetic isolated white dwarf (T ≳ 10.000 K), while a cold isolated white dwarf is still compatible with our limits.

Detection of Ubiquitous Weak and Impulsive Nonthermal Emissions from the Solar Corona

A ubiquitous presence of weak energy releases is one of the most promising hypotheses to explain coronal heating, referred to as the nanoflare hypothesis. The accelerated electrons associated with such weak heating events are also expected to give rise to coherent impulsive emission via plasma instabilities in the meterwave radio band, making this a promising spectral window to look for their presence. Recently Mondal et al. reported the presence of weak and impulsive emissions from quiet Sun regions which seem to meet the requirements of being radio counterparts of the hypothesized nanoflares. Detection of such low-contrast weak emission from the quiet Sun is challenging and, given their implications, it is important to confirm their presence. In this work, using data from the Murchison Widefield Array, we explore the use of an independent robust approach for their detection by separating the dominant, slowly varying component of emission from the weak impulsive one in the visibility domain. We detect milli-Solar Flux Unit-level bursts taking place all over the Sun and characterize their brightness temperatures, distributions, morphologies, durations, and associations with features seen in extreme-UV images. We also attempt to constrain the energies of the nonthermal particles using inputs from the FORWARD coronal model along with some reasonable assumptions, and find them to lie in the subpico flare (∼1019–1021 erg) range. In the process, we also discover perhaps the weakest type III radio burst and another that shows clear signatures of the weakest quasi-periodic pulsations.

Multi-frequency angular power spectrum of the 21 cm signal from the Epoch of Reionisation using the Murchison Widefield Array

Context. The Multi-frequency Angular Power Spectrum (MAPS) is an alternative to spherically averaged power spectra, and computes local fluctuations in the angular power spectrum without need for line-of-sight spectral transform. Aims. We aimed to test different approaches to MAPS and treatment of the foreground contamination, and compare with the spherically averaged power spectrum, and the single-frequency angular power spectrum. Methods. We applied the MAPS to 110 h of data in z = 6.2 − 7.5 obtained for the Murchison Widefield Array Epoch of Reionisation experiment to compute the statistical power of 21 cm brightness temperature fluctuations. In the presence of bright foregrounds, a filter was applied to remove large-scale modes prior to MAPS application, significantly reducing MAPS power due to systematics. Results. The MAPS showed a contrast of 102–103 to a simulated 21 cm cosmological signal for spectral separations of 0−4 MHz after application of the filter, reflecting results for the spherically averaged power spectrum. The single-frequency angular power spectrum was also computed. At z = 7.5 and l = 200, we found an angular power of 53 mK2, exceeding a simulated cosmological signal power by a factor of one thousand. Residual spectral structure, inherent to the calibrated data, and not spectral leakage from large-scale modes, was the dominant source of systematic power bias. The single-frequency angular power spectrum yielded slightly poorer results compared with the spherically averaged power spectrum, having applied a spectral filter to reduce foregrounds. Exploration of other filters may improve this result, along with consideration of wider bandwidths.

Diffuse non-thermal emission in the disks of the Magellanic Clouds

Context. The Magellanic Clouds, two dwarf galaxy companions to the Milky Way, are among the Fermi Large Area Telescope (LAT) brightest γ-ray sources. Aims. We present comprehensive modeling of the non-thermal electromagnetic and neutrino emission in both Clouds. Methods. We self-consistently model the radio and γ-ray spectral energy distribution from their disks based on recently published Murchison Widefield Array and Fermi/LAT data. All relevant radiative processes involving relativistic and thermal electrons (synchrotron, Compton scattering, and bremsstrahlung) and relativistic protons (neutral pion decay following interaction with thermal protons) are considered, using exact emission formulae. Results. Joint spectral analyses indicate that radio emission in the Clouds has both primary and secondary electron synchrotron and thermal bremsstrahlung origin, whereas γ rays originate mostly from π0 decay with some contributions from relativistic bremsstrahlung and Comptonized starlight. The proton spectra in both galaxies are modeled as power laws in energy with similar spectral indices, ~2.4, and energy densities, ~1 eV cm−3. The predicted 0.1–10 GeV neutrino flux is too low for detection by current and upcoming experiments. Conclusions. We confirm earlier suggestions of a largely hadronic origin of the γ-ray emission in both Magellanic Clouds.

Collimation of the kiloparsec-scale radio jets in NGC 2663

ABSTRACT We present the discovery of highly collimated radio jets spanning a total of 355 kpc around the nearby elliptical galaxy NGC 2663, and the possible first detection of recollimation on kiloparsec scales. The small distance to the galaxy (∼28.5 Mpc) allows us to resolve portions of the jets to examine their structure. We combine multiwavelength data: radio observations by the Murchison Widefield Array (MWA), the Australian Square Kilometre Array Pathfinder (ASKAP) and the Australia Telescope Compact Array (ATCA), and X-ray data from Chandra, Swift, and SRG/eROSITA. We present intensity, rotation measure, polarization, spectral index, and X-ray environment maps. Regions of the southern jet show simultaneous narrowing and brightening, which can be interpreted as a signature of the recollimation of the jet by external, environmental pressure, though it is also consistent with intermittent active galactic nuclei or complex internal jet structure. X-ray data suggest that the environment is extremely poor; if the jet is indeed recollimating, the large recollimation scale (40 kpc) is consistent with a slow jet in a low-density environment.

Generating extremely large-volume reionization simulations

ABSTRACT Preparing for a first detection of the 21-cm signal during reionization by large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. Validating that these do not erroneously add/remove features mimicking the signal (e.g. from side lobes or large-scale power leakage) requires simulations extending beyond the primary field of view. However, the Murchison Wide Field Array (MWA) with a field of view of ∼252 deg2 would require simulations spanning several Gpcs, which are currently infeasible. To address this, we developed a simplified version of the seminumerical reionization simulation code 21cmfast, sacrificing some physical accuracy (linear structure formation) in favour of extremely large volumes. We then constructed a 7.5 Gpc comoving volume specifically tailored to the binned spectral resolution of the MWA (∼1.17 cMpc), required for validating the pipeline used in the 2020 MWA 21-cm power spectrum (PS) upper limits. With this large-volume simulation, we then explored: (i) whether smaller volume simulations are biased by missing large-scale modes, (ii) non-Gaussianity in the cosmic variance uncertainty, (iii) biases in the recovered 21-cm PS following foreground wedge avoidance, and (iv) the impact of tiling smaller simulations to achieve large volumes. We found (i) no biases from missing large-scale power, (ii) significant contribution from non-Gaussianity, as expected, (iii) a 10–20 per cent overestimate of the 21-cm PS following wedge mode excision, and (iv) tiling smaller simulations underestimates the large-scale power and cosmic variance.

Working Principle of the Calibration Algorithm for High Dynamic Range Solar Imaging with the Square Kilometre Array Precursor

Imaging the low-frequency radio Sun is an intrinsically challenging problem. Meter-wavelength solar emission spans angular scales from a few arcminutes to a few degrees. These emissions show temporal and spectral variability on sub-second and sub-MHz scales. The brightness temperature of these emissions also varies by many orders of magnitude, which requires high-dynamic-range spectroscopic snapshot imaging. With the unique array configuration of the Murchison Widefield Array (MWA) and the robust calibration and imaging pipeline, Automated Imaging Routine for the Compact Arrays for the Radio Sun (AIRCARS) produces the best spectroscopic snapshot solar images available to date. The working principle and the strength of this algorithm are demonstrated using statistical analysis and simulation. AIRCARS uses the partial phase stability of the MWA, which has a compact core with many antenna elements distributed over a small array footprint. The strength of this algorithm makes it a state-of-the-art calibration and imaging pipeline for low-frequency solar imaging, which is expected to be highly suitable for the upcoming Square Kilometre Array (SKA) and other future radio interferometers for producing high-dynamic-range and high-fidelity images of the Sun.

First measurement of interplanetary scintillation with the ASKAP radio telescope: Implications for space weather

We report on a measurement of interplanetary scintillation (IPS) using the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope. Although this proof-of-concept observation utilised just 3 s of data on a single source, this is nonetheless a significant result, since the exceptional wide field of view of ASKAP, and this validation of its ability to observe within 10 degrees of the Sun, mean that ASKAP has the potential to observe an interplanetary coronal mass ejection (CME) after it has expanded beyond the field of view of white light coronagraphs, but long before it has reached the Earth. We describe our proof of concept observation and extrapolate from the measured noise parameters to determine what information could be gleaned from a longer observation using the full field of view. We demonstrate that, by adopting a ‘Target Of Opportunity’ (TOO) approach, where the telescope is triggered by the detection of a CME in white-light coronagraphs, the majority of interplanetary CMEs could be observed by ASKAP while in an elongation range <30°. It is therefore highly complementary to the colocated Murchison Widefield Array, a lower-frequency instrument which is better suited to observing at elongations >20°.

Understanding Weak Impulsive Narrowband Quiet Sun Emissions (WINQSEs)

AbstractThe confluence of data from the Murchison Widefield Array and an imaging pipeline tailored for spectroscopic snapshot images of the Sun at low radio frequencies have led to enormous improvements in the imaging quality of the Sun. These developments have lowered the detection thresholds by up to two orders of magnitude as compared to earlier studies, and have enabled the discovery of Weak Impulsive Narrowband Quiet Sun Emissions (WINQSEs). Their spatial distribution and various other properties are consistent with being the radio signatures of coronal nanoflares hypothesized by Parker (1988) to explain coronal heating in the quiet Sun emissions. We present the status of the multiple projects we have been pursuing to improve the detection and characterisation of WINQSEs, ranging from looking for them in multiple independent datasets using independent detection techniques to looking for their counter parts to estimate the energy associated with them and understanding their morphologies.

System equivalent flux density of Stokes I, Q, U, and V of a polarimetric interferometer

Aims. We present the system equivalent flux density (SEFD) expressions for all four Stokes parameters: I, Q, U, and V. Methods. The expressions were derived based on our derivation of SEFD I (for Stokes I) and subsequent extensions of that work to phased array and multipole interferometers. The key to the derivation of the SEFD Q, U, and V expressions is to recognize that the noisy estimates of Q, U, and V can be written as the trace of a matrix product. This shows that the SEFD I is a special case, where the general case involves a diagonal or anti-diagonal 2 × 2 matrix interposed in the matrix multiplication. Following this step, the relation between the SEFD for I as well as Q, U, and V immediately becomes evident. Results. We present example calculations for a crossed dipole based on the formulas derived and the comparison between simulation and observation using the Murchison Widefield Array (MWA).

Puzzling large-scale polarization in the galaxy cluster Abell 523

ABSTRACT Large-scale magnetic fields reveal themselves through diffuse synchrotron sources observed in galaxy clusters such as radio halos. Total intensity filaments of these sources have been observed in polarization as well, but only in three radio halos out of about 100 currently known. In this paper we analyse new polarimetric Very Large Array data of the diffuse emission in the galaxy cluster Abell 523 in the frequency range 1–2 GHz. We find for the first time evidence of polarized emission on scales of ∼2.5 Mpc. The total intensity emission is observed only in the central part of the source, likely due to observational limitations. To look for total intensity emission beyond the central region, we combine these data with single-dish observations from the Sardinia Radio Telescope and we compare them with multifrequency total intensity observations obtained with different instruments, including the LOw Frequency ARray and the Murchison Widefield Array. By analysing the rotation measure properties of the system and utilizing numerical simulations, we infer that this polarized emission is associated with filaments of the radio halo located in the outskirts of the system, in the peripheral region closest to the observer.

Independent Discovery of a Nulling Pulsar with Unusual Subpulse Drifting Properties with the Murchison Widefield Array

We report the independent discovery of PSR J0026-1955 with the Murchison Widefield Array (MWA) in the ongoing Southern-sky MWA Rapid Two-metre pulsar survey. J0026-1955 has a period of ∼1.306 s, a dispersion measure of ∼20.869 pc cm−3, and a nulling fraction of ∼77%. This pulsar highlights the advantages of the survey's long dwell times (∼80 minutes), which, when fully searched, will be sensitive to the expected population of similarly bright, intermittent pulsars with long nulls. A single-pulse analysis in the MWA's 140–170 MHz band also reveals a complex subpulse drifting behavior, including both rapid changes of the drift rate characteristic of mode switching pulsars, as well as a slow, consistent evolution of the drift rate within modes. In some longer drift sequences, interruptions in the otherwise smooth drift rate evolution occur preferentially at a particular phase, typically lasting a few pulses. These properties make this pulsar an ideal test bed for prevailing models of drifting behavior such as the carousel model.

Tackling the Unique Challenges of Low-frequency Solar Polarimetry with the Square Kilometre Array Low Precursor: The Algorithm

Coronal magnetic fields are well known to be one of the crucial parameters defining coronal physics and space weather. However, measuring the global coronal magnetic fields remains challenging. The polarization properties of coronal radio emissions are sensitive to coronal magnetic fields. While they can prove to be useful probes of coronal and heliospheric magnetic fields, their usage has been limited by technical and algorithmic challenges. We present a robust algorithm for precise polarization calibration and imaging of low-radio frequency solar observations and demonstrate it on data from the Murchison Widefield Array, a Square Kilometre Array (SKA) precursor. This algorithm is based on the Measurement Equation framework, which forms the basis of all modern radio interferometric calibration and imaging. It delivers high-dynamic-range and high-fidelity full-Stokes solar radio images with instrumental polarization leakages <1%, on par with general astronomical radio imaging, and represents the state of the art. Opening up this rewarding, yet unexplored, phase space will enable multiple novel science investigations and offer considerable discovery potential. Examples include detection of low-level circular polarization from thermal coronal emission to estimate large-scale quiescent coronal fields; polarization of faint gyrosynchrotron emissions from coronal mass ejections for robust estimation of plasma parameters; and detection of the first-ever linear polarization at these frequencies. This method has been developed with the SKA in mind and will enable a new era of high-fidelity spectropolarimetric snapshot solar imaging at low radio frequencies.

High time resolution search for prompt radio emission from the long GRB 210419A with the Murchison Widefield Array

ABSTRACT We present a low-frequency (170–200 MHz) search for prompt radio emission associated with the long GRB 210419A using the rapid-response mode of the Murchison Widefield Array (MWA), triggering observations with the Voltage Capture System for the first time. The MWA began observing GRB 210419A within 89 s of its detection by Swift, enabling us to capture any dispersion delayed signal emitted by this gamma-ray burst (GRB) for a typical range of redshifts. We conducted a standard single pulse search with a temporal and spectral resolution of $100\, \mu$s and 10 kHz over a broad range of dispersion measures from 1 to $5000\, \text{pc}\, \text{cm}^{-3}$, but none were detected. However, fluence upper limits of 77–224 Jy ms derived over a pulse width of 0.5–10 ms and a redshift of 0.6 &amp;lt; z &amp;lt; 4 are some of the most stringent at low radio frequencies. We compared these fluence limits to the GRB jet–interstellar medium interaction model, placing constraints on the fraction of magnetic energy (ϵB ≲ [0.05–0.1]). We also searched for signals during the X-ray flaring activity of GRB 210419A on minute time-scales in the image domain and found no emission, resulting in an intensity upper limit of $0.57\, \text{Jy}\, \text{beam}^{-1}$, corresponding to a constraint of ϵB ≲ 10−3. Our non-detection could imply that GRB 210419A was at a high redshift, there was not enough magnetic energy for low-frequency emission, or the radio waves did not escape from the GRB environment.

Establishing the Capabilities of the Murchison Widefield Array as a Passive Radar for the Surveillance of Space

This paper describes the use of the Murchison Widefield Array, a low-frequency radio telescope at a radio-quiet Western Australian site, as a radar receiver forming part of a continent-spanning multistatic radar network for the surveillance of space. This paper details the system geometry employed, the orbit-specific radar signal processing, and the orbit determination algorithms necessary to ensure resident space objects are detected, tracked, and propagated. Finally, the paper includes the results processed after a short collection campaign utilising several FM radio transmitters across the country, up to a maximum baseline distance of over 2500 km. The results demonstrate the Murchison Widefield Array is able to provide widefield and persistent coverage of objects in low Earth orbit.