We present a catalog of Faraday rotation measures (RMs) and redshifts for 4003 extragalactic radio sources detected at 1.4 GHz, derived by identifying optical counterparts and spectroscopic redshifts for linearly polarized radio sources from the NRAO VLA Sky Survey. This catalog is more than an order of magnitude larger than any previous sample of RM vs. redshift, and covers the redshift range 0 < z < 5.3 ; the median redshift of the catalog is z = 0.70, and there are more than 1500 sources at redshifts z > 1. For 3650 of these sources at Galactic latitudes |b| >= 20 degrees, we present a second catalog in which we have corrected for the foreground Faraday rotation of the Milky Way, resulting in an estimate of the residual rotation measure (RRM) that aims to isolate the contribution from extragalactic magnetic fields. We find no significant evolution of RRM with redshift, but observe a strong anti-correlation between RRM and fractional polarization, p, that we argue is the result of beam depolarization from small-scale fluctuations in the foreground magnetic field or electron density. We suggest that the observed variance in RRM and the anti-correlation of RRM with p both require a population of magnetized intervening objects that lie outside the Milky Way but in the foreground to the emitting sources.

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.