Euclid preparation: XVII. Cosmic Dawn Survey. Spitzer Space Telescope observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer Space Telescope’s Infrared Array Camera (IRAC). We combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11 % of the total Spitzer Space Telescope mission time. The resulting mosaics cover a total of approximately 71.5 deg2 in the 3.6 and 4.5 μm bands, and approximately 21.8 deg2 in the 5.8 and 8 μm bands. They reach at least 24 AB magnitude (measured to 5σ, in a 2′′ .5 aperture) in the 3.6 μm band and up to ∼ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16 < [3.6] < 19 is . 0′′ .15. The photometric calibration is in excellent agreement with previous WISE measurements. We extracted source number counts from the 3.6 μm band mosaics, and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned, these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields, which are set to become some of the premier fields for extragalactic astronomy in the 2020s.


Introduction
The Euclid mission will survey 15 000 deg 2 of the extragalactic sky to investigate the nature of dark energy and dark matter and to study the formation and evolution of galaxies (Laureijs et al. 2011). To this end, Euclid will obtain high resolution and high signal-to-noise imaging of a billion galaxies in a broad optical filter to measure their shapes, and in three near-infrared (NIR) filters to measure their colours. It will also obtain high signal-tonoise NIR spectroscopy of about thirty million of these galaxies to measure abundances and redshifts. Additionally, photometric redshifts will be determined by combining the Euclid data with optical photometry from external surveys.
To reach the required precision on cosmological parameters and satisfy the stringent mission requirements on completeness and spectroscopic purity and to shape noise bias, Euclid must also obtain observations with a 40 times longer exposure per pixel than the main survey over regions covering at least 40 deg 2 . To this end, three 'deep' fields were selected by the Euclid Consortium. They are described in detail in Scaramella et al. (in prep.), which gives references to prior studies; we give just a very brief description here. They are as follows: (1) Euclid Deep Field North (EDF-N), a roughly circular, 10 deg 2 region centred on the well-studied north Ecliptic pole; (2) Euclid Deep Field Fornax (EDF-F), also roughly circular and of 10 deg 2 , centred on the Chandra Deep Field South (CDF-S, Mauduit et al. 2012;Lacy et al. 2021), and including the GOODS-S (Giavalisco et al. 2004) and the Hubble Ultra Deep Field (HUDF, Beckwith et al. 2006); and (3) the Euclid Deep Field South (EDF-S), a pill-shaped area of 20 deg 2 with no previous dedicated observations. In addition to these three deep fields, Euclid will observe several fields for the calibration of photometric redshifts (photo-z). These fields need to be observed to a level five times deeper than the main survey, and they are centred on some of the best studied extragalactic survey fields that already have extensive spectroscopic data: (1) the COSMOS field (Scoville et al. 2007;Sanders et al. 2007), (2) the Extended Groth Strip (EGS, Grogin et al. 2011), (3) the Hubble Deep Field North (HDF, also GOODS-N), and (4) the XMM-Large Scale Structure Survey field, which includes the Subaru XMM Deep Survey field (SXDS, Lonsdale et al. 2003;Mauduit et al. 2012;Lacy et al. 2021), and the VIMOS VLT Deep Survey (VVDS) 1 .
While Euclid will observe these fields primarily for calibration purposes, those observations will provide an unprecedented data set to study galaxies at faint magnitudes and high redshifts. The survey efficiency of Euclid in the NIR bands is orders of magnitude greater than that of ground-based telescopes (e.g. VISTA). The Euclid deep fields alone will be 30 times larger and one magnitude deeper than the latest UltraVISTA data release (Moneti et al. 2019) covering the COSMOS field and will reach a depth of 26 mag in the Y, J, and H filters (5 σ). In addition, Euclid carries a wide-field near-infrared grism spectrograph, the Near Infrared Spectrometer and Photometer (NISP), covering the 0.92 < λ < 1.85 µm region, which will provide multiple spectra at numerous grism orientations for more than one million sources to a line flux limit similar to 3D-HST (Brammer et al. 2012) and over an area 200 times larger than the COS-MOS field (depending on the scheduling of the blue grism observations). The observations of the deep fields will result in the most complete and deepest spectroscopic coverage produced by Euclid. Such a spectroscopic data set will be unique for the ree-mail: moneti@iap.fr 1 These are now known as the "Euclid Auxiliary Fields" in Euclid terminology construction of the galaxy environment at cosmic noon and for measuring the star formation rate from the Hα emission line intensity.
The deep and wide NIR data from Euclid are also ideal for detecting significant numbers of high-redshift (7 < z < 10) galaxies, as the Lyman α line is redshifted out of the optical into the NIR. However, in order to distinguish galaxy candidates from stars (primarily brown dwarfs), faint Balmer-break galaxies, and dusty star-forming galaxies at lower redshifts, which can all have similar NIR magnitudes and colours, deep optical, and mid-infrared (MIR) data are also needed (Bouwens et al. 2019;Bridge et al. 2019;Bowler et al. 2020).
The Cosmic Dawn Survey (Toft et al. in prep) aims to obtain uniform, multi-wavelength imaging of the Euclid deep and calibration fields to limits matching the Euclid data for the characterisation of high-redshift galaxies. The optical data will be provided by the Hawaii-Two-0 Subaru telescope/Hyper-SuprimeCam (HSC) survey (McPartland et al. in prep.) for the EDF-N and EDF-F, and likely by the Vera C. Rubin Observatory for EDF-S and EDF-F. For the COSMOS and SXDS fields, optical data are provided by the Subaru HSC Strategic programme (HSC-SSP Aihara et al. 2011).
In this paper, we present the Spitzer Space Telescope (Werner et al. 2004) component of the Cosmic Dawn Survey, which consists of (1) observing the full extent of the deep fields at 3.6 and 4.5 µm (only small parts of the EDF-N and EDF-F had already been observed with IRAC, see above), and (2) processing the new observations together with all relevant archival IRAC data, thus including data at 5.8 and 8.0 µm obtained during the cryogenic mission, and tying them to the latest Gaia astrometric reference system. In this way, we strive to produce the deepest possible and most modern MIR images (mosaics) of the deep and calibration fields.
In addition to being essential for the identification of highredshift galaxies, MIR data are crucial to revealing the stellar mass content of the high-redshift Universe (which is outside the scope of Euclid core science). The Euclid data alone are not sufficient to characterise the stellar masses at z > 3.5, as the Balmer break is redshifted out of the reddest band of the NISP. Without MIR data, the interpretation of spectral energy distributions (SEDs) would rely on rest-frame ultraviolet emission which is strongly affected by dust attenuation and dominated by stellar light of new-born stars. Therefore, integrated quantities like the stellar mass would be highly unreliable (Bell & de Jong 2001). Moreover, photometric redshifts would be prone to catastrophic failures resulting from the misidentification of the Lyman and Balmer breaks (e.g. Le Fèvre et al. 2015;Kauffmann et al. 2020). In summary, Spitzer Space Telescope/IRAC data are crucial for identifying the most distant objects (e.g. Bridge et al. 2019), for improving the accuracy of their photometric redshifts, and deriving their physical properties such as stellar masses, dust content, age, and star-formation rate from population synthesis models (e.g. Pérez-González et al. 2008;Caputi et al. 2015;Davidzon et al. 2017). The build-up of stellar mass, especially when confronted with the amount of matter residing in dark matter halos at high redshifts can be a highly discriminating test for galaxy formation models (Legrand et al. 2019). The extrapolation of recent work in the COSMOS field (Bowler et al. 2020) suggests that hundreds of the rarest, brightest z > 7 galaxies are expected to be discovered in the Euclid deep fields. These provide unique constraints on cosmic reionisation, as the brightest galaxies form in the highest density regions of the Universe, which are expected to be the sites of the first generation of stars and galaxies, and thus of reionisation bubbles (Trac et al. 2008).

Observations
All observations described here were made with IRAC. In brief, IRAC is a four-channel array camera on the Spitzer Space Telescope, simultaneously observing four fields slightly separated on the sky at 3.6, 4.5, 5.8, and 8.0 µm, known as channels 1-4, respectively. Spitzer Space Telescope science observations began in August 2003, but observations in channels 3 and 4 ceased once the on-board cryogen was exhausted (May 15, 2009). During the following 'warm mission' phase, channels 1 and 2 continued to operate until the end of operations in late January 2020, albeit with somewhat lower, but still comparable, performance. The earliest observations presented here are archival observations that were obtained in September 2003; the observations of the dedicated Capak programme began in 2017, and the ones of the dedicated Scarlata programme began in 2019. The dedicated observations continued until January 2020, shortly before the shutdown of the satellite. Figure 1 shows a histogram of the integration time accumulated in bins of 30 days over the observing period. These observations account for almost 1.5 million frames, a total integration time of 34 000 hr (all channels combined), and a total on-target time, omitting overheads, of just over 15 600 hr, or nearly 1.8 yr, which is approximately 11 % of the Spitzer Space Telescope mission time.
For our dedicated observations of EDF-N, EDF-S, and EDF-F we adopted a consistent observing strategy that comprises blocks of 3 × 3 maps with a step size of 310 , a large three-point dither pattern, and four repeats per position. Each block covers a 15 .1 × 15 .1 region with a coverage of 3 × 4 × 100 s exposures per pixel. The block centres are offset between passes in order to ensure uniform coverage and enable self-calibration. Each block forms an astronomical observation request (AOR), in IRAC jargon. All other data included in our processing are archival data. They were obtained with a variety of observing strategies that we did not investigate in detail and that we do not attempt to summarise here. In Appendix C, we list the programme IDs of all the observations processed; the ones of our dedicated observations in bold. The combination of the archival data and our own dedicated data produces a spatially variable depth in most fields; this is discussed further in Sect. 4.
All observations are summarised in Table 1, which gives, for each field and channel, the number of frames (data collection events -DCEs -in IRAC terminology) used to produce the mosaics (we note that this can be lower than the number of frames downloaded as some were discarded; see Sect. 3) together with the total observing time. For channels 1 and 2, on the left side of the table, the information is sub-divided into the cryogenic part and the warm part of the mission.

Pre-processing and calibration
Processing begins with the Level 1 data products generated by the Spitzer Science Center via their 'basic calibrated cata' pipeline (Lowrance et al. 2016), which were downloaded from the NASA/IPAC Infrared Science Archive (IRSA 2 ). They have had all well-understood instrumental signatures removed, have been flux-calibrated in units of MJy sr −1 , and are delivered with an uncertainty image and a mask image; they are described in detail in the IRAC Instrument Handbook 3 . More precisely, we begin from the 'corrected basic calibration data' products, which have file extensions .cbcd for the image, .cbunc for the uncertainty, and .bimsk for the mask. The files are grouped by AORs, namely sets of a few to several hundred DCEs obtained sequentially. All frames are 256 × 256 pixels, the pixels are 1 . 2 wide, and the image file header contains the photometric solution and an initial astrometric solution.
The processing was done region by region. A first pass over the files is used to check the headers for completeness and to discard a few incomplete AORs, which accounts for most of the differences in the number of frames listed in Table 1 between channels 1 and 2 or 3 and 4. This is followed by the correction of the 'first frame' bias effect 4 . Next, the positions and magnitudes of WISE (Wright et al. 2010;Mainzer et al. 2011) andGaia DR2 (Gaia Collaboration et al. 2018) sources falling within the field are downloaded. The Gaia sources are first 'projected' to their location at the time of the observations using the Gaia proper motions. Next, they are identified on each IRAC frame, their observed fluxes and positions are determined in each frame using the APEX software (the point-source extractor in MOPEX 5 ) in forced-photometry mode, and the positions are used to update the astrometric solution of each frame. There are typically 30-40 Gaia DR2 sources available for each frame. In channels 1 and 2, most of them are detected and used for the astrometric correction. In the longer-wavelength channels, 3 and 4, only a few sources in total are detected and usable, but that is still sufficient to determine an astrometric solution with negligible distortion as shown in Sect. 4.2.
An attempt was made to subtract bright stars in order to recover faint sources in their wings. For each AOR, a model star built from the template PSFs described in the IRAC Instrument Handbook 6 (see Fig. 4.9 there) is scaled to the median of the fluxes of the star measured in that AOR, and it is subtracted from each frame (of the AOR). Different templates are available for each filter and separately for the cryogenic and the warm missions. While this procedure worked quite well for moderately bright stars (which are of course the vast majority and which represent only a small loss in area), it introduced significant artefacts around the (few) very bright stars in the final mosaics. These artefacts included diffraction spikes corrected only out to a certain distance (out to where the template extends beyond the frame), other edge effects, and the subtraction of the core of bright galaxies. For these reasons, the bright-star subtraction was not performed and the bright stars are left as they are.

Stacking and image combination
In the next step, we computed a median image for all frames within an AOR, which corrects for persistence in the detectors and also for any residual first-frame pattern that introduces structure in the background. In parallel, a background map is also created by iteratively clipping objects and masking them, and finally that background is subtracted from each frame of the AOR.
The final processing steps consist of resampling the background-subtracted frames onto a common grid with a scale of 0 . 6 pix −1 , that is, half the instrument pixel size, which covers all data in all channels and which is the same in all channels.   Notes. Here, 'Num' is the number of frames used, and 'Time' is the total integration time, in hours, they contribute. The left part of the table is for channels 1 and 2, split between cryogenic and warm mission, the right part is for channels 3 and 4, which were used during the cryogenic mission only. We note that the EDF-S field was observed only during the warm mission.
We experimented with two MOPEX interpolation schemes to produce our final mosaics. We first tried the 'drizzling' (Fruchter & Hook 2002) scheme in which the final value of the output pixels is computed by considering the contribution of each input pixel in a smaller pixel grid in the output image. This procedure has excellent noise properties (it does not suffer from correlated pixels) when many input frames are available, but with few input frames it can produce artefacts in the output images. The second, simpler approach is to compute the value of each output pixel as a linear combination of the input pixel values. Although this procedure produces correlated noise, it works reliably for all the fields considered in this work, which can have widely varying numbers of input images. Noise correlations can be estimated through simulations or by comparing sources in our drizzled and non-drizzled images. These comparisons show that the linear interpolation procedure leads to an underestimation of aperture magnitude errors by 30 − 40%, while the magnitudes themselves are unaffected. Next, we used MOPEX to produce an average-combined image while rejecting outliers and excluding masked regions. Outlier rejection was done using the multi-frame temporal, dual outlier, and box outlier modules from MOPEX. We used the same parameters used by the Spitzer super mosaics, which were tuned to produce balanced outlier rejection for a wide range of data types and uses. The outlier rejection methods were then combined using the MosaicRMASK module. The specific parameters used can A&A proofs: manuscript no. output 3.0 be found in Table 2. Frames with different exposure times are combined using the exposure time weighting feature in MOPEX. This appropriately scales the uncertainties and coverage maps to match the median exposure time before carrying out the combination.
The stacking pipeline also produces the following ancillary characterisation maps: (1) an uncertainty map produced by stacking the input uncertainty maps using the same shifts as for the signal stack, (2) a coverage map giving the number of frames contributing to each pixel, and (3) an exposure time map giving the total exposure time per pixel. As the exposure times are not the same for all the observing programmes, these last two maps are not simply scaled versions of each other.

Spatial variation of the PSF in the stacks
The observations described here were made at many different satellite position angles (PAs), and thus when the images are stacked they must be rotated back to north upwards. This has the effect of rotating the PSF, which is fixed in the satellite's reference frame. Since the PSF is not rotationally symmetric, due in particular to the diffraction spikes, the stacked image of a star will depend on when it was observed. As all parts of the stack were not observed at the same time (or at the same PA), the PSF varies spatially in the stack.
The COSMOS field, which is near the equator, was observable only at specific times and therefore with a very restricted range of PAs; the PSF in the COSMOS stacks is thus quite homogeneous. However, in the EDF-N, which was in a continuous viewing zone, observations were obtained at many different PAs, yielding more complicated and more spatially variable PSF. This effect is very important for PSF-based photometry: the PSF at each position of the stack has to be reconstructed by stacking the nominal PSF at the PAs of the observations at that position, as did, for example, Labbe et al. (2015) for the GOODS-South and HUDF fields, and also Weaver et al. (Weaver et al. (2021), henceforth COSMOS2020) for the production of the COSMOS2020 catalogue. The latter used the PRFmap 7 code by Andreas Faisst. While such photometry is beyond the scope of this paper, we nevertheless provide, for each stack, a table of the PAs of each frame used in the stack. For completeness, those tables also contain the frame coordinates, the MJD of the observation, and the exposure time (see Appendix A for more details).

Products
As an example of the data quality, Fig. 2 shows a zoomed-in image of a section of the EDF-F mosaic in the four channels near the region of maximum coverage. We do not provide figures of the full mosaics here as they would be physically too small to show anything informative other than the overall coverage. Maps of the integration time per pixel for channels 1 and 3 of all the fields are presented in Appendix B. Since channel 2 is observed together with channel 1, and similarly for channels 4 and 3, the paired channels have very similar coverage, albeit slightly shifted in position. The 10 deg 2 circular area of EDF-N and EDF-F and the 20 deg 2 pill-shaped area of EDF-S are easily seen in those figures. Also, and with the exception of EDF-S, for which there are only observations done specifically for this programme and no archival data, the integration time per pixel, and consequently the depth reached, is far from uniform, with only a small part of the total area of each field having been observed for more than a few hours. In fact, the median integration time per pixel is greater than 1 hr for only two fields. Table 3 gives the median and maximum pixel integration time for each field and each channel.    The variation of area covered as a function of exposure time for channels 1 and 3 and for all fields is shown graphically in Fig. 3, which presents a cumulative histogram of the area covered versus exposure time. The intersection of the curve with the vertical axis thus gives the total area covered for that field and these areas are also listed in Table 4. EDF-S is the most uniformly observed field and it covers the largest area, but it is also the shallowest, with only 0.1 hr per pixel on average, and it is also the only field with no channel 3 and 4 data. EDF-F and EDF-N reach the target coverage of 10 deg 2 with about 1 hr of exposure time, with the latter showing deeper coverage over smaller zones. The other fields were covered by many observing programmes with different objectives and which covered specific areas to different depths. The combination of these programmes with our own yields a curve with many plateaus. Finally, there are a few small parts of the EDF-F and HDFN fields that have more than 100 hr of exposure time.

Final sensitivities
We estimated the sensitivities of the stacked images by measuring the flux in circular 2 . 5 diameter apertures randomly placed across each image after masking the regions with detected objects using the SExtractor (Bertin & Arnouts 1996) segmentation map. The sensitivity is then computed as the standard deviation of these fluxes (3σ clipped). This procedure is done in 200 × 200 pixel cells (4 arcmin 2 ). Figure 4 shows the cumulative area covered as a function of sensitivity for the channel 1 mosaics. We note the similarity between this figure and the top panel of Fig. 3 once the latter is rotated by 90 degrees. The solid line shows our total depth, summed over all our survey fields. Also shown in the figure are the published sensitivities of the surveys that are included in our data and analyses. Generally, our measured sensitivities are consistent with literature measurements for surveys of equivalent exposure time.

Validation and quality control
As part of our validation process, we compare photometry and astrometry of sources in our stacks with reference catalogues and also extract number counts that can be compared to previous works.

Catalogue extraction
We began by extracting source catalogues from channel 1 and 2 stacks of all fields using SExtractor. We adopted the usual approach of searching for objects that contain a minimum number of connected pixels above a specified noise threshold (in this case 2 σ) and measuring their aperture magnitudes. In the case of our moderately deep IRAC data, where many sources are blended due to the large IRAC PSF, this approach is known to miss faint sources. However, these faint sources are not required for our quality assessment purposes and a shallower catalogue is entirely sufficient. SExtractor estimates a global background on a grid with mesh size of 32 × 32 pixels (we remind the reader that pixels are 0 . 6 wide). This background is smoothed with a 5 × 5 pixel Gaussian kernel with a FWHM = 1 . 5. For each source, the flux is measured within a circular aperture of 7 diameter, and a local background is estimated within an annulus of width 32 pixels around the isophotal limits. The measured fluxes were converted from MJy/sr to AB magnitude using a zero-point of 21.58 (which accounts for a zero-magnitude flux of 3631 Jy and a pixel size of 0 . 6 8 ), and the latter were converted to total magnitude using the aperture corrections given in the IRAC Instrument Handbook for the warm mission (−0.1164 and −0.1158 for channel 1 and channel 2, respectively), which covers the vast majority of the data, while the correction for the cryogenic mission differs at only the 1-2% level 9 . A list of relevant SExtractor parameters used for the catalogue extraction can be found in Table 5.

Astrometric and photometric validation
Using the catalogues extracted above, we evaluated the astrometric accuracy of our stacked images. For each field, we crossmatched sources with magnitude 16 < [3.6] < 19 within 1 of their counterparts in the Gaia DR2 catalogue. This magnitude range was adopted to ensure only bright, non-blended sources were chosen. We now present a detailed analysis for EDF-N, but other fields are similar. Figure 5 shows the difference between reference and measured coordinates (for clarity, only one point in ten is shown). The heavy blue dashed line gives the size of one pixel in the stacked image (which is half the size of the instrument pixel). Similarly (again showing only one in ten points), Fig. 6 shows, for each coordinate, the difference between the reference and the measured value as a function of position along the other coordinate. The thick red dashed line shows a running median computed over a bin containing 20 points. The flatness of this line indicates that there is no significant spatial variation in astrometric precision. Considering all fields, we find that the 1 σ precision (measured as the RMS of the difference between positions in our catalogue and those in Gaia DR2) is 0 . 15. Furthermore, the median value is always < ∼ 0 . 1, with the exception of the sparselycovered HDF-N field where it is < ∼ 0 . 2. These measurements demonstrate that the astrometric solutions have been correctly applied to the individual images and that the combined images are free of residuals on a scale much smaller than an individual mosaic pixel, which is more than sufficient to measure precise infrared and optical-infrared colours. Finally, we perform a simple check on the photometric calibration of our mosaics. As described previously, individual images are photometrically calibrated by the Spitzer Science Center (SSC). Following the validation procedures outlined by the SSC, we compare magnitudes of objects in our catalogues with

Magnitude number counts
We computed the differential number counts in channel 1 in each field using the corrected 7 aperture magnitudes. Since the IRAC PSF is too large to perform morphological source classification, we simply included all objects detected. These are shown in Fig. 8, where the red circles with uncertainties present our measurements, and the lines show the number counts from the literature; the bottom-right panel shows the mean of all fields. We compared our number counts with those presented in Ashby et al. (2013), which also surveyed many of our fields and with those computed using the new COSMOS2020 photometric catalogue (COSMOS2020) that we used as the reference.
There is a general agreement in the number counts in all the fields with Ashby et al. (2013) and COSMOS2020 for 16 < [3.6] < 22. At brighter magnitudes, the COSMOS2020 counts drop off as bright sources were not included. At fainter magnitudes, our aperture-based catalogues are confusion-limited and thus incomplete. Conversely, the COSMOS2020 catalogue, which uses a high-resolution prior for the detection and a profilefitting method for the measurement, is complete up to significantly fainter magnitudes.
Counts for EDF-N are slightly higher than the other fields at bright magnitudes. To investigate this difference, we simulated a stellar catalogue of 1 deg 2 centred on EDF-N using TRILEGAL (Girardi et al. 2005) and compared counts from this simulated catalogue with our observations, shown in Fig. 9. At bright magnitudes, where stars are expected to outnumber galaxies, our counts are in reasonable agreement with TRILEGAL predictions, and in excellent agreement with the number counts extracted from the AllWISE (Wright et al. 2010) catalogue for this field. These comparisons indicate that the difference between EDF-N and other fields is largely due to the higher density of stellar sources in there, which is consistent with its lower Galactic latitude.

Summary
We have presented the Spitzer Space Telescope/IRAC midinfrared component of the Cosmic Dawn Survey: an effort to complement the Euclid mission's observations of deep and calibration fields with deep longer-wavelength data to enable high redshift legacy science.
The survey consists of two major new programmes covering the three Euclid deep fields (EDF-N, EDF-F and EDF-S) and a homogeneous reprocessing of all existing data in Euclid's four calibration fields (COSMOS, XMM, EGS and HDFN). We have processed new data together with all relevant archival data to produce mosaics of these fields covering a total of ∼ 71 deg 2 in IRAC channels 1 and 2. Furthermore, the new mosaics are tied to the Gaia astrometric reference system. The MIR data will be essential for a wide range of legacy science with Euclid, including improved star/galaxy separation, more accurate photometric redshifts, determination of stellar masses of galaxies, and the construction of complete galaxy samples at z > 2 with well understood selection effects.
We validated our final products by comparing catalogues extracted from channels 1 and 2 to external catalogues. In all fields, comparing with Gaia DR2, the residual astrometric uncertainty for sources with total magnitudes 16 < [3.6] < 19 is around 0 . 15 (1σ). Our photometric measurements are in excellent agreement with WISE photometry and our number counts are consistent with previous determinations. The Cosmic Dawn Survey Spitzer Space Telescope survey presented here represents the first essential step in assembling the required multi-wavelength coverage in the Euclid deep fields which are set to become some of the most important fields in extragalactic astronomy for the coming decade. Since the Spitzer Space Telescope mission has finished, and all available data in these fields have been processed with the latest reduction pipeline, the resulting mosaics will remain the deepest and widest MIR imaging survey for the foreseeable future. No existing or approved future observatories are capable of obtaining such data. While JWST is more sensitive and has higher spatial resolution at these wavelengths, its mapping speed is too slow to cover comparable degree-scale areas.
In the context of the Cosmic Dawn Survey, several programmes are currently underway to add data at other wavelengths to the Euclid deep fields and calibration fields. In particular deep optical data in the EDF-N and EDF-F are currently being obtained with the Subaru's Hyper-Suprime-Cam instrument as part of the Hawaii-Two-0 observing programme (McPartland et al. in prep). These fields are also being targeted with high spatial resolution millimeter observations as part of the planned Large-scale Structure Survey with the Toltech Camera 10 on the Large Millimeter Telescope (LMT Pope et al. 2019). A deep Uband survey is also underway with the CFHT (Zalesky et al. in prep). EDF-S is being covered with K-band observations from the VISTA telescope (Nonino et al. private communication), and planning is ongoing to obtain optical data with the Vera C. Rubin Observatory.
The Cosmic Dawn Survey Spitzer Space Telescope mosaics and associated products described here can be downloaded from the IRSA web site, Appendix A gives the details of the download site and the naming convention used. The community is encouraged to make use of them for their science. Integration time maps for channel 3. A blank field is included for EDF-S, which was not observed in that channel, and in order to have the same structure as B.1.