The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

Context. The magnetic field in spiral galaxies is known to have a large-scale spiral structure along the galactic disk and is observed as X-shaped in the halo of some galaxies. While the disk field can be well explained by dynamo action, the 3-dimensional structure of the halo field and its physical nature is still unclear. Aims. As first steps towards understanding the halo fields, we want to clarify whether the observed X-shaped field is a wide-spread pattern in the halos of spiral galaxies and whether these halo fields are just turbulent fields ordered by compression or shear (anisotropic turbulent fields), or have a large-scale regular structure. Methods. The analysis of the Faraday rotation in the halo is the tool to discern anisotropic turbulent fields from large-scale magnetic fields. This, however, has been challenging until recently because of the faint halo emission in linear polarization. Our sensitive VLA broadband observations C-band and L-band of 35 spiral galaxies seen edge-on (called CHANG-ES) allowed us to perform RM-synthesis in their halos and to analyze the results. We further accomplished a stacking of the observed polarization maps of 28 CHANG-ES galaxies at C-band. Results. Though the stacked edge-on galaxies were of different Hubble types, star formation and interaction activities, the stacked image clearly reveals an X-shaped structure of the apparent magnetic field. We detected a large-scale (coherent) halo field in all 16 galaxies that have extended polarized intensity in their halos. We detected large-scale field reversals in all of their halos. In six galaxies they are along lines about vertical to the galactic midplane (vertical RMTL) with about 2 kpc separation. Only in NGC 3044 and possibly in NGC 3448 we observed vertical giant magnetic ropes (GMRs) similar to those detected recently in NGC 4631.

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$n\alpha$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $\tau$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

Magnetic fields are involved in every astrophysical process on every scale: from planetary and stellar interiors to neutron stars, stellar wind bubbles and supernova remnants; from the interstellar medium in galactic disks, nuclei, spiral arms and halos to the intracluster and intergalactic media. They are involved in essentially every particle acceleration process and are thus fundamental to non-thermal physics in the Universe. Key questions include the origin of magnetic fields, their evolution over cosmic time, the amplification and decay processes that modify their strength, and their impact on other processes such as star formation and galaxy evolution. Astrophysical plasmas provide a unique laboratory for testing magnetic dynamo theory. The study of magnetic fields requires observations that span the wavelength range from radio through infrared, optical, UV, X-ray, and gamma-ray. Canada has an extremely strong record of research in cosmic magnetism, and has a significant leadership role in several ongoing and upcoming global programs. This white paper will review the science questions to be addressed in the study of cosmic magnetic fields and will describe the observational and theoretical opportunities and challenges afforded by the telescopes and modelling capabilities of today and tomorrow.

We present new narrow-band H-alpha imaging for 24 nearby edge-on galaxies in the CHANG-ES survey. We use the images in conjunction with WISE 22 micron imaging of the sample to estimate improved star formation rates (SFRs) using the updated recipe from Vargas et al. (2018). We explore correlations between the updated star formation properties and radio continuum scale heights, scale lengths, and diameters, measured in Krause et al. (2018). We find a newly discovered correlation between SFR and radio scale height that did not exist using mid-IR only SFR calibrations. This implies that a mid-IR extinction correction should be applied to SFR calibrations when used in edge-on galaxies, due to attenuation by dust. The updated SFR values also show newly discovered correlations with radio scale length and radio diameter, implying that the previously-measured relationship between radio scale height and radio diameter originates from star formation within the disk. We also identify a region of star formation located at extreme distance from the disk of NGC 4157, possibly ionized by a single O5.5 V star. This region is spatially coincident with an XUV disk feature, as traced by GALEX NUV imaging. We theorize that the star formation feature arose due to gravitational instability within gas from an accretion event. New H-alpha images from this work can be found at the CHANG-ES data release web site, https://www.queensu.ca/changes.

The NRAO VLA Sky Survey (NVSS) Rotation Measure (RM) catalogue has enabled numerous studies in cosmic magnetism, and will continue being a unique dataset complementing future polarisation surveys. Robust comparisons with these new surveys will however require further understandings in the systematic effects present in the NVSS RM catalogue. In this paper, we make careful comparisons between our new on-axis broadband observations with the Karl G. Jansky Very Large Array and the NVSS RM results for 23 sources. We found that two unpolarised sources were reported as polarised at about 0.5% level in the RM catalogue, and noted significant differences between our newly derived RM values and the catalogue values for the remaining 21 sources. These discrepancies are attributed to off-axis instrumental polarisation in the NVSS RM catalogue. By adopting the 0.5% above as the typical off-axis instrumental polarisation amplitude, we quantified its effect on the reported RMs with a simulation, and found that on average the RM uncertainties in the catalogue have to be increased by $\approx$ 10% to account for the off-axis instrumental polarisation effect. This effect is more substantial for sources with lower fractional polarisation, and is a function of the source's true RM. Moreover, the distribution of the resulting RM uncertainty is highly non-Gaussian. With the extra RM uncertainty incorporated, we found that the RM values from the two observations for most (18 out of 21) of our polarised targets can be reconciled. The remaining three are interpreted as showing hints of time variabilities in RM.

The NRAO VLA Sky Survey (NVSS) Rotation Measure (RM) catalogue is invaluable for the study of cosmic magnetism. However, the RM values reported in it can be affected by n$\pi$-ambiguity, resulting in deviations of the reported RM from the true values by multiples of +-652.9 rad m-2. We therefore set off to observationally constrain the fraction of sources in the RM catalogue affected by this ambiguity. New broadband spectro-polarimetric observations were performed with the Karl G. Jansky Very Large Array (VLA) at 1--2 GHz, with 23 n$\pi$-ambiguity candidates selected by their peculiarly high |RM| values. We identified nine sources with erroneous RM values due to n$\pi$-ambiguity and 11 with reliable RM values. In addition, we found two sources to be unpolarised and one source to be inconsistent with neither n$\pi$-ambiguity nor reliable RM cases. By comparing the statistical distributions of the above two main classes, we devised a measure of how much a source's RM deviates from that of its neighbours: $\Delta/\sigma$, which we found to be a good diagnostic of n$\pi$-ambiguity. With this, we estimate that there are at least 50 sources affected by n$\pi$-ambiguity among the 37,543 sources in the catalogue. Finally, we explored the Faraday complexities of our sources revealed by our broadband observations.

With rotation measure (RM) towards 37,543 polarised sources, the Taylor et al. (2009) RM catalogue has been widely exploited in studies of the foreground magneto-ionic media. However, due to limitations imposed by observations in survey mode in the narrowband era, the listed RM values are inevitably affected by various systematic effects. With new Karl G. Jansky Very Large Array (VLA) broadband spectro-polarimetric observations at L-band, we set off to observationally examine the robustness of the Taylor catalogue. This would facilitate combinations and comparisons of it with results from current and future polarisation surveys such as Polarization Sky Survey of the Universe's Magnetism (POSSUM), VLA Sky Survey (VLASS), and the eventual Square Kilometre Array (SKA). Our on-axis pointed observations, in conjunction with simulations, allowed us to estimate the impact of off-axis polarisation leakage on the measured RM values. This demonstrates the importance to properly calibrate for the off-axis leakage terms in future all-sky polarisation surveys, in order to obtain high fidelity polarisation information from sources down to low fractional polarisation.

Anderson et al. (2017) identified seventy six candidate supernova remnants (SNRs) using data from The HI, OH, Recombination line survey of the Milky Way (THOR). The spectral index and polarization properties can help distinguish between SNRs and H II regions, which are often confused. We confirm two SNR candidates using spectral index data and morphology. However, we observe that the fractional linear polarization cannot distinguish between SNRs and H II regions, likely due to contamination by diffuse Galactic synchrotron emission. We also comment on the association of SNR candidates with pulsars through geometric and age considerations.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.

Observations of Galactic H I gas for seven targeted regions at intermediate Galactic latitude are presented at 1' angular resolution using data from the DRAO Synthesis Telescope (ST) and the Green Bank Telescope (GBT). The DHIGLS data are the most extensive arcminute resolution measurements of the diffuse atomic interstellar medium beyond those in the Galactic plane. The acquisition, reduction, calibration, and mosaicking of the DRAO ST data and the cross calibration and incorporation of the short-spacing information from the GBT are described. The high quality of the resulting DHIGLS products enables a variety of new studies in directions of low Galactic column density. We analyze the angular power spectra of maps of the integrated H I emission (column density) from the data cubes for several distinct velocity ranges. Fitting power spectrum models based on a power law, but including the effects of the synthesized beam and noise at high spatial frequencies, we find exponents ranging from -2.5 to -3.0.. Power spectra of maps of the centroid velocity for these components give similar results. These exponents are interpreted as being representative of the 3D density and 3D velocity fields of the atomic gas, respectively. We find evidence for dramatic changes in the H I structures in channel maps over even small changes in velocity. This narrow line emission has counterparts in absorption spectra against bright background radio sources, quantifying that the gas is cold and dense and can be identified as the cold neutral medium phase. Fully reduced DHIGLS H I data cubes and other data products are available at www.cita.utoronto.ca/DHIGLS.

The Cygnus Loop is a well-studied supernova remnant (SNR) that has been observed across the electromagnetic spectrum. Although widely believed to be an SNR shell with a blow- out region in the south, we consider the possibility that this object is two SNRs projected along the same line-of-sight by using multi-wavelength images and modelling. Our results show that a model of two objects including some overlap region/interaction between the two objects has the best match to the observed data.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

We investigate the magnetic field in an intermediate-velocity filament for which the H$\alpha$ intensity in the WHAM survey correlates with excess Faraday rotation of extragalactic radio sources over the length of the filament from b ~ 20 degr. to b ~ 55 degr. The density-weighted mean magnetic field is 2.8 +/- 0.8 microgauss, derived from rotation measures and an empirical relation between H-alpha emission measure and dispersion measure from Berkhuijsen et al. (2006). In view of the uncertainties in the derived magnetic field strength, we propose an alternative use of the available data, rotation measure and emission measure, to derive a lower limit to the Alfven speed, weighted by electron density n_e^1.5. We find lower limits to the Alfven speed that are comparable to, or larger than the sound speed in a 10^4 K plasma, and conclude that the magnetic field is dynamically important. We discuss the role of intermediate-velocity gas as a locus of Faraday rotation in the interstellar medium, and propose this lower limit to the Alfven speed may also be applicable to Faraday rotation by galaxy clusters.

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

We have mapped cold atomic gas in 21cm line HI self-absorption (HISA) at arcminute resolution over more than 90% of the Milky Way's disk. To probe the formation of H2 clouds, we have compared our HISA distribution with CO J=1-0 line emission. Few HISA features in the outer Galaxy have CO at the same position and velocity, while most inner-Galaxy HISA has overlapping CO. But many apparent inner-Galaxy HISA-CO associations can be explained as chance superpositions, so most inner-Galaxy HISA may also be CO-free. Since standard equilibrium cloud models cannot explain the very cold HI in many HISA features without molecules being present, these clouds may instead have significant CO-dark H2.

We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

(abridged) Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground reconstruction is not a good estimate for the extragalactic contributions. We point out a degeneracy between the contributions to the observed values due to extragalactic magnetic fields and observational noise and comment on the dangers of over-interpreting an estimate without taking into account its uncertainty information. To overcome these difficulties, we develop an extended reconstruction algorithm based on the assumption that the observational uncertainties are accurately described for a subset of the data, which can overcome the degeneracy with the extragalactic contributions. We present a probabilistic derivation of the algorithm and demonstrate its performance using a simulation, yielding a high quality reconstruction of the Galactic Faraday rotation foreground, a precise estimate of the typical extragalactic contribution, and a well-defined probabilistic description of the extragalactic contribution for each data point. We then apply this reconstruction technique to a catalog of Faraday rotation observations. We vary our assumptions about the data, showing that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m^2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.

Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking.

Stacking polarized radio emission in SKA surveys provides statistical information on large samples that is not accessible otherwise due to limitations in sensitivity, source statistics in small fields, and averaging over frequency (including Faraday synthesis). Polarization is a special case because one obvious source of stacking targets is the Stokes I source catalog, possibly in combination with external catalogs, for example an SKA HI survey or a non-radio survey. We point out the significance of stacking sub-samples selected by additional observable parameters to investigate relations that reveal more about the physics of the source. Applications of stacking polarization include, but are not limited to, obtaining in a statistical sense polarization information to the detection limit in total intensity, depolarization as a function of cosmic time at consistent source-frame wavelengths, magnetic field properties in objects with a low radio luminosity such as dwarf and low-surface-brightness galaxies, and investigating potential correlations of observable parameters with the average magnetic field direction in a sample. We also point out the potential use of stacking in validating the polarization calibration of a survey. While stacking is flexible in terms of survey definition, we discuss optimal survey parameters for the science experiments presented, as well as computing and archiving requirements.

Magnetic fields play an important role in shaping the structure and evolution of the interstellar medium (ISM) of galaxies, but the details of this relationship remain unclear. With SKA1, the 3D structure of galactic magnetic fields and its connection to star formation will be revealed. A highly sensitive probe of the internal structure of the magnetoionized ISM is the partial depolarization of synchrotron radiation from inside the volume. Different configurations of magnetic field and ionized gas within the resolution element of the telescope lead to frequency-dependent changes in the observed degree of polarization. The results of spectro-polarimetric observations are tied to physical structure in the ISM through comparison with detailed modeling, supplemented with the use of new analysis techniques that are being actively developed and studied within the community such as Rotation Measure Synthesis. The SKA will enable this field to come into its own and begin the study of the detailed structure of the magnetized ISM in a sample of nearby galaxies, thanks to its extraordinary wideband capabilities coupled with the combination of excellent surface brightness sensitivity and angular resolution.

(abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying science: a) an average RM weighted by polarized intensity, RMwtd, b) the separation $\Delta\phi$ of two Faraday components and c) the reduced chi-squared. Based on the current test data of signal to noise ratio of about 32, we find that: (1) When only one Faraday thin component is present, most methods perform as expected, with occasional failures where two components are incorrectly found; (2) For two Faraday thin components, QU-fitting routines perform the best, with errors close to the theoretical ones for RMwtd, but with significantly higher errors for $\Delta\phi$. All other methods including standard Faraday synthesis frequently identify only one component when $\Delta\phi$ is below or near the width of the Faraday point spread function; (3) No methods, as currently implemented, work well for Faraday thick components due to the narrow bandwidth; (4) There exist combinations of two Faraday components which produce a large range of acceptable fits and hence large uncertainties in the derived single RMs; in these cases, different RMs lead to the same Q, U behavior, so no method can recover a unique input model.

We present polarisation properties at $1.4\,$GHz of two separate extragalactic source populations: passive quiescent galaxies and luminous quasar-like galaxies. We use data from the {\it Wide-Field Infrared Survey Explorer} data to determine the host galaxy population of the polarised extragalactic radio sources. The quiescent galaxies have higher percentage polarisation, smaller radio linear size, and $1.4\,$GHz luminosity of $6\times10^{21}<L_{\rm 1.4}<7\times10^{25}\,$W Hz$^{-1}$, while the quasar-like galaxies have smaller percentage polarisation, larger radio linear size at radio wavelengths, and a $1.4\,$GHz luminosity of $9\times10^{23}<L_{\rm 1.4}<7\times10^{28}\,$W Hz$^{-1}$, suggesting that the environment of the quasar-like galaxies is responsible for the lower percentage polarisation. Our results confirm previous studies that found an inverse correlation between percentage polarisation and total flux density at $1.4\,$GHz. We suggest that the population change between the polarised extragalactic radio sources is the origin of this inverse correlation and suggest a cosmic evolution of the space density of quiescent galaxies. Finally, we find that the extragalactic contributions to the rotation measures (RMs) of the nearby passive galaxies and the distant quasar-like galaxies are different. After accounting for the RM contributions by cosmological large-scale structure and intervening Mg\,{II} absorbers we show that the distribution of intrinsic RMs of the distant quasar-like sources is at most four times as wide as the RM distribution of the nearby quiescent galaxies, if the distribution of intrinsic RMs of the WISE-Star sources itself is at least several rad m$^{-2}$ wide.

The new broad-band capabilities of large radio interferometers such as the GMRT and JVLA allow for long-integration mosaic imaging observations to create ultra-deep full-polarization images of the sky over wide frequency ranges. Achieving rms sensitivities of order 1 $\mu$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio galaxies and Active Galactic Nuclei. We present initial results from radio sources revealed with deep mosaicking observations with the GMRT and JVLA at respectively 0.6 and 5 GHz, and evidence that the $\mu$Jy sensitivity level marks the transition to detection of polarized emission from a population of sources dominated by emission from magnetic fields in the disks of starburst and normal galaxies.

We present stacking polarized intensity as a means to study the polarization of sources that are too faint to be detected individually in surveys of polarized radio sources. Stacking offers not only high sensitivity to the median signal of a class of radio sources, but also avoids a detection threshold in polarized intensity, and therefore an arbitrary exclusion of source with a low percentage of polarization. Correction for polarization bias is done through a Monte Carlo analysis and tested on a simulated survey. We show that the non-linear relation between the real polarized signal and the detected signal requires knowledge of the shape of the distribution of fractional polarization, which we constrain using the ratio of the upper quartile to the lower quartile of the distribution of stacked polarized intensities. Stacking polarized intensity for NVSS sources down to the detection limit in Stokes I, we find a gradual increase in median fractional polarization that is consistent with a trend that was noticed before for bright NVSS sources, but is much more gradual than found by previous deep surveys of radio polarization. Consequently, the polarized radio source counts derived from our stacking experiment predict fewer polarized radio sources for future surveys with the Square Kilometre Array and its pathfinders.

(Abridged)One of the fundamental challenges for astrophysics in the 21st century is finding a way to untangle the physical processes that govern galaxy formation and evolution. Given the importance and scope of this problem, the multi-wavelength astronomical community has used the past decade to build up a wealth of information over specific extragalactic deep fields to address key questions in galaxy formation and evolution. These fields generally cover at least 10square degrees to facilitate the investigation of the rarest, typically most massive, galaxies and AGN. Furthermore, such areal coverage allows the environments to be fully accounted for, thereby linking the single halo to the two-halo terms in the halo occupation distribution. Surveys at radio wavelengths have begun to lag behind those at other wavelengths, especially in this medium-deep survey tier. However, the survey speed offered by the JVLA means that we can now reach a point where we can begin to obtain commensurate data at radio wavelengths to those which already exists from the X-ray through to the far-infrared over ~10 square degrees. We therefore present the case for a 10 square degree survey to 1.5uJy at L-band in A or B Array, requiring ~4000 hours to provide census of star-formation and AGN-accretion activity in the Universe. For example, the observations will allow galaxies forming stars at 10Msolar/yr to be detected out to z~1 and luminous infrared galaxies (1000Msolar/yr to be found out to z~6. Furthermore, the survey area ensures that we will have enough cosmic volume to find these rare sources at all epochs. The bandwidth will allow us to determine the polarisation properties galaxies in the high-redshift Universe as a function of stellar mass, morphology and redshift.

A VLA Sky Survey of the extragalactic sky at S band (2-4 GHz) with polarization information can uniquely probe the magneto-ionic medium in a wide range of astrophysical environments over cosmic time. For a shallow all-sky survey, we expect to detect over 4 million sources in total intensity $>$ 0.45 mJy beam$^{-1}$ and over 2.2$\times$10$^5$ sources in polarized intensity. With these new observations, we expect to discover new classes of polarized radio sources in very turbulent astrophysical environments and those with extreme values of Faraday depth. Moreover, by determining reliable Faraday depths and by modeling depolarization effects, we can derive properties of the magneto-ionic medium associated with AGNs, absorption line systems and galaxies, addressing the following unresolved questions: (1) What is the covering fraction, the degree of turbulence and the origin of absorption line systems? (2) What is the thermal content in AGNs and radio galaxies? (3) How do AGNs and galaxies evolve over cosmic time? (4) What causes the increase in percentage polarization with decreasing flux densities at the low flux density end of the polarized source count? (5) What is the growth rate of large-scale magnetic fields in galaxies?

We use previously published high-resolution synchrotron polarization data to perform an angular dispersion analysis with the aim of charactering magnetized turbulence in M51. We first analyze three distinct regions (the center of the galaxy, and the northwest and southwest spiral arms) and can clearly discern the turbulent correlation length scale from the width of the magnetized turbulent correlation function for two regions and detect the imprint of anisotropy in the turbulence for all three. Furthermore, analyzing the galaxy as a whole allows us to determine a two-dimensional Gaussian model for the magnetized turbulence in M51. We measure the turbulent correlation scales parallel and perpendicular to the local mean magnetic field to be, respectively, delta_{para} = 98 +/- 5 pc and delta_{perp} = 54 +/- 3 pc, while the turbulent to ordered magnetic field strength ratio is found to be Bt/B0 = 1.01 +/- 0.04. These results are consistent with those of Fletcher et al. (2011), who performed a Faraday rotation dispersion analysis of the same data, and our detection of anisotropy is consistent with current magnetized turbulence theories.

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

We present a study of the Milky Way disk and halo magnetic field, determined from observations of Faraday rotation measure (RM) towards 641 polarized extragalactic radio sources in the Galactic longitude range 100-117 degs, within 30 degs of the Galactic plane. For |b| < 15 degs, we observe a symmetric RM distribution about the Galactic plane. This is consistent with a disk field in the Perseus arm of even parity across the Galactic mid-plane. In the range 15<|b|<30 degs, we find median rotation measures of -15+/-4 rad/m^2 and -62+/-5 rad/m^2 in the northern and southern Galactic hemispheres, respectively. If the RM distribution is a signature of the large-scale field parallel to the Galactic plane, this suggests that the halo magnetic field toward the outer Galaxy does not reverse direction across the mid-plane. The variation of RM as a function of Galactic latitude in this longitude range is such that RMs become more negative at larger |b|. This is consistent with an azimuthal magnetic field of strength 2 microGauss (7 microGauss) at a height 0.8-2 kpc above (below) the Galactic plane between the local and the Perseus spiral arm. We propose that the Milky Way could possess spiral-like halo magnetic fields similar to those observed in M51.

We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of pi. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior with a power law index of -2.14 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.

GALFACTS is a large-area spectro-polarimetric survey on the Arecibo Radio telescope. It uses the seven-beam focal plane feed array receiver system (ALFA) to carry out an imaging survey project of the 12,700 square degrees of sky visible from Arecibo at 1.4 GHz with 8192 spectral channels over a bandwidth of 300 MHz sampled at 1 millisecond. The aggregate data rate is 875 MB/s. GALFACTS observations will create full-Stokes image cubes at an angular resolution of 3.5' with a band-averaged sensitivity of 90 $\mu$Jy, allowing sensitive imaging of polarized radiation and Faraday Rotation Measure from both diffuse emission and extragalactic sources. GALFACTS is a scientific pathfinder to the SKA in the area of cosmic magnetism. Key to magnetism science with the SKA is the technique of RM synthesis. The technique of RM synthesis is introduced and we discuss practical aspects of RM synthesis including efficient computational techniques and detection thresholds in the resulting Faraday spectrum. We illustrate the use of the technique by presenting the current development of the RM synthesis pipeline for GALFACTS and present early results.

Faraday rotation measure synthesis is a method for analyzing multichannel polarized radio emissions, and it has emerged as an important tool in the study of galactic and extra-galactic magnetic fields. The method requires the recovery of the Faraday dispersion function from measurements restricted to limited wavelength ranges, which is an ill-conditioned deconvolution problem. Here, we discuss a recovery method, which assumes a sparse approximation of the Faraday dispersion function in an over-complete dictionary of functions. We discuss the general case, when both thin and thick components are included in the model, and we present the implementation of a greedy deconvolution algorithm. We illustrate the method with several numerical simulations that emphasize the effect of the covered range and sampling resolution in the Faraday depth space, and the effect of noise on the observed data.

Detection thresholds in polarized intensity and polarization bias correction are investigated for surveys where the polarization information is obtained from RM synthesis. Considering unresolved sources with a single rotation measure, a detection threshold of $8 \sigma_{QU}$ applied to the Faraday spectrum will retrieve the RM with a false detection rate less than $10^{-4}$, but polarized intensity is more strongly biased than Ricean statistics suggest. For a detection threshold of $5 \sigma_{QU}$, the false detection rate increases to ~4%, depending also on $\lambda^2$ coverage and the extent of the Faraday spectrum. Non-Gaussian noise in Stokes Q and U due to imperfect imaging and calibration can be represented by a distribution that is the sum of a Gaussian and an exponential. The non-Gaussian wings of the noise distribution increase the false detection rate in polarized intensity by orders of magnitude. Monte-Carlo simulations assuming non-Gaussian noise in Q and U, give false detection rates at $8 \sigma_{QU}$ similar to Ricean false detection rates at $4.9 \sigma_{QU}$.

The Dominion Radio Astrophysical Observatory Deep Field polarization study has been matched with the Spitzer Wide-Area Infrared Extragalactic survey of the European Large Area Infrared Space Observatory Survey North 1 field. We have used VLA observations with a total intensity rms of 87 microJy beam^-1 to match SWIRE counterparts to the radio sources. Infrared color analysis of our radio sample shows that the majority of polarized sources are elliptical galaxies with an embedded active galactic nucleus. Using available redshift catalogs, we found 429 radio sources of which 69 are polarized with redshifts in the range of 0.04 < z <3.2. We find no correlation between redshift and percentage polarization for our sample. However, for polarized radio sources, we find a weak correlation between increasing percentage polarization and decreasing luminosity.

We have determined 194 Faraday rotation measures (RMs) of polarized extragalactic radio sources using new, multi-channel polarization observations at frequencies around 1.4 GHz from the Very Large Array (VLA) in the Galactic plane at $17^\circ \leq l \leq 63^\circ$ and $205^\circ \leq l \leq 253^\circ$. This catalog fills in gaps in the RM coverage of the Galactic plane between the Canadian Galactic Plane Survey and Southern Galactic Plane Survey. Using this catalog we have tested the validity of recently-proposed axisymmetric and bisymmetric models of the large-scale (or regular) Galactic magnetic field, and found that of the existing models we tested, an axisymmetric spiral model with reversals occurring in rings (as opposed to along spiral arms) best matched our observations. Building on this, we have performed our own modeling, using RMs from both extragalactic sources and pulsars. By developing independent models for the magnetic field in the outer and inner Galaxy, we conclude that in the inner Galaxy, the magnetic field closely follows the spiral arms, while in the outer Galaxy, the field is consistent with being purely azimuthal.Furthermore, the models contain no reversals in the outer Galaxy, and together seem to suggest the existence of a single reversed region that spirals out from the Galactic center.

An analysis of structure in rotation measure (RM) across the sky based on the RM catalog of Taylor et al. (2009) is presented. Several resolved RM structures are identified with structure in the local ISM, including radio loops I, II, and III, the Gum nebula, and the Orion-Eridanus super bubble. Structure functions (SFs) of RM are presented for selected areas, and maps of SF amplitude and slope across the sky are compared with H-alpha intensity and diffuse polarized intensity. RM variance on an angular scale of 1 degree is correlated with length of the line of sight through the Galaxy, with a contribution from local structures. The slope of the SFs is less concentrated to the Galactic plane and less correlated with length of the line of sight through the Galaxy, suggesting a more local origin for RM structure on angular scales ~ 10 degrees. The RM variance is a factor ~2 higher towards the SGP than towards the NGP, reflecting a more wide-spread asymmetry between the northern and southern Galactic hemispheres. Depolarization of diffuse Galactic synchrotron emission at latitudes < 30 degrees can be explained largely by Faraday dispersion related to small-scale variance in RM, but the errors allow a significant contribution from differential Faraday rotation along the line of sight.

Observational evidence for the existence of cosmic magnetic fields in intergalactic space coherent over Mpc scales is presented. Using an unprecedentedly large sample of Faraday rotation measures of radio sources from the NRAO VLA Sky Survey data and the photometric redshift galaxy catalog from the sixth Data Release of the Sloan Digital Sky Survey (SDSS), we measure the cross-correlations between the rotation measures and the galaxy density field distributed along the source sightlines. It is shown that the rotation measures for sightlines passing through high density regions at separation r >= 1 Mpc/h from the locations of background radio sources are significantly enhanced. We discuss possible generators of this enhancement and interpret it to be intergalactic magnetic fields coherent over 1 Mpc/h with mean field strength B~30 nG.

An electromagnetic wave can be uniquely characterized by the four Stokes parameters: I, Q, U, and V. Typical observations in astronomy rely solely on total intensity measurements of the incoming radiation (Stokes I). However, a significant amount of information both about the emitting region and the propagation path is carried in the remaining Stokes parameters. These data provide a means to observe parts of the interstellar medium which remain invisible in Stokes I, at any wavelength. For example, when an electromagnetic wave propagates through a region containing free electrons and a magnetic field, the plane of polarisation of the wave will rotate - an effect recorded only in Stokes Q and U. The interstellar medium of the Galaxy is such a region, containing free electrons (observed as HII) and a magnetic field of a few microgauss. By imaging in Stokes Q and U we are able to observe signatures of magnetic field perturbations from the small scale (tens of pc) to the large scale (kpc). In this paper, we review the status of Canadian polarisation studies of cosmic magnetic fields and discuss the leading role Canada is playing in polarsation studies around the world.

We present results of deep polarization imaging at 1.4 GHz with the Dominion Radio Astrophysical Observatory as part of the DRAO Planck Deep Fields project. This deep extragalactic field covers 15.16 square degrees centered at RA = 16h 14m and DEC = 54d 56', has an angular resolution of 42" x 62" at the field center, and reaches a sensitivity of 55 microJy/beam in Stokes I and 45 microJy/beam in Stokes Q and U. We detect 958 radio sources in Stokes I of which 136 are detected in polarization. We present the Euclidean-normalized polarized differential source counts down to 400 microJy. These counts indicate that sources have a higher degree of fractional polarization at fainter Stokes I flux density levels than for brighter sources, confirming an earlier result. We find that the majority of our polarized sources are steep-spectrum objects with a mean spectral index of -0.77, and there is no correlation between fractional polarization and spectral index. We also matched deep field sources to counterparts in the Faint Images of the Radio Sky at Twenty Centimeters catalogue. Of the polarized sources, 77% show structure at the arc-second scale whereas only 38% of the sources with no detectable polarization show such structure. The median fractional polarization is for resolved sources is 6.8%, while it is 4.4% for compact objects. The polarized radio sources in our deep field are predominantly those sources which are resolved and show the highest degrees of fractional polarization, indicating that the lobe dominated structure may be the source of the highly polarized sources. These resolved radio galaxies dominate the polarized source counts at P_0 = sqrt(Q^2 + U^2) < 3 mJy.

We present three-dimensional magnetohydrodynamic (MHD) simulations of superbubbles, to study the importance of MHD effects in the interpretation of images from recent surveys of the Galactic plane. These simulations focus mainly on atmospheres defined by an exponential density distribution and the Dickey & Lockman (1990) density distribution. In each case, the magnetic field is parallel to the Galactic plane and we investigate cases with either infinite scale height (constant magnetic field) or a constant ratio of gas pressure to magnetic pressure. The three-dimensional structure of superbubbles in these simulations is discussed with emphasis on the axial ratio of the cavity as a function of magnetic field strength and the age of the bubble. We investigate systematic errors in the age of the bubble and scale height of the surrounding medium that may be introduced by modeling the data with purely hydrodynamic models. Age estimates derived with symmetric hydrodynamic models fitted to an asymmetric magnetized superbubble can differ by up to a factor of four, depending on the direction of the line of sight. The scale height of the surrounding medium based on the Kompaneets model may be up to 50% lower than the actual scale height. We also present the first ever predictions of Faraday rotation by a magnetized superbubble based on three-dimensional MHD simulations. We emphasize the importance of MHD effects in the interpretation of observations of superbubbles.

We present results from simulations of the extragalactic polarized sky at 1.4 GHz. As the basis for our polarization models, we use a semi-empirical simulation of the extragalactic total intensity (Stokes I) continuum sky developed at the University of Oxford (http://scubed.physics.ox.ac.uk) under the European SKA Design Study (SKADS) initiative, and polarization distributions derived from analysis of polarization observations. By considering a luminosity dependence for the polarization of AGN, we are able to fit the 1.4 GHz polarized source counts derived from the NVSS and the DRAO ELAIS N1 deep field survey down to approximately 1 mJy. This trend is confirmed by analysis of the polarization of a complete sample of bright AGN. We are unable to fit the additional flattening of the polarized source counts from the deepest observations of the ELAIS N1 survey, which go down to ~0.5 mJy. Below 1 mJy in Stokes I at 1.4 GHz, starforming galaxies become an increasingly important fraction of all radio sources. We use a spiral galaxy integrated polarization model to make realistic predictions of the number of polarized sources at microJy levels in polarized flux density and hence, realistic predictions of what the next generation radio telescopes such as ASKAP, other SKA pathfinders and the SKA itself will see.

Three recent surveys of 21-cm line emission in the Galactic plane, combining single dish and interferometer observations to achieve resolution of 1 arcmin to 2 arcmin, 1 km/s, and good brightness sensitivity, have provided some 650 absorption spectra with corresponding emission spectra for study of the distribution of warm and cool phase H I in the interstellar medium. These emission-absorption spectrum pairs are used to study the temperature of the interstellar neutral hydrogen in the outer disk of the Milky Way, outside the solar circle, to a radius of 25 kpc. The cool neutral medium is distributed in radius and height above the plane with very similar parameters to the warm neutral medium. In particular, the ratio of the emission to the absorption, which gives the mean spin temperature of the gas, stays nearly constant with radius to 25 kpc radius. This suggests that the mixture of cool and warm phases is a robust quantity, and that the changes in the interstellar environment do not force the H I into a regime where there is only one temperature allowed. The mixture of atomic gas phases in the outer disk is roughly 15% to 20% cool (40 K to 60 K), the rest warm, corresponding to mean spin temperature 250 to 400 K. The Galactic warp appears clearly in the absorption data, and other features on the familiar longitude-velocity diagram have analogs in absorption with even higher contrast than for 21-cm emission. In the third and fourth Galactic quadrants the plane is quite flat, in absorption as in emission, in contrast to the strong warp in the first and second quadrants. The scale height of the cool gas is similar to that of the warm gas, and both increase with Galactic radius in the outer disk.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

We present integrated polarization properties of nearby spiral galaxies at 4.8 GHz, and models for the integrated polarization of spiral galaxy disks as a function of inclination. Spiral galaxies in our sample have observed integrated fractional polarization in the range < 1% to 17.6%. At inclinations less than 50 degrees, the fractional polarization depends mostly on the ratio of random to regular magnetic field strength. At higher inclinations, Faraday depolarization associated with the regular magnetic field becomes more important. The observed degree of polarization is lower (<4%) for more luminous galaxies, in particular those with L_{4.8} > 2 x 10^{21} W/Hz. The polarization angle of the integrated emission is aligned with the apparent minor axis of the disk for galaxies without a bar. In our axially symmetric models, the polarization angle of the integrated emission is independent of wavelength. Simulated distributions of fractional polarization for randomly oriented spiral galaxies at 4.8 GHz and 1.4 GHz are presented. We conclude that polarization measurements, e.g. with the SKA, of unresolved spiral galaxies allow statistical studies of the magnetic field in disk galaxies using large samples in the local universe and at high redshift. As these galaxies behave as idealized background sources without internal Faraday rotation, they can be used to detect large-scale magnetic fields in the intergalactic medium.

Predictions of the number of faint polarised radio sources that can be detected by SKA pathfinder telescopes and the SKA depend on the polarisation properties of radio sources with a total flux density around 1 mJy. Total intensity source counts suggest a transition in the dominant population from AGN to galaxies around this flux density, and the properties of brighter radio sources may not be representative for this fainter population. We show that unresolved spiral galaxies can be highly polarised radio sources, up to ~ 20% polarised at 4.8 GHz. This result is partly based on observations of nearby galaxies, including galaxies with significant deviations from axial symmetry and other peculiarities. A first analysis of polarised source counts divided into steep-spectrum AGN, flat-spectrum AGN and star forming galaxies is presented, including a prediction of polarised source counts to microjansky levels.

We report depolarization of extragalactic sources in the NRAO VLA Sky Survey (NVSS) by local structures in the interstellar medium. The sky density of polarized sources drops by a factor 2 - 4 in regions with angular scales approx. 10 degrees, implying up to 40% depolarization on average per source. Some of these polarization shadows are associated with HII regions, but three are associated with regions of depolarized diffuse Galactic emission. The absence of a correlation between the depth of polarization shadows and H-alpha intensity suggests that some shadows are related to structure in the magnetic field. At least some polarization shadows are caused by partial bandwidth depolarization in the NVSS. Alternatively, some may be caused by regions with small scale (< 1 arcsec) variations in rotation measure.

We present deep polarimetric observations at 1420 MHz of the European Large Area ISO Survey North 1 region (ELAIS N1) as part of the Dominion Radio Astrophysical Observatory Planck Deep Fields project. By combining closely spaced aperture synthesis fields, we image a region of 7.43 square degrees to a maximum sensitivity in Stokes Q and U of 78 microJy/beam, and detect 786 compact sources in Stokes I. Of these, 83 exhibit polarized emission. We find that the differential source counts (log N - log p) for polarized sources are nearly constant down to p > 500 microJy, and that these faint polarized radio sources are more highly polarized than the strong source population. The median fractional polarization is (4.8 +/- 0.7)% for polarized sources with Stokes I flux density between 1 and 30 mJy; approximately three times larger than sources with I > 100 mJy. The majority of the polarized sources have been identified with galaxies in the Spitzer Wide Area Infrared Extragalactic Survey (SWIRE) image of ELAIS N1. Most of the galaxies occupy regions in the IRAC 5.8/3.6 micron vs. 8.0/4.5 micron color-color diagram associated with dusty AGNs, or with ellipticals with an aging stellar population. A few host galaxies have colors that suggests significant PAH emission in the near-infrared. A small fraction, 12%, of the polarized sources are not detected in the SWIRE data. None of the polarized sources in our sample appears to be associated with an actively star-forming galaxy.

HI 21cm-line self-absorption (HISA) reveals the shape and distribution of cold atomic clouds in the Galactic disk. Many of these clouds lack corresponding CO emission, despite being colder than purely atomic gas in equilibrium models. HISA requires background line emission at the same velocity, hence mechanisms that can produce such backgrounds. Weak, small-scale, and widespread absorption is likely to arise from turbulent eddies, while strong, large-scale absorption appears organized in cloud complexes along spiral arm shocks. In the latter, the gas may be evolving from an atomic to a molecular state prior to star formation, which would account for the incomplete HISA-CO agreement.

The VLA Galactic Plane Survey (VGPS) is a survey of HI and 21-cm continuum emission in the Galactic plane between longitude 18 degrees 67 degr. with latitude coverage from |b| < 1.3 degr. to |b| < 2.3 degr. The survey area was observed with the Very Large Array (VLA) in 990 pointings. Short-spacing information for the HI line emission was obtained by additional observations with the Green Bank Telescope (GBT). HI spectral line images are presented with a resolution of 1 arcmin x 1 arcmin x 1.56 km/s (FWHM) and rms noise of 2 K per 0.824 km/s channel. Continuum images made from channels without HI line emission have 1 arcmin (FWHM) resolution. VGPS images are compared with images from the Canadian Galactic Plane Survey (CGPS) and the Southern Galactic Plane Survey (SGPS). In general, the agreement between these surveys is impressive, considering the differences in instrumentation and image processing techniques used for each survey. The differences between VGPS and CGPS images are small, < 6 K (rms) in channels where the mean HI brightness temperature in the field exceeds 80 K. A similar degree of consistency is found between the VGPS and SGPS. The agreement we find between arcminute resolution surveys of the Galactic plane is a crucial step towards combining these surveys into a single uniform dataset which covers 90% of the Galactic disk: the International Galactic Plane Survey (IGPS). The VGPS data will be made available on the World Wide Web through the Canadian Astronomy Data Centre (CADC).

The VLA Galactic Plane Survey (VGPS) of the first Galactic quadrant was searched for HI emission with velocities well above the maximum velocity allowed by Galactic rotation. A sample of 17 small fast-moving clouds was identified. The distribution of the ensemble of clouds in longitude and velocity indicates that the clouds are part of the Galactic disk, despite their large forbidden velocity. The median angular diameter of the clouds detected in the VGPS is 3.4 arcminutes. These clouds would not be noticed in previous low resolution surveys because of strong beam dilution. Assuming each cloud is located at the tangent point, a median cloud has a diameter of 10 pc, HI mass of 60 M_sun, and a velocity more than 25 km/s beyond the local terminal velocity derived from 12CO observations. Three clouds in the sample have a velocity between 50 and 60 km/s in excess of the local terminal velocity. The longitude distribution of the sample peaks near l = 30 degrees, while the latitude distribution of the clouds is nearly flat. The observed longitude and latitude distributions are compared with simulated distributions taking into account the selection criteria of the cloud search. It is found that the number of clouds declines with distance from the Galactic center, with an exponential scale length 2.8 - 8 kpc at the 99% confidence level. We find a lower limit to the scale height of the clouds of 180 pc (HWHM), but the true value is likely significantly higher.

We discuss the interpretation of the dust SED from the mid-infrared to the millimeter range of NGC 1569. The model developed by D\'esert et al. (1990) including three dust components (Polyaromatic Hydrocarbons, Very Small Grains and big grains) can explain the data using a realistic interstellar radiation field and adopting an enhanced abundance of VSGs. A simple three-temperature model is also able to reproduce the data but requires a very low dust temperature which is considered to be unlikely in this low-metallicity starburst galaxy. The high abundance of Very Small Grains might be due to large grain destruction in supernova shocks. This possibility is supported by ISO data showing that the emission at 14.3 $\mu$m, tracing VSGs, is enhanced with respect to the emission at 6.7 $\mu$m and 850 $\mu$m in regions of high star formation.

We present 21-cm HI line and 408 MHz and 1.4 GHz continuum observations of the dwarf galaxy UGCA 86 with the DRAO synthesis telescope. UGCA 86 is detected in the continuum at 408 MHz (S_{408} = 120 +/- 30 mJy) and 1.4 GHz (S_{1400} = 79 +/- 3 mJy). The HI structure of UGCA 86 is complex, with two separate components: a rotating disk and a highly elongated spur that is kinematically disjunct from the disk. The HI disk is centered on the optical galaxy with similar axial ratio and orientation of the major axis. An area of the disk with a peculiar velocity of 25 km/s relative to the regular rotation of the disk is found on the southern side, where most of the star formation is concentrated. The spur is seen along the minor axis of UGCA 86 and overlaps in part with the disk. Towards the optical center of UGCA 86, the velocity difference between the spur and the disk is 40 km/s, about one third of the rotation velocity of the HI disk at 6 kpc from the center. This implies a large radial component of the orbital velocity of the spur, and therefore a significantly non-circular orbit. The velocity dispersion of the disk is 8.8 km/s, whereas the velocity dispersion of the spur varies from 10 km/s to 30 km/s. A possible tidal origin of the spur is considered in view of the proximity of the large Scd galaxy IC 342. However, evidence that the spur is located far outside the plane of the HI disk, and the absence of evidence for a warp in the outer HI disk, poses a problem for the interpretation of the spur as a tidal tail induced by IC 342. Detailed modelling of the IC 342/UGCA 86 system will be required before a tidal origin of the spur can be dismissed conclusively. The possibility that the spur is part of the nascent cloud of UGCA 86 or the remains of a small dwarf galaxy is presented as an alternative interpretation.

GSH23.0-0.7+117 is a well-defined neutral hydrogen shell discovered in the VLA Galactic Plane Survey (VGPS). Only the blueshifted side of the shell was detected. The expansion velocity and systemic velocity were determined through the systematic behavior of the HI emission with velocity. The center of the shell is at (l,b,v)=(23.05,-0.77,+117 km/s). The angular radius of the shell is 6.8', or 15 pc at a distance of 7.8 kpc. The HI mass divided by the volume of the half-shell implies an average density n_H = 11 +/- 4 cm^{-3} for the medium in which the shell expanded. The estimated age of GSH23.0-0.7+117 is 1 Myr, with an upper limit of 2 Myr. The modest expansion energy of 2 * 10^{48} erg can be provided by the stellar wind of a single O4 to O8 star over the age of the shell. The 3 sigma upper limit to the 1.4 GHz continuum flux density (S_{1.4} < 248 mJy) is used to derive an upper limit to the Lyman continuum luminosity generated inside the shell. This upper limit implies a maximum of one O9 star (O8 to O9.5 taking into account the error in the distance) inside the HI shell, unless most of the incident ionizing flux leaks through the HI shell. To allow this, the shell should be fragmented on scales smaller than the beam (2.3 pc). If the stellar wind bubble is not adiabatic, or the bubble has burst (as suggested by the HI channel maps), agreement between the energy and ionization requirements is even less likely. The limit set by the non-detection in the continuum provides a significant challenge for the interpretation of GSH23.0-0.7+117 as a stellar wind bubble. A similar analysis may be applicable to other Galactic HI shells that have not been detected in the continuum.

The VLA Galactic Plane Survey (VGPS) consists of measurements of the 21cm HI line and 1420 MHz continuum at 1' angular resolution over much of the first quadrant of Galactic longitude within a few degrees of the Galactic plane. In combination with similar surveys of the fourth longitude quadrant and the outer Galaxy made with other instruments, the VGPS will provide comprehensive data on the interstellar medium in most of the regions of active star formation in the Galaxy. This paper describes the parameters of the VGPS, reports on its status, and discusses some early scientific results from the survey.

We present WSRT observations of high sensitivity and resolution of the neutral hydrogen in the starburst dwarf galaxy NGC 1569. Assuming a distance of 2.2 Mpc, we find a total HI mass of 1.3 * 10^8 M_sun to be distributed in the form of a dense, clumpy ridge surrounded by more extended diffuse HI containing a few additional discrete features, such as a Western HI Arm and an HI bridge reaching out to a small counterrotating companion cloud. About 10% by mass of all HI in NGC 1569 is at unusually high velocities. Some of this HI may be associated with the mass outflow evident from H-alpha measurements, but some may also be associated with NGC 1569's HI companion and intervening HI bridge, in which case, infall rather than outflow might be the cause of the discrepant velocities. No indication of a large bubble structure was found in position-velocity maps of the high-velocity HI. The galaxy as a whole is in modest overall rotation, but the HI gas lacks any sign of rotation within 60'' (0.6 kpc) from the center, i.e. over most of the optical galaxy. Here, turbulent motions resulting from the starburst appear to dominate over rotation. In the outer disk, the rotational velocities reach a maximum of 35 \pm 6 km/s, but turbulent motion remains significant. Thus, starburst effects are still noticeable in the outer HI disk, although they are no longer dominant beyond 0.6 kpc. Even excluding the most extreme high-velocity HI clouds, NGC 1569 still has an unusually high mean HI velocity dispersion of sigma_v=21.3 km/s, more than double that of other dwarf galaxies.

This paper is the first in a series presenting a sample of 29 late-type dwarf galaxies, observed with the Westerbork Synthesis Radio Telescope (WSRT) in the 21-cm line of neutral atomic hydrogen (HI). The sample itself, the HI content of and the HI distribution in the sample galaxies are briefly discussed. Four sample galaxies were also detected in the continuum.

This paper is the second in a series presenting a sample of 29 late-type dwarf galaxies observed with the Westerbork Synthesis Radio Telescope in the 21-cm line of neutral hydrogen (HI). Here we present rotation curves, maps of the velocity field and maps of the velocity dispersion across the sample galaxies.

We present new data of the dwarf galaxy NGC 1569 at 450 $\mu$m, 850 $\mu$m and 1200$\mu$m taken with SCUBA at the JCMT and the bolometer array at the IRAM 30m telescope. After including data from IRAS at 12, 25, 60 and 100 $\mu$m, we have successfully fitted the dust grain population model of D\'esert et al. (1990) to the observed midinfrared-to-millimeter spectrum. The fit requires a combination of both large and very small grains exposed to a strong radiation field as well as an enhancement of the number of very small grains relative to the number of large grains. We interpret this as the consequence of large grain destruction due to shocks in the turbulent interstellar medium of NGC 1569. The contribution of polyaromatic hydrocarbons (PAH's) is found to be negligible. Comparison of the dust emission maps with an HI map of similar resolution shows that both dust and molecular gas distributions peak close to the radio continuum maximum and at a minimum in the HI distribution. From a comparison of these three maps and assuming that the gas-to-dust mass ratio is the same everywhere, we estimate the ratio of molecular hydrogen column density to integrated CO intensity to be about 25 -- 30 times the local Galactic value. The gas-to-dust ratio is 1500 -- 2900, about an order of magnitude higher than in the Solar Neighbourhood.

The properties of the Galactic HI shell GSH 138-01-94 are derived from data of the Canadian Galactic Plane Survey. The basic parameters of GSH 138-01-94 were determined by fitting the expansion of a thin shell to the expansion velocity field on the sky. The kinematic distance is 16.6 kpc for v_LSR=-94.2 $\pm$ 0.5 km/s. The radius is 180 $\pm$ 10 pc, the expansion velocity v_exp = 11.8 $\pm$ 0.9 km/s, and the mass 2 * 10^5 M_sun. No radio continuum counterpart of the shell was detected at 21 cm or at 74 cm. Absorption of a background continuum source constrains the spin temperature of \HI in the shell to T_s=230^{367}_{173} K. The expansion age of GSH 138-01-94 is 4.3 Myr. These observables are in excellent agreement with predictions from hydrodynamic models for a supernova remnant in a low-density low-metallicity environment such as the outer Galaxy. GSH 138-01-94 is then the largest and the oldest supernova remnant known. It provides direct evidence for the release of mechanical energy in the interstellar medium by stars in the outer galaxy. It is argued that such old supernova remnants be found in low-density, low-metallicity environments such as the outer Galaxy, dwarf galaxies and low surface brightness galaxies.

Deep SCUBA observations of NGC 1275 at 450 micron and 850 micron along with the application of deconvolution algorithms have permitted us to separate the strong core emission in this galaxy from the fainter extended emission around it. The core has a steep spectral index and is likely due primarily to the AGN. The faint emission has a positive spectral index and is clearly due to extended dust in a patchy distribution out to a radius of $\sim$ 20 kpc from the nucleus. These observations have now revealed that a large quantity of dust, $\sim$ 6 $\times$ 10$^7$ $M_\odot$, 2 orders of magnitude larger than that inferred from previous optical absorption measurements, exists in this galaxy. We estimate the temperature of this dust to be $\sim$ 20 K (using an emissivity index of $\beta$ = 1.3) and the gas/dust ratio to be 360. These values are typical of spiral galaxies. The dust emission correlates spatially with the hot X-ray emitting gas which may be due to collisional heating of broadly distributed dust by electrons. Since the destruction timescale is short, the dust cannot be replenished by stellar mass loss and must be externally supplied, either via the infalling galaxy or the cooling flow itself.

We have studied the relation between HI linewidth and optical luminosity in two samples of dwarf galaxies, with special attention to the accuracy of the inclination correction. The first sample consists of rotationally supported galaxies, and has DDO 154 as a prototype. They have HI linewidths significantly larger than predicted by an extrapolation of the Tully-Fisher relation. The deviation from the extrapolated Tully-Fisher relation correlates with the distance-independent mass-to-luminosity ratio M_HI/L_B. The baryonic correction for these gas rich galaxies is too small to explain the deviation from the Tully-Fisher relation, unless a small mass-to-light ratio (M/L_B)_* <= 0.25 for the stellar population is assumed. The second sample contains dwarf galaxies with comparable rotational and random velocities. This sample, having Leo A as a prototype, is consistent with an extrapolation of the classical Tully-Fisher relation to lower luminosities. These dwarf galaxies are not particularly rich in gas, but their mass is dominated by baryonic matter within their HI radius. We suggest that any evolution of the deviating galaxies towards the Tully-Fisher relation will combine a decrease in linewidth by radial transport of HI with an increase in luminosity due to a larger accumulated stellar mass.

High-sensitivity maps of the large-scale structure of atomic hydrogen in the starburst dwarf galaxy NGC 1569 show evidence for an HI cloud with a mass of 7*10**6 M_sun, at a projected distance of 5 kpc from the parent galaxy. This cloud may be a condensation in a low-column-density HI halo or a companion galaxy/HI-cloud. NGC 1569 and its companion are connected by a low surface brightness HI bridge. At the edge of NGC1569, the HI bridge coincides with H_alpha arcs, also detected in soft X-rays.