4.3. Source confusion

Source confusion is an additional noise factor closely related to the astronomical background, described in Section 4.1, “Background radiation”. The sensitivity limit due to confusion is determined by the telescope aperture, observation wavelength and the position on the sky. The sensitivity cannot be improved by increasing the integration time after reaching the confusion limit. The most important contributions to source confusion are:

The confusion noise is usually defined as the (stochastic) fluctuations of the background sky brightness below which sources cannot be detected individually. In addition to the diffuse Galactic foreground cirrus component, these fluctuations are caused by intrinsically discrete extragalactic sources in the beam. Due to the limited telescope diameter compared to the wavelength, these fluctuations play an important, if not dominant, role in the total noise budget in extragalactic surveys carried out in the MIR, FIR and sub-mm range. Moreover, the noise due to extragalactic sources depends strongly on the shape of the source counts at a given wavelength.

There are two different criteria to derive the confusion noise, and thus the detectability of a point-like or compact source:

Generally, we should compare the confusion noise derived from both criteria, in order not to underestimate it artificially. The confusion noise, σc, and confusion limit, Slim are defined as follows:

σc² = ƒ²(θ,φ)dθdφ0S_limS²(dN/dS)dS

where ƒ(θ,φ) is the instrumental 2D beam profile, that can be approximated by a Gaussian profile with the same FWHM as the expected PSF, or by an Airy function, S is the source flux density (in Jy) and dN/dS is the differential source number counts (in Jy-1sr-1).

Then, the total noise is computed by adding in quadrature the different noise contributions, in this case the photon (and instrumental) noise and the confusion noise, i.e. σtotal = (σp2 + σc2)½

The photometric criterion is defined by choosing the S/N ratio qphot between the faintest source (of flux Slim and the noise σc due to fluctuations from beam to beam caused by sources fainter than Slim, as given by the implicit equation:

qphot = Slim / σc(Slim)

q is usually chosen between 3 and 5, depending on the specific objectives.

The source density criterion is defined by setting the minimum degree of completeness of the detection of sources above the limiting flux Slim, which is driven by the fraction of sources lost in the detection process due to a nearest neighbour source with flux above Slim too close to be separated given an instrumental beam size. For a given Poissonian source density N(>S), the probability P of finding a nearest neighbour with S Slim at a distance closer than the minimum angular separation θmin is given by:

P(< θmin) = 1 - exp(-πNθ²min)

An acceptable probability limit is Pmax = 0.1. The minimum distance is usually parameterised using the FWHM of the beam profile θmin = kθFWHM, and 0.8 k 1. Fixing the probability we obtain the corresponding "source density criterion" limiting density of sources:

NSDC = -ln(1-P(<θmin)) / πNk²θ²FWHM

The instrumental beam area, is given by Ω 1.14θ²FWHM. Therefore, for P = 0.1 and k = 0.8, the density is 1/16.7 sources/beam. The limiting source flux, SSDC is thus determined by using existing number counts results and a suitable model for infrared galaxy evolution extrapolating the data to the appropriate wavelengths and (faint) flux levels. The confusion noise, σSDC is computed using the same relation as for the photometric criterion, as the S/N ratio qSDC = SSDC / σSDC.

The Herschel confusion noise levels due to extragalactic sources in the different instruments/bands are being computed from deep maps on 'blank' fields. The preliminary results ([RD11] and [RD12]) are shown in Table 4.1 along with the values predicted by by Lagache et al. 2003 ([RD8]) using number counts derived from a phenomenological model based on template spectra of starburst and normal galaxies, and on the local infrared luminosity function. This model has been found to be in very good overall agreement with ISOCAM at 15 μm, IRAS at 60 and 170 μm and SCUBA at 850 μm (see references within [RD8]). Confusion level predictions for Herschel/PACS have been also computed by Dole et al. 2004 ([RD9]) based on recent Spitzer/MIPS number counts from Papovich et al. 2004 ([RD10]) shown in Figure 4.4. They obtain SSDC(70 μm) = 0.16 mJy and SSDC(160 μm) = 10.0 mJy. Please refer to the specific instruments' Observers' Manual for an up-to-date information regarding confusion noise.

Table 4.1. PACS and SPIRE measured confusion noise, compared to predictions computed according to photometric and source density criteria. From [RD9]. I.

 σobserved (mJy) σ (mJy)Slim (mJy)
PACS 75 μm N/Aqphot = 5.02.26 × 10-31.12 × 10-2
  qSDC = 8.91.42 × 10-21.26 × 10-1
PACS 110 μm0.1-0.3qphot = 5.01.98 × 10-21.0 × 10-1
  qSDC = 8.71.02 × 10-18.91 × 10-1
PACS 170 μm0.7-0.9qphot = 5.03.97 × 10-12.00
  qSDC = 7.139.93 × 10-17.08
SPIRE 250 μm 5.8qphot = 5.02.5112.6
  qSDC = 5.22.7014.1
SPIRE 360 μm6.3qphot = 5.04.422.4
  qSDC = 3.63.5212.6
SPIRE 550 μm6.8qphot = 5.03.6917.8
  qSDC = 2.53.187.94

Cumulative (left) and differential (right) 24 μm number counts from RD10. The differential counts have been normalised to an Euclidean slope, dN/dSν ∼ Sν-2.5. The curves show predictions from different recent models, including that from Lagache et al. 2003.

Figure 4.4. Cumulative (left) and differential (right) 24 μm number counts from [RD10]. The differential counts have been normalised to an Euclidean slope, dN/dSν Sν-2.5. The curves show predictions from different recent models, including that from Lagache et al. 2003.