Hostname: page-component-cd9895bd7-gxg78 Total loading time: 0 Render date: 2024-12-27T06:12:20.511Z Has data issue: false hasContentIssue false

The Phase II Murchison Widefield Array: Design overview

Published online by Cambridge University Press:  23 November 2018

Randall B. Wayth*
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
Steven J. Tingay
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
Cathryn M. Trott
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia
David Emrich
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
Melanie Johnston-Hollitt
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
Ben McKinley
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia
B. M. Gaensler
Affiliation:
ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St George St, Toronto, ON M5S 3H4, Canada
A. P. Beardsley
Affiliation:
School of Earth and Space Exploration, Arizona State University, Tempe, AZ 85287, USA
T. Booler
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
B. Crosse
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
T. M. O. Franzen
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
L. Horsley
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
D. L. Kaplan
Affiliation:
Department of Physics, University of Wisconsin–Milwaukee, Milwaukee, WI 53201, USA
D. Kenney
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
M. F. Morales
Affiliation:
Department of Physics, University of Washington, Seattle, WA 98195, USA
D. Pallot
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), University of Western Australia, Crawley, WA 6009, Australia
G. Sleap
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
K. Steele
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
M. Walker
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
A. Williams
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), Curtin University, Bentley 6845, Australia
C. Wu
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), University of Western Australia, Crawley, WA 6009, Australia
Iver. H. Cairns
Affiliation:
Sydney Institute for Astronomy (SIfA), School of Physics, The University of Sydney, NSW 2006, Australia
M. D. Filipovic
Affiliation:
Western Sydney University, Locked Bag 1797, Penrith South DC, NSW, 1797, Australia
S. Johnston
Affiliation:
CSIRO Astronomy & Space Science, Australia Telescope National Facility, P.O. Box 76, Epping, NSW 1710, Australia
T. Murphy
Affiliation:
Sydney Institute for Astronomy (SIfA), School of Physics, The University of Sydney, NSW 2006, Australia
P. Quinn
Affiliation:
International Centre for Radio Astronomy Research (ICRAR), University of Western Australia, Crawley, WA 6009, Australia
L. Staveley-Smith
Affiliation:
ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia International Centre for Radio Astronomy Research (ICRAR), University of Western Australia, Crawley, WA 6009, Australia
R. Webster
Affiliation:
ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia School of Physics, The University of Melbourne, Parkville, VIC 3010, Australia
J. S. B. Wyithe
Affiliation:
ARC Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), Curtin University, Bentley, 6102, Australia School of Physics, The University of Melbourne, Parkville, VIC 3010, Australia
*
Author for correspondence: Randall B.Wayth, Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

We describe the motivation and design details of the ‘Phase II’ upgrade of the Murchison Widefield Array radio telescope. The expansion doubles to 256 the number of antenna tiles deployed in the array. The new antenna tiles enhance the capabilities of the Murchison Widefield Array in several key science areas. Seventy-two of the new tiles are deployed in a regular configuration near the existing array core. These new tiles enhance the surface brightness sensitivity of the array and will improve the ability of the Murchison Widefield Array to estimate the slope of the Epoch of Reionisation power spectrum by a factor of ∼3.5. The remaining 56 tiles are deployed on long baselines, doubling the maximum baseline of the array and improving the array u, v coverage. The improved imaging capabilities will provide an order of magnitude improvement in the noise floor of Murchison Widefield Array continuum images. The upgrade retains all of the features that have underpinned the Murchison Widefield Array’s success (large field of view, snapshot image quality, and pointing agility) and boosts the scientific potential with enhanced imaging capabilities and by enabling new calibration strategies.

Type
Research Article
Copyright
Copyright © Astronomical Society of Australia 2018 

1. Introduction

Radio astronomy at low radio frequencies (RFs) (<500 MHz) has undergone a major renaissance over the last decade, with the establishment and operation of a new generation of large-scale observational facilities. The renaissance has been enabled by advances in signal processing and computing and has been motivated by challenging problems in astrophysics and cosmology.

Perhaps the most visible science challenge at low RFs is the detection of the redshifted H I signal from the early universe, the so-called Cosmic Dawn and Epoch of Reionisation (EoR) (e.g. Furlanetto, Oh, & Briggs Reference Furlanetto, Oh and Briggs2006; Koopmans et al. Reference Koopmans2015). A number of responses to this challenge have emerged, in the form of new radio telescopes and the refurbishment of existing radio telescopes. For example, the Giant Metrewave Radio Telescope (GMRT; Gupta Reference Gupta2014) has undergone a number of upgrades. New telescopes such as the Low Frequency ARray (LOFAR; van Haarlem et al. Reference van Haarlem2013) and the Murchison Widefield Array (MWA; Tingay et al. Reference Tingay2013a), while strongly motivated by the EoR power spectrum experiment, are general-purpose instruments capable of a wide range of investigations. Other new instruments such as the Precision Array for Probing the Epoch of Reionisation (PAPER; Parsons et al. Reference Parsons2010), the Hydrogen Epoch of Reionisation Array (HERA; DeBoer et al. Reference DeBoer2017), and the 21 Centimeter Array (21CMA; Zheng et al. Reference Zheng, Wu, Johnston-Hollitt, Gu and Xu2016) are entirely motivated by, and designed for, the EoR power spectrum experiment. Collectively, these telescopes are making steady progress towards an EoR power spectrum detection over a range of redshifts (e.g. Paciga et al. Reference Paciga2011; Jacobs et al. Reference Jacobs2015; Ali et al. Reference Ali2015; Beardsley et al. Reference Beardsley2016; Patil et al. Reference Patil2017).

The nature of these instruments also provides us with a number of new practical challenges. The large data challenge is one. Other challenges include imaging and calibration over wide fields of view, wide fractional bandwidths, and taking account of complex propagation effects (e.g. ionospheric and interplanetary scintillation).

For those instruments designed to be more general purpose than just an ‘EoR machine’, these challenges are more than compensated for by the vast new areas of parameter space opened up by the unique characteristics of these instruments: (1) large fields of view on the sky; (2) sensitive, high-fidelity, and high-resolution imaging at low frequencies; (3) flexibility in choice of output data products; and (4) large fractional bandwidths at low frequencies. Already, LOFAR and the MWA have proven the value of this new parameter space, with important new fields of research strongly emerging well before all the technical challenges have been completely met.

The future Square Kilometre Array (SKA) will build from the current generation of low-frequency instruments, in particular the SKA precursors and pathfinders MWA, HERA, and LOFAR. In the vein of the MWA and LOFAR, the low-frequency SKA (SKA Low) will have a wide range of science goals and will open up an even wider parameter space for discovery.

The MWA is located at the same remote and radio quiet site as SKA Low will be, the Murchison Radio-astronomy Observatory (MRO). The design, construction, commissioning, and operation of the MWA have generated deep experience that is directly relevant to SKA Low. In addition, the MWA has been directly utilised as a development platform to prototype and test various subsystems for SKA Low (e.g. Sutinjo et al. Reference Sutinjo2015; Benthem et al. Reference Benthem2017; Wayth et al. Reference Wayth2017).

Low-frequency radio telescopes such as the MWA are naturally suited to upgrades, since the front-end infrastructure and antenna systems are relatively simple and the back-end signal processing and data processing are driven by the ever increasing capabilities of the commodity electronics market. The design philosophy of the MWA, incorporating commercial off-the-shelf components where possible (e.g. Ord et al. Reference Ord2015; Tremblay et al. Reference Tremblay2015), has led to a flexible and intrinsically extensible instrument that is responsive to the dynamic scientific environment created when a new and large parameter space is entered.

Here we describe the next steps in the evolution of the MWA, with a doubling of the number of antennas in the array, a doubling of the maximum baseline in the array, and an exploration of calibration techniques based on alternative array architectures. These front-end upgrades address a number of different science areas, in keeping with a general-purpose and flexible instrument, and pave the way for future upgrades that will be more focused on the signal and data processing in the back-end. The material described in this paper demonstrates that once the investment in the basic infrastructure has been made, substantial upgrades to the scientific capabilities of instruments like the MWA can be made in a highly efficient and cost-effective manner. We expect similar benefits will apply to the SKA Low telescope.

In this paper, we refer to the original system that was deployed and commissioned in 2012/2013 as the ‘Phase I’ MWA, which is described by Tingay et al. (Reference Tingay2013a).

2. Scientific motivation for the upgrade

The Phase I MWA was motivated by a vast array of science goals spanning initial observations of the EoR power spectrum, Galactic and Extragalactic continuum and polarimetry, solar and ionospheric science, and detection of transient sources (Bowman et al. Reference Bowman2013). In addition to covering the initial science goals, the Phase I MWA also produced a number of unforeseen results including the first detection of plasma tubes in the ionosphere (Loi et al. Reference Loi2015a, Reference Loi2015b), the potential to utilise the MWA for space debris detection and tracking (Tingay et al. Reference Tingay2013b; Zhang et al. Reference Zhang2018), tentative detections of new molecules (Tremblay et al. Reference Tremblay, Hurley-Walker, Cunningham, Jones, Hancock, Wayth and Jordan2017), and use of interplanetary scintillation to identify and measure sub-arcsecond compact components in low-resolution, low-frequency radio surveys without the need for long baseline interferometry or ionospheric calibration (Kaplan et al. Reference Kaplan2015a; Morgan et al. Reference Morgan2018; Chhetri et al. Reference Chhetri, Morgan, Ekers, Macquart, Sadler, Giroletti, Callingham and Tingay2018). The wide field of view and pointing agility of the MWA have also proven invaluable for both targeted and archival searches (from serendipitous observations) for radio emission associated with astronomical events such as gamma-ray bursts (Kaplan et al. Reference Kaplan2015b), neutrino detections (Croft et al. Reference Croft2016), gravitational wave events (Abbott et al. Reference Abbott2016, Reference Abbott2017), and interstellar visitors (Tingay et al. Reference Tingay2018).

Given the successes of the Phase I MWA, when considering an upgrade to the facility, scientists generally desire improved capability in all key performance metrics of the telescope (sensitivity, angular resolution, frequency coverage, frequency resolution, time resolution, number of beams, etc). Budgetary, timeline, and practical constraints for the Phase II upgrade meant that the focus was put on expanding the physical array, with future upgrades to consider receivers and downstream digital systems.

The design of the Phase II array thus focuses on adding antenna tiles to improve capability for MWA Key Science and to enhance capabilities in science areas where the MWA is a particularly well-suited instrument. Broadly, this comes down to (1) improving the sensitivity (and calibratability) for the EoR power spectrum measurement, and (2) improving the imaging properties of the array through improved angular resolution, u, v coverage, and reduced confusion.

2.1. Epoch of Reionisation

When the Phase I MWA was deployed, it was one of the first generation of telescopes aiming to detect the EoR power spectrum. It was designed as a general-purpose array with emphasis on imaging capability and brightness sensitivity on the angular scales relevant to the EoR. In the years following the initial design and deployment, it has become increasingly apparent that the dynamic range requirements (e.g. accuracy of calibration, accuracy of sky model, etc.) of the EoR detection experiment are just as challenging as the sensitivity requirements (Jacobs et al. Reference Jacobs2016; Trott et al. Reference Trott2016; Beardsley et al. Reference Beardsley2016; Barry et al. Reference Barry, Hazelton, Sullivan, Morales and Pober2016; Patil et al. Reference Patil2016, Reference Patil2017).

In contrast to the MWA, the PAPER array (Parsons et al. Reference Parsons2010) pursued a more targeted EoR power spectrum experiment, using a regular grid of simple dipole antennas. The regular antenna arrangement is especially useful for the EoR power spectrum, because many identical (redundant) baselines formed by a regular grid boost the sensitivity for a power spectrum measurement (Parsons et al. Reference Parsons, Pober, McQuinn, Jacobs and Aguirre2012). A regular array configuration also enables one to exploit the benefits of redundant calibration. This removes some dependence on having an accurate sky model to calibrate the antenna gains. The under-construction HERA telescope (DeBoer et al. Reference DeBoer2017) also uses a regular configuration.

Motivated by this experience with existing EoR telescopes and experiments, 72 of the new antenna tiles were deployed in a regular hexagonal configuration. Section 3 provides details of the deployment and Section 4 discusses the performance implications of the upgrade.

2.2. Image-based science

Aside from raw sensitivity, the capability of a radio interferometer is usually constrained by its angular resolution and by the spatial scales that are measured given the antenna layout. The Phase I MWA layout is a hybrid configuration having baselines ranging in length from approximately 7–2 800 m, with many tiles in the core region. The synthesised beam size of the Phase I array is approximately 2 arcmin at 150 MHz. Improving the angular resolution of the MWA for Phase II was identified for its potential to impact a broad range of science areas. In addition to the ability to better resolve objects and structures within sources, improved angular resolution directly impacts the classical and sidelobe confusion in continuum images.

For Phase II, 56 new antenna tiles were deployed beyond the circumference of the Phase I array. The new tiles double the maximum baseline of the array and generate more uniform u, v coverage. The new tiles require additional hardware to transmit signals over long distances as described in 3.2. The performance implications of the new tiles are discussed in 4.3.

3. The Phase II upgrade

The existing MWA consists of 16 in-field receiver units, each of which services 8 antenna tiles (see Tingay et al. Reference Tingay2013a). The Phase II upgrade adds an additional 128 tiles, taking the total number of deployed tiles to 256. The 128 new tiles comprise:

  • 72 tiles arranged in two regular hexagonal configurations near the existing MWA core (Figure 1); and

  • 56 long baseline tiles placed beyond the limit of the existing array (Figure 2), taking the maximum baseline for the Phase II array to approximately 5.3 km.

The existing receivers (Prabu et al. Reference Prabu2015) and correlator (Ord et al. Reference Ord2015) were retained for this expansion, which limits the total number of tiles that can operate at any time to 128.Footnote a As such, the Phase II MWA will be periodically reconfigured between compact and extended configurations, each of which contains 128 operating tiles. Reconfiguration involves physically disconnecting tiles from some receivers, relocating the receivers, then connecting new tiles.Footnote b

Figure 1. The compact configuration including the 72 new tiles arranged as two regular hexagonal arrays (filled squares). Note the size of the squares is not to scale.

Figure 2. The extended configuration including the 56 new long baseline tiles (filled squares). The inset rectangle shows the area bounded by Figure 1.

The overall goal was to have the maximum baseline of the Phase II array be as large as possible while retaining the excellent snapshot u, v coverage that has underpinned much of the MWA’s success. Placement of the long baseline tiles was constrained by a number of practical and logistical considerations.

A feature of the original MWA deployment that has proven invaluable for the Phase II expansion is the excess power and communications infrastructure that was initially installed. The MWA’s in-field receivers are placed on concrete pads next to service points for buried power and communication infrastructure. During the infrastructure work for the initial MWA deployment, it was decided to install extra buried infrastructure beyond what was needed for Phase I. The reason was that the additional cost of service infrastructure was small compared to the overall infrastructure costs, and this excess capacity would be useful for future expansion. During Phase II construction, extra concrete pads were placed near the central region of the MWA, thereby creating sufficient space to utilise the extra buried capacity to host all 16 receivers near the core, which is required for the new compact configuration.

When choosing which existing tiles were to be used in either of the new configurations, we opted to include or exclude all eight tiles on a receiver (hence, we retain or change all tiles connected to a receiver) for practical reasons. Hence, the two new configurations consist of a subset of the receivers from the Phase I array (and all their associated tiles), plus new tiles, which are connected to relocated receivers.

3.1. Compact hexagonal configuration antenna tiles

The compact configuration consists of 7 Phase I receivers and their corresponding 56 tiles, which comprise the bulk of the existing MWA core tiles, plus the 72 new hexagonal configuration tiles. The nine receivers required for the new tiles are relocated from elsewhere in the array, mostly from the existing Phase I outermost tiles. The layout of the compact configuration is shown in Figure 1.

The vast majority of the baselines in the compact configuration are shorter than 200 m, consistent with requirements for an EoR power spectrum detection being the primary science driver for this configuration. Clearly the compact configuration is a hybrid between the regular hexagonal configuration tiles and the pseudo-random core of the existing array, with the regular hexagonal configurations providing many redundant baselines.

The hexagonal configuration of new tiles was chosen over other possibilities (such as a rectangular grid which would have better theoretical sensitivity for the power spectrum) for several reasons. Firstly, the hexagonal configuration has less pronounced grating features compared to a rectangular grid; the configuration matches the HERA configuration, hence will allow comparison of results for the same modes; the size of the hexagons is a good match to the size of the existing Phase I core region; finally, the hexagonal arrangement avoids potential issues with a rectangular grid’s grating features beating against the primary beam of the MWA’s square tiles.

3.2. Long baseline tiles

The extended configuration consists of 9 existing receivers (essentially the Phase I array minus the core) and their corresponding 72 tiles, plus the 56 new long baseline tiles. The layout of the extended configuration is shown in Figure 2. The seven receivers required to service the new long baseline tiles are relocated from the core region.

All existing MWA tiles (and the new hex configuration tiles) are serviced by twin coaxial cables, which carry power and communications to the MWA beamformers, and the RF signal back from the tiles. The new long baseline tiles are substantially further from the nearest receiver pad that can be feasibly reached via coaxial cable, hence the new long baseline tiles use self-contained solar power units and RF-over-fibre (RFoF) technology to transmit the RF signal back to receivers.

Significant drivers for the cost and schedule of the expansion project were environmental and local heritage clearances, which are required if the land is disturbed by clearing and/or digging. As such, the additional infrastructure work to support the long baselines was limited to extending the existing access tracks; no new buried cables were installed. All new equipments for the long baseline tiles (including the tiles themselves, power supplies, and fibre cables) were placed on the surface.

The locations of the long baseline tiles were determined by initially selecting target locations within existing boundaries of the site, such that the baselines were as long as possible and the resulting u, v coverage was sensible. Each target location was visited and inspected such that a suitably clear, flat site near the target location could be chosen. For all but two of the long baseline tiles, a suitable site was found within approximately 50 m of the original target location, and the remaining two were moved by several hundred metres to avoid locally sensitive areas. This process of on-the-ground micro-siting was critical in meeting the cost, schedule, and compliance requirements of the upgrade project.

3.3. Power and signal transport for long baseline tiles

All Phase I tiles (and the new Phase II hex configuration tiles) have beamformers that connect to MWA receivers via two coaxial cables that range in length up to a few hundred metres. The cables supply power to beamformers, carry RF signals from the beamformers, and provide two-way communications via a simple serial protocol. This system is supported by custom-designed ‘data over coax’ (DoC) circuit boards, one of which is in the beamformer (the beamformer-side DoC), and the other is in the receiver (the receiver-side DoC).

The long baseline tiles use standard MWA beamformers. In order to operate these beamformers, supporting infrastructure that mimics the functionality provided by the receiver was required. A bespoke system was designed to do this, as illustrated in Figure 3.

Figure 3. Block diagram of long baseline tile power, control, and signal transport system.

The long baseline tiles are connected to the rest of the MWA via a 6-core optical fibre cable. This cable provides standard two-way ethernet communications (using two fibres) and carries the RF signals from the tiles (on another two fibres, hence there are two spares). Power and communications with the tiles are provided via the beamformer interface (BFIF) unit, which includes a Raspberry Pi single-board computer (SBC). The SBC controls a receiver-side DoC board, which in turn connects to the beamformer via short lengths of coaxial cable. The dual-polarisation RF output from the DoC board is fed into two RF-over-fibre drivers (ASTRON model 350).

Power is provided by two 315 W, 24 V solar panels (Suntech STP315/24/VEM), which are mounted on a fixed frame nearby the tile. Four 150 Ah, 12 V lead-acid batteries provide energy storage, arranged as two sets of 24-V batteries. The power subsystem is controlled by a solar power interface unit, which provides a low-level monitor and control serial interface that is connected to the SBC.

A key requirement for the power subsystem is to be radio quiet, so a typical off-the-shelf solution (which commonly uses rapid switching) was not feasible. Instead, a custom low-power low-frequency switching supply was designed to trickle charge the batteries while also powering the tile. The trickle charge regulators can be bypassed under software control, allowing one or both panels to be directly connected to the batteries to provide a more rapid charge, depending on battery charge state and solar panel output. Software on the SBC monitors and controls this subsystem closely to ensure optimum battery lifetime and additionally will shut down most of the loads if required to protect batteries from over-discharge. Finally, a carefully selected and tested off-the-shelf medium-frequency switch-mode power supply is used internally to convert the 24 V from the panel/battery subsystem to the 48 V required for the beamformer.

Given the proximity of this solar power subsystem to its nearby antenna tile, we undertook careful electromagnetic compatibility (EMC) design to ensure that any radio emission generated by the charging electronics would not interfere with the astronomical signals. We therefore designed the enclosure to be very similar to the beamformer enclosure, which is known to provide good shielding from radiated emissions, and built a custom EMC filter to prevent any emissions from leaving the enclosure on power cables. Figure 4 shows an image of a Phase II long baseline tile.

Figure 4. A Phase II long baseline tile. The batteries and BFIF systems are located under the solar panels.

Modifications were also required to the receivers connected to long baseline tiles. An 8-port network switch with optical ports was fitted in each modified receiver to distribute the control and telemetry signals for the tiles. Also, the original 8 receiver-side DoC boards were removed and replaced with a 16-channel RFoF receiver board, which converts the 16 optical RF signals back to electrical signals. These were then connected to the existing Analogue Signal Conditioner modules. The remaining components of the receivers were unchanged.

3.4. Correlator and archive

The digital systems of the MWA remain the same for this upgrade, although the new longer baselines do have an effect on the operation of the correlator and onsite archive. Phase I observations were typically performed with time and frequency averaging in the correlator set to keep the output data rate within the capabilities of the archive. For example, observations requiring 0.5-s time resolution were required to average the native 10-kHz frequency resolution down to 40 kHz. Likewise, observations requiring the full 10-kHz frequency resolution were required to time-average the data to 2 s.

The Phase I correlator (Ord et al. Reference Ord2015) does not fringe track. While this greatly simplifies the requirements of the correlator itself, the natural fringe frequency of the longest baselines sets an upper limit to the time averaging that can be performed in the correlator. Similarly, cable length differences that are uncompensated before the correlator set an upper limit to the frequency averaging that can be performed in the correlator. The net effect of this is that for the extended array, the correlator must be run with higher time and frequency resolution as compared to Phase I. Typically this is 0.5-s time resolution and 10-kHz frequency resolution. Compared to Phase I, this increases the maximum data rate generated by the correlator by a factor of 4.

To handle the larger data rate, two enhancements were made to the correlator and archive: expansion of the onsite archive system and enabling in-situ compression of visibility data in the correlator. The onsite archive (Wu et al. Reference Wu, Wicenec, Pallot and Checcucci2013) was upgraded by replacing the two existing data servers with six Dell PowerEdge R730xd servers (each with approximately 66 TB of usable storage) and upgrading the switching network between the correlator machines and archive servers.

The in-situ compression of data takes advantage of the fact that the raw correlation coefficients generated by the correlator are integers, since the data going into the correlator are 4-bit integers (Ord et al. Reference Ord2015). Although the xGPU library performs cross-multiplication and accumulation as floating-point numbers, the resulting correlation coefficients are simply integers with very small floating-point accumulation/rounding errors. Thus, instead of storing the coefficients as floating-point numbers, they are rounded and stored as 32-bit integers using the in-built Rice compression functionality provided by the CFITSIO library. This transformation library.Footnote c This transformation is lossless. For typical observations, the correlation coefficients are less than 10 000 or so, and hence the 32-bit integers are highly compressible, whereas the same number stored in floating-point format would not be. For typical observations, this scheme reduces the data rate out of the correlator by approximately 2.5:1 compared to uncompressed data.

The downstream post-processing tool Cotter (Offringa et al. Reference Offringa2015) (which applies delay corrections and creates a phase centre for the data, in addition to flagging) provides functionality for time and/or frequency averaging such that end users may not have to deal with the larger data volumes when the data are being converted from raw format (stored in FITS files) to a standard radio astronomy format. In addition, the compression feature specified in the FITS standard is handled transparently by all tools; hence, no modifications to downstream tools were required.

4. Capability of the upgraded MWA

4.1. Calibration using redundant baselines

Redundant calibration, whereby the observation of the sky using many identical baselines recovers the antenna gains without requiring a sky model (Wieringa Reference Wieringa1992; Liu et al. Reference Liu, Tegmark, Morrison, Lutomirski and Zaldarriaga2010; Zheng et al. Reference Zheng2014; Dillon & Parsons Reference Dillon and Parsons2016), was used by the PAPER telescope and is central to the proposed design for HERA. Despite its promise of more accurate calibration by removing the dependence on a sky model, redundant calibration has its own set of challenges for real instruments with differing antenna responses and suffers from degeneracies in some situations. The MWA Phase II compact configuration addresses these challenges by combining information from redundant and non-redundant baselines, allowing a more holistic approach to calibration while also providing information about the variations between different antenna tiles (Li et al. Reference Li2018).

4.2. Power spectrum sensitivity

Figure 5 shows the theoretical noise levels for a 1 000-h observation with the original MWA and the upgraded MWA, with the extra sensitivity at small k-modes, for a typical 21-cm signal model at 150 MHz (using 21cmFAST; Mesinger, Furlanetto, & Cen Reference Mesinger, Furlanetto and Cen2011). Over the range of scales of relevance to the EoR power spectrum at z∼8 (u = [2, 200]), the raw sensitivity increase is a factor of 4, as can be observed in the figure.

Figure 5. Typical EoR power spectrum model at 150 MHz with associated noise levels available to the Phase I and Phase II arrays with 1 000-h observation. ‘Phase II 256’ represents the result from a future MWA upgrade where all 256 tiles are used simultaneously.

To provide a realistic assessment of the performance improvement for Phase II, a realistic foreground model is included in the analysis. For both Phase I and II, compact sources are assumed to be subtracted to five times the classical confusion level, with the lower confusion level for Phase II accounted for. Phase II therefore offers both improved sensitivity on large scales and lower foreground contamination. (Figure 5 includes a future ‘Phase II 256’ array where all 256 antenna tiles are correlated.) Note that Phase II arrays with both the compact configuration of 128 tiles and a full extended 256-tile array yield comparable performance on the large angular scales of importance for EoR science. Finally, we perform a Fisher Analysis of the theoretical signal-to-noise ratio attainable on the slope and offset of the power spectrum (in the presence of realistic foregrounds), yielding an expected improvement in parameter estimation performance by a factor of 3.5 for Phase II compared with Phase I.

4.3. Imaging

Since both the compact and extended arrays contain 128 tiles, the theoretical radiometric sensitivity of MWA Phase II is the same as the Phase I array. In practice, however, the spatial resolution and u, v coverage of the Phase II arrays mean the two configurations perform very differently depending on the application. The snapshot monochromatic u, v coverage of the compact and extended configurations are shown in Figure 6.

Figure 6. Snapshot monochromatic u, v coverage of the Phase II configurations at 154 MHz. The left panel shows the compact array (zenith pointed) and the right panel the extended array (zenith pointed).

The extended configuration doubles the maximum baseline length compared to the Phase I array and also does not contain the many short baselines in the core. This means the u, v plane is more uniformly filled compared to the Phase I array, with an associated improvement in the synthesised beam. Figure 7 shows example 1D cuts through the MWA’s synthesised beam for the compact and extended configurations and the Phase I array for two different visibility weighting schemes, roughly corresponding to uniform weights (robust −1) and natural weights (robust 1). The example uses only 1 MHz of bandwidth, so corresponds closely to the ideal monochromatic synthesised beam.

Figure 7. Example cuts of the synthesised beam (zenith pointed) at 154 MHz for the Phase I and Phase II arrays with different robust weights. In this example, 1 MHz of bandwidth is used.

For continuum imaging compact sources with the Phase I array, it was usually required to use visibility weighting that was closer to uniform weighting (e.g. Hurley-Walker et al. Reference Hurley-Walker2017), with the associated loss in sensitivity. This is because the naturally weighted Phase I synthesised beam has a large ‘halo’ around the central peak due to the compact core, as can be seen in Figure 7. Conversely, the Phase II extended array has an almost identical beam for both visibility weighting schemes and very low sidelobes away from the phase centre.

Previous deep 154-MHz images made with the MWA Phase I (Offringa et al. Reference Offringa2016; Franzen et al. Reference Franzen2016) were limited by confusion, with sidelobe confusion contributing to the overall image noise. The improved angular resolution of the extended array will reduce the classical confusion in continuum images, with the precise factor depending on the slope of the differential source count distribution around 1 mJy (Condon Reference Condon1974). Using data from the Phase I MWA and other low-frequency telescopes, Franzen et al. (Reference Franzen2016) estimated the classical confusion for Phase I MWA at 154 MHz to be 1.7 mJy and the expected reduction in classical confusion to be a factor of 5–10 for Phase II. The Phase II extended array’s improved synthesised beam will also reduce sidelobe confusion. The combination of improved sensitivity, improved classical confusion, and reduced sidelobe confusion means an order of magnitude improvement in the noise floor of continuum images from the extended array is expected.

An example of the improvement in resolution achieved between the Phase I and Phase II MWA is shown in Figure 8 which shows the 185-MHz Phase I and Phase II images of Fornax A. These images were made with approximately 6 min of data for each, with the same frequency, pointing, and Local Sidereal Time for each data set. The same procedure to calibrate and image was used for both data sets using robust 0 visibility weighting.

Figure 8. 185 MHz images of Fornax A taken with the Phase I (left) and Phase II (right) MWA. The ellipse in the lower left of each image shows the synthesised beam sizes, which are 2.6×2.3 and 1.0×0.89 arcmin.

By design, the compact configuration has very many identical baselines, which are favourable for the EoR power spectrum experiment but are sub-optimal for imaging applications. The synthesised beam contains regularly spaced grating-like features separated by approximately 9° at 154 MHz, which are due to the 14-m separation between tiles in the hexagons. Using uniform-like weighting can suppress these grating sidelobes at the expense of approximately 50% loss in theoretical sensitivity. Imaging with the compact array will require careful selection of the visibility weighting scheme depending on the application.

4.4. The link to SKA Low

It is worth noting that the trade-offs associated with sensitivity and synthesised beam quality of the MWA (especially the Phase I MWA) are directly relevant to SKA Low. Like Phase I MWA, the proposed layout of SKA1 Low (Dewdney & Braun Reference Dewdney and Braun2016) contains a densely packed core of approximately 45% of the antennas with the remaining antennas spread out for angular resolution and u, v coverage. The naturally weighted SKA Low synthesised beam will thus also have a significant halo around the central peak, making deconvolution difficult. High-resolution continuum imaging with SKA1 Low will thus likely require significant down-weighting of the shortest baselines with associated reduction in sensitivity.

A future upgrade to the MWA’s correlator will allow all 256 antennas to be processed simultaneously. This upgrade will make the MWA’s signal processing requirements comparable to SKA Low. As a precursor to SKA Low, the MWA is an ideal system for proving calibration and imaging pipelines. The ratio of the maximum baseline to station size (hence number of pixels in images) is almost identical, as is the fraction of antennas in the core region. Expertise and lessons from MWA data processing are thus very relevant for SKA Low.

5. Conclusion

This paper describes the Phase II upgrade of the MWA—the first major expansion of the radio telescope since it was constructed in 2012. The MWA continues to be a general-purpose facility underpinned by its key science programmes (Bowman et al. Reference Bowman2013). In the first 5 years of telescope operations, the capabilities of the MWA have been demonstrated across a wide variety of science domains, producing over 100 refereed publications that have generated over 3 000 citations.

The Phase II upgrade enhances the capabilities of the MWA in angular resolution, brightness sensitivity, and sensitivity to the EoR power spectrum while retaining (and even improving) some of the unique features that have made the originally deployed array so successful: the very large field of view combined with the excellent snapshot imaging capability. The extended configuration, which includes new tiles that double the maximum baseline length, doubles the angular resolution of the array and improves the snapshot u, v coverage. The improved image resolution and synthesised beam will translate to an order of magnitude improvement in continuum image sensitivity due to the reduction in classical and sidelobe confusion.

The compact configuration, which includes the new regularly spaced hexagonal arrays, will improve the sensitivity to the EoR power spectrum by a factor of 3.5. The compact configuration also allows improved and novel array calibration techniques that make use of the hybrid nature of the array layout, which contains redundant and non-redundant baselines.

The improvement to the MWA does not stop with the Phase II upgrade; work is ongoing to upgrade several of the digital systems of the MWA. These projects include an upgraded correlator that will correlate all 256 tiles with increased bandwidth, and a new digital receiver that removes the limitations of the existing system. Work is also ongoing to incorporate commensal data processing capability for the Breakthrough Listen project.

Finally, future upgrades of the MWA may include (but are not limited to) retrofitting RFoF transmitters to core tiles to improve the bandpass spectral smoothness; adding more tiles to further extend the maximum baseline; changing the low-frequency cutoff of the analogue components to 50 MHz (since the dipoles have been shown to perform down to 50 MHz with modified Low-Noise Amplifiers in the Engineering Development Array (EDA, Wayth et al. Reference Wayth2017); and fully integrating the signal paths from the EDA, SKA Low Aperture Array Verification System (Benthem et al. Reference Benthem2017) and other SKA Low prototype systems.

Acknowledgements

The MWA Phase II upgrade project was supported by Australian Research Council LIEF grant LE160100031 and the Dunlap Institute for Astronomy and Astrophysics at the University of Toronto. Parts of this research were supported by the Australian Research Council Centre of Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project number CE170100013. BM is supported by an ARC Discovery Early-Career Researcher Award, under project number DE160100849. Parts of this research were supported by the Australian Research Council Centre of Excellence for All-sky Astrophysics (CAASTRO), through project number CE110001020. This scientific work makes use of the Murchison Radio-astronomy Observatory, operated by CSIRO. We acknowledge the Wajarri Yamatji people as the traditional owners of the Observatory site. Support for the operation of the MWA is provided by the Australian Government (NCRIS), under a contract to Curtin University administered by Astronomy Australia Limited. We acknowledge the Pawsey Supercomputing Centre which is supported by the Western Australian and Australian Governments.

Footnotes

a No new receiver hardware is available for Phase II; however, potential receiver upgrades are under development as noted in Section 5.

b The reconfiguration takes approximately 2 weeks and happened approximately once per year between mid-2016 and mid-2018. Reconfigurations are anticipated to be twice yearly moving forward.

References

Abbott, B. P., et al. 2016, ApJ, 826, L13CrossRefGoogle Scholar
Abbott, B. P., et al. 2017, ApJ, 848, L12CrossRefGoogle Scholar
Ali, Z. S., et al. 2015, ApJ, 809, 61CrossRefGoogle Scholar
Barry, N., Hazelton, B., Sullivan, I., Morales, M. F., & Pober, J. C. 2016, MNRAS, 461, 3135CrossRefGoogle Scholar
Beardsley, A. P., et al. 2016, ApJ, 833, 102CrossRefGoogle Scholar
Benthem, P., et al. 2017, in 2017 XXXIInd General Assembly and Scientific Symposium of the International Union of Radio Science (URSI GASS), 1, doi: 10.23919/URSIGASS.2017.8104992CrossRefGoogle Scholar
Bowman, J. D., et al. 2013, PASA, 30, e031CrossRefGoogle Scholar
Chhetri, R., Morgan, J., Ekers, R. D., Macquart, J.-P., Sadler, E. M., Giroletti, M., Callingham, J. R., & Tingay, S. J. 2018, MNRAS, 474, 4937CrossRefGoogle Scholar
Condon, J. J. 1974, ApJ, 188, 279CrossRefGoogle Scholar
Croft, S., et al. 2016, ApJ, 820, L24CrossRefGoogle Scholar
DeBoer, D. R., et al. 2017, PASP, 129, 045001CrossRefGoogle Scholar
Dewdney, P. E., & Braun, R. 2016, Technical Report SKA-TEL-SKO-0000422, SKA1-LOW CONFIGURATION COORDINATES - COMPLETE SET. SKA OfficeGoogle Scholar
Dillon, J. S., & Parsons, A. R. 2016, ApJ, 826, 181CrossRefGoogle Scholar
Franzen, T. M. O., et al. 2016, MNRAS, 459, 3314CrossRefGoogle Scholar
Furlanetto, S. R., Oh, S. P., & Briggs, F. H. 2006, PhR, 433, 181Google Scholar
Gupta, Y. 2014, in Astronomical Society of India Conference Series, 441Google Scholar
Hurley-Walker, N., et al. 2017, MNRAS, 464, 1146CrossRefGoogle Scholar
Jacobs, D. C., etal. 2015, ApJ, 801, 51CrossRefGoogle Scholar
Jacobs, D. C., etal. 2016, ApJ, 825, 114CrossRefGoogle Scholar
Kaplan, D. L., et al. 2015a, ApJ, 809, L12CrossRefGoogle Scholar
Kaplan, D. L., et al. 2015b, ApJ, 814, L25CrossRefGoogle Scholar
Koopmans, L., et al. 2015, in Advancing Astrophysics with the Square Kilometre Array (AASKA14), 1Google Scholar
Li, W., et al. 2018, ApJ, 863, 170CrossRefGoogle Scholar
Liu, A., Tegmark, M., Morrison, S., Lutomirski, A., & Zaldarriaga, M. 2010, MNRAS, 408, 1029CrossRefGoogle Scholar
Loi, S. T., et al. 2015a, GeoRL, 42, 3707Google Scholar
Loi, S. T., et al. 2015b, RaSc, 50, 574Google Scholar
Mesinger, A., Furlanetto, S., & Cen, R. 2011, MNRAS, 411, 955CrossRefGoogle Scholar
Morgan, J. S., et al. 2018, MNRAS, 473, 2965CrossRefGoogle Scholar
Offringa, A. R., et al. 2015, PASA, 32, e008CrossRefGoogle Scholar
Offringa, A. R., et al. 2016, MNRAS, 458, 1057CrossRefGoogle Scholar
Ord, S. M., et al. 2015, PASA, 32, e006CrossRefGoogle Scholar
Paciga, G., et al. 2011, MNRAS, 413, 1174CrossRefGoogle Scholar
Parsons, A. R., et al. 2010, AJ, 139, 1468CrossRefGoogle Scholar
Parsons, A., Pober, J., McQuinn, M., Jacobs, D., & Aguirre, J. 2012, ApJ, 753, 81CrossRefGoogle Scholar
Patil, A. H., et al. 2016, MNRAS, 463, 4317CrossRefGoogle Scholar
Patil, A. H., et al. 2017, ApJ, 838, 65CrossRefGoogle Scholar
Prabu, T., et al. 2015, ExA, 39, 73Google Scholar
Sutinjo, A. T., et al. 2015, ITAP, 63, 5433Google Scholar
Tingay, S. J., et al. 2013a, PASA, 30, e007CrossRefGoogle Scholar
Tingay, S. J., et al. 2013b, AJ, 146, 103CrossRefGoogle Scholar
Tingay, S. J., etal. 2018, ApJ, 857, 11CrossRefGoogle Scholar
Tremblay, S. E., et al. 2015, PASA, 32, e005CrossRefGoogle Scholar
Tremblay, C. D., Hurley-Walker, N., Cunningham, M., Jones, P. A., Hancock, P. J., Wayth, R., & Jordan, C. H. 2017, MNRAS, 471, 4144CrossRefGoogle Scholar
Trott, C. M., et al. 2016, ApJ, 818, 139Google Scholar
van Haarlem, M. P., et al. 2013, A&A, 556, A2Google Scholar
Wayth, R., et al. 2017, PASA, 34, e034CrossRefGoogle Scholar
Wieringa, M. H. 1992, ExA, 2, 203Google Scholar
Wu, C., Wicenec, A., Pallot, D., & Checcucci, A. 2013, ExA, 36, 679Google Scholar
Zhang, X., et al. 2018, MNRAS, 477, 5167CrossRefGoogle Scholar
Zheng, H., et al. 2014, MNRAS, 445, 1084CrossRefGoogle Scholar
Zheng, Q., Wu, X.-P., Johnston-Hollitt, M., Gu, J.-h., & Xu, H. 2016, ApJ, 832, 190CrossRefGoogle Scholar
Figure 0

Figure 1. The compact configuration including the 72 new tiles arranged as two regular hexagonal arrays (filled squares). Note the size of the squares is not to scale.

Figure 1

Figure 2. The extended configuration including the 56 new long baseline tiles (filled squares). The inset rectangle shows the area bounded by Figure 1.

Figure 2

Figure 3. Block diagram of long baseline tile power, control, and signal transport system.

Figure 3

Figure 4. A Phase II long baseline tile. The batteries and BFIF systems are located under the solar panels.

Figure 4

Figure 5. Typical EoR power spectrum model at 150 MHz with associated noise levels available to the Phase I and Phase II arrays with 1 000-h observation. ‘Phase II 256’ represents the result from a future MWA upgrade where all 256 tiles are used simultaneously.

Figure 5

Figure 6. Snapshot monochromatic u, v coverage of the Phase II configurations at 154 MHz. The left panel shows the compact array (zenith pointed) and the right panel the extended array (zenith pointed).

Figure 6

Figure 7. Example cuts of the synthesised beam (zenith pointed) at 154 MHz for the Phase I and Phase II arrays with different robust weights. In this example, 1 MHz of bandwidth is used.

Figure 7

Figure 8. 185 MHz images of Fornax A taken with the Phase I (left) and Phase II (right) MWA. The ellipse in the lower left of each image shows the synthesised beam sizes, which are 2.6×2.3 and 1.0×0.89 arcmin.