We investigate the statistical properties of the polarized emission of extragalactic radio sources and estimate their contribution to the power spectrum of polarization fluctuations in the microwave region. The basic ingredients of our analysis are the NVSS polarization data, the multifrequency study of polarization properties of the B3-VLA sample (Mack et al. 2002) which has allowed us to quantify Faraday depolarization effects, and the 15 GHz survey by Taylor et al. (2001), which has provided strong constraints on the high-frequency spectral indices of sources. The polarization degree of both steep- and flat-spectrum at 1.4 GHz is found to be anti-correlated with the flux density. The median polarization degree at 1.4 GHz of both steep- and flat-spectrum sources brighter than S(1.4GHz)=80mJy is ≃2.2. The data by Mack et al. (2002) indicate a substantial mean Faraday depolarization at 1.4 GHz for steep spectrum sources, while the depolarization is undetermined for most flat/inverted-spectrum sources. Exploiting this complex of information we have estimated the power spectrum of polarization fluctuations due to extragalactic radio sources at microwave frequencies. We confirm that extragalactic sources are expected to be the main contaminant of Cosmic Microwave Background (CMB) polarization maps on small angular scales. At frequencies <30GHz the amplitude of their power spectrum is expected to be comparable to that of the E-mode of the CMB. At higher frequencies, however, the CMB dominates.
