The extraction of cosmological parameters from microwave background observations relies on specific assumptions about the statistical properties of the data, in particular that the p-point distributions of temperature fluctuations are jointly-normal. Using a battery of statistical tests, we assess the multivariate Gaussian nature of the Wilkinson Microwave Anisotropy Probe (WMAP) 1st year data. The statistics we use fall into three classes which test different aspects of joint-normality: the first set assess the normality of marginal (one-point) distributions using familiar univariate methods; the second involves statistics that directly assess joint-normality; and the third explores the evidence of non-linearity in the relationship between variates. We applied these tests to frequency maps, `foreground-cleaned' assembly maps and all-sky CMB-only maps. The assembly maps are of particular interest as when combined with the kp2 mask, we recreate the region used in the computation of the angular power spectrum. Significant departures from normality were found in all the maps. In particular, the kurtosis coefficient, D'Agostino's statistic and bivariate kurtosis calculated from temperature pairs extracted from all the assembly maps were found to be non-normal at 99% confidence level. We found that the results were unaffected by the size of the Galactic cut and were evident on either hemisphere of the CMB sky. The latter suggests that the non-Gaussianity is not simply related to previous claims of north-south asymmetry or localized abnormalities detected through wavelet techniques.

Using a set of compilations of measurements for extragalactic radio sources we construct all-sky maps of the Faraday Rotation produced by the Galactic magnetic field. In order to generate the maps we treat the radio source positions as a kind of "mask" and construct combinations of spherical harmonic modes that are orthogonal on the masked sky. As long as relatively small multipoles are used the resulting maps are quite stable to changes in selection criteria for the sources, and show clearly the structure of the local Galactic magnetic field. We also suggest the use of these maps as templates for CMB foreground analysis, illustrating the idea with a cross-correlation analysis between the Wilkinson Microwave Anisotropy Probe (WMAP) data and our maps. We find a significant cross-correlation, indicating the presence of significant residual contamination.

We present a new technique for constraining the topology of the universe. The method exploits the existence of correlations in the phases of the spherical harmonic coefficients of the CMB temperature pattern associated with matched pairs of circles seen in the sky in universes with non-trivial topology. The method is computationally faster than all other statistics developed to hunt for these matched circles. We applied the method to a range of simulations with topologies of various forms and on different scales. A characteristic form of phase correlation is found in the simulations. We also applied the method to preliminary CMB maps derived from WMAP, but the separation of topological effects from, e.g., foregrounds is not straightforward.

We study the statistical properties of spherical harmonic modes of temperature maps of the cosmic microwave background. Unlike other studies, which focus mainly on properties of the amplitudes of these modes, we look instead at their phases. In particular, we present a simple measure of phase correlation that can be diagnostic of departures from the standard assumption that primordial density fluctuations constitute a statistically homogeneous and isotropic Gaussian random field, which should possess phases that are uniformly random on the unit circle. The method we discuss checks for the uniformity of the distribution of phase angles using a non-parametric descriptor based on the use order statistics, which is known as Kuiper's statistic. The particular advantage of the method we present is that, when coupled to the judicious use of Monte Carlo simulations, it can deliver very interesting results from small data samples. In particular, it is useful for studying the properties of spherical harmonics at low l for which there are only small number of independent values of m and which therefore furnish only a small number of phases for analysis. We apply the method to the COBE-DMR and WMAP sky maps, and find departures from uniformity in both. In the case of WMAP, our results probably reflect Galactic contamination or the known variation of signal-to-noise across the sky rather than primordial non-Gaussianity.

We present a diagnostic test of possible Galactic contamination of cosmic microwave background sky maps designed to provide an independent check on the methods used to compile these maps. The method involves a non-parametric measurement of cross-correlation between the Faraday rotation measure (RM) of extragalactic sources and the measured microwave signal at the same angular position. We argue that statistical properties of the observed distribution of rotation measures are consistent with a Galactic origin, an argument reinforced by a direct measurement of cross-correlation between dust, free-free and synchrotron foreground maps and RM values with the strongest correlation being for dust and free-free. We do not find any statistically compelling evidence for correlations between the RM values and the COBE DMR maps at any frequency, so there is no evidence of residual contamination in these CMB maps. On the other hand, there is a statistically significant correlation of RM with the preliminary WMAP individual frequency maps which remains significant in the Tegmark et al. Wiener-filtered map but not in the Internal Linear Combination map produced by the WMAP team.