Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer 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's Infrared Array Camera (IRAC). We have 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 mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um 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 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um 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, spectroscopic purity and shape noise bias, Euclid must also obtain observations with 40 times longer exposure per pixel than the main survey over regions covering at least 40 deg 2 .To this end, three 'deep' fields have been selected by the Euclid Consortium.They are described in detail in Scaramella et al. (2021) and we give just a very brief description here.They are: (1) the Euclid Deep Field North (EDF-N), a roughly circular, 10 deg 2 region centred on the well-studied north Ecliptic pole, (2) the Euclid Deep Field Fornax (EDF-F), also roughly circular and of 10 deg 2 , centred on the Chandra Deep Field South and including the GOODS-S (Giavalisco et al. 2004) and the Hubble Ultra Deep Field (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 5 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 (2) the Extended Groth Strip (EGS), (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), and 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 to faint magnitudes and high redshifts.The survey efficiency of Euclid in the NIR bands is orders of magnitudes 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 covering the COSMOS field, and will reach a depth of 26 mag in the Y, J, 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 COSMOS 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 reconstruction 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 all 1 These are now known as the "Euclid Auxiliary Fields" in Euclid terminology can 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 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 program (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, consisting primarily of 3.6 and 4.5 µm observations of the three deep fields and parts of the calibration fields acquired with Spitzer's IRAC camera (Fazio et al. 2004).Two dedicated programs were submitted for this purpose: the Euclid/WFIRST Spitzer Legacy Survey (SLS, requesting 5 286 h, PI: Capak) covering the EDF-N and EDF-F fields, and the EDF-S survey (requesting 687 h, PI: Scarlata).These programmes were established based on the Euclid plans for the deep fields that were available at the time.All the fields had been observed, at least in part, before our new observations, and we processed our new data together with all relevant archival IRAC data, thus including data obtained during the cryogenic mission, i.e., data at 5.8 and 8.0 µm.In this way we strive to produce the deepest possible MIR images (mosaics) of these fields to date.A significant improvement in our processing is that our pipeline ties the astrometry to the Gaia reference system which, given its higher precision, will greatly facilitate cross-identification with other data, which will of course also have to be tied to Gaia.
In addition to being essential for the identification of highredshift galaxies, MIR data are crucial to reveal the stellar mass content of the high-redshift Universe (which is outside the scope of the 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 mis-identification of the Lyman and Balmer breaks (e.g.Le Fèvre et al. 2015;Kauffmann et al. 2020).In summary, Spitzer/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 for 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).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).
Article number, page 2 of 15 The layout of this paper is as follows: Sect. 2 describes the observations, Sect. 3 presents our data processing techniques, and Sect. 4 compares our results to previous ones.

Observations
All observations described here were made with IRAC.In brief, IRAC is a four-channel array camera on the Spitzer Space Telescope, observing simultaneously 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 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 program began in 2017 and the ones of the dedicated Scarlata program 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 ontarget 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 AOR, or Astronomical Observation Request, in IRAC jargon.All other data included in our processing is archival data.It was obtained with a variety of observing strategies which we did not investigate in detail and which we do not attempt to summarise here.In Appendix C we list the Program IDs of all the observations processed; in bold the ones of our dedicated observations.The combination of the archival data with our own dedicated data produces a spatially variable depth in most fields; this is discussed further in Sect. 4.
A total of 292 IRAC observing programs are used in this work.Table 1 lists the ten largest programs in terms in terms of observing time together with the PI of the program, the field concerned and program's total integration time.
All observations are summarised in Table 2 which gives, for each field and channel, the number of frames (Data Collection Events or DCEs in IRAC terminology) used to produce the mosaics (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 subdivided 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 Data' pipeline (Lowrance et al. 2016), which were downloaded from the NASA/IPAC Infrared Science Archive (IRSA2 ).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 Handbook3 .More precisely, we begin from the 'corrected basic calibration data' products, which have file extensions .cbcdfor the image, .cbuncfor the uncertainty, and .bimskfor 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 is 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 2 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)) and Gaia 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 determined in each frame using the APEX software (the point-source extractor in MOPEX5 ) 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 Handbook6 (see Fig. 4.9 there) is scaled to the median of the fluxes of the star measured in that AOR, and 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 compute 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 final processing steps consist of resampling the background-subtracted frames onto a common grid with a scale of 0 .6 pix −1 , i.e., half the instrument pixel size, that covers all data in all channels and which is the same in all channels.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 use MOPEX to produce an average-combined image while rejecting outliers and excluding masked regions.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.

On the 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 Notes.Here 'Num' is the number of frames used, and 'Time' is the 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 & 4 which were used during the cryogenic mission only.Note that the EDF-S field was observed only during the warm mission.
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.But 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 e.g., Labbe et al. (2015) for the GOODS-South and HUDF fields, and also Weaver et al. (submitted) for the production of the COSMOS2020 catalogue.The latter used the PRFmap code by Andreas Faisst, available at https://github.com/cosmic-dawn/prfmapfor that purpose.While doing 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 section of the EDF-F mosaic in the four channels near the region of maximum coverage.We do not provide here figures of the full mosaics 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  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 larger than 1 hr for only two fields.That 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 vs. 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 programs with different objectives and which covered specific areas to different depths.The combination of these programs 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 estimate 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.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 begin by extracting source catalogues from the channels 1 and 2 stacks of all fields using SExtractor.We adopt 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 (recall that pixels are 0 .6 wide).This background is smoothed with a 5 × 5 pixel Gaussian kernel with 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 .67 ), 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 8 .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 evaluate the astrometric accuracy of our stacked images.For each field we crossmatch 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 Article number, page 7 of 15

Magnitude number counts
We compute 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 include 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 compare our number counts with those presented in Ashby et al. (2013) who also surveyed many of our fields and also with those computed using the new COSMOS2020 photometric catalogue (Weaver et al., submitted) which we use as a 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.
EDF-N counts 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, consistent with its lower Galactic latitude.

Summary
We have presented the Spitzer/IRAC mid-infrared 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 programs 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.become some of the most important fields in extragalactic astronomy for the coming decade.Since the Spitzer mission has finished, and all available data in these fields have been pro-cessed 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 programs 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 program (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 9 on the Large Millimeter Telescope (LMT Pope et al. 2019).A deep U-band 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 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.
Acknowledgements.We thank the MOPEX support team for fixing issues that appeared when combining large numbers of files.This publication is based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA, and has made use of the NASA/IPAC Infrared Science Archive, which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology.This publication has also made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia),processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium).Funding for the Gaia Data Processing and Analysis Consortium (DPAC) has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.This publication makes use of data products from the Wide-field Infrared Survey Explorer (WISE), which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.H.J.McC.acknowledges support from the PNCG.This work used the CANDIDE computer system at the IAP supported by grants from the PNCG and the DIM-ACAV and maintained by S. Rouberol.S.T. and J.W. acknowledge support from the European Research Council (ERC) Consolidator Grant funding scheme (project ConTExt, grant No. 648179).I.D. has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 896225.The Cosmic Dawn Center is funded by the Danish National Research Foundation under grant No. 140.H.Hildebrandt is supported by a Heisenberg grant of the Deutsche Forschungsgemeinschaft (Hi 1495/5-1) as well as an ERC Consolidator Grant (No. 770935).The Euclid Consortium acknowledges the European Space Agency and a number of agencies and institutes that have supported the development of Euclid, in particular the Academy of Finland, the Agenzia Spaziale Italiana, the Belgian Science Policy, the Canadian Euclid Consortium, the Centre National d'Etudes Spatiales, the Deutsches Zentrum für Luft-und Raumfahrt, the Danish Space Research Institute, the Fundação para a Ciência e a Tecnologia, the Ministerio de Economia y Competitividad, the National Aeronautics and Space Administration, the National Astronomical Observatory of Japan, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Romanian Space Agency, the State Secretariat for Education, Research and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency.A complete and detailed list is available on the Euclid web site (http://www.euclid-ec.org).Notes.Longitude and latitude are Equatorial and J2000, for the image tangent point.These values are valid for all four channels of each field and for their ancillary images.-1700 1910-1949 1951 1953-1961 1963-1983

Fig. 1 .
Fig.1.Histogram of the exposure time of the data analysed here (including the few discarded observations) using bins of 30 days.Our dedicated observations began in November 2016 and comprise most of the data after that date.The red part of each bar accounts for observations in channels 3 and 4, the blue part for those in channels 1 and 2; the vertical dotted line at 2009.37 indicates the end of the cryogenic mission.No observations are made in channel 3 and 4 after the end of the cryogenic mission.

Fig. 2 .
Fig. 2. Detail of EDF-F mosaic in the region near that of maximum exposure time, which here is the same for all four channels.Images are 200 × 200 pixels, or 2 × 2 .Display levels are −σ to +8σ, where σ is the standard deviation of the sky pixels, which is ∼ 0.005 MJy sr −1 for channels 1 & 2, and ∼ 0.013 MJy sr −1 for channels 3 & 4.

Fig. 4 .
Fig.4.Sensitivity of the Spitzer/IRAC channel 1 data as a function of cumulative area coverage.The coloured lines illustrate 1 σ depths measured in empty 2 .5 diameter apertures in each field.The grey solid line is the total area observed to a given depth summed over different surveys.The data points indicate point-source sensitivities at 1σ compiled inAshby et al. (2018) (note that some of these data are included in our stacks).The circles and squares represent surveys executed during cryogenic and warm missions, respectively.

Fig. 5 .
Fig. 5.The difference between the reference and the measured position, in arcseconds, of Gaia DR2 catalogue sources with 16 < [3.6] < 19 total magnitudes extracted from the EDF-N channel 1 mosaic.The blue dashed lines indicate the size of one mosaic pixel.The blue dotted lines go through the origin.The shaded regions are ellipses containing 68 % and 99 % of all sources respectively.For clarity, only one in ten sources is plotted.

Fig. 6 .Fig. 7 .
Fig. 6.The difference between the reference Gaia DR2 catalogue and the measured RA (top panel) and DEC (bottom panel) of sources in the EDF-N channel 1 mosaic with 16 < [3.6] < 19 total magnitudes as a function of the coordinate.The solid red line shows a running median computed in bins of 20 points, and the shaded areas indicate the regions containing 68 % and 99 % of all sources respectively.

Fig. 8 .Fig. 9 .
Fig. 8. Magnitude number counts in channel 1 (red circles) together with COSMOS2020 (long dashed lines) and Ashby et al. (2013, hereafter A13; short dotted green lines).The bottom right panel shows the mean of all fields compared to , and the legend there applies to all panels.

Table 1 .
The ten largest programmes, by Program ID 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.

Table 3 .
Cumulative area coverage as a function of exposure time for channels 1 and 3, for all fields.The figures for channels 2 and 4 are similar to the ones above, as explained in the text.Median and maximum pixel integration time in hours Table 3 gives the median and maximum pixel integration time for each field and each channel.Article number, page 5 of 15 A&A proofs: manuscript no.Cosmic_dawn_spitzer

Table 4 .
Location and area, in deg 2 , covered in each field

Table 5 .
SExtractor parameters used for detection and photometry.
The Cosmic Dawn Survey Spitzer survey presented here represents the first essential step in assembling the required multiwavelength coverage in the Euclid deep fields which are set to Article number, page 8 of 15 Moneti et al.: The Cosmic Dawn Survey

Table A .
1. Data products information