Herschel and Planck are surveying the sky at unprecedented angular scales and sensitivities over large areas. But both experiments are limited by source confusion in the submillimeter. The high confusion noise in particular restricts the study of the clustering properties of the sources that dominate the cosmic infrared background. At these wavelengths, it is more appropriate to consider the statistics of the unresolved component. In particular, high clustering will contribute in excess of Poisson noise in the power spectra of CIB anisotropies. These power spectra contain contributions from sources at all redshift. We show how the stacking technique can be used to separate the different redshift contributions to the power spectra. We use simulations of CIB representative of realistic Spitzer, Herschel, Planck, and SCUBA-2 observations. We stack the 24um sources in longer wavelengths maps to measure mean colors per redshift and flux bins. The information retrieved on the mean spectral energy distribution obtained with the stacking technique is then used to clean the maps, in particular to remove the contribution of low-redshift undetected sources to the anisotropies. Using the stacking, we measure the mean flux of populations 4 to 6 times fainter than the total noise at 350um at redshifts z=1 and z=2, respectively, and as faint as 6 to 10 times fainter than the total noise at 850um at the same redshifts. In the deep Spitzer fields, the detected 24um sources up to z~2 contribute significantly to the submillimeter anisotropies. We show that the method provides excellent (using COSMOS 24um data) to good (using SWIRE 24um data) removal of the z<2 (COSMOS) and z<1 (SWIRE) anisotropies. Using this cleaning method, we then hope to have a set of large maps dominated by high redshift galaxies for galaxy evolution study (e.g., clustering, luminosity density).

The coming Planck and Herschel missions will survey the sky at unprecedented angular scales and sensitivities. Simulations are needed for better interpretating the results of the surveys and for testing new methods of, e.g., source extraction and component separation. We present new simulations of the infrared and submillimeter cosmic background, including the correlation between infrared galaxies. The simulations were used to quantify the source-detection thresholds for Herschel/SPIRE and Planck/HFI, as well as to study the detectability of the cosmic infrared background correlated fluctuations. The simulations are based on an empirical model of IR galaxy evolution. For the correlations, we only included the linear clustering, assuming that infrared galaxies are biased tracers of the dark-matter fluctuation density field. We used the simulations with different bias parameters to predict the confusion noise for Herschel/SPIRE and Planck/HFI and the completeness levels. We also discuss the detectability of the linear clustering in Planck/HFI power spectra, including the foreground and backgrounds components. Simulated maps and catalogs are publicly available online at http://www.ias.u-psud.fr/irgalaxies/simulations.php

We report the detection of correlated anisotropies in the Cosmic Far-Infrared Background at 160 microns. We measure the power spectrum in the Spitzer/SWIRE Lockman Hole field. It reveals unambiguously a strong excess above cirrus and Poisson contributions, at spatial scales between 5 and 30 arcminutes, interpreted as the signature of infrared galaxy clustering. Using our model of infrared galaxy evolution we derive a linear bias b=1.74 \pm 0.16. It is a factor 2 higher than the bias measured for the local IRAS galaxies. Our model indicates that galaxies dominating the 160 microns correlated anisotropies are at z~1. This implies that infrared galaxies at high redshifts are biased tracers of mass, unlike in the local Universe.

Aims. We quantify the contributions of 24um galaxies to the Far-Infrared (FIR) Background at 70 and 160um. We provide new estimates of the Cosmic Infrared Background (CIB), and compare it with the Cosmic Optical Background (COB). Methods. Using Spitzer data at 24, 70 and 160um in three deep fields, we stacked more than 19000 MIPS 24um sources with S24>60uJy at 70 and 160um, and measured the resulting FIR flux densities. Results. This method allows a gain up to one order of magnitude in depth in the FIR. We find that the Mid-Infrared (MIR) 24um selected sources contribute to more than 70% of the CIB at 70 and 160um. This is the first direct measurement of the contribution of MIR-selected galaxies to the FIR CIB. Galaxies contributing the most to the total CIB are thus z~1 luminous infrared galaxies, which have intermediate stellar masses. We estimate that the CIB will be resolved at 0.9 mJy at 70 and 3 mJy at 160um. By combining the extrapolation of the 24um source counts below 60uJy, with 160/24 and 70/24 colors as measured with the stacking analysis, we obtain lower limits of 7.1+/-1.0 and 13.4+/-1.7 nW/m2/sr for the CIB at 70 and 160um, respectively. Conclusions. The MIPS surveys have resolved more than three quarters of the MIR and FIR CIB. By carefully integrating the Extragalactic Background Light (EBL) SED, we also find that the CIB has the same brightness as the COB, around 24 nW/m2/sr. The EBL is produced on average by 115 infrared photons for one visible photon. Finally, the galaxy formation and evolution processes emitted a brightness equivalent to 5% of the primordial electromagnetic background (CMB).