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.

Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds.

We present early results from a project to measure the sky-averaged (global), redshifted $21\,$cm signal from the Epoch of Reionisation (EoR), using the Murchison Widefield Array (MWA) telescope. Because interferometers are not sensitive to a spatially-invariant global average, they cannot be used to detect this signal using standard techniques. However, lunar occultation of the radio sky imprints a spatial structure on the global signal, allowing us to measure the average brightness temperature of the patch of sky immediately surrounding the Moon. In this paper we present one night of Moon observations with the MWA between 72 - 230 MHz and verify our techniques to extract the background sky temperature from measurements of the Moon's flux density. We improve upon previous work using the lunar occultation technique by using a more sophisticated model for reflected `earthshine' and by employing image differencing to remove imaging artefacts. We leave the Moon's (constant) radio brightness temperature as a free parameter in our fit to the data and as a result, measure $T_{\rm{moon}} = 180 \pm 12 $ K and a Galactic synchrotron spectral index of $-2.64\pm0.14$, at the position of the Moon. Finally, we evaluate the prospects of the lunar occultation technique for a global EoR detection and map out a way forward for future work with the MWA.

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV-coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL, and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results. (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD, and demonstrate these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EoR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes.

We use observations made with the Giant Metrewave Radio Telescope (GMRT) to probe the neutral hydrogen (HI) gas content of field galaxies in the VIMOS VLT Deep Survey (VVDS) 14h field at $z \approx 0.32$. Because the HI emission from individual galaxies is too faint to detect at this redshift, we use an HI spectral stacking technique using the known optical positions and redshifts of the 165 galaxies in our sample to co-add their HI spectra and thus obtain the average HI mass of the galaxies. Stacked HI measurements of 165 galaxies show that 95 per cent of the neutral gas is found in blue, star-forming galaxies. Among these galaxies, those having lower stellar mass are more gas-rich than more massive ones. We apply a volume correction to our HI measurement to evaluate the HI gas density at $z \approx 0.32$ as $\Omega_{HI}=(0.50\pm0.18) \times 10^{-3}$ in units of the cosmic critical density. This value is in good agreement with previous results at z < 0.4, suggesting no evolution in the neutral hydrogen gas density over the last $\sim 4$ Gyr. However the $z \approx 0.32$ gas density is lower than that at $z \sim 5$ by at least a factor of two.

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

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

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.

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.

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

The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $\Delta^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications.

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.

Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.

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

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

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window.

We present the results of HI spectral stacking analysis of Giant Metrewave Radio Telescope (GMRT) observations targeting the COSMOS field. The GMRT data cube contains 474 field galaxies with redshifts known from the zCOSMOS-bright 10k catalogue. Spectra for the galaxies are co-added and the stacked spectrum allows us to make a $\sim 3\sigma$ measurement of the average HI mass. Using this average HI mass along with the integral optical $B$-band luminosity of the galaxies and the luminosity density of the COSMOS field, a volume normalisation is applied to obtain the cosmic HI mass density ($\Omega_{\rm HI}$). We find a cosmic HI mass density of $\Omega_{\rm HI}$ = (0.42 $\pm$ 0.16) $\times$ 10$^{-3}$ at $z \sim 0.37$, which is the highest-redshift measurement of $\Omega_{\rm HI}$ ever made using HI spectral stacking. The value we obtained for $\Omega_{\rm HI}$ at $z \sim 0.37$ is consistent with that measured from large blind 21-cm surveys at $z = 0$ as well as measurements from other HI stacking experiments at lower redshifts. Our measurement in conjunction with earlier measurements indicates that there has been no significant evolution of HI gas abundance over the last 4 Gyr. A weighted mean of $\Omega_{\rm HI}$ from all 21-cm measurements at redshifts $z \lesssim 0.4$ gives $\Omega_{\rm HI}$ = (0.35 $\pm$ 0.01) $\times$ 10$^{-3}$. The $\Omega_{\rm HI}$ measured (from HI 21-cm emission measurements) at $z \lesssim 0.4$ is however approximately half that measured from Damped Lyman-$\alpha$ Absorption (DLA) systems at $z \gtrsim 2$. Deeper surveys with existing and upcoming instruments will be critical to understand the evolution of $\Omega_{\rm HI}$ in the redshift range intermediate between $z \sim 0.4$ and the range probed by DLA observations.

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

Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle \alpha_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population.

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.

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

Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by PSF sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source.

In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum.

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

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.

Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz.

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 report on observations of giant pulses from the Crab pulsar performed simultaneously with the Parkes radio telescope and the incoherent combination of the Murchison Widefield Array (MWA) antenna tiles. The observations were performed over a duration of approximately one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5$\sigma$ and 6.5$\sigma$ respectively. We detected 51$\%$ of the MWA giant pulses at the Parkes radio telescope, with spectral indices in the range of $-3.6>\alpha> -4.9$ ($S_{\rm \nu} \propto \nu^\alpha$). We present a Monte Carlo analysis supporting the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than $100\%$ correlation being due to the relative sensitivities of the two instruments and the width of the spectral index distribution. Our observations are consistent with the hypothesis that the spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 MHz.

We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known spectral width of any GPS source. Spectral coverage of the source spans from 0.118 to 22 GHz, which includes data from the Murchison Widefield Array and the wide bandpass receivers on the Australia Telescope Compact Array. We have implemented a Bayesian inference model fitting routine to fit the data with various absorption models. We find that without the inclusion of a high-frequency exponential break the absorption models can not accurately fit the data, with significant deviations above and below the peak in the radio spectrum. The addition of a high-frequency break provides acceptable spectral fits for the inhomogeneous free-free absorption and double-component synchrotron self-absorption models, with the inhomogeneous free-free absorption model statistically favored. The requirement of a high-frequency spectral break implies that the source has ceased injecting fresh particles. Additional support for the inhomogeneous free-free absorption model as being responsible for the turnover in the spectrum is given by the consistency between the physical parameters derived from the model fit and the implications of the exponential spectral break, such as the necessity of the source being surrounded by a dense ambient medium to maintain the peak frequency near the gigahertz region. The discovery of PKS B0008-421 suggests that the next generation of low radio frequency surveys could reveal a large population of GPS sources that have ceased activity, and that a portion of the ultra-steep spectrum source population could be composed of these GPS sources in a relic phase.

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.

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

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.

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

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.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

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

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.

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

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.

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

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

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

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

Using the new wideband capabilities of the Australia Telescope Compact Array (ATCA), we obtain spectra for PKS 1718-649, a well-known gigahertz-peaked spectrum radio source. The observations, between approximately 1 and 10 GHz over three epochs spanning approximately 21 months, reveal variability both above the spectral peak at ~3 GHz and below the peak. The combination of the low and high frequency variability cannot be easily explained using a single absorption mechanism, such as free-free absorption or synchrotron self-absorption. We find that the PKS 1718-649 spectrum and its variability are best explained by variations in the free-free optical depth on our line-of-sight to the radio source at low frequencies (below the spectral peak) and the adiabatic expansion of the radio source itself at high frequencies (above the spectral peak). The optical depth variations are found to be plausible when X-ray continuum absorption variability seen in samples of Active Galactic Nuclei is considered. We find that the cause of the peaked spectrum in PKS 1718-649 is most likely due to free-free absorption. In agreement with previous studies, we find that the spectrum at each epoch of observation is best fit by a free-free absorption model characterised by a power-law distribution of free-free absorbing clouds. This agreement is extended to frequencies below the 1 GHz lower limit of the ATCA by considering new observations with Parkes at 725 MHz and 199 MHz observations with the newly operational Murchison Widefield Array. These lower frequency observations argue against families of absorption models (both free-free and synchrotron self-absorption) that are based on simple homogenous structures.

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 new low-frequency observations of the nearby radio galaxy Fornax A at 154 MHz with the Murchison Widefield Array, microwave flux-density measurements obtained from WMAP and Planck data, and gamma-ray flux densities obtained from Fermi data. We also compile a comprehensive list of previously published images and flux-density measurements at radio, microwave and X-ray energies. A detailed analysis of the spectrum of Fornax A between 154 MHz and 1510 MHz reveals that both radio lobes have a similar spatially-averaged spectral index, and that there exists a steep-spectrum bridge of diffuse emission between the lobes. Taking the spectral index of both lobes to be the same, we model the spectral energy distribution of Fornax A across an energy range spanning eighteen orders of magnitude, to investigate the origin of the X-ray and gamma-ray emission. A standard leptonic model for the production of both the X-rays and gamma-rays by inverse-Compton scattering does not fit the multi-wavelength observations. Our results best support a scenario where the X-rays are produced by inverse-Compton scattering and the gamma-rays are produced primarily by hadronic processes confined to the filamentary structures of the Fornax A lobes.

We present the results of a survey for low frequency radio emission from 17 known exoplanetary systems with the Murchison Widefield Array. This sample includes 13 systems that have not previously been targeted with radio observations. We detected no radio emission at 154 MHz, and put 3 sigma upper limits in the range 15.2-112.5 mJy on this emission. We also searched for circularly polarised emission and made no detections, obtaining 3 sigma upper limits in the range 3.4-49.9 mJy. These are comparable with the best low frequency radio limits in the existing literature and translate to luminosity limits of between 1.2 x 10^14 W and 1.4 x 10^17 W if the emission is assumed to be 100% circularly polarised. These are the first results from a larger program to systematically search for exoplanetary emission with the MWA.

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.

The Murchison Widefield Array (MWA) is a new low frequency interferometric radio telescope, operating in the remote Murchison Radio Observatory in Western Australia. In this paper we present the first MWA observations of the well known radio relics in Abell 3667 (A3667) between 120 and 226 MHz. We clearly detect the radio relics in A3667 and present flux estimates and spectral indices for these features. The average spectral index of the north-west (NW) and south-east (SE) relics is -0.9 +/- 0.1 between 120 and 1400 MHz. We are able to resolve spatial variation in the spectral index of the NW relic from -1.7 to -0.4, which is consistent with results found at higher frequencies.

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

We report on the detection of the millisecond pulsar PSR J0437-4715 with the Murchison Widefield Array (MWA) at a frequency of 192 MHz. Our observations show rapid modulations of pulse intensity in time and frequency that arise from diffractive scintillation effects in the interstellar medium (ISM), as well as prominent drifts of intensity maxima in the time-frequency plane that arise from refractive effects. Our analysis suggests that the scattering screen is located at a distance of $\sim$80-120 pc from the Sun, in disagreement with a recent claim that the screen is closer ($\sim$10 pc). Comparisons with higher frequency data from Parkes reveals a dramatic evolution of the pulse profile with frequency, with the outer conal emission becoming comparable in strength to that from the core and inner conal regions. As well as demonstrating high time resolution science capabilities currently possible with the MWA, our observations underscore the potential to conduct low-frequency investigations of timing-array millisecond pulsars, which may lead to increased sensitivity for the detection of nanoHertz gravitational waves via the accurate characterisation of ISM effects.

The detection of the Epoch of Reionization (EoR) in the redshifted 21-cm line is a challenging task. Here we formulate the detection of the EoR signal using the drift scan strategy. This method potentially has better instrumental stability as compared to the case where a single patch of sky is tracked. We demonstrate that the correlation time between measured visibilities could extend up to 1-2 hr for an interferometer array such as the Murchison Widefield Array (MWA), which has a wide primary beam. We estimate the EoR power based on cross-correlation of visibilities across time and show that the drift scan strategy is capable of the detection of the EoR signal with comparable/better signal-to-noise as compared to the tracking case. We also estimate the visibility correlation for a set of bright point sources and argue that the statistical inhomogeneity of bright point sources might allow their separation from the EoR signal.

Astronomical widefield imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new widefield interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependencies of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarisation correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2-12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square-Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.

The Sun has remained a difficult source to image for radio telescopes, especially at the low radio frequencies. Its morphologically complex emission features span a large range of angular scales, emission mechanisms involved and brightness temperatures. In addition, time and frequency synthesis, the key tool used by most radio interferometers to build up information about the source being imaged is not effective for solar imaging, because many of the features of interest are short lived and change dramatically over small fractional bandwidths. Building on the advances in radio frequency technology, digital signal processing and computing, the kind of instruments needed to simultaneously capture the evolution of solar emission in time, frequency, morphology and polarization over a large spectral span with the requisite imaging fidelity, and time and frequency resolution have only recently begun to appear. Of this class of instruments, the Murchison Widefield Array (MWA) is best suited for solar observations. The MWA has now entered a routine observing phase and here we present some early examples from MWA observations.

We present techniques for bridging the gap between idealized inverse covariance weighted quadratic estimation of 21 cm power spectra and the real-world challenges presented universally by interferometric observation. By carefully evaluating various estimators and adapting our techniques for large but incomplete data sets, we develop a robust power spectrum estimation framework that preserves the so-called "EoR window" and keeps track of estimator errors and covariances. We apply our method to observations from the 32-tile prototype of the Murchinson Widefield Array to demonstrate the importance of a judicious analysis technique. Lastly, we apply our method to investigate the dependence of the clean EoR window on frequency--especially the frequency dependence of the so-called "wedge" feature--and establish upper limits on the power spectrum from z = 6.2 to z = 11.7. Our lowest limit is Delta(k) < 0.3 Kelvin at 95% confidence at a comoving scale k = 0.046 Mpc^-1 and z = 9.5.

We report results from a programme aimed at investigating the temperature of neutral gas in high-redshift damped Lyman-$\alpha$ absorbers (DLAs). This involved (1) HI 21cm absorption studies of a large DLA sample, (2) VLBI studies to measure the low-frequency quasar core fractions, and (3) optical/ultraviolet spectroscopy to determine DLA metallicities and velocity widths. Including literature data, our sample consists of 37 DLAs with estimates of the spin temperature $T_s$ and the covering factor. We find a strong $4\sigma$) difference between the $T_s$ distributions in high-z (z>2.4) and low-z (z<2.4) DLA samples. The high-z sample contains more systems with high $T_s$ values, $\gtrsim 1000$ K. The $T_s$ distributions in DLAs and the Galaxy are also clearly (~$6\sigma$) different, with more high-$T_s$ sightlines in DLAs than in the Milky Way. The high $T_s$ values in the high-z DLAs of our sample arise due to low fractions of the cold neutral medium. For 29 DLAs with metallicity [Z/H] estimates, we confirm the presence of an anti-correlation between $T_s$ and [Z/H], at $3.5\sigma$ significance via a non-parametric Kendall-tau test. This result was obtained with the assumption that the DLA covering factor is equal to the core fraction. Monte Carlo simulations show that the significance of the result is only marginally decreased if the covering factor and the core fraction are uncorrelated, or if there is a random error in the inferred covering factor. We also find evidence for redshift evolution in DLA $T_s$ values even for the z>1 sub-sample. Since z>1 DLAs have angular diameter distances comparable to or larger than those of the background quasars, they have similar efficiency in covering the quasars. Low covering factors in high-z DLAs thus cannot account for the observed redshift evolution in spin temperatures. (Abstract abridged.)

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.

In this paper we explore for the first time the relative magnitudes of three fundamental sources of uncertainty, namely, foreground contamination, thermal noise and sample variance in detecting the HI power spectrum from the Epoch of Reionization (EoR). We derive limits on the sensitivity of a Fourier synthesis telescope to detect EoR based on its array configuration and a statistical representation of images made by the instrument. We use the Murchison Widefield Array (MWA) configuration for our studies. Using a unified framework for estimating signal and noise components in the HI power spectrum, we derive an expression for and estimate the contamination from extragalactic point-like sources in three-dimensional k-space. Sensitivity for EoR HI power spectrum detection is estimated for different observing modes with MWA. With 1000 hours of observing on a single field using the 128-tile MWA, EoR detection is feasible (S/N > 1 for $k\lesssim 0.8$ Mpc$^{-1}$). Bandpass shaping and refinements to the EoR window are found to be effective in containing foreground contamination, which makes the instrument tolerant to imaging errors. We find that for a given observing time, observing many independent fields of view does not offer an advantage over a single field observation when thermal noise dominates over other uncertainties in the derived power spectrum.

We present new wide-field observations of Centaurus A (Cen A) and the surrounding region at 118 MHz with the Murchison Widefield Array (MWA) 32-tile prototype, with which we investigate the spectral-index distribution of Cen A's giant radio lobes. We compare our images to 1.4 GHz maps of Cen A and compute spectral indices using temperature-temperature plots and spectral tomography. We find that the morphologies at 118 MHz and 1.4 GHz match very closely apart from an extra peak in the southern lobe at 118 MHz, which provides tentative evidence for the existence of a southern counterpart to the northern middle lobe of Cen A. Our spatially-averaged spectral indices for both the northern and southern lobes are consistent with previous analyses, however we find significant spatial variation of the spectra across the extent of each lobe. Both the spectral-index distribution and the morphology at low radio frequencies support a scenario of multiple outbursts of activity from the central engine. Our results are consistent with inverse-Compton modelling of radio and gamma-ray data that supports a value for the lobe age of between 10 and 80 Myr.

The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the possibility that the MWA can be used for the purposes of Space Situational Awareness (SSA). In particular we propose that the MWA can be used as an element of a passive radar facility operating in the frequency range 87.5 - 108 MHz (the commercial FM broadcast band). In this scenario the MWA can be considered the receiving element in a bi-static radar configuration, with FM broadcast stations serving as non-cooperative transmitters. The FM broadcasts propagate into space, are reflected off debris in Earth orbit, and are received at the MWA. The imaging capabilities of the MWA can be used to simultaneously detect multiple pieces of space debris, image their positions on the sky as a function of time, and provide tracking data that can be used to determine orbital parameters. Such a capability would be a valuable addition to Australian and global SSA assets, in terms of southern and eastern hemispheric coverage. We provide a feasibility assessment of this proposal, based on simple calculations and electromagnetic simulations that shows the detection of sub-metre size debris should be possible (debris radius of >0.5 m to ~1000 km altitude). We also present a proof-of-concept set of observations that demonstrate the feasibility of the proposal, based on the detection and tracking of the International Space Station via reflected FM broadcast signals originating in south-west Western Australia. These observations broadly validate our calculations and simulations.

We measure the neutral atomic hydrogen (HI) gas content of field galaxies at intermediate redshifts of z ~ 0.1 and z ~ 0.2 using hydrogen 21-cm emission lines observed with the Westerbork Synthesis Radio Telescope (WSRT). In order to make high signal-to-noise ratio detections, an HI signal stacking technique is applied: HI emission spectra from multiple galaxies, optically selected by the CNOC2 redshift survey project, are co-added to measure the average HI mass of galaxies in the two redshift bins. We calculate the cosmic HI gas densities ({\Omega}_{HI}) at the two redshift regimes and compare those with measurements at other redshifts to investigate the global evolution of the HI gas density over cosmic time. From a total of 59 galaxies at z ~ 0.1 we find {\Omega}_{HI} = (0.33 $\pm$ 0.05) ~ $\times$ 10$^{-3}$, and at z ~ 0.2 we find {\Omega}_{HI} = (0.34 $\pm$ 0.09) ~ $\times$ 10$^{-3}$, based on 96 galaxies. These measurements help bridge the gap between high-z damped Lyman-$\alpha$ observations and blind 21-cm surveys at $z=$ 0. We find that our measurements of {\Omega}_{HI} at z ~ 0.1 and 0.2 are consistent with the HI gas density at z ~ 0 and that all measurements of {\Omega}_{HI} from 21-cm emission observations at $z \la$ ~ 0.2 are in agreement with no evolution of the HI gas content in galaxies during the last 2.4 Gyr.

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements (known as tiles) distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper concentrates on the capabilities of the MWA for solar science and summarises some of the solar science results to date, in advance of the initial operation of the final instrument in 2013.

The Murchison Widefield Array is a low frequency (80 - 300 MHz) SKA Precursor, comprising 128 aperture array elements distributed over an area of 3 km diameter. The MWA is located at the extraordinarily radio quiet Murchison Radioastronomy Observatory in the mid-west of Western Australia, the selected home for the Phase 1 and Phase 2 SKA low frequency arrays. The MWA science goals include: 1) detection of fluctuations in the brightness temperature of the diffuse redshifted 21 cm line of neutral hydrogen from the epoch of reionisation; 2) studies of Galactic and extragalactic processes based on deep, confusion-limited surveys of the full sky visible to the array; 3) time domain astrophysics through exploration of the variable radio sky; and 4) solar imaging and characterisation of the heliosphere and ionosphere via propagation effects on background radio source emission. This paper will focus on a brief discussion of the as-built MWA system, highlighting several novel characteristics of the instrument, and a brief progress report (as of June 2012) on the final construction phase. Practical completion of the MWA is expected in November 2012, with commissioning commencing from approximately August 2012 and operations commencing near mid 2013. A brief description of recent science results from the MWA prototype instrument is given.

A new generation of low frequency radio telescopes is seeking to observe the redshifted 21 cm signal from the Epoch of Reionization (EoR), requiring innovative methods of calibration and imaging to overcome the difficulties of widefield low frequency radio interferometry. Precise calibration will be required to separate the small expected EoR signal from the strong foreground emission at the frequencies of interest between 80 and 300 MHz. The Moon may be useful as a calibration source for detection of the EoR signature, as it should have a smooth and predictable thermal spectrum across the frequency band of interest. Initial observations of the Moon with the Murchison Widefield Array 32 tile prototype show that the Moon does exhibit a similar trend to that expected for a cool thermally emitting body in the observed frequency range, but that the spectrum is corrupted by reflected radio emission from Earth. In particular, there is an abrupt increase in the observed flux density of the Moon within the internationally recognised Frequency Modulated (FM) radio band. The observations have implications for future low frequency surveys and EoR detection experiments that will need to take this reflected emission from the Moon into account. The results also allow us to estimate the equivalent isotropic power emitted by the Earth in the FM band and to determine how bright the Earth might appear at metre wavelengths to an observer beyond our own solar system.

Using the final 128 antenna locations of the Murchison Widefield Array (MWA), we calculate its sensitivity to the Epoch of Reionization (EoR) power spectrum of red- shifted 21 cm emission for a fiducial model and provide the tools to calculate the sensitivity for any model. Our calculation takes into account synthesis rotation, chro- matic and asymmetrical baseline effects, and excludes modes that will be contaminated by foreground subtraction. For the fiducial model, the MWA will be capable of a 14{\sigma} detection of the EoR signal with one full season of observation on two fields (900 and 700 hours).

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. Our new method is an extension of A-projection/software-holography/forward modeling analysis techniques and shares their precision deconvolution and widefield polarimetry, while being significantly faster than current implementations that use full direction-dependent antenna gains. Using data from the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology observations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately related to the precision of the data reduction.

The Murchison Widefield Array (MWA) is one of three Square Kilometre Array Precursor telescopes and is located at the Murchison Radio-astronomy Observatory in the Murchison Shire of the mid-west of Western Australia, a location chosen for its extremely low levels of radio frequency interference. The MWA operates at low radio frequencies, 80-300 MHz, with a processed bandwidth of 30.72 MHz for both linear polarisations, and consists of 128 aperture arrays (known as tiles) distributed over a ~3 km diameter area. Novel hybrid hardware/software correlation and a real-time imaging and calibration systems comprise the MWA signal processing backend. In this paper the as-built MWA is described both at a system and sub-system level, the expected performance of the array is presented, and the science goals of the instrument are summarised.

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferometer under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50-degree diameter fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, and used it to produce ~15' angular resolution maps of the two fields. We perform a blind source extraction using these confusion-limited images, and detect 655 sources at high significance with an additional 871 lower significance source candidates. We compare these sources with existing low-frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency surveys in these regions of the sky and with statistical distributions of point sources derived from Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 110 MHz to 200 MHz band.

We have selected a sample of local E+A galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 7 for follow up integral field spectroscopy with the Wide Field Spectrograph (WiFeS) on the ANU 2.3-m telescope. The sample was selected using the Halpha line in place of the [OII]3727 line as the indicator of on-going star formation (or lack thereof). This allowed us to select a lower redshift sample of galaxies than available in the literature since the [OII]3727 falls off the blue end of the wavelength coverage in the SDSS for the very lowest redshift objects. This low redshift selection means that the galaxies have a large angular to physical scale which allows us to resolve the central ~1kpc region of the galaxies; the region where stellar population gradients are expected. Such observations have been difficult to make using other higher redshift samples because even at redshifts z~0.1 the angular to physical scale is similar to the resolution provided by ground based seeing. Our integral field spectroscopy has enabled us to make the first robust detections of Balmer line gradients in the centres of E+A galaxies. Six out of our sample of seven, and all the galaxies with regular morphologies, are observed to have compact and centrally-concentrated Balmer line absorption. This is evidence for compact young cores and stellar population gradients which are predicted from models of mergers and tidal interactions which funnel gas into the galaxy core. Given the generally isolated nature of our sample this argues for the galaxies being seen in the late stage of a merger where the progenitors have already coalesced.

We present the first spectroscopic images of solar radio transients from the prototype for the Murchison Widefield Array (MWA), observed on 2010 March 27. Our observations span the instantaneous frequency band 170.9-201.6 MHz. Though our observing period is characterized as a period of `low' to `medium' activity, one broadband emission feature and numerous short-lived, narrowband, non-thermal emission features are evident. Our data represent a significant advance in low radio frequency solar imaging, enabling us to follow the spatial, spectral, and temporal evolution of events simultaneously and in unprecedented detail. The rich variety of features seen here reaffirms the coronal diagnostic capability of low radio frequency emission and provides an early glimpse of the nature of radio observations that will become available as the next generation of low frequency radio interferometers come on-line over the next few years.

The Murchison Wide-field Array (MWA) is a low frequency radio telescope, currently under construction, intended to search for the spectral signature of the epoch of re-ionisation (EOR) and to probe the structure of the solar corona. Sited in Western Australia, the full MWA will comprise 8192 dipoles grouped into 512 tiles, and be capable of imaging the sky south of 40 degree declination, from 80 MHz to 300 MHz with an instantaneous field of view that is tens of degrees wide and a resolution of a few arcminutes. A 32-station prototype of the MWA has been recently commissioned and a set of observations taken that exercise the whole acquisition and processing pipeline. We present Stokes I, Q, and U images from two ~4 hour integrations of a field 20 degrees wide centered on Pictoris A. These images demonstrate the capacity and stability of a real-time calibration and imaging technique employing the weighted addition of warped snapshots to counter extreme wide field imaging distortions.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.

We report evidence for an anti-correlation between spin temperature $T_s$ and metallicity [Z/H], detected at $3.6 \sigma$ significance in a sample of 26 damped Lyman-$\alpha$ absorbers (DLAs) at redshifts $0.09 < z < 3.45$. The anti-correlation is detected at $3 \sigma$ significance in a sub-sample of 20 DLAs with measured covering factors, implying that it does not stem from low covering factors. We obtain $T_s = (-0.68 \pm 0.17) \times {\rm [Z/H]} + (2.13 \pm 0.21)$ from a linear regression analysis. Our results indicate that the high $T_s$ values found in DLAs do not arise from differences between the optical and radio sightlines, but are likely to reflect the underlying gas temperature distribution. The trend between $T_s$ and [Z/H] can be explained by the larger number of radiation pathways for gas cooling in galaxies with high metal abundances, resulting in a high cold gas fraction, and hence, a low spin temperature. Conversely, low-metallicity galaxies have fewer cooling routes, yielding a larger warm gas fraction and a high $T_s$. Most DLAs at $z>1.7$ have low metallicities, [Z/H] $< -1$, implying that the HI in high-$z$ DLAs is predominantly warm. The anti-correlation between $T_s$ and [Z/H] is consistent with the presence of a mass-metallicity relation in DLAs, suggested by the tight correlation between DLA metallicity and the kinematic widths of metal lines. Most high-$z$ DLAs are likely to arise in galaxies with low masses ($M_{\rm vir} < 10^{10.5} M_\odot$), low metallicities ([Z/H]$< -1$, and low cold gas fractions.

We used observations from the Giant Metrewave Radio Telescope to measure the atomic hydrogen gas content of 324 galaxies around the galaxy cluster Abell 370 at a redshift of z = 0.37 (a look-back time of ~4 billion years). The HI 21-cm emission from these galaxies was measured by coadding their signals using precise optical redshifts obtained with the Anglo-Australian Telescope. The average HI mass measured for all 324 galaxies is (6.6 +- 3.5)x10^9 solar masses, while the average HI mass measured for the 105 optically blue galaxies is (19.0 +- 6.5)x10^9 solar masses. The significant quantities of gas found around Abell 370, suggest that there has been substantial evolution in the gas content of galaxy clusters since redshift z = 0.37. The total amount of HI gas found around Abell 370 is up to ~8 times more than that seen around the Coma cluster, a nearby galaxy cluster of similar size. Despite this higher gas content, Abell 370 shows the same trend as nearby clusters, that galaxies close to the cluster core have lower HI gas content than galaxies further away. The Abell 370 galaxies have HI mass to optical light ratios similar to local galaxy samples and have the same correlation between their star formation rate and HI mass as found in nearby galaxies. The average star formation rate derived from [OII] emission and from de-redshifted 1.4 GHz radio continuum for the Abell 370 galaxies also follows the correlation found in the local universe. The large amounts of HI gas found around the cluster can easily be consumed by the observed star formation rate in the galaxies over the ~4 billion years (from z = 0.37) to the present day.

Modern graphics processing units (GPUs) are inexpensive commodity hardware that offer Tflop/s theoretical computing capacity. GPUs are well suited to many compute-intensive tasks including digital signal processing. We describe the implementation and performance of a GPU-based digital correlator for radio astronomy. The correlator is implemented using the NVIDIA CUDA development environment. We evaluate three design options on two generations of NVIDIA hardware. The different designs utilize the internal registers, shared memory and multiprocessors in different ways. We find that optimal performance is achieved with the design that minimizes global memory reads on recent generations of hardware. The GPU-based correlator outperforms a single-threaded CPU equivalent by a factor of 60 for a 32 antenna array, and runs on commodity PC hardware. The extra compute capability provided by the GPU maximises the correlation capability of a PC while retaining the fast development time associated with using standard hardware, networking and programming languages. In this way, a GPU-based correlation system represents a middle ground in design space between high performance, custom built hardware and pure CPU-based software correlation. The correlator was deployed at the Murchison Widefield Array 32 antenna prototype system where it ran in real-time for extended periods. We briefly describe the data capture, streaming and correlation system for the prototype array.

We have used the GMOS instrument on the Gemini-South telescope to obtain spatially-resolved two-colour imaging and IFU spectroscopy of a sample of ten nearby E+A galaxies. Surface brightness profiles measured using our imaging data show the isophotal profiles of our sample are generally r^{1/4}-like, consistent with a sample dominated by early-type galaxies. This is further underscored by all galaxies having early-type morphological classifications, and showing a behaviour in the central velocity dispersion-absolute magnitude plane that is consistent with the Faber-Jackson relation, once the transitory brightening that occurs in the E+A phase is corrected for. In addition, two-thirds of our sample shows clear evidence of either ongoing or recent tidal interactions/mergers. While all the galaxies in our sample have total integrated colours that are relatively blue, they show a diversity of colour gradients, possessing central core regions that are either redder, bluer, or indistinct in colour relative to their outer regions. Kinematically, the most striking property is the significant and unambiguous rotation that is seen in all our E+A galaxies, with it being generally aligned close to the photometric major axis. This is contrary to the findings of Norton et al. 2001, who found little or no evidence for rotation in a very similar sample of nearby E+A galaxies. We also clearly demonstrate that our E+A galaxies are, in all but one case, consistent with being "fast rotators" (Emsellem et al. 2007). We argue that the combination of disturbed morphologies and significant rotation in these galaxies supports their production via gas-rich galaxy mergers and interactions. The large fraction of fast rotators argues against equal mass mergers being the dominant progenitor to the E+A population.

We have used the Very Long Baseline Array to image 18 quasars with foreground damped Lyman-$\alpha$ systems (DLAs) at 327, 610 or 1420 MHz, to measure the covering factor $f$ of each DLA at or near its redshifted HI 21cm line frequency. Including six systems from the literature, we find that none of 24 DLAs at $0.09 < z < 3.45$ has an exceptionally low covering factor, with $f \sim 0.45 - 1$ for the 14 DLAs at $z > 1.5$, $f \sim 0.41 - 1$ for the 10 systems at $z < 1$, and consistent covering factor distributions in the two sub-samples. The observed paucity of detections of HI 21cm absorption in high-$z$ DLAs thus cannot be explained by low covering factors and is instead likely to arise due to a larger fraction of warm HI in these absorbers.

The Murchison Widefield Array (MWA) is a dipole-based aperture array synthesis telescope designed to operate in the 80-300 MHz frequency range. It is capable of a wide range of science investigations, but is initially focused on three key science projects. These are detection and characterization of 3-dimensional brightness temperature fluctuations in the 21cm line of neutral hydrogen during the Epoch of Reionization (EoR) at redshifts from 6 to 10, solar imaging and remote sensing of the inner heliosphere via propagation effects on signals from distant background sources,and high-sensitivity exploration of the variable radio sky. The array design features 8192 dual-polarization broad-band active dipoles, arranged into 512 tiles comprising 16 dipoles each. The tiles are quasi-randomly distributed over an aperture 1.5km in diameter, with a small number of outliers extending to 3km. All tile-tile baselines are correlated in custom FPGA-based hardware, yielding a Nyquist-sampled instantaneous monochromatic uv coverage and unprecedented point spread function (PSF) quality. The correlated data are calibrated in real time using novel position-dependent self-calibration algorithms. The array is located in the Murchison region of outback Western Australia. This region is characterized by extremely low population density and a superbly radio-quiet environment,allowing full exploitation of the instrumental capabilities.

We report on the detection of giant pulses from the Crab Nebula pulsar at a frequency of 200 MHz using the field deployment system designed for the Mileura Widefield Array's Low Frequency Demonstrator (MWA-LFD). Our observations are among the first high-quality detections at such low frequencies. The measured pulse shapes are deconvolved for interstellar pulse broadening, yielding a pulse-broadening time of 670$\pm$100 $\mu$s, and the implied strength of scattering (scattering measure) is the lowest that is estimated towards the Crab nebula from observations made so far. The sensitivity of the system is largely dictated by the sky background, and our simple equipment is capable of detecting pulses that are brighter than $\sim$9 kJy in amplitude. The brightest giant pulse detected in our data has a peak amplitude of $\sim$50 kJy, and the implied brightness temperature is $10^{31.6}$ K. We discuss the giant pulse detection prospects with the full MWA-LFD system. With a sensitivity over two orders of magnitude larger than the prototype equipment, the full system will be capable of detecting such bright giant pulses out to a wide range of Galactic distances; from $\sim$8 to $\sim$30 kpc depending on the frequency. The MWA-LFD will thus be a highly promising instrument for the studies of giant pulses and other fast radio transients at low frequencies.

We use observations from the Giant Metrewave Radio Telescope (GMRT) to measure the atomic hydrogen gas content of star-forming galaxies at z = 0.24 (i.e. a look-backtime of ~3 Gyr). The sample of galaxies studied were selected from Halpha-emitting field galaxies detected in a narrow-band imaging survey with the Subaru Telescope. The Anglo-Australian Telescope was used to obtain precise optical redshifts for these galaxies. We then coadded the HI 21 cm emission signal for all the galaxies within the GMRT spectral line data cube. From the coadded signal of 121 galaxies, we measure an average atomic hydrogen gas mass of (2.26 +- 0.90)*10^9 solar masses. We translate this HI signal into a cosmic density of neutral gas at z = 0.24 of Omega_gas = (0.91 +- 0.42)*10^-3. This is the current highest redshift at which Omega_gas has been constrained from 21 cm emission and our value is consistent with that estimated from damped Lyman-alpha systems around this redshift. We also find that the correlations between the Halpha luminosity and the radio continuum luminosity and between the star formation rate and the HI gas content in star-forming galaxies at z = 0.24 are consistent with the correlations found at z = 0. These two results suggest that the star formation mechanisms in field galaxies ~3 Gyr ago were not substantially different from the present, even though the star formation rate is 3 times higher.

Experiments were performed with prototype antenna tiles for the Mileura Widefield Array--Low Frequency Demonstrator (MWA-LFD) to better understand the widefield, wideband properties of their design and to characterize the radio frequency interference (RFI) between 80 and 300 MHz at the site in Western Australia. Observations acquired during the six month deployment confirmed the predicted sensitivity of the antennas, sky-noise dominated system temperatures, and phase-coherent interferometric measurements. The radio spectrum is remarkably free of strong terrestrial signals, with the exception of two narrow frequency bands allocated to satellite downlinks and rare bursts due to ground-based transmissions being scattered from aircraft and meteor trails. Results indicate the potential of the MWA-LFD to make significant achievements in its three key science objectives: epoch of reionziation science, heliospheric science, and radio transient detection.

We present the results of global VLBI observations of the BL Lac object PKS 0003-066 at 2.3 and 8.6 GHz, over a 10 year period using data from new and archived VLBA, geodetic VLBI, and global VLBI (including the use of new disk-based recording and software correlation systems) observations. Inter-continental baselines resolve this source into a compact core and three jet components that move away from the core with proper motions ranging between \~0.1 mas/yr for the outer jet component and ~1.4 mas/yr for the fastest inner jet component, i.e. apparent speeds of ~2.0 c to ~21.6 c respectively. This result is in contrast to previous studies of this source which did not reveal significant jet motions from limited data sets. We find that the fast inner jet components catch up to the more slowly moving outer jet components and interact strongly, causing local brightening of the jet in the interaction. We document the body of VLBI data for this source, as well as other supporting radio observations, as a basis for future work that will investigate the underlying physics of the pc-scale jet in PKS 0003-066.

Observations of the high-redshift Universe with the 21 cm hyperfine line of neutral hydrogen promise to open an entirely new window onto the early phases of cosmic structure formation. Here we review the physics of the 21 cm transition, focusing on processes relevant at high redshifts, and describe the insights to be gained from such observations. These include measuring the matter power spectrum at z~50, observing the formation of the cosmic web and the first luminous sources, and mapping the reionization of the intergalactic medium. The epoch of reionization is of particular interest, because large HII regions will seed substantial fluctuations in the 21 cm background. We also discuss the experimental challenges involved in detecting this signal, with an emphasis on the Galactic and extragalactic foregrounds. These increase rapidly toward low frequencies and are especially severe for the highest redshift applications. Assuming that these difficulties can be overcome, the redshifted 21 cm line will offer unique insight into the high-redshift Universe, complementing other probes but providing the only direct, three-dimensional view of structure formation from z~200 to z~6.

The Westerbork Synthesis Radio Telescope finds a weak 21cm line emission feature at the coordinates (RA-Dec-velocity) of the sub-Damped Lyman-alpha absorber observed at z_abs=0.00632 in the spectrum of PG1216+069. The emission feature, WSRT-J121921+0639, lies within 30" of the quasar sightline, is detected at 99.8% (3 sigma) confidence level, has M_HI between 5 and 15x10^6 M_solar, and has velocity spread between 20 and 60 km/s. Other HI emitters in the field include VCC297 at a projected distance of 86/h_75 kpc and a previously unreported HI cloud, WSRT-J121919+0624 at 112/h_75 kpc with M_HI ~ 3x10^8 M_solar. The optically identified, foreground galaxy that is closest to the quasar sightline appears to be VCC339 (~L*/25) at 29/h_75 kpc and velocity offset 292 km/s . A low surface brightness galaxy with the HI mass of the sub-DLA absorber WSRT-J121921+0639 would likely have m_B ~ 17, and its diffuse optical emission would need to compete with the light of both the background QSO and a brighter foreground star ~10" from the QSO sight line.

We have detected the four 18cm OH lines from the $z \sim 0.765$ gravitational lens toward PMN J0134-0931. The 1612 and 1720 MHz lines are in conjugate absorption and emission, providing a laboratory to test the evolution of fundamental constants over a large lookback time. We compare the HI and OH main line absorption redshifts of the different components in the $z \sim 0.765$ absorber and the $z \sim 0.685$ lens toward B0218+357 to place stringent constraints on changes in $F \equiv g_p [\alpha^2/\mu]^{1.57}$. We obtain $[\Delta F/F] = (0.44 \pm 0.36^{\rm stat} \pm 1.0^{\rm syst}) \times 10^{-5}$, consistent with no evolution over the redshift range $0 < z < 0.7$. The measurements have a $2 \sigma$ sensitivity of $[\Delta \alpha/\alpha] < 6.7 \times 10^{-6}$ or $[\Delta \mu/\mu] < 1.4 \times 10^{-5}$ to fractional changes in $\alpha$ and $\mu$ over a period of $\sim 6.5$ Gyr, half the age of the Universe. These are among the most sensitive current constraints on changes in $\mu$.

A comprehensive analysis of 355 high-quality WSRT HI 21-cm line maps of nearby galaxies shows that the properties and incident rate of Damped Lyman-alpha (DLA) absorption systems observed in the spectra of high redshift QSOs are in good agreement with DLAs originating in gas disks of galaxies like those in the z~0 population. Comparison of low-z DLA statistics with the HI incidence rate and column density distribution f(N) for the local galaxy sample shows no evidence for evolution in the integral "cross section density" below z~1.5, implying that there is no need for a hidden population of galaxies or HI clouds to contribute significantly to the DLA cross section. Compared with z~4, our data indicates evolution of a factor of two in the comoving density along a line of sight. We find that dN/dz(z=0)=0.045 +/- 0.006. The idea that the local galaxy population can explain the DLAs is further strengthened by comparing the properties of DLAs and DLA galaxies with the expectations based on our analysis of local galaxies. The distribution of luminosities of DLA host galaxies, and of impact parameters between QSOs and the centres of DLA galaxies, are in good agreement with what is expected from local galaxies. Approximately 87% of low z DLA galaxies are expected to be fainter than L* and 37 per cent have impact parameters less than 1'' at z=0.5. The analysis shows that some host galaxies with very low impact parameters and low luminosities are expected to be missed in optical follow up surveys. The well-known metallicity-luminosity relation in galaxies, in combination with metallicity gradients in galaxy disks, cause the expected median metallicity of low redshift DLAs to be low (~1/7 solar), which is also in good agreement with observations of low z DLAs. (Abridged)

We calculate in detail the expected properties of low redshift DLA systems under the assumption that they arise in the gaseous disks of galaxies like those in the z=0 population. A sample of 355 nearby galaxies is analysed, for which high quality HI 21-cm emission line maps are available as part of an extensive survey with the Westerbork telescope (WHISP). We find that expected luminosities, impact parameters between quasars and DLA host galaxies, and metal abundances are in good agreement with the observed properties of DLAs and DLA galaxies. The measured redshift number density of z=0 gas above the DLA limit is dN/dz=0.045 +/- 0.006, which compared to higher redshift measurements implies that there is no evolution in the comoving density of DLAs along a line of sight between z=1.5 and z=0, and a decrease of only a factor of two from z=4 to the present time. We conclude that the local galaxy population can explain all properties of low redshift DLAs.

This overview provides an interface between lines of thought on rfi mitigation in the fields of radio astronomy and signal processing. The goal is to explore the commonality of different approaches to help researchers in both fields interpret each others' concepts and jargon. The paper elaborates on the astronomers' concept of gain closure relations and how they may be used in a "self-calibrating" system of rfi cancellation. Further discussion of the eigen decomposition method in terms of rfi power and antenna gains introduces adaptive nulling and rfi cancellation through matrix partitioning. While multipath scattering appears at first glance to be fatal to methods of rfi cancellation, its effect is easily incorporated in frequency dependent gain coefficients under many circumstances.

Neutral intergalactic clouds are so greatly out numbered by galaxies that their integral HI content is negligible in comparison to that contained in optically luminous galaxies. In fact, no HI cloud that is not associated with a galaxy or grouping of galaxies has yet been identified. This points to a causal relationship that relies on gravitational potentials that bind galaxies also being responsible for confining HI clouds to sufficient density that they can become self-shielding to the ionizing background radiation. Unconfined clouds of low density become ionized, but confined clouds find themselves vulnerable to instability and collapse, leading to star formation.

We report on our experiences using the Metsahovi Radio Observatory's (MRO) VLBI Standard Interface (VSI, Whitney 2002) recorder in a number of astronomical applications. The PC-EVN device is a direct memory access (DMA) interface which allows 512 megabit per second (Mbps) or better recording to "off the shelf" PC components. We have used this setup to record at 640 Mbps for a pulsar coherent dispersion system and at 256 Mpbs for a global VLBI session. We have also demonstrated recording at 512 Mbps and will shortly form cross correlations between the CPSR-II and the PC-EVN systems.

We discuss the prospects for ``tomography'' of the intergalactic medium (IGM) at high redshifts using the 21 cm transition of neutral hydrogen. Existing observational constraints on the epoch of reionization imply a complex ionization history that may require multiple generations of sources. The 21 cm transition provides a unique tool to probe this era in detail, because it does not suffer from saturation effects, retains full redshift information, and directly probes the IGM gas. Observations in the redshifted 21cm line will allow one to study the history and morphology of reionization in detail. Depending on the characteristics of the first sources, they may also allow us to probe the era before reionization, when the first structures and luminous sources were forming. The construction of high signal-to-noise ratio maps on arcminute scales will require approximately one square kilometer of collecting area.

The epoch of reionization (EoR) sets a fundamental benchmark in cosmic structure formation, corresponding to the formation of the first luminous objects that act to ionize the neutral intergalactic medium (IGM). Recent observations at near-IR and radio wavelengths imply that we are finally probing into this key epoch of galaxy formation at $z \ge 6$. The Square Kilometer Array (SKA) will provide critical insight into the EoR in a number of ways. First, the ability of the SKA to image the neutral IGM in 21cm emission is a truly unique probe of the process of reionization, and is recognized as the next necessary and fundamental step in our study of the evolution of large scale structure and cosmic reionization. Second, study of HI 21cm absorption toward the first radio loud objects probes small to intermediate scale structure in the neutral 'cosmic web', as well as HI in the first collapsed structures (proto-disks and mini-halos). And third, the incomparable sensitivity of the SKA allows for the study of the molecular gas, dust, and star formation activity in the first galaxies, as well as the radio continuum emission from the first accreting massive black holes. Such objects will be obscured at optical wavelengths due to absorption by the neutral IGM.

HI 21~cm absorption spectroscopy provides an excellent probe of the neutral gas content of absorbing galaxies, yielding information on their kinematics, mass, physical size and ISM conditions. The high sensitivity, unrivaled frequency coverage and RFI suppression techniques of the SKA will enable it to use HI absorption to study the ISM of high column density intervening systems along thousands of lines of sight out to high redshifts. Blind SKA 21~cm surveys will yield large, unbiased absorber samples, tracing the evolution of normal galaxies and active galactic nuclei from $z \gtrsim 6$ to the present epoch. It will thus be possible to directly measure the physical size and mass of typical galaxies as a function of redshift and, hence, to test hierarchical models of structure formation.

We have used the Westerbork Synthesis Radio Telescope to detect HI 21cm absorption at $z \sim 0.7645$ in the gravitational lens system towards PMN J0134-0931. The 21cm profile has two broad components, with peak optical depths of $0.047 \pm 0.007$ and $0.039 \pm 0.007$, at heliocentric redshifts $0.76470 \pm 0.00006$ and $0.76348 \pm 0.00006$, respectively. The redshift of the stronger component matches that of CaII H and K absorption detected earlier. The absorption has a total velocity width of $\sim 500$ km/s (between nulls) and an equivalent width of $7.1 \pm 0.08$ km/s. This would imply a total HI column density of $2.6 \pm 0.3 \times 10^{21}$ per cm$^2$, for a spin temperature of 200 K and a covering factor of unity. The high estimated HI column density is consistent with the presence of large amounts of dust at the lens redshift; the intervening dust could be responsible for the extremely red colour of the background quasar.

We present results of two 21-cm HI surveys performed with the Australia Telescope Compact Array in the nearby Centaurus A and Sculptor galaxy groups. These surveys are sensitive to compact HI clouds and galaxies with HI masses as low as 3E+06 Msun, and are therefore among the most sensitive extragalactic HI surveys to date. The surveys consist of sparsely spaced pointings that sample approximately 2% of the groups' area on the sky. We detected previously known group members, but we found no new HI clouds or galaxies down to the sensitivity limit of the surveys. If the HI mass function had a faint end slope of alpha = 1.5 below M_{HI} = 10^{7.5} Msun in these groups, we would have expected ~3 new objects. Cold dark matter theories of galaxy formation predict the existence of a large number low mass DM sub-halos that might appear as tiny satellites in galaxy groups. Our results support and extend similar conclusions derived from previous HI surveys that a HI rich population of these satellites does not exist.

Westerbork Radio Synthesis Telescope observation of the galaxy VV29=Arp188=UGC10214 shows that there are at least three distinct dynamical components whose kinematics can be traced in 21cm line emission. The system appears to be the result of a galaxy-galaxy interaction. We identify a sufficient number of dynamical elements containing baryons (stars and neutral gas) that there is no compelling reason to postulate the presence of an additional dark matter halo that is devoid of detectable baryons. The central galaxy VV29a is massive (V_rot = 330 km/s) and gas rich (M_HI} = 6x10^9 Msolar). The distinctive optical plume (VV29b), which extends eastward from the main galaxy, is also gas rich (M_HI = 3x10^9 Msolar) and has a very low gradient in line of sight velocity (<30 km/s) over 70kpc. On the western side, there is an HI feature of M_HI = 4x10^8 Msolar that participates strongly in orbital motion about the host in the same sense of rotation as the VV29a itself. A blue, less massive, gas-rich galaxy "VV29c" (M_HI = 9x10^8 Msolar) appears clearly in the HI maps as an 170 km/s wide spectral feature, seen in projection against or, more likely, behind the west side of the host disk. Its high recessional velocity is counter to the host rotation direction. The optical images of Trentham et al (2001) show signs of this blue dwarf against the redder VV29a disk. The companion galaxy CGCG27-021=MGC09-26-54 (at projected distance 115 kpc) is not detected in 21cm line emission (M_HI<10^9Msolar).

As part of a larger program to identify low redshift radio analogues of the damped Lyman-alpha (DLA) absorbers seen in the spectra of high redshift quasars, Westerbork Synthesis Radio Telescope (WSRT) observations have discovered two new HI 21cm absorption lines at z=0.394 and z=0.437 in the spectra of the radio sources B 0248+430 and B 1243-072 respectively. These sightlines and redshifts were selected for study on the basis of the previously known low ionization absorption lines of MgII, and neither has been observed in the Lyman-alpha line. The 21cm line observations provide information on column densities, temperatures and kinematics of the thickest cold neutral clouds in the absorbers.

Blind 21cm surveys in the local universe have shown that the local HI mass density, Omega_HI, is dominated by luminous, high surface brightness, spiral galaxies. On the other hand, surveys for host galaxies of damped Ly-alpha systems have not always been successful in finding bright spiral galaxies. From an analysis of 21cm aperture synthesis maps of nearby galaxies we show that this apparent contradiction can be resolved by realizing that the HI mass density is dominated by L* galaxies, but the HI cross section near the DLA threshold is more evenly distributed over galaxies with a large range in luminosity, gas mass, and surface brightness. The distributions of column densities and impact parameters of optically identified and non-identified DLA host galaxies in the literature and the HI maps are qualitatively in agreement. Due to poor number statistics of low redshift DLA systems, there is no firm indication that the redshift number count of low redshift DLA systems is inconsistent with that calculated from the nearby galaxy population.

We measure the z=0 B-band optical luminosity function (LF) for galaxies selected in a blind HI survey. The total LF of the HI selected sample is flat, with Schechter parameters M*=-19.38_{-0.62}^{+1.02} + 5 log h mag and alpha=-1.03_{-0.15}^{+0.25}, in good agreement with LFs of optically selected late-type galaxies. Bivariate distribution functions of several galaxy parameters show that the HI density in the local Universe is more widely spread over galaxies of different size, central surface brightness, and luminosity than is the optical luminosity density. The number density of very low surface brightness (>24.0 mag/arcsec^2) gas-rich galaxies is considerably lower than that found in optical surveys designed to detect dim galaxies. This suggests that only a part of the population of LSB galaxies is gas rich and that the rest must be gas poor. However, we show that this gas-poor population must be cosmologically insignificant in baryon content. The contribution of gas-rich LSB galaxies (>23.0 mag/arcsec^2) to the local cosmological gas and luminosity density is modest (18_{-5}^{+6} and 5_{-2}^{+2} per cent respectively); their contribution to Omega_matter is not well-determined, but probably < 11 per cent. These values are in excellent agreement with the low redshift results from the Hubble Deep Field.

Westerbork Radio Synthesis Telescope observations of the redshifted 21cm line absorber against the z_em = 0.871 double lobed quasar 3C196 show that the intervening absorber is an L approximately equal L* spiral galaxy (3C196-G1) and that the absorbing layer of cold gas extends to radii of at least 30 h_50^-1 kpc. The new data solve several long standing puzzles about this system by (1) discovering a second 21cm absorption feature, corresponding to absorption against the NE lobe of the background radio source and (2) spatially ``resolving'' the two absorption features to isolate the absorption along the two lines of sight to the opposing radio lobes. These findings resolve the disagreement in redshift between the UV metal and 21cm lines, and as well as demonstrating that the neutral layer does absorb both lobes of the background radio source. Simple kinematic models with an inclined, rotating gas disk match the observed 21cm profile and are also compatible with both the redshift and velocity spread of the absorption measured in UV resonance lines along a third, independent line of sight to the quasar nucleus and with the lack of 21cm absorption in as earlier VLBI experiment that was sensitive to opacity against the hot spot in the northern lobe. The inferred rotation speed and luminosity for the galaxy are compatible with the z approximately 0 Tully-Fisher Relation. This system illustrates well how 21cm absorption against extended background radio sources is a powerful tool in determining the nature of the damped Lyman-alpha class of QSO absorption line system.

We present measurements of absorption by the 21cm hyperfine transition of neutral hydrogen toward radio sources at substantial look-back times. These data are used in combination with observations of rotational transitions of common interstellar molecules, to set limits on the evolution of the fine structure constant: (dot alpha)/alpha < 3.5e-15/yr to a look-back time of 4.8 Gyr. The neutral hydrogen observations employed Very Long Baseline Interferometry in order to mitigate the substantial uncertainty arising from the fact that observations at very different wavelengths may probe different lines-of-site due to frequency dependent structure of the background source. We discuss the implication of these results on theories unifying natural forces based on compact quantum dimensions. In the context of string theory, the limit on the secular evolution of the scale factor of these compact dimensions, R, is (dot R)/R < 1e-15/yr. Including terrestrial and other astronomical measurements places limits (2sigma) on slow oscillations of R from the present to the epoch of cosmic nucleosynthesis, just seconds after the big bang, of (Delta R)/R < 1e-5.

NGC 2146 is a peculiar spiral galaxy which is currently undergoing a major burst of star formation and is immersed in a extended HI structure that has morphological and kinematical resemblence to a strong tidal interaction. This paper reports aperture synthesis observations carried out in the 21cm line with the Very Large Array (VLA - The National Radio Astronomy Observatory (NRAO) is operated by Associated Universities, Inc. under cooperative agreement with the National Science Foundation.) of two fields positioned to optimally cover the HI streams to the north and south of the galaxy, along with a 300 ft total power spectral mapping program to recover the low surface brightness extended emission. The observations reveal elongated streams of neutral hydrogen towards both the north and the south of the optical galaxy extending out up to 6 Holmberg radii. The streams are not in the principle plane of rotation of the galaxy, but instead are suggestive of a tidal interaction between NGC 2146 and a LSB companion that was destroyed by the encounter and remains undetected at optical wavelengths. Part of the southern stream is turning back to fall into the main galaxy, where it will create a long-lived warp in the HI disk of NGC 2146. Analysis of the trajectory of the outlying gas suggests that the closest encounter took place about 0.8 billion years ago and that infall of debris will continue for a similar time span.

The growing level of radio frequency interference (RFI) is a recognized problem for research in radio astronomy. This paper describes an intuitive but powerful RFI cancellation technique that is suitable for radio spectroscopy where time-averages are recorded. An RFI ``reference signal,'' is constructed from the cross power spectrum of the signals from the two polarizations of a reference horn pointed at the source of the RFI signal. The RFI signal paths obey simple phase and amplitude closure relations, which allows computation of the RFI contamination in the astronomical data and the corrections to be applied to the astronomical spectra. Since the method is immune to the effects of multipath scattering in both the astronomy and reference signal channels, ``clean copies'' of the RFI signal are not required. The method could be generalized (1) to interferometer arrays, (2) to correct for scattered solar radiation that causes spectral ``standing waves'' in single-dish spectroscopy, and (3) to pulsar survey and timing applications where a digital correlator plays an important role in broadband pulse dedispersion. Future large radio telescopes, such as the proposed LOFAR and SKA arrays, will require a high degree of RFI suppression and could implement the technique proposed here with the benefit of faster electronics, greater digital precision and higher data rates.

Neutral hydrogen traces gravitational potentials. In the nearby universe, 21cm emission-line surveys show that the bulk of the HI resides in well-formed, optically-luminous galaxies. At high redshift, 21cm line absorption against background radio quasars occurs in gas-rich systems identified with the highest HI column densities -- the ``damped Lyman alpha'' quasar absorption-line systems. High spatial-resolution observations of the redshifted 21cm line absorbers measure sizes and kinematics of the neutral absorbers.

The Arecibo HI Strip Survey probed the halos of ~300 cataloged galaxies and the environments of ~14 groups with sensitivity to neutral hydrogen masses > 10^7 M_sun. The survey detected no objects with properties resembling the High Velocity Clouds (HVCs) associated with the Milky Way or Local Group. If the HVCs were typically M_HI=10^{7.5} M_sun objects distributed throughout groups and galaxy halos at distances of 1 Mpc, the survey should have made ~70 HVC detections in groups and ~250 detections around galaxies. The null detection implies that HVCs are deployed at typical distances of <200 kpc from the galaxies or group barycenters. If the clouds are in virial equilibrium, their average dark matter fraction must be 98% or higher.

We present a detailed study of the HI 21cm absorption system at z=0.0912 towards the radio quasar B0738+313. The uncommonly narrow main absorption line and weak secondary line are resolved for the first time. In addition we find it necessary to add a third, broader shallow component to obtain a good fit to the spectrum. Although the harmonic mean spin temperature calculated by comparison of the 21cm lines to the damped Ly-alpha line is T_s = 775 K, the thermal kinetic temperatures of the two narrow components, calculated from their widths, are much lower: T_k \leq 297 and \leq 103 K respectively. This is the first case of a redshifted absorption system for which T_k is measured to be less than T_s. We discuss this result in the context of a two phase gas model, in which the damped Ly-alpha gas is sensitive to a significant neutral column density of warm phase gas as well as the cold phase gas of the narrow 21cm lines. The third component is interpreted as representing the warm phase gas with with T_k \leq 5050 K. The combined column density of the three 21cm components is approximately equal to that derived from fits to the damped Ly-alpha line.

A Square Kilometer Array radio telescope will detect tens of thousands of galaxies per square degree in the 21cm emission line of neutral hydrogen. The telescope will be sensitive to ordinary galaxy populations at redshifts z>3 when the mass density in neutral gas greatly exceeded the mass in luminous stars. Spectroscopy in the HI line will trace the kinematics of evolving galactic potentials, as well as monitoring consumption of the fuel for star formation. The OH megamaser population may be a significant indicator of merger and interaction rates for galaxies over cosmic time. Spectroscopy of the redshifted HI and OH lines provides a major incentive for designing a large radio telescope with continuous frequency coverage below 2 GHz.

We present a measurement of the HI column density distribution function, f(N), at the present epoch for column densities log N > 20 cm^-2. These high column densities compare to those measured in damped Ly alpha lines seen in absorption against background quasars. Although observationally rare, it appears that the bulk of the neutral gas in the Universe is associated with these damped Ly alpha systems. In order to obtain a good anchor point at z=0 we determine f(N) in the local Universe by using 21cm synthesis observations of a complete sample of spiral galaxies. We show that f(N) for damped Ly alpha systems has changed significantly from high z to the present and that change is greatest for the highest column densities. The measurements indicate that low surface brightness galaxies make a minor contribution to the cross section for HI, especially for log N > 21^-2.

Surveys in the 21cm line of neutral hydrogen are testing the completeness of the catalogs of nearby galaxies. The remarkable observational fact is that the potential wells that confine gas to sufficient density that it can remain neutral in the face of ionizing radiation also provide sites for star formation, so that there are no known cases of neutral intergalactic clouds without associated star light.

A radio spectroscopic driftscan survey in the 21cm line with the Nancay Radio Telescope of 0.08 steradians of sky in the direction of the constellation Canes Venatici covering a heliocentric velocity range of -350 < V_hel < 2350 km/s produced 53 spectral features, which was further reduced to a sample of 33 reliably detected galaxies by extensive follow-up observations. With a typical noise level of rms = 10 mJy after Hanning smoothing, the survey is -- depending on where the detections are located with regard to the center of the beam -- sensitive to M(HI)=1-2x10^8 Msun at 23 Mpc and to M(HI)=4-8x10^7 Msun throughout the CVn Groups (Ho=100). The survey region had been previously examined on deep optical plates by Binggeli et al. 1990 and contains loose groups with many gas-rich galaxies as well as voids. No galaxies that had not been previously identified in these deep optical surveys were uncovered in our HI survey, neither in the groups nor the voids. The implication is that no substantial quantity of neutral hydrogen contained in gas-rich galaxies has been missed in these well-studied groups. All late-type members of our sample are listed in the Fisher and Tully (1981b) optically selected sample of nearby late-type galaxies; the only system not contained in Fisher and Tullys' Catalog is the S0 galaxy NGC 4203. Within the well--sampled CVn group volume with distances corrected for flow motions, the HI mass function is best fitted with the Zwaan et al. (1997) HI mass function (alpha=-1.2) scaled by a factor of f=4.5 in account of the locally overdense region.

Bright OH megamaser galaxies, radiating 1667/1665 MHz lines, could be detected at redshifts from z from 1 to 3 in moderate integration times with existing radio telescopes. The superluminous FIR galaxies that host the megamasers are relatively rare at z equal 0, but they may have been more common at high redshift, if the galaxy merger rate increases steeply with redshift. Therefore, blind radio spectroscopic surveys at frequencies of 400 to 1000 Mhz can form an independent test of the galaxy merger rate as a function of time over the redshift interval z = 4 to 0.7. The redshift range z = 0.17 to 0.4 will be difficult to survey for OH masers, since spectroscopic survey signals will be confused with HI emission from normal galaxies at redshifts less than 0.3. In fact, the signals from OH masers are likely to dominate over 21cm line emission from normal galaxies at frequencies below 1200 MHz (i.e. large redshifts z_HI > 0.18 and z_OH > 0.4). Surveyors of nearby galaxies in the 21cm line may find that OH masers form a contaminant to deep, blind HI surveys for redshift velocities less than a few hundred kilometers per second. At frequencies just above 1420 MHz, sensitive sky surveys might detect OH masers, which could be mistaken for a population of ``infalling, compact High Velocity Clouds'' but would ultimately be traced to luminous FIR background galaxies at z equal 0.17 once optical and IR follow-up has been performed.

The present generation of radio telescopes, combined with powerful new spectrometers, is opening a new age of redshifted radio absorption-line studies. Outfitting of arrays of antennas, such as the European VLBI Network and the upgraded VLA, with flexibly tuned receivers, will measure sizes and kinematics of intervening galaxies as a function of cosmic time.

We report the discovery of two low redshift HI 21cm absorbers, one at z = 0.2212 towards the z_{em} = 0.630 quasar OI 363 (B0738+313), and the other at z = 0.3127 towards PKS B1127-145 (z_{em} = 1.187). Both were found during a survey of MgII selected systems at redshifts 0.2 < z < 1 using the new UHF-high system at the Westerbork Synthesis Radio Telescope (WSRT). New HST/FOS observations also identify both systems as damped Ly-alpha (DLa) absorbers. By comparing the column density from the DLa line with that from the HI 21cm line, we calculate the spin temperature, and find that T_s is approximately 1000 K for both of these low redshift absorbers. We briefly discuss some implications of this result.

The HI mass function (HIMF) for galaxies in the local universe is constructed from the results of the Arecibo HI Strip Survey, a blind extragalactic survey in the 21cm line. The survey consists of two strips covering in total 65 square degrees of sky, with a depth of cz = 7400 km/s and was optimized to detect column densities of neutral gas N_HI > 10^18 cm^-2 (5 sigma). The survey yielded 66 significant extragalactic signals of which approximately 50% are cataloged galaxies. No free floating HI clouds without stars are found. VLA follow-up observations of all signals have been used to obtain better measurements of the positions and fluxes and allow an alternate determination of the achieved survey sensitivity. The resulting HIMF has a shallow faint end slope (alpha ~ 1.2), and is consistent with earlier estimates computed for the population of optically selected gas rich galaxies. This implies that there is not a large population of gas rich low luminosity or low surface brightness galaxies that has gone unnoticed by optical surveys. The cosmological mass density of HI at the present time determined from the survey, Omega_HI = (2.0 +/- 0.5) x 10^-4, is in good agreement with earlier estimates. We determine lower limits to the average column densities <N_HI> of the galaxies detected in the survey and find that none of the galaxies have <N_HI> below 10^19.7 cm^-2, although there are no observational selection criteria against finding lower density systems.

The Fisher-Tully `10 Mpc Catalog of Late-Type Galaxies' (FT, 1981) is remarkably complete. Despite the considerable effort that has been spent searching for and cataloging low surface brightness (LSB) galaxies, almost no new HI-rich galaxies have been added to the volume to which the FT observations were both sensitive (d < 10 (M_{HI}/10^{8.2} Msolar)^{5/12} Mpc) and intended to be complete (d < 10 Mpc). It has not yet been demonstrated by the new surveys that HI-rich LSB galaxies are not properly represented in this catalog of nearby galaxies. Although new optical surveys are discovering heretofore uncatalogued galaxies at greater distances than the depth of the FT Catalog, they do not alter the completeness of the FT Catalog for HI-rich LSB objects. Thus, unless the galaxy population of the local volume is atypical, there is at present no evidence for a significant, new population of gas-rich galaxies that has escaped recognition until now.

The space density of low surface brightness and tiny gas-rich dwarf galaxies are estimated for two recent catalogs: The Arecibo Survey of Northern Dwarf and Low Surface Brightness Galaxies (Schneider, Thuan, Magri & Wadiak 1990) and The Catalog of Low Surface Brightness Galaxy, List II (Schombert, Bothun, Schneider & McGaugh 1992). The goals are (1) to evaluate the additions to the completeness of the Fisher and Tully (1981) 10 Mpc Sample and (2) to estimate whether the density of galaxies contained in the new catalogs adds a significant amount of neutral gas mass to the the inventory of HI already identified in the nearby, present-epoch universe. Although tiny dwarf galaxies (M_HI < ~10^7 solar masses) may be the most abundant type of extragalactic stellar system in the nearby Universe, if the new catalogs are representative, the LSB and dwarf populations they contain make only a small addition (<10%) to the total HI content of the local Universe and probably constitute even smaller fractions of its luminous and dynamical mass.

Radio spectroscopy at 323 MHz using the Arecibo Telescope and the VLA has produced a tentative detection of the 21cm line of neutral hydrogen at Z_abs=3.38716+-0.00007 in absorption against the radio continuum of the QSO MG0201+113 (Z_em = 3.61). This redshift roughly agrees with one determined by the Westerbork Synthesis Radio Telescope and reported earlier; however the observations at different telescopes produce different results for width and optical depth of the line. If the detection holds, it provides evidence for a high column density of neutral gas that is confined to a dynamically cold layer with velocity dispersion ~10 km/s. Although the interpretation is uncertain due to a lack of detailed knowledge of the extended radio structure of the background quasar and the relative quantities of neutral gas in the cold and turbulent components, the observations specify high spin temperatures, T_s > 1000 K for both this 21cm line absorbing cloud and the turbulent component, which together have N_HI ~ 10^21.4 cm^-2. Results of optical spectroscopy require the additional presence of metal-enriched clouds of still broader velocity dispersion than the 21cm line, creating a picture which is consistent with this system being a young disk galaxy that is sufficiently evolved by z=3.4 to have collapsed to a flattened system and produced a population of stars that have aged to pollute a turbulent halo. The observations constrain the neutral gas mass of a possible 21cm line emitter associated with the intervening absorber to be less than ~(Delta V/200 km/s) 10^13 Msolar for velocity widths Delta V (FWHM) in the range 200 to 1200 km/s (H_o=100 km/s/Mpc, Omega_o=1).

Analyses of QSO absorption lines are showing that HI content has evolved over the redshift range z=5 to z= 0. The 21cm line measurements of the z=0 HI content avoid several biases inherent in the absorption line technique, such as the influence of evolving dust content in the absorbers, and will produce a reliable measure to anchor theories of galaxy evolution. Examples of important questions to be addressed by local HI surveys are: (1) is there a significant population of gas-rich galaxies or intergalactic clouds that are missing from the census of optically selected galaxies? (2) is there an adequate reservoir of neutral gas to substantially prolong star formation at its present rate? and (3) are there massive objects of such low HI column density that they can have escaped detection in the ``unbiased'' HI surveys that have been conducted so far?

Driftscan methods are highly efficient, stable techniques for conducting extragalactic surveys in the 21cm line of neutral hydrogen. Holding the telescope still while the beam scans the sky at the sidereal rate produces exceptionally stable spectral baselines, increased stability for RFI signals, and excellent diagnostic information about system performance. Data can be processed naturally and efficiently by grouping long sequences of spectra into an image format, thereby allowing thousands of individual spectra to be calibrated, inspected and manipulated as a single data structure with standard tools that already exist in astronomical software. The behavior of spectral standing waves (multi-path effects) can be appraised and excised in this environment, making observations possible while the Sun is up. The method is illustrated with survey data from Arecibo and Nancay.

Results from the Arecibo HI Strip Survey, an unbiased extragalactic HI survey, combined with optical and 21cm follow-up observations, determine the HI Mass Function and the cosmological mass density of HI at the present epoch. Both are consistent with earlier estimates, computed for the population of optically selected galaxies. This consistency occurs because, although the distribution of optical central surface brightnesses among galaxies is flat, we fail to find a population of galaxies with central surface brightnesses fainter than 24 B-mag/arcsec^2, even though there is no observational selection against them.