New determinations are presented of the cosmic infrared background monopole brightness in the Planck HFI bands from 100 GHz to 857 GHz. Planck was not designed to measure the monopole component of sky brightness, so cross-correlation of the 2015 HFI maps with COBE/FIRAS data is used to recalibrate the zero level of the HFI maps. For the HFI 545 and 857 GHz maps, the brightness scale is also recalibrated. Correlation of the recalibrated HFI maps with a linear combination of Galactic H I and H alpha data is used to separate the Galactic foreground emission and determine the cosmic infrared background brightness in each of the HFI bands. We obtain CIB values of 0.007 +- 0.014, 0.010 +- 0.019, 0.060 +- 0.023, 0.149 +- 0.017, 0.371 +- 0.018, and 0.576 +- 0.034 MJy/sr at 100, 143, 217, 353, 545, and 857 GHz, respectively. The estimated uncertainties for the 353 to 857 GHz bands are about 3 to 6 times smaller than those of previous direct CIB determinations at these frequencies. Our results are compared with integrated source brightness results from selected recent submillimeter and millimeter wavelength imaging surveys.

This is an overview of the data products and other resources available through NASA's LAMBDA site https://lambda.gsfc.nasa.gov/. An up-to-date version of this document, along with code tools actively maintained and developed by LAMBDA staff, can be found on the LAMBDA GitHub page at https://github.com/nasa-lambda/lambda_overview. New data products and other updates are announced on LAMBDA's twitter account at https://twitter.com/NASA_LAMBDA. If you have questions or suggestions relating to LAMBDA, or are interested in joining a LAMBDA advisory group, please contact us using the form here: https://lambda.gsfc.nasa.gov/contact/contact.cfm.

While the upcoming telescopes will reveal correspondingly fainter, more distant galaxies, a question will persist: what more is there that these telescopes cannot see? One answer is the source-subtracted Cosmic Infrared Background (CIB). The CIB is comprised of the collective light from all sources remaining after known, resolved sources are accounted for. Ever-more-sensitive surveys will identify the brightest of these, allowing them to be removed, and - like peeling layers off an onion - reveal deeper layers of the CIB. In this way it is possible to measure the contributions from populations not accessible to direct telescopic observation. Measurement of fluctuations in the source-subtracted CIB, i.e., the spatial power spectrum of the CIB after subtracting resolved sources, provides a robust means of characterizing its faint, and potentially new, populations. Studies over the past 15 years have revealed source-subtracted CIB fluctuations on scales out to ~100' which cannot be explained by extrapolating from known galaxy populations. Moreover, they appear highly coherent with the unresolved Cosmic X-ray Background, hinting at a significant population of accreting black holes among the CIB sources. Characterizing the source-subtracted CIB with high accuracy, and thereby constraining the nature of the new populations, is feasible with upcoming instruments and would produce critically important cosmological information in the next decade. New coextensive deep and wide-area near-infrared, X-ray, and microwave surveys will bring decisive opportunities to examine, with high fidelity, the spatial spectrum and origin of the CIB fluctuations and their cross-correlations with cosmic microwave and X-ray backgrounds, and determine the formation epochs and the nature of the new sources (stellar nucleosynthetic or accreting black holes).

Accurate modeling of the spectrum of thermal dust emission at millimeter wavelengths is important for improving the accuracy of foreground subtraction for CMB measurements, for improving the accuracy with which the contributions of different foreground emission components can be determined, and for improving our understanding of dust composition and dust physics. We fit four models of dust emission to high Galactic latitude COBE/FIRAS and COBE/DIRBE observations from 3 millimeters to 100 microns and compare the quality of the fits. We consider the two-level systems model because it provides a physically motivated explanation for the observed long wavelength flattening of the dust spectrum and the anticorrelation between emissivity index and dust temperature. We consider the model of Finkbeiner, Davis, and Schlegel because it has been widely used for CMB studies, and the generalized version of this model recently applied to Planck data by Meisner and Finkbeiner. For comparison we have also fit a phenomenological model consisting of the sum of two graybody components. We find that the two-graybody model gives the best fit and the FDS model gives a significantly poorer fit than the other models. The Meisner and Finkbeiner model and the two-level systems model remain viable for use in Galactic foreground subtraction, but the FIRAS data do not have sufficient signal-to-noise ratio to provide a strong test of the predicted spectrum at millimeter wavelengths.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

[Abridged] We present updated estimates of Galactic foreground emission using seven years of WMAP data. Using the power spectrum of differences between multi-frequency template-cleaned maps, we find no evidence for foreground contamination outside of the updated (KQ85y7) foreground mask. We place a 15 microKelvin upper bound on rms foreground contamination in the cleaned maps used for cosmological analysis. We find no indication in the polarization data of an extra "haze" of hard synchrotron emission from energetic electrons near the Galactic center. We provide an updated map of the cosmic microwave background (CMB) using the internal linear combination (ILC) method, updated foreground masks, and updates to point source catalogs with 62 newly detected sources. Also new are tests of the Markov chain Monte Carlo (MCMC) foreground fitting procedure against systematics in the time-stream data, and tests against the observed beam asymmetry. Within a few degrees of the Galactic plane, WMAP total intensity data show a rapidly steepening spectrum from 20-40 GHz, which may be due to emission from spinning dust grains, steepening synchrotron, or other effects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepening synchrotron models fit the combination of WMAP and 408 MHz data equally well. ARCADE data appear inconsistent with the steepening synchrotron model, and consistent with the spinning dust model, though some discrepancies remain regarding the relative strength of spinning dust emission. More high-resolution data in the 10-40 GHz range would shed much light on these issues.

(Abridged) The 7-year WMAP data and improved astrophysical data rigorously test the standard cosmological model and its extensions. By combining WMAP with the latest distance measurements from BAO and H0 measurement, we determine the parameters of the simplest LCDM model. The power-law index of the primordial power spectrum is n_s=0.968+-0.012, a measurement that excludes the scale-invariant spectrum by 99.5%CL. The other parameters are also improved from the 5-year results. Notable examples of improved parameters are the total mass of neutrinos, sum(m_nu)<0.58eV, and the effective number of neutrino species, N_eff=4.34+0.86-0.88. We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis. We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z=1090 and the dominance of adiabatic scalar fluctuations. With the 7-year TB power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved to Delta(alpha)=-1.1+-1.4(stat)+-1.5(syst) degrees. We report significant detections of the SZ effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data. However, it is a factor of 0.5 to 0.7 times the predictions from "universal profile" of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically-expected SZ power spectrum recently measured by the South Pole Telescope collaboration.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

We describe a sampling method to estimate the polarized CMB signal from observed maps of the sky. We use a Metropolis-within-Gibbs algorithm to estimate the polarized CMB map, containing Q and U Stokes parameters at each pixel, and its covariance matrix. These can be used as inputs for cosmological analyses. The polarized sky signal is parameterized as the sum of three components: CMB, synchrotron emission, and thermal dust emission. The polarized Galactic components are modeled with spatially varying power law spectral indices for the synchrotron, and a fixed power law for the dust, and their component maps are estimated as by-products. We apply the method to simulated low resolution maps with pixels of side 7.2 degrees, using diagonal and full noise realizations drawn from the WMAP noise matrices. The CMB maps are recovered with goodness of fit consistent with errors. Computing the likelihood of the E-mode power in the maps as a function of optical depth to reionization, tau, for fixed temperature anisotropy power, we recover tau=0.091+-0.019 for a simulation with input tau=0.1, and mean tau=0.098 averaged over 10 simulations. A `null' simulation with no polarized CMB signal has maximum likelihood consistent with tau=0. The method is applied to the five-year WMAP data, using the K, Ka, Q and V channels. We find tau=0.090+-0.019, compared to tau=0.086+-0.016 from the template-cleaned maps used in the primary WMAP analysis. The synchrotron spectral index, beta, averaged over high signal-to-noise pixels with standard deviation sigma(beta)<0.25, but excluding ~6% of the sky masked in the Galactic plane, is -3.03+-0.04. This estimate does not vary significantly with Galactic latitude, although includes an informative prior.

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced chi-squared of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

We present the list of point sources found in the WMAP 5-year maps. The technique used in the first-year and three-year analysis now finds 390 point sources, and the five-year source catalog is complete for regions of the sky away from the galactic plane to a 2 Jy limit, with SNR > 4.7 in all bands in the least covered parts of the sky. The noise at high frequencies is still mainly radiometer noise, but at low frequencies the CMB anisotropy is the largest uncertainty. A separate search of CMB-free V-W maps finds 99 sources of which all but one can be identified with known radio sources. The sources seen by WMAP are not strongly polarized. Many of the WMAP sources show significant variability from year to year, with more than a 2:1 range between the minimum and maximum fluxes.

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that Ka band data (33 GHz) is suitable for use in cosmological analysis, after foreground cleaning. This significantly reduces the overall polarization uncertainty. With the 5 year WMAP data, we detect no convincing deviations from the minimal 6-parameter LCDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations, we find (68% CL uncertainties): Omega_bh^2 = 0.02267 \pm 0.00059, Omega_ch^2 = 0.1131 \pm 0.0034, Omega_Lambda = 0.726 \pm 0.015, n_s = 0.960 \pm 0.013, tau = 0.084 \pm 0.016, and Delta_R^2 = (2.445 \pm 0.096) x 10^-9. From these we derive: sigma_8 = 0.812 \pm 0.026, H_0 = 70.5 \pm 1.3 km/s/Mpc, z_{reion} = 10.9 \pm 1.4, and t_0 = 13.72 \pm 0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.22 (95% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and spatial curvature: -0.14 < 1+w < 0.12 and -0.0179 < Omega_k < 0.0081 (both 95% CL). The number of relativistic degrees of freedom (e.g. neutrinos) is found to be N_{eff} = 4.4 \pm 1.5, consistent with the standard value of 3.04. Models with N_{eff} = 0 are disfavored at >99.5% confidence.

Cosmology and other scientific results from the WMAP mission require an accurate knowledge of the beam patterns in flight. While the degree of beam knowledge for the WMAP one-year and three-year results was unprecedented for a CMB experiment, we have significantly improved the beam determination as part of the five-year data release. Physical optics fits are done on both the A and the B sides for the first time. The cutoff scale of the fitted distortions on the primary mirror is reduced by a factor of ~2 from previous analyses. These changes enable an improvement in the hybridization of Jupiter data with beam models, which is optimized with respect to error in the main beam solid angle. An increase in main-beam solid angle of ~1% is found for the V2 and W1-W4 differencing assemblies. Although the five-year results are statistically consistent with previous ones, the errors in the five-year beam transfer functions are reduced by a factor of ~2 as compared to the three-year analysis. We present radiometry of the planet Jupiter as a test of the beam consistency and as a calibration standard; for an individual differencing assembly, errors in the measured disk temperature are ~0.5%.

We present the temperature and polarization angular power spectra of the cosmic microwave background (CMB) derived from the first 5 years of WMAP data. The 5-year temperature (TT) spectrum is cosmic variance limited up to multipole l=530, and individual l-modes have S/N>1 for l<920. The best fitting six-parameter LambdaCDM model has a reduced chi^2 for l=33-1000 of chi^2/nu=1.06, with a probability to exceed of 9.3%. There is now significantly improved data near the third peak which leads to improved cosmological constraints. The temperature-polarization correlation (TE) is seen with high significance. After accounting for foreground emission, the low-l reionization feature in the EE power spectrum is preferred by \Delta\chi^2=19.6 for optical depth tau=0.089 by the EE data alone, and is now largely cosmic variance limited for l=2-6. There is no evidence for cosmic signal in the BB, TB, or EB spectra after accounting for foreground emission. We find that, when averaged over l=2-6, l(l+1)C^{BB}_l/2\pi < 0.15 uK^2 (95% CL).

Determination of the cosmic infrared background (CIB) at far infrared wavelengths using COBE/DIRBE data is limited by the accuracy to which foreground interplanetary and Galactic dust emission can be modeled and subtracted. Previous determinations of the far infrared CIB (e.g., Hauser et al. 1998) were based on the detection of residual isotropic emission in skymaps from which the emission from interplanetary dust and the neutral interstellar medium were removed. In this paper we use the Wisconsin H-alpha Mapper (WHAM) Northern Sky Survey as a tracer of the ionized medium to examine the effect of this foreground component on determination of the CIB. We decompose the DIRBE far infrared data for five high Galactic latitude regions into H I and H-alpha correlated components and a residual component. We find the H-alpha correlated component to be consistent with zero for each region, and we find that addition of an H-alpha correlated component in modeling the foreground emission has negligible effect on derived CIB results. Our CIB detections and 2 sigma upper limits are essentially the same as those derived by Hauser et al. and are given by nu I_nu (nW m-2 sr-1) < 75, < 32, 25 +- 8, and 13 +- 3 at 60, 100, 140, and 240 microns, respectively. Our residuals have not been subjected to a detailed anisotropy test, so our CIB results do not supersede those of Hauser et al. We derive upper limits on the 100 micron emissivity of the ionized medium that are typically about 40% of the 100 micron emissivity of the neutral atomic medium. This low value may be caused in part by a lower dust-to-gas mass ratio in the ionized medium than in the neutral medium, and in part by a shortcoming of using H-alpha intensity as a tracer of far infrared emission.

We present a full-sky model of polarized Galactic microwave emission based on three years of observations by the Wilkinson Microwave Anisotropy Probe (WMAP) at frequencies from 23 to 94 GHz. The model compares maps of the Stokes Q and U components from each of the 5 WMAP frequency bands in order to separate synchrotron from dust emission, taking into account the spatial and frequency dependence of the synchrotron and dust components. This simple two-component model of the interstellar medium accounts for at least 97% of the polarized emission in the WMAP maps of the microwave sky. Synchrotron emission dominates the polarized foregrounds at frequencies below 50 GHz, and is comparable to the dust contribution at 65 GHz. The spectral index of the synchrotron component, derived solely from polarization data, is -3.2 averaged over the full sky, with a modestly flatter index on the Galactic plane. The synchrotron emission has mean polarization fraction 2--4% in the Galactic plane and rising to over 20% at high latitude, with prominent features such as the North Galactic Spur more polarized than the diffuse component. Thermal dust emission has polarization fraction 1% near the Galactic center, rising to 6% at the anti-center. Diffuse emission from high-latitude dust is also polarized with mean fractional polarization 0.036 +/- 0.011.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

The WMAP satellite has completed 3 years of observations of the cosmic microwave background radiation. The 3-year data products include several sets of full sky maps of the Stokes I, Q and U parameters in 5 frequency bands, spanning 23 to 94 GHz, and supporting items, such as beam window functions and noise covariance matrices. The processing used to produce the current sky maps and supporting products represents a significant advancement over the first year analysis, and is described herein. Improvements to the pointing reconstruction, radiometer gain modeling, window function determination and radiometer spectral noise parametrization are presented. A detailed description of the updated data processing that produces maximum likelihood sky map estimates is presented, along with the methods used to produce reduced resolution maps and corresponding noise covariance matrices. Finally two methods used to evaluate the noise of the full resolution sky maps are presented along with several representative year-to-year null tests, demonstrating that sky maps produced from data from different observational epochs are consistent.

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

We present full sky microwave maps in five bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 2<l<900 anisotropy power spectrum is cosmic variance limited for l<354 with a signal-to-noise ratio >1 per mode to l=658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is 0.17 +/- 0.04, which implies a reionization epoch of 180+220-80 Myr (95% CL) after the Big Bang at a redshift of 20+10-9 (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. The age of the best-fit universe is 13.7 +/- 0.2 Gyr old. Decoupling was 379+8-7 kyr after the Big Bang at a redshift of 1089 +/- 1. The thickness of the decoupling surface was dz=195 +/- 2. The matter density is Omega_m h^2 = 0.135 +0.008 -0.009, the baryon density is Omega_b h^2 = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega_tot = 1.02 +/- 0.02. The spectral index of scalar fluctuations is fit as n_s = 0.93 +/- 0.03 at wavenumber k_0 = 0.05 Mpc^-1, with a running index slope of dn_s/d ln k = -0.031 +0.016 -0.018 in the best-fit model. This flat universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy equation of state is limited to w<-0.78 (95% CL). Inflation theory is supported with n_s~1, Omega_tot~1, Gaussian random phases of the CMB anisotropy, and superhorizon fluctuations. An admixture of isocurvature modes does not improve the fit. The tensor-to-scalar ratio is r(k_0=0.002 Mpc^-1)<0.90 (95% CL).

Full sky maps are made in five microwave frequency bands to separate the temperature anisotropy of the CMB from foreground emission. We define masks that excise regions of high foreground emission. The effectiveness of template fits to remove foreground emission from the WMAP data is examined. These efforts result in a CMB map with minimal contamination and a demonstration that the WMAP CMB power spectrum is insensitive to residual foreground emission. We construct a model of the Galactic emission components. We find that the Milky Way resembles other normal spiral galaxies between 408 MHz and 23 GHz, with a synchrotron spectral index that is flattest (beta ~ -2.5) near star-forming regions, especially in the plane, and steepest (beta ~ -3) in the halo. The significant synchrotron index steepening out of the plane suggests a diffusion process in which the halo electrons are trapped in the Galactic potential long enough to suffer synchrotron and inverse Compton energy losses and hence a spectral steepening. The synchrotron index is steeper in the WMAP bands than in lower frequency radio surveys, with a spectral break near 20 GHz to beta < -3. The modeled thermal dust spectral index is also steep in the WMAP bands, with beta ~ 2.2. Microwave and H alpha measurements of the ionized gas agree. Spinning dust emission is limited to < ~5% of the Ka-band foreground emission. A catalog of 208 point sources is presented. Derived source counts suggest a contribution to the anisotropy power from unresolved sources of (15.0 +- 1.4) 10^{-3} microK^2 sr at Q-band and negligible levels at V-band and W-band.

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of 1/f noise. We infer and correct for a small ~1% transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

The COBE Diffuse Infrared Background Experiment (DIRBE) was designed to search for the cosmic infrared background (CIB) radiation. Scattered light and thermal emission from the interplanetary dust (IPD) are major contributors to the diffuse sky brightness at most infrared wavelengths. Accurate removal of this zodiacal light foreground is a necessary step toward a direct measurement of the CIB. The zodiacal light foreground contribution in each of the 10 DIRBE wavelength bands ranging from 1.25 to 240 microns is distinguished by its apparent seasonal variation over the whole sky. This contribution has been extracted by fitting the brightness calculated from a parameterized physical model to the time variation of the all-sky DIRBE measurements over 10 months of observations. The model brightness is evaluated as the integral along the line of sight of the product of a source function and a three-dimensional dust density distribution function. The dust density distribution is composed of multiple components: a smooth cloud, three asteroidal dust bands, and a circumsolar ring near 1 A.U. By using a directly measurable quantity which relates only to the IPD cloud, we exclude other contributors to the sky brightness from the IPD model. Using the IPD model described here, high-quality maps of the infrared sky with the zodiacal foreground removed have been generated. Imperfections in the model reveal themselves as low-level systematic artifacts in the residual maps which correlate with components of the IPD. The most evident of these artifacts are located near the ecliptic plane in the mid-infrared, and are less than 2% of the zodiacal foreground brightness. Uncertainties associated with the model are discussed, including implications for the CIB search.

The DIRBE on the COBE spacecraft was designed primarily to conduct systematic search for an isotropic CIB in ten photometric bands from 1.25 to 240 microns. The results of that search are presented here. Conservative limits on the CIB are obtained from the minimum observed brightness in all-sky maps at each wavelength, with the faintest limits in the DIRBE spectral range being at 3.5 microns (\nu I_\nu < 64 nW/m^2/sr, 95% CL) and at 240 microns (\nu I_\nu < 28 nW/m^2/sr, 95% CL). The bright foregrounds from interplanetary dust scattering and emission, stars, and interstellar dust emission are the principal impediments to the DIRBE measurements of the CIB. These foregrounds have been modeled and removed from the sky maps. Assessment of the random and systematic uncertainties in the residuals and tests for isotropy show that only the 140 and 240 microns data provide candidate detections of the CIB. The residuals and their uncertainties provide CIB upper limits more restrictive than the dark sky limits at wavelengths from 1.25 to 100 microns. No plausible solar system or Galactic source of the observed 140 and 240 microns residuals can be identified, leading to the conclusion that the CIB has been detected at levels of \nu I_\nu = 25+-7 and 14+-3 nW/m^2/sr at 140 and 240 microns respectively. The integrated energy from 140 to 240 microns, 10.3 nW/m^2/sr, is about twice the integrated optical light from the galaxies in the Hubble Deep Field, suggesting that star formation might have been heavily enshrouded by dust at high redshift. The detections and upper limits reported here provide new constraints on models of the history of energy-releasing processes and dust production since the decoupling of the cosmic microwave background from matter.

The Cosmic Infrared Background (CIB) is hidden behind veils of foreground emission from our own solar system and Galaxy. This paper describes procedures for removing the Galactic IR emission from the 1.25 - 240 micron COBE DIRBE maps as steps toward the ultimate goal of detecting the CIB. The Galactic emission models are carefully chosen and constructed so that the isotropic CIB is completely retained in the residual sky maps. We start with DIRBE data from which the scattered light and thermal emission of the interplanetary dust (IPD) cloud have already been removed. Locations affected by the emission from bright compact and stellar sources are excluded from the analysis. The unresolved emission of faint stars at near- and mid-IR wavelengths is represented by a model based on Galactic source counts. The 100 micron DIRBE observations are used as the spatial template for the interstellar medium (ISM) emission at high latitudes. Correlation of the 100 micron data with H I column density allows us to isolate the component of the observed emission that is associated with the ISM. Limits are established on the far-IR emissivity of the diffuse ionized medium, which indicate a lower emissivity per H nucleus than in the neutral medium. At 240 micron, we find that adding a second spatial template to the ISM model can greatly improve the accuracy of the model at low latitudes. The crucial product of this analysis is a set of all-sky IR maps from which the Galactic (and IPD) emission has been removed. We discuss systematic uncertainties and potential errors in the foreground subtraction process that may have an impact on studies seeking to detect the CIB in the residual maps.

In this paper we examine the cosmological constraints of the recent DIRBE and FIRAS detection of the extragalactic background light between 125-5000 microns on the metal and star formation histories of the universe.

We have created a map of the large-scale infrared surface brightness in excess of that associated with the atomic interstellar medium, using region-by-region correlations between the far-infrared and 21-cm line surface brightness. Our study updates and extends a previous attempt with IRAS and Berkeley/Parkes H I surveys. The far-infrared observations used here are from DIRBE, which extends far-infrared wavelength coverage to 240 um, so that we are reliably sampling the emission of large, thermal-equilibrium grains that dominate the dust mass. The H I data are from the combined Leiden-Dwingeloo and Parkes 21-cm line surveys. Using the maps of excess infrared emission at 100, 140, and 240 um, we created an atlas and identified the coherent structures. These infrared excess clouds can be caused both by dust that is warmer than average, or by dust associated with gas other than the atomic interstellar medium. We find very few warm clouds, such as the H II region around Spica. The majority of the infrared excess clouds are colder than the average atomic interstellar medium. These clouds are peaks of column density, and their excess infrared emission is due to dust associated with molecular gas. We identify essentially all known high-latitude molecular clouds in the infrared excess maps, and further identify a sample of new clouds with similar infrared properties. The infrared excess was correlated with CO line brightness, allowing us to measure the ratio of N(H2)/W(CO) for high-latitude clouds. The atlas of infrared excess may be a useful guide to regions of relatively high column density, which might cause high extinction toward extragalactic objects at optical and ultraviolet wavelengths and confusion for cosmic infrared and microwave backgrounds.