The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.

The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of $\Delta^2 \leq 2.7 \times 10^4$ mK$^2$ at $k=0.27$ h~Mpc$^{-1}$ and $z=7.1$, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications.

We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. The sensitivity available in our deep observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ~54' resolution with a sensitivity of 5.9 mJy beam^-1 and compact polarized sources at ~2.4' resolution with a sensitivity of 2.3 mJy beam^-1 for a subset (400 deg^2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (~1 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from -2 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+1 rad m^-2. We estimate a distance of 51+/-20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.

Detection of the Epoch of Reionization HI signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a catalogue of 7394 extragalactic sources at 182 MHz detected in the RA=0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of -0.85. Spectral flattening is seen toward lower frequencies with a median of -0.71 predicted at 182 MHz. The astrometric error is 7 arcsec. compared to a 2.3 arcmin. beam FWHM. The resulting catalogue covers approximately 1400 sq. deg. and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.

We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple, independent, data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds.

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). Three hours of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1-3 hours of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 18$ but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window.

In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum.

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

A new generation of low frequency radio telescopes is seeking to observe the redshifted 21 cm signal from the Epoch of Reionization (EoR), requiring innovative methods of calibration and imaging to overcome the difficulties of widefield low frequency radio interferometry. Precise calibration will be required to separate the small expected EoR signal from the strong foreground emission at the frequencies of interest between 80 and 300 MHz. The Moon may be useful as a calibration source for detection of the EoR signature, as it should have a smooth and predictable thermal spectrum across the frequency band of interest. Initial observations of the Moon with the Murchison Widefield Array 32 tile prototype show that the Moon does exhibit a similar trend to that expected for a cool thermally emitting body in the observed frequency range, but that the spectrum is corrupted by reflected radio emission from Earth. In particular, there is an abrupt increase in the observed flux density of the Moon within the internationally recognised Frequency Modulated (FM) radio band. The observations have implications for future low frequency surveys and EoR detection experiments that will need to take this reflected emission from the Moon into account. The results also allow us to estimate the equivalent isotropic power emitted by the Earth in the FM band and to determine how bright the Earth might appear at metre wavelengths to an observer beyond our own solar system.

The Murchison Widefield Array (MWA) is a new low-frequency, wide field-of-view radio interferometer under development at the Murchison Radio-astronomy Observatory (MRO) in Western Australia. We have used a 32-element MWA prototype interferometer (MWA-32T) to observe two 50-degree diameter fields in the southern sky in the 110 MHz to 200 MHz band in order to evaluate the performance of the MWA-32T, to develop techniques for epoch of reionization experiments, and to make measurements of astronomical foregrounds. We developed a calibration and imaging pipeline for the MWA-32T, and used it to produce ~15' angular resolution maps of the two fields. We perform a blind source extraction using these confusion-limited images, and detect 655 sources at high significance with an additional 871 lower significance source candidates. We compare these sources with existing low-frequency radio surveys in order to assess the MWA-32T system performance, wide field analysis algorithms, and catalog quality. Our source catalog is found to agree well with existing low-frequency surveys in these regions of the sky and with statistical distributions of point sources derived from Northern Hemisphere surveys; it represents one of the deepest surveys to date of this sky field in the 110 MHz to 200 MHz band.

It is shown that the excellent Murchison Radio-astronomy Observatory site allows the Murchison Widefield Array to employ a simple RFI blanking scheme and still calibrate visibilities and form images in the FM radio band. The techniques described are running autonomously in our calibration and imaging software, which is currently being used to process an FM-band survey of the entire southern sky.

The physical properties of galactic cirrus emission are not well characterized. BOOMERanG is a balloon-borne experiment designed to study the cosmic microwave background at high angular resolution in the millimeter range. The BOOMERanG 245 and 345GHz channels are sensitive to interstellar signals, in a spectral range intermediate between FIR and microwave frequencies. We look for physical characteristics of cirrus structures in a region at high galactic latitudes (b~-40{\deg}) where BOOMERanG performed its deepest integration, combining the BOOMERanG data with other available datasets at different wavelengths. We have detected eight emission patches in the 345 GHz map, consistent with cirrus dust in the Infrared Astronomical Satellite maps. The analysis technique we have developed allows to identify the location and the shape of cirrus clouds, and to extract the flux from observations with different instruments at different wavelengths and angular resolutions. We study the integrated flux emitted from these cirrus clouds using data from Infrared Astronomical Satellite (IRAS), DIRBE, BOOMERanG and Wilkinson Microwave Anisotropy Probe (WMAP) in the frequency range 23-3000 GHz (13 mm 100 micron wavelength). We fit the measured spectral energy distributions with a combination of a grey body and a power-law spectra considering two models for the thermal emission. The temperature of the thermal dust component varies in the 7 - 20 K range and its emissivity spectral index is in the 1 - 5 range. We identified a physical relation between temperature and spectral index as had been proposed in previous works. This technique can be proficiently used for the forthcoming Planck and Herschel missions data.

In order to construct accurate point sources simulations at the frequencies relevant to 21-cm experiments, the angular correlation of radio sources must be taken into account. This paper presents a measurement of angular two-point correlation function, w(\theta), at 232 MHz from the MIYUN survey - tentative measurements of w(\theta) are also performed at 151 MHz. It is found that double power law with shape w(\theta) = A \theta^{-\gamma} fits the 232 MHz data well. For the angular lenght of 0.2 degrees < \theta < 0.6 degrees, \gamma ~ -1.12, and this value of slope is independent of the flux-density threshold; while for angular lenghts much greater than 0.6 degrees, \gamma has a shallower value of about -0.16. By comparing the results of this paper with previous measurements of w(\theta), it is discussed how w(\theta) changes with the change of frequency and completness limit.

In order to construct accurate point sources simulations at the frequencies relevant to 21 cm experiments, the angular correlation of radio sources must be taken into account. Using the 74 MHz VLSS survey, we measured the angular 2-point correlation function, w(\theta). We obtain the first measurement of clustering at the low frequencies relevant to 21 cm tomography. We find that a single power law with shape w(\theta) = A \theta^{-\gamma} fits well the data. For a galactic cut of 10 degrees, with a data cut of \delta less than -10 degrees, and a flux limit of S = 770 mJy, we obtain a slope of \gamma = (-1.2 +/- 0.35). This value of \gamma is consistent with that measured from other radio catalogues at the millimeter wavelengths. The amplitude of clustering has a length of 0.2 degrees - 0.6 degrees, and it is independent of the flux-density threshold.

The Sunyaev-Zel'dovich (SZ) effect is the inverse Compton-scattering of cosmic microwave background (CMB) photons by hot electrons in the intervening gas throughout the universe. The effect has a distinct spectral signature that allows its separation from other signals in multifrequency CMB datasets. Using CMB anisotropies measured at three frequencies by the BOOMERanG 2003 flight we constrain SZ fluctuations in the 10 arcmin to 1 deg angular range. Propagating errors and potential systematic effects through simulations, we obtain an overall upper limit of 15.3 uK (2 sigma) for rms SZ fluctuations in a broad bin between multipoles of of 250 and 1200 at the Rayleigh-Jeans (RJ) end of the spectrum. When combined with other CMB anisotropy and SZ measurements, we find that the local universe normalization of the density perturbations is sigma-8(SZ) < 0.96 at the 95% confidence level, consistent with sigma-8 determined from primordial perturbations.

We use Minkowski Functionals (MF) to constrain a primordial non-Gaussian contribution to the CMB intensity field as observed in the 150 GHz and 145 GHz BOOMERanG maps from the 1998 and 2003 flights, respectively, performing for the first time a joint analysis of the two datasets. A perturbative expansion of the MF formulae in the limit of a weakly non-Gaussian field yields analytical formulae, derived by Hikage et al. (2006), which can be used to constrain the coupling parameter f_NL without the need for non-Gaussian simulations. We find -1020<f_NL<390 at 95% CL, significantly improving the previous constraints by De Troia et al. (2007) on the BOOMERanG 2003 dataset. These are the best f_NL limits to date for suborbital probes.

Understanding diffuse Galactic radio emission is interesting both in its own right and for minimizing foreground contamination of cosmological measurements. Cosmic Microwave Background experiments have focused on frequencies > 10 GHz, whereas 21 cm tomography of the high redshift universe will mainly focus on < 0.2 GHz, for which less is currently known about Galactic emission. Motivated by this, we present a global sky model derived from all publicly available total power large-area radio surveys, digitized with optical character recognition when necessary and compiled into a uniform format, as well as the new Villa Elisa data extending the 1.4 GHz map to the entire sky. We quantify statistical and systematic uncertainties in these surveys by comparing them with various global multi-frequency model fits. We find that a principal component based model with only three components can fit the 11 most accurate data sets (at 10, 22, 45 & 408 MHz and 1.4, 2.3, 23, 33, 41, 61, 94 GHz) to an accuracy around 1%-10% depending on frequency and sky region. Both our data compilation and our software returning a predicted all-sky map at any frequency from 10 MHz to 100 GHz are publicly available at http://space.mit.edu/home/angelica/gsm .

We analyze the BOOMERanG 2003 (B03) 145 GHz temperature map to constrain the amplitude of a non Gaussian, primordial contribution to CMB fluctuations. We perform a pixel space analysis restricted to a portion of the map chosen in view of high sensitivity, very low foreground contamination and tight control of systematic effects. We set up an estimator based on the three Minkowski functionals which relies on high quality simulated data, including non Gaussian CMB maps. We find good agreement with the Gaussian hypothesis and derive the first limits based on BOOMERanG data for the non linear coupling parameter f_NL as -300<f_NL<650 at 68% CL and -800<f_NL<1050 at 95% CL.

Most analysis of Cosmic Microwave Background spherical harmonic coefficients a_lm has focused on estimating the power spectrum C_l=<|a_lm|^2> rather than the coefficients themselves. We present a minimum-variance method for measuring a_lm given anisotropic noise, incomplete sky coverage and foreground contamination, and apply it to the WMAP data. Our method is shown to constitute lossless data compression in the sense that the widely used quadratic estimators of the power spectrum C_l can be computed directly from our a_lm-estimators. As the Galactic cut is increased, the error bars Delta-a_lm on low multipoles go from being dominated by foregrounds to being dominated by sample variance from other multipoles, with the intervening minimum defining the optimal cut. Applying our method to the WMAP quadrupole and octopole, we find that their previously reported "axis of evil" alignment appears to be rather robust to Galactic cut and foreground contamination.

We present the BOOMERanG-03 experiment and maps of the Stokes parameters I, Q, U of the microwave sky obtained during a 14 day balloon flight in 2003. Three regions of the southern sky were surveyed: a deep survey (~ 90 square degrees) and a shallow survey (~ 750 square degrees) at high Galactic latitudes (both centered at RA ~ 5.5 h, dec ~ -45 deg) and a survey of ~ 300 square degrees across the Galactic plane at RA ~ 9.1 h, dec ~ -47 deg. All three surveys were carried out in three wide frequency bands centered at 145, 245 and 345 GHz, with an angular resolution of ~ 10'. The 145 GHz maps of Stokes I are dominated by Cosmic Microwave Background (CMB) temperature anisotropy, which is mapped with high signal to noise ratio. The measured anisotropy pattern is consistent with the pattern measured in the same region by BOOMERanG-98 and by WMAP. The 145 GHz maps of Stokes Q and U provide a robust statistical detection of polarization of the CMB when subjected to a power spectrum analysis. This amplitude of the polarization is consistent with that of the CMB in the $\Lambda$CDM cosmological scenario. At 145 GHz, in the CMB surveys, the intensity and polarization of the astrophysical foregrounds are found to be negligible with respect to the cosmological signal. At 245 and 345 GHz we detect ISD emission correlated to the 3000 GHz IRAS/DIRBE maps, and give upper limits for any other non-CMB component. We also present intensity maps of the surveyed section of the Galactic plane. These are compared to monitors of different interstellar components, showing that a variety of emission mechanisms is present in that region.

We report measurements of the CMB polarization power spectra from the January 2003 Antarctic flight of BOOMERANG. The primary results come from six days of observation of a patch covering 0.22% of the sky centered near R.A. = 82.5 deg., Dec= -45 deg. The observations were made using four pairs of polarization sensitive bolometers operating in bands centered at 145 GHz. Using two independent analysis pipelines, we measure a non-zero <EE> signal in the range 100< l <1000 with a significance 4.8-sigma, a 2-sigma upper limit of 8.6 uK^2 for any <BB> contribution, and a 2-sigma upper limit of 7.0 uK^2 for the <EB> spectrum. Estimates of foreground intensity fluctuations and the non-detection of <BB> and <EB> signals rule out any significant contribution from galactic foregrounds. The results are consistent with a Lambda-CDM cosmology seeded by adiabatic perturbations. We note that this is the first detection of CMB polarization with bolometric detectors.

We present a measurement of the temperature-polarization angular cross power spectrum, <TE>, of the Cosmic Microwave Background. The result is based on $\sim 200$ hours of data from 8 polarization sensitive bolometers operating at 145 GHz during the 2003 flight of BOOMERANG. We detect a significant <TE> correlation in the $\ell$-range between 50 and 950 with a statistical significance > 3.5 sigma. Contamination by polarized foreground emission and systematic effects are negligible in comparison with statistical uncertainty. The spectrum is consistent with previous detections and with the "concordance model" that assumes adiabatic initial conditions. This is the first measurement of <TE> using bolometric detectors.

We present the cosmological parameters from the CMB intensity and polarization power spectra of the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains the parameters of the $\Lambda$CDM model remarkably well and is consistent with constraints from a multi-experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the combined CMB data set and test several extensions to the standard model including: running of the spectral index, curvature, tensor modes, the effect of massive neutrinos, and an effective equation of state for dark energy. We also include an analysis of constraints to a model which allows a CDM isocurvature admixture.

We report on observations of the Cosmic Microwave Background (CMB) obtained during the January 2003 flight of Boomerang . These results are derived from 195 hours of observation with four 145 GHz Polarization Sensitive Bolometer (PSB) pairs, identical in design to the four 143 GHz Planck HFI polarized pixels. The data include 75 hours of observations distributed over 1.84% of the sky with an additional 120 hours concentrated on the central portion of the field, itself representing 0.22% of the full sky. From these data we derive an estimate of the angular power spectrum of temperature fluctuations of the CMB in 24 bands over the multipole range (50 < l < 1500). A series of features, consistent with those expected from acoustic oscillations in the primordial photon-baryon fluid, are clearly evident in the power spectrum, as is the exponential damping of power on scales smaller than the photon mean free path at the epoch of last scattering (l > 900). As a consistency check, the collaboration has performed two fully independent analyses of the time ordered data, which are found to be in excellent agreement.

The DASI discovery of CMB polarization has opened a new chapter in cosmology. Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and remove polarized foregrounds. We use the POLAR experiment, COBE/DMR and radio surveys to provide the strongest limits to date on the TE cross power spectrum of the CMB on large angular scales and to quantify the polarized synchrotron radiation, which is likely to be the most challenging polarized contaminant for the MAP satellite. We find that the synchrotron E- and B-contributions are equal to within 10% from 408-820 MHz with a hint of E-domination at higher frequencies. We quantify Faraday Rotation and Depolarization effects in the two-dimensional (l,nu)-plane and show that they cause the synchrotron polarization percentage to drop both towards lower frequencies and towards lower multipoles.

Great efforts are currently being channeled into ground- and balloon-based CMB experiments, mainly to explore anisotropy on small angular scales and polarization. To optimize instrumental design and assess experimental prospects, it is important to understand in detail the atmosphere-related systematic errors that limit the science achievable with new suborbital instruments. For this purpose, we spatially compare the 648 square degree ground- and balloon-based QMASK map with the atmosphere-free WMAP map, finding beautiful agreement on all angular scales where both are sensitive. This is a reassuring quantitative assessment of the power of the state-of-the-art FFT- and matrix-based mapmaking techniques that have been used for QMASK and virtually all subsequent experiments.

The DASI discovery of CMB polarization, confirmed by WMAP, has opened a new chapter in cosmology. Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst. The goal of the present review is to provide the state-of-the-art of the CMB polarization from a practical point of view, connecting real-world data to physical models. We present the physics of this polarized phenomena and illustrate how it depends of various cosmological parameters for standard adiabatic models. We also present all observational constraints to date and discuss how much we have learned about polarized foregrounds so far from the CMB studies. Finally, we comment on future prospects for the measurement of CMB polarization.

The WMAP team has produced a foreground map that can account for most of the low-frequency Galactic microwave emission in the WMAP maps, tentatively interpreting it as synchrotron emission. Finkbeiner and collaborators have challenged these conclusions, arguing that the WMAP team "synchrotron" template is in fact not dominated by synchrotron radiation, but by some dust-related Galactic emission process, perhaps spinning dust grains, making dramatically different predictions for its behavior at lower frequencies. By cross-correlating this "synchrotron" template with 10 and 15 GHz CMB observations, we find that its spectrum turns over in a manner consistent with spinning dust emission, falling about an order of magnitude below what the synchrotron interpretation would predict.

We investigate anomalies reported in the Cosmic Microwave Background maps from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite on very large angular scales and discuss possible interpretations. Three independent anomalies involve the quadrupole and octopole: 1. The cosmic quadrupole on its own is anomalous at the 1-in-20 level by being low (the cut-sky quadrupole measured by the WMAP team is more strikingly low, apparently due to a coincidence in the orientation of our Galaxy of no cosmological significance); 2. The cosmic octopole on its own is anomalous at the 1-in-20 level by being very planar; 3. The alignment between the quadrupole and octopole is anomalous at the 1-in-60 level. Although the a priori chance of all three occurring is 1 in 24000, the multitude of alternative anomalies one could have looked for dilutes the significance of such a posteriori statistics. The simplest small universe model where the universe has toroidal topology with one small dimension of order half the horizon scale, in the direction towards Virgo, could explain the three items above. However, we rule this model out using two topological tests: the S-statistic and the matched circle test.

BOOMERANG is a balloon-borne telescope designed for long duration (LDB) flights around Antarctica. The second LDB Flight of BOOMERANG took place in January 2003. The primary goal of this flight was to measure the polarization of the CMB. The receiver uses polarization sensitive bolometers at 145 GHz. Polarizing grids provide polarization sensitivity at 245 and 345 GHz. We describe the BOOMERANG telescope noting changes made for 2003 LDB flight, and discuss some of the issues involved in the measurement of polarization with bolometers. Lastly, we report on the 2003 flight and provide an estimate of the expected results.

We perform an independent foreground analysis of the WMAP maps to produce a cleaned CMB map (available online) useful for cross-correlation with, e.g., galaxy and X-ray maps. We use a variant of the Tegmark & Efstathiou (1996) technique that is completely blind, making no assumptions about the CMB power spectrum, the foregrounds, WMAP detector noise or external templates. Compared with the foreground-cleaned internal linear combination map produced by the WMAP team, our map has the advantage of containing less non-CMB power (from foregrounds and detector noise) outside the Galactic plane. The difference is most important on the the angular scale of the first acoustic peak and below, since our cleaned map is at the highest (13') rather than lowest (49') WMAP resolution. We also produce a Wiener filtered CMB map, representing our best guess as to what the CMB sky actually looks like, as well as CMB-free maps at the five WMAP frequencies useful for foreground studies. We argue that our CMB map is clean enough that the lowest multipoles can be measured without any galaxy cut, and obtain a quadrupole value that is slightly less low than that from the cut-sky WMAP team analysis. This can be understood from a map of the CMB quadrupole, which shows much of its power falling within the Galaxy cut region, seemingly coincidentally. Intriguingly, both the quadrupole and the octopole are seen to have power suppressed along a particular spatial axis, which lines up between the two, roughly towards (l,b) \~ (-110,60) in Virgo.

The coming flood of CMB polarization experiments, spurred by the recent detection of CMB polarization by DASI and WMAP, will be confronted by many new analysis tasks specific to polarization. For the analysis of CMB polarization data sets, the devil is truly in the details. With this in mind, we present the details of the data analysis for the POLAR experiment, which recently led to the tightest upper limits on the polarization of the Cosmic Microwave Background Radiation at large angular scales. We discuss the data selection process, mapmaking algorithms, offset removal, and the likelihood analysis which were used to find upper limits on the polarization. Stated using the modern convention for reporting CMB Stokes parameters, these limits are 5.0 \uK on both E-type and B-type polarization at 95% confidence. Finally, we discuss simulations used to test our analysis techniques and to probe the fundamental limitations of the experiment.

Most of the useful information about inflationary gravitational waves and reionization is on large angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, quantify and remove polarized foregrounds. We use the Leiden radio surveys to quantify the polarized synchrotron radiation at large angular scales, which is likely to be the most challenging polarized contaminant for the WMAP satellite. We find that the synchrotron E- and B-contributions are equal to within 10% from 408-820MHz with a hint of E-domination at higher frequencies. We quantify Faraday Rotation & Depolarization effects and show that they cause the synchrotron polarization percentage to drop both towards lower frequencies and towards lower multipoles.

Separation of the E and B components of a microwave background polarization map or a weak lensing map is an essential step in extracting science from it, but when the map covers only part of the sky and/or is pixelized, this decomposition cannot be done perfectly. We present a method for decomposing an arbitrary sky map into a sum of three orthogonal components that we term ``pure E'', ``pure B'' and ``ambiguous''. This method is useful both for providing intuition for experimental design and for analyzing data sets in practice. We show how to find orthonormal bases for all three components in terms of bilaplacian eigenfunctions. The number of ambiguous modes is proportional to the length of the map boundary so fairly round maps are preferred. For real-world data analysis, we present a simple matrix eigenvalue method for calculating nearly pure E and B modes in pixelized maps. We find that the dominant source of leakage between E and B is aliasing of small-scale power caused by the pixelization. This problem can be eliminated by heavily oversampling the map, but is exacerbated by the fact that the E power spectrum is expected to be much larger than the B power spectrum and extremely blue. We found that a factor of 2 to 3 more pixels are needed in a polarization map to achieve the same level of contamination by aliased power than in a temperature map.

We compute upper limits on CMB cross-polarization by cross-correlating the PIQUE and Saskatoon experiments. We also discuss theoretical and practical issues relevant to measuring cross-polarization and illustrate them with simulations of the upcoming BOOMERanG 2002 experiment. We present a method that separates all six polarization power spectra (TT, EE, BB, TE, TB, EB) without any other "leakage" than the familiar EE-BB mixing caused by incomplete sky coverage. Since E and B get mixed, one might expect leakage between TE and TB, between EE and EB and between BB and EB - our method eliminates this by preserving the parity symmetry under which TB and EB are odd and the other four power spectra are even.

We measure the cosmic microwave background (CMB) power spectrum on angular scales l~30-200 (1-6 degrees) from the QMASK map, which combines the data from the QMAP and Saskatoon experiments. Since the accuracy of recent measurements leftward of the first acoustic peak is limited by sample-variance, the large area of the QMASK map (648 square degrees) allows us to place among the sharpest constraints to date in this range, in good agreement with BOOMERanG and (on the largest scales) COBE/DMR. By band-pass-filtering the QMAP and Saskatoon maps, we are able to spatially compare them scale-by-scale to check for beam- and pointing-related systematic errors.

We present a new puzzle involving Galactic microwave emission and attempt to resolve it. On one hand, a cross-correlation analysis of the WHAM H-alpha map with the Tenerife 10 and 15 GHz maps shows that the well-known DIRBE correlated microwave emission cannot be dominated by free-free emission. On the other hand, recent high resolution observations in the 8-10 GHz range with the Green Bank 140 ft telescope by Finkbeiner et al. failed to find the corresponding 8 sigma signal that would be expected in the simplest spinning dust models. So what physical mechanism is causing this ubiquitous dust-correlated emission? We argue for a model predicting that spinning dust is the culprit after all, but that the corresponding small grains are well correlated with the larger grains seen at 100 micron only on large angular scales. In support of this grain segregation model, we find the best spinning dust template to involve higher frequency maps in the range 12-60 micron, where emission from transiently heated small grains is important. Upcoming CMB experiments such as ground-based interferometers, MAP and Planck LFI with high resolution at low frequencies should allow a definitive test of this model.

We present a method for measuring CMB polarization power spectra given incomplete sky coverage and test it with simulated examples such as Boomerang 2001 and MAP. By augmenting the quadratic estimator method with an additional step, we find that the E and B power spectra can be effectively disentangled on angular scales substantially smaller than the width of the sky patch in the narrowest direction. We find that the basic quadratic and maximum-likelihood methods display a unneccesary sensitivity to systematic errors when T-E cross-correlation is involved, and show how this problem can be eliminated at negligible cost in increased error bars. We also test numerically the widely used approximation that sample variance scales inversely with sky coverage, and find it to be an excellent approximation on scales substantially smaller than the sky patch.

We present an upper limit on the polarization of the Cosmic Microwave Background at 7 degree angular scales in the frequency band between 26 and 36 GHz, produced by the POLAR experiment. The campaign produced a map of linear polarization over the R.A. range 112 degrees - 275 degrees at declination 43degrees. The model-independent upper limit on the E-mode polarization component of the CMB at angular scales l = 2 - 20 is 10 microKelvin (95% confidence). The corresponding limit for the B-mode is also 10 microKelvin. Constraining the B-mode power to be zero, the 95% confidence limit on E-mode power alone is 8 microKelvin.

We present a method for comparing and combining maps with different resolutions and beam shapes, and apply it to the Saskatoon, QMAP and COBE/DMR data sets. Although the Saskatoon and QMAP maps detect signal at the 21 sigma and 40 sigma levels, respectively, their difference is consistent with pure noise, placing strong limits on possible systematic errors. In particular, we obtain quantitative upper limits on relative calibration and pointing errors. Splitting the combined data by frequency shows similar consistency between the Ka- and Q-bands, placing limits on foreground contamination. The visual agreement between the maps is equally striking. Our combined QMAP+Saskatoon map, nicknamed QMASK, is publicly available at www.hep.upenn.edu/~xuyz/qmask.html together with its 6495x6495 noise covariance matrix. This thoroughly tested data set covers a large enough area (648 square degrees -- currently the largest degree-scale map available) to allow a statistical comparison with COBE/DMR, showing good agreement.

We present Cosmic Microwave Background (CMB) maps from the Santa Barbara HACME balloon experiment (Staren etal 2000), covering about 1150 square degrees split between two regions in the northern sky, near the stars gamma Ursae Minoris and alpha Leonis, respectively. The FWHM of the beam is about 0.77 degrees in three frequency bands centered on 39, 41 and 43 GHz. The results demonstrate that the thoroughly interconnected scan strategy employed allows efficient removal of 1/f-noise and slightly variable scan-synchronous offsets. The maps display no striping, and the noise correlations are found to be virtually isotropic, decaying on an angular scale around one degree. The noise performance of the experiment resulted in an upper limit on CMB anisotropy. However, our results demonstrate that atmospheric contamination and other systematics resulting from the circular scanning strategy can be accurately controlled, and bodes well for the planned follow-up experiments BEAST and ACE, since they show that even with the overly cautious assumption that 1/f-noise and offsets will be as dominant as for HACME, the problems they pose can be readily overcome with the mapmaking algorithm discussed. Our prewhitened notch-filter algorithm for destriping and offset removal should be useful also for other balloon- and ground-based experiments whose scan strategies involve substantial interleaving.

We quantify the level of foreground contamination in the QMAP Cosmic Microwave Background (CMB) data with two objectives: (a) measuring the level to which the QMAP power spectrum measurements need to be corrected for foregrounds and (b) using this data set to further refine current foreground models. We cross-correlate the QMAP data with a variety of foreground templates. The 30 GHz Ka-band data is found to be significantly correlated with the Haslam 408 MHz and Reich and Reich 1420 MHz synchrotron maps, but not with the Diffuse Infrared Background Experiment (DIRBE) 240, 140 and 100 micron maps or the Wisconsin H-Alpha Mapper (WHAM) survey. The 40 GHz Q-band has no significant template correlations. We discuss the constraints that this places on synchrotron, free-free and dust emission. We also reanalyze the foreground-cleaned Ka-band data and find that the two band power measurements are lowered by 2.3% and 1.3%, respectively.

The recent discovery of dust-correlated diffuse microwave emission has prompted two rival explanations: free-free emission and spinning dust grains. We present new detections of this component at 10 and 15 GHz by the switched-beam Tenerife experiment. The data show a turnover in the spectrum and thereby supports the spinning dust hypothesis. We also present a significant detection of synchrotron radiation at 10 GHz, useful for normalizing foreground contamination of CMB experiments at high-galactic latitudes.

We describe the HEMT Advanced Cosmic Microwave Explorer (HACME), a balloon borne experiment designed to measure sub-degree scale Cosmic Microwave Background anisotropy over hundreds of square degrees, using a unique two dimensional scanning strategy. A spinning flat mirror that is canted relative to its spin axis modulates the direction of beam response in a nearly elliptical path on the sky. The experiment was successfully flown in February of 1996, achieving near laboratory performance for several hours at float altitude. A map free of instrumental systematic effects is produced for a 3.5 hour observation of 630 square degrees, resulting in a flat band power upper limit of (l(l+1)C_l/2 pi)^0.5 < 77 microK at l = 38 (95% confidence). The experiment design, flight operations and data, including atmospheric effects and noise performance, are discussed.

One of the main challenges facing upcoming CMB experiments will be to distinguish the cosmological signal from foreground contamination. We present a comprehensive treatment of this problem and study how foregrounds degrade the accuracy with which the Boomerang, MAP and Planck experiments can measure cosmological parameters. Our foreground model includes not only the normalization, frequency dependence and scale dependence for each physical component, but also variations in frequency dependence across the sky. When estimating how accurately cosmological parameter can be measured, we include the important complication that foreground model parameters (we use about 500) must be simultaneously measured from the data as well. Our results are quite encouraging: despite all these complications, precision measurements of most cosmological parameters are degraded by less than a factor of 2 for our main foreground model and by less than a factor of 5 in our most pessimistic scenario. Parameters measured though large-angle polarization signals suffer more degradation: up to 5 in the main model and 25 in the pessimistic case. The foregrounds that are potentially most damaging and therefore most in need of further study are vibrating dust emission and point sources, especially those in the radio frequencies. It is well-known that E and B polarization contain valuable information about reionization and gravity waves, respectively. However, the cross-correlation between polarized and unpolarized foregrounds also deserves further study, as we find that it carries the bulk of the polarization information about most other cosmological parameters.

We cross-correlate a 19 GHz full sky Cosmic Microwave Background (CMB) survey with other maps to quantify the foreground contribution. Correlations are detected with the Diffuse Infrared Background Experiment (DIRBE) 240, 140 and 100 micron maps at high latitudes (|b|>30degrees), and marginal correlations are detected with the Haslam 408 MHz and the Reich & Reich 1420 MHz synchrotron maps. The former agree well with extrapolations from higher frequencies probed by the COBE DMR and Saskatoon experiments and are consistent with both free-free and rotating dust grain emission.

We present results from the QMAP balloon experiment, which maps the Cosmic Microwave Background (CMB) and probes its angular power spectrum on degree scales. In two separate flights, data were taken in six channels at two frequency bands between 26 to 46 GHz. We describe our method for mapmaking (removal of 1/f-noise and scan-synchronous offsets) and power spectrum estimation, as well as the results of a joint analysis of the data from both flights. This produces a 527 square degree map of the CMB around the North Celestial Pole, allowing a wide variety of systematic cross-checks. The frequency dependence of the fluctuations is consistent with CMB and inconsistent with Galactic foreground emission. The anisotropy is measured in three multipole bands from l~40 to l~200, and the angular power spectrum shows a distinct rise which is consistent with the Saskatoon results.

We report on the first flight of the balloon-borne QMAP experiment. The experiment is designed to make a map of the cosmic microwave background anisotropy on angular scales from 0.7 to several degrees. Using the map we determine the angular power spectrum of the anisotropy in multipole bands from l~40 to l~140. The results are consistent with the Saskatoon (SK) measurements. The frequency spectral index (measured at low l) is consistent with that of CMB and inconsistent with either Galactic synchrotron or free-free emission. The instrument, measurement, analysis of the angular power spectrum, and possible systematic errors are discussed.

We report the results from the second flight of QMAP, an experiment to map the cosmic microwave background near the North Celestial Pole. We present maps of the sky at 31 and 42 GHz as well as a measurement of the angular power spectrum covering the l-range 40-200. Anisotropy is detected at about 20 sigma and is in agreement with previous results at these angular scales. We also report details of the data reduction and analysis techniques which were used for both flights of QMAP.

For high-precision cosmic microwave background (CMB) experiments, contamination from extragalactic point sources is a major concern. It is therefore useful to be able to detect and discard point source contaminated pixels using the map itself. We show that the sensitivity with which this can be done can often be greatly improved (by factors between 2.5 and 18 for the upcoming Planck mission) by a customized hi-pass filtering that suppresses fluctuations due to CMB and diffuse galactic foregrounds. This means that point source contamination will not severely degrade the cleanest Planck channels unless current source count estimates are off by more than an order of magnitude. A catalog of around 40,000 far infra-red sources at 857 GHz may be a useful by-product of Planck.

We cross-correlate the Saskatoon Ka and Q-Band Cosmic Microwave Background (CMB) data with different maps to quantify possible foreground contamination. We detect a marginal correlation (2 sigma) with the Diffuse Infrared Background Experiment (DIRBE) 240, 140 and 100 microm maps, but we find no significant correlation with point sources, with the Haslam 408 MHz map or with the Reich and Reich 1420 MHz map. The rms amplitude of the component correlated with DIRBE is about 20% of the CMB signal. Interpreting this component as free-free emission, this normalization agrees with that of Kogut et al. (1996a; 1996b) and supports the hypothesis that the spatial correlation between dust and warm ionized gas observed on large angular scales persists to smaller angular scales. Subtracting this contribution from the CMB data reduces the normalization of the Saskatoon power spectrum by only a few percent.

We present a new data analysis method to study rectangular T^3 ``small universes" with one or two of its dimensions significantly smaller than the present horizon (which we refer to as T^1- and T^2-models, respectively). We find that the 4 year COBE/DMR data set a lower limit on the smallest cell size for T^1- and T^2-models of 3000h^{-1} Mpc, at 95% confidence, for a scale invariant power spectrum (n=1).

The cosmic microwave background (CMB) is a unique probe of cosmological parameters and conditions. There is a connection between anisotropy in the CMB and the topology of the Universe. Adopting a universe with the topology of a 3-Torus, or a universe where only harmonics of the fundamental mode are allowed, and using 2-years of COBE/DMR data, we obtain constraints on the topology of the Universe. Previous work constrained the topology using the slope information and the correlation function of the CMB. We obtain more accurate results by using all multipole moments, avoiding approximations by computing their full covariance matrix. We obtain the best fit for a cubic toroidal universe of scale 7200h^{-1} Mpc for n=1. The data set a lower limit on the cell size of 4320h^{-1} Mpc at 95% confidence and 5880h^{-1} Mpc at 68% confidence. These results show that the most probable cell size would be around 1.2 times larger than the horizon scale, implying that the 3-Torus topology is no longer an interesting cosmological model.

Although the cubic T^3 "small universe" has been ruled out by COBE/DMR results as an interesting cosmological model, we still have the possibility of living in a universe with a more anisotropic topology such as a rectangular T^3 "small universe" with one or two of its dimensions significantly smaller than the present horizon (which we refer to as T^1- and T^2-models, respectively). In order to rule out these anisotropic topologies as well, we apply a new data analysis method that searches for the specific kind of symmetries that these models should produce. We find that the 4 year COBE/DMR data set a lower limit on the smallest cell size for T^1- and T^2-models of 3000h^{-1} Mpc, at 95% confidence, for a scale invariant power spectrum (n=1). These results imply that ALL toroidal universes (cubes and rectangles) are ruled out as interesting cosmological models.

We cross-correlate the Saskatoon Q-Band data with different spatial template maps to quantify possible foreground contamination. We detect a correlation with the Diffuse Infrared Background Experiment (DIRBE) 100 microm map, which we interpret as being due to Galactic free-free emission. Subtracting this foreground power reduces the Saskatoon normalization of the Cosmic Microwave Background (CMB) power spectrum by roughly 2%.

We present a Wiener filtered map of the Cosmic Microwave Background (CMB) fluctuations in a disk with 15 degree diameter, centered at the North Celestial Pole. The map is based on the 1993-1995 data from the Saskatoon (SK) experiment, with an angular resolution around 1 degree in the frequency range 27.6-44.1 GHz. The signal-to-noise ratio in the map is of order two, and some individual hot and cold spots are significant at the 5 sigma level. The spatial features are found to be consistent from year to year, reenforcing the conclusion that the SK results are not dominated by residual atmospheric contamination or other non-celestial signals.