Final analysis of ELAIS 15-μm observations: method, reduction and catalogue

We present the final analysis of the European Large Area Infrared Space Observatory ( ISO ) Survey (ELAIS) 15-μm observations, carried out with the ISO Camera (ISOCAM) instrument on board the ISO .
The data-reduction method, known as the Lari Method, is based on a mathematical model of the behaviour of the detector and was specifically designed for the detection of faint sources in ISOCAM/ISO Photopolarimeter (ISOPHOT) data. The method is fully interactive and leads to very reliable and complete source lists.
The resulting catalogue includes 1923 sources detected with signal-to-noise ratio of  > 5  in the 0.5–100 mJy flux range and over an area of 10.85 deg² split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from the ISO data.
This paper presents the concepts underlying the data-reduction method together with its latest enhancements. The data-reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.     

The European Large Area ISO Survey – I. Goals, definition and observations

We describe the European Large Area ISO Survey (ELAIS). ELAIS was the largest single Open Time project conducted by ISO , mapping an area of 12 deg² at 15 μm with ISOCAM and at 90 μm with ISOPHOT. Secondary surveys in other ISO bands were undertaken by the ELAIS team within the fields of the primary survey, with 6 deg² being covered at 6.7 μm and 1 deg² at 175 μm.
This paper discusses the goals of the project and the techniques employed in its construction, as well as presenting details of the observations carried out, the data from which are now in the public domain. We outline the ELAIS 'preliminary analysis' which led to the detection of over 1000 sources from the 15 and 90-μm surveys (the majority selected at 15 μm with a flux limit of ∼3 mJy), to be fed into a ground-based follow-up campaign, as well as a programme of photometric observations of detected sources using both ISOCAM and ISOPHOT.
We detail how the ELAIS survey complements other ISO surveys in terms of depth and areal coverage, and show that the extensive multi-wavelength coverage of the ELAIS fields resulting from our concerted and on-going follow-up programme has made these regions amongst the best studied areas of their size in the entire sky, and, therefore, natural targets for future surveys. This paper accompanies the release of extremely reliable subsets of the 'preliminary analysis' products. Subsequent papers in this series will give further details of our data reduction techniques, reliability and completeness estimates and present the 15‐ and 90-μm number counts from the 'preliminary analysis', while a further series of papers will discuss in detail the results from the ELAIS 'final analysis', as well as from the follow-up programme.     

The European Large Area ISO Survey: optical identifications of 15-μm and 1.4-GHz sources in N1 and N2

We present the multiwavelength properties and catalogue of the 15 μm and 1.4 GHz radio sources detected in the European Large Area ISO Survey ( ELAIS ) areas N1 and N2. Using the optical data from the Wide Field Survey we use a likelihood ratio method to search for the counterparts of the 1056 and 691 sources detected at 15 μm and 1.4 GHz, respectively, down to flux limits of   S ₁₅= 0.5 mJy  and   S _1.4 GHz= 0.135 mJy  . We find that ∼92 per cent of the 15 μm ELAIS sources have an optical counterpart down to   r '= 24  . All mid-infrared (IR) sources with fluxes   S ₁₅≥ 3 mJy  have an optical counterpart. The magnitude distribution of the sources shows a well-defined peak at relatively bright magnitudes   r '∼ 18  . The mid-IR-to-optical and radio-to-optical flux diagrams are presented and discussed in terms of actual galaxy models. About 15 per cent of the sources are bright galactic stars; of the extragalactic objects ∼65 per cent are compatible with being normal or starburst galaxies and ∼25 per cent active galactic nuclei (AGNs). Objects with mid-IR-to-optical fluxes larger than 100 are found, comprising ∼20 per cent of the sample. We suggest that that these sources are highly obscured luminous and ultraluminous starburst galaxies and AGNs.     

The European Large Area ISO Survey – IV. The preliminary 90-μm luminosity function

We present the luminosity function of 90-μm-selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z =0.3. Their luminosities are in the range 10⁹65⁻² L /L_⊙<10¹², i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. S ₉₀100 mJy data base, we have found optical, 15-μm or 1.4-GHz identifications for 24 (65 per cent). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89 per cent of identifications to rigid multivariate flux limits. We construct a luminosity function assuming that (i) our spectroscopic subset is an unbiased sparse sample, and (ii) there are no galaxies that would not be represented in our spectroscopic sample at any redshift. We argue that we can be confident of both assumptions. We find that the luminosity function is well described by the local 100-μm luminosity function of Rowan-Robinson, Helou & Walker. Assuming this local normalization, we derive luminosity evolution of (1+ z )^2.45±0.85 (95 per cent confidence). We argue that star formation dominates the bolometric luminosities of these galaxies, and we derive comoving star formation rates in broad agreement with the Flores et al. and Rowan-Robinson et al. mid-infrared-based estimates.     

Radio source stacking and the infrared/radio correlation at μJy flux densities

We investigate the infrared/radio correlation using the technique of source stacking, in order to probe the average properties of radio sources that are too faint to be detected individually. We compare the two methods used in the literature to stack sources and demonstrate that the creation of stacked images leads to a loss of information. We stack infrared sources in the Spitzer Extragalactic First Look Survey (xFLS) field, and the three northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, using radio surveys created at 610 MHz and 1.4 GHz, and find a variation in the absolute strength of the correlation between the xFLS and SWIRE regions, but no evidence for significant evolution in the correlation over the 24-μm flux density range 150 μJy to 2 mJy. We carry out the first radio source stacking experiment using 70-μm-selected galaxies, and find no evidence for significant evolution over the 70-μm flux density range 10–100 mJy.     

A far-infrared view of the Lockman hole from ISO 95-μm observations – II. Optical identifications and insights into the nature of the far-infrared sources

We present the optical identifications of a 95-μm ISOPHOT sample in the Lockman hole over an area of approximately half a deg². The Rodighiero et al. catalogue includes 36 sources, making up a complete flux-limited sample for   S _95 μm≥ 100 mJy  . Reliable sources were detected, with decreasing but well-controlled completeness, down to   S _95 μm≃ 20 mJy  . We have combined mid-infrared (IR) and radio catalogues in this area to identify the potential optical counterparts of the far-IR sources. We found 14 radio and 13 15-μm associations, 10 of which have both associations. For the 11 sources with spectroscopic redshift, we have performed a spectrophotometric analysis of the observed spectral energy distributions (SEDs). Four of these 95-μm sources have been classified as faint IR (FIR) galaxies  ( L _FIR < 1. e 11 L_⊙)  , six as luminous IR galaxies (LIRGs) and only one as an ultraluminous IR galaxy (ULIRG). We have discussed the redshift distribution of these objects, comparing our results with evolutionary model predictions 95 and 175 μm. Given their moderate distances (the bulk of the closest spectroscopically identified objects lying at   z < 0.2  ), their luminosities and star formation rates (SFR; median value  ∼ 10 M_⊙ yr⁻¹  ), the sources unveiled by ISOPHOT at 95 μm seem to correspond to the low redshift  ( z < 0.3)  FIRBACK 175-μm population, composed of dusty, star-forming galaxies with moderate SFRs. We computed and compared different SFR estimators, and found that the SF derived from the bolometric IR luminosity is well correlated with that computed from the radio and mid-IR fluxes.     

Observations of the Hubble Deep Field South with the Infrared Space Observatory– I. Observations, data reduction and mid-infrared source counts

We present results from a deep mid-infrared survey of the Hubble Deep Field South (HDF-S) region performed at 6.7 and 15 μm with the ISOCAM instrument on board the Infrared Space Observatory ( ISO ). The final map in each band was constructed by the co-addition of four independent rasters, registered using bright sources securely detected in all rasters, with the absolute astrometry being defined by a radio source detected at both 6.7 and 15 μm. We sought detections of bright sources in a circular region of radius 2.5 arcmin at the centre of each map, in a manner that simulations indicated would produce highly reliable and complete source catalogues using simple selection criteria. Merging source lists in the two bands yielded a catalogue of 35 distinct sources, which we calibrated photometrically using photospheric models of late-type stars detected in our data. We present extragalactic source count results in both bands, and discuss the constraints that they impose on models of galaxy evolution, given the volume of space sampled by this galaxy population.     

The European Large Area ISO Survey – III. 90-μm extragalactic source counts

We present results and source counts at 90 μm extracted from the preliminary analysis of the European Large Area ISO Survey (ELAIS). The survey covered about 12 deg² of the sky in four main areas and was carried out with the ISOPHOT instrument onboard the Infrared Space Observatory ( ISO ). The survey is at least an order of magnitude deeper than the IRAS 100-μm survey and is expected to provide constraints on the formation and evolution of galaxies. The majority of the detected sources are associated with galaxies on optical images. In some cases the optical associations are interacting pairs or small groups of galaxies, suggesting that the sample may include a significant fraction of luminous infrared galaxies. The source counts extracted from a reliable subset of the detected sources are in agreement with strongly evolving models of the starburst galaxy population.     

SCUBA observations of the sources detected in the MAMBO 1200-μm survey

We have observed 23 sources from the Max-Planck Millimetre Bolometer (MAMBO) array 1200-μm survey with SCUBA at 850 μm, detecting 19 of the sources. The sources generally have low values for the ratio of 850- to 1200-μm flux. Two possible explanations for the low values are either that the sources are at very high redshifts or that the global properties of the dust in the MAMBO sources are different from the global properties of dust in low-redshift galaxies. If the former explanation is correct, we estimate that 15 of the MAMBO sources lie at   z > 3  .     

Starburst galaxies and structure in the submillimetre background towards the Hubble Deep Field

We use an 850-μm SCUBA map of the Hubble Deep Field (HDF) to study the dust properties of optically-selected starburst galaxies at high redshift. The optical/infrared (IR) data in the HDF allow a photometric redshift to be estimated for each galaxy, together with an estimate of the visible star-formation rate. The 850-μm flux density of each source provides the complementary information: the amount of hidden, dust-enshrouded star formation activity. Although the 850-μm map does not allow detection of the majority of individual sources, we show that the galaxies with the highest UV star-formation rates are detected statistically, with a flux density of about S ₈₅₀=0.2 mJy for an apparent UV star-formation rate of 1  h ⁻² M_⊙ yr⁻¹. This level of submillimetre output indicates that the total star-forming activity is on average a factor of approximately 6 times larger than the rate inferred from the UV output of these galaxies. The general population of optical starbursts is then predicted to contribute at least 25 per cent of the 850-μm background. We carry out a power-spectrum analysis of the map, which yields some evidence for angular clustering of the background source population, but at a level lower than that seen in Lyman-break galaxies. Together with other lines of argument, particularly from the NICMOS HDF data, this suggests that the 850-μm background originates over an extremely wide range of redshifts – perhaps 1≲ z ≲6.     

Hot dust in two hard Chandra X-ray sources

The two brightest hard X-ray sources discovered serendipitously by Chandra in the field of the lensing cluster A2390 are found to have ISOCAM counterparts at 6.7 and 15 μm. We use this fact, together with their non-detection by SCUBA at 850 μm, as the basis for dusty radiative transfer modelling of their infrared spectral energy distributions. For the best-fitting models, we find that the dust that reprocesses the optical–ultraviolet light in these Compton-thin active galactic nuclei (AGN) is heated to near its sublimation temperature (above 1000 K), with an inner radius within a parsec of the nucleus. Some warm-dust models with inner temperatures of 200 K are also formally acceptable. These findings strongly support the obscured AGN hypothesis for the new hard X-ray Chandra sources, which lack both strong emission lines and 850-μm SCUBA detections.     

The SCUBA HAlf Degree Extragalactic Survey – VI. 350-μm mapping of submillimetre galaxies

A follow-up survey using the Submillimetre High-Angular Resolution Camera (SHARC-II) at 350 μm has been carried out to map the regions around several 850-μm-selected sources from the Submillimetre HAlf Degree Extragalactic Survey (SHADES). These observations probe the infrared (IR) luminosities and hence star formation rates in the largest existing, most robust sample of submillimetre galaxies (SMGs). We measure 350-μm flux densities for 24 850-μm sources, seven of which are detected at ≥2.5σ within a 10 arcsec search radius of the 850-μm positions. When results from the literature are included the total number of 350-μm flux density constraints of SHADES SMGs is 31, with 15 detections. We fit a modified blackbody to the far-IR (FIR) photometry of each SMG, and confirm that typical SMGs are dust-rich  ( M _dust≃ 9 × 10⁸ M_⊙)  , luminous  ( L _FIR≃ 2 × 10¹² L_⊙)  star-forming galaxies with intrinsic dust temperatures of ≃35 K and star formation rates of  ≃400 M_⊙ yr⁻¹  . We have measured the temperature distribution of SMGs and find that the underlying distribution is slightly broader than implied by the error bars, and that most SMGs are at 28 K with a few hotter. We also place new constraints on the 350-μm source counts, N ₃₅₀(>25 mJy) ∼ 200–500 deg⁻².     

Variations in 24-μm morphologies among galaxies in the Spitzer Infrared Nearby Galaxies Survey: new insights into the Hubble sequence

To study the distribution of star formation and dust emission within nearby galaxies, we measured five morphological parameters in the 3.6- and 24-μm wavebands for 65 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS) and eight galaxies that were serendipitously observed by SINGS. The morphological parameters demonstrate strong variations along the Hubble sequence, including statistically significant differences between S0/a-Sab and Sc-Sd galaxies. Early-type galaxies are generally found to be compact, centralized, symmetric sources in the 24-μm band, while late-type galaxies are generally found to be extended, asymmetric sources. These results suggest that the processes that increase the real or apparent sizes of galaxies' bulges also lead to more centralized 24-μm dust emission. Several phenomena, such as strong nuclear star formation, Seyfert activity, or outer ring structures, may cause galaxies to deviate from the general morphological trends observed at 24 μm. We also note that the 24-μm morphologies of Sdm-Im galaxies are quite varied, with some objects appearing very compact and symmetric but others appearing diffuse and asymmetric. These variations reflect the wide variation in star formation in irregular galaxies as observed at other wavelengths. The variations in the 24-μm morphological parameters across the Hubble sequence mirror many of the morphological trends seen in other tracers of the ISM and in stellar emission. However, the 24-μm morphological parameters for the galaxies in this sample do not match the morphological parameters measured in the stellar wavebands. This implies that the distribution of dust emission is related to but not equivalent to the distribution of stellar emission.     

The SCUBA Local Universe Galaxy Survey – II. 450-μm data: evidence for cold dust in bright IRAS galaxies

This is the second in a series of papers presenting results from the SCUBA Local Universe Galaxy Survey. In our first paper we provided 850-μm flux densities for 104 galaxies selected from the IRAS Bright Galaxy Sample and we found that the 60-, 100-μm ( IRAS ) and 850-μm (SCUBA) fluxes could be adequately fitted by emission from dust at a single temperature. In this paper we present 450-μm data for the galaxies. With the new data, the spectral energy distributions of the galaxies can no longer be fitted with an isothermal dust model – two temperature components are now required. Using our 450-μm data and fluxes from the literature, we find that the 450/850-μm flux ratio for the galaxies is remarkably constant, and this holds from objects in which the star formation rate is similar to our own Galaxy, to ultraluminous infrared galaxies (ULIRGs) such as Arp 220. The only possible explanation for this is if the dust emissivity index for all of the galaxies is ∼2 and the cold dust component has a similar temperature in all galaxies     . The 60-μm luminosities of the galaxies were found to depend on both the dust mass and the relative amount of energy in the warm component, with a tendency for the temperature effects to dominate at the highest L ₆₀. The dust masses estimated using the new temperatures are higher by a factor of ∼2 than those determined previously using a single temperature. This brings the gas-to-dust ratios of the IRAS galaxies into agreement with those of the Milky Way and other spiral galaxies which have been intensively studied in the submm.     

An ultradeep submillimetre map: beneath the SCUBA confusion limit with lensing and robust source extraction

Extracting sources with low signal-to-noise ratio (S/N) from maps with structured background is a non-trivial task which has become important in studying the faint end of the submillimetre (submm) number counts. In this paper, we study the source extraction from submm jiggle-maps from the Submillimetre Common-User Bolometer Array (SCUBA) using the Mexican hat wavelet (MHW), an isotropic wavelet technique. As a case study, we use a large (11.8-arcmin²) jiggle-map of the galaxy cluster Abell 2218 (A2218), with a 850-μm 1σ rms sensitivity of 0.6–1 mJy. We show via simulations that MHW is a powerful tool for the reliable extraction of low-S/N sources from the SCUBA jiggle-maps and nine sources are detected in the A2218 850-μm image. Three of these sources are identified as images of a single background source with an unlensed flux of 0.8 mJy. Further, two single-imaged sources also have unlensed fluxes <2 mJy, below the blank-field confusion limit. In this ultradeep map, the individual sources detected resolve nearly all of the extragalactic background light at 850 μm, and the deep data allow to put an upper limit of 44 sources arcmin⁻² to 0.2 mJy at 850 μm.     

A SCUBA scan-map of the Hubble Deep Field: measuring the bright end of the submillimetre source counts

Using the 850-μm SCUBA camera on the JCMT and a scanning technique different from other sub-mm surveys, we have obtained a 125 arcmin² map centred on the Hubble Deep Field . The 1 σ sensitivity to point sources is roughly 3 mJy and thus our map probes the brighter end of the sub-mm source counts. We find six sources with a flux greater than about 12 mJy (>4 σ ) and, after a careful accounting of incompleteness and flux bias, estimate the integrated density of bright sources     (68 per cent confidence bounds).     

Timeline analysis and wavelet multiscale analysis of the AKARI All-Sky Survey at 90 μm

We present a careful analysis of the point-source detection limit of the AKARI All-Sky Survey in the WIDE-S 90-μm band near the North Ecliptic Pole (NEP). Timeline analysis is used to detect IRAS ( Infrared Astronomy Satellite ) sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90-μm flux of a source. Combined with a robust noise measurement, the point-source flux detection limit at signal-to-noise ratio  (S/N) > 5  for a single detector row is  1.1 ± 0.1 Jy  which corresponds to a point-source detection limit of the survey of ∼0.4 Jy.
Wavelet transform offers a multiscale representation of the Time Series Data ( tsd ). We calculate the continuous wavelet transform of the tsd and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above 4σ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non- IRAS sources detected in timeline analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non- IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.     

Confusion limit resulting from galaxies: using the Infrared Array Camera on board SIRTF

Recent ISO data have allowed, for the first time, observationally based estimates for source confusion in mid-infrared surveys. We use the extragalactic source counts from ISOCAM in conjunction with K -band counts to predict the confusion resulting from galaxies in deep mid-infrared observations. We specifically concentrate on the near-future Space Infrared Telescope Facility ( SIRTF ) mission, and calculate expected confusion for the Infrared Array Camera (IRAC) on board SIRTF . A defining scientific goal of the IRAC instrument will be the study of high-redshift galaxies using a deep, confusion-limited wide-field survey at 3–10 μm . A deep survey can reach 3-μJy sources with reasonable confidence in the shorter wavelength IRAC bands. Truly confusion-limited images with the 8 μm will be difficult to obtain because of practical time constraints, unless infrared galaxies exhibit very strong evolution beyond the deepest current observations. We find L * galaxies to be detectable to z =3–3.5 at 8 μm, which is slightly more pessimistic than found in 1999 by Simpson & Eisenhardt.     

Tracing spiral density waves in M81

We use Spitzer IRAC 3.6 and 4.5 μm near-infrared data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), optical B, V and I and Two-Micron All-Sky Survey K _s-band data to produce mass surface density maps of M81. The IRAC 3.6- and 4.5-μm data, whilst dominated by emission from old stellar populations, are corrected for small-scale contamination by young stars and polycyclic aromatic hydrocarbon emission. The I -band data are used to produce a mass surface density map by a   B − V   colour correction, following the method of Bell and de Jong. We fit a bulge and exponential disc to each mass map, and subtract these components to reveal the non-axisymmetric mass surface density. From the residual mass maps, we are able to extract the amplitude and phase of the density wave, using azimuthal profiles. The response of the gas is observed via dust emission in the 8-μm IRAC band, allowing a comparison between the phase of the stellar density wave and gas shock. The relationship between this angular offset and radius suggests that the spiral structure is reasonably long-lived and allows the position of corotation to be determined.     

The radio–mid-infrared correlation and the contribution of 15-μm galaxies to the 1.4-GHz source counts

The radio counterparts to the 15-μm sources in the European Large Area ISO Survey southern fields are identified in 1.4-GHz maps down to ∼80 μJy. The radio–mid-infrared correlation is investigated and derived for the first time at these flux densities for a sample of this size. Our results show that radio and mid-infrared (MIR) luminosities correlate almost as well as radio and far-infrared (FIR), at least up to   z ≃ 0.6  . Using the derived relation and its spread together with the observed 15-μm counts, we have estimated the expected contribution of the 15-μm extragalactic populations to the radio source counts and the role of MIR starburst galaxies in the well-known 1.4-GHz source excess observed at sub-mJy levels. Our analysis demonstrates that IR emitting starburst galaxies do not contribute significantly to the 1.4-GHz counts for strong sources, but start to become a significant fraction of the radio source population at flux densities ≲0.5–0.8 mJy. They are expected to be responsible for more than 60 per cent of the observed radio counts at ≲0.05 mJy. These results are in agreement with the existing results on optical identifications of faint radio sources.