Search

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year., 19 pages, 11 figures

The radio-wavelength detection of extensive air showers (EAS) initiated by cosmic-ray interactions in the Earth's atmosphere is a promising technique for investigating the origin of these particles and the physics of their interactions. The Low Frequency Array (LOFAR) and the Owens Valley Long Wavelength Array (OVRO-LWA) have both demonstrated that the dense cores of low frequency radio telescope arrays yield detailed information on the radiation ground pattern, which can be used to reconstruct key EAS properties and infer the primary cosmic-ray composition. Here, we demonstrate a new observation mode of the Murchison Widefield Array (MWA), tailored to the observation of the sub-microsecond coherent bursts of radiation produced by EAS. We first show how an aggregate 30.72 MHz bandwidth (3072x 10 kHz frequency channels) recorded at 0.1 ms resolution with the MWA's voltage capture system (VCS) can be synthesised back to the full bandwidth Nyquist resolution of 16.3 ns. This process, which involves `inverting' two sets of polyphase filterbanks, retains 90.5% of the signal-to-noise of a cosmic ray signal. We then demonstrate the timing and positional accuracy of this mode by resolving the location of a calibrator pulse to within 5 m. Finally, preliminary observations show that the rate of nanosecond radio-frequency interference (RFI) events is 0.1 Hz, much lower than that found at the sites of other radio telescopes that study cosmic rays. We conclude that the identification of cosmic rays at the MWA, and hence with the low-frequency component of the Square Kilometre Array, is feasible with minimal loss of efficiency due to RFI., 9 pages, 8 figures. To be published in the Journal of Astronomical Instrumentation

We report the discovery of the first new pulsar with the Murchison Widefield Array (MWA), PSR J0036$-$1033, a long-period (0.9 s) nonrecycled pulsar with a dispersion measure (DM) of 23.1 ${\rm pc\,cm^{-3}}$. It was found after processing only a small fraction ($\sim$1%) of data from an ongoing all-sky pulsar survey. Follow-up observations have been made with the MWA, the upgraded Giant Metrewave Radio Telescope (uGMRT), and the Parkes 64 m telescopes, spanning a frequency range from $\sim$150 MHz to 4 GHz. The pulsar is faint, with an estimated flux density ($S$) of $\sim$1 mJy at 400 MHz and a spectrum $S(\nu)\,\propto\,\nu^{-2.0 \pm 0.2}$, where $\nu$ is frequency. The DM-derived distance implies that it is also a low-luminosity source ($\sim$ 0.1 ${\rm mJy\,kpc^2}$ at 1400 MHz). The analysis of archival MWA observations reveals that the pulsar's mean flux density varies by up to a factor of $\sim$5-6 on timescales of several weeks to months. By combining MWA and uGMRT data, the pulsar position was determined to arcsecond precision. We also report on polarization properties detected in the MWA and Parkes bands. The pulsar's nondetection in previous pulsar and continuum imaging surveys, the observed high variability, and its detection in a small fraction of the survey data searched to date, all hint at a larger population of pulsars that await discovery in the southern hemisphere, with the MWA and the future low-frequency Square Kilometre Array., Comment: 12 pages, 6 figures, 2 tables, Accepted for publication in The ApJ Letters

One of the major challenges for pulsar timing array (PTA) experiments is the mitigation of the effects of the turbulent interstellar medium (ISM) from timing data. These can potentially lead to measurable delays and/or distortions in the pulse profiles and scale strongly with the inverse of the radio frequency. Low-frequency observations are therefore highly appealing for characterizing them. However, in order to achieve the necessary time resolution to resolve profile features of short-period millisecond pulsars, phase-coherent de-dispersion is essential, especially at frequencies below $300$ MHz. We present the lowest-frequency ($80$-$220$ MHz), coherently de-dispersed detections of one of the most promising pulsars for current and future PTAs, PSR J2241$-$5236, using our new beam-former software for the MWA's voltage capture system (VCS), which reconstructs the time series at a much higher time resolution of $\sim 1 \mu$s by re-synthesizing the recorded voltage data at $10$-kHz/$100$-$\mu$s native resolutions. Our data reveal a dual-precursor type feature in the pulse profile that is either faint or absent in high-frequency observations from Parkes. The resultant high-fidelity detections have enabled dispersion measure (DM) determinations with very high precision, of the order of $(2$-$6)\times10^{-6}$ $\rm pc\,cm^{-3}$, owing to the microsecond level timing achievable for this pulsar at the MWA's low frequencies. This underscores the usefulness of low-frequency observations for probing the ISM toward PTA pulsars and informing optimal observing strategies for PTA experiments., Comment: 09 pages, 4 figures, 2 tables

We have observed the relativistic binary pulsar PSR J1141$-$6545 over a period of $\sim$6 years using the Parkes 64 m radio telescope, with a focus on modelling the diffractive intensity scintillations to improve the accuracy of the astrometric timing model. The long-term scintillation, which shows orbital and annual variations, allows us to measure parameters that are difficult to measure with pulsar timing alone. These include: the orbital inclination $i$; the longitude of the ascending node $\Omega$; and the pulsar system transverse velocity. We use the annual variations to resolve the previous ambiguity in the sense of the inclination angle. Using the correct sense, and a prior probability distribution given by a constraint from pulsar timing ($i=73\pm3^\circ$), we find $\Omega=24.8\pm1.8^\circ$ and we estimate the pulsar distance to be $D=10^{+4}_{-3}$ kpc. This then gives us an estimate of this pulsar's proper motion of $\mu_{\alpha}\cos{\delta}=2.9\pm1.0$ mas yr$^{-1}$ in right ascension and $\mu_{\delta}=1.8\pm0.6$ mas yr$^{-1}$ in declination. Finally, we obtain measurements of the spatial structure of the interstellar electron density fluctuations, including: the spatial scale and anisotropy of the diffraction pattern; the distribution of scattering material along the line of sight; and spatial variation in the strength of turbulence from epoch to epoch. We find that the scattering is dominated by a thin screen at a distance of $(0.724\pm0.008)D$, with an anisotropy axial ratio $A_{\rm r} = 2.14\pm0.11$., Comment: 17 pages, 8 figures, 2 tables. Accepted for publication in MNRAS

The high sensitivity and wide frequency coverage of the Murchison Widefield Array allow for the measurement of the spectral scaling of the pulsar scattering timescale, $\alpha$, from a single observation. Here we present three case studies targeted at bright, strongly scattered pulsars J0534+2200 (the Crab pulsar), J0835-4510 (the Vela pulsar) and J0742-2822. We measure the scattering spectral indices to be $-3.8\pm0.2$, $-4.0\pm1.5$, and $-2.5\pm0.6$ for the Crab, Vela, and J0742-2822, respectively. We find that the scattered profiles of both Vela and J0742-2822 are best described by a thin screen model where the Gum Nebula likely contributes most of the observed scattering delay. For the Crab pulsar we see characteristically different pulse shapes compared to higher frequencies, for which none of the scattering screen models we explore are found to be optimal. The presence of a finite inner scale to the turbulence can possibly explain some of the discrepancies., Comment: 14 pages, 8 figures, accepted for publication in ApJ

The Australian Square Kilometre Array Pathfinder (ASKAP) is an interferometric array of 36 fully steerable dishes operating at low radio frequencies (700–1800 MHz) with wide instantaneous bandwidths of 300MHz. Located at the radio-quiet Murchison Radio-astronomy Observatory (MRO) in the deserts of Western Australia, ASKAP is unique in several ways. Each of these 12-metre dishes is equipped with an array of 188 Phased Array Feeds (PAFs) at their focal planes. The PAFs provide enormous flexibility in reliably measuring wide regions of the sky. While conventional telescopes observe the sky with a single instantaneous beam, the ASKAP PAFs, with its current hardware, can electronically form up to 36 simultaneous dual-polarised beams on the sky. The “Roll axis” allows rotation of the dishes about their optic axes, thus ensuring the multiple beams formed by the ASKAP PAFs remain fixed with respect to the sky throughout the observation duration thereby facilitating high dynamic range imaging. While these make ASKAP a very powerful wide field-of-view instrument (with instantaneous $\mathrm {f}-\mathrm {o}- \mathrm {v}\sim 30$square-degrees at the lower end of its operating frequency band), the resulting high data volumes present new challenges in data processing (calibration, imaging, archiving etc.,) and require the algorithms and the processing pipelines to adopt high performance computing solutions.

The Australian Square Kilometre Array Pathfinder (ASKAP) is a new generation low frequency (700-1800 MHz) interferometer capable of observing an instantaneous field of view of 30 square degrees. Such a wide field of view is realized using the Phased Array Feed (PAF) technology that allow forming multiple beams on the sky. The telescope is already in the early science phase with 28 of the 36 antennas routinely available for observations at the time of writing. The hardware has been installed in all antennas and the full array is nearing completion. It is capable of generating about 200 Tb/day of raw data and producing spectral cubes up to 70 Tb in size per observation (if all the data are used; the current plan is to sacrifice the spatial resolution). This fact makes ASKAP an excellent pathfinder and test-bed for the big data challenge which the Square Kilometre Array (SKA) project is facing.

We present a Wide-field ASKAP L-Band Legacy All-sky Blind surveY (WALLABY) study of the nearby ($v_{\rm sys}$ = 915 km s$^{-1}$) spiral galaxy IC 5201 using the Australian Square Kilometre Array Pathfinder (ASKAP). IC 5201 is a blue, barred spiral galaxy that follows the known scaling relations between stellar mass, SFR, HI mass and diameter. We create a four-beam mosaicked HI image cube, from 175 hours of observations made with a 12-antenna sub-array. The RMS noise level of the cube is 1.7 mJy beam$^{-1}$ per channel, equivalent to a column density of $N_{\rm HI}$ = 1.4 $\times$ 10$^{20}$ cm$^{-2}$ over 25 km s$^{-1}$. We report 9 extragalactic HI detections $-$ 5 new HI detections including the first velocity measurements for 2 galaxies. These sources are IC 5201, 3 dwarf satellite galaxies, 2 galaxies and a tidal feature belonging to the NGC 7232/3 triplet and 2 potential infalling galaxies to the triplet. There is evidence of a previous tidal interaction between IC 5201 and the irregular satellite AM 2220$-$460. A close fly-by is likely responsible for the asymmetric optical morphology of IC 5201 and warping its disc, resulting in the irregular morphology of AM 2220$-$460. We quantify the HI kinematics of IC 5201, presenting its rotation curve as well as showing that the warp starts at 14 kpc along the major axis, increasing as a function of radius with a maximum difference in position angle of 20$^\circ$. There is no evidence of stripped HI, triggered or quenched star formation in the system as measured using DECam optical and $GALEX$ UV photometry., Comment: 19 pages, 11 figures, accepted in MNRAS

In this paper we build upon recent work in the radio astronomy community to experimentally demonstrate the viability of passive radar for Space Situational Awareness. Furthermore, we show that the six state parameters of objects in orbit may be measured and used to perform orbit characterisation/estimation., Comment: Published in: 2017 IEEE Radar Conference (RadarConf) URL: http://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=7944483&isnumber=7944108

The Murchison Widefield Array is a low-frequency Square Kilometre Array precursor located at the Murchison Radio-astronomy Observatory in Western Australia. Primarily designed as an imaging telescope, but with a flexible signal path, the capabilities of this telescope have recently been extended to include off-line incoherent and tied-array beam formation using recorded antenna voltages. This has provided the capability for high-time and frequency resolution observations, including a pulsar science program. This paper describes the algorithms and pipeline that we have developed to form the tied array beam products from the summation of calibrated signals of the antenna elements, and presents example polarimetric profiles for PSRs J0437-4715 and J1900-2600 at 185 MHz., Comment: 14 pages, 3 figures, accepted for publication in PASA

Rotating Radio Transients (RRATs) represent a relatively new class of pulsar, primarily characterised by their sporadic bursting emission of single pulses on time scales of minutes to hours. In addition to the difficulty involved in detecting these objects, low-frequency ($, Comment: 18 pages, 8 figures, 5 tables, accepted for publication in PASA

Polarimetric studies of pulsars at low radio frequencies provide important observational insights into the pulsar emission mechanism and beam models, and probe the properties of the magneto-ionic interstellar medium (ISM). Aperture arrays are the main form of next-generation low-frequency telescopes, including the Murchison Widefield Array (MWA). These require a distinctly different approach to data processing (e.g. calibration and beamforming) compared to traditional dish antennas. As the second paper of this series, we present a verification of the MWA's pulsar polarimetry capability, using two bright southern pulsars, PSRs J0742-2822 and J1752-2806. Our observations simultaneously cover multiple frequencies (76-313 MHz) and were taken at multiple zenith angles during a single night for each pulsar. We show that the MWA can be reliably calibrated for zenith angles < 45 degree and frequencies < 270 MHz. We present the polarimetric profiles for PSRs J0742-2822 and J1752-2806 at frequencies lower than 300 MHz for the first time, along with an analysis of the linear polarisation degree and pulse profile evolution with frequency. For PSR J0742-2822, the measured degree of linear polarisation shows a rapid decrease at low frequencies, in contrast with the generally expected trend, which can be attributed to depolarisation effects from small-scale, turbulent, magneto-ionic ISM components. This effect has not been widely explored for pulsars in general, and will be further investigated in future work., accepted for publication in PASA

We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $, Comment: MNRAS Accepted, 17 pages, 10 figures

We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg^2, a resolution of 2.3 arcmin and contains 13,458 sources above 5 sigma. The rms noise in the centre of the image is 4-5 mJy/beam. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30-200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ~3.5 mJy/beam level, due to incompletely-peeled and out-of-image sources, and classical confusion becomes apparent at ~1.7 mJy/beam. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ~150 MHz surveys as it cannot be directly inferred from higher frequency data., 13 pages, 13 figures, accepted for publication in MNRAS

The rare intermittent pulsars pose some of the most challenging questions surrounding the pulsar emission mechanism, but typically have relatively minimal low-frequency ($\lesssim$ 300 MHz) coverage. We present the first low-frequency detection of the intermittent pulsar J1107-5907 with the Murchison Widefield Array (MWA) at 154 MHz and the simultaneous detection from the recently upgraded Molonglo Observatory Synthesis Telescope (UTMOST) at 835 MHz, as part of an on-going observing campaign. During a 30-minute simultaneous observation, we detected the pulsar in its bright emission state for approximately 15 minutes, where 86 and 283 pulses were detected above a signal-to-noise threshold of 6 with the MWA and UTMOST, respectively. Of the detected pulses, 51 had counterparts at both frequencies and exhibited steep spectral indices for both the bright main pulse component and the precursor component. We find that the bright state pulse energy distribution is best parameterised by a log-normal distribution at both frequencies, contrary to previous results which suggested a power law distribution. Further low-frequency observations are required in order to explore in detail aspects such as pulse-to-pulse variability, intensity modulations and to better constrain the signal propagation effects due to the interstellar medium and intermittency characteristics at these frequencies. The spectral index, extended profile emission covering a large fraction of pulse longitude, and the broadband intermittency of PSR J1107-5907 suggests that future low-frequency pulsar searches, for instance those planned with SKA-Low, will be in an excellent position to find and investigate new pulsars of this type., 19 pages, 7 figures, 4 tables, accepted for publication in ApJ; minor updates to keep in-line with published version

TPID5179907. These provide a wide, 30 square degree instantaneous field-of-view by forming up to 36 separate dual-polarisation beams at once. This results in a high data rate: 70 TB of correlated visibilities in an 8 hour observation, requiring high-performance calibration & imaging pipelines. The ASKAPsoft package has been developed specifically for meeting this challenging problem. Particular design decisions in the imaging algorithms, and the framework under which they run, have been made to meet processing requirements on memory and run-time, and scientific requirements.

We present the first low-frequency Murchison Widefield Array (MWA) radio continuum maps of the Magellanic Clouds (MCs), using mosaics from the GaLactic Extragalactic All-Sky MWA (GLEAM) survey. In this paper, we discuss the overall radio continuum morphology between 76 and 227 MHz and compare them with neutral hydrogen maps, 1.4 GHz continuum maps and optical images. Variation of diffuse emission is noticeable across the Large Magellanic Cloud (LMC) but absent across the bar of the Small Magellanic Cloud (SMC). We also measure the integrated flux densities and derive the spectral indices for the MCs. A double power-law model with fixed α1 = -0.1 fit between 19.7 MHz and 8.55 GHz yields α0 = -0.66 ± 0.08 for the LMC. A power-law model yields α _{8.55 GHz}^{85.5 MHz} = -0.82± 0.03 for the SMC. The radio spectral index maps reveal distinctive flat and steep spectral indices for the H II regions and supernova remnants, respectively. We find strong correlation between H II regions and Hα emission. Using a new 150 MHz-Hα relation as a star formation rate indicator, we estimate global star formation rates of 0.068-0.161 M☉ yr-1 and 0.021-0.050 M☉ yr-1 for the LMC and SMC, respectively. Images in 20 frequency bands, and wideband averages are made available via the GLEAM virtual observatory server.

Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field.

Detection of the fluctuations in 21 cm line emission from neutral hydrogen during the Epoch of Reionization in thousand hour integrations poses stringent requirements on calibration and image quality, both of which necessitate accurate primary beam models. The Murchison Widefield Array (MWA) uses phased array antenna elements which maximize collecting area at the cost of complexity. To quantify their performance, we have developed a novel beam measurement system using the 137 MHz ORBCOMM satellite constellation and a reference dipole antenna. Using power ratio measurements, we measure the {\it in situ} beampattern of the MWA antenna tile relative to that of the reference antenna, canceling the variation of satellite flux or polarization with time. We employ angular averaging to mitigate multipath effects (ground scattering), and assess environmental systematics with a null experiment in which the MWA tile is replaced with a second reference dipole. We achieve beam measurements over 30 dB dynamic range in beam sensitivity over a large field of view (65\% of the visible sky), far wider and deeper than drift scans through astronomical sources allow. We verify an analytic model of the MWA tile at this frequency within a few percent statistical scatter within the full width at half maximum. Towards the edges of the main lobe and in the sidelobes, we measure tens of percent systematic deviations. We compare these errors with those expected from known beamforming errors.

Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures., Accepted for publication in Geophysical Research Letters

An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements., Comment: 14 pages, 7 figures

Recent observations with the Murchison Widefield Array at 185~MHz have serendipitously unveiled a heretofore unknown giant and relatively nearby ($z = 0.0178$) radio galaxy associated with NGC\,1534. The diffuse emission presented here is the first indication that NGC\,1534 is one of a rare class of objects (along with NGC\,5128 and NGC\,612) in which a galaxy with a prominent dust lane hosts radio emission on scales of $\sim$700\,kpc. We present details of the radio emission along with a detailed comparison with other radio galaxies with disks. NGC1534 is the lowest surface brightness radio galaxy known with an estimated scaled 1.4-GHz surface brightness of just 0.2\,mJy\,arcmin$^{-2}$. The radio lobes have one of the steepest spectral indices yet observed: $\alpha=-2.1\pm0.1$, and the core to lobe luminosity ratio is $, Comment: 12 pages, 7 figures, accepted MNRAS

Low radio frequency solar observations using the Murchison Widefield Array have recently revealed the presence of numerous weak, short-lived and narrow-band emission features, even during moderately quiet solar conditions. These non-thermal features occur at rates of many thousands per hour in the 30.72 MHz observing bandwidth, and hence, necessarily require an automated approach for their detection and characterization. Here, we employ continuous wavelet transform using a mother Ricker wavelet for feature detection from the dynamic spectrum. We establish the efficacy of this approach and present the first statistically robust characterization of the properties of these features. In particular, we examine distributions of their peak flux densities, spectral spans, temporal spans and peak frequencies. We can reliably detect features weaker than 1 SFU, making them, to the best of our knowledge, the weakest bursts reported in literature. The distribution of their peak flux densities follows a power law with an index of -2.23 in the 12-155 SFU range, implying that they can provide an energetically significant contribution to coronal and chromospheric heating. These features typically last for 1-2 seconds and possess bandwidths of about 4-5 MHz. Their occurrence rate remains fairly flat in the 140-210 MHz frequency range. At the time resolution of the data, they appear as stationary bursts, exhibiting no perceptible frequency drift. These features also appear to ride on a broadband background continuum, hinting at the likelihood of them being weak type-I bursts., 15 pages, 14 figures, Accepted for publication in ApJ

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). 3 h of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1–3 h of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of 104 mK on comoving scales k ≲ 0.5 h Mpc[superscript −1]. This represents the first upper limits on the 21 cm power spectrum fluctuations at redshifts 12 ≲ z ≲ 18 but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window., National Science Foundation (U.S.) (AST-0457585), National Science Foundation (U.S.) (AST-0821321), National Science Foundation (U.S.) (AST-1105835), National Science Foundation (U.S.) (AST-1410719), National Science Foundation (U.S.) (AST-1410484), National Science Foundation (U.S.) (AST-1411622), National Science Foundation (U.S.) (AST-1440343), National Science Foundation (U.S.) (1122374), National Science Foundation (U.S.) (PHY-0835713), National Science Foundation (U.S.) (CAREER-0847753), National Science Foundation (U.S.) (AST-0908884), United States. Air Force Office of Scientific Research (FA9550-0510247)

Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (, 16 pages, 14 figures, accepted to AJ

The Murchison Widefield Array, and its recently developed Voltage Capture System, facilitates extending the low-frequency range of pulsar observations at high-time and -frequency resolution in the Southern Hemisphere, providing further information about pulsars and the ISM. We present the results of an initial time-resolved census of known pulsars using the Murchison Widefield Array. To significantly reduce the processing load, we incoherently sum the detected powers from the 128 Murchison Widefield Array tiles, which yields ~10% of the attainable sensitivity of the coherent sum. This preserves the large field-of-view (~450 deg2 at 185 MHz), allowing multiple pulsars to be observed simultaneously. We developed a WIde-field Pulsar Pipeline that processes the data from each observation and automatically folds every known pulsar located within the beam. We have detected 50 pulsars to date, 6 of which are millisecond pulsars. This is consistent with our expectation, given the telescope sensitivity and the sky coverage of the processed data (~17 000 deg2). For 10 pulsars, we present the lowest frequency detections published. For a subset of the pulsars, we present multi-frequency pulse profiles by combining our data with published profiles from other telescopes. Since the Murchison Widefield Array is a low-frequency precursor to the Square Kilometre Array, we use our census results to forecast that a survey using the low-frequency component of the Square Kilometre Array Phase 1 can potentially detect around 9 400 pulsars.

The phenomenon of subpulse drifting may hold the key to understanding the pulsar emission mechanism. Here, we report on new observations of PSR J0034-0721 (B0031-07), carried out with the Murchison Widefield Array at 185 MHz. We observe three distinct drift modes whose "vertical" drift band separations ($P_3$) and relative abundances are consistent with previous studies at similar and higher frequencies. The driftbands, however, are observed to change their slopes over the course of individual drift modes, which can be interpreted as a continuously changing drift rate. The implied acceleration of the intrinsic carousel rotation cannot easily be explained by plasma models based on ExB drift. Furthermore, we find that methods of classifying the drift modes by means of $P_3$ measurements can sometimes produce erroneous identifications in the presence of a changing drift rate. The "horizontal" separation between driftbands ($P_2$) is found to be larger at later rotation phases within the pulse window, which is inconsistent with the established effects of retardation, aberration, and the motion of the visible point. Longer observations spanning at least ~10,000 pulses are required to determine how the carousel rotation parameters change from one drift sequence to the next., Comment: Accepted for publication in ApJ, 11 pages, 6 figures

We present the first Murchison Widefield Array observations of the well-known cluster of galaxies Abell 3667 (A3667) between 105 and 241 MHz. A3667 is one of the best known examples of a galaxy cluster hosting a double radio relic and has been reported to contain a faint radio halo and bridge. The origins of radio halos, relics and bridges is still unclear, however galaxy cluster mergers seems to be an important factor. We clearly detect the North-West (NW) and South-East (SE) radio relics in A3667 and find an integrated flux density at 149 MHz of 28.1 +/- 1.7 and 2.4 +/- 0.1 Jy, respectively, with an average spectral index, between 120 and 1400 MHz, of -0.9 +/- 0.1 for both relics. We find evidence of a spatial variation in the spectral index across the NW relic steepening towards the centre of the cluster, which indicates an ageing electron population. These properties are consistent with higher frequency observations. We detect emission that could be associated with a radio halo and bridge. How- ever, due to the presence of poorly sampled large-scale Galactic emission and blended point sources we are unable to verify the exact nature of these features., 17 pages, 11 figures

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < �� < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.

We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII regions transition from the optically thin to thick regime we derive the physical properties including the electron density, ionised gas mass and ionising photon flux, towards 61 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date., Comment: 20 pages, 9 figures, 3 tables

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307,455 radio sources with 20 separate flux density measurements across 72--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5\sigma$, which itself is typically 50mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low frequency flux density scale of the southern sky to better than 10%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results., Comment: 30 pages, 18 figures, 6 tables, published in Monthly Notices of the Royal Astronomical Society

The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data., Comment: 17 pages, 10 figures, accepted for publication in The Astrophysical Journal

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2, Comment: 34 pages, accepted for publication in ApJ

We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds., Comment: accepted to ApJ

We estimate spatial gradients in the ionosphere using the Global Positioning System (GPS) and GLONASS (Russian global navigation system) observations, utilising data from multiple GPS stations in the vicinity of Murchison Radio-astronomy Observatory (MRO). In previous work the ionosphere was characterised using a single-station to model the ionosphere as a single layer of fixed height and this was compared with ionospheric data derived from radio astronomy observations obtained from the Murchison Widefield Array (MWA). Having made improvements to our data quality (via cycle slip detection and repair) and incorporating data from the GLONASS system, we now present a multi-station approach. These two developments significantly improve our modelling of the ionosphere. We also explore the effects of a variable-height model. We conclude that modelling the small-scale features in the ionosphere that have been observed with the MWA will require a much denser network of Global Navigation Satellite System (GNSS) stations than is currently available at the MRO., 15 pages, 12 Figures, paper accepted for publication in Publications of the Astronomical Society of Australia (PASA)

Accurate antenna beam models are critical for radio observations aiming to isolate the redshifted 21cm spectral line emission from the Dark Ages and the Epoch of Reionization and unlock the scientific potential of 21cm cosmology. Past work has focused on characterizing mean antenna beam models using either satellite signals or astronomical sources as calibrators, but antenna-to-antenna variation due to imperfect instrumentation has remained unexplored. We characterize this variation for the Murchison Widefield Array (MWA) through laboratory measurements and simulations, finding typical deviations of order +/- 10-20% near the edges of the main lobe and in the sidelobes. We consider the ramifications of these results for image- and power spectrum-based science. In particular, we simulate visibilities measured by a 100m baseline and find that using an otherwise perfect foreground model, unmodeled beamforming errors severely limit foreground subtraction accuracy within the region of Fourier space contaminated by foreground emission (the "wedge"). This region likely contains much of the cosmological signal, and accessing it will require measurement of per-antenna beam patterns. However, unmodeled beamforming errors do not contaminate the Fourier space region expected to be free of foreground contamination (the "EOR window"), showing that foreground avoidance remains a viable strategy., Comment: Accepted to ApJ

Using the final 128 antenna locations of the Murchison Widefield Array (MWA), we calculate its sensitivity to the epoch of reionization (EoR) power spectrum of redshifted 21 cm emission for a fiducial model and provide the tools to calculate the sensitivity for any model. Our calculation takes into account synthesis rotation, chromatic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. For the fiducial model, the MWA will be capable of a 14σ detection of the EoR signal with one full season of observation on two fields (900 and 700 h)., National Science Foundation (U.S.) (Grant AST CAREER-0847753), National Science Foundation (U.S.) (Grant AST-0457585), National Science Foundation (U.S.) (Grant AST-0908884), National Science Foundation (U.S.) (Grant PHY-0835713), United States. Air Force Office of Scientific Research (Grant FA9550-0510247), Smithsonian Astrophysical Observatory, MIT School of Science

Antenna layout is an important design consideration for radio interferometers because it determines the quality of the snapshot point spread function (PSF, or array beam). This is particularly true for experiments targeting the 21 cm Epoch of Reionization signal as the quality of the foreground subtraction depends directly on the spatial dynamic range and thus the smoothness of the baseline distribution. Nearly all sites have constraints on where antennas can be placed---even at the remote Australian location of the MWA (Murchison Widefield Array) there are rock outcrops, flood zones, heritages areas, emergency runways and trees. These exclusion areas can introduce spatial structure into the baseline distribution that enhance the PSF sidelobes and reduce the angular dynamic range. In this paper we present a new method of constrained antenna placement that reduces the spatial structure in the baseline distribution. This method not only outperforms random placement algorithms that avoid exclusion zones, but surprisingly outperforms random placement algorithms without constraints to provide what we believe are the smoothest constrained baseline distributions developed to date. We use our new algorithm to determine antenna placements for the originally planned MWA, and present the antenna locations, baseline distribution, and snapshot PSF for this array choice.

We present a method for subtracting point sources from interferometric radio images via forward modeling of the instrument response and involving an algebraic nonlinear 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%. 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.

The Murchison Widefield Array (MWA) is a next-generation radio telescope currently under construction in the remote Western Australia Outback. Raw data will be generated continuously at 5GiB/s, grouped into 8s cadences. This high throughput motivates the development of on-site, real time processing and reduction in preference to archiving, transport and off-line processing. Each batch of 8s data must be completely reduced before the next batch arrives. Maintaining real time operation will require a sustained performance of around 2.5TFLOP/s (including convolutions, FFTs, interpolations and matrix multiplications). We describe a scalable heterogeneous computing pipeline implementation, exploiting both the high computing density and FLOP-per-Watt ratio of modern GPUs. The architecture is highly parallel within and across nodes, with all major processing elements performed by GPUs. Necessary scatter-gather operations along the pipeline are loosely synchronized between the nodes hosting the GPUs. The MWA will be a frontier scientific instrument and a pathfinder for planned peta- and exascale facilities., Comment: Version accepted by Comp. Phys. Comm

We report on simultaneous wideband observations of Crab giant pulses with the Parkes radio telescope and the Murchison Widefield Array (MWA). The observations were conducted simultaneously at 732 and 3100 MHz with Parkes, and at 120.96, 165.76 and 210.56 MHz with the MWA. Flux density calibration of the MWA data was accomplished using a novel technique based on tied-array beam simulations. We detected between 90-648 giant pulses in the 120.96-210.56 MHz MWA subbands above a $ 5.5\sigma $ threshold while in the Parkes subbands we detected 6344 and 231 giant pulses above a threshold of $ 6\sigma $ at 732 and 3100 MHz, respectively. We show, for the first time over a wide frequency range, that the average spectrum of Crab giant pulses exhibits a significant flattening at low frequencies. The spectral index, $ \alpha $, for giant pulses evolves from a steep, narrow distribution with a mean $ \alpha=-2.6 $ and width $\sigma_\alpha=0.5 $ between 732 and 3100 MHz, to a wide, flat distribution of spectral indices with a mean $ \alpha=-0.7 $ and width $ \sigma_\alpha=1.4 $ between 120.96 and 165.76 MHz. We also comment on the plausibility of giant pulse models for Fast Radio Bursts based on this spectral information., Comment: 21 pages, 12 figures, 8 tables, accepted for publication in ApJ

We present a novel experiment with the capacity to independently measure both the electron density and the magnetic field of the solar corona. We achieve this through measurement of the excess Faraday rotation due to propagation of the polarised emission from a number of pulsars through the magnetic field of the solar corona. This method yields independent measures of the integrated electron density, via dispersion of the pulsed signal and the magnetic field, via the amount of Faraday rotation. In principle this allows the determination of the integrated magnetic field through the solar corona along many lines of sight without any assumptions regarding the electron density distribution. We present a detection of an increase in the rotation measure of the pulsar J1801$-$2304 of approximately 160 \rad at an elongation of 0.95$^\circ$ from the centre of the solar disk. This corresponds to a lower limit of the magnetic field strength along this line of sight of $> 393��\mathrm{G}$. The lack of precision in the integrated electron density measurement restricts this result to a limit, but application of coronal plasma models can further constrain this to approximately 20mG, along a path passing 2.5 solar radii from the solar limb. Which is consistent with predictions obtained using extensions to the Source Surface models published by Wilcox Solar Observatory, 16 pages, 4 figures (1 colour): Submitted to Solar Physics

We present an analysis of the variations seen in the dispersion measures (DMs) of 20 millisecond pulsars observed as part of the Parkes Pulsar Timing Array project. We carry out a statistically rigorous structure function analysis for each pulsar and show that the variations seen for most pulsars are consistent with those expected for an interstellar medium characterised by a Kolmogorov turbulence spectrum. The structure functions for PSRs J1045-4509 and J1909-3744 provide the first clear evidence for a large inner scale, possibly due to ion-neutral damping. We also show the effect of the solar wind on the DMs and show that the simple models presently implemented into pulsar timing packages cannot reliably correct for this effect. For the first time we clearly show how DM variations affect pulsar timing residuals and how they can be corrected in order to obtain the highest possible timing precision. Even with our presently limited data span, the residuals (and all parameters derived from the timing) for six of our pulsars have been significantly improved by correcting for the DM variations., 16 pages, 6 figures, 5 tables, accepted by MNRAS

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (MWA) using a drift scan method that makes efficient use of the MWA's very large field-of-view. We present the observation details, imaging strategies and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40 kHz frequency resolution and 0.5 s time resolution; the second year using 10 kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing and image weighting scheme. At 154 MHz the image resolution is approximately 2.5 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA's primary beam., 17 pages, 7 Figures, 4 tables. Accepted for publication in the Publications of the Astronomical Society of Australia (PASA)

Low-frequency, wide-field radio telescopes such as the Murchison Widefield Array (MWA) enable the dense spatial sampling of the ionosphere and plasmasphere on regional scales. For a physically compact array such as the MWA, the refractive shifts in the positions of celestial sources in the synthesised radio images are proportional to spatial gradients in the total electron content (TEC) transverse to the line of sight. By measuring the angular position shifts of celestial radio sources, one can probe waves and disturbances in the intervening plasma. Radio telescopes differ fundamentally from other techniques for measuring plasma fluctuations in that they are sensitive to TEC gradients/differences rather than absolute TEC. This makes them sensitive specifically to fluctuations about the ambient density, and therefore powerful probes of plasma density waves and irregularities.

The Murchison Widefield Array (MWA) is a Square Kilometre Array (SKA) Precursor. The telescope is located at the Murchison Radio--astronomy Observatory (MRO) in Western Australia (WA). The MWA consists of 4096 dipoles arranged into 128 dual polarisation aperture arrays forming a connected element interferometer that cross-correlates signals from all 256 inputs. A hybrid approach to the correlation task is employed, with some processing stages being performed by bespoke hardware, based on Field Programmable Gate Arrays (FPGAs), and others by Graphics Processing Units (GPUs) housed in general purpose rack mounted servers. The correlation capability required is approximately 8 TFLOPS (Tera FLoating point Operations Per Second). The MWA has commenced operations and the correlator is generating 8.3 TB/day of correlation products, that are subsequently transferred 700 km from the MRO to Perth (WA) in real-time for storage and offline processing. In this paper we outline the correlator design, signal path, and processing elements and present the data format for the internal and external interfaces., Comment: 17 pages, 9 figures. Accepted for publication in PASA. Some figures altered to meet astro-ph submission requirements

The science cases for incorporating high time resolution capabilities into modern radio telescopes are as numerous as they are compelling. Science targets range from exotic sources such as pulsars, to our Sun, to recently detected possible extragalactic bursts of radio emission, the so-called fast radio bursts (FRBs). Originally conceived purely as an imaging telescope, the initial design of the Murchison Widefield Array (MWA) did not include the ability to access high time and frequency resolution voltage data. However, the flexibility of the MWA's software correlator allowed an off-the-shelf solution for adding this capability. This paper describes the system that records the 100 micro-second and 10 kHz resolution voltage data from the MWA. Example science applications, where this capability is critical, are presented, as well as accompanying commissioning results from this mode to demonstrate verification., Comment: 9 pages, 6 figures, 1 table. Accepted for publication in PASA

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively)., Comment: Accepted for publication in PASA