We have extended our previous work to use the Murchison Widefield Array (MWA) as a non-coherent passive radar system in the FM frequency band, using terrestrial FM transmitters to illuminate objects in Low Earth Orbit LEO) and the MWA as the sensitive receiving element for the radar return. We have implemented a blind detection algorithm that searches for these reflected signals in difference images constructed using standard interferometric imaging techniques. From 20 hours of archived MWA observations, we conduct a survey of LEO, detecting 74 unique objects over multiple passes and demonstrating the MWA to be a valuable addition to the global Space Domain Awareness network. We detected objects with ranges up to 977 km and as small as 0.03 m^2 radar cross section. We found that 30 objects were either non-operational satellites or upper-stage rocket body debris. Additionally, we also detected FM reflections from Geminid meteors and aircraft flying over the MWA. Most of the detections of objects in LEO were found to lie within the parameter space predicted by previous feasibility studies, verifying the performance of the MWA for this application.

The entire southern sky (Declination, $\delta <$ 30 deg) has been observed using the Murchison Widefield Array (MWA), which provides radio imaging of $\sim$2-arcmin resolution at low frequencies (72-231 MHz). This is the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we have previously used a combination of visual inspection, cross-checks against the literature, and internal matching to identify the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy) in the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg). We refer to these 1,863 sources as the GLEAM 4-Jy (G4Jy) Sample, and use radio images (of $\leq$ 45-arcsec resolution), and multi-wavelength information, to assess their morphology and identify the galaxy that is hosting the radio emission (where appropriate). Details of how to access all of the overlays used for this work are available at https://github.com/svw26/G4Jy. Alongside this we conduct further checks against the literature, which we document in this paper for individual sources. Whilst the vast majority of the G4Jy Sample are active galactic nuclei with powerful radio-jets, we highlight that it also contains a nebula, two nearby, star-forming galaxies, a cluster relic, and a cluster halo. There are also three extended sources for which we are unable to infer the mechanism that gives rise to the low-frequency emission. In the G4Jy catalogue we provide mid-infrared identifications for 86% of the sources, and flag the remainder as: having an uncertain identification (129 sources), having a faint/uncharacterised mid-infrared host (126 sources), or it being inappropriate to specify a host (2 sources). For the subset of 129 sources, there is ambiguity concerning candidate host-galaxies, and this includes four sources (B0424$-$728, B0703$-$451, 3C 198, and 3C 403.1) where we question the existing identification.

The Murchison Widefield Array (MWA) has observed the entire southern sky (Declination, $\delta <$ 30 deg) at low radio-frequencies, over the range 72-231 MHz. These observations constitute the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we use the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg) to define the GLEAM 4-Jy (G4Jy) Sample. This is a complete sample of the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy), the majority of which are active galactic nuclei with powerful radio-jets. Crucially, low-frequency observations allow the selection of such sources in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). We then use higher-resolution radio images, and information at other wavelengths, to morphologically classify the brightest components in GLEAM. We also conduct cross-checks against the literature, and perform internal matching, in order to improve sample completeness (which is estimated to be $>$ 95.5%). This results in a catalogue of 1,863 sources, making the G4Jy Sample over 10 times larger than that of the revised Third Cambridge Catalogue of Radio Sources (3CRR; $S_{\mathrm{178MHz}} >$ 10.9 Jy). Of these G4Jy sources, 78 are resolved by the MWA (Phase-I) synthesised beam ($\sim$2 arcmin at 200 MHz), and we label 67% of the sample as 'single', 26% as 'double', 4% as 'triple', and 3% as having 'complex' morphology at $\sim$1 GHz (45-arcsec resolution). Alongside this, our value-added catalogue provides mid-infrared source associations (subject to 6-arcsec resolution at 3.4 micron) for the radio emission, as identified through visual inspection and thorough checks against the literature. As such, the G4Jy Sample can be used as a reliable training set for cross-identification via machine-learning algorithms. [Abstract abridged for arXiv submission.]

The low-frequency linearly-polarised radio source population is largely unexplored. However, a renaissance in low-frequency polarimetry has been enabled by pathfinder and precursor instruments for the Square Kilometre Array. In this second paper from the POlarised GaLactic and Extragalactic All-Sky Murchison Widefield Array (MWA) Survey -- the POlarised GLEAM Survey, or POGS -- we present the results from our all-sky MWA Phase I Faraday Rotation Measure survey. Our survey covers nearly the entire Southern sky in the Declination range $-82^{\circ}$ to $+30^{\circ}$ at a resolution between around three and seven arcminutes (depending on Declination) using data in the frequency range 169$-$231 MHz. We have performed two targeted searches: the first covering 25,489 square degrees of sky, searching for extragalactic polarised sources; the second covering the entire sky South of Declination $+30^{\circ}$, searching for known pulsars. We detect a total of 517 sources with 200 MHz linearly-polarised flux densities between 9.9 mJy and 1.7 Jy, of which 33 are known radio pulsars. All sources in our catalogues have Faraday rotation measures in the range $-328.07$ rad m$^{-2}$ to $+279.62$ rad m$^{-2}$. The Faraday rotation measures are broadly consistent with results from higher-frequency surveys, but with typically more than an order of magnitude improvement in the precision, highlighting the power of low-frequency polarisation surveys to accurately study Galactic and extragalactic magnetic fields. We discuss the properties of our extragalactic and known-pulsar source population, how the sky distribution relates to Galactic features, and identify a handful of new pulsar candidates among our nominally extragalactic source population.

The search for emission from weakly interacting massive particle (WIMP) dark matter annihilation and decay has become a multi-pronged area of research not only targeting a diverse selection of astrophysical objects, but also taking advantage of the entire electromagnetic spectrum. The decay of WIMP particles into standard model particles has been suggested as a possible channel for synchrotron emission to be detected at low radio frequencies. Here, we present the stacking analysis of a sample of 33 dwarf spheroidal (dSph) galaxies with low-frequency (72 - 231 MHz) radio images from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. We produce radial surface brightness profiles of images centred upon each dSph galaxy with background radio sources masked. We remove ten fields from the stacking due to contamination from either poorly subtracted, bright radio sources or strong background gradients across the field. The remaining 23 dSph galaxies are stacked in an attempt to obtain a statistical detection of any WIMP-induced synchrotron emission in these systems. We find that the stacked radial brightness profile does not exhibit a statistically significant detection above the 95% confidence level of $\sim$1.5 mJy beam$^{-1}$. This novel technique shows the potential of using low-frequency radio images to constrain fundamental properties of particle dark matter.

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

The number of active and non active satellites in Earth orbit has dramatically increased in recent decades, requiring the development of novel surveillance techniques to monitor and track them. In this paper, we build upon previous non-coherent passive radar space surveillance demonstrations undertaken using the Murchison Widefield Array (MWA). We develop the concept of the Dynamic Signal to Noise Ratio Spectrum (DSNRS) in order to isolate signals of interest (reflections of FM transmissions of objects in orbit) and efficiently differentiate them from direct path reception events. We detect and track Alouette-2, ALOS, UKube-1, the International Space Station, and Duchifat-1 in this manner. We also identified out-of-band transmissions from Duchifat-1 and UKube-1 using these techniques, demonstrating the MWA's capability to look for spurious transmissions from satellites. We identify an offset from the locations predicted by the cataloged orbital parameters for some of the satellites, demonstrating the potential of using MWA for satellite catalog maintenance. These results demonstrate the capability of the MWA for Space Situational Awareness and we describe future work in this area.

This work makes available a further 2,860deg$^2$ of the GLEAM survey, covering half of the accessible Galactic Plane, across twenty frequency bands sampling $72-231$MHz, with resolution $4'-2'$. Unlike previous GLEAM data releases, we used multi-scale clean to better deconvolve large-scale Galactic structure. For the Galactic longitude ranges $345^\circ < l < 67^\circ$, $180^\circ < l < 240^\circ$, we provide a compact source catalogue of 22,037 components selected from a 60-MHz bandwidth image centred at 200-MHz, with RMS noise $\approx10-20$mJy beam$^{-1}$ and position accuracy better than $2"$. The catalogue has a completeness of 50% at $\approx120$mJy, and a reliability of 99.86%. It covers Galactic latitudes $1^\circ\leq|b|\leq10^\circ$ toward the Galactic Centre and $|b|\leq10^\circ$ for other regions, and is available from Vizier; images covering $|b|\leq10^\circ$ for all longitudes are made available on the GLEAM VO server and SkyView.

We have detected 27 new supernova remnants (SNRs) using a new data release of the GLEAM survey from the Murchison Widefield Array (MWA) telescope, including the lowest surface-brightness SNR ever detected, G0.1-9.7. Our method uses spectral fitting to the radio continuum to derive spectral indices for 26/27 candidates, and our low-frequency observations probe a steeper-spectrum population than previously discovered. None of the candidates have coincident Wide-field Infrared Survey Explorer mid-IR emission, further showing that the emission is non-thermal. Using pulsar associations we derive physical properties for six candidate SNRs, finding G0.1-9.7 may be younger than 10kyr. 60% of the candidates subtend areas larger than 0.2deg$^{2}$ on the sky, compared to $<25$% of previously-detected SNRs. We also make the first detection of two SNRs in the Galactic longitude range $220^\circ-240^\circ$.

We examined the latest data release from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey covering $345^\circ < l < 60^\circ$, $180^\circ < l < 240^\circ$, using these data and that of the Widefield Infrared Survey Explorer to follow up proposed candidate Supernova Remnants from other sources. Of the 101 candidates proposed in the region, we are able to definitively confirm ten as SNRs, tentatively confirm two as SNRs, and reclassify five as Hii regions. A further two are detectable in our images but difficult to classify; the remaining 82 are undetectable in these data. We also investigated the 18 unclassified Multi-Array Galactic Plane Imaging Survey (MAGPIS) candidate SNRs, newly confirming three as SNRs, reclassifying two as Hii regions, and exploring the unusual spectra and morphology of two others.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

The Murchison Widefield Array (MWA) is an electronically steered low frequency ($<300$\,MHz) radio interferometer, with a `slew' time less than 8seconds. Low frequency ($\sim 100$\,MHz) radio telescopes are ideally suited for rapid-response follow-up of transients due to their large field of view, the inverted spectrum of coherent emission, and the fact that the dispersion delay between a 1GHz and 100MHz pulse is on the order of $1-10$\,min for dispersion measures of $100-2000$\,pc/cm$^3$. The MWA has previously been used to provide fast follow up for transient events including gamma-ray bursts, fast radio bursts, and gravitational waves, using systems that respond to gamma-ray coordinates network (GCN) packet-based notifications. We describe a system for automatically triggering MWA observations of such events, based on VOEvent triggers, which is more flexible, capable, and accurate than previous systems. The system can respond to external multi-messenger triggers, which makes it well-suited to searching for prompt coherent radio emission from gamma-ray bursts, the study of fast radio bursts and gravitational waves, single pulse studies of pulsars, and rapid follow-up of high-energy superflares from flare stars. The new triggering system has the capability to trigger observations in both the regular correlator mode (limited to $\geq 0.5$\,s integrations) or using the Voltage Capture System (VCS, $0.1$\,ms integration) of the MWA, and represents a new mode of operation for the MWA. The upgraded standard correlator triggering capability has been in use since MWA observing semester 2018B (July-Dec 2018), and the VCS and buffered mode triggers will become available for observing in a future semester.

We report the discovery of a radio transient in the spiral galaxy M81. The transient was detected in early 2015 as part of a two-year survey of M81 made up of 12 epochs using the Karl G. Jansky Very Large Array. While undetected on 2014 September 12, the source was first detected on 2015 January 2, from which point it remained visible at an approximately constant luminosity of $L_{R} = 1.5 \pm 0.1 \times 10^{24}$erg s$^{-1}$ Hz$^{-1}$ at the observing frequency of 6 GHz for at least 2 months. Assuming this is a synchrotron event with a rise-time between 2.6 and 112 days, the peak luminosity (at equipartition) corresponds to a minimum energy of $10^{44} \lesssim E_{min} \lesssim 10^{46}$erg and jet power of $P_{min} \sim 10^{39}$erg s$^{-1}$, which are higher than most known X-ray binaries. Given its longevity, lack of short-term radio variability, and the absence of any multi-wavelength counterpart (X-ray luminosity $L_{X} \lesssim 10^{36}$erg s$^{-1}$), it does not behave like known Galactic or extragalactic X-ray binaries. The M81 transient radio properties more closely resemble the unidentified radio transient 43.78+59.3 discovered in M82, which has been suggested to be a radio nebula associated with an accreting source similar to SS 433. One possibility is that both the new M81 transient and the M82 transient may be the birth of a short-lived radio bubble associated with a discrete accretion event similar to those observed from the ULX Holmberg II X-1. However, it is not possible to rule out other identifications including long-term supernova shockwave interactions with the surrounding medium from a faint supernova or a background active galaxy.

We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0.18 is accelerating particles beyond TeV energies, thus making this object a promising new cosmic ray hadron source candidate. The remnant cannot be seen in current optical, infrared and X-ray data-sets. We do find, however, a dip in CO-traced molecular gas at a line-of-sight velocity of ~85 km/s, suggesting the existence of a G23.11+0.18 progenitor wind-blown bubble. Furthermore, the discovery of molecular gas clumps at a neighbouring velocity towards HESS J1832-085 adheres to the notion that a hadronic gamma-ray production mechanism is plausible towards the north of the remnant. Based on these morphological arguments, we propose an interstellar medium association for G23.11+0.18 at a kinematic distance of 4.6+/-0.8 kpc.

Surveys for transient and variable phenomena can be confounded by the presence of extrinsic variability such as refractive interstellar scintillation (RISS). We have developed an all-sky model for RISS which can predict variability on a variety of timescales, survey locations, and observing frequencies. The model makes use of Halpha intensity maps to probe the emission measure along the line of sight, convert this to a scattering measure, and finally a scintillation strength. The model uses previously developed and long understood physics along with (indirect) measurements of the electron content and distribution within the Milky Way. We develop a set of expectations that are useful in the planning of future surveys for transient and radio variability, and demonstrate that the 1-GHz sky is a poor predictor of the variable nature of the $100$-MHz sky. Interestingly, the correlation between the incidence of variability and Galactic latitude which has been seen at 1GHz, is reversed at 100MHz. We compare the predictions of our model to a low-frequency radio survey that was conducted with the Murchison Widefield Array, and find good qualitative agreement. We discuss the implications, current limitations, and future development of the model. The model has been implemented in a Python code and is available on GitHub/Zenodo.

In this paper we demonstrate improved techniques to extend coherent processing intervals for passive radar processing, with the Murchison Widefield Array. Specifically, we apply a two stage linear range and Doppler migration compensation by utilising Keystone Formatting and a recent dechirping method. These methods are used to further demonstrate the potential for the surveillance of space with the Murchison Widefield Array using passive radar, by detecting objects orders of magnitude smaller than previous work. This paper also demonstrates how the linear Doppler migration methods can be extended to higher order compensation to further increase potential processing intervals.

We present Murchison Widefield Array (MWA) monitoring of the Kepler K2 mission Fields 3, 4, and 5 at frequencies of 155 and 186 MHz, from observations contemporaneous with the K2 observations. This work follows from previous MWA and GMRT surveys of Field 1, with the current work benefiting from a range of improvements in the data processing and analysis. We continue to build a body of systematic low frequency blind surveys overlapping with transient/variable survey fields at other wavelengths, providing multi-wavelength data for object classes such as flare stars. From the current work, we detect no variable objects at a surface density above 2e-4 per square degree, at flux densities of ~500 mJy, and observation cadence of days to weeks, representing almost an order of magnitude decrease in measured upper limits compared to previous results in this part of observational parameter space. This continues to show that radio transients at metre and centimetre wavelengths are rare.

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.

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$.

We present Robbie: a general work-flow for the detection and characterization of radio variability and transient events in the image domain. Robbie is designed to work in a batch processing paradigm with a modular design so that components can be swapped out or upgraded to adapt to different input data, whilst retaining a consistent and coherent methodological approach. Robbie is based on commonly used and open software, and is encapsulated in a Makefile to aid portability and reproducibility. In this paper wedescribe the methodology behind Robbie, and demonstrate its use on real and simulated data. Robbie is available on GitHub.

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.

Powerful radio-galaxies feature heavily in our understanding of galaxy evolution. However, when it comes to studying their properties as a function of redshift and/or environment, the most-detailed studies tend to be limited by small-number statistics. During Focus Meeting 3, on "Radio Galaxies: Resolving the AGN phenomenon", SVW presented a new sample of nearly 2,000 of the brightest radio-sources in the southern sky (Dec. $<$ 30 deg). These were observed at low radio-frequencies as part of the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, which is a continuum survey conducted using the Murchison Widefield Array (MWA). This instrument is the precursor telescope for the low-frequency component of the Square Kilometre Array, and allows us to select radio galaxies in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). Being brighter than 4 Jy at 151 MHz, we refer to these objects as the GLEAM 4-Jy (G4Jy) Sample. The G4Jy catalogue is close to being finalised, with SVW describing how multi-wavelength data have been used to determine the morphology of the radio emission, and identify the host galaxy. In addition, the MWA's excellent spectral-coverage and sensitivity to extended/diffuse emission were highlighted. Both of these aspects are important for understanding the physical mechanisms that take place within active galaxies, and how they interact with their environment.

The low-frequency polarisation properties of radio sources are poorly studied, particularly in statistical samples. However, the new generation of low-frequency telescopes, such as the Murchison Widefield Array (MWA; the precursor for the low-frequency component of the Square Kilometre Array) offers an opportunity to probe the physics of radio sources at very low radio frequencies. In this paper, we present a catalogue of linearly-polarised sources detected at 216 MHz, using data from the Galactic and Extragalactic All-sky MWA (GLEAM) survey. Our catalogue covers the Declination range $-17^{\circ}$ to $-37^{\circ}$ and 24 hours in Right Ascension, at a resolution of around 3 arcminutes. We detect 81 sources (including both a known pulsar and new pulsar candidate) with linearly-polarised flux densities in excess of 18 mJy across a survey area of approximately 6400 square degrees, corresponding to a surface density of 1 source per 79 square degrees. The level of Faraday rotation measured for our sources is broadly consistent with those recovered at higher frequencies, with typically more than an order of magnitude improvement in the uncertainty compared to higher-frequency measurements. However, our catalogue is likely incomplete at low Faraday rotation measures, due to our practice of excluding sources in the region where instrumental leakage appears. The majority of sources exhibit significant depolarisation compared to higher frequencies; however, a small sub-sample repolarise at 216 MHz. We also discuss the polarisation properties of four nearby, large-angular-scale radio galaxies, with a particular focus on the giant radio galaxy ESO 422$-$G028, in order to explain the striking differences in polarised morphology between 216 MHz and 1.4 GHz.

The Earth's ionosphere refracts radio waves incident on an interferometer, resulting in shifts to the measured positions of radio sources. We present a method to smoothly remove these shifts and restore sources to their reference positions, in both the catalogue and image domains. The method is applicable to instruments and ionospheric weather such that all antennas see the same ionosphere. The method is generalisable to repairing any sparsely-sampled vector field distortion to some input data. The code is available under the Academic Free License (https://opensource.org/licenses/AFL-3.0) from https://github.com/nhurleywalker/fits_warp

We derive the synchrotron distribution in the Milky Way disk from HII region absorption observations over -40{\deg} < l < 40{\deg} at six frequencies of 76.2, 83.8, 91.5, 99.2, 106.9, and 114.6 MHz with the GaLactic and Extragalactic All-sky Murchison widefield array survey (GLEAM). We develop a new method of emissivity calculation by taking advantage of the Haslam et al., (1981) map and known spectral indices, which enable us to simultaneously derive the emissivity and the optical depth of HII regions at each frequency. We show our derived synchrotron emissivities based on 152 absorption features of HII regions using both the method previously adopted in the literature and our improved method. We derive the synchrotron emissivity from HII regions to the Galactic edge along the line of sight and, for the first time, derive the emissivity from HII regions to the Sun. These results provide direct information on the distribution of the Galactic magnetic field and cosmic-ray electrons for future modelling.

We present a pilot search for variable and transient sources at 1.4 GHz with the Australian Square Kilometre Array Pathfinder (ASKAP). The search was performed in a 30 deg$^{2}$ area centred on the NGC 7232 galaxy group over 8 epochs and observed with a near-daily cadence. The search yielded nine potential variable sources, rejecting the null hypothesis that the flux densities of these sources do not change with 99.9% confidence. These nine sources displayed flux density variations with modulation indices m $\geq 0.1$ above our flux density limit of 1.5 mJy. They are identified to be compact AGN/quasars or galaxies hosting an AGN, whose variability is consistent with refractive interstellar scintillation. We also detect a highly variable source with modulation index m $ > 0.5$ over a time interval of a decade between the Sydney University Molonglo Sky Survey (SUMSS) and our latest ASKAP observations. We find the source to be consistent with the properties of long-term variability of a quasar. No transients were detected on timescales of days and we place an upper limit $\rho < 0.01$ deg$^{2}$ with 95% confidence for non-detections on near-daily timescales. The future VAST-Wide survey with 36-ASKAP dishes will probe the transient phase space with a similar cadence to our pilot survey, but better sensitivity, and will detect and monitor rarer brighter events.

Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_\nu \propto \nu^\alpha$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteor emission and other transient events using the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down to a $5\sigma$ detection threshold of 3.5 Jy/beam, no transient candidates were identified as intrinsic emission from meteors. We derived an upper limit of -3.7 (95% confidence limit) on the spectral index in our frequency range. We also report detections of other transient events, like reflected FM broadcast signals from small satellites, conclusively demonstrating the ability of the MWA to detect and track space debris on scales as small as 0.1 m in low Earth orbits.

We report the discovery of an unusual stellar system, SMSS J130522.47-293113.0. The optical spectrum is dominated by a blue continuum together with emission lines of hydrogen, neutral and ionized helium, and the N III, C III blend at around 4640-4650 Angstrom. The emission line profiles vary in strength and position on timescales as short as 1 day, while optical photometry reveals fluctuations of as much as ~0.2 mag in g on timescales as short as 10-15 min. The system is a weak X-ray source (f_{0.3-10} = (1.2 +/- 0.1) x 10^{-13} ergs cm^{-2} s^{-1} in the 0.3-10 keV band) but is not detected at radio wavelengths (3-sigma upper limit of 50 microJy at 5.5 GHz). The most intriguing property of the system, however, is the existence of two "blobs", a few arcsec in size, that are symmetrically located 3.8 arcmin (2.2 pc for our preferred system distance of ~2 kpc) each side of the central object. The blobs are detected in optical and near-IR broadband images but do not show any excess emission in Halpha images. We discuss the interpretation of the system, suggesting that the central object is most likely a nova-like CV, and that the blobs are relics of a pc-scale accretion-powered collimated outflow.

We present Very Long Baseline Interferometry observations of a faint and low-luminosity ($L_{\rm 1.4 GHz} < 10^{27}~\mbox{W Hz}^{-1}$) Gigahertz-Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sample. We select eight sources from deep radio observations that have radio spectra characteristic of a GPS or CSS source and an angular size of $\theta \lesssim 2$ arcsec, and detect six of them with the Australian Long Baseline Array. We determine their linear sizes, and model their radio spectra using Synchrotron Self Absorption (SSA) and Free Free Absorption (FFA) models. We derive statistical model ages, based on a fitted scaling relation, and spectral ages, based on the radio spectrum, which are generally consistent with the hypothesis that GPS and CSS sources are young and evolving. We resolve the morphology of one CSS source with a radio luminosity of $10^{25}~\mbox{W Hz}^{-1}$, and find what appear to be two hotspots spanning 1.7 kpc. We find that our sources follow the turnover-linear size relation, and that both homogenous SSA and an inhomogeneous FFA model can account for the spectra with observable turnovers. All but one of the FFA models do not require a spectral break to account for the radio spectrum, while all but one of the alternative SSA and power law models do require a spectral break to account for the radio spectrum. We conclude that our low-luminosity sample is similar to brighter samples in terms of their spectral shape, turnover frequencies, linear sizes, and ages, but cannot test for a difference in morphology.

We present new criteria for selecting HII regions from the Infrared Astronomical Satellite (IRAS) Point Source catalogue (PSC), based on an HII region catalogue derived manually from the all-sky Wide-field Infrared Survey Explorer (WISE). The criteria are used to augment the number of HII region candidates in the Milky Way. The criteria are defined by the linear decision boundary of two samples: IRAS point sources associated with known HII regions, which serve as the HII region sample, and IRAS point sources at high Galactic latitudes, which serve as the non-HII region sample. A machine learning classifier, specifically a support vector machine (SVM), is used to determine the decision boundary. We investigate all combinations of four IRAS bands and suggest that the optimal criterion is log(F$_{\rm 60}$/F$_{\rm 12}$)$\ge$(-0.19$\times$log(F$_{\rm 100}$/F$_{\rm 25}$)+ 1.52), with detections at 60 and 100 micron. This selects 3041 HII region candidates from the IRAS PSC. We find that IRAS HII region candidates show evidence of evolution on the two-colour diagram. Merging the WISE HII catalogue with IRAS HII region candidates, we estimate a lower limit of approximately 10200 for the number of HII regions in the Milky Way.

In the era of the SKA precursors, telescopes are producing deeper, larger images of the sky on increasingly small time-scales. The greater size and volume of images place an increased demand on the software that we use to create catalogues, and so our source finding algorithms need to evolve accordingly. In this paper we discuss some of the logistical and technical challenges that result from the increased size and volume of images that are to be analysed, and demonstrate how the Aegean source finding package has evolved to address these challenges. In particular we address the issues of source finding on spatially correlated data, and on images in which the background, noise, and point spread function, vary across the sky. We also introduce the concept of forced or priorized fitting.

We describe an image-based method that uses two radio criteria, compactness and spectral index, to identify promising pulsar candidates among Fermi Large Area Telescope (LAT) unassociated sources. These criteria are applied to those radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150 MHz (TGSS ADR1) found within the error ellipses of unassociated sources from the 3FGL catalog and a preliminary source list based on 7 years of LAT data. After follow-up interferometric observations to identify extended or variable sources, a list of 16 compact, steep-spectrum candidates is generated. An ongoing search for pulsations in these candidates, in gamma rays and radio, has found six millisecond pulsars and one normal pulsar. A comparison of this method with existing selection criteria based on gamma-ray spectral and variability properties suggests that the pulsar discovery space using Fermi may be larger than previously thought. Radio imaging is a hitherto underutilized source selection method that can be used, as with other multi-wavelength techniques, in the search for Fermi pulsars.

The Murchison Widefield Array (MWA), located in Western Australia, is one of the low-frequency precursors of the international Square Kilometre Array (SKA) project. In addition to pursuing its own ambitious science program, it is also a testbed for wide range of future SKA activities ranging from hardware, software to data analysis. The key science programs for the MWA and SKA require very high dynamic ranges, which challenges calibration and imaging systems. Correct calibration of the instrument and accurate measurements of source flux densities and polarisations require precise characterisation of the telescope's primary beam. Recent results from the MWA GaLactic Extragalactic All-sky MWA (GLEAM) survey show that the previously implemented Average Embedded Element (AEE) model still leaves residual polarisations errors of up to 10-20 % in Stokes Q. We present a new simulation-based Full Embedded Element (FEE) model which is the most rigorous realisation yet of the MWA's primary beam model. It enables efficient calculation of the MWA beam response in arbitrary directions without necessity of spatial interpolation. In the new model, every dipole in the MWA tile (4 x 4 bow-tie dipoles) is simulated separately, taking into account all mutual coupling, ground screen and soil effects, and therefore accounts for the different properties of the individual dipoles within a tile. We have applied the FEE beam model to GLEAM observations at 200 - 231 MHz and used false Stokes parameter leakage as a metric to compare the models. We have determined that the FEE model reduced the magnitude and declination-dependent behaviour of false polarisation in Stokes Q and V while retaining low levels of false polarisation in Stokes U.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

We have acquired radio continuum data between 70\,MHz and 48\,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500\,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500\,MHz the radio-continuum at low frequency ($\nu < 200$\,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $\alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.

Synchrotron emission pervades the Galactic plane at low radio frequencies, originating from cosmic ray electrons interacting with the Galactic magnetic field. Using a low-frequency radio telescope, the Murchison Widefield Array (MWA), we measure the free-free absorption of this Galactic synchrotron emission by intervening HII regions along the line of sight. These absorption measurements allow us to calculate the Galactic cosmic-ray electron emissivity behind and in front of 47 detected HII regions in the region $250^\circ < l < 355^\circ$, $|b| < 2^\circ$. We find that all average emissivities between the HII regions and the Galactic edge along the line of sight ($\epsilon_b$) are in the range of 0.24$\,\,\sim\,\,$0.70$\,\,$K$\,\,$pc$^{-1}$ with a mean of 0.40$\,\,$K$\,\,$pc$^{-1}$ and a variance of 0.10$\,\,$K$\,\,$pc$^{-1}$ at 88$\,\,$MHz. Our best model, the Two-circle model, divides the Galactic disk into three regions using two circles centring on the Galactic centre. It shows a high emissivity region near the Galactic centre, a low emissivity region near the Galactic edge, and a medium emissivity region between these two regions, contrary to the trend found by previous studies.

We detail and present results from a pilot study to assess the feasibility of detecting molecular lines at low radio frequencies. We observed a 400 square degree region centred on the Galactic Centre with the Murchison Widefield Array (MWA) between 103 and 133\,MHz targeting 28 known molecular species that have significant transitions. The results of this survey yield tentative detections of nitric oxide (NO) and the mercapto radical (SH). Both of these molecules appear to be associated with evolved stars.

We present a study of the three-dimensional structure of the molecular clouds in the Galactic Centre (GC) using CO emission and OH absorption lines. Two CO isotopologue lines, $^{12}$CO ($J$=1$\rightarrow$0) and $^{13}$CO ($J$=1$\rightarrow$0), and four OH ground-state transitions, surveyed by the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH), contribute to this study. We develop a novel method to calculate the OH column density, excitation temperature, and optical depth precisely using all four OH lines, and we employ it to derive a three-dimensional model for the distribution of molecular clouds in the GC for six slices in Galactic latitude. The angular resolution of the data is 15.5 arcmin, which at the distance of the GC (8.34 kpc) is equivalent to 38 pc. We find that the total mass of OH in the GC is in the range 2400-5100 Solar mass . The face-on view at a Galactic latitude of b = 0{\deg} displays a bar-like structure with an inclination angle of 67.5 $\pm$ 2.1{\deg} with respect to the line of sight. No ring-like structure in the GC is evident in our data, likely due to the low spatial resolution of the CO and OH maps.

The current generation of experiments aiming to detect the neutral hydrogen signal from the Epoch of Reionisation (EoR) is likely to be limited by systematic effects associated with removing foreground sources from target fields. In this paper we develop a model for the compact foreground sources in one of the target fields of the MWA's EoR key science experiment: the `EoR1' field. The model is based on both the MWA's GLEAM survey and GMRT 150 MHz data from the TGSS survey, the latter providing higher angular resolution and better astrometric accuracy for compact sources than is available from the MWA alone. The model contains 5049 sources, some of which have complicated morphology in MWA data, Fornax A being the most complex. The higher resolution data show that 13% of sources that appear point-like to the MWA have complicated morphology such as double and quad structure, with a typical separation of 33~arcsec. We derive an analytic expression for the error introduced into the EoR two-dimensional power spectrum due to peeling close double sources as single point sources and show that for the measured source properties, the error in the power spectrum is confined to high $k_\bot$ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest ten mis-modelled sources in the field contribute 90% of the power bias in the data, suggesting that it is most critical to improve the models of the brightest sources. With this hybrid model we reprocess data from the EoR1 field and show a maximum of 8% improved calibration accuracy and a factor of two reduction in residual power in $k$-space from peeling these sources. Implications for future EoR experiments including the SKA are discussed in relation to the improvements obtained.

We present low-frequency spectral energy distributions of 60 known radio pulsars observed with the Murchison Widefield Array (MWA) telescope. We searched the GaLactic and Extragalactic All-sky MWA (GLEAM) survey images for 200-MHz continuum radio emission at the position of all pulsars in the ATNF pulsar catalogue. For the 60 confirmed detections we have measured flux densities in 20 x 8 MHz bands between 72 and 231 MHz. We compare our results to existing measurements and show that the MWA flux densities are in good agreement.

The S-band Polarisation All-Sky Survey (S-PASS) has observed the entire southern sky using the 64-metre Parkes radio telescope at 2.3GHz with an effective bandwidth of 184MHz. The surveyed sky area covers all declinations $\delta\leq 0^\circ$. To analyse compact sources the survey data have been re-processed to produce a set of 107 Stokes $I$ maps with 10.75arcmin resolution and the large scale emission contribution filtered out. In this paper we use these Stokes $I$ images to create a total intensity southern-sky extragalactic source catalogue at 2.3GHz. The source catalogue contains 23,389 sources and covers a sky area of 16,600deg$^2$, excluding the Galactic plane for latitudes $|b|<10^\circ$. Approximately 8% of catalogued sources are resolved. S-PASS source positions are typically accurate to within 35arcsec. At a flux density of 225mJy the S-PASS source catalogue is more than 95% complete, and $\sim$94% of S-PASS sources brighter than 500mJy beam$^{-1}$ have a counterpart at lower frequencies.

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).

We have studied radio haloes and relics in nine merging galaxy clusters using the Murchison Widefield Array (MWA). The images used for this study were obtained from the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey which was carried out at 5 frequencies, viz. 88, 118, 154, 188 and 215 MHz. We detect diffuse radio emission in 8 of these clusters. We have estimated the spectra of haloes and relics in these clusters over the frequency range 80-1400 MHz; the first such attempt to estimate their spectra at low frequencies. The spectra follow a power law with a mean value of $\alpha = -1.13\pm0.21$ for haloes and $\alpha = -1.2\pm0.19$ for relics where, $S \propto \nu^{\alpha}$. We reclassify two of the cluster sources as radio galaxies. The low frequency spectra are thus an independent means of confirming the nature of cluster sources. Five of the nine clusters host radio haloes. For the remaining four clusters, we place upper limits on the radio powers of possible haloes in them. These upper limits are a factor of 2-20 below those expected from the $L_{\rm X}-P_{\rm 1.4}$ relation. These limits are the lowest ever obtained and the implications of these limits to the hadronic model of halo emission are discussed.

Many astronomical sources produce transient phenomena at radio frequencies, but the transient sky at low frequencies (<300 MHz) remains relatively unexplored. Blind surveys with new widefield radio instruments are setting increasingly stringent limits on the transient surface density on various timescales. Although many of these instruments are limited by classical confusion noise from an ensemble of faint, unresolved sources, one can in principle detect transients below the classical confusion limit to the extent that the classical confusion noise is independent of time. We develop a technique for detecting radio transients that is based on temporal matched filters applied directly to time series of images rather than relying on source-finding algorithms applied to individual images. This technique has well-defined statistical properties and is applicable to variable and transient searches for both confusion-limited and non-confusion-limited instruments. Using the Murchison Widefield Array as an example, we demonstrate that the technique works well on real data despite the presence of classical confusion noise, sidelobe confusion noise, and other systematic errors. We searched for transients lasting between 2 minutes and 3 months. We found no transients and set improved upper limits on the transient surface density at 182 MHz for flux densities between ~20--200 mJy, providing the best limits to date for hour- and month-long transients.

We present a sample of 1,483 sources that display spectral peaks between 72 MHz and 1.4 GHz, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The GLEAM survey is the widest fractional bandwidth all-sky survey to date, ideal for identifying peaked-spectrum sources at low radio frequencies. Our peaked-spectrum sources are the low frequency analogues of gigahertz-peaked spectrum (GPS) and compact-steep spectrum (CSS) sources, which have been hypothesized to be the precursors to massive radio galaxies. Our sample more than doubles the number of known peaked-spectrum candidates, and 95% of our sample have a newly characterized spectral peak. We highlight that some GPS sources peaking above 5 GHz have had multiple epochs of nuclear activity, and demonstrate the possibility of identifying high redshift ($z > 2$) galaxies via steep optically thin spectral indices and low observed peak frequencies. The distribution of the optically thick spectral indices of our sample is consistent with past GPS/CSS samples but with a large dispersion, suggesting that the spectral peak is a product of an inhomogeneous environment that is individualistic. We find no dependence of observed peak frequency with redshift, consistent with the peaked-spectrum sample comprising both local CSS sources and high-redshift GPS sources. The 5 GHz luminosity distribution lacks the brightest GPS and CSS sources of previous samples, implying that a convolution of source evolution and redshift influences the type of peaked-spectrum sources identified below 1 GHz. Finally, we discuss sources with optically thick spectral indices that exceed the synchrotron self-absorption limit.

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN.

We present a search for transient and highly variable sources at low radio frequencies (150-200 MHz) that explores long timescales of 1-3 years. We conducted this search by comparing the TIFR GMRT Sky Survey Alternative Data Release 1 (TGSS ADR1) and the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey catalogues. To account for the different completeness thresholds in the individual surveys, we searched for compact GLEAM sources above a flux density limit of 100 mJy that were not present in the TGSS ADR1; and also for compact TGSS ADR1 sources above a flux density limit of 200 mJy that had no counterpart in GLEAM. From a total sample of 234 333 GLEAM sources and 275 612 TGSS ADR1 sources in the overlap region between the two surveys, there were 99658 GLEAM sources and 38 978 TGSS ADR sources that passed our flux density cutoff and compactness criteria. Analysis of these sources resulted in three candidate transient sources. Further analysis ruled out two candidates as imaging artefacts. We analyse the third candidate and show it is likely to be real, with a flux density of 182 +/- 26 mJy at 147.5 MHz. This gives a transient surface density of rho = (6.2 +/- 6) x 10-5 deg-2 . We present initial follow-up observations and discuss possible causes for this candidate. The small number of spurious sources from this search demonstrates the high reliability of these two new low-frequency radio catalogues.

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.

We present the extensive follow-up campaign on the afterglow of GRB 110715A at 17 different wavelengths, from X-ray to radio bands, starting 81 seconds after the burst and extending up to 74 days later. We performed for the first time a GRB afterglow observation with the ALMA observatory. We find that the afterglow of GRB 110715A is very bright at optical and radio wavelengths. We use optical and near infrared spectroscopy to provide further information about the progenitor's environment and its host galaxy. The spectrum shows weak absorption features at a redshift z = 0.8225, which reveal a host galaxy environment with low ionization, column density and dynamical activity. Late deep imaging shows a very faint galaxy, consistent with the spectroscopic results. The broadband afterglow emission is modelled with synchrotron radiation using a numerical algorithm and we determine the best fit parameters using Bayesian inference in order to constrain the physical parameters of the jet and the medium in which the relativistic shock propagates. We fitted our data with a variety of models, including different density profiles and energy injections. Although the general behaviour can be roughly described by these models, none of them are able to fully explain all data points simultaneously. GRB 110715A shows the complexity of reproducing extensive multi-wavelength broadband afterglow observations, and the need of good sampling in wavelength and time and more complex models to accurately constrain the physics of GRB afterglows.

We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.

Using the narrowband all-sky imager mode of the LWA1 we have now detected 30 transients at 25.6 MHz, 1 at 34 MHz, and 93 at 38.0 MHz. While we have only optically confirmed that 37 of these events are radio afterglows from meteors, evidence suggests that most, if not all, are. Using the beam-forming mode of the LWA1 we have also captured the broadband spectra between 22.0 and 55.0 MHz of four events. We compare the smooth, spectral components of these four events and fit the frequency dependent flux density to a power law, and find that the spectral index is time variable, with the spectrum steepening over time for each meteor afterglow. Using these spectral indices along with the narrow band flux density measurements of the 123 events at 25.6 and 38 MHz, we predict the expected flux densities and rates for meteor afterglows potentially observable by other low frequency radio telescopes.

Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.

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.

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.

We present Murchison Widefield Array observations of the supernova remnant (SNR) 1987A between 72 and 230 MHz, representing the lowest frequency observations of the source to date. This large lever arm in frequency space constrains the properties of the circumstellar medium created by the progenitor of SNR 1987A when it was in its red supergiant phase. As of late-2013, the radio spectrum of SNR 1987A between 72 MHz and 8.64 GHz does not show any deviation from a non-thermal power-law with a spectral index of $-0.74 \pm 0.02$. This spectral index is consistent with that derived at higher frequencies, beneath 100 GHz, and with a shock in its adiabatic phase. A spectral turnover due to free-free absorption by the circumstellar medium has to occur below 72 MHz, which places upper limits on the optical depth of $\leq$ 0.1 at a reference frequency of 72 MHz, emission measure of $\lesssim$ 13,000 cm$^{-6}$ pc, and an electron density of $\lesssim$ 110 cm$^{-3}$. This upper limit on the electron density is consistent with the detection of prompt radio emission and models of the X-ray emission from the supernova. The electron density upper limit implies that some hydrodynamic simulations derived a red supergiant mass loss rate that is too high, or a wind velocity that is too low. The mass loss rate of $\sim 5 \times 10^{-6}$ $M_{\odot}$ yr$^{-1}$ and wind velocity of 10 km s$^{-1}$ obtained from optical observations are consistent with our upper limits, predicting a current turnover frequency due to free-free absorption between 5 and 60 MHz.

We use archival data from the Phoenix Deep Survey to investigate the variable radio source population above 1mJy/beam at 1.4GHz. Given the similarity of this survey to other such surveys we take the opportunity to investigate the conflicting results which have appeared in the literature. Two previous surveys for variability conducted with the Very Large Array (VLA) achieved a sensitivity of 1mJy/beam. However, one survey found an areal density of radio variables on timescales of decades that is a factor of ~4 times greater than a second survey which was conducted on timescales of less than a few years. In the Phoenix deep field we measure the density of variable radio sources to be $\rho =0.98\mathrm{deg}^{-2}$ on timescales of 6 months to 8 years. We make use of WISE infrared cross-ids, and identify all variable sources as an AGN of some description. We suggest that the discrepancy between previous VLA results is due to the different time scales probed by each of the surveys, and that radio variability at 1.4 GHz is greatest on timescales of 2 - 5 years.

We present the first dedicated radio continuum survey of a Kepler K2 mission field, Field 1 covering the North Galactic Cap. The survey is wide field, contemporaneous, multi-epoch, and multi-resolution in nature and was conducted at low radio frequencies between 140 and 200 MHz. The multi-epoch and ultra wide field (but relatively low resolution) part of the survey was provided by 15 nights of observation with the Murchison Widefield Array (MWA) over a period of approximately a month, contemporaneous with K2 observations of the field. The multi-resolution aspect of the survey was provided by the low resolution (4') MWA imaging, complemented by non-contemporaneous but much higher resolution (20") observations using the Giant Metrewave Radio Telescope (GMRT). The survey is therefore sensitive to the details of radio structures across a wide range of angular scales. Consistent with other recent low radio frequency surveys, no significant radio transients or variables were detected in the survey. The resulting source catalogs consist of 1,085 and 1,468 detections in the two MWA observation bands (centered at 154 and 185 MHz, respectively) and 7,445 detections in the GMRT observation band (centered at 148 MHz), over 314 square degrees. The survey is presented as a significant resource for multi-wavelength investigations of the more than 21,000 target objects in the K2 field. We briefly examine our survey data against K2 target lists for dwarf star types (stellar types M and L) that have been known to produce radio flares.

We present 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 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.

Electron density irregularities in the ionosphere are known to be magnetically anisotropic, preferentially elongated along the lines of force. While many studies of their morphology have been undertaken by topside sounding and whistler measurements, it is only recently that detailed regional-scale reconstructions have become possible, enabled by the advent of widefield radio telescopes. Here we present a new approach for visualising and studying field-aligned irregularities (FAIs), which involves transforming interferometric measurements of TEC gradients onto a magnetic shell tangent plane. This removes the perspective distortion associated with the oblique viewing angle of the irregularities from the ground, facilitating the decomposition of dynamics along and across magnetic field lines. We apply this transformation to the dataset of Loi et al. [2015a], obtained on 15 October 2013 by the Murchison Widefield Array (MWA) radio telescope and displaying prominent FAIs. We study these FAIs in the new reference frame, quantifying field-aligned and field-transverse behaviour, examining time and altitude dependencies, and extending the analysis to FAIs on sub-array scales. We show that the inclination of the plane can be derived solely from the data, and verify that the best-fit value is consistent with the known magnetic inclination. The ability of the model to concentrate the fluctuations along a single spatial direction may find practical application to future calibration strategies for widefield interferometry, by providing a compact representation of FAI-induced distortions.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a new relatively low resolution, contiguous 72-231 MHz survey of the entire sky south of declination +25 deg. In this paper, we outline one approach to determine the relative contribution of system noise, classical confusion and sidelobe confusion in GLEAM images. An understanding of the noise and confusion properties of GLEAM is essential if we are to fully exploit GLEAM data and improve the design of future low-frequency surveys. Our early results indicate that sidelobe confusion dominates over the entire frequency range, implying that enhancements in data processing have the potential to further reduce the noise.

This paper outlines how the new GaLactic and Extragalactic All-sky MWA Survey (GLEAM, Wayth et al. 2015), observed by the Murchison Widefield Array covering the frequency range 72 - 231 MHz, allows identification of a new large, complete, sample of more than 2000 bright extragalactic radio sources selected at 151 MHz. With a flux density limit of 4 Jy this sample is significantly larger than the canonical fully-complete sample, 3CRR (Laing, Riley & Longair 1983). In analysing this small bright subset of the GLEAM survey we are also providing a first user check of the GLEAM catalogue ahead of its public release (Hurley-Walker et al. in prep). Whilst significant work remains to fully characterise our new bright source sample, in time it will provide important constraints to evolutionary behaviour, across a wide redshift and intrinsic radio power range, as well as being highly complementary to results from targeted, small area surveys.

Current models suggest gamma-ray bursts could be used as a way of probing Population III stars - the first stars in the early Universe. In this paper we use numerical simulations to demonstrate that late time radio observations of gamma-ray burst afterglows could provide a means of identifying bursts that originate from Population III stars, if these were highly massive, independently from their redshift. We then present the results from a pilot study using the Australia Telescope Compact Array at 17 GHz, designed to test the hypothesis that there may be Population III gamma-ray bursts amongst the current sample of known events. We observed three candidates plus a control gamma-ray burst, and make no detections with upper limits of 20-40 uJy at 500-1300 days post explosion.

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 $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field.

Compact binary mergers may have already been observed as they are the leading model for short gamma-ray bursts (sGRBs). Radioactive decay within the ejecta from these mergers is expected to produce an infra-red flare, dubbed macronova (or kilonova), on a time scale of a week. Recently two such macronova candidates were identified in followup observations of sGRBs, strengthening the possibility that those indeed arise from mergers. The same ejecta will also produce a long term (months to years) radio emission due to its interaction with the surrounding ISM. In search for this emission, we observed the two macronova candidates, GRB 130603B and GRB 060614 with the Jansky very large array (VLA) and the Australia Telescope Compact Array (ATCA). Our observations resulted in null-detections, putting strong upper limits on the kinetic energy and mass of the ejecta. A possible outcome of a merger is a highly magnetized neutron star (a magnetar), which has been suggested as the central engine for GRBs. Such a magnetar will deposit a significant fraction of its energy into the ejecta leading to a brighter radio flare. Our results, therefore, rule out magnetars in these two events.

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

We present the results of a pilot study search for Fast Radio Bursts (FRBs) using the Murchison Widefield Array (MWA) at low frequencies (139 - 170 MHz). We utilised MWA data obtained in a routine imaging mode from observations where the primary target was a field being studied for Epoch of Reionisation detection. We formed images with 2 second time resolution and 1.28~MHz frequency resolution for 10.5 hours of observations, over 400 square degrees of the sky. We de-dispersed the dynamic spectrum in each of 372,100 resolution elements of 2$\times$2 arcmin$^{2}$, between dispersion measures of 170 and 675~pc~cm$^{-3}$. Based on the event rate calculations in Trott, Tingay & Wayth (2013), which assumes a standard candle luminosity of $8\times10^{37}$ Js$^{-1}$, we predict that with this choice of observational parameters, the MWA should detect ($\sim10$,$\sim2$,$\sim0$) FRBs with spectral indices corresponding to ($-$2, $-$1, 0), based on a 7$\sigma$ detection threshold. We find no FRB candidates above this threshold from our search, placing an event rate limit of $<700$ above 700 Jy.ms per day per sky and providing evidence against spectral indices $\alpha<-1.2$ ($S\propto\nu^{\alpha}$). We compare our event rate and spectral index limits with others from the literature. We briefly discuss these limits in light of recent suggestions that supergiant pulses from young neutron stars could explain FRBs. We find that such supergiant pulses would have to have much flatter spectra between 150 and 1400 MHz than have been observed from Crab giant pulses to be consistent with the FRB spectral index limit we derive.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

We present observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array.

We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with the ionosphere as inferred from Global Positioning System (GPS) data. The first order ionosphere manifests itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be computed from global ionosphere maps provided by GPS analysis centres, namely the Center for Orbit Determination in Europe (CODE), using data from globally distributed GPS receivers. However, for the more accurate local ionosphere estimates required for precision radio astronomy applications, data from local GPS networks needs to be incorporated into ionospheric modelling. For GPS observations, the ionospheric parameters are biased by GPS receiver instrument delays, among other effects, also known as receiver Differential Code Biases (DCBs). The receiver DCBs need to be estimated for any non-CODE GPS station used for ionosphere modelling, a requirement for establishing dense GPS networks in arbitrary locations in the vicinity of the MWA. In this work, single GPS station-based ionospheric modelling is performed at a time resolution of 10 minutes. Also the receiver DCBs are estimated for selected Geoscience Australia (GA) GPS receivers, located at Murchison Radio Observatory (MRO1), Yarragadee (YAR3), Mount Magnet (MTMA) and Wiluna (WILU). The ionospheric gradients estimated from GPS are compared with the ionospheric gradients inferred from radio source position shifts observed with the MWA. The ionospheric gradients at all the GPS stations show a correlation with the gradients observed with the MWA. The ionosphere estimates obtained using GPS measurements show promise in terms of providing calibration information for the MWA.

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.

We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the frequencies. The relics, with a linear extent of $\sim$ 1 Mpc each, are separated by $\sim$ 2 Mpc. Combining the current observations with those in the literature, we have obtained the spectra of these relics over the frequency range 80 -- 1400 MHz. The spectra follow power laws, with $\alpha$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto \nu^{\alpha}$). Assuming the break frequency to be near the lower end of the spectrum we estimate the age of the relics to be $\sim$ 0.4 Gyr. No diffuse radio emission from the central regions of the cluster (halo) was detected. The upper limit on the radio power of any possible halo that might be present in the cluster is a factor of 35 lower than that expected from the radio power and X-ray luminosity correlation for cluster halos. From this we conclude that the cluster halo is very extended ($>$ 500 kpc) and/or most of the radio emission from the halo has decayed. The current limit on the halo radio power is a factor of ten lower than the existing upper limits with possible implications for models of halo formation.

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.

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.

We present a three epoch survey for transient and variables in the extended Chandra Deep Field South at 5.5 GHz with the Australia Telescope Compact Array. A region covering $\sim$0.3 deg$^{2}$ was observed on timescales of 2.5 months and 2.5 years and typical sensitivities 12.1 $-$ 17.1 $\mu$Jy beam$^{-1}$ (1$\sigma$) were achieved. This survey represents the deepest search for transient and variable radio sources at 5.5 GHz. In total 124 sources were detected above the 5.5$\sigma$ level. One highly variable radio source was found with $\Delta S > 50%$ implying a surface density of $\sim$3 deg$^{-2}$. A further three radio sources were found with lower levels of variability equating to a surface density of $\sim$13 deg$^{-2}$ above a detection threshold of 82.3 $\mu$Jy. All of the variable sources have inverted radio spectra (between 1.4 and 5.5 GHz) and are associated with active galactic nuclei. We conclude that these variables are young gigahertz peaked-spectrum sources with active and self-absorbed radio jets. We explore the variability completeness of this sample and conclude that the fairly low levels of variability would only be detectable in 3$-$25% of all sources within the field. No radio transients were detected in this survey and we place an upper limit on the surface density of transient events $ < 7.5$ deg$^{-2}$ above a detection threshold of 68.8~$\mu$Jy.

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).

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 $<0.1$\%. We estimate the space density of this low brightness (dying) phase of radio galaxy evolution as $7\times10^{-7}$\,Mpc$^{-3}$ and argue that normal AGN cannot spend more than 6\% of their lifetime in this phase if they all go through the same cycle.

We present the results of an approximately 6,100 square degree 104--196MHz radio sky survey performed with the Murchison Widefield Array during instrument commissioning between 2012 September and 2012 December: the Murchison Widefield Array Commissioning Survey (MWACS). The data were taken as meridian drift scans with two different 32-antenna sub-arrays that were available during the commissioning period. The survey covers approximately 20.5 h < Right Ascension (RA) < 8.5 h, -58 deg < Declination (Dec) < -14 deg over three frequency bands centred on 119, 150 and 180 MHz, with image resolutions of 6--3 arcmin. The catalogue has 3-arcmin angular resolution and a typical noise level of 40 mJy/beam, with reduced sensitivity near the field boundaries and bright sources. We describe the data reduction strategy, based upon mosaiced snapshots, flux density calibration and source-finding method. We present a catalogue of flux density and spectral index measurements for 14,110 sources, extracted from the mosaic, 1,247 of which are sub-components of complexes of sources.

We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. Fifty-one images were obtained that cover a field of view of 1430 deg^2 centred on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 minutes and there was additional temporal information on day and year timescales. We explore the variability of a sample of 105 low frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short timescale variability (26 minutes). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the timescale of one year we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources rho < 7.5 x 10^-5 deg^-2 with flux densities > 5.5 Jy, and characteristic timescales of both 26 minutes and one year.

The intermediate Palomar Transient Factory reports our discovery of a young supernova, iPTF13bvn, in the nearby galaxy, NGC5806 (22.5Mpc). Our spectral sequence in the optical and infrared suggests a likely Type Ib classification. We identify a single, blue progenitor candidate in deep pre-explosion imaging within a 2{\sigma} error circle of 80 mas (8.7 pc). The candidate has a MB luminosity of -5.2 +/- 0.4 mag and a B-I color of 0.1+/-0.3 mag. If confirmed by future observations, this would be the first direct detection for a progenitor of a Type Ib. Fitting a power law to the early light curve, we find an extrapolated explosion date around 1.1 days before our first detection. We see no evidence of shock cooling. The pre-explosion detection limits constrain the radius of the progenitor to be smaller than a few solar radii. iPTF13bvn is also detected in cm and mm-wavelengths. Fitting a synchrotron self-absorption model to our radio data, we find a mass loading parameter of 1.3*10^12 g/cm. Assuming a wind velocity of 10^3km/s, we derive a progenitor mass loss rate of 3*10^-5Msun/yr. Our observations, taken as a whole, are consistent with a Wolf Rayet progenitor of the supernova iPTF13bvn.

The detection rate of gamma-ray burst (GRB) afterglows is only ~30% at radio wavelengths, much lower than in the X-ray (~95%) or optical (~70%) bands. The cause of this low radio detection rate has previously been attributed to limited observing sensitivity. We use visibility stacking to test this idea, and conclude that the low detection rate is instead due to two intrinsically different populations of GRBs, radio bright and radio faint. We calculate that no more than 70% of GRB afterglows are truly radio bright, leaving a significant population of GRBs that lack a radio afterglow. These radio bright GRBs have higher gamma-ray fluence, isotropic energies, X-ray fluxes and optical fluxes than the radio faint GRBs, confirming the existence of two physically distinct populations. We suggest that the gamma-ray efficiency of the prompt emission is responsible for the difference between the two populations. We also discuss the implications for future radio and optical surveys.

The Australian Square Kilometre Array Pathfinder (ASKAP) will give us an unprecedented opportunity to investigate the transient sky at radio wavelengths. In this paper we present VAST, an ASKAP survey for Variables and Slow Transients. VAST will exploit the wide-field survey capabilities of ASKAP to enable the discovery and investigation of variable and transient phenomena from the local to the cosmological, including flare stars, intermittent pulsars, X-ray binaries, magnetars, extreme scattering events, interstellar scintillation, radio supernovae and orphan afterglows of gamma ray bursts. In addition, it will allow us to probe unexplored regions of parameter space where new classes of transient sources may be detected. In this paper we review the known radio transient and variable populations and the current results from blind radio surveys. We outline a comprehensive program based on a multi-tiered survey strategy to characterise the radio transient sky through detection and monitoring of transient and variable sources on the ASKAP imaging timescales of five seconds and greater. We also present an analysis of the expected source populations that we will be able to detect with VAST.

We present a detailed analysis of four of the most widely used radio source finding packages in radio astronomy, and a program being developed for the Australian Square Kilometer Array Pathfinder (ASKAP) telescope. The four packages; SExtractor, SFind, IMSAD and Selavy are shown to produce source catalogues with high completeness and reliability. In this paper we analyse the small fraction (~1%) of cases in which these packages do not perform well. This small fraction of sources will be of concern for the next generation of radio surveys which will produce many thousands of sources on a daily basis, in particular for blind radio transients surveys. From our analysis we identify the ways in which the underlying source finding algorithms fail. We demonstrate a new source finding algorithm Aegean, based on the application of a Laplacian kernel, which can avoid these problems and can produce complete and reliable source catalogues for the next generation of radio surveys.

The process of determining the number and characteristics of sources in astronomical images is so fundamental to a large range of astronomical problems that it is perhaps surprising that no standard procedure has ever been defined that has well understood properties with a high degree of statistical rigour on completeness and reliability. There are now a large number of commonly used software tools for accomplishing this task, typically with different tools being used for images acquired using different technologies. Despite this, there have been relatively few quantitative analyses of the robustness or reliability of individual tools, or the details of the techniques they implement. We have an opportunity to redress this omission in the context of surveys planned with the Australian Square Kilometre Array Pathfinder (ASKAP). The Evolutionary Map of the Universe (EMU) survey with ASKAP, a continuum survey of the Southern Hemisphere up to declination +30 deg, aims to utilise an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. A key stage in source extraction methods is the background estimation (background level and noise level) and the choice of a threshold high enough to reject false sources yet not so high that the catalogues are significantly incomplete. In this analysis we present results from testing such algorithms as implemented in the SExtractor, Selavy (Duchamp), and sfind tools on simulated data. In particular the effects of background estimation, threshold and false-discovery rate settings are explored. For parameters that give similar completeness, the false-discovery rate method employed by sfind results in a more reliable catalogue compared to the peak threshold methods of SExtractor and Selavy.

The Australia Telescope 20GHz (AT20G) survey is a large area (2{\pi} sr), sensitive (40mJy), high frequency (20GHz) survey of the southern sky. The survey was conducted in two parts: an initial fast scanning survey, and a series of more accurate follow-up observations. The follow-up survey catalog has been presented by Murphy et al. 2010. In this paper we discuss the hardware setup and scanning survey strategy as well as the production of the scanning survey catalog.

Our current understanding of radio-loud AGN comes predominantly from studies at frequencies of 5 GHz and below. With the recent completion of the Australia Telescope 20 GHz (AT20G) survey, we can now gain insight into the high-frequency radio properties of AGN. This paper presents supplementary information on the AT20G sources in the form of optical counterparts and redshifts. Optical counterparts were identified using the SuperCOSMOS database and redshifts were found from either the 6dF Galaxy survey or the literature. We also report 144 new redshifts. For AT20G sources outside the Galactic plane, 78.5% have optical identifications and 30.9% have redshift information. The optical identification rate also increases with increasing flux density. Targets which had optical spectra available were examined to obtain a spectral classification. There appear to be two distinct AT20G populations; the high luminosity quasars that are generally associated with point-source optical counterparts and exhibit strong emission lines in the optical spectrum, and the lower luminosity radio galaxies that are generally associated with passive galaxies in both the optical images and spectroscopic properties. It is suggested that these different populations can be associated with different accretion modes (cold-mode or hot-mode). We find that the cold-mode sources have a steeper spectral index and produce more luminous radio lobes, but generally reside in smaller host galaxies than their hot-mode counterparts. This can be attributed to the fact that they are accreting material more efficiently. Lastly, we compare the AT20G survey with the S-cubed semi-empirical (S3-SEX) models and conclude that the S3-SEX models need refining to correctly model the compact cores of AGN. The AT20G survey provides the ideal sample to do this.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

We describe the process of stacking radio interferometry visibilities to form a deep composite image and its application to the observation of transient phenomena. We apply "visibility stacking" to 46 archival Very Large Array observations of nearby type Ia supernovae (SNeIa). This new approach provides an upper limit on the SNIa ensemble peak radio luminosity of 1.2x10^{25}erg/s/Hz at 5GHz, which is 5-10 times lower than previously measured. This luminosity implies an upper limit on the average companion stellar wind mass loss rate of 1.3x10^{-7}M_o/yr. This mass loss rate is consistent with the double degenerate scenario for SNeIa and rules out intermediate and high mass companions in the single degenerate scenario. In the era of time domain astronomy, techniques such as visibility stacking will be important in extracting the maximum amount of information from observations of populations of short lived events.

The Australia Telescope 20 GHz (AT20G) survey is a blind survey of the whole Southern sky at 20 GHz with follow-up observations at 4.8, 8.6, and 20 GHz carried out with the Australia Telescope Compact Array (ATCA). In this paper we present an analysis of radio spectral properties in total intensity and polarisation, sizes, optical identifications, and redshifts of the sample of the 5808 extragalactic sources in the survey catalogue of confirmed sources over the whole Southern sky excluding the strip at Galactic latitude |b|<1.5deg. The sample has a flux density limit of 40 mJy. Completeness has been measured as a function of scan region and flux density. Averaging over the whole survey area the follow-up survey is 78% complete above 50mJy and 93% complete above 100mJy. 3332 sources with declination <-15deg have good quality almost simultaneous observations at 4.8, 8.6, and 20GHz. The spectral analysis shows that the sample is dominated by flat-spectrum sources. The fraction of flat-spectrum sources decreases from 81% for 20GHz flux densities S>500mJy, to 60% for S<100mJy. There is also a clear spectral steepening at higher frequencies with the median spectral index decreasing from -0.16 between 4.8 and 8.6GHz to -0.28 between 8.6 and 20GHz. Simultaneous observations in polarisation are available for all the sources at all the frequencies. 768 sources have a good quality detection of polarised flux density at 20GHz; 467 of them were also detected in polarisation at 4.8 and/or at 8.6GHz so that it has been possible to compare the spectral behaviour in total intensity and polarisation. We have found that the polarised fraction increases slightly with frequency and decreases with flux density. Cross matches and comparisons have been made with other catalogues at lower radio frequencies, and in the optical, X-ray and gamma-ray bands. Redshift estimates are available for 825 sources.

We used the Australia Telescope Compact Array (ATCA) to obtain 40 GHz and 95 GHz observations of a number of sources that were selected from the Australia Telescope Compact Array 20 GHz (AT20G) survey . The aim of the observations was to improve the spectral coverage for sources with spectral peaks near 20 GHz or inverted (rising) radio spectra between 8.6 GHz and 20 GHz. We present the radio observations of a sample of 21 such sources along with optical spectra taken from the ANU Siding Spring Observatory 2.3m telescope and the ESO-New Technology Telescope (NTT). We find that as a group the sources show the same level of variability as typical GPS sources, and that of the 21 candidate GPS sources roughly 60% appear to be genuinely young radio galaxies. Three of the 21 sources studied show evidence of being restarted radio galaxies. If these numbers are indicative of the larger population of AT20G radio sources then as many as 400 genuine GPS sources could be contained within the AT20G with up to 25% of them being restarted radio galaxies.

We present the full source catalogue from the Australia Telescope 20 GHz (AT20G) Survey. The AT20G is a blind radio survey carried out at 20 GHz with the Australia Telescope Compact Array (ATCA) from 2004 to 2008, and covers the whole sky south of declination 0 deg. The AT20G source catalogue presented here is an order of magnitude larger than any previous catalogue of high-frequency radio sources, and includes 5890 sources above a 20 GHz flux-density limit of 40 mJy. All AT20G sources have total intensity and polarisation measured at 20 GHz, and most sources south of declination -15 deg also have near-simultaneous flux-density measurements at 5 and 8 GHz. A total of 1559 sources were detected in polarised total intensity at one or more of the three frequencies. We detect a small but significant population of non-thermal sources that are either undetected or have only weak detections in low-frequency catalogues. We introduce the term Ultra-Inverted Spectrum (UIS) to describe these radio sources, which have a spectral index alpha(5, 20) > +0.7 and which constitute roughly 1.2 per cent of the AT20G sample. The 20 GHz flux densities measured for the strongest AT20G sources are in excellent agreement with the WMAP 5-year source catalogue of Wright et al. (2009), and we find that the WMAP source catalogue is close to complete for sources stronger than 1.5 Jy at 23 GHz.

GHz-peaked spectrum (GPS) radio sources are thought to be young objects which later evolve into FR-I and FR-II radio galaxies. We have used the Australia Telescope 20GHz (AT20G) survey catalogue to select a uniform sample of GPS sources with spectral peaks above 5GHz, which should represent the youngest members of this class. In this paper, we present e-VLBI observations of ten such objects which are associated with nearby (z<0.15) galaxies and so represent a new population of local, low--power GPS sources. Our e-VLBI observations were carried out at 4.8GHz with the Australia Telescope Long Baseline Array (LBA) using a real--time software correlator. All ten sources were detected, and were unresolved on scales of ~100mas, implying that they are typically less than 100pc in linear size.

The Australia Telescope 20GHz (AT20G) survey was used to select a complete sample of 656 Gigahertz Peaked Spectrum (GPS) sources with spectral turnovers above 5GHz. The AT20G has near simultaneous observations at 4.8, 8.6 and 20GHz, which makes it possible to exclude flat spectrum variability as a cause of a source's peaked spectrum. Optical identification of the sample results in 361 QSOs and 104 galaxies and 191 blank fields. Redshifts are known for 104 of the GPS sources. The GPS sources from the AT20G are discussed and compared to previously known samples. The new sample of high frequency peaking GPS sources is found at a lower redshift than previous samples and to also have a lower 5GHz radio power. Evidence is found to support the idea that the origin of the GPS spectral shape are intrinsically different for galaxies and QSOs. This paper is an elaboration and extension of the talk given at the $4^{th}$ CSS/GPS conference in Riccione in May 2008.

The Australia Telescope 20 GHz (AT20G) Survey is a blind survey of the whole Southern sky at 20 GHz (with follow-up observations at 4.8 and 8.6 GHz) carried out with the Australia Telescope Compact Array (ATCA) from 2004 to 2007. The Bright Source Sample (BSS) is a complete flux-limited subsample of the AT20G Survey catalogue comprising 320 extragalactic (|b|>1.5 deg) radio sources south of dec = -15 deg with S(20 GHz) > 0.50 Jy. Of these, 218 have near simultaneous observations at 8 and 5 GHz. In this paper we present an analysis of radio spectral properties in total intensity and polarisation, size, optical identifications and redshift distribution of the BSS sources. The analysis of the spectral behaviour shows spectral curvature in most sources with spectral steepening that increases at higher frequencies (the median spectral index \alpha, assuming S\propto \nu^\alpha, decreases from \alpha_{4.8}^{8.6}=0.11 between 4.8 and 8.6 GHz to \alpha_{8.6}^{20}=-0.16 between 8.6 and 20 GHz), even if the sample is dominated by flat spectra sources (85 per cent of the sample has \alpha_{8.6}^{20}>-0.5). The almost simultaneous spectra in total intensity and polarisation allowed us a comparison of the polarised and total intensity spectra: polarised fraction slightly increases with frequency, but the shapes of the spectra have little correlation. Optical identifications provided an estimation of redshift for 186 sources with a median value of 1.20 and 0.13 respectively for QSO and galaxies.

We have combined optical data from the 2dF-SDSS Luminous Red Galaxy and QSO (2SLAQ) redshift survey with radio measurements from the 1.4 GHz VLA FIRST and NVSS surveys to identify a volume-limited sample of 391 radio galaxies at redshift 0.4<z<0.7. By determining an accurate radio luminosity function for early-type galaxies in this redshift range, we can investigate the cosmic evolution of the radio-galaxy population over a wide range in radio luminosity. The low-power radio galaxies in our LRG sample (those with 1.4 GHz radio luminosities in the range 10^{24} to 10^{25} W/Hz, corresponding to FR I radio galaxies in the local universe) undergo significant cosmic evolution over the redshift range 0<z<0.7, consistent with pure luminosity evolution of the form (1+z)^k where k=2.0+/-0.3. Our results appear to rule out (at the 6-7 sigma level) models in which low-power radio galaxies undergo no cosmic evolution. The most powerful radio galaxies in our sample (with radio luminosities above 10^{26} W/Hz) may undergo more rapid evolution over the same redshift range. The evolution seen in the low-power radio-galaxy population implies that the total energy input into massive early-type galaxies from AGN heating increases with redshift, and was roughly 50% higher at z~0.55 (the median redshift of the 2SLAQ LRG sample) than in the local universe.

We present some first results on the variability, polarization and general properties of radio sources selected in a blind survey at 20 GHz, the highest frequency at which a sensitive radio survey has been carried out over a large area of sky. Sources with flux densities above 100 mJy in the AT20G Pilot Survey at declination -60 to -70 were observed at up to three epochs during 2002-4, including near-simultaneous measurements at 5, 8 and 18 GHz in 2003. Of the 173 sources detected, 65% are candidate QSOs, BL Lac objects or blazars, 20% galaxies and 15% faint (b > 22 mag) optical objects or blank fields. On a 1-2 year timescale, the general level of variability at 20 GHz appears to be low. For the 108 sources with good-quality measurements in both 2003 and 2004, the median variability index at 20 GHz was 6.9% and only five sources varied by more than 30% in flux density. Most sources in our sample show low levels of linear polarization (typically 1-5%), with a median fractional polarization of 2.3% at 20 GHz. There is a trend for fainter sources to show higher fractional polarization. At least 40% of sources selected at 20GHz have strong spectral curvature over the frequency range 1-20 GHz. We use a radio `two-colour diagram' to characterize the radio spectra of our sample, and confirm that the radio-source population at 20 GHz (which is also the foreground point-source population for CMB anisotropy experiments like WMAP and Planck) cannot be reliably predicted by extrapolating the results of surveys at lower frequencies. As a result, direct selection at 20 GHz appears to be a more efficient way of identifying 90 GHz phase calibrators for ALMA than the currently-proposed technique of extrapolation from all-sky surveys at 1-5 GHz.