# Diffuse polarized emission in the LOFAR Two-meter Sky Survey

C. van Eck, M. Haverkorn, M. Alves, R. Beck, P. Best, E. Carretti, K. Chyży, T. Ensslin, J. Farnes, K. Ferrière, et al.

► To cite this version:

C. van Eck, M. Haverkorn, M. Alves, R. Beck, P. Best, et al. Diffuse polarized emission in the LOFAR Two-meter Sky Survey. Astronomy and Astrophysics - A&A, 2019, 623, pp.A71. 10.1051/0004-6361/201834777. hal-02377436

HAL Id: hal-02377436

https://hal.science/hal-02377436

Submitted on 1 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Diffuse polarized emission in the LOFAR Two-meter Sky Survey*

C. L. Van Eck<fnref n="1" />,<fnref n="2" />, M. Haverkorn<fnref n="2" />, M. I. R. Alves<fnref n="2" />, R. Beck<fnref n="3" />, P. Best<fnref n="4" />, E. Carretti<fnref n="5" />, K. T. Chyży<fnref n="3" />,<fnref n="6" />, T. Enßlin<fnref n="7" />,<fnref n="8" />, J. S. Farnes<fnref n="2" />,<fnref n="9" />, K. Ferrière<fnref n="10" />, G. Heald<fnref n="11" />, M. Iacobelli<fnref n="12" />, V. Jelić<fnref n="13" />, W. Reich<fnref n="3" />, H. J. A. Röttgering<fnref n="14" />, and D. H. F. M. Schnitzeler<fnref n="3" />

Received 4 December 2018 / Accepted 25 January 2019

## ABSTRACT

Faraday tomography allows us to map diffuse polarized synchrotron emission from our Galaxy and use it to interpret the magnetic field in the interstellar medium (ISM). We have applied Faraday tomography to 60 observations from the LOFAR Two-meter Sky Survey (LOTSS) and produced a Faraday depth cube mosaic covering 568 square degrees at high Galactic latitudes, at 4'3 angular resolution and 1 rad m⁻<fnref n="2" /> Faraday depth resolution, with a typical noise level of 50–100 µJy per point spread function (PSF) per rotation measure spread function (RMSF; 40–80 mK RMSF⁻<fnref n="1" />). While parts of the images are strongly affected by instrumental polarization, we observed diffuse polarized emission throughout most of the field, with typical brightness between 1 and 6 K RMSF⁻<fnref n="1" />, and Faraday depths between −7 and +25 rad m⁻<fnref n="2" />. We observed many new polarization features, some up to 15° in length. These include two regions with very uniformly structured, linear gradients in the Faraday depth; we measured the steepness of these gradients as 2.6 and 13 rad m⁻<fnref n="2" /> deg⁻<fnref n="1" />. We also observed a relationship between one of the gradients and an H I filament in the local ISM. Other ISM tracers were also checked for correlations with our polarization data and none were found, but very little signal was seen in most tracers in this region. We conclude that the LOTSS data are very well suited for Faraday tomography, and that a full-scale survey with all the LOTSS data has the potential to reveal many new Galactic polarization features and map out diffuse Faraday depth structure across the entire northern hemisphere.

**Key words.** ISM: magnetic fields – polarization

## 1. Introduction

Magnetic fields are present throughout interstellar space and play an important role in many aspects of the interstellar medium (ISM), such as cloud collapse during star formation (van Loo et al. 2012), energy transports and cascades in magnetohydrodynamic turbulence (Beresnyak & Lazarian 2015), and pressure balance between different gas phases (Boulares & Cox 1990).

Interstellar magnetic fields can be measured using radio polarization through two processes: synchrotron emission and Faraday rotation. Synchrotron emission is produced throughout interstellar space by cosmic-ray electrons as they are accelerated by interstellar magnetic fields, resulting in polarized radio emission. Faraday rotation occurs when polarized emission passes

through magnetized plasma (which fills most of the volume of the ISM), which causes a frequency-dependent rotation of the polarization angle. The change in the polarization angle (Δθ) is given by

$$ \Delta\theta = \lambda^2 \phi(d) = \lambda^2 \left[ 0.812 \text{ rad m}^{-2} \int_{\text{source}}^{\text{observer}} \left( \frac{n_e}{\text{cm}^{-3}} \right) \left( \frac{\mathbf{B}}{\mu\text{G}} \right) \cdot \left( \frac{\text{d}l}{\text{pc}} \right) \right], \quad (1) $$

where φ(d) is the Faraday depth which depends on the free electron density (n_e) and magnetic field (B) along the line of sight integrated from the emission source at distance d to the observer.

Polarized synchrotron emission is produced throughout the Galaxy, and then undergoes Faraday rotation as it propagates through the ISM; therefore along any line of sight we see the superposition of the emission at all distances and with correspondingly different Faraday depths. The rotation measure (RM)

<footnote n="1">¹ Dunlap Institute for Astronomy and Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON M5S 3H4, Canada e-mail: cameron.van.eck@dunlap.utoronto.ca</footnote>

<footnote n="2">² Department of Astrophysics/IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands</footnote>

<footnote n="3">³ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany</footnote>

<footnote n="4">⁴ SUPA, Institute for Astronomy, Royal Observatory, Blackford Hill, Edinburgh, EH9 3HJ, UK</footnote>

<footnote n="5">⁵ INAF – Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047 Selargius (CA), Italy</footnote>

<footnote n="6">⁶ Astronomical Observatory, Jagiellonian University, ul. Orla 171, 30-244 Kraków, Poland</footnote>

<footnote n="7">⁷ Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85748 Garching, Germany</footnote>

<footnote n="8">⁸ Ludwig-Maximilians-Universität München, Geschwister-Scholl-Platz 1, 80539 München, Germany</footnote>

<footnote n="9">⁹ Oxford e-Research Centre (OeRC), Keble Road, Oxford OX1 3QG, UK</footnote>

<footnote n="10">¹⁰ IRAP, Université de Toulouse, CNRS, 9 avenue du Colonel Roche, BP 44346, 31028 Toulouse, Cedex 4, France</footnote>

<footnote n="11">¹¹ CSIRO Astronomy and Space Science, PO Box 1130, Bentley, WA 6102, Australia</footnote>

<footnote n="12">¹² ASTRON, the Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA Dwingeloo, The Netherlands</footnote>

<footnote n="13">¹³ Ruder Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia</footnote>

<footnote n="14">¹⁴ Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands</footnote>

<footnote>* The data shown in Figs. 2–7 are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/623/A71</footnote>
synthesis technique (Brentjens & de Bruyn 2005) can be used to transform the observed wavelength-dependent polarization into the distribution of polarized emission as a function of Faraday depth. When applied to 3D image-frequency data this is called Faraday tomography, which produces Faraday depth cubes that map out the diffuse polarized emission as a function of position on the sky and Faraday depth. Such observations can be used to constrain magnetic fields in the ISM and study their properties (e.g., Schnitzeler et al. 2007; Jelić et al. 2015; Lenc et al. 2016; Van Eck et al. 2017).

The resolution in Faraday depth of Faraday tomography depends on the range of wavelength-squared sampled by the observations, so low-frequency radio telescopes (which can produce very large wavelength-squared coverage) are capable of achieving much finer Faraday depth resolution than higher-frequency instruments. The newest generation of very low-frequency radio telescopes such as the Low Frequency Array (LOFAR, van Haarlem et al. 2013) and the Murchison Wide-field Array (MWA, Tingay et al. 2013), which operate in the 100–200 MHz range, are capable of reaching Faraday depth resolutions around 1 rad m⁻<fnref n="2" />, approximately two orders of magnitude smaller than the resolution that can be achieved at 1.4 GHz or higher frequencies.

In recent years, there have been many studies using Faraday tomography to study diffuse polarization, each with different observations balancing tradeoffs between field-of-view, angular resolution, and Faraday depth resolution: Jelić et al. (2014, 2015), and Van Eck et al. (2017) used LOFAR observations, Bernardi et al. (2013) and Lenc et al. (2016) used MWA observations, Hill et al. (2017) used data from the Global Magneto-Ionic Medium Survey (GMIMS, Wolleben et al. 2008), while Iacobelli et al. (2013) used observations from the Westerbork Synthesis Radio Telescope (WSRT). The observations with high angular resolution tend to have smaller fields of view, and vice versa, making it difficult to observe the full size of large (several degree) polarization features (such as the filament in Jelić et al. 2015 that extends outside the field of view) while still minimizing the effects of beam depolarization (which reduces the signal present in low-resolution observations).

In this paper we report on the polarization processing of 60 LOFAR observations from the LOFAR Two-meter Sky Survey (LOTSS, Shimwell et al. 2017) and present a Faraday depth mosaic covering an area of 568 square degrees. This is the first work to combine multiple LOFAR observations into a single Faraday depth cube. This allows us to take advantage of the high angular and Faraday depth resolution of LOFAR while still looking at large-scale features in the diffuse polarization. In Sect. 2 we present our data reduction and the production of Faraday depth cubes for each observation. In Sect. 3 we present a Faraday depth cube mosaic produced from these observations, and highlight the diffuse polarization features seen. In Sect. 4 we compare our mosaic against other tracers for components of the ISM, and in Sect. 5 we discuss possible origins for the polarized features we observe. In Sect. 6 we consider the depolarization effects of Faraday rotation gradients and the implications for our observations and others. In Sect. 7 we summarize our results.

## 2. Data processing

We analyzed 63 calibrated visibility datasets generated with observations from the LOFAR Two-meter Sky Survey; full details of the observational parameters and calibration methods can be found in Shimwell et al. (2017). The observations

are from the LOTSS test region, right ascension from 10h30m to 15h30m and declination from 45° to 57°, which covers the HETDEX Spring field (Hill et al. 2008), a region near the Galactic north pole. Each observation has a nominal duration of 8 h and covers the frequency range 120–168 MHz with 488 channels. We received the observations after direction-independent amplitude and phase calibration with the LOTSS pipeline. Figure 1 shows the coverage of these observations in equatorial and Galactic coordinates. The polarization calibration and imaging closely follows that of Van Eck et al. (2017) and is described below.

Polarization calibration in the form of a correction for ionospheric Faraday rotation was performed prior to imaging. This was done using the RMextract package<fnref n="1" /> written by Maaijke Mevius, combined with maps of the ionospheric total electron content from the Center for Orbit Determination in Europe (CODE)<fnref n="2" />. This package produced predictions for the ionospheric Faraday rotation that were used by the Black Board Selfcal software (BBS, Pandey et al. 2009) to derotate the polarization of the visibilities and thereby remove the ionospheric Faraday rotation. The estimated error in Faraday depth correction is approximately 0.1 rad m⁻<fnref n="2" /> (Van Eck et al. 2018). No correction for polarization leakage was performed, as it was deemed too computationally expensive using the currently available methods. The effects of polarization leakage on the results are discussed in Sect. 3.1.

Polarization imaging was performed independently for each channel using AWImageR (Tasse et al. 2013). A baseline length upper limit of 800λ was applied, as this made the synthesized beam (hereafter, point spread function or PSF) more uniform across the band, giving a FWHM of 4.3. The short baselines between the two halves of each high-band antenna (HBA) station (i.e., substations CS(X)HBA0 and CS(X)HBA1, iterating (X) over all HBA stations; see van Haarlem et al. 2013, for a description of the substation layout and naming) were removed, as these were known to often suffer from significant mutual interference. A robust weighting of 1.0 was used for imaging, and no CLEANing was performed, as the signal to noise ratio in individual channels was expected to be too low for CLEAN to be effective. By default AWimageR produces images both with and without correction for the LOFAR primary beam; the images without correction were used for the quality control step described next, while the images with correction were used for producing Faraday depth cubes.

We found that a small fraction of the images were strongly affected by interference or calibration problems, which usually manifested as very strong patterns throughout the image. To identify these channels in an automated way, we calculated the standard deviation over all pixels in each of the images without primary beam correction (separately for Stokes Q and U). The channels affected by these problems stood out as having abnormally large standard deviations, but variations in the background noise across the band and between different observations made it difficult to assign a single threshold value for classifying images as bad. We chose to calculate, for each channel, the median of the standard deviations of the images within 50 channels above or below in frequency; if the standard deviation (for either Q or U separately) in a channel was more than 1.5 times the median standard deviation of these neighbouring channels then that channel was flagged as bad. These bad channels were removed before RM synthesis. During this process, three observations were found to

<footnote n="1">¹ https://github.com/maaijke/RMextract/</footnote>

<footnote n="2">² http://aiuws.unibe.ch/ionosphere/</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Locations of the LOTSS HETDEX observations, in both Equatorial (dotted lines) and Galactic (dashed lines) coordinates. Each gray circle covers the area inside the primary beam FWHM (3.4° diameter) for one observation, at the highest frequency (168 MHz). This image, and all that follow, are in an orthographic projection. The gaps are caused by the three fields that were removed for being very strongly affected by instrumental leakage.</figcaption></figure>

have much higher standard deviations across the full bandwidth than the others, which was due to the presence of extremely bright radio sources; these observations were removed and not included in the following steps. After this step, the typical noise level in a single channel was 2–5 mJy PSF⁻¹.

RM synthesis was performed using pyrmsynth<fnref n="3" />. Channel weights were calculated using the inverse square of the standard deviation of each channel that was calculated in the previous step. For most observations, the weights were fairly constant across the band, so this weighting choice was not expected to produce any significant changes from uniform weighting. Weights were applied independently per channel; no adjustments for sampling density in λ² were applied. From the frequency coverage and weighting used, the resulting rotation measure spread function (RMSF) had a typical full-width at half-max (FWHM) of 1.2 rad m⁻². From Eq. (61) of Brentjens & de Bruyn (2005), the theoretical FWHM (for uniform weighting) is 1.15 rad m⁻², so these values are consistent and show that the natural weighting has not significantly degraded the Faraday depth resolution. From Eqs. (62) and (63) of Brentjens & de Bruyn (2005), our observations are not sensitive to extended structures in Faraday depth wider than about 1.0 rad m⁻², and we have limited sensitivity to Faraday depths larger in absolute value than 170 rad m⁻². Faraday depth measurements beyond this range may also be unreliable (Schnitzeler & Lee 2015). As is typical for LOFAR observations, we are unable to resolve any features with Faraday thickness greater than the resolution, therefore we are only able to pick up unresolved features (or sharp edges of Faraday thick features, as described in Van Eck et al. 2017). No spectral index correction was included in the RM synthesis step; this may produce minor errors in the polarized intensity and small secondary peaks in the RM spectrum, but these were not expected to affect our analysis. Each RM cube covered Faraday depths from -100 to +100 rad m⁻² in steps of 0.25 rad m⁻². The RM-CLEAN

algorithm (Heald et al. 2009) was applied to each Faraday depth cube, to a depth of 2 mJy PSF⁻¹ RMSF⁻¹.

In order to combine the Faraday depth cubes from the different observations, we estimated the position dependent noise in the cubes. We did this on a per-pixel basis by taking the distribution of polarized intensity values taken from the regions of the spectrum that were expected to be dominated by noise (|φ| > 20 rad m⁻²) and fitting a Rayleigh distribution (which is the expected distribution for polarized intensity in the absence of signal, Macquart et al. 2012; Hales et al. 2012). The resulting Rayleigh σ parameter (which is equivalent to the Gaussian σ of the underlying Stokes Q and U distributions) was taken to be the noise in that pixel.

To merge the separate Faraday depth cubes we averaged the overlapping regions using inverse-variance weighting to produce a final mosaic covering the entire HETDEX region. Alternative weighting schemes, where each pixel used only the observation with the lowest noise or the observation with the nearest pointing center, were also tried but were found to produce significant artifacts in the overlap regions of adjacent pointings. However, these were useful for assessing the quality of the data and the effectiveness of the mosaicing; we observed that many features in the cubes, such as depolarization canals, could be traced continuously across the boundaries between different observations.

### 3. The HETDEX mosaic

The combined Faraday depth cube resulting from the processing described previously covers a total area of sky of 568 square degrees, with a typical noise level between 50 and 100 µJy PSF⁻¹ RMSF⁻¹, with higher noise at the edges and near very bright Stokes *I* sources. Figures 2–7 show selected slices from the cube. These images give the polarized intensity in mJy PSF⁻¹ RMSF⁻¹; the conversion to brightness temperature units is 0.8 K (mJy PSF⁻¹)⁻¹. Figures 8–11 show collapsed versions of

<footnote n="3">³ https://github.com/mrbell/pyrmsynth</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between −7 and −4 rad m⁻². **Top panel:** typical quiescent slice, with strong artifacts around 3C295. **Lower two panels:** diffuse polarized emission appearing in the top right.</figcaption></figure>

sections of the cube where, per pixel, the highest polarized intensity was located and the polarized intensity and Faraday depth of the peak were used to determine the brightness and color, respectively. Figure 8 shows the full mosaic, while Figs. 9–11 focus on the specific diffuse polarized features described in Sects. 3.2.1–3.2.3. Figure 12 shows some typical Faraday depth spectra from each region. In the cube we see many diffuse polarized emission features, polarized point sources, and instrumental polarization leakage, which are discussed in detail below.

While this mosaic covers a large area of sky and allows large features to be followed across several pointings, it is not sensitive to very large polarized features. The shortest baselines present in the observations are approximately 15λ in (projected) length, corresponding to an angular scale of about 3.8°; uniform polarization features larger than this scale will be missed and this should be considered when interpreting the mosaic. The distribution of baselines is shown in Fig. 13. However, at such low frequencies this may not be as significant a
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between −3 and −1 rad m⁻². The feature in the top right travels south-westward, and another feature appears in the top center and spreads outwards. The instrumental leakage from bright Stokes *I* sources also appear strongly at −2 and −1 rad m⁻².</figcaption></figure>

problem as at higher frequencies, due to Faraday rotation gradients driving emission to smaller angular scales (Schnitzeler et al. 2009). We can observe polarization features with a change in polarization angle greater than $\pi$ radians over a 3.8° angular scale; if we assume that all of this change is due to Faraday rotation fluctuations, we can calculate the amplitude of Faraday

rotation variation necessary to cause this polarization change. For a typical frequency of 150 MHz, this corresponds to a Faraday depth change of 0.8 rad m⁻² per 3.8° or 0.21 rad m⁻² deg⁻¹. Polarization features with greater variation in Faraday depth than this value should be detected by our observations.
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between 0 and 2 rad m⁻². The top-right feature continues to travel south-westward, becoming very filamentary-looking, while the top-center feature continues to expand outwards. Two more diffuse features appear in the lower right corner and lower left-of-center.</figcaption></figure>

**3.1. Instrumental polarization leakage**

The instrumental polarization leakage for LOFAR, which causes emission from Stokes *I* to incorrectly appear in the other Stokes parameters, has the convenient property of being effectively

independent of frequency. The result of this is that after RM synthesis, the leakage should appear only at Faraday depth 0 rad m⁻², which would allow astrophysical emission at other Faraday depths to be identified. However, the ionospheric Faraday rotation correction applies a frequency-dependent
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between 3 and 5 rad m⁻². The top right feature now appears as a very long filamentary structure spanning the height of the image, before fading away. The lower right feature also fades away at higher Faraday depths. The top center and lower left-of-center features now appear as a series of patches running south-westward just left of center. A faint, very straight filamentary feature appears in the bottom mid-left, slowly traveling upwards.</figcaption></figure>

polarization angle rotation to the data, which shifts the leak-
age to the opposite Faraday depth as the correction (e.g., if the
ionospheric Faraday rotation is +2 rad m⁻², the leakage will
appear at −2 rad m⁻² after the correction). Since the ionospheric

correction varies with time, this also has the result of caus-
ing partial depolarization of the leakage, as the post-correction
instrumental leakage will have different polarization angles at
different times.
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between 7 and 12 rad m⁻². The intensity scale has been adjusted as the emission at these depths is much fainter. The top right and left-of-center features continue to fade out, while the faint very straight filament continues to travel upwards with new, thin features appearing above it.</figcaption></figure>

The net effect of the instrumental polarization leakage is to cause all Stokes *I* sources to appear in the Faraday depth cubes between approximately −2 and 0 rad m⁻² (as the ionospheric correction was typically between 0 and 2 rad m⁻²). Since the different observations had different ionospheric

conditions, the leakage sources appear at slightly different Faraday depths at different locations in the mosaic, and some regions also appear to have stronger or weaker leakage depending on how much depolarization the time-variability produced.
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Selected slices from the Faraday depth cube mosaic, showing the polarized intensity between 15 and 21 rad m⁻². The faint, straight filamentary features slowly converge and disappear. **Top panel:** a very bright polarized source appears to the left of 3C295. Some emission appears to be present at the bottom of the image, just left of center, but it is not clear if this is real or caused by enhanced noise at the edge of the mosaic.</figcaption></figure>

In addition there is significant off-source polarization leakage, which is the dominant cause of the strong artifacts appearing widely distributed in Faraday depth around very bright sources such as 3C295 (α = 14h 11m 20.59s, δ = +52° 12' 9.6"', 84 Jy in Stokes *I* at 178 MHz; Spinrad et al. 1985). The spurious

polarization caused by the polarization leakage does not appear just at the location of the source, but is convolved with the PSF, just as the Stokes *I* and real polarized emission is. Since the locations and intensities of the sidelobes of the PSF vary with frequency, pixels near sources can pass in and out of
<figure><img src="image_9.png" /></figure>

**Fig. 8.** Collapsed view of the full mosaic, where the polarized intensity and Faraday depth of the brightest feature in the Faraday depth spectrum were used to determine the brightness and color, respectively. Instrumental leakage from bright point sources, particularly 3C295, dominate many areas with features at Faraday depths around -2 rad m⁻², but diffuse emission can be seen through most of the mosaic.

sidelobes as a function of frequency. The result is a very complicated frequency-dependent spurious polarization signal that, when RM synthesis is applied, results in structure broadly distributed in Faraday depth. Thus, while the on-source polarization leakage is usually confined to a single Faraday depth (determined by the ionosphere), the off-source leakage is broadly distributed in Faraday depth and covers larger areas around brighter sources.

In the HETDEX region mosaic, the off-source leakage of the two brightest sources occupies a large portion of the field: 3C295 (84 Jy) is in the mid-east (left) part of the mosaic, and 3C280 ($α = 12^h56^m59^s, δ = +47°18′48″, 24 Jy at 178 MHz; Spinrad et al. 1985) is in the bottom center of the mosaic. As a result the region around 3C295 is unusable for polarization analysis with these data.

## 3.2. Diffuse emission

Diffuse polarized emission is seen through much of the mosaic, with polarized intensity levels ranging from 6–8 mJy PSF⁻¹ RMSF⁻¹ for the brighter features to 1 mJy PSF⁻¹ RMSF⁻¹ for the faintest identified features. The presence of polarization leakage has made quantitative analysis of the diffuse emission difficult, as it is not clear how to separate the diffuse emission from the off-source leakage. We have divided the diffuse emission into three regions, based on their morphologies and positions, and discuss each in turn below.

### 3.2.1. Bright northwest gradient and southwest patchy emission

This region, shown in Fig. 9, occupies the field west (right) of 12^h right ascension, and is dominated by two bright features that appear connected. The most striking feature is a gradient in Faraday depth that occurs in the northwest part of the field, from approximately 11^h to 12^h in right ascension and +53° to +57°30′ in declination<fnref n="4" />. This emission feature appears at Faraday depths as low as −7 rad m⁻² (Fig. 2, top panel) and can be traced as a single continuous sheet in the cube up to +7 rad m⁻² (Fig. 6,

top panel), with typical brightness of 6–8 mJy PSF⁻¹ RMSF⁻¹. Between −4 rad m⁻² and +3 rad m⁻² (Fig. 2, bottom panel, to Fig. 5, top panel), this feature has the appearance of a “traveling filament”, caused by a very linear gradient in Faraday depth. We measured locations along the “filament” at −3.5 rad m⁻² and at +3 rad m⁻², and found that the separation between corresponding points was approximately 2.5°, resulting in a Faraday depth gradient of 2.6 rad m⁻² deg⁻¹.

At Faraday depths around −1 rad m⁻², additional emission emerges from many locations in the west-most third of the field and expands with increasing Faraday depth. While this emission looks patchy and disconnected in individual slices (particularly at 0, +1, and +5 rad m⁻²) and in Fig. 9, careful inspection of the full cube shows that the emission is a continuous sheet “wrinkled” in Faraday depth, with many local minima and maxima in Faraday depth surrounded by a complex web of emission at intermediate Faraday depths. This is best seen at +3 rad m⁻² (Fig. 5, top panel) where many “loops” of emission can be seen (surrounding the locations of local Faraday depth maxima or minima); this slice also best shows the extent of the emission, which extends the full height of the field.

Near the region occupied by the linear gradient, a second emission feature is also present at higher Faraday depths. It can be seen at +3 rad m⁻² at the same location occupied by the first feature at −5 rad m⁻² (as a result, it cannot be seen in Fig. 9, which only shows the brightest Faraday-depth peak found in each pixel). It is seen over a much smaller area than the first feature, but this may be due to the much lower brightness of this feature, 1–2 mJy PSF⁻¹ RMSF⁻¹. It can be seen at Faraday depths up to +14 rad m⁻², and while it has some morphological similarities to the first feature (gradients and local maxima at similar locations), it also has some significant differences (a notable absence of emission in certain regions that are bright in the first feature).

### 3.2.2. Central sheet

The region between 12^h and 14^h, shown in Fig. 10, is also dominated by a bright feature. This feature appears to be a continuous sheet of emission distributed in Faraday depth between −3.5 and +12 rad m⁻², with several local maxima and minima in Faraday depth at various locations and filamentary-looking

<footnote n="4">⁴ When referencing coordinates to locations in the images, it is important to recall that this is not a rectilinear projection; the projection of equatorial coordinates into our images is shown in Fig. 1.</footnote>
<figure><img src="image_10.png" /></figure>

**Fig. 9.** Collapsed view of the west region of the mosaic, where the polarized intensity and Faraday depth of the brightest feature in the Faraday depth spectrum were used to determine the brightness and color respectively. The features in the north west and north center, marked by the gray box and described in Sect. 3.2.1, can be clearly seen as sheets of polarized emission with Faraday rotation gradients, while the feature in the southwest shows much more patchy structure in polarized intensity. The instrumental leakage causes bright Stokes *I* sources to appear as small blue circular features. The three crosses mark the locations of the spectra shown in Fig. 12.

emission. The typical brightness is 3–4 mJy PSF⁻¹ RMSF⁻¹. The brightest region occurs at the northern edge of the field, where emission appears at −3.5 rad m⁻² and spreads outwards in all directions with increasing Faraday depth, with a typical brightness of 6–8 mJy PSF⁻¹ RMSF⁻¹. A second bright spot occurs to the southeast, around (13h15′,47°), with Faraday depths of 5–7 rad m⁻².

### 3.2.3. Southeast “filament”

The third region, shown in Fig. 11 is between 14^h^ and 15^h^ right ascension, and between +45° and +50° declination, and contains another continuous feature with a very linear Faraday depth gradient. The feature appears as a very long thin filament in Faraday depth slices from +1 to +11 rad m⁻², which in some slices can be seen to extend over 12° in length. As with the other features, this can be interpreted as a continuous sheet of polarized emission.

Around +9 rad m⁻², additional emission can be seen further north, and as Faraday depth increases these merge into the filaments toward several local maxima at Faraday depths between +19 and +25 rad m⁻². The typical brightness of this feature is 1–2 mJy PSF⁻¹ RMSF⁻¹. The change in position of the filament was measured between +3 and +11 rad m⁻², and was found to be 0.6°, indicating a Faraday depth gradient of 13 rad m⁻² deg⁻¹.

### 3.3. Point sources

In addition to diffuse Galactic polarized emission, polarized emission is also observed from some of the unresolved point sources. The identification, measurement, and analysis of these sources is done in a parallel paper (Van Eck et al. 2018), and will not be discussed here. The majority of sources are found at Faraday depths between +10 and +25 rad m⁻², higher than most of the diffuse emission. This implies that the observed diffuse
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** As Fig. 9, but for the central region of the mosaic described in Sect. 3.2.2. The brightest region in the top, right of center and marked by the gray box, shows a minimum in Faraday depth of –3 rad m⁻², with increasing Faraday depth in the surrounding regions. The three crosses mark the locations of the spectra shown in Fig. 12.</figcaption></figure>

emission is from only a portion of the line of sight through the Galaxy, with the more distant parts contributing additional Faraday rotation to the extragalactic sources.

**4. Comparison with other tracers**

Previously published studies of Faraday depth cubes have looked at relating the observed diffuse polarization to tracers of other ISM components. Zaroubi et al. (2015) related Faraday depth structure observed in one LOFAR field to high-frequency (353 GHz) polarization maps, which trace dust emission, from the Planck mission, and found that the observed filaments in Faraday depth followed the magnetic field orientation inferred from the Planck polarization. Van Eck et al. (2017) suggested that neutral clouds could act as sources of Faraday-thin polarized emission (which would not be strongly depolarized at low frequencies), and associated two observed polarization features with neutral clouds in the local ISM.

We investigated several ISM tracers for features corresponding to those we see in our Faraday depth cube. For each tracer,

we regridded the data onto the same projection as our field, for ease of comparison. These tracers and their sources are:

- $H\alpha$ (integrated over velocity), from Finkbeiner (2003), with a resolution of 1°;

- 408 MHz radio continuum, from Remazeilles et al. (2015); Haslam et al. (1982), with a resolution of 56';

- Planck thermal dust emission, with a resolution of 60', at a reference frequency of 545 GHz; all Planck maps taken from the Commander component separation in Planck Collaboration X (2016);

- Planck thermal dust polarization, with a resolution of 10', at a reference frequency of 353 GHz;

- Planck synchrotron emission, with a resolution of 60', at the reference frequency of 408 MHz;

- Planck synchrotron polarization, with a resolution of 40', at the reference frequency of 30 GHz;

- Planck CO(1–0), with a resolution of 60';

- Planck free–free emission, with a resolution of 60';

- integrated HI line emission from the Effelsberg-Bonn HI survey (EBHIS, Winkel et al. 2016), with a resolution of 11';
<figure><img src="image_12.png" /></figure>

**Fig. 11.** As Fig. 9, focusing on the region around the southeast “filament” described in Sect. 3.2.3 and marked by the gray box. The contours are the H I intensity integrated over the velocity range –46 to –40 km s⁻¹ and were selected to highlight the H I filament present in this region; the levels are 3, 6, 9, and 12 K km s⁻¹ transitioning from white to dark gray. The smooth gradient in the Faraday depth is located below the edge of the H I filament, while inside the filament the Faraday depth behavior is more complex. Leakage from the bright source 3C295 dominate the top of the image. The three crosses mark the locations of the spectra shown in Fig. 12.

– dust extinction, calculated from Green et al. (2015) as described below, with a resolution of 42';

– 1.4 GHz polarized intensity, from Wolleben et al. (2006), with a resolution of 36'.

The maps of each of these tracers are shown in Figs. 14 and 15. The Planck polarization maps were converted from Stokes Q and U to polarized intensity. The extinction maps were made by using the MWDUST package (Bovy et al. 2016) to get the optical extinction as a function of distance, taking the numerical derivative with respect to distance (which should serve as a proxy for dust density). This produced a 3D data cube of extinction per unit distance, which was then integrated over selected distance ranges to produce extinction maps for different distances. The distance ranges were chosen to be 0–75, 75–250, and 250–1000 pc. The boundary at 75 pc was chosen as the Green et al. (2015) extinction model is often poorly constrained for distances less than this; the boundary at 1000 pc was chosen as no significant extinction contribution was seen beyond this distance; the boundary at 250 pc was chosen as the extinction showed different spatial structure on both sides of this boundary (i.e., the morphological differences between the two were maximized by placing the boundary at 250 pc). The extinction maps were observed to have a large scatter between adjacent pixels, so they were smoothed using a Gaussian kernel with a FWHM of 42' ($\sigma = 18'$) to produce the maps shown in Fig. 15.

The synchrotron emission and polarization, 408 MHz radio continuum, $H\alpha$, free-free emission, and CO emission maps are

all largely featureless, so no comparisons with the low-frequency radio polarization were possible. The H I column density, thermal dust emission and polarization, and extinction between 75 and 250 pc all show similar structure, and the location of this structure does not appear to strongly correlate with the position of the low-frequency polarized emission.

However, there is a very interesting relationship between the southeast Faraday depth gradient feature and the H I filament in the southeast part of our field. Figure 11 shows this region of the field, comparing our Faraday depth cube with contours made by integrating the H I data over the velocity range of the filament (–46 to –40 km s⁻¹). The H I filament has a sharp boundary with minimal emission (at any velocity) south of the lowest contour. The linear gradient in Faraday depth occurs south of the H I filament, and at the location of the filament the Faraday depth structure is more complex, with several local maxima, and the polarized intensity is more patchy. This transition from linear gradient to more complicated structure is very closely aligned with the boundary of the filament, including the slight bend around $\alpha = 14^h30^m$, $\delta = +47^\circ20'$, which is strong evidence that this is not coincidental positioning. We discuss interpretations of this relationship in Sect. 5.3. We note that while the southeast gradient shows this clear alignment with an H I filament, this is not the case for the northwest gradient, which does also occur in a region of low H I emission but has no H I features in the area where the gradient transitions to a more complex morphology.
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Selected Faraday depth spectra from the three regions described in Sect. 3.2, covering points in the northwest gradient (top panel), the central sheet (middle panel), and southeast filament (bottom panel) regions. The locations corresponding to each spectrum are shown in Figs. 9–11.</figcaption></figure>

The comparison of our 150 MHz polarization data to the 1.4 GHz polarization data from Wolleben et al. (2006) is interesting, but complicated by a number of factors. First is that we are comparing a Faraday depth cube, which can separate out multiple Faraday depth components, to a frequency-averaged polarization map, which will show a combination of all Faraday depth components present. Second, and related to the first, is that our low-frequency data will not be sensitive to Faraday-thick components while the higher frequency data may include emission from such components (albeit possibly partially depolarized). Third, the 1.4 GHz data is at a much lower resolution, 36', so it will likely be more strongly affected by beam depolarization.

With those caveats, there are some interesting relationships between the 1.4 GHz polarized intensity and some of the features

in our Faraday depth cube. Several of the brighter regions in the 1.4 GHz data correspond to brighter regions in the Faraday depth cube: the brighter region in the northwest region, around (11$^{h}$30′, +57°), overlaps with the bright, low-Faraday depth region in the northwest gradient (the green region in Fig. 9); the bright region in the north-center, around (13$^{h}$, +57°), overlaps the brightest part of the central sheet (the bright green region at the top of Fig. 10); and the bright region slightly east and south of the center, around (13$^{h}$30′, +49°), overlaps with the second bright region in the central sheet (red region in the lower left of Fig. 10). It is particularly noteworthy that each of these regions corresponds to a maximum or a minimum in Faraday depth, and that the regions with intermediate, more spatially variable Faraday depths show weaker polarized intensity at 1.4 GHz. This is not
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Distribution of baseline lengths in a typical LOTSS observation (after applying an upper limit of 800λ), and corresponding angular scales. The shortest baselines set the largest angular scale at approximately 3.8°, while the resolution of 4.3 is set by the imposed limit in baseline length.</figcaption></figure>

a perfect correlation: the bright 1.4 GHz region in the southeast region of the field, around (14$^h$30′, +47°), overlaps with the steep Faraday depth gradient (the “filament” in Fig. 11).

One possible explanation is that this might be depolarization caused by the Faraday depth gradients. However, this is not supported by the results of Appendix A, as the gradient of 2.5 rad m$^{-2}$ deg$^{-1}$ would produce only a 4° change in polarization angle across the 1.4 GHz beam, so minimal depolarization by the gradient is expected. Based on this, we rule out this explanation.

Finding an explanation becomes more difficult when the polarized intensities are compared. We estimate the typical total-intensity spectral index in this region is -2.6 (using the 150 MHz data from Landecker & Wielebinski 1970 and the 1420 MHz data from Reich 1982), so under the simplest assumption of a constant polarized fraction (perfectly Faraday-thin emission) the expected polarized intensity at 150 MHz can be estimated. In the areas of faintest 1.4 GHz polarization, the polarized brightness temperature at 1.4 GHz is typically 20–40 mK, which would be equivalent to 7–14 K at 150 MHz. Correspondingly, the brighter areas with polarized brightness temperatures of 100–150 mK would be equivalent to 35–50 K. In the fainter regions, where the northwest gradient is, the typical polarized intensity we observe at 150 MHz is about 3–5 K; correspondingly in the brightest regions of the central sheet we observe maximal intensities of about 8 K.

These values imply we are only detecting roughly 40% of the polarized intensity in the faint regions and about 20% in the brighter regions. The two most straightforward explanations for this are that the “missing” polarized intensity is either in Faraday-thick components and correspondingly depolarized at low-frequencies, or exists on the very large angular scales not observed by LOFAR. However, it’s not clear why either of these effects would cause a higher degree of missing emission in regions of higher polarized intensity.

In retrospect, it may not be surprising that many of the ISM tracers show very little signal or structure. The HETDEX field was originally selected for high-redshift cosmology observations in the optical, and was chosen for minimal contamination from Galactic foregrounds. As a result, many of the tracers we

investigated were barely detected in this region. In this respect, the choice of initial region for LOTSS is somewhat unfortunate for Galactic foreground science, but it also suggests that future investigations using observations in other regions, particularly at lower Galactic latitudes, will probably have more visible ISM structure against which to compare the Faraday tomography observations.

# 5. Interpretation of diffuse emission

The lack of correlations between the polarized emission and the ISM tracers provides constraints on the source of the polarized emission. Below we consider separately the source of the Faraday thin polarized emission, and the Faraday rotation of that emission.

## 5.1. Origin of Faraday thin emission

Due to the frequency coverage of our observations, we are insensitive to polarized features with Faraday thickness greater than about 1.0 rad m$^{-2}$. This places strong constraints on the ISM structures that we can observe; Van Eck et al. (2017) discussed which ISM conditions could cause Faraday-thin features. Here we consider these conditions in relation to the polarization features were observed.

A localized enhancement in the magnetic field perpendicular to the line of sight, such as in a shock, could produce a Faraday-thin feature, but would also produce a synchrotron excess which should be seen in total intensity. An enhancement in the degree of order in the magnetic field could also produce such features. There are also no obvious sources of such enhancements, such as supernova remnants, in this field, so both of these possibilities seem unlikely.

Strong localized enhancements or diminishments in the free electron density could also produce Faraday thin features. Enhancements would be associated with objects like shocks or HII regions, neither of which are observed in this field. Typical locations of low free electron density would be neutral clouds (associated with the warm neutral or cold neutral ISM phases) or regions of the hot ionized ISM phase. The tracers of neutral material (H1, thermal dust, extinction) do not show any correlation to the observed polarized emission, and the local ISM models of Lallement et al. (2014) do not show any clouds within the nearest few 100 pc in this direction. The Local Bubble is a volume of hot ionized medium (HIM) surrounding the Sun, and may (depending on the path length in this direction) have internal Faraday rotation of less than 1 rad m$^{-2}$ which would result in a Faraday thin feature (Van Eck et al. 2017). However, this feature would not have any foreground to produce Faraday rotation, so the feature would appear at 0 ± 0.5 rad m$^{-2}$ and could not be the source of the emission features seen at other Faraday depths. More distant regions of HIM could be the source of that emission, provided they are not so large as to become Faraday thick.

Another possibility is a region where the component of the magnetic field parallel to the line of sight is very small. If this occurs as a result of a reversal in the sign of the parallel magnetic field this is called a Faraday caustic (Bell et al. 2011), but lines of sight where the parallel component is very small in some volume but has the same sign in front and behind this volume are also possible. Such regions are possible in any phase of the ISM. Bell et al. (2011) showed that Faraday caustics can form sheets of Faraday-thin emission, covering the area of sky with the sign-reversal. While Bell et al. (2011) show
<figure><img src="image_15.png" /><figcaption>**Fig. 14.** Selected ISM tracers, in the same coordinates and projection as the previous images. The H I filament in the southeast region of the field is well correlated, but slightly offset in position, with the steep linear gradient in Faraday depth described in Sect. 3.2.3; this is shown in more detail in Fig. 11.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Fig. 15.** Continuation of Fig. 14, with additional tracers. The DRAO 21-cm polarized intensity shows stronger emission in several regions where we observe bright features and maxima/minima in Faraday depth in our data. To guide the eye, the gray boxes from Figs. 9–11 have been overlaid.</figcaption></figure>
that Faraday caustics should have a distinctive one-sided tail in the Faraday depth spectrum, they also state that to resolve this structure requires that the ratio between the highest and lowest frequencies must be at least 1.5 (their Eq. (19)), which is not satisfied for our observations. As a result we cannot identify this characteristic feature of Faraday caustics, but future observations with larger bandwidth may be able to do so.

Finally, these features could represent sharp “edges” in the intrinsic Faraday spectra, caused by sharp transitions in the Faraday spectrum amplitude at the edges of Faraday thick features. Such transitions would have similar causes as the Faraday-thin features described above: a sharp transition in the free electron density, the strength of the parallel component of the magnetic field, or the synchrotron emissivity. However, such sharp edges tend to be strongly depolarized (Van Eck et al. 2017), which in turn requires a very high intrinsic polarization to produce a detectable feature in the Faraday depth cubes.

We are not able to determine which of these possible causes are the source of the Faraday-thin features we observe. The potential signatures in the other ISM tracers of these conditions are difficult to measure due to limited sensitivity of the tracer observations, or to the fact that most of the tracer data are line-of-sight integrated rather than tomographic measurements.

## 5.2. Gradients in Faraday rotation

The Faraday depth structure seen in the diffuse emission shows a wealth of complex morphologies, such as the linear gradients and local minima/maxima previously described. A detailed analysis of all these individual features is beyond the scope of this work, but we do consider some possible causes of the Faraday depth gradients.

A gradient in Faraday depth can be caused by a gradient in the three factors that determine the Faraday depth: the free electron density, the parallel magnetic field, and the (physical) depth of the Faraday-rotating volume. A gradient in the depth of the Faraday-rotating region would naturally result if the emitting region was a sheet that was not perpendicular to the line of sight and the rotation occurred directly in front of the emitting region; if such an emitting sheet was reasonably flat the distance to the region, and the corresponding Faraday depth, could increase linearly with position. Gradients in the electron density might be expected on larger (kpc) scales, due to effects like the scale height of free electrons, and also on the scales of individual clouds (e.g., a transition layer between a neutral cloud and an ionized skin or exterior).

The linear gradient in the northwest feature (Sect. 3.2.1) contains both negative and positive Faraday depths. The simplest explanation for this is a gradient in the parallel component of the magnetic field that causes the sign to change from negative to positive. Alternative explanations would require two regions at different distances with oppositely directed parallel magnetic fields and an electron density or distance gradient in one of the regions. However, this would produce a Faraday caustic associated with the reversal in parallel magnetic field between these two regions and a corresponding Faraday-thin emission feature (which we do not observe). While there are many possible magnetic field configurations that would produce a gradient in the parallel component of the magnetic field, one possible configuration that would be consistent both with the observations and theoretical models (as discussed in Beck et al. 1996; Shukurov 2004) is a magnetic filament bent into a loop such as shown in Fig. 16.

<figure><img src="image_17.png" /><figcaption>**Fig. 16.** A schematic diagram of a possible source for Faraday depth gradients, where a magnetic loop (black lines) would possess a gradient in the parallel component of the magnetic field and in turn produce a gradient in Faraday depth. An (unrelated) Faraday-thin emission feature in the background would provide the observed source of polarized emission (gray lines).</figcaption></figure>

## 5.3. Southeast Faraday depth gradient-HI correlation

The transition from a smooth linear gradient in Faraday depth outside of the H I filament, shown in Fig. 11, to a more complicated morphology overlapping the filament is very intriguing. Since we are able to observe smooth transitions in the Faraday depth and polarized intensity in the filament region, we conclude that the H I filament is not associated with the source of the polarized emission, and must lie in the foreground with the rest of the Faraday-rotating structure. Furthermore, since the orientation of the gradient so closely aligns with that of the filament boundary, we conclude that the gradient is most likely related to, and co-distant with, the H I filament.

One possible interpretation of the source of the gradient is in terms of an envelope of ionized material surrounding the (predominantly neutral) H I filament. If this is the case, then we expect to see a similar gradient on the other side of the filament; unfortunately that corresponds to the region where instrumental leakage from 3C295 overwhelms the real signal. Improvements in the leakage calibration for LOFAR data should allow for future investigation of this possibility.

# 6. Depolarization by Faraday depth gradients

In addition to constraining possible magnetic field configurations, the Faraday depth gradients seen in this field are noteworthy in that they can be difficult to observe, requiring specific observational parameters. At higher frequencies, it is often the case that the Faraday depth resolution is so low that gradients in Faraday depth can only be measured in cases with very high signal-to-noise ratios. However, at lower frequencies the synthesized beam is typically larger than at higher frequencies, leading to stronger beam depolarization masking the polarized signal and making the gradients undetectable. Appendix A describes this beam depolarization, and gives an approximate threshold for the conditions under which the depolarization becomes significant: polarization angle gradients greater than 180° over the FWHM of the beam will be strongly depolarized, and gradients shallower than this may also be significantly depolarized depending on the exact beam parameters.
We can use this threshold to determine the depolarization caused by the Faraday depth gradients in our observations. For our observations, with a FWHM of 4.3 and wavelengths between 1.9 and 2.5 m, a polarization angle gradient of 180° $FWHM^{-1}$ would be created by a Faraday depth gradient of 7.5–13 rad $m^{-2} deg^{-1}$. In the HETDEX field, we observed two very linear gradients in Faraday depth, with values of 2.6 rad $m^{-2} deg^{-1}$ (Sect. 3.2.1), which is well inside the predicted range of minimal depolarization, and 13 rad $m^{-2} deg^{-1}$ (Sect. 3.2.3), at the upper limit of this range. Beam depolarization by the gradient may explain the lower polarized flux density of the steeper gradient compared to the other polarized features. For these gradients, the Gaussian model predicts that the shallower gradient would be detected at 87% of its intrinsic intensity and the steeper gradient at 3%. Our synthesized beam likely falls somewhere between the two models, implying that the measured brightness of the shallower gradient is probably accurate within approximately 10% while the steeper gradient is probably somewhere between 1–30 times brighter than our measurement.

We can also consider the detectability of Faraday depth gradients with other low-frequency radio observations. Lenc et al. (2016) observed faint diffuse emission with the MWA, at 154 MHz ($\lambda$=1.95 m) with a resolution of 54' and similar sensitivity as our observations. With the simple beam model of Appendix A, the depolarization threshold corresponds to a Faraday depth gradient of 0.9 rad $m^{-2} deg^{-1}$. However, they also applied a Gaussian ($u, v$) taper to their data, which probably has the effect of making their beam more Gaussian-like and increasing the depolarization present. These observations are probably strongly affected by beam depolarization due to Faraday depth gradients; the majority of the structure seen in the HETDEX field could not be detected with these observational parameters. Observing such features may be possible with the MWA by including the longer baselines and minimizing the effects of the Gaussian taper.

There are plans to use the data from the Canadian H I Mapping Experiment (CHIME) to measure polarized foregrounds. CHIME will operate at 400–800 MHz, with a resolution of approximately 30'. The Faraday depth gradient threshold for these parameters varies from 11–45 rad $m^{-2} deg^{-1}$ from the bottom of the frequency band to the top, indicating that these observations should be less affected by gradient beam depolarization than our LOTSS observations.

The sensitivity of observations to steep Faraday depth gradients may be a consideration in planning future observations, as it may have strong consequences on the interpretation of such observations. If polarized features with steep gradients are present in an observed region of sky, they may strongly depolarized and thus not detected. This would lead to an underestimate in the amount of polarized flux present, and could introduce a bias into any inferred properties of the diffuse polarization and ISM magnetic field. Also, polarized features that include regions of both shallow and steep gradients would have missing regions, which could cause large connected Faraday depth structures to be interpreted as separate features.

We have not used the full resolution available with the LOTSS data, so a reprocessing of these data at higher resolution could improve the sensitivity to steeper Faraday depth gradients. At the same resolution as the initial LOTSS Stokes *I* data products, 25", the threshold for depolarization increases to 72–125 rad $m^{-2} deg^{-1}$ at the bottom and top ends of the frequency band. However, the higher resolution results in a correspondingly lower sensitivity to diffuse emission, making it difficult to select a compromise between maximizing sensitivity

to diffuse emission while preventing beam depolarization due to unresolved polarization angle gradients. Higher resolution also significantly increases the computational and data storage requirements of such reprocessing, so future surveys will need to balance computational expense against sensitivity to steeper Faraday depth gradients.

## 7. Conclusions

We have used 60 8-h observations from the LOFAR Two-meter Sky Survey to perform Faraday tomography covering the HETDEX Spring field (right ascension from $10^h30^m$ to $15^h30^m$ and declination from 45° to 57°). We have produced a mosaic Faraday depth cube of this region at 4.3 resolution, which shows polarized emission as a function of Faraday depth. The low-frequency nature of our data gives us a Faraday depth resolution of 1 rad $m^{-2}$, allowing us to probe very small variations in Faraday depth. We achieve a typical sensitivity of 50–100 µJy $PSF^{-1}$ RMSF$^{-1}$.

In our Faraday depth cube we see diffuse polarized emission across most of the region, at Faraday depths between -7 and +25 rad $m^{-2}$. This diffuse emission mostly takes the form of “sheets”, where the emission appears to be filamentary at any single Faraday depth but can be seen as a continuous feature distributed smoothly over Faraday depth. We are able to map out several of these sheets in different positions in the region. A few of these show very linear features, where emission over several square degrees has a smooth linear gradient in Faraday depth.

We compared our Faraday depth cube with several other tracers of different ISM components, and found two interesting relationships. One of the large linear Faraday depth gradients follows the edge of an H I filament, with the gradient occurring just outside the observed edge of the filament with the Faraday depth decreasing with distance from the filament. We also found that the locations of maxima and minima in the Faraday depth of the features we observe tend to overlap with regions of high 21-cm polarized intensity.

Motivated by the Faraday depth gradients seen in our observations, we considered beam depolarization caused by the polarization angle gradients produced by such Faraday depth gradients. Many previous authors have considered the depolarization effects of a Gaussian PSF, which produces very strong depolarization. We have shown that an idealized PSF for an interferometer with uniform ($u, v$) coverage causes much weaker depolarization for some cases. We expect that most realistic observations will fall between these two cases, and suggest that Faraday depth gradients may be more likely to be detected than previously expected.

We have shown that the LOTSS data are well suited for Faraday tomography. However, the quality of the Faraday depth cubes could be improved significantly by performing additional processing to remove the instrumental polarization leakage, for example through the use of complex CLEAN (Pratley & Johnston-Hollitt 2016) to remove image-plane sidelobes. It may also be worthwhile to explore higher (angular) resolution, as this could increase sensitivity to emission with stronger Faraday depth gradients. At the time of this work, LOTSS data with direction-dependent calibration was not available for most observations. Early investigations into the effects of direction-dependent polarization on the polarization properties suggest that there is a modest improvement due to fewer image-plane sidelobes of the instrumental leakage around bright sources (S. O'Sullivan, priv. comm.), which is encouraging for future
polarization work with LOTSS data. LOTSS will observe the
entire sky north of declination zero, so a Faraday tomography
survey, with excellent sensitivity and resolution (both image-
plane and in Faraday depth), will soon be possible.

Acknowledgements. C.V.E. would like to thank Justin Bray and Jeroen Stil for their helpful comments on the derivation in the appendix. This work is part of the research program 639.042.915, which is (partly) financed by the Netherlands Organisation for Scientific Research (NWO). LOFAR, the Low Frequency Array designed and constructed by ASTRON, has facilities in several countries, that are owned by various parties (each with their own funding sources), and that are collectively operated by the International LOFAR Telescope (ILT) foundation under a joint scientific policy. This research used ionospheric TEC maps produced by the Centre for Orbit Determination in Europe (CODE, http://aiuws.unibe.ch/ionosphere/). Some of the ISM tracer maps were downloaded from the Legacy Archive for Microwave Background Data Analysis (LAMBDA), part of the High Energy Astrophysics Science Archive Center (HEASARC). HEASARC/LAMBDA is a service of the Astrophysics Science Division at the NASA Goddard Space Flight Center<fnref n="5" />. This research made extensive use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013); SciPy (Jones et al. 2001); NumPy (van der Walt et al. 2011); IPython (Pérez & Granger 2007); matplotlib (Hunter 2007); the Common Astronomy Software Applications (CASA, McMullin et al. 2007); and the Karma visualization tools (Gooch 1996).

References

Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33

Beck, R., Brandenburg, A., Moss, D., Shukurov, A., & Sokoloff, D. 1996, ARA&A, 34, 155

Bell, M. R., Junklewitz, H., & Enßlin, T. A. 2011, A&A, 535, A85

Beresnyak, A., & Lazarian, A. 2015, Astrophys. Space Sci. Lib., 407, 163

Bernardi, G., Greenhill, L. J., Mitchell, D. A., et al. 2013, ApJ, 771, 105

Boulares, A., & Cox, D. P. 1990, ApJ, 365, 544

Bovy, J., Rix, H.-W., Green, G. M., Schlafly, E. F., & Finkbeiner, D. P. 2016, ApJ, 818, 130

Brentjens, M. A., & de Bruyn A. G. 2005, A&A, 441, 1217

Finkbeiner, D. P. 2003, ApJS, 146, 407

Gooch, R. 1996, in Astronomical Data Analysis Software and Systems V, eds. G. H. Jacoby & J. Barnes, ASP Conf. Ser., 101, 80

Green, G. M., Schlafly, E. F., Finkbeiner, D. P., et al. 2015, ApJ, 810, 25

Hales, C. A., Gaensler, B. M., Norris, R. P., & Middelberg, E. 2012, MNRAS, 424, 2160

Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 1

Heald, G., Braun, R., & Edmonds, R. 2009, A&A, 503, 409

Hill, G. J., Gebhardt, K., Komatsu, E., et al. 2008, in Panoramic Views of Galaxy Formation and Evolution, eds. T. Kodama, T. Yamada, & K. Aoki, ASP Conf. Ser., 399,115

Hill, A. S., Landecker, T. L., Carretti, E., et al. 2017, MNRAS, 467, 4631

Hunter, J. D. 2007, Comput. Sci. & Eng., 9, 90

Iacobelli, M., Haverkorn, M., & Katgert, P. 2013, A&A, 549, A56

Jelić, V., de Bruyn, A. G., Mevius, M., et al. 2014, A&A, 568, A101

Jelić, V., de Bruyn, A. G., Pandey, V. N., et al. 2015, A&A, 583, A137

Jones, E., Oliphant, T., Peterson, P., et al. 2001, SciPy: Open source scientific tools for Python, www.scipy.org

Lallement, R., Vergely, J.-L., Valette, B., et al. 2014, A&A, 561, A91

Landecker, T. L., & Wielebinski, R. 1970, Aust. J. Phys. Astrophys. Suppl., 16, 1

Lenc, E., Gaensler, B. M., Sun, X. H., et al. 2016, ApJ, 830, 38

Macquart, J.-P., Ekers, R. D., Feain, I., & Johnston-Hollitt, M. 2012, ApJ, 750, 139

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astronomical Data Analysis Software and Systems XVI, eds. R. A. Shaw, F. Hill, & D. J. Bell, ASP Conf. Ser., 376, 127

Pandey, V. N., van Zwieten, J. E., de Bruyn, A. G., & Nijboer, R. 2009, in The Low-Frequency Radio Universe, eds. D. J. Saikia, D. A. Green, Y. Gupta, & T. Venturi, ASP Conf. Ser., 407, 384

Pérez, F., & Granger, B. E. 2007, Comput. Sci. Eng., 9, 21

Planck Collaboration X. 2016, A&A, 594, A10

Patley, L., & Johnston-Hollitt, M. 2016, MNRAS, 462, 3483

Reich, W. 1982, A&AS, 48, 219

Remazeilles, M., Dickinson, C., Banday, A. J., Bigot-Sazy, M.-A., & Ghosh, T. 2015, MNRAS, 451, 4311

Schnitzeler, D. H. F. M., & Lee, K. J. 2015, MNRAS, 447, L26

Schnitzeler, D. H. F. M., Katgert, P., & de Bruyn A. G. 2007, A&A, 471, L21

Schnitzeler, D. H. F. M., Katgert, P., & de Bruyn A. G. 2009, A&A, 494, 611

Shimwell, T. W., Röttgering, H. J. A., Best, P. N., et al. 2017, A&A, 598, A104

Shukurov, A. 2004, ArXiv Astrophysics e-prints [arXiv:astro-ph/0411739]

Sokoloff, D. D., Bykov, A. A., Shukurov, A., et al. 1998, MNRAS, 299, 189

Spinrad, H., Marr, J., Aguilar, L., & Djorgovski, S. 1985, PASP, 97, 932

Tasse, C., van der Tol, S., van Zwieten, J., van Diepen, G., & Bhatnagar, S. 2013, A&A, 553, A105

Tingay, S. J., Goeke, R., Bowman, J. D., et al. 2013, PASA, 30, e007

Tribble, P. C. 1991, MNRAS, 250, 726

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Comput. Sci. Eng., 13, 22

Van Eck, C. L., Haverkorn, M., Alves, M. I. R., et al. 2017, A&A, 597, A98

Van Eck, C. L., Haverkorn, M., Alves, M. I. R., et al. 2018, A&A, 613, A58

van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2

van Loo, S., Hartquist, T. W., & Falle, S. A. E. G. 2012, Astron. Geophys., 53, 5,31

Winkel, B., Kerp, J., Flöer, L., et al. 2016, A&A, 585, A41

Wolleben, M., Landecker, T. L., Reich, W., & Wielebinski, R. 2006, A&A, 448, 411

Wolleben, M., Landecker, T. L., Carretti, E., et al. 2008, Proc. Int. Astron. Union, 4, 89

Zaroubi, S., Jelić, V., de Bruyn, A. G., et al. 2015, MNRAS, 454, L46

<footnote n="5">⁵ https://lambda.gsfc.nasa.gov/</footnote>
## Appendix A: Beam depolarization in a linear Faraday depth gradient

Gradients in Faraday depth with respect to position on the sky result in corresponding gradients in the polarization angle (if the intrinsic emitted angles are uniform across the source), which in turn produce beam depolarization when observed with a finite resolution. Here we explore the depolarization caused by Faraday depth gradients by considering the case of an infinite, uniform background emission source, with a (foreground) gradient in Faraday depth, corresponding to the sheets of emission and Faraday depth gradients we observe in our data. When a Gaussian beam is assumed, an analytical solution can be found, as described in Sokoloff et al. (1998). However, more realistic synthesized beams have not been considered in the literature, and below we show that this can make a significant difference in the expected depolarization.

Consider a gradient in Faraday depth with respect to an angular position variable $x$, $\frac{d\phi}{dx}$. The resulting intrinsic (sky) polarization, in the complex-polarization notation<fnref n="6" />, is

$$ \tilde{P}_{\text{sky}}(x, y, \lambda) = \tilde{P}_0 \exp \left( 2i\lambda^2 \frac{d\phi}{dx} x \right) \quad (A.1) $$

where $\lambda$ is the observing wavelength and $\tilde{P}_0$ is the intrinsic emitted polarization prior to Faraday rotation<fnref n="7" />. The observed polarization is the convolution of the intrinsic polarization with the telescope beam. In the case of a circular Gaussian beam ($\frac{1}{2\pi\sigma^2} \exp(-\frac{x^2+y^2}{2\sigma^2})$), the convolution integral can be solved analytically, resulting in the previously known result

$$ \tilde{P}_{\text{obs}}(x, y, \lambda) = \tilde{P}_0 \exp \left( 2i\lambda^2 \frac{d\phi}{dx} x \right) \exp \left( -2 \left( \lambda^2 \frac{d\phi}{dx} \sigma \right)^2 \right). \quad (A.2) $$

The first part of this equation is identical to the intrinsic polarization and contains all the position dependence and the only complex terms, which means that the correct polarization angle, Faraday depth, and Faraday depth gradient are recovered, but the polarized intensity is modified by the later parts of the equation. The second part of the equation gives the depolarization, and shows that the depolarization is a very strong function of the wavelength, gradient, and beam size.

However, a Gaussian is often not an accurate representation of the synthesized beam, depending on the $(u, v)$ coverage of the instrument used. Previous calculations such as Tribble (1991), while not for the exact same model, have shown that Gaussian functions tend to cause much stronger depolarization than other functions. With this motivation, we have developed a more general formulation for beam depolarization of gradients.

The general case, and the more natural formulation for interferometric observations, can be found simply by exploiting properties of Fourier transforms. Since the observed polarization is the convolution of the sky polarization with the telescope beam, we can move to the Fourier domain where the convolution becomes a multiplication and the beam becomes the $(u, v)$ coverage of the observation:

$$ \mathcal{F}\{\tilde{P}_{\text{obs}}(x, y)\} = \mathcal{F}\{\tilde{P}_{\text{sky}}(x, y) * \text{Beam}(x, y)\}, \quad (A.3) $$

$$ = \mathcal{F}\{\tilde{P}_{\text{sky}}(x, y)\} \mathcal{F}\{\text{Beam}(x, y)\}, \quad (A.4) $$

$$ = \tilde{P}_0 \delta \left( u - \frac{2\lambda^2}{2\pi} \frac{d\phi}{dx} \right) \delta(v) W(u, v), \quad (A.5) $$

where $W(u, v)$ is the sampling weight function in the $(u, v)$ plane, and in the last step we have used the definition of $\tilde{P}_{\text{sky}}(x)$ from Eq. (A.1) and recalled that the Fourier transform of a complex exponential is a suitably shifted Dirac delta function. Inverting the Fourier transform to recover the observed sky distribution gives:

$$ \tilde{P}_{\text{obs}}(x, y) = \mathcal{F}^{-1} \left\{ \tilde{P}_0 \delta \left( u - \frac{2\lambda^2}{2\pi} \frac{d\phi}{dx} \right) \delta(v) W(u, v) \right\}, \quad (A.6) $$

$$ = \tilde{P}_0 \exp \left( 2i\lambda^2 \frac{d\phi}{dx} x \right) W \left( \frac{\lambda^2}{\pi} \frac{d\phi}{dx}, 0 \right). \quad (A.7) $$

To this point, the result is general and works for any sampling function; using Gaussian weights gives the result from Eq. (A.2). The correct Faraday depth, and by extension the gradient in Faraday depth, is recovered at all positions, but the polarized intensity is modified by the weighting function at the location of the gradient in the $(u, v)$ plane. As a result, gradients at positions in the $(u, v)$ plane that are not sampled by the instrument or that have been down-weighted significantly will have zero or strongly reduced polarized intensity and may not be detected; for such an idealized gradient the beam depolarization depends only on the weight function.

The steepest gradient that can be measured for a given observation can be determined from three parameters: the observing wavelength, $\lambda$, the largest sampled $(u, v)$ coordinate, $UV_{\text{max}}$ (normalized by wavelength), and the Faraday depth gradient, $\frac{d\phi}{dx}$. We found that we could reduce this to one parameter by converting our model into scale-free units: we replaced the $UV_{\text{max}}$ parameter with the synthesized beam full-width half max (FWHM), which is defined as $FWHM = UV_{\text{max}}^{-1}$; replaced the $\lambda^2 \frac{d\phi}{dx}$ term in the intrinsic polarization with a single variable, the polarization angle gradient (with units of radians of rotation per beam FWHM); and expressed the position variable $x$ in terms of the FWHM. The resulting condition for being observable is $\frac{d(\phi, \lambda^2)}{dx} FWHM < \pi$, or equivalently less than a 180° change in polarization angle across the FWHM of the synthesized beam.

To derive Eq. (A.7), we assumed that the source of the polarized emission has an infinite extent, so that it shows up in the $(u, v)$ plane at a single point. Real sources, which have a finite extent, cover an area in the $uv$-plane, not a single point. A more general argument can be constructed from Fourier properties, using the ideal case as a guide. Consider a polarized intensity feature of arbitrary morphology and angular size (and its corresponding Fourier transform to the $(u, v)$ plane), and apply a linear Faraday depth gradient of the form in Eq. (A.1) which has the effect of adding a position dependent phase modulation. From the Fourier shift-modulation theorem, the net effect is to shift the Fourier transform of the polarized feature in the $(u, v)$ plane, with the length and orientation of the shift determined by the steepness and orientation of the Faraday depth gradient. This shift can cause some or all of the polarized feature to fall outside of the $(u, v)$ coverage of the observation, causing depolarization.

This result mirrors that of Schnitzeler et al. (2009), who used the same gradient model and Fourier properties to derive the corresponding result for the shallowest gradient that can be detected, given the shortest baseline present in an observation. These can be interpreted together: the longest baseline sets the steepest Faraday depth gradient that can be observed, while

<footnote n="6">⁶ We use a tilde to denote complex quantities, specifically phasors describing linear polarization.</footnote>

<footnote n="7">⁷ For convenience we have aligned the gradient with the x-axis and set $\phi=0$ at $x=0$, but these are arbitrary and do not affect the final result.</footnote>
the shortest baseline sets the shallowest gradient that can be
observed. Single-dish observations, which measure around the
center of the (u, v) plane, are sensitive to large polarized features
with small Faraday depth gradients.

It should also be noted that this derivation describes beam
depolarization in a single frequency channel. This depolariza-
tion may have a very strong frequency dependence; as frequency
decreases the location of the gradient in the (u, v) plane will
move outwards to larger (u, v) distances, scaling as the square
of the wavelength, while a baseline of a fixed length moves
inwards proportionally to the wavelength. For a given obser-
vation there will be a range of gradients that can be detected
in a fraction of the observed bandwidth, producing a complex
frequency-dependence which may introduce artificial structure
into the corresponding Faraday depth spectrum. Determining
the specific effects of a “partial detection” of a gradient is beyond
the scope of this derivation.

Overall, we find three important results from this derivation.
First, we have identified a key figure of merit for characterizing
Faraday depth gradients for the purpose of evaluating their depo-
larization within a given observation: the change in polarization

angle across the FWHM of the synthesized beam, with gradi-
ents steeper than 180° per beam FWHM being beyond the (u, v)
coverage of the observation

Second, beam depolarization may often not be as strong
as may have been previously predicted using a Gaussian beam
model. Just below the threshold of 180° per FWHM the depolar-
ization of an ideal gradient by the Gaussian beam would leave
only 2.9% of the original polarized flux, whereas a uniform-
weighted observation with the same beam FWHM would not be
depolarized. The effects will be less extreme for a more realistic
gradient, depending on the weighting of the (u, v) coverage of the
observation and the angular extent of the gradient.

Third, and perhaps most importantly, this suggests that the
choice of (u, v) weighting in the imaging process can have very
strong effects on the measured polarization. A full study of this
effect is beyond the scope of this work, but these simple mod-
els imply that the use of Gaussian tapers and similar weighting
schemes that highly down-weight longer baselines may cause
significant depolarization; the use of such tapers should be con-
sidered cautiously to avoid removing possible polarized features
of interest.