Free Access
Issue
A&A
Volume 662, June 2022
Article Number A112
Number of page(s) 41
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/202141938
Published online 28 June 2022

© ESO 2022

1. Introduction

Observations of distant type Ia supernovae (e.g. Riess et al. 1998; Perlmutter et al. 1999) together with those of the cosmic microwave background (CMB; e.g. de Bernardis et al. 2000, 2002; Hanany et al. 2000; Pryke et al. 2002; Hinshaw et al. 2013; Planck Collaboration I 2020; Planck Collaboration VI 2020) suggest that the spatial curvature of the Universe is close to zero. Despite the large contribution of dark matter (DM), however, the total matter density is still much lower than the critical matter density. As a consequence, a non-zero value for the cosmological constant, Λ, is usually introduced to complete the cosmological model. Although a cosmological constant can describe the data, it is not generally appealing (Weinberg 1989), and alternative solutions have been investigated, such as an evolving quantum field – dark energy (DE) – and a modification to general relativity on cosmological scales. We refer to Amendola et al. (2018) for an extensive review of theoretical models.

To learn more about the nature of DE and DM, we need to quantify their impact on cosmological observations. In particular, we need to determine H(z), the expansion history of the Universe as a function of redshift, z, using geometrical tests and to measure the growth of large-scale structures through gravitational instability. The latter can be captured using the time derivative of the matter density contrast δ ≡ δρ/ρ,

d ln ( δ ) d ln ( a ) Ω m γ , $$ \begin{aligned} \frac{\mathrm{d}\ln (\delta )}{\mathrm{d}\ln (a)} \, \simeq \, \Omega _{\rm m}^\gamma , \end{aligned} $$(1)

where a ≡ 1/(1  +  z) is the cosmic scale factor and Ωm is the mean density divided by the critical density. For the canonical Λ cold dark matter (CDM) model in linear theory, γ ≃ 0.55, whereas it differs for other models of DE (Amendola et al. 2018). We note that the growth of cosmic structures is also influenced by the DM characteristics.

The exact nature of DE can be tested via its equation-of-state, w = p/(ρc2), which directly influences the expansion history. In the general case, w is a function of the scale factor, w = w(a), and simplifies to w = −1 in the case of a cosmological constant. Using a truncated Taylor expansion, one can write w  =  w0  +  wa (1 − a) and seek constraints on the possible values in the w0 − wa plane. Ideally, we should have a model that describes the redshift dependence of w(a) and predicts how the structure growth is affected by a modification of gravity. Nonetheless, γ, w0, and wa provide convenient generic parameterisations that can be used to compare the expected performance of various cosmological probes. When Fisher matrix techniques are used (Euclid Collaboration 2020a), the confidence areas in two-dimensional parameter spaces are ellipses. The inverse of the area of the 2σ ellipse in the w0 − wa plane, after marginalisation over all other cosmological and nuisance parameters, defines the DE figure of merit (FoM; Albrecht et al. 2006; Laureijs et al. 2011). Hence, the larger the FoM, the better (more informative) the experiment is.

Two of the best cosmological probes are galaxy clustering (GC) and weak gravitational lensing (WL), especially once combined in the so-called 3 × 2 pt statistics. These are the two-point correlation of galaxies positions (GC uses galaxies as test particles in the expanding space-time to map the mass density contrast, δ, over time), the shear two-point correlations (WL exploits the cumulative distortion effect of the tidal gravitational fields along the line of sight on the shapes of the galaxy images; see e.g. Kilbinger 2015, for a review), and the cross-correlation of the lens positions with the shear of the source galaxies, known as galaxy-galaxy lensing. Despite tremendous progress in GC and WL experiments in recent years (e.g. BOSS Collaboration 2017; Dark Energy Survey Collaboration 2022; Gil-Marin et al. 2020; Hildebrandt et al. 2020), much larger cosmological volumes need to be surveyed.

Even though a clear ‘target precision’ is lacking, Laureijs et al. (2011) argues that a FoM ≥ 400 provides constraints on w(a) and γ that can test key aspects of our current cosmological model. This target FoM ≥ 400 has driven the design of Euclid, a medium class mission of the European Space Agency (ESA) that combines GC and WL.

To meet its primary science goal, Euclid (Laureijs et al. 2011; Racca et al. 2016) has to observe a large fraction of the extra-galactic sky with both multi-band imaging and slitless spectroscopy. The sky area and mean number density of galaxies are specified by the scientific requirements of the GC and WL experiments (Laureijs et al. 2011; Rassat et al. 2008; Cropper et al. 2013; Massey et al. 2013):

(i) a 15 000 deg2 survey of the extra-galactic sky, jointly for WL and GC, to be completed in six years with all the necessary calibrations; (ii) an average galaxy number density of 30 arcmin−2 that is useful for WL in the optical imaging data; and (iii) an average galaxy number density of 1700 deg2 with reliable redshifts from the Hα emission line spectroscopic data, useful for GC1.

Moreover, we want to minimise systematic residuals in the error budget by maximising the uniformity in covering the observed sky by using the same observation sequence (see Sect. 4.1 and Fig. 8) as much as possible (Laureijs et al. 2011; Scaramella et al. 2015). The resulting survey is the Euclid Wide Survey (EWS), which will cover 15 000 deg2 to a minimum depth of mAB = 24.5 mag in the visible band IE with a signal-to-noise ratio (S/N) of 10 for extended sources such as the z ∼ 1 galaxies (details in Sect. 5.2.2 and Laureijs et al. 2011; Cropper et al. 2016). In the near-infrared (NIR) YE, JE and HE bands, a depth of mAB = 24.0 mag will be reached with a minimum S/N of 5 for point sources (Laureijs et al. 2011). This is sufficient to complement ground-based multi-band observations that will be used to determine photometric redshifts (photo zs) for the WL sources (Euclid Collaboration 2020b). Using slitless spectroscopy, Euclid will detect line emission with a sensitivity of fHα ≥ 2  ×  10−16 erg s−1 cm−2 and a S/N of 3.5 for a typical source of size 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $ (Maciaszek et al. 2014, 2016). Space-based observations provide excellent and consistent image quality at visible wavelengths for WL and sufficient depth in the NIR bands, both of which are unattainable from the ground for this type of survey.

With this design, Fisher matrix analyses forecast that EWS 3 × 2 pt datasets will obtain a DE FoM of 500 for a non-flat w0 − wa CDM model in an optimistic setting (defined by the range of scales used; see Euclid Collaboration 2020a). It is interesting to compare Euclid 3 × 2 pt forecasts with the constraints obtained by DES, which can be considered as representative of ongoing Stage-III surveys. In DES Year 1 constraints for a w0 − wa CDM model (Dark Energy Survey Collaboration 2019), the uncertainties on w0 and wa are a factor of 10 larger than the Euclid forecasts in the pessimistic setting (roughly translating to a factor of 100 in the FoM). The DES DE FoM from the DES Year 3 analysis of a more complete dataset is not yet available, but we can use instead the results on S8 = σ8m/0.3)1/2 to compare Stage-III constraints with Stage-IV Euclid forecasts. Using a 3 × 2 pt dataset and marginalising over 25 nuisance parameters and 7 cosmological parameters of a flat wCDM model, DES Year 3 gets 3.2% errors on S8 (Dark Energy Survey Collaboration 2022). Using the Fisher matrix for 3 × 2 pt obtained in Euclid Collaboration (2020a) for a w0 − wa CDM model, and adding an additional ten nuisance parameters to represent the uncertainty on the shear multiplicative bias of each redshift bin (in order to increase the number of nuisance parameters to have a fair comparison with DES), we marginalise over a total of 25 nuisance and cosmological parameters, getting 1.25% errors on S8 in a pessimistic setting and 0.68% errors in an optimistic setting. Although the comparison is not completely fair – because of the different intrinsic alignments and bias models, the fewer nuisance parameters considered in the Euclid forecast, and a slightly different set of cosmological parameters – these numbers are roughly in agreement with what one would obtain by simply scaling the DES Year 3 results for the Euclid increase in area and source number density.

The EWS has to be carried out within the six-year mission baseline and will start 3 months after the launch, following a commissioning period (1 month) and a performance verification (PV) period (2 months). Figure 1 shows the Euclid timeline with the data release planning. The first major data release (DR1), corresponding to 2500 deg2 of the EWS, is planned to take place one year after T1 (T1 = 14 months after launch), the second data release (DR2) is expected to release 7500 deg2 three years after T1, and the final one (DR3) will release the full survey (15 000 deg2) six years after T1. In between, there will be other ‘quick data releases’: Q1 of 50 deg2 is planned at T1, and Q2, Q3, and Q4 will take place two, four, and five years after Q1, respectively.

thumbnail Fig. 1.

Euclid timeline. The EWS has to be carried out within the six year mission baseline and will start 3 months after the launch, following a commissioning period of 1 month and a PV period of 2 months. An extension of operations beyond the six years is possible and will be decided in due time. During the first six months of the main mission (the early operations phase) a faster replanning of the survey is allowed. The three main data releases (DR#) are shown. The plan is to have 2500 deg2 made public in DR1, to grow to 7500 in DR2, and be complete at DR3 for 15 000 deg2. In addition, four quick data releases (Q#) are foreseen, each of ∼50 deg2.

In addition to the main survey, a significant fraction of time will be spent on calibrating the instruments and characterising the target galaxies. As a result, some fields will be observed to greater depths than with the wide survey (typically two magnitudes deeper). These deep fields have a great legacy value beyond the cosmological core science. While aspects of non-core science did not influence the design of the spacecraft and instruments, they are taken into account in the design of the EWS to maximise the Euclid scientific return. In fact, it must be noticed that the large decrease in the background with the wavelength dramatically increases the S/N in the NIR bands when compared to Earth-based observations affected by airglow, which instead increases with wavelength. This makes even a small space telescope competitive with a large ground telescope that suffers from a background dominated by atmospheric emission in the NIR bands. The relative gain is such that, in order to cover the same areas planned for Euclid and at the same depths, a ground-based NIR survey on existing facilities would need to observe for several centuries. Regarding other space-based facilities, we notice that the James Webb Space Telescope (JWST) will be in orbit as well and, with its diameter of 6.5 m, will go much deeper and faster than Euclid, although only on very small areas (the JWST field of view is 75 times smaller than the Euclid one). Hence, the two facilities are complementary and, moreover, JWST will likely benefit from targets selected from the Euclid surveys.

The challenge is to fit all these observations into a finite time allocation set by the limitation of the mission, which is six years, whilst fulfilling a wide range of constraints, which are reviewed in detail in this paper. Part of the survey optimisation involves selecting the best areas of the sky to use, which in turn relies on a good model of the properties of the observable sky, such as Galactic extinction and the zodiacal background. We also need to model the distribution of (bright) stars as their stray light lowers the observed galaxy number density.

This paper focuses on the design of the EWS, and the deep fields will be described in a companion paper (Scaramella et al., in prep., hereafter [Sc23]). The EWS design takes into account the main backgrounds that impact any large area survey, the sequence of operations, and the many limitations to the pointing of the telescope. The EWS is at an advanced stage, fulfilling the key survey requirements over the full mission. Survey scenarios at this stage therefore show the detailed feasibility of the mission but are subject to further optimisation. Nevertheless, the results we present and their discussion are instructive and useful for any future large area survey from space or the ground that aims to combine imaging and spectroscopy.

The paper is organised as follows. The spacecraft is described in Sect. 2, followed by a summary of Euclid’s instruments in Sect. 3. In Sect. 4 the reference observation sequence (ROS) is introduced, including a study of dithering scenarios. Models of zodiacal light, stray light effects, and other environmental properties define the region of interest (RoI) used as input for the implementation of the Euclid reference survey definition (RSD). These effects and the properties of the resulting RoI are presented in Sect. 5, where we also discuss complementary ground-based observations. Section 6 describes the implementation of the calibration programme. Observations of sample characterisation fields and the Euclid Deep Survey (EDS) are briefly mentioned in this context. The construction of the EWS is presented in Sect. 7. We present the most recent outcome of the survey optimisation (mid 2021) in Sect. 8. This solution is a good proxy for the actual survey. We conclude in Sect. 9.

In the appendix we provide a list of the acronyms used in this paper.

2. The spacecraft and telescope

2.1. The spacecraft

The spacecraft comprises a service module (SVM) and a payload module (PLM), connected by an interface structure designed to maximise thermal decoupling. The PLM includes the main instruments, the folded beam optical components of the telescope, the radiators, and the fine guidance system (FGS). The SVM provides the main spacecraft services – power generation, conditioning and distribution, Sun shield and solar array, telecommunication with the ground (low and high gain antenna), and attitude and orbit control system (AOCS; including the FGS) – and support the instruments’ warm electronics. Details are given in Laureijs et al. (2011) and Racca et al. (2016).

Euclid has severe constraints in pointing to ensure maximal thermal stability, which are described in this paper and limit the standard operations. Therefore, it is important to describe in detail the attitude of the spacecraft. The Euclid spacecraft reference frame (OSC, XSC, YSC, ZSC) is defined as follows (see Fig. 2 for a graphical representation): (i) OSC, the origin is at the point of intersection of the longitudinal launcher axis with the launcher adapter interface plane (the plane of separation of the spacecraft from the launcher); (ii) +ZSC in the direction perpendicular to the launcher interface plane, positive in the direction of the launch; (iii) +XSC in the launcher interface plane, directed to a physical mark on the interface ring nominally aligned with the solar array such that the +XSC vector is perpendicular to the solar array and pointing towards the Sun; and (iv) +YSC in the remaining direction of the right-handed orthogonal triad.

thumbnail Fig. 2.

Euclid reference frames. Left panel: Euclid spacecraft reference frame. XSC points towards the Sun disk centre. The edge and centre of the joint FoV are offset by 0 . ° 47 $ 0{{\overset{\circ}{.}}}47 $ and 0 . ° 82 $ 0{{\overset{\circ}{.}}}82 $, respectively. The longer side of the FoV is typically aligned with ecliptic meridians during observations. Right panel: Euclid optical reference frame. We note that the XoptYopt plane is defined as looking onto the sky with z ̂ $ \hat z $ towards the spacecraft. There are four additional chips used as fine guidance sensors that are placed on each side of the VIS FoV.

The orientation of the telescope optical reference frame, projected onto the sky, is also specified in Fig. 2. The field of view (FoV) reference frame is centred on the centre of the FoV itself and is such that XFoV = −XOPT and Y FoV = Y OPT 0 . ° 82 $ Y_{\mathrm{FoV}}=Y_{\mathrm{OPT}}-0{{\overset{\circ}{.}}}82 $, taking into account the shift of the edge of the FoV of 0 . ° 47 $ 0{{\overset{\circ}{.}}}47 $, and its half size of 0 . ° 35 $ 0{{\overset{\circ}{.}}}35 $ (see Fig. 2).

2.2. The telescope

Euclid’s PLM (Racca et al. 2016) is designed around a three-mirror anastigmat Korsch design telescope with silicon carbide (SiC) mirrors and truss (Korsch 1972; Pamplona et al. 2016). The sizes of the telescope components are: primary pupil R1 = 0.6 m, primary mirror (M1) stopper R2 = 0.1975 m, spider arm mean length R3 = 0.44 m, spider arm thickness L = 0.012 m. This provides a total collecting area of A=π( R 1 2 R 2 2 )3 R 3 L=0.99 $ A=\pi \left(R_1^2-R_2^2\right)\,{-}\,3 R_3 L = 0.99 $ m2. Euclid has two instruments on board, the visible imager (VIS; Sect. 3.1) and the near-infrared spectrometer and photometer (NISP; Sect. 3.2). The wavelength separation at ∼920 nm between the two instruments is performed by a dichroic plate located at the exit pupil of the telescope. The two focal planes image the same part of the sky, allowing multiple data acquisition with a single telescope pointing (see Sect. 4.1). The coordinates of the focal plane as projected on the sky are shown in Fig. 2.

2.3. Pointing angles

The main reference frames are shown in Fig. 2. The following angles (see Fig. 3) have operational ranges that constrain a pointing and therefore target visibility: (i) solar aspect angle (SAA), which is the angle between the spacecraft’s +ZSC axis (telescope pointing direction) and the direction to the centre of the Solar disk; (ii) alpha angle (AA), which is the angle between the Sun vector projected onto the XSC − YSC plane and the +XSC axis, increasing as the spacecraft rotates clockwise about its +ZSC axis; and (iii) solar panel solar aspect angle (SPSAA), which is the angle between the spacecraft +XSC axis and the direction to the centre of the Solar disk.

thumbnail Fig. 3.

Illustration of the relevant pointing angles defined in the text. In the top right image the minimal SAA (pointing towards the Sun) is indicated, while the maximum SAA (pointing away from the Sun) adopted in the survey is 110° (see also Sect. 2.5). In the lower right image, AA is shown; the maximum range allowed for the survey is ±5°, which corresponds to a large margin with respect to the spacecraft capabilities.

2.4. Orbit and operation mode

Euclid will operate at the Sun-Earth Lagrangian point L2, following a yearly orbit with a libration within ± 0 . ° 41 $ {\pm}0{{\overset{\circ}{.}}}41 $ across the ecliptic plane (Fig. 4). The Lissajous orbit is dynamically unstable and requires regular orbital maintenance, currently planned to last one day every four weeks (i.e. ∼3% of the total mission time). This orbit offers a very stable thermal environment and maximises the visible sky at any time.

thumbnail Fig. 4.

Euclid’s orbit and operation mode. Left panel: Euclid’s transfer orbit following launch at time T0, and the subsequent Lissajous orbit at the Sun-Earth Lagrangian point L2 (SEL2). Right panel: Euclid’s main step-and-stare operating mode, showing the steps as rotations around XSC.

Euclid employs a step-and-stare mode, acquiring data on a fixed sky field and then slewing to the next pointing. Slews come in two types, depending on the value of the eigenslew ϵ defined as the angle between one field quaternion2 and the next, which can be decomposed in an arc connecting two separate pointings on the sky plus a rotation around ZSC). Eigenslews ϵ 3 . ° 6 $ \epsilon\leq 3{{\overset{\circ}{.}}}6 $ are considered ‘small slews’, and ϵ > 3 . ° 6 $ \epsilon > 3{{\overset{\circ}{.}}}6 $ are ‘large slews’. Euclid has adopted a specific hybrid on-board AOCS architecture, where slews are performed using four reaction wheels in stop and go mode, and pointing stabilisation is achieved using low-noise cold gas micro-thrusters. With this solution, slew time is reduced thanks to the high torque provided by the reaction wheels. Slews only consume cold gas for the tranquilisation transient phase, large slews requiring longer tranquilisation periods. Therefore, Euclid’s lifetime slew budget is limited to 950 large and 2.5 × 105 small slews. The latter are weighted in the budget: for ϵ 1 . ° 2 $ \epsilon\leq 1{{\overset{\circ}{.}}}2 $ they count as a single slew, and for 1 . ° 2 < ϵ 3 . ° 6 $ 1{{\overset{\circ}{.}}}2 < \epsilon\leq 3{{\overset{\circ}{.}}}6 $ a penalty occurs proportional to ϵ. The slew constraints imply that the EWS must be implemented mostly with small slews (preferably with ϵ 1 . ° 2 $ \epsilon\leq 1{{\overset{\circ}{.}}}2 $), and that fields observed consecutively in time must be spatially adjacent.

2.5. Transits and visibility

The SAA and AA ranges define how much the spacecraft can deviate from observing at ‘transit’ meridian, which are the two ecliptic meridians defined by the perpendicular to the spacecraft’s XSC axis (for a transit SAA = 90°). The size and geometry of the Sun shield limit the SAA and AA ranges that can be used for observations. The ranges (and variations in) SAA and AA are constrained further by the fact that Euclid needs great thermal stability to minimise temporal point spread function (PSF) variations. The SAA limits allow the telescope to ‘de-point’ (i.e. to rotate around YSC) from transit to a maximum of 3° towards the Sun (SAA = 87°), and up to 20° away from the Sun (SAA = 110°), while the AA limits allow the telescope to rotate around ZSC up to |AA|≤5°.

In addition, the orbit libration mentioned in Sect. 2.4 imposes an additional 0 . ° 41 $ 0{{\overset{\circ}{.}}}41 $ buffer for the spacecraft orientation angles with respect to the Sun. This decreases the allowed ranges of both SAA and AA by 0 . ° 41 $ 0{{\overset{\circ}{.}}}41 $ on each side of their range intervals.

The allowed ranges of SAA and AA define the instantaneous sky visibility, shaped along the full circle defined by the two meridian transits (see Fig. 5). The orbit progresses with Euclid’s revolution around the Sun, continuously changing the visible sky enabling a full sky survey. Given the symmetry of the transit meridians, the spacecraft has access to the same region of the sky every six months by pointing in the ‘leading’ direction (towards the direction of motion around the Sun) or, flipping the telescope, six months later pointing in the ‘trailing’ direction.

thumbnail Fig. 5.

Strip of sky observable by Euclid at a given epoch (in light red). The dashed blue line is the great circle orthogonal to the Sun with the satellite observing orthogonal to the Sun for a transit at 0 deg. The asymmetry of the red region with respect to the transit blue lines is due to the asymmetry of the usable SAA range with respect to 90° (wider when pointing away from the Sun).

Figure 5 shows an example of the instantaneous sky visibility for a generic transit. It is evident that two small regions located at the ecliptic poles have perennial (continuous) visibility, whereas the lowest ecliptic latitudes can be observed only when crossed by a transit meridian. In practice, at any given time (or from a given position in the orbit) Euclid can scan an annulus on the sky, and consequently most of the sky must be observed at or close to transit. Observations that require long and regular visibility (such as for the EDS) can only be fulfilled in a very limited area on the sky at high ecliptic latitudes.

In general, a de-pointing (SAA ≠ 90°) induces a rotation of the focal plane with respect to the transit ecliptic meridian. To counterbalance this, the spacecraft must rotate around the +ZSC axis to keep the alignment with the transit meridian. The amplitude of this rotation must stay within | AA | 4 . ° 59 $ |\mathrm{AA}|\leq 4{{\overset{\circ}{.}}}59 $ to fulfil the thermal and orbit libration constraints. This effect becomes larger with increasing ecliptic latitudes and requires ad hoc solutions for the scheduling (see Sect. 7.4.2).

3. The instruments

3.1. VIS instrument

The VIS instrument contains a focal plane array (FPA) that consists of 6 × 6 Teledyne e2v charge-coupled devices (CCDs; 4k × 4k pixels each) with a Nyquist driven plate scale of 0 . 1 $ 0{{\overset{\prime\prime}{.}}}1 $ pixel−1, yielding a field-of-view (FoV) of 0.56 deg2, including detector gaps (for details, see Table 1 and the left panel of Fig. 6). In addition to the gaps, the central four rows in each detector serve as charge injection lines (which cannot be read out).

thumbnail Fig. 6.

VIS and NISP focal plane arrays. Left panel: VIS FPA, illustrating the detector placement. The dashed line shows the joint FoV of both instruments. Two narrow strips at the extremes of the Y axis are outside the joint FoV. Right panel: NISP FPA, with two narrow strips at the extremes of the X axis outside the FoV (the reference frame is XFoV − YFoV; see also Fig. 2 and the definition at the end of Sect. 2.1).

Table 1.

Sizes of the VIS and NISP focal planes and their corresponding FoVs, as well as the angular sizes of detectors and detector gaps and the width of the VIS charge injection lines (see Fig. 2).

The VIS is optimised to detect spatially resolved images of galaxies in the 550−900 nm passband (hereafter referred to as the ‘VIS band’; Cropper et al. 2014, 2016). The VIS nominal survey images (dithered; see Sect. 4) will have at least a S/N of 10σ for extended sources at the detection limit of IE = 24.5 AB mag, the average S/N being 15.8σ (see Sect. 5.2.2). We note that an extended source has a full width at half maximum (FWHM) such that F W H M gal > 1.25 F W H M PSF 0 . 225 $ FWHM_{\mathrm{gal}} > 1.25\,FWHM_{\mathrm{PSF}} \simeq 0{{\overset{\prime\prime}{.}}}225 $. This feature will enable accurate galaxy shape measurements for an average of 30 arcmin−2 galaxies over the survey area (Laureijs et al. 2011; Massey et al. 2013). Besides WL shape measurement, VIS data are also used to improve photo-z estimation, by enabling optimal photometric extraction of the less resolved, complementary ground observations thanks to its diffraction-limited image quality. To maximise the S/N for the shape measurements, the VIS band IE is rather broad (see Fig. 7), encompassing the Sloan Digital Sky Survey (SDSS) r and i bands, and the bluer half of the z band.

thumbnail Fig. 7.

Total sensitivity of Euclid’s photometric IE, YE, JE, and HE, and spectroscopic BGE, RGE (solid and dashed lines, respectively). The sensitivity (in electrons per photon) includes all optical surfaces as well as the detectors’ average quantum efficiency, all considered at their expected EOL performance after six years at L2 (degraded by radiation damage and contamination; see e.g. Venancio et al. 2020).

The VIS central data processing unit constructs the images from the pixel data and compresses them in a lossless manner in approximately 250 s No additional image processing will be done on board to maintain full control over systematic errors. The data will be transferred to the ground with a rate of approximately 520 Gbit/day (Racca et al. 2016).

3.2. NISP instrument

NISP (Maciaszek et al. 2014, 2016) contains an array of 4 × 4 HAWAII-2RGs detectors (2k × 2k pixels each) with a plate scale of 0 . 3 $ 0{{\overset{\prime\prime}{.}}}3 $ pixel−1, under-sampling its diffraction limited PSF (Fig. 6, right panel). Table 1 shows the size of the FPA, FoV, and gaps between the detectors. It should be noted that in the Y direction of the focal plane, the central gap ( 86 . 1 $ 86{{\overset{\prime\prime}{.}}}1 $ wide) is narrower than the two outer gaps ( 101 . 4 $ 101{{\overset{\prime\prime}{.}}}4 $ wide).

NISP is designed to carry out slitless spectroscopy (NISP-S) and imaging photometry (NISP-P) at NIR wavelengths (see Fig. 7). By using its grism and filter wheel assemblies (GWA and FWA, respectively), NISP can switch between slitless spectroscopy and imaging modes, which are detailed in the following. NISP will transfer data to the ground with a rate of approximately 290 Gbit/day, for a total of 810 Gbit/day, smaller than the spacecraft allocation of 850 Gbit/day.

3.2.1. Slitless spectroscopy

Euclid has a ‘blue’ grism (BGS000) covering the 0.92−1.25 μm wavelength range, and three ‘red’ grisms (RGS000, RGS180, and RGS270) covering RGE band, 1.25−1.85 μm (Costille et al. 2016). The blue grism, covering the BGE band, is not used for the EWS observations, and will be only employed for part of the EDS. The numeric labels indicate the dispersion directions, offset by 90° for the red grisms. Different dispersion directions are required to disentangle the spectra of various objects in the slitless spectroscopic exposures of the EWS. Due to a non-conformity discovered in 2020 (Laureijs et al. 2020), the RGS270 will not be used in the survey observations. Instead, the RGS000 and RGS180 will be rotated in the ROS by −4° and +4°, respectively (see Sect. 4.1 for details).

The red grisms disperse the light with nearly constant spectral resolution of 1.354 nm px−1, which gives ℛ = λλ ∼ 450 for an object of diameter 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $. This is larger than the minimum required value of 380 to achieve an error on the measured redshift of σ(z) < 0.001(1 + z). The spectroscopic observations support the GC probe and are optimised to detect the redshifted Hα emission of galaxies at z = 0.9 − 1.8. With a detection limit of 2 × 10−16 erg s−1 cm−2 (3.5σ) for a typical source of size 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $ at 1600 nm (see Sect. 5.2.2), NISP should be able to determine spectroscopic redshifts for at least 1700 galaxies deg−2 on average in the corresponding wavelength range 1250−1850 nm. This estimate, however, strongly depends on the intrinsic luminosity function of Hα emitters, which is still uncertain in the redshift range observed by Euclid (Pozzetti et al. 2016). Because the redshift is based on an emission line, passive galaxies will be underrepresented in the spectroscopic sample, with a bias against dense environments.

3.2.2. Photometry

Photometry will be measured for objects down to a minimum of mAB = 24 AB mag for 5σ point-like source in the JE, YE, and HE passbands. The photometric data support the GC experiment by providing the reference images needed to extract the spectra in the (slitless) dispersed images. The NIR photometric data, however, critically complement the ground-based observations (Sect. 5.5) in getting accurate photometric redshift estimates, at the primary probe level mainly needed for the WL experiment and essential for many other astronomical science aspects.

3.3. The Euclid joint field of view

The intersection of the VIS and NISP FoVs defines the Euclid joint FoV, with the X and Y dimensions defined by VIS and NISP, respectively. The Euclid FoV is 0.53 deg2. Its borders are shown by the dashed lines in Fig. 6, resulting from the overlap of the VIS and NISP FoVs aligned on an edge. The left and right edges of the NISP FoV and the top and bottom edges of the VIS FoV are outside the joint Euclid FoV.

4. Observing and dithering sequences

4.1. The reference observation sequence (ROS)

Euclid executes a highly optimised ROS (see Fig. 8) at every survey field, exploiting the instruments’ inter-operability. The ROS visits four nearby pointings at every field, covering an area of 0.53 deg2 (see Sect. 3.3) common to both instruments and fulfilling the galaxy number density and S/N requirements detailed in Laureijs et al. (2011). Small ‘dither slews’ are performed between the pointings, taking 66 s.

thumbnail Fig. 8.

(Simplified) ROS executed for each field, observing four dithered pointings. Each pointing results in a IE, YE, JE, and HE image and one spectral exposure, RGE. Filter and grism wheel movements are shown, together with stabilisation times (‘Stab’). Not shown are wheel movements during the dither slew, readout and processing done simultaneously with other operations and slews, overheads, and various inline calibrations.

At each pointing, VIS takes an image and NISP a simultaneous spectral exposure with the red grism, both lasting about 570 s (we note that these times are not frozen yet). Once the VIS shutter is closed, the GWA and FWA move for the three NISP images of 112 s each. A 2 s margin is allocated between the end of a NISP exposure and the wheel actuation, ensuring the NISP exposure is completed before the wheel is moving to avoid compromising the last frame. Moreover, a stabilisation time of 10 s is considered between a wheel movement and the following exposure. During the NISP imaging, VIS takes biases, flats, and other calibration frames; these ‘parallel’ observations are highly synchronised and optimised so that the instrument operations do not interfere with each other. In addition, VIS also takes a shorter science exposure of 108 s during the HE exposure in the first pointing, in order to help with the PSF dynamic range on relatively bright stars that saturate during the standard, longer exposures. Details of these sequences are given in, for example, Cropper et al. (2016) and Maciaszek et al. (2016).

After each pointing a dither step is applied and a new grism position is selected. The ROS uses the RGS000 and RGS180 at two angles each, offset by four degrees, to allow for sufficient decontamination of the overlapping slitless spectra.

The total duration of the ROS, including dither slews and overheads is 4214 s. At the end of the ROS a slew towards the next field is performed. Most of these slews are small ( ϵ 3 . ° 6 $ \epsilon\leq 3{{\overset{\circ}{.}}}6 $) and referred to as ‘field slews’. The slew duration is a function of the (eigen-) angular rotation. On average it is 182 s, implying a total length of 4396 s for the ROS (slews included). This is less than the upper limit of 4400 s as defined at mission system level during budget allocation. On occasion, the Euclid Reference Survey (ERS) requires larger slews that are limited to a maximum number of 950 over the full mission. In the most recent EWS solution (see Sect. 8), the ‘large slews’ comprise 1% of all non-dither slews applied.

4.2. The dithering strategy

For each field, the ROS obtains multiple exposures with dithered pointings to mitigate detector defects and cosmic rays, and to meet the required depth. The depth will vary across the field, not only because of masked defects, but predominantly because the NISP and VIS focal planes have different detector and gap sizes (see Fig. 6, and Table 1). The dithering strategy between pointings, used for the ERS, must meet the following requirements: (i) 95% of the survey area shall be covered with at least three exposures in IE. 90% of the survey area shall be covered by at least three exposures in each of the NISP bands YE, JE, HE; (ii) 90% of each survey field shall be covered by three or more RGE spectroscopic exposures, and 50% by four or more RGE spectroscopic exposures (using different grism orientations); and (iii) the NISP imaging of the fields covered by the NISP spectroscopic channel in the EWS shall be acquired over the whole image with depth on average fainter than mAB = 24 mag for 5σ point-like source.

To analyse the sky coverage of a given dither pattern, we simulated the ROS observations at the pixel level for a nearly square sky region. The latter covers several FoVs in each of the two dimensions (N × N joint FoVs, where 3 ≤ N ≤ 8) to avoid boundary effects. The pixel count statistics (number of exposures per sky area) are then computed with an integration time map calculator.

In the various dither patterns, the minimum dither step in each of the two directions is sized so as to ensure that the slew is larger than the largest of the detector gaps in that direction (which, according to Table 1 are the NISP gaps). The maximum dither step is sized to prevent gaps of one line or row of detectors to overlap with the next gaps after the three dithers. In our simulations, we used values for the detector gaps reflecting the ones present in the as-built instruments (see Table 1). For VIS, we used 13 . 6 $ 13{{\overset{\prime\prime}{.}}}6 $ in the X direction, 67 . 6 $ 67{{\overset{\prime\prime}{.}}}6 $ in the Y direction and 0 . 4 $ 0{{\overset{\prime\prime}{.}}}4 $ for the width of the charge injection lines, while for NISP we used 50″ in X and 100″ in Y. Subsequent fields are shifted by about 0.7°, setting in the simulation an overlapping of 1% (in area) between contiguous FoVs. This constitutes the basic set-up of our simulations.

The dither step size is affected by errors, namely after a dither slew there is a required absolute pointing error (APE) smaller than 7 . 5 $ 7{{\overset{\prime\prime}{.}}}5 $ at 3σ in both XFoV and YFoV, for a total displacement error of 11″. Also the overlap between contiguous FoVs is also affected by APE as well. The line of sight has an uncertainty off-set introduced by the rotation of the filter and grism wheels (that is compensated by the spacecraft AOCS).

We produce coverage maps in two ways. The first method is deterministic (method D). In this case, we consider the basic set-up and add a deterministic small shift (of 11″) to each step of the dithering pattern. The implementation of this further displacement is necessary to be safe in filling the gaps when considering the estimate of the uncertainties in the dithering step. No random errors are considered in this method. In the second method, the Monte Carlo (MC) method, we consider uncertainties on dithering steps and directions (a pointing error of 11″ at 3σ), on off-set repeatability, and other relevant parameters. In this method, the actual overlapping between contiguous FoVs is affected by the APE uncertainty; we simulate an APE with a Gaussian random amplitude of 3 . 5 $ 3{{\overset{\prime\prime}{.}}}5 $ at 1σ and uniform random orientation. Monte Carlo iterations are then run on a representative patch to extract the coverage maps and the statistics of pixel numbers for the desired area inside the simulated patch for each of the dither patterns presented in Sects. 4.2.1 and 4.2.2.

4.2.1. Minimal pattern

The dither pattern defined in Laureijs et al. (2011) was minimal in the sense that it only respected the stringent constraints on the size of a dither step. The latter is constrained to a minimum of 100″ by the reaction wheels to prevent mechanical damage to wheel ball bearings in small rotation regime, and to a maximum of 396″ by the size of the star catalogue available to the AOCS. The theoretical pattern (referred to as a ‘J’ pattern given its shape; see Fig. 9) is defined as follows, starting from the pointing of the first exposure: (i) dither step 1: ΔXSC = 50″, ΔYSC = 100″; (ii) dither step 2: ΔXSC = 0″, ΔYSC = 100″; (iii) dither step 3: ΔXSC = 0″, ΔYSC = 100″.

thumbnail Fig. 9.

Seven dither patterns that we have analysed shown in the spacecraft reference frame (see Fig. 2). For clarity, a shift of 5″ (or of −5″) in the Xsc direction is applied to the ‘J’ (or to the ‘N’) pattern in the figure so as to distinguish it from the ‘S’ pattern. Only the ‘J’ and the ‘S’ patterns, the first two steps of which are identical, are fully compliant with the stringent survey requirements.

These are the minimum values that need to be ensured by the AOCS. As discussed earlier, in order to account for the error in the dithering step, the dither is commanded to be 11″ larger than the minimum step.

In method D, uncertainties in the telescope pointing were not considered. In total, we simulated 3 × 3 joint FoVs (i.e. nine adjacent executions of the ROS). For a simple visualisation of the integration time map, we simulated the sky coverage with the resolution of the VIS and NISP FPAs degraded to 1″. The width of the VIS charge injection lines were increased to 1″ to be included in the simulation. The top panels in Fig. 10 show, for VIS and NISP, the central 50 . 0 × 56 . 7 $ {\simeq}50{{\overset{\prime}{.}}}0\,{\times}\,56{{\overset{\prime}{.}}}7 $ of the ‘J’ integration time map, that is, a single survey field and its boundaries (containing ≃107 simulated pixels).

thumbnail Fig. 10.

Illustration of the four stacked integration time maps for ‘J’ (top) and ‘S’ (bottom) dithering patterns. The left half displays the maps for VIS and the right half those for NISP. The colour scaling gives the number of stacked exposures per map pixel. Details are given in the text.

In the MC method, the statistics on the number of exposures per map pixel were evaluated with the same tool, but now considering random errors on some pointing parameters. We simulated a larger sky area (7 × 7 joint FoVs) to allow a corresponding increase in the number of realisations and a better characterisation of the overlapping regions between the FoVs. Here, we adopted a resolution of 0 . 2 $ 0{{\overset{\prime\prime}{.}}}2 $ to accurately quantify the effect of the finer geometric characteristics of the VIS and NISP FoVs. Statistics extracted for areas of different sizes (for example of 4 . ° 0 × 4 . ° 0 $ 4{{\overset{\circ}{.}}}0 \times 4{{\overset{\circ}{.}}}0 $) around the central joint FoV show that the coverage requirements are fulfilled at a global level. On a field-to-field basis, the coverage can vary within the associated standard deviations.

4.2.2. The baseline ‘S’ dither pattern

Looking for a potentially better dither pattern, we analysed six other patterns with four vertices (Fig. 9). Two cases were proposed by Arendt, Kashlinski, and Mosley based on previous experience (priv. comm.; labelled as AKM2 and AKM3), and four cases (‘S’, ‘R’, ‘N’, and ‘X’) were taken from Markovič et al. (2017), who focused on optimising the performance of the spectroscopic survey only, without taking into account the concurrent constraints required by the imaging part. The two AKM patterns exceed the maximum step size, with only AKM3 satisfying the constraint on the joint visibility.

Among the four possibilities suggested in Markovič et al. (2017), only the ‘S’ pattern, which is the one closest to the reference ‘J’, meets all constraints. Its statistics on the number of exposures per pixel derived from simulations performed considering the most relevant sources of uncertainty are given in Table 2. The ‘S’ pattern improves, particularly for NISP, upon the ‘J’ pattern, decreasing the fraction of pixels with a single exposure (from 0.59% to 0.23% in the case of VIS and from 3.47% to 0.99% in the case of NISP) while increasing the fraction of pixels with two exposures (from 3.91% to 4.34% in the case of VIS and from 4.79% to 8.13% in the case of NISP). For the ‘S’ pattern, the percentages of X ≥ 3 covering are 95.43 ± 0.03 and 90.88 ± 0.05 for VIS and NISP, respectively, while the percentage of X ≥ 4 covering for NISP is 50.42 ± 0.06 (the statistics reported in Table 2 are extracted from 500 realisations). Thus, at least under the considered uncertainty in specifications, the requirement 1 of Sect. 4.2 for VIS and NISP imaging is satisfied at ∼14.3σ level and at ∼17.6σ level, respectively, while the requirement 2 of Sect. 4.2 for NISP spectroscopy is satisfied at ∼17.6σ level and at ∼7σ level for the X ≥ 3 covering and X ≥ 4 covering, respectively. We note that the standard deviation of the cumulative case (X ≥ 3 or X ≥ 4) covering cannot be derived by a simple analytical propagation of the standard deviations of simpler cases because of the presence of mutual correlations. These, however, are taken into account in the numerical simulation.

Table 2.

Statistics (mean and standard deviation) on the number, X, of exposures per pixel for the ‘S’ dither pattern: percentages of covering for individual (X = 0, 1, 2, 3, 4) and cumulative (X > 4, X ≥ 3, X ≥ 4) bins.

The ‘S’ pattern was therefore chosen as the baseline dither pattern for the ERS (see also Markovič et al. 2017). The coverage map is shown in the bottom panels of Fig. 10, and in Fig. 11 we display the ‘S’ pattern as it appears on sky. It is defined as follows: (i) dither step 1: ΔXsc = 50″, ΔYsc = 100″; (ii) dither step 2: ΔXsc = 0″, ΔYsc = 100″; (iii) dither step 3: ΔXsc = 50″, ΔYsc = 100″. However, the random errors present in pointing could cause the actual shifts to be smaller than nominal and so be insufficient to cover the gaps. Therefore, in order to be on the conservative side, to each move of this nominal ‘S’ pattern a buffer of 11″ is added to each ΔXsc and ΔYsc, resulting in the operational ‘S’ which has: (i) dither step 1: ΔXsc = 61″, ΔYsc = 111″; (ii) dither step 2: ΔXsc = 11″, ΔYsc = 111″; (iii) dither step 3: ΔXsc = 61″, ΔYsc = 111″. We note that while the latter sequence has been used in the simulations to get the statistics presented in Table 2, the ERS presented in Sect. 8 was built with the nominal ‘S’ pattern.

thumbnail Fig. 11.

Adopted ‘S’ dither pattern projected on the sky for a leading (left) and a trailing (right) pointing. See Fig. 2 for reference frames. The red frame is the first to be observed, and the dots mark the centre of each frame. We note how the pattern on the sky is invariant if the telescope is flipped.

We note that by increasing the number of dithers and their step size, a more uniform coverage can be achieved (see e.g. Rowe et al. 2011). However, this needs to be balanced against the total survey area, the mission duration, readout overheads etc. Other patterns can still be studied and implemented in case of a revision of the current survey and hardware limitations.

5. The EWS region of interest

The constraints that drive the implementation of the Euclid survey areas fall into three main categories (see Sect. 1): (i) environmental; (ii) calibration needs; and (iii) spacecraft constraints.

Environmental properties are external physical constraints of the observable sky, namely zodiacal background, Galactic extinction, foreground stellar density and blinding stars. The zodiacal background and extinction affect the S/N, as does the stray light from bright stars and Galactic plane. The star density also affects the number of useful objects, because their light is dispersed and overlaps with the galaxy spectra of interest, or affects the ability to measure shapes reliably. Moreover, an increased star density also increase the multiplicative bias in weak lensing shear estimates, if unaccounted for (Hoekstra et al. 2017).

Together, these constraints lead to the definition of the EWS RoI, a collection of four contours enclosing two larger ‘mainlands’ and two smaller ‘islands’. Each contour is defined by a series of joint segments derived from the zodiacal light (ecliptic latitude segment), dust extinction (Galactic caps segment), or the stellar density (Galactic latitude segment).

5.1. Environment models

The cosmological measurements are strongly mediated by the nature of the areas of the sky selected for the survey. In the following, we discuss the models of the environment used in the survey planning.

5.1.1. Zodiacal background

A large amount of literature is available on this subject and several models and estimates have been proposed over the years. Updating previous work done by the DUNE consortium (Refregier & Douspis 2008), we used a lean and conservative model of the zodiacal background obtained by combining the spectral dependence proposed by Aldering (2001) for the proposal for the SNAP satellite (Aldering et al. 2002) and the angular dependence found in Leinert et al. (1998). This ‘basic’ model assumes a cylindrical symmetry with respect to the Sun. In this time-invariant and symmetrical model, the zodiacal background flux density ζ depends solely on wavelength λ and ecliptic latitude β,

ζ ( λ , β ) = f λ ( λ ) g ( | β | ) . $$ \begin{aligned} \zeta (\lambda , \beta ) = f_\lambda (\lambda ) \, g(|\beta |). \end{aligned} $$(2)

For the spectral dependence, Aldering (2001) suggests

f λ ( λ ) = { κ , λ [ 0.41 μ m ; λ ref ] κ 10 0.730 ( λ λ ref ) / μ m , λ [ λ ref , 2.2 μ m ] , $$ \begin{aligned} f_\lambda (\lambda ) = {\left\{ \begin{array}{ll} \kappa ,&\lambda \in [0.41\,\upmu \mathrm{m}; \lambda _{\rm ref}] \\ \kappa \,10^{-0.730\,(\lambda -\lambda _{\rm ref})/\upmu \mathrm{m}},&\lambda \in [\lambda _{\rm ref},2.2\,\upmu \mathrm{m}] \end{array}\right.}, \end{aligned} $$(3)

with λref = 0.61 μm. The normalisation constant κ fixes the flux density to the north ecliptic pole (NEP; β = 90°), such that κ = 1.76 × 10−18 erg cm−2 Å−1 arcsec−2, corresponding to mAB = 23.05 mag arcsec−2. For comparison, this is just 3% lower than the value given by Leinert et al. (1998) for the NEP at 0.50 μm.

For the dependence on ecliptic latitude β, we have g(|β|) as a dimensionless, monotonically decreasing function over the interval β ∈ [0° ,90° ], normalised to g(90° ) = 1. Leinert et al. (1998) report values for ζ (their Table 17) as a function of β and elongation from the Sun, for a wavelength of 0.50 μm. We reproduce their values for an elongation of 90° (applicable to Euclid) in Table 3, showing that at β = 20° (10°) the zodiacal background is 2 times (3 times) higher than at the NEP. This dependence on latitude is in good agreement with values measured by the Solar Mass Ejection Imager (SMEI) satellite (Buffington et al. 2016). Euclid uses fixed integration times and thus does not compensate for increased background. We therefore limit the EWS to |β|≥10°, corresponding to a reduction of 20/30% in VIS/J S/N compared to the NEP when taking into account also the stray light (see Sect. 5.1.3 and Fig. 19). This still allows for a suitable number density of detected galaxies for WL and GC averaged over the survey.

Table 3.

Ecliptic latitude dependence of the simple zodiacal background model of Leinert et al. (1998).

This basic model results in a background that is constant in time, depends on the ecliptic latitude as in Leinert et al. (1998), and follows an exponential decay for λ ≥ 0.61 μm. In Table 4 we report corresponding numeric values and magnitudes for the various Euclid bands. Here we used the simple average wavelength of a band, defined as

λ = λ T ( λ ) d λ / T ( λ ) d λ , $$ \begin{aligned} \lambda _* = \int \lambda \,T(\lambda )\,\mathrm{d}\lambda /\int T(\lambda )\,\mathrm{d}\lambda , \end{aligned} $$(4)

Table 4.

Normalisation of the basic zodiacal model at the average wavelength of Euclid bands (see Fig. 7).

where T(λ) is the end-of-life (EOL) throughput (see Fig. 7). For background-limited observations, we could instead weight the integrals in Eq. (4) by the background flux density; this would slightly change the values of λ*, as we will discuss in a separate paper.

Better fits to the dust emission detected by the COsmic Background Explorer (COBE) Differential Microwave Radiometer, can be obtained by considering more sophisticated models, such as the ones of Kelsall et al. (1998). These models include a slab model that is not centred on the ecliptic plane and add inhomogeneous clouds. These result in a zodiacal background that depends not only on the wavelength and direction of observation, but also on the observation’s epoch. In such models the minimum background no longer coincides with the ecliptic poles. Instead, it circles the poles with a yearly period, and the value at the NEP has a corresponding sinusoidal variation with a 20% peak-to-peak variation (see also Pyo et al. 2012).

We implemented such a model (Maris et al. 2006, 2016), which includes time dependence, differences in trailing and leading directions and possible deviations of the pointing directions from the orthogonal direction to the Sun. The model is evaluated for 103 random fields of the reference survey binned in ecliptic latitude, and shows scattering due to different observation epochs. It predicts for all bands on average 15 − 20% lower background values than the basic model, and 30% less flux at the NEP (see Fig. 12). We use this more complex model for more precise estimates done a posteriori once a reference pointing solution is obtained using the basic model as an input.

thumbnail Fig. 12.

Comparison of the zodiacal emission at 1.25 μm. For a random subset of the pointings in the reference survey, we show results from the conservative basic model (red triangles) and the estimates by the more complex model (blue diamonds). The latter (Maris et al. 2016) depends on the observation time and angular direction, and it no longer has cylindrical symmetry. This causes the scatter in the blue points. The EuclidYJH bands show a similar behaviour.

We adopted the basic model as our baseline to define the ecliptic latitude exclusion zone presented in Sect. 5.2, because it is more conservative, whilst providing reasonable margins. This model is also the reference model adopted for Mission Performance Evaluations, and hence Euclid might detect a slightly larger number density of usable galaxies than our current predictions. We note that among the different models in the literature some could yield background values larger than the adopted basic model, because of a different normalisation (see e.g. Wright 1998).

5.1.2. Galactic extinction

Extinction by interstellar dust is estimated from the E(B − V) reddening maps3 produced by the Planck Collaboration XI (2014). The map’s resolution of 5′ is high compared to the linear size of Euclid’s FoV of 0 . ° 7 $ {\sim}0{{\overset{\circ}{.}}}7 $ ∼44′. We smoothed the map with a 2 deg wide Gaussian kernel, such that the segments of the RoI boundary that are determined by E(B − V) have a comparable smoothness as the segments determined by Galactic and ecliptic latitude.

To apply rigidly the original E(B − V) = 0.08 mag limit of Laureijs et al. (2011) would cause highly convoluted RoI (see Sect. 5.2) contours and holes inside the contiguous survey areas. To achieve a larger RoI with compact regions, a slightly adjusted upper limit of E(B − V) = 0.09 mag was chosen, while allowing local excursions up to 0.17 mag to simplify contours and avoid local holes. These settings define a first version of the Galactic exclusion zone. The introduction of the ecliptic exclusion zone (Sect. 5.1.1) in the two resulting Galactic caps, divides them into two larger mainlands and two smaller islands, as discussed in Sect. 5.2. The median value of E(B − V) over the RoI is 0.037 mag (more statistics are given in Sect. 5.2.2). This approach meets the performance requirement on mean galaxy number density, while preserving a connected survey that optimally complements ground-based data (Sect. 5.5).

We also use the smoothed E(B − V) to estimate the S/N in the RoI. To this end we must compute the total extinction in magnitude for an Euclid band of central wavelength λ:

2.5 log 10 ( F obs / F 0 ) = A λ 3.1 E ( B V ) . $$ \begin{aligned} -2.5\,\mathrm{log}_{10}{(F_{\rm obs}/F_0)} = A_\lambda \; 3.1 \; E(B-V). \end{aligned} $$(5)

Here, Fobs is the observed flux and F0 is the flux in the absence of extinction. The total extinction in the V band is quantified by AV = RVE(B − V) where RV = 3.1 parameterises the dust extinction in our Galaxy. We infer the extinction scaling coefficient Aλ with respect to the V band from the dust extinction curves of Gordon et al. (2003) for the Euclid channels (IE, YE, JE, HE, RGE) based on their central wavelength. Results are given in Table 5 for the NEP with E(B − V) = 0.07 mag.

Table 5.

Extinction scaling coefficients Aλ for the Euclid passband central wavelengths and total extinction A at the NEP.

5.1.3. Bright sources and stray light

The EWS allows for 15% of its area to be lost due to various effects such as dead pixels, cosmic rays, etc. This masking budget allocates 2% for bright stars and an additional 2% for their ghost images. Stars brighter than mAB ∼ 17.5 mag for VIS and 16.0 mag for NISP will saturate the detectors for the baseline integration times, and we refer to them as ‘bright’ stars.

Due to image persistence constraints for NISP, stars with mAB ≤ 4 mag must never be present within the FoV for all three photometric bands. We also apply this rule to the VIS instrument due to stray light considerations as discussed below. In the following, we refer to both cases as ‘blinding’ stars. We selected three catalogues in the literature that include all the brightest stars in the sky and sample the wavelength domain of the two instruments: i band for VIS from the ATLAS All-Sky Stellar Reference Catalog (Tonry et al. 2018), Y band for NISP from the spectrally matched magnitudes of stars from the Tycho2/2MASS catalogues (Pickles & Depagne 2010), and J and H band for NISP from the 2MASS Point Source Catalog (Skrutskie et al. 2006). In total, there are 1034 blinding stars within the RoI (Fig. 13). 1003 of those stars were first selected through the Y < 4 or J < 4 or H < 4 AB mag NISP criteria, and 31 stars with Y, J, H > 4 mag were added through the i < 4 mag VIS criteria (there are 275 i < 4 mag stars within the RoI, the majority already selected through the NISP criteria). The blinding star catalog will be refined when color terms for conversion to the Euclid bands will be available. Survey fields with blinding stars will be skipped on all instruments, amounting to 3% of the RoI area. Bright (non-blinding) stars will be observed but locally masked during data reduction. The observed area lost in this way is accounted for in the pre-allocated masking budget, which does not incorporate areas skipped on purpose because of the presence of blinding stars.

thumbnail Fig. 13.

1034 blinding stars present in Euclid’s RoI. There are 1003 stars over the RoI that have Y, or J, or H < 4 AB mag (in red), and 31 stars with i < 4 AB mag and all three bands Y, J, H > 4 AB mag (in blue). The size of the dots scales down with the magnitude from −2 to +4. Those nearby bright stars are evenly distributed over the RoI, uncorrelated with the Galactic latitude. The three EDFs (see Sect. 6.2), free of blinding stars, are shown in green.

In the following we show how stray light from stellar sources inside and outside the FoV affects the observations and dictates the RoI Galactic latitude threshold. Stray light is generated in various ways, such as reflections on mechanical structures, scattering on contaminated optical surfaces and their intrinsic surface roughness, by diffraction on the edges of mechanical parts such as baffles, spider arms, or by multiple internal reflections between optical elements. Stray light creates an additional background contribution that can be diffuse or structured (ghosts).

Up to 2014, the stray light was expected to be a minor fraction of the total diffuse background (this was initially specified to be less than 20% of the zodiacal background at the NEP). However, with lessons learned from the Gaia mission, a study on stray light contamination was carried out by ESA and the industry for Euclid. Venancio et al. (2020) have studied the stray light aspects extending the analysis from pure in-field (dominated by the mirrors particulate contamination) to the far out-of-field domain (dominated by the internal structural multiple reflections). Both can contribute significantly to the stray light level, as out-of-field stray light, though largely attenuated by diffuse scattering process, integrates over the full sky and becomes dominant when getting close to the Galactic plane. They found that particulate contamination on the mirrors will be the main contributor: stray light in some sky areas can become comparable and even exceed the local zodiacal background. Consequently, Euclid must stay clear of the Galactic plane, since sources both inside and outside the FoV (in-field and out-field stray light, respectively) contribute.

The normalised diffusion irradiance profile (NDI; Venancio et al. 2016) describes the profile of the scattered light in the telescope focal plane for a point source at a given position either within the FoV or up to 20° away from the optical axis. Euclid’s Korsch optical design (Korsch 1972) effectively suppresses scattered light. For Euclid and its enhanced baffling (Venancio et al. 2020), in-field stray light will impact the S/N of faint galaxies from our science goal defined limit of IE = 24.5 AB mag, causing them to fall below the minimal value of S/N = 10. To ensure such a S/N level is realised through the mission science pipeline involving all the steps of data processing and signal extraction, the system team in charge of scaling the mission design built margins by adopting a conservative goal of S/N = 14 using their own internal radiometric S/N metric. The S/N of a IE = 24.5 AB mag galaxy degrades from the mission median value of 17.5 to 14 if the galaxy is at an angular separation of 6.8, 4.1, 2.5, 1.5, and 1.0 arcmin from a bright star of IE = 4, 5, 6, 7 and 8 AB mag, respectively. The left panel in Fig. 14 illustrates the case for a IE = 4 star, the impacted area representing just 9% of the entire FoV, a truly remarkable performance made possible by the Korsch optical design. For fainter stars, of magnitude 9 and 10, the radius of S/N degradation to a level of 14 is reached at a 8 and 1.5 arcsec radius, respectively. As shown in the right panel of Fig. 14 the impact of these stars is negligible beyond the core of the PSF. The average density of the 8, 9, and 10 mag stars (Zakharov et al. 2013) over the RoI amounts to 0.6, 1.6, and 3.8 in-field stars per FoV for the VIS: their collected impact will be limited to a tenth of a percent of area loss of the Euclid FoV on average across the RoI. Accounting for the sparse 4 to 7 AB mag stars does not alter these statistics. A comparable performance is expected from the NISP as the mission design drove an NDI dominated by the telescope, not by the instruments.

thumbnail Fig. 14.

Impact from bright stars on the S/N achievable by a target galaxy IE = 24.5 AB mag (extended source) when placed over a full Euclid FoV ( 0 . ° 7 × 0 . ° 7 $ {\sim}0{{\overset{\circ}{.}}}7 \times 0{{\overset{\circ}{.}}}7 $, with a scheme of the 6 × 6 CCD mosaic layout superimposed) in a typical part of the Euclid RoI. The colour scale is the S/N of a target source placed in the field at that position. The degradation in S/N when the target galaxy is placed closer and closer to a bright star is evident. Left panel: IE = 4 AB mag star with effects on almost 50% of the FoV, but science impacted on less than 10% of the FoV. The solid yellow area is where the target galaxy is at S/N < 14 (see text). Right panel: common stars of AB magnitude IE = 8, 9, and 10, which have a very limited impact overall with respect to the FoV size: less than 0.1% of the FoV lost on average across the RoI.

In summary, for individual bright stars Euclid will skip observing tiles in which at least one of the four exposures will contain a mAB ≤ 4 mag star in any of the Euclid bands. For a handful of extremely bright stars also nearby tiles will be skipped, according to an avoidance radius for the tile centre set at a level of stray light yielding a 15% degradation in S/N (NDI model). Areas affected by in-field stars fainter than mAB = 4 mag will be masked during the data reduction phase, as well as ghosts originating from the dichroic.

All the stars outside the telescope FoV will also contribute globally to the background level of stray light. Their combined effect is to add a diffuse cumulative component that depends on the pointing direction of the telescope. In consequence, this effects scales with the Galactic latitude, the NDI defining an intensity ratio of the collected brightness of the Galaxy. The out-field stray light two-dimensional map adopted in this paper is a Besançon model of the Galaxy flux (Robin et al. 2012, 2014). The model is scaled at the relevant wavelength to match at the one percent level the out-field stray light level computed at 12 selected points across the whole sky. The levels were estimated by the system team for the 2018 mission critical design review (MCDR) effort.

A spacecraft stray light model based on the estimation of the NDI was then used by the system team to describe the stray light due to the diffusion in the telescope and gauge its effects on the PSF and local background level on the focal plane (Venancio et al. 2016, 2020). Assuming the entrance of the telescope is illuminated by a distant point source (collimated light), then the NDI is defined as the ratio of light irradiance (power per unit area) on the image plane to the source irradiance in object space at the entrance of the telescope. The NDI is computed for both VIS and NISP using the ASAP optical software4, a ray-tracing program that uses a statistical MC approach. The computation is done with telescope and instrument optical and mechanical models and associated contamination assumptions. Then the NDI is applied on the sky for a mesh of pointing directions on the sky (with sampling equal to the Euclid FoV) over the full sky in order to estimate for each possible pointing direction the cumulative out-of-field stray light maps.

Conservative estimates of these contaminants established by the mission system team are adopted for the background and associated noise computations presented in this paper. This drives in particular the RoI definition with respect to Galactic latitude (from a minimum of |β|≥23° to nominal |β|≥25°), and how close the EWS can get to the Galactic bulge, refining the Galactic exclusion zone.

Finally, an additional concern is that due to stray light: Euclid must avoid pointing within a circle centred on the position of Solar system planets, Mars and Jupiter having the largest radius, of 13 deg.

5.2. The Euclid region of interest

5.2.1. The four main areas of the RoI

The dominant factors that determine the RoI are the zodiacal background, the Galactic extinction, and stray light due to the Galactic stellar density. Minor contributors such as emission from Galactic cirrus were ignored, being at least five magnitudes fainter than the total background over most of the RoI (Sect. 5.3). The EWS will also skip fields containing blinding stars (Sect. 5.1.3), leaving only their faint effect imprinted on the out-field stray light.

The main outline of the RoI is defined by the extinction limits, an ecliptic latitude threshold of |β|≤10°, and a Galactic latitude threshold of |b|≤23°. We note that for declination δ ≥ +30° we set |b|≤25°, since the Euclid complementary ground surveys were designed and started in 2017 with this value from Laureijs et al. (2011). Section 5.5 describes how the RoI is affected by these ground-based surveys.

The RoI fragments into four quadrants delimited by the yellow contours in Fig. 15 and detailed in Table 6. The RoI is best presented on this equi-rectangular projection. We use elliptical projections when highlighting aspects of area conservation. Our sky projections were produced with IPAC’s Montage package5 and the University of Groningen’s Kapteyn package for Python6.

thumbnail Fig. 15.

RoI outline (17 354 deg2) with the accepted ranges of the stellar density, dust extinction, and zodiacal light.

Table 6.

Four separate quadrants of the RoI (see Fig. 17).

The environment limits and their impact on the RoI are shown in Fig. 15, while Fig. 16 shows the distribution of extinction, stellar counts from Gaia (limited at mag = 20) and zodiacal background within the RoI. In total, the present RoI encompasses 17 354 deg2 that are compliant for Euclid’s core cosmology science. The EWS can be constructed from any 15 000 deg2 within. Some parts of the EWS will inevitably be of lower quality for cosmology, yet their legacy value is high. For example, the Small Magellanic Cloud (SMC) is inside the RoI, although clearly at odds with the survey constraints. Pushing further into the Galactic plane, though, would rapidly reach extinction levels unacceptable for Euclid’s core science, as is evident from Fig. 17 that shows the combination of all constraints, highlighting the best parts of the Euclid sky.

thumbnail Fig. 16.

Distribution of reddening values (left), stellar density (middle), and zodiacal background (right) across the RoI. Below each graph, the box-and-whiskers plot marks the mean (red line), median (black line), and interquartile range (empty box), plus the minimum and maximum values.

From Fig. 17 it is clear that the area of the RoI changes with ecliptic longitude. Figure 18 shows the area of the RoI as a function of ecliptic longitude (in bins of 1° of ecliptic longitude). The plot contains two global maxima, which coincide with the longitudes that cross both a mainland and an island, and two global minima, which coincide with the intersection of the ecliptic and Galactic planes.

thumbnail Fig. 17.

Sky constraints as experienced by Euclid from L2, and the RoI avoiding the worst regions. The four quadrants of the RoI are in blue. The solid yellow lines trace the ecliptic latitude threshold of |β| = 10°. The solid red lines trace the galactic latitude threshold of |b| = 23°. For declination δ ≥ +30°, the threshold is |b| = 25° to match the northern ground surveys (Sect. 5.2.1).

Given that one degree of longitude is scanned by one day of orbit and that in one day ∼10 deg2 of EWS are observed, in a six year mission it is only possible to observe a maximum of 60 deg2 of EWS sky, per degree of longitude. In practice, this time must be shared with calibrations and Euclid deep field (EDF) observations that collectively take ∼20% of the total time. This lowers the EWS allocation to a maximum average value of 48 deg2 per degree of longitude. From the analysis of Figs. 17 and 18, it is clear that the sky in the RoI is not uniform enough to fill this quota. Given the limited pointing range of the telescope, with observations at or close to transit, this inevitably leads to the depletion of the available (i.e. yet unobserved) sky in some ecliptic longitudes. This reveals an intrinsic limitation to the maximum efficiency attainable by the EWS, in which in some parts of the year there will be unallocated time periods that increase in duration towards the end of the mission. This is an important feature of EWS solutions, as discussed in Sect. 8.2.

thumbnail Fig. 18.

Available RoI area as a function of the ecliptic longitude. The troughs are produced by the intersection of the Galactic and ecliptic planes and lead to unallocated time.

5.2.2. S/N and survey depth in the RoI

Figure 17 highlights the best parts of the Euclid sky. In the following, we compute the corresponding S/N maps, which provide the quantitative context.

For our S/N computations we take into account the following aspects at the hardware level: telescope and instruments’ internal backgrounds, photometric zero points (encoding the total throughput), read noise and dark current. These are independent of sky position and were taken from the latest available ground characterisation measurements. At the environmental level, we include all-sky maps for the zodiacal background, extinction, and stray light from the Galaxy as detailed in Sect. 5.1.

At the operational level, we allow for three exposures (VIS and NISP imaging) and four exposures for NISP spectroscopy, the FPA geometries, integration times, and the size of the measurement apertures. This is motivated by the fact that 90% (50%) of the survey area is covered with at least three (four) imaging exposures (see Table 1). The S/N measurement metrics are evaluated as follows: for VIS, we consider an extended source with a total magnitude of IE = 24.5 AB mag in a 1 . 3 $ 1{{\overset{\prime\prime}{.}}}3 $ diameter aperture, capturing 94% of the flux. For NISP photometry, we consider a point-like source with a total magnitude of mAB = 24.0 mag in the YE, JE, HE bands in a 0 . 9 × 0 . 9 $ 0{{\overset{\prime\prime}{.}}}9 \times 0{{\overset{\prime\prime}{.}}}9 $ (3 × 3 pixel) aperture, capturing ∼80% of the flux for YE and JE, 70% for HE. For NISP spectroscopy, we consider an emission line with a flux of 2 × 10−16 erg s−1 cm−2 at an observed wavelength of 1600 nm, measured in a 4  ×  4 pixel wide aperture in the dispersed images.

In this way we verify that the scientific requirements of the Euclid project are met. Global statistics of the S/N are summarised in Table 7. The median survey depths converted and scaled to a 5σ point-like source (5σ point-like source) performance metric for imaging are listed.

Table 7.

S/N statistics for the RoI, for each channel: IE band (boldface values refer to extended objects), NIR bands (JE, YE, and HE values refer to point-like objects), and red grism band, RGE (italic values in the last column refer to 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $ diameter sources).

The resulting S/N maps for VIS and NISP are shown on Fig. 19. All four quadrants are fully green, within specifications, for all channels for their respective depth metrics.

thumbnail Fig. 19.

Euclid S/N maps at detection limit in equi-rectangular celestial projection (same referential as Fig. 17). Top: VIS IE band S/N; the entire RoI is within specifications (S/N ≥ 10, the first level of light green). Middle: NISP S/N for the JE band (the YE and HE bands are similar; see Table 7 for overall statistics). Bottom: NISP RGE band S/N. The more convoluted isocontours of the IE S/N are caused by the stronger dust extinction versus the NIR. Greyed regions outside the RoI are excluded due to extreme reddening, high stellar density, and/or high zodiacal background (cf. RoI definition and Fig. 17).

We note that the S/N computations do not consider the contamination of galaxy samples by stars; to this end we introduced the thresholds to Galactic latitude. The greyed areas in Fig. 19 illustrate where a certain component (such as extinction) is out of range. These may appear inside the RoI (e.g. at the location of the SMC). Non greyed areas outside the RoI reflect an evolution of the criteria that led to the RoI definition, for example by tightening the Galactic latitude threshold from |b|≥23° to |b|≥25° after the northern ground-based surveys had been defined for |b|≥25°.

Our more complex zodiacal model (Sect. 5.1.1) predicts a lower background that varies with time and position along the orbit. This modulation happens at a level far below the typical range of zodiacal background within the RoI (Fig. 12), and hence we do not expect the median performance to change with this model.

In summary, Fig. 19 shows that the S/N in the VIS band exceeds the requirement of S/N ≥ 10 over the whole RoI, with a median value of nearly 16. This gain is mostly related to longer than required integration times, driven by the needs of the spectroscopic channel. Likewise, the YJH photometric data are well above the S/N ≥ 5 requirement. A negligible area (less than 50 deg2) of the NISP spectroscopy is below the S/N ≥ 3.5 requirement (Fig. 19 bottom). The median S/N for spectroscopy is 4.5, a comfortable margin. Hence all specifications are exceeded, and on average Euclid will go deeper than initially planned.

5.2.3. Best S/N areas of the Euclid sky

The areas of sky where the largest S/N can be achieved on average are offset from the ecliptic poles due to the out of field stray light from the Galaxy and the Large Magellanic Cloud (LMC). The S/N in these areas is close to the maximum values listed in Table 7. Figure 20 shows yellow and orange filled areas that were derived from an average of the VIS, Y, J and H S/N maps. The boundaries have been smoothed in this representation. The areas shown in yellow represent the best 1300 deg2 in each Galactic cap (or celestial hemisphere), and the orange area the best 2600 deg2 (including the yellow area). The EWS seeks to cover these best areas first.

thumbnail Fig. 20.

Best S/N areas in each Galactic cap. The best 1300 deg2 in each Galactic cap is shown in yellow, and the best 2600 deg2 in red (including the yellow area). These are similar to the areas planned for the first and second data releases, 2500 deg2 and 5000 deg2, respectively (cf. Fig. 46).

5.3. Limiting surface brightness

As described above, we have a complete knowledge of the background reaching Euclid’s focal planes. Similarly to what we did for faint compact sources, we can derive the detection performance for diffuse emission such as non-resolved stellar populations in tidal streams of galaxies, intracluster light, and the cosmic infrared background (CIB). In the following, we do this considering the same noise properties as before, assuming a stack of three exposures for both VIS and NISP, as this applies to the 90% level coverage of the imaging survey (see Sect. 4).

We use the limiting surface brightness metric adopted by Mihos et al. (2013), which is based on the asinh magnitude introduced by Lupton et al. (1999). This conservative metric is a good description of the actual signal properties at very low σ-levels relevant for such science. It has the merit of reflecting an actual science performance: the determination of the light profile of Messier 101 down to a surface brightness of μB = 29.5 mag arcsec−2 corresponds to the 1σ limit in the Mihos et al. (2013) study. At 1σ, the asinh magnitude is 0.5 mag smaller than the corresponding standard magnitude. This limiting surface brightness is computed directly from the photometric zero point of the system and the background noise property from our S/N study. Because of the very limited contamination of bright stars (see Sect. 5.1.3) here we only need to consider the diffuse background. The limiting surface brightness expressed at the pixel scale is

μ = ZP 2.5 log 10 ( b ) a asinh ( f / 2 b ) , $$ \begin{aligned} \upmu = \mathrm{ZP} - 2.5\,\mathrm{log}_{10}(b) - a\,\mathrm{asinh}(f/2b), \end{aligned} $$(6)

where, at the native resolution, ZP is the photometric zero point in electrons per second; b is the noise per pixel in the image, assumed to be the Poisson standard deviation of the background counts per pixel, B (i.e. b = B $ b=\sqrt{B} $; a = 2.5 log10(e) = 1.08574); and f is the level in electrons per pixel of the extended astronomical source.

Given Euclid’s small plate scales ( 0 . 1 $ 0{{\overset{\prime\prime}{.}}}1 $ and 0 . 3 $ 0{{\overset{\prime\prime}{.}}}3 $ pixel−1 for VIS and NISP, respectively), the depth metric relation must be brought to the physical scale of common features encountered in the near-field (galaxies, streams, shells, dwarfs, etc): we adopt a generic 10″ × 10″ scale while conforming to the standard performance unit for extended emissions in magnitude per square arcsecond. We shift from the pixel scale to our scale of interest considering a square area of n native pixels on the side. By averaging over the larger area, the estimate for the noise is scaled down by a factor n × n $ \sqrt{n\times n} $ while the zero point gets shifted by −2.5 log10(n × n) for flux conservation. Scaling to the magnitude per square arcsecond unit adds −2.5 log10(p), with p the area of the native pixel in square arcsecond (0.01 for VIS, 0.09 for NISP). For the adopted 10″ × 10″ scale (n = 100 for VIS, 33.33 for NISP) the combined effects on the limiting surface brightness amount to −2.5 log10(1/(0.01 × 100)) = 0 for VIS, and −2.5 log10(1/(0.09 × 33.33)) = +1.193 for NISP (the larger the physical scale, the lower the noise and the greater the performance).

We can now explore the Euclid RoI at the asinh 1σ level (f = b). The result is a map of limiting surface brightness for each band. The VIS and NISP maps show essentially the same structures (shape, amplitude, location), and in Fig. 21 they are combined in a single map (see the colour bar for the amplitudes in IE and JE band). The maximum range from the best area to the worst, at the ecliptic plane limit, is ∼0.4 mag arcsec−2, corresponding to the total background level ratio of 2.25 between these best and worst areas. The median limiting surface brightness across the four bands over the RoI is: IE = 29.8, YE = 28.2, JE = 28.4, HE = 28.4 AB mag arcsec−2, 1σ asinh magnitude at the 10″ × 10″ scale (−0.25 mag for the minimum performance, +0.15 mag for the maximum over the RoI). We note that our various background components are conservative estimates and these levels can be considered a safe performance. Our more complex zodiacal background model presented in Sect. 5.1.1 indicates how intensity varies with time and position along the orbit. This will in consequence modulate the total background, hence the depth, although at a level lower than the 0.4 mag range depth seen here over the RoI (Fig. 12). The median performance is not expected to change. We also note that a depth metric based on the integrated Sersic radial profile over a whole galaxy, digging deep into the noise, typically adds at least two magnitudes with respect to this contrast-oriented metric. This is taking particularly into account Euclid’s pristine image quality that will enable an effective masking of the foreground and background compact sources.

thumbnail Fig. 21.

EuclidIE and JE band 1σ asinh limiting surface brightness at the 10″ × 10″ scale (low surface brightness performance) and the nearby Universe: Local Group galaxies are represented by large cyan dots and galaxies up to z = 0.03 by red dots, while the purple ellipses indicate the actual sky area covered by the nearest clusters of galaxies (Virgo, Coma, Fornax, and Hydra are within the RoI). The RoI exhibits a 0.4 mag range in sensitivity for the four Euclid bands from the best area down to the ecliptic plane.

For an illustration of the scientific potential, the map in Fig. 21 features the nearby extra-galactic Universe up to a redshift of z = 0.03: more than 10 000 bright (K-band magnitude < 12) galaxies (2MRS catalogue; Huchra et al. 2012), including several members from the Local Group, and four nearby clusters of galaxies all falling within the Euclid RoI (some additional targets that are located outside the RoI might be observed during the unallocated time; see Sect. 8.2). The three EDFs are shown in yellow on the map of Fig. 21. They are designed to be +2 mag deeper than the EWS for compact sources, but the gain in depth is comparable for the diffuse emission (for specific details, see our companion paper on the EDFs, [Sc23]).

Such a capacity for detecting faint nebulous objects will make the diffuse galactic light (DGL; cirrus) an ubiquitous component of the background over the entire RoI as it averages to a level of 27.1 mag arcsec−2 in IE (derived from multi-band dedicated CFHT-MegaCam observations to help characterise the RoI). Based on the scaling of the DGL’s albedo (Gordon 2004), we derive an average of 27.9 mag arcsec−2 over the RoI in the NISP JE band. This extra background is shown in the map in Fig. 22. We note that it is included in our derivation of the limiting surface brightness, but it has no impact since it is more than five magnitudes fainter than the combination of zodiacal background and out-field stray light. We know that structures exist in the DGL down to the arcsecond scale even at high Galactic latitude (Miville-Deschênes et al. 2016). The DGL is, however, truly diffuse over the great majority of the RoI, although as shown in Fig. 22 some parts of the Euclid sky are inevitably worse than others. The two islands of the RoI (regions II and III; see Sect. 5.2.1) correspond to the worst parts of the Euclid sky due to their proximity to the Galactic plane and a considerable presence of cirrus.

thumbnail Fig. 22.

E(B − V) over the RoI at the native 5′ resolution of the Planck/IRAS map, and surface brightness of the DGL (cirrus) in the VIS and J bands. The two islands (smaller parts of the RoI) collect the worst areas of the Euclid sky due to the combined presence of cirrus and the proximity of the Galactic plane. The small green ellipses show the three EDFs.

5.4. Expected number of sources

The expected counts of galaxies that satisfy the WL requirements in terms of S/N > 10 (as measured by SourceExtractor, Bertin & Arnouts 1996) and size (FWHMgal > 1.25 FWHMPSF) can be inferred using realistic image simulations. In our analysis, we used sky simulations produced with SkyLens7, in which the input galaxies magnitude, spectro-morphological, and redshift distributions were drawn from the HST Ultra Deep Field (HUDF). Since the HUDF is a very small field (11 square arcmin), the galaxies therein may not be representative of the mean properties of galaxies on the whole sky. Given the larger size of the Cosmic Evolution Survey (COSMOS) field, its counts are more robust. We thus corrected the HUDF counts so as to reproduce the magnitude distribution in the COSMOS field. We also estimated the S/N based on the specified throughput of the VIS instrument and including Galactic extinction and the zodiacal background. We note that these estimates are conservative because the latest throughput estimates for the VIS S/N are larger than the initially specified ones that we used.

By simulating images under various conditions, and extracting the sources from them, the dependences of the galaxy density with extinction and background were determined, as shown in Fig. 23. The results were interpolated from each pointing to derive an estimate of the spatial map of number counts, shown in Fig. 24. These estimates are currently being updated to include stray light, which is quite dependent on the local environment of bright stars and their spectral energy distribution (SED; cf. Sect. 5.1.3), and to include the latest instrument parameters and data reduction methods.

thumbnail Fig. 23.

Expected average number density of IE galaxies that satisfy the WL requirements (we notice that the local impact of foreground galaxies or clusters is not considered here). Upper panel: counts as a function of the local background (in mag arcsec−2) for different levels of E(B − V) reddening. Lower panel: counts as a function of reddening for different levels of the background. Lines are linear fits to the simulations.

thumbnail Fig. 24.

Expected counts of WL galaxies in arcmin−2 (see text; we note that the local impact of foreground galaxies is not considered here) in ecliptic coordinates with the Galactic plane removed (black). The decrease in the number counts is evident going towards the ecliptic plane, where the zodiacal background is higher. See also Scaramella et al. (2015).

A similar approach was carried out for GC, where spectra were simulated, then extracted and finally measured (Zoubian et al. 2014), yielding a similar sky map of expected number of reliable redshifts over the sky (see Jamal et al. 2018, other methods are being developed). Figure 25 shows an updated version from the science performance verification #2 (SPV2) GC simulation (SPV is an end-to-end simulated exercise of the whole data reduction chain8). As for the WL case, the GC simulations are currently being updated to include the latest estimated effects of stray light, instrument characteristics and data reduction procedures. We expect that the overall fraction of recovered redshifts in simulations will increase with updated models and throughput, coupled with better simulations and data reduction, since the simple but up to date estimates for spectra now yield S/N > 3.5 almost everywhere (see Fig. 19). The indicative expected number of reliable redshifts of each field (used as a relative weight) as a function of time for a typical EWS (see Sect. 7) is shown in Fig. 26. Only 2% of the fields have a weight below 1000, the consequence of an excellent survey efficiency and a good RoI selection. The overall normalisation of these numbers is based on model estimates of the intrinsic number of Hα emitters plus current estimates of contamination rate and a preliminary data reduction success rate. The latter is expected to be improved. At present, the number of H-alpha emitters observable with Euclid remains uncertain, but recently Bagley et al. (2020) forecasted this number to be ∼3300 deg−2 in the redshift range measurable by Euclid. According to this forecast, to meet the initial Euclid specification of 1700 deg−2 would require the redshifts of half of the underlying population to be reliably measured.

thumbnail Fig. 25.

Expected counts of GC galaxies (with reliable spectroscopic redshifts) in deg−2 (see text) in ecliptic coordinates. Here the number counts decrease faster with latitude than for the WL counts (cf. Fig. 24) since the impact of stars (spectra) from our Galaxy on the GC one-dimensional spectra is more severe than the impact of stars (point-like) on the two-dimensional WL data.

thumbnail Fig. 26.

Weight (defined as the GC expected number of good redshifts per deg2) of observed fields along a typical EWS progression (see Sect. 7). The moving median of the weights in the previous year is superimposed (red line, starting on year 1). The approximately regular cusps every 6 months are due to the paucity of observable fields close to the Galactic and ecliptic planes and their increased background and extinction. The decrease in the median and the upper envelope in the second half of the mission is due to a relatively low quality of the sky areas still available for observations in those periods. In fact, the best areas are observed in the first half of the mission (see Figs. 20 and 46).

5.5. The EWS coverage by ground-based telescopes

Weak lensing tomography and the need to account for the contamination by galaxy intrinsic alignments require solid estimates of the redshifts of the galaxies used as sources in the weak lensing analysis. Euclid will exploit galaxies up to a redshift of z ∼ 2 with the majority of galaxies, at the lensing goal of mVIS = 24.5, lying at a redshift z ∼ 1. A typical SED of such a galaxy is shown on the simulated spectrum of Fig. 27. The NISP depth goals (YE, JE, HE = 24.0 AB mag; see Sect. 3.2.2) were scaled to capture the flux of this galaxy population at the required S/N for proper photometric redshift derivation. However, since the weak lensing imaging band through the broad VIS band IE does not sample key features of a z = 1 galaxy energy distribution, in particular the 4000 Å break that falls within the i band (Fig. 27), complementary bands are needed to reach the required redshift precision. The g, r, i, z bands are critical in particular, as introduced by Laureijs et al. (2011).

thumbnail Fig. 27.

Simulated spectrum of a typical galaxy (in terms of mass and star formation rate) from the Euclid weak lensing population at z = 1 (D. Masters, priv. comm.), overlapping the IE, YE, JE and HE bands. The Vera C. Rubin Observatory photometric bands (u, g, r, i, z, y) are also shown, highlighting the importance of the ground-based data to the photometric redshift derivation as the 4000 Å break falls within the i band. The total sensitivity for the ground is given at 1.2 airmass.

Large projects aiming at obtaining photometry through the Sloan bands over large parts of the sky were on the rise at the time of the Euclid mission definition. Since photometry in those bands does not require the observing conditions of a space observatory, this critical part of the mission was left for external up-and-coming photometric surveys by ground-based facilities located across the two celestial hemispheres in order to reach the entire Euclid sky. The minimal depths needed to derive photometric redshifts for the WL probe are 25.7, 25.1, 24.8, 24.6 AB mag in the g, r, i, z bands, respectively, for 5σ point-like source. These levels were first introduced in Laureijs et al. (2011) and later on fine tuned to optimally match the SED of the z = 1 IE = 24.5 galaxy populations anchored on the YE, JE, HE depths of 24.0 that will be achieved by Euclid.

At the time of the mission selection in 2012, Laureijs et al. (2011) commented on ground surveys that were still speculative: only the Dark Energy Survey (DES) was about to start on its broadband imaging effort, with nearly 4500 deg2 of its of 5000 deg2 goal overlapping the EWS. As of the end of 2020, the now completed DES has secured in the g, r, i, z bands coverage of nearly a third of the EWS area over the south Galactic cap while the rest of the Euclid RoI is an on-going effort. Together, the following six most powerful ground-based wide-field telescopes will eventually deliver the photometry needed by Euclid across the u, g, r, i, z bands from the northern and southern hemispheres (Fig. 28), the u band being a solid bonus for photometric redshifts at any depth.

thumbnail Fig. 28.

Current plan (2021) of the ground-based coverage of the Euclid sky based on the six most powerful ground-based wide-field imagers and the Euclid RoI (in solid yellow). Vera C. Rubin Observatory (LSST and suggested northern and southern extensions: purple, green, and pink rectangles) and Blanco (DES, red) in the south. Subaru (WISHES), JST (JEDIS-g), Pan-STARRS, and CFHT (CFIS) in the north (orange). The best S/N 2600 deg2 areas for Euclid in each Galactic cap (black dotted contours) are key to steering the ground efforts with respect to the mission’s first data release (DR1). The Galactic referential is shown in light green.

Extensive community-based actions led to the Canada-France Imaging Survey (CFIS; Ibata et al. 2017), which ought to cover by 2025 the northernmost 4800 deg2 of the Euclid RoI in the u and r band, using 314 MegaCam nights on the 3.6 m Canada-France-Hawaii Telescope (CFHT). Spain’s 2.6 m Javalambre Survey Telescope (Cenarro et al. 2018) should start in 2021 covering that area with the Javalambre-Euclid Deep Imaging Survey in the g band (JEDIS-g; 100 nights). Pan-STARRS (USA, 2 × 1.8 m telescopes, Chambers et al. 2016) joined in 2018 to provide the i band by 2025 as a result of their on-going near-Earth object (NEO) search. Finally, a group of Japanese scientists joined the Euclid Consortium in 2020 through the contribution of Subaru Hyper Suprime-Cam (HSC; Miyazaki et al. 2018) time (40 nights). WISHES (Wide Imaging with Subaru HSC of the Euclid Sky) will cover the northern area in the z band, the most demanding band in terms of depth, hence requiring an 8 m class telescope. The telescopes actively collecting data are now working in concert as part of the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS), an independent consortium motivated by the shared effort for Euclid, to cover the northernmost sky over the complete set of photometric bands, with a completion date around 2025. Canadian and University of Hawaii UNIONS members launched an effort in 2019 to gather g-band data with Subaru-HSC to complement the Spanish effort. We note that since CFHT and Subaru cannot effectively observe from Hawai’i at declinations δ ≥ +80°, the EWS RoI has been trimmed around the equatorial pole by a few tens of degrees, a minor hit since the area was already mostly rejected due to high dust extinction.

Meanwhile the Vera C. Rubin Observatory (USA, Ivezić et al. 2019) is approaching first light and the start of the Legacy Survey of Space and Time (LSST) should be in phase with Euclid. When on the sky, Rubin will be the most powerful wide-field imager ever built and the Euclid minimal depths in the ugriz bands will be reached within one year of normal Rubin LSST operations over the 8000 deg2 overlapping the Euclid RoI in the southern sky; this will supersede the DES dataset. The Rubin Observatory being such a powerful machine, the Euclid Consortium is investigating with the Rubin community a northern survey extension serving various strategic Rubin scientific niches (Rhodes et al. 2017). Such an extension (3000 deg2 of Euclid RoI area) would fill the +2 to +30 deg declination gap between the main component of the LSST and the on-going Euclid northernmost sky effort (Fig. 28).

6. Calibration and deep field observations

Euclid’s three major modes of observation (VIS, NISP-P, and NISP-S) and their tight scientific requirements imply a thorough and extensive calibration programme throughout the mission, which serves two main purposes, namely the calibration of the flight hardware (instrument calibrations; see Sect. 6.1) and the characterisation of the target galaxies and quantification of any biases that may arise in the WL and GC experiments (sample characterisation; see Sect. 6.2).

6.1. Telescope and instrument calibrations

The hardware calibration focuses on the properties of the optics, detectors, electronics, and the opto-mechanical aspects of structural components. The associated performance regarding PSF, throughput, quantum efficiency, noise, and bias is subject to change due either to variations in operational conditions (orbit, de-pointing, on-board power dissipation) or to long-term ageing (micro meteorite pitting, particles and UV radiation damage) The optical performance is furthermore affected by contamination from material outgassing. Monitoring and accurate correction of these effects is paramount to the scientific success of Euclid and requires repeated execution of calibration observations with varying cadences.

The hardware calibrations can be divided into on-sky (e.g. transmission) and off-sky calibrations (e.g. flat fields). The latter have the least constraints as they can be executed independently of the spacecraft’s pointing; the scheduling must merely respect the required cadence within some tolerance. On-sky calibrations have the additional constraint that they must minimally disturb the thermal equilibrium of the spacecraft. To this end, a selection of targets is available, from which we choose those that optimally merge with the scientific observations.

In the following we provide a summary of the main aspects of the calibrations that impact the building of the EWS (excluding the additional calibration data taken during the PV phase, and shorter instrumental calibrations that are integrated in the ROS; see Sect. 4.1).

6.1.1. Self-calibration and VIS non-linearity

Quite important for hardware calibrations are the self-calibration observations, a block of about 18 h observing a field near the NEP with perennial visibility (for details about this field see [Sc23]). Besides monitoring the total system transmission, these observations provide the data for a large range of additional calibration products. The self-calibration observations are scheduled eleven times per year, approximately on a monthly basis, and back-to-back with a VIS non-linearity sequence of about 9 h duration.

6.1.2. VIS PSF calibration

Another large calibration block is a sequence of VIS PSF observations lasting 20.5 h, targeting one of about a dozen stellar fields featuring a suitable range of magnitudes and SEDs, while minimising Galactic extinction and polarisation effects on the PSF ellipticity. The Euclid VIS PSF has in fact a weak dependence on polarisation of the incident light. Most fields in the EWS and EDS have low levels of Galactic polarisation at VIS wavelengths, but we must ensure that the PSF calibration fields are also selected to have low polarisation. Additional dedicated observations of polarised regions are separately planned to measure the PSF polarisation dependence, in orbit. These observations are used to routinely update and validate the VIS PSF model, and must be taken with the spacecraft in thermal equilibrium. This condition is met after about one week without large changes in the SAA and AA attitude angles (i.e. as close as possible to the values used in the preceding days of EWS observations). The VIS PSF calibration data hence need to be embedded in the currently executed EWS patch. In the global schedule (Fig. 29), the longitudes of the PSF calibration fields are marked on the top row, and the related observations are shown by the red strips within the EWS observation blocks.

thumbnail Fig. 29.

Example of the full schedule, including the stage-1 schedule discussed in Sect. 6.3.3. Time goes from left to right for each row, which spans one year of the survey. Subsequent years go from top to bottom. See text for a description.

6.1.3. NISP calibrations

Two types of NISP calibrations are relevant in the description of the EWS. First, the NISP nonlinearity calibration requires a data-intensive special readout mode and must be done one detector at a time, requiring a total of 49 h. In the current implementation of the EWS, these calibrations are scheduled approximately every six months. In a future version of the EWS, these observations might be partially executed in parallel with one of the VIS on-sky calibrations, pending a confirmation of the instrument inter-operability.

The second set of NISP calibrations is a one-time check of the NISP-S wavelength dispersion solution (NISP-S-PN-1 in the 2nd year of Fig. 29). This is initially obtained during the PV phase prior to the beginning of the survey, and repeated once, about one year after the start of the EWS. The dispersion solution is obtained from compact planetary nebulae (PNe) with strong emission lines, being stepped across a larger number of positions in the NISP focal plane. The dispersion solution is transferred to the self-cal field (Sect. 6.1.1) to establish a set of secondary standards to monitor the stability of the dispersion. Currently, the PN can be chosen from a list of 24 suitably compact PNe, which will be down-selected to some extent pending further ground-based spectroscopy.

6.2. Euclid deep fields and Euclid auxiliary fields

The Euclid schedule devotes much time to deep observations for galaxy sample characterisation. For GC, one needs to quantify biases in redshift measurements due to contamination and emission line misclassification (completeness-purity calibrations). For WL one needs to quantify biases in shear estimation due to noise (Viola et al. 2014), colour gradients (Semboloni et al. 2013; Er et al. 2018), and the calibration of photo zs. To this end, Euclid will observe three types of fields: (i) deep observations of six well-known fields for photometric redshift calibration and colour gradient calibration purposes (CDFS, COSMOS-Wide, SXDS, VVDS-Deep, CANDELS/AEGIS, and CANDELS/GOODS), hereafter known as the Euclid auxiliary fields (EAFs), which have extensive ground- and space-based multi-wavelength photometric and spectroscopic coverage, and are covered with 1−4 Euclid FoVs (i.e. spanning 0.5−2.0 deg2); (ii) repeated observations of two 20 deg2 fields at different times to obtain different dispersion angles to calibrate spectral confusion, named the completeness-purity-calibration fields (CPC); and (iii) deep observations of large (10−20 deg2) fields, two magnitudes deeper than the EWS, for calibration of the noise bias, that constitute the three EDFs (EDF-North, EDF-South, and EDF-Fornax).

Defining the exact locations and footprints of the three EDFs required substantial effort9 EDF-North, located at the NEP, is visited 40 times with the ROS to reach a depth two magnitudes deeper than the EWS, while the other two EDFs, at lower latitudes and larger local background, need more visits to reach the required depth. To maximise synergy, EDF-North and EDF-South are chosen to have the same centres of CPC-North and CPC-South, respectively.

The EDFs and EAFs will have great scientific and legacy value due to the depth of the observations. A detailed description of their planned observations is presented in [Sc23].

6.3. Computing the stage-1 schedule

ECTile is the software purposely developed to generate the scheduling of the EWS. It has two major stages. In the preparatory stage-1, ECTile computes the schedule of the various calibrations, EDFs and EAFs. We review stage-1 in this section. The core of ECTile is stage-2, which computes the EWS; it is described in Sect. 7.

The observations of calibration fields, EDFs, EAFs, and also ecliptic ‘polar caps’ (see below), hereafter called ‘targets’, share the same traits: most are only visible during a short interval at a given time (apart from the poles themselves), must observe a specific region (some with a specific orientation), and some must be re-observed with a fixed cadence. Given these restrictive properties, their schedule is computed in stage-1, followed by the scheduling of the EWS in stage-2 in the remaining time available. When scheduling the targets we need to make sure to leave enough time for EWS in each year, in order to fulfil the planned public delivery of reduced data to the community (Sect. 1).

We note that for the reasons explained in Sect. 7.4.2, the polar caps at high ecliptic latitudes (|β|≥79°) are also part of the stage-1 schedule. Each of the two polar caps is covered with a fixed number of patches (thirteen in the northern polar cap, six in the southern polar cap), defining targets that are scheduled with a locally fixed patch area and sequence (i.e. these two regions of the EWS are observed with the same strategy as the EDFs and EAFs).

The resulting schedule is called the ‘stage-1 schedule’, and its computation consists of three steps: (i) analysis of each target (strategy and visibility); (ii) placement of each target in longitude and year; and (iii) assignment of timestamps to each placement.

6.3.1. Target analysis

In general, a target has four types of constraints: a fixed location on the sky, a specific observation sequence, a minimum depth, and a cadence (or other time constraints). From these we compute the target’s ‘window of visibility’, namely, the range of ecliptic longitudes (of the Sun) in which the target may be observed (on the leading or on the trailing side of the orbit). Typically, each target has two windows of visibility per year.

The computation of each target window of visibility reveals its flexibility in placement and, more important, possible conflicts (of placement) with other targets. The procedure is essentially manual (aided by software tools). Three cases are of particular interest (see [Sc23] for details).

The first case is the EDF-Fornax, which is planned to be observed gradually along the mission for a total of 52 times, taking into account the larger local background. It can be observed twice per year. However, the window of visibility of the EDF-Fornax partially collides with that of the EDF-South, and with some orientations of the CPC-South. This means that scheduling one of these targets strongly restricts the placement of the others. To overcome this conflict, the visits to CPC-South are all scheduled in the first year, leaving those longitudes free for EDF-Fornax, in the following years (5 times once every six months, plus a short 2 times visit on the first year, as required). Likewise, the EDF-South is observed from the second year onwards (but offset from EDF-Fornax).

The second case is the scheduling of the COSMOS and SXDS photo-z calibration targets. These require an observation of a 2 × 2 pattern that, if observed in one go, would take 15 days each. This would pose a great difficulty for the scheduling of the EWS, by adding a long interruption. In general, when scheduling the EWS, it is possible to bridge over interruptions (such as calibrations), if shorter than five days. Long interruptions are not bridgeable, forcing a interruption of the EWS buildup, something that must be compensated in the next passage (at least a half-year later). It also reduces the opportunity to place PSF calibrations (that require a week of EWS observations prior to calibration). The solution is to observe these fields in four visits, observing half of it each time, first to an intermediate depth, and then a second time to the final depth.

The third case are the self-calibration and VIS non-linearity calibrations, which both require approximately a monthly cadence of observations. In order to decrease the number of interruptions to the EWS observations and the number of large slews used, these observations are always scheduled in sequence and executed at the same sky field, defining a ‘recurring calibration block’.

Finally, every four weeks starting on Mondays at noon (UTC) ±1.2 h, there is one block of 12 h that is reserved for spacecraft orbit maintenance operations (SOP). During SOP time survey data cannot be taken.

6.3.2. Target placement

The second part of the computation of the stage-1 schedule consists on placing each target in a table (named the ‘design schedule’), at a given longitude and year, striving to avoid collisions with other targets. In that table, choices have to be made such as what target goes into each year, and at what longitude to place it. There is some freedom in this procedure. For instance, many targets do not need to be observed in a specific year. However, some targets are more stringent than others (i.e. they have shorter visibility windows or have a single annual visibility). Therefore, for simplicity, it is preferable to place the targets by decreasing order of perceived difficulty, filling the table year by year. This process is performed manually.

While being manual, the process does not need to be very precise: some overlap in longitude between targets is allowed. Any overlaps are rectified in the next step and the longitudes are converted to timestamps, transforming the design schedule into the stage-1 schedule.

The design schedule is filled according to a strategy that prioritises the placement of the targets with the most constraining observing conditions. In the first year the order of placement is as follows: (i) place all ten CPC-South visits at transit longitude (we note these are the most demanding targets); (ii) place recurring calibration blocks (self-cal+VIS non-linearity) at an approximate step of 360 ° / 11 = 32 . ° 73 $ 360\circ/11 = 32{{\overset{\circ}{.}}}73 $, which promotes a synergy with CPC-North; (iii) place the ten CPC-North visits (at orientations differing by steps of 32 . ° 73 $ 32{{\overset{\circ}{.}}}73 $); (iv) place eleven of the thirteen targets of the northern polar cap, so that the width of these patches and the angular offset matches the placement of the CPC-North targets (i.e. they are observed immediately before or after), we note that the remaining two targets are placed in the second year; (v) place all eight targets of the southern polar cap; (vi) place the first two targets of the COSMOS-wide, covering half of the area needed on this target in the first year; (vii) place a ×2 pass visit to EDF-Fornax, slightly offset from its ideal placement, to avoid collision with one of the CPC-South targets; (viii) place the AEGIS target (required for the first year); (ix) place the two NISP non-linearity targets (we note that because these are freely placeable, the choice is to place them on the slot of longitudes where there is less area of EWS within the RoI); and (x) place a double visit to EDF-North to avoid a small gap of time between other targets (we note that since it can be placed all year round, the visits to EDF-North may be used as a filler).

At this point, the first year is filled. It would be possible to pack more targets into it, but that would reduce the size and number of EWS windows, which in turn reduces the opportunity to place VIS PSF calibration fields. The latter need to be scheduled within EWS windows larger than seven days (see Sect. 6.1.2).

With all CPC observations placed on the design schedule, the most demanding target left to be scheduled is EDF-Fornax. Unlike the EDF-South, the EDF-Fornax cannot be reached far from transit through a large de-pointing because of its low ecliptic latitude. The remaining years follow a common strategy, with most of the targets (EDF-Fornax, EDF-South, recurring calibration blocks, EDF-North, NISP wavelength dispersion) being scheduled almost exactly the same way. The exception are the EAF targets that vary from year to year: In the current schedule, the SXDS is placed in the third and fourth years, and the VVDS and GOODS-North in the fifth year.

The strategy for the remaining years is as follows: (i) place two ×5 passes to EDF-Fornax at their transits; (ii) place two pairs of visits to EDF-South, clustered around the two EDF-Fornax pointings, consisting of a first pair of a one pass visit and a ×2 pass visit, separated by 5°, and a second pair of two ×2 pass visits, also separated by 5° (so to ensure survey windows not smaller than 5 days, given an average orbital progression of 1° per day). Each group is clustered around each EDF-Fornax pointing; (iii) place the 11 recurring calibration blocks, slightly offset from a nominal cadence to avoid conflicts with the EDF-Fornax and EDF-South observations; (iv) place the two remaining northern polar targets (second year only); (v) place the large block of NISP wavelength dispersion target (second year only); (vi) place EAF targets (these vary from year to year); (vii) place two visits to EDF-North, one ×4 pass and another ×2 pass, near the centre of the range of EWS RoI scarcity, matching surrounding recurring calibration blocks; and (viii) place the two NISP non-linearity targets.

6.3.3. Timestamp assignment

The last stage of the computation of the intermediate schedule is the automated conversion of the year-by-year longitudes into timestamps. Given a starting date to the survey routine phase (currently but not frozen yet, this is expected to start on 8392.5 modified Julian date, or 2022-12-23T11:59:23Z), it is easy to compute the corresponding Sun longitude ( 271 . ° 30 $ 271{{\overset{\circ}{.}}}30 $). Then, it is a simple matter of ‘reading’ the design schedule, year-wise.

The first target following the starting longitude is a recurring calibration block at 290°, approximately two days after the beginning. This timestamp is assigned to the first target found. Practically, this involves only computing the next timestamp when the Sun is at a given longitude. The process continues, assigning timestamps to targets of the first year, by order of longitude (wrapping around at longitude 360°). Once the traversing reaches the initial longitude, the process continues in the second year, starting at the same longitude where the first year ended. Once the second year is completed, the process continues on the third year and so on, until the intermediate schedule is completed. There is a possibility that the initial longitude for traversing any year coincides with the middle of a target. In those cases, the traversing continues after that target.

Besides converting longitudes to timestamps, an algorithm disentangles overlaps between targets and prevents, when possible, the occurrence of too small EWS windows. The overlap of two or more targets (allowed in the previous stage), is fixed by offsetting those targets from their initial placement, minimising the overall offset, within the range allowed by their window of visibility. In most cases, this process is sufficient to resolve overlaps. In case of failure, the solution is to go back to the previous stage and fix the overlap manually. The same process is applied to eliminate the occurrence of small windows. Except that, now, the offsetting is in the opposite direction, pushing targets closer to each other.

6.3.4. Results

Figure 29 shows an example of the full Euclid schedule. The result of the stage-1 procedure, in this diagram, is the sequence of coloured boxes. The time allocated for each observation of a target is represented by a labelled box of a unique colour. The pink and white boxes represent the periods available to observe the EWS (at the end of the stage-1 scheduling they are all still unallocated). During the stage-2 scheduling (Sect. 7) the EWS is scheduled in part of the available time (shown by the pink boxes), while some time periods remain unallocated (shown by the white boxes that are increasingly longer towards the final years of the survey, cf. Sect. 8.2).

The alignment of boxes with the same colour across years indicates that the respective targets are scheduled at the same time every year. The labels (and widths) of the boxes indicate when consecutive visits are made to the same field. The EDFs are often scheduled with more than one visit in a row, in particular EDF-Fornax is usually visited five times to efficiently use its short time visibility that occurs twice per year.

Differently from the other calibrations, PSF calibrations are scheduled within the pink boxes: the EWS observation is interrupted to point to a PSF field (for ∼20.5 h) and then returns to the same position. This allows for a much better stability of the SAA and AA values used (that will match the ones used on that observation of the EWS), than if they were scheduled in the recurring calibration block, always together with the other targets of approximately monthly cadence. The visibility windows of twenty preliminary PSF fields are shown as horizontal red bars within the pink EWS boxes, while the vertical red bars show the actual scheduled time of the PSF observations (the corresponding Sun longitude for each of the fields is indicated by the numbers 1 to 20 on the axes at the edges of the Fig. 29).

Table 8 summarises the time allocated to make all calibration, EDF and EAF observations, which is 427 days. We note that the time allocated for the self-cal field is here included in the EAFs budget and not in the instrument calibrations, since its repeated observations will make it the deepest of the Euclid fields. We also note that due to the synergy between CPC and EDFs, 36 of the 48 days needed for CPC calibrations also contribute for the completion of the EDS. Details are given in [Sc23].

Table 8.

Time budget for calibration, EDF, and EAF observations.

The breakdown of the time-allocation by observing type is: (i) instrument calibrations, 25%; (ii) auxiliary fields, 21%; and (iii) deep fields, 54%.

7. Computing the Euclid Wide Survey

The computation of the EWS is a complex optimisation problem, for which we developed the scheduling tool ECTile. After the preliminary stage-1 described in Sect. 6.3 that determines the schedule of the various calibrations, EDFs, EAFs, and polar caps targets, ECTile proceeds with the computation of the EWS in its stage-2. Before turning to a detailed description of stage-2 of ECTile from Sect. 7.2 onwards, we give a brief description of pre-ECTile explorations.

7.1. Early explorations

The derivation of the optimal survey is a complex process, and ECTile is one of many possible solutions. It is, however, worth stressing that it is the outcome of a lengthy process, in which alternatives have been explored, but ultimately rejected. For instance, the first solutions of the EWS were delivered by industry to demonstrate the feasibility of the survey, but ignored some important additional considerations, such as observing areas with low zodiacal background first (see Sect. 5.1.1). Other solutions (Amiaux et al. 2012) were produced using ESA’s Euclid Sky Survey Planning Tool (ESSPT; Gómez-Alvarez et al. 2018), which allows the user to manually place patches on a sky map and fill them with Euclid FoVs. In this section we provide a brief overview of the prior investigations that led to ECTile as it is today.

To cover a maximum area with minimum overlap between single observations, a pre-determined tiling is almost unavoidable. The early approaches (see Tereno et al. 2015) therefore defined rigid FoVs placed parallel to the ecliptic meridians. The fields were scheduled starting at high latitudes and moving up and down along ecliptic meridians, observing in transit, that is, with SAA = 90° (and AA = 0°) when the local meridian coincides with the YSC direction defined in Sect. 2.1 (see also Sect. 2.5). The extent to which we move across latitudes before moving to the next longitude, defines a latitude band, to be observed in one year. Due to the convergence towards the poles, the number of fields per band decreases with latitude, and conversely the height of the bands increase with latitude.

Whenever there is an interruption in the EWS schedule due to an observation of calibration or EDF targets, a corresponding gap is left in the band. After one revolution in the orbit, a latitude band is finished and the scheduling of the next band, on a lower latitude, starts. The gaps can be recovered in the following year when in transit again. For this, the height of the next latitude band needs to be smaller, in order to create a time buffer that allows one to cover the gaps while not creating new gaps in the lower band. This way, the missing area can be observed slightly off-transit, tilting the telescope arouns XSC. The tilt must be accompanied by a rotation around the ZSC axis (hence changing AA) to compensate and keep the field aligned with the tiling. However, given the very stringent constraints on AA (|AA|< 1° at the time of the early explorations) the time buffers are necessarily small, and it is only possible to partially cover the gaps. This process is illustrated in Fig. 30. It is clear that gaps in the survey are unavoidable, and attempts to fill them only leads to more gaps elsewhere. It was thus concluded that the use of a fixed tesselation and fixed latitude bands was not viable.

thumbnail Fig. 30.

Example EWS scheduling with a rigid tiling and requiring |AA|< 1°, covering parts of the northern ecliptic hemisphere. Different colours indicate different observing epochs. Left panel: empty dark regions caused by interruptions of the EWS. Right panel: empty regions only partially recovered at a later stage, and new gaps formed at lower latitudes.

The AA limit was eventually relaxed to 3° at the preliminary design review (PDR) and later to 6° at the critical design review (CDR), being currently fixed at 5°. This makes the use of a fixed tessellation a viable approach, enabling an efficient coverage of the sky. However, the use of fixed latitude bands remains not viable and more complex strategies needed to be developed, as described in the next sections.

7.2. Introducing ECTile’s stage-2

The computation of the EWS revolves around the concept of a patch, which may be loosely defined as a compact set of fields with a valid observation sequence. This is the basic building block for this computation. This section describes the steps that precede the computation of a patch, summarised in Fig. 31. The first step is the computation of the tessellated RoI, a set of tiles covering the designed RoI. In parallel, survey windows are computed from the stage-1 schedule, defining the intervals of time left to EWS observations. These two elements are then combined into patch sources, compact sub-sets of each quadrant of the tessellated RoI within reach of a given survey window. A patch source is, in general, further divided into one or more patch segments, which are sets of tiles guaranteed to be not only compact but also having their contour shaped like a lat-long rectangle (an essential property for the computation that follows it). Patch segments from the same patch source are then transformed into a patch by establishing upon them an ordered sequence of observations.

thumbnail Fig. 31.

Steps and concepts required for the computation of the EWS, linking the RoI and stage-1 schedule to patches (defined later).

7.3. Major constraints for the EWS

In this section we present the main inputs and associated constraints for the EWS optimisation algorithm, and briefly review their impact on the current solution.

7.3.1. Constraints due to overlap

We define a tile as the largest rectangle in latitude–longitude that is completely contained in a single FoV (see Fig. 32). The survey area must then be observed through geometrically contiguous tiles, with an overlap of boundaries between adjacent tiles of a 0.5% wide strip (1% overlap, overall), to cope with the non-null pointing error. The goal of this requirement is to enable efficient coverage of the sky, whilst ensuring a minimum overlap between adjacent fields. This can be achieved with a tessellation of tiles laid out along parallels of latitude, with adjacent tiles on the same row and tiles between adjacent rows overlapping by 1%. Due to the convergence towards the poles, the number of tiles per row decreases with latitude at the same time that the overlap between adjacent fields increases, as shown in Fig. 32. This source of overlap is inevitable and a consequence of tessellating a sphere with a rectangle (i.e. a FoV). Nevertheless, this overlap is only significant at latitudes closer to the poles (above ±79°). A second source of overlap comes from the placement of tiles per row. Each row requires an integer number of tiles. However, in general, this integer number of tiles covers slightly more than the 360° of a full circle. This excess contributes to an increase in the overlap between adjacent tiles (shared evenly along the row) but, being no more than one tile per row, its amortised cost is relatively small.

thumbnail Fig. 32.

Example of a few rows of a tile-driven tessellation drawn in a cylindrical projection. Individual tiles are indicated by the red rectangles, whereas the (distorted) footprints of the FoV for non-zero latitudes are shown in blue (sizes are exaggerated for clarity). The number of tiles per row varies with latitude.

7.3.2. Constraints due to SAA and AA

As described in Sect. 2.5, the limited range of the pointing angles implies that most of the sky must be observed at, or at least close to transit. As we motivate later, most of the EWS is observed with fields aligned with the ecliptic meridians. In general, these fields are almost never observed at transit, thus requiring a rotation around ZSC to realign the FoV with the local meridian. However, the constraints on AA and SAA severely limit the extent to which a field may be observed away from transit. This is highlighted in Fig. 33, which shows the region of the sky that is observable at a given transit. It is mostly constrained by the AA range, except at lower latitudes, where the constraint on SAA dominates.

thumbnail Fig. 33.

Subset of the tessellation observable at a given epoch, with the Sun pointing at 0 deg and the transit meridians at 90 deg (trailing side) and 270 deg (leading side). This shows that the need to observe fields aligned with meridians strongly restricts visibility (see Fig. 5).

7.3.3. Constraints due to slews

The limitations associated with the cost of a slew (see Sect. 2.4), imply that EWS fields observed consecutively in time must also be, as much as possible, spatially adjacent to each other. In this way, large slews are mostly reserved for moving between EWS fields and deep or calibration fields or between patches of the EWS.

Figure 34 shows the reach from a given field, when slewing with a small-slew within a limit of 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $ and 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $, for a field placed at two different latitudes. The example demonstrates that, at low latitude, it is possible to slew to all adjacent tiles; if considering slews up to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $, it is possible to slew to tiles two rows away (we recall that the slew can be depicted as the arc separating two different directions on the sky plus a rotation around the latter). However, at high latitude, the slew between adjacent tiles is much more limited by the size of the not small change in AA needed to keep the alignment with the local meridian of the tessellation. There, adjacent tiles on the same row are farther apart, strongly limiting the field-to-field slewing, which, in this case, is performed by a rotation around the ZSC-axis. For example, given the FoV width of 0 . ° 701 $ 0{{\overset{\circ}{.}}}701 $, tiles placed at latitude of 79 . ° 08 $ 79{{\overset{\circ}{.}}}08 $ are separated exactly by 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $ of longitude. Thus, above 79°, it is not possible to slew sideways between two adjacent fields (aligned with the tessellation).

thumbnail Fig. 34.

Example of adjacent FoV slew ability at high latitudes (left panel) and low latitudes (right panel). The plots show the slew from a central tile (blue) to its neighbours: up to 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $ in light red, and from 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $ to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $ in dark red. At high latitudes, slew movements between FoVs of a tessellation are severely restricted.

7.4. From EWS tessellation to patches

7.4.1. The tessellated RoI

The generation of the tessellated RoI begins with the computation of a global tessellation on the sphere, covering the sphere without polar caps from latitudes ±79° towards the equator, with non-overlapping tiles aligned with the meridians. This tessellation is then filtered by selecting tiles that have at least one of their corners inside the RoI. This represents the tessellated RoI and specifies the FoVs eligible to be observed (see Fig. 35). Figure 36 shows two zoomed views, highlighting the dependence of overlap with latitude. The EWS solution schedules a large subset of the tessellated RoI, which then becomes the Euclid ‘footprint’.

thumbnail Fig. 35.

Tessellated RoI (in green), covering the four quadrants. The dashed red line represents the RoI.

thumbnail Fig. 36.

Zoomed-in view of the tessellation for the EWS at low latitudes (left panel) and high latitudes (right panel). The increased overlap of neighbouring tiles at high latitudes is evident. Orange tiles contain a blinding star (black dots) and are excluded from the observation.

7.4.2. Covering the polar caps

As discussed in Sect. 7.3.3, above 79° it is not possible to slew sideways, when traversing fields aligned with the tessellation. Indeed, at high latitudes, the centres of adjacent fields on the same row of the tessellation are separated by an eigenslew larger than 3 . ° 7 $ 3{{\overset{\circ}{.}}}7 $. As the amplitude of an eigenslew includes the rotations needed to align the fields and not only the angular separation between the fields centres, the part of the EWS above +79° and below −79°, the polar caps, are not part of the global tessellation.

Instead, each polar cap is covered by a fixed number of patches of fixed area, as shown in Fig. 37. The northern cap is covered with 13 patches, while there are six patches covering the southern cap. The southern polar cap is smaller due to the presence of the LMC, which lies outside the EWS RoI. Each patch of the polar caps is a target field scheduled during stage-1 with the strategy used for observing the EDFs and calibration fields.

thumbnail Fig. 37.

Coverage of the northern and southern ecliptic polar caps (left and right panels, respectively), with the design of polar caps drawn in red and patches drawn in blue. Extra patches below 79° are also added to polar caps (two in the north and one in the south) to observe corners of the RoI that would be difficult to observe with the regular wide strategy. The presence of the Galactic plane affects both caps. In addition, the southern boundary avoids the LMC. Blinding stars are also shown; in particular, the small gap in the southern cap is due to the presence of the extremely bright star R-Doradus.

At high latitudes, overlap between neighbouring fields cannot be avoided, because there is longer room for a rotation of the FoV. The average FoV overlap on the polar caps is ∼18%, while the average overlap on the regular wide is under 3%. This leads to a small loss in survey efficiency.

7.4.3. Survey windows and patch sources

The next step in the computation of the EWS is to match the global tessellation with the survey windows defined by the stage-1 schedule. A ‘survey window’ is the span of time between consecutive calibration blocks. It defines uninterrupted time intervals available for observing the RoI. Currently, there are approximately 100 survey windows. As explained later, these are processed in chronological order, one at a time. But we first show how a single survey window intersects with the RoI and how it is populated with EWS observations.

A survey window begins at the end of the last pointing of a calibration block, and lasts until the first pointing of the following calibration block. Let a and b be the two pointings delimiting a survey window, and Δt a nominal observation time, defined as ROS time, including the typical field-to-field slew time. The approximate number of FoVs possible to observe within a survey window (i.e. its capacity) is then given by n = (b − a)/Δt.

During the slot of time defined by a survey window, spanning from a to b, two transit meridians scan two opposite sectors of the sphere, representing the areas within reach of the survey window. In turn, both of these sectors intersect with two or more of the four quadrants of the RoI, identifying the eligible FoVs within reach of a given survey window. The intersection of each of these sectors with a single quadrant of the RoI defines a ‘patch source’ (i.e. a contiguous subset of the tessellated RoI within reach of a survey window, as exemplified in Fig. 38). The number of patch sources per survey window varies from two, intersecting only the mainlands, to four, intersecting all quadrants. In rare configurations, a survey window intersects the same mainland twice, defining two separate patch sources.

thumbnail Fig. 38.

Example of four patch sources (in red) generated by a survey window. We notice how the patch source in the northern island intersects that quadrant only partially.

A patch source provides a simplification of the wide survey time window (‘survey window’). Given the need to observe fields in sequence, ones near each other (within the slew constraints), the observations within a wide survey window must necessarily form a contiguous compact set. As described below, these compact sets of tiles defines a patch.

A ‘patch’ is the building block of the survey. It constitutes a unit of observation of the wide, using some of the time budget for the wide (which may all or part of a survey window) and covering some of the RoI. By construction, different patches do not intersect, neither in time nor space. Moreover, the RoI is covered as much as possible in an orderly fashion, with patches stacking one on top of each other. In that respect, the EWS may be seen as a long sequence of patches (scheduled around the calibration blocks).

Patch sources of the same survey window compete for the same time, with each one (potentially) generating a patch. Hence, in general, several candidate patches are available, with the actual choice of one among the possible ones determined by several aspects. In the simplest scenario, one or more patches fill all the time available and, in that case, it suffices to choose one of them (and discard the rest). In more complex scenarios, a survey window only intersects a quadrant partially, producing a patch source that does not consume all the available time. In these cases, the solution is to fill the survey window with patches from several patch sources. One patch is selected, reducing the extent of the survey window. Then, the process is restarted, recomputing the patch sources and generating a new set of patches. The process is iterated until all available time is exhausted.

Inevitably, due to the cyclic nature of the scanning of the sky, the selection of a patch reduces the RoI available for the generation of later patches. This does not cause any problem and it is easily coped with by flagging observed tiles as they are scheduled, thus avoiding selecting them again in subsequent compilations of patch sources. However, cyclic placement of patches around the sphere creates a ‘dented’ boundary of observed tiles. Over time, this leads some of the observed regions in the RoI to acquire a boundary shaped like a polyline in latitude and longitude (see Fig. 39).

thumbnail Fig. 39.

Example of a patch source reduced by previously scheduled patches (plot of region IV only, no star skipping for clarity). Previously scheduled patches are shown in blue, with the red region depicting the patch source of some windows. We notice how the bottom side of this quadrant acquired a lat-long shape after the stacking of some patches.

7.4.4. Patch-source partitioning

In the previous section, the process of extracting a patch from a patch source was simplified, for the sake of clarity. Actually, the process is slightly more complex, requiring the definition of the concept of a patch segment. First, tiles are selected on the condition that tiles observed in the same time slot should also be close to each other. However, because the contour line of the unscheduled part of the RoI may become irregular, the condition on proximity might not be feasible to meet. The solution is to partition the tiles of a patch source into patch segments, where a patch segment is simply a group of tiles amenable to be visited with a sequence of small-slews.

When creating a patch segment, tiles must be selected evenly across longitude, matching the rate of fields observable per unit of time (approximately, 20 fields per day ≡ per degree of longitude). The first time a patch is extracted from a patch source, the base of the corresponding RoI is bounded by a straight line (the side of the mainlands close to the poles), resulting in a single, possibly large, patch segment of tiles. But, after a few iterations, the base becomes a polyline (see Fig. 39). Then, the process of selecting tiles evenly across longitude may, potentially, produce a fragmented selection, made of two or more separate patch segments (Fig. 40 shows an example of this). This fragmentation is inevitable, but it is not an obstacle. In general, a patch source produces several patch segments in order to fill its span of time. In the case of multiple patch segments, the solution is to schedule them separately, followed by a merge into a continuous single schedule.

thumbnail Fig. 40.

Example of sibling patch segments from a fragmented patch source (of Fig. 39). New patch segments are displayed in shades of light blue, while previously observed patches are in dark blue. The discontinuities in longitude of the previously observed patches induce the fragmentation of this patch source. Tiles containing a blinding star are excluded from the observation, here displayed in orange. Tiles in the vicinity of an extremely bright star are also excluded (even if not containing it), as is the case of the star marked with a circle at ∼(323, −36).

7.5. Scheduling patches

This section describes the core scheduling functionality for the computation of the EWS. We describe two algorithms that we used. Each algorithm takes a patch segment as input and produces an ordering of the tiles (i.e. a sequence) that is geometrically appropriate for observation (within the constraints of the spacecraft). We then describe how several patch segments, properly ordered, are combined in a single schedule.

We describe first the ‘look-ahead’ algorithm and later on the ‘diffusion’ algorithm that replaced the former.

7.5.1. The look-ahead algorithm

The look-ahead algorithm was the first successful attempt at scheduling a patch segment of tiles, within mission constraints, allowing the generation of a compliant survey. While superseded by the diffusion algorithm (Sect. 7.5.2), it gave much insight into the key factors at play, and paved the way for the design of the latest algorithm.

The look-ahead algorithm was designed around the idea that a scheduling sequence must traverse a patch segment of tiles following a ‘natural’ zigzag scheme, monotonously across ecliptic longitudes. At its core, the algorithm traverses the patch segment from right to left (longitude), going up and down (latitude), with minimal reversing of direction. The natural ordering is computed iteratively, moving from a given position and given direction of traversal (going-up or going-down) to the next. If the current direction is going-up, the next tile in the sequence is the first unvisited tile of the row above found by scanning the patch segment from right to left. Should that tile not exist, then the tile on the same row immediately to the left is selected. In the latter case, the direction is reversed from going-up to going-down, setting a flag that the top border of the patch segment was reached. If the current direction is going-down, the choices are reversed; the next tile in the sequence is the first unvisited tile of the row below found by scanning the patch segment from right to left or, if this does not exists, it is the tile on the same row immediately to the left. Likewise, in the later case, the direction is reversed from going-down to going-up, setting a flag that the bottom border was reached.

With this algorithm one can then define the full process. At the beginning, a starting tile is chosen (i.e. a tile on the rightmost side of the patch segment), as well as an initial direction, and a starting timestamp. Typically, there are many tiles close to or at the same longitude as the rightmost tile; all are suitable as starting tiles. Afterwards, the algorithm computes a path to traverse the patch segment, propagating along timestamps for the observations of each tile.

In general, this algorithm does not cover the patch segment completely. It may get stuck in one of three ways: (i) a dead-end is reached, with no unvisited tiles to jump to; (ii) an unvisited adjacent tile is identified, but it cannot be observed within the slew constraints; and (iii) a non-valid observation is encountered because an observation went outside the valid range of AA or SAA when assigning timestamps.

This can be understood as follows. Due to convergence towards the ecliptic poles, patch segments at high latitudes are very asymmetric, having most of their tiles either at the top or the bottom. In this case, a sequence of observations should spend more time at latitudes with many tiles, making only occasional excursions to less populated latitudes. However, the algorithm is designed for full vertical excursions (whenever possible). Also, some patch segments have unique geometric features; or corners odd enough to trap the single path of the traversing strategy.

The solution is to extend the algorithm in two ways: making it explore more paths and allowing changes of direction in mid-excursion. The first should promote sequences that adapt to odd shapes, while the second should be able to cope with the asymmetry of patch segments at high latitudes. To these ends a ‘probing step’ is implemented. Instead of blindly progressing up and down, inverting direction only at the boundaries, the algorithm first runs two probes to determine if, in the next step, the sequence should advance by going-up or going-down (i.e. whether it should continue or invert direction).

The probing step is simple. Before advancing, the algorithm first computes the natural sequence for the remaining unvisited tiles for two scenarios: first, for a natural sequence that continues going-up, and then for a natural sequence that continues going-down. The direction to take for the next step is then given by the length of the two natural sequences just computed. If the lengths are different, it takes the direction of the longest sequence (possibly inverting direction), otherwise, it just keeps going in the same direction as before. The resulting sequence is obtained by applying this probing process iteratively, until all tiles are visited (returning success) or until the algorithm gets stuck (returning failure). As an optimisation, if one of the probing sequences traverses all of the remaining tiles, that sequence is taken and completes the schedule. If the above run fails, the process is repeated by trying other starting tiles (from the subset of tiles with longitude close to the rightmost tile). Varying the starting point explores different configurations, greatly improving the chances of success. Figure 41 illustrates this process.

thumbnail Fig. 41.

Example output of the look-ahead algorithm. The sequence starts at mid-height, going downwards first. It traverses all fields running up and down, from right to left. We notice how it sometimes shows an irregular behaviour, zigzagging to the sides instead of smoothly following the expected up and down directions.

This algorithm was used to generate the EWS from the period of time between the mission PDR to the CDR. However, following the CDR it was realised that the resulting surveys exhibited uncontrolled excursions over the full AA range (visible in Fig. 41 and the left panel of Fig. 43), degrading the thermal stability of the telescope and thus affecting the PSF estimation. This was particularly acute for patch segments at high latitudes. Part of this failure is due to the restrictiveness of small slews, limited (at the time this was developed) to a maximum of 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $. The slew range was then relaxed, allowing now for a small number of slews up to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $. However, the look-ahead algorithm is intrinsically limited by its simplicity and lack of flexibility (low number of parameters), making it difficult to accommodate new constraints such as skipping tiles with bright stars. These disadvantages prompted the development of a more capable algorithm, as explained in the following section.

7.5.2. The diffusion algorithm

Following analysis of the full satellite structural thermal optical performance (STOP) made by Euclid’s industrial prime contractor, we analysed the impact of the spacecraft attitude on the PSF stability. It was found that that, in addition to the applicable limitations on SAA and AA, a further minimisation of the field to field variation in these angles was desirable.

The new diffusion algorithm, which is presented here, achieves this goal, whilst also facilitating the avoidance of bright stars. To minimise the angle variations, the patch segment must be traversed in columns of alternating directions. This strategy avoids lateral slews on the same row that would cause spikes in δAA and δSAA, detrimental to PSF stability. Moreover, since the slew reach is longer in the vertical direction (see Sect. 2.5), moving mostly in the vertical direction increases the chances of successfully skipping across rows. This feature allows the algorithm to skip adjacent tiles (if already observed), and jump over holes created by bright stars.

To illustrate the diffusion algorithm, we consider an input patch segment in the northern hemisphere close to the polar cap. The proximity to the pole highlights the patch segment convergence, which is relevant in this context. The same strategy is easily adapted to other latitudes (with less convergence) and the southern hemisphere (by swapping up and down). The algorithm is divided in two steps. The first step computes parts of the final sequence, called ‘threads’, which connects tiles along columns. The second step obtains the final sequence by tying adjacent threads together.

The computation of the threads begins by selecting the widest row of the patch segment (not necessarily the one at the bottom), assigning a thread to each of its tiles. Then, the threads are computed in parallel, sequentially joining tiles of the current row with tiles of the row immediately above, until the top row is reached. The algorithm then returns to the starting point (the widest row), extending the threads downwards and thus completing the threads.

The double step approach is needed because, in general, the widest row is not at the bottom or at the top. Typically, a patch is bounded by a lat-long rectangle, like the examples of Figs. 41 and 42. However, patches alongside the RoI boundary may get asymmetrical shapes, acquiring some of the shape of the adjoining boundary. In those cases, the widest row is usually some row in the middle.

thumbnail Fig. 42.

Example run of the diffusion algorithm. The green and blue squares represent the fields in a patch. The blue ones mark the base row. The red line segments depict the progress of the algorithm as it traverses all the fields. Panel a: diffusion after fours iterations, extending the threads upwards. The remaining panels illustrate several points of the diffusion process: in panel b two threads skip across the hole left by a bright star, while two other threads merge, in panel c two threads merge while the remaining ones reach the top, in panel d the threads are extended two rows downwards, and panel e shows the final sequence, after the linking of adjacent threads. The plots are in elliptic coordinates.

The computation of threads is based on two parameters of the rows: row length and row capacity. The length of a row i, len(i), is the number of tiles of that row. The capacity of a row i, cap(i), is the maximum length of a row or of any row above it, defining the minimum number of threads that must cross a given row. If cap(i) = len(i), then all threads visit all tiles in row i. If cap(i) > len(i), then cap(i)−len(i) threads do not touch row i, and must skip it, but participate in some row above it. Due to the convergence towards the poles, the patch segments in the northern hemisphere funnel on the upper part. Similarly, the rows capacity reduces with increase in latitude.

Let k be the row (or one of the rows) of largest length, being also the starting row for computing the threads. Each thread is initialised with exactly one tile, of row k. Let nthread be the initial number of threads. The process then begins by extending the threads, now at row k, to the tiles of row k + 1, guided by the following criteria that depend on len(k + 1) and cap(k + 1): (i) if len(k + 1) = nthread, then we have the simplest case with a one to one correspondence. It suffices to extend each thread to row k  +  1 by making an ordered assignment, from right to left (ensuring the threads do not cross); (ii) if cap(k + 1) < nthread  ∧  cap(k + 1) > nthread − 2, then it is not possible to extend all threads, since the number of tiles of the row above is smaller than the number of threads. In this case, we extend the threads closely aligned to some tile of the row above it. This is performed by trying all combinations of ordered (non-crossing) assignments between a subset of len(k + 1) threads and all of the len(k + 1) tiles of row k + 1. For each trial combination that was accepted, we compute the cumulative variation in longitude (‘verticality’) of the thread. The combination with the lowest variation in longitude is selected, extending len(k + 1) threads (and leaving the rest unchanged); (iii) if cap(k + 2) < nthread − 2, it is also not possible to extend all the threads but, unlike the previous case, it is possible to eliminate a pair of threads (leaving nthread − 2 threads). Actually, it is necessary to eliminate threads to avoid the risk of having threads stalled (i.e. not reaching the top of the patch). A pair of threads may be eliminated if they are adjacent in the list of threads. To eliminate them, it suffices to merge the two threads together, short-circuiting their paths, by connecting the two adjacent threads to the same tile (of row k + 1, creating an inverted ‘v’). This case is slightly more complex than the previous one. Now, it needs to try all combinations of ordered (non-crossing) assignments between a subset of len(k  +  1)+1 threads and all of the len(k + 1) tiles of row k + 1, considering that pairs of consecutive threads extend to the same tile (merging those threads). Like before, all combinations are scored against verticality, selecting the one with lowest cumulative variation in longitude. It extends len(k  +  1)−2 threads, leaving nthread − 2 active threads, and merges two threads (ending their progress).

This process is illustrated in Fig. 42. Figure 42a shows the first step, extending threads from the initial row (in blue) to the row immediately above. Given the difference in row length, one thread is held up. After a few iterations, Fig. 42b shows a case where a thread on the left-side jumps over four threads (approximately at longitude 210°, latitude 60 − 68°), and a thread in the middle (approximately at longitude 209°, latitude 66 − 68°) is kept straightly vertical by skipping a row. This is an example of the capability of the diffusion process to adjust to a varying row length, by ‘squeezing’ more threads than the length of each individual row, promoting straighter threads. After a few more iterations, Fig. 42c shows two mergers of two neighbouring threads (rows at latitudes 70° and 71°). This reduces the number of threads from seven to five (which is enough to cover the tiles of the rows above). This mechanism ensures a monotonous decrease in the number of threads, keeping it close to the row’s capacity. Figure 42d shows the threads fully extended upwards (with two pair of threads ended being merged in the process). In this case, there is no need to extend the threads downwards also, since the starting row is also the bottom row. Lastly, Fig. 42e shows the final sequence, obtained after connecting adjacent threads (arrows show the temporal sequence of the covering).

The net effect of the above strategy is to grow threads upwards, striving to be as vertical as possible, and merging adjacent threads when necessary, to cope with the reduction of rows capacity (with latitude). In a way, this growth resembles a diffusion process, hence its name. It should be noticed that the number of threads crossing each row may be larger than the actual number of tiles at that row. This is desirable. It is the mechanism that enables the threads to accommodate to the slight irregularities of the patch segment, while still be vertically aligned. In addition, this same mechanism allows the threads to skip tiles containing blinding stars (to be implemented in a future release).

The input patch segment is bounded, by construction, to a lat-long rectangle. Therefore, in general, all the threads traverse the full extent of the patch segment, from top to bottom (the exception being the cases truncated by the RoI boundary). So, the thread endpoints (top and bottom) are expected to be close to each other meridian-wise. The final step of the computation is to pairwise connect the threads from right to left. This last step creates larger moves in longitude, detrimental for PSF stability. However, these moves are limited to the number of threads, which is much lower than the length of the patch segment).

The diffusion algorithm is by far the most expensive part of the EWS computation, with the combinatorial exploration taking most of its cost. However, it has so far proven to be stable and robust, solving a large number of diverse patch segment configurations while optimising thermal stability, thus yielding a much more stable PSF. Figure 43 shows a comparison of the resulting AA field to field variation in the same patch, when scheduling it with the look-ahead algorithm versus diffusion algorithm. The time behaviour of AA in the latter case is much smoother and with fewer spikes than the previous results. This improvement translates into a better thermal stability and an overall decrease in time variations in the PSF, which therefore can be better modelled.

thumbnail Fig. 43.

Variation in AA when observing the same EWS patch with the two different algorithms: look-ahead (left panel) and diffusion (right panel). The latter has a comparatively more monotonic progression with much smaller field-to-field AA variations and spikes. Bottom panels: are the time derivative of the top panels.

7.5.3. Linking patch segments into a patch

The diffusion algorithm proposes a sequence of observations for a pattern of tiles, which is a segment of a patch. This is the building block for computing a patch. First, the diffusion algorithm is applied to all patch segments of a patch source. Then, if all are successful, the resulting patch segments are linked together in a single observation sequence, by assigning timestamps for observation, generating a single sequence that covers the time slot of a survey window (i.e. a patch; see Fig. 44).

thumbnail Fig. 44.

Example of three sibling patch segments (cf. Fig. 40). Patch segments are displayed in alternating colours. Previously observed tiles are displayed in dark blue. Red arrows represent large slews (moving from patch segment to patch segment).

Timestamps are assigned adding observation and slew times along the order of the sequence, checking compliance with constraints. This is performed first forwards in time, starting at the pointing (quaternion and timestamp) that defines the start of the survey window. Next, the timestamp of the first observation is computed by adding the slew time from the start pointing. Then, the timestamp of the second time observation is computed by adding the ROS observation time (a fixed value) plus the slew time from the previous observation (a variable value). The process is repeated throughout the sequence, assigning timestamps sequentially. In this process, slew-times are computed according to the slew-time estimator. By construction, slews within a patch segment are expected to be small slews. In contrast, jumps between patch segments are considered to be large slews, adding to the large-slew budget. However, this number of large slews is relatively small. Typically, survey windows are filled with patches generated from one or two patch sources, with each patch being split in a low number of patch segments, at most. Hence, the initial ∼100 survey windows produce no more than a few hundred patch segments in the end.

As explained above, the algorithm first attempts to assign timestamps from the beginning of the survey window forwards. If this succeeds, it creates a patch flushed backwards in time, leaving some idle time at the end of the window. This is because in general, an integer number of observations with varying slew-times does not fit perfectly into the slot of time previously defined by a survey window. If this succeeds, this is the preferable solution. If it fails, then a reverse assignment is attempted, assigning timestamps from the end of the survey window backwards (and reversing the computations of slew-time). If this succeeds, it may leave some idle-time at the beginning of the window.

Usually, there is a sufficient number of tiles distributed along the range of longitudes (covered by a patch source) to generate a patch covering all the time slot of a survey window. However, in some cases, such as when a patch source intersects only a corner of a quadrant, it is not possible to have a path extending the full width of the respective survey window. Trying to flush a patch to both end-sides of a survey window enables patching those odd cases, promoting also the generation of patches that adapt to the boundary of the RoI.

In the process of assigning timestamps, an observation of a field with a given timestamp may fail to comply with the constraints of SAA and/or AA. If, at any point, a failure of compliance is encountered, the computation is terminated and all segments are discarded. The rationale for this strategy is simple: if a sequence, computed as parallel to meridians as possible, fails the timestamp assignment, then it might not be viable in the first place. Most likely, in those cases, the patch has some geometry feature that stretches the scheduling flexibility too far. The result of this stage is a patch that is flushed, time-wise, to the beginning or to the end of the interval of time allotted to it. The path is thus a composite structure, made by a string of segments.

7.6. Placing patches

The EWS is computed sequentially, covering the RoI by placing one patch at a time. In some way, it is like solving a jigsaw puzzle on a surface of a sphere, but where the shapes of the pieces are not fixed from the start, instead being computed as the scheduling progresses. Two objectives guide the buildup of the EWS, which can be summarised as: ‘observe as much as possible, as early as possible’. More specifically: (i) we want a compact footprint (per quadrant), so that the footprint of observations forms a single, continuous region, implying that patches should match perfectly next to each other, with no holes in between10; and (ii) we want to observe fields as early as possible, which means that all slots in survey windows are assigned to observations, provided there is an unobserved part of the RoI within reach. We note that, given the cyclic nature of the process, after a few iterations, the RoI within reach of some of the survey windows is typically scheduled.

Eventually, either no area is left to schedule (within a range of latitudes of the RoI) or there is no time left (within the six years of the mission). When no suitable area can be scheduled for a given period, the schedule just leaves it unassigned, giving rise to unallocated time. This highlights that presently it is not possible to assign it to particular observations, but that it can be done at a later time.

The EWS is computed by filling the RoI one survey window at the time, in chronological order. At each step, a survey window is processed by assigning observations to its time slot. As explained below, this may require a few iterations. Only when all possible assignments are handled, the computation moves on to the next survey window.

As discussed in Sect. 7.4.3, each survey window generates one or more patch sources, which in turn may generate one patch (made of a single patch segment or of a string of several patch segments), flushed to the beginning or to the end of its survey window. All the patches are continuously linked to a particular stage-1 observation (e.g. a calibration, deep-field, or polar cap patch). All these patches compete for the same slot of time, and we need to decide which patch to add to the survey. In general, the choice is to build the EWS layer by layer of ecliptic latitude, striving to go from high to low S/N regions. When this criteria is not decisive, the choice is to select the patch that better matches (horizontally, along ecliptic latitudes) some previously selected adjacent patch.

The last step in the processing of a survey window is to place PSF calibrations and SOPs. At most, one PSF calibrations is inserted per window. This is achieved by analysing the patches just selected, identifying the timestamps where a PSF calibration target is within reach. A jump to a PSF calibration takes place at the end of some wide observation. Typically, a single patch generates several such candidates. In parallel, candidate timestamps for the required SOPs are identified. Again, these are timestamps of the end of some wide observation. There are one to two SOP candidates per window, at most. Then, the two types of candidates are considered, inserting them in chronological order; one PSF calibration, if available, and one or two required SOPs. Among all PSF candidates, preference is given to the ones occurring between patch segments or between a preceding SOP and a following patch segment, in order to save large slews. The insertion is performed by cutting a wide patch at the required timestamp, adding the PSF calibration or the SOP, and pushing the remaining of the patch forwards in time. This process may require also an adjustment of the following calibration block, pushing it forwards in time (if now overlapped by the preceding patch). The process of SOP insertion is also applied to the following calibration block before processing the next survey window.

The islands are located at lower latitudes and comprise less area. Consequently, most of the survey buildup is shared between the two mainlands (I and IV; see Fig. 17) with the islands (II and III) becoming relevant in the last year of the survey. This can be seen in Fig. 45, which shows an example of the buildup of EWS patches over the duration of the mission, coloured from deep blue to light blue as time progresses. This strategy does not only forces the two mainlands to grow at an equal rate but, more importantly, it is crucial to guarantee a compact survey footprint.

thumbnail Fig. 45.

Example of the ECTile output, showing the buildup of the coverage of the EWS in RSD_2021A (cf. Sect. 8) over time in ecliptic coordinates. In the main panel the EWS patches are coloured according to different years, in correspondence to the time bars below. The part of the RoI that is not observed is in light green (we recall that the RoI covers an area larger than the one needed to be observed). Under the main panel, the progression over the six years is unfolded, with two rows for each year indicating the trailing and leading directions of the telescope. Different colours represent the progression in time, from deep blue to light blue, whereas the thin line segments represent time intervals reserved for calibrations and deep fields. The light grey boxes for the leading or trailing directions mark when the telescope is actually observing in the opposite direction (for each grey box there is a corresponding coloured box). Towards the end of the mission, some of the survey windows are completely or partially light red, showing that that particular time slot has run out of sky areas to observe.

The height of a patch is dictated by the speed of the orbit, which varies slightly with latitude. It roughly corresponds to the time the orbit takes to scan the width of a FoV divided by the nominal observation time. To some degree, the height (and width) may vary slightly from its natural size (approximately, ±2 rows). Hence, in practice the height of a patch is more or less fixed, making the layer approach optimal. It guarantees a maximum of free RoI above previously computed patches, giving ample space (i.e. height) for the generation of each new patch. The sole exception is when the top layer of the survey reaches the boundary of the RoI. If the available height is less than the minimum patch height, the patch generation fails. Due to the layering approach, this obstacle arises only for a small part of the RoI, namely the top layer of each quadrant. For the rest of the RoI, the stacking of patches ensures a compact filling.

8. The Euclid reference survey definition

In this section we present the latest version of the Euclid RSD, RSD_2021A. This is the result of the stage-1 and stage-2 scheduling procedures (cf. Sects. 6.3 and 7), and the corresponding ECTile outputs are the ones shown in Figs. 29 and 45.

We recall that the RSD observations start three months after launch. They are preceded by a one-month commissioning phase, followed by a two-month PV phase. During the PV phase the first survey-like data will be obtained, which are used to verify the data processing, and to validate and eventually tune or adjust the nominal sequence of operations (possible minor changes to the RSD can be implemented in less than one week time, producing a new version of the RSD; more complex changes of course require more time). Moreover, during PV some survey specific observations will be carried out. The latter will, for instance, provide direct measurements of the zodiacal light and stray light to verify and refine our models.

8.1. Survey performance

Figure 46 shows the footprint of RSD_2021A The different colours indicate different observing epochs of the EWS. The three EDFs and the six EAFs (cf. Sect. 6.2) are tied into the EWS. In the two mainlands of the EWS, the observations start from the ecliptic poles and progress towards the equator. The best sky areas around both Galactic caps are covered within the first three years of the mission. The observations of the two EWS islands only take place in the final two years of the mission. We note that some of the worst parts of the RoI (totalling an area of ∼2696 deg2) are left unobserved. These correspond to the uncoloured areas in the islands and at low latitudes in the mainlands. The areas of the sky with longitudes between 150° and 225°, and between 330° and 45°, are observable at the same time since they are separated by 180°. They contain much area within the RoI, and moreover the EDF-F and EDF-S are also located there. This means that there is not enough time for the EWS to observe all that area in the six years of the mission and hence their worst-quality regions are not observed.

thumbnail Fig. 46.

RSD 2021A (14 514 deg2) chronology shown in celestial coordinates. RoI boundaries are shown as solid red lines. Blinding stars cause 809 avoidance areas within the reference survey, with an average of 0.785 deg2 per avoidance area, totalling 635 deg2. Dashed lines (1300 deg2 in white per Galactic cap) delimit the highest S/N areas. The ecliptic referential is over-plotted in red. The three EDFs (bright green) and the six EAFs (red diamonds, not to scale) are shown.

The RSD contains 44 065 fields (28 080 to build the EWS and 15 985 for EDFs, EAFs and calibration targets observations). The EWS fields are contained in 256 patches (seen in Figs. 45 and 46). The vast majority of the field slews, used to point the telescope, are below 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $, as shown in the right panel of Fig. 47. This is the most efficient slew regime in terms of propellant usage. As shown in Fig. 47, all telescope rotations are done within the allowed SAA and AA limits. Most of the observations are done close to transit, with 90% of the SAA values used between 88° and 94°. The statistics of AA usage shows that 71% of the telescope rotations are done with |AA|< 1°. Even though SAA and AA values spread over the full range allowed, the field-to-field variations (between consecutive observations) of SAA and AA are very small throughout the survey: smaller than 1° in 97.4% (SAA) and 98.6% (AA) of the field-to-field transitions over the full mission. This feature is extremely important for the thermal stability, which ensures a stable PSF for WL shape measurements. It was possible to achieve this performance thanks to the implementation of the diffusion algorithm described in Sect. 7.5.2.

thumbnail Fig. 47.

Distributions of angles and slews. Left panel: AA angle. Central panel: SAA angle. Right panel: field slews up to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $. All quantities comply with the allowed ranges.

8.2. Unallocated time

The existence of a deficit of area on some longitudes (see Sect. 5.2.1), compared to the available observing time, is evident from Fig. 48. The blue curve is the area available in the RoI at a given ecliptic longitude (in bins of 1°). The RoI areas in longitudes separated by 180° are added, since that pair of longitudes can be observed at the same time, from the trailing or the leading direction. Due to this six-month periodicity, the x-axis range only extends to 180°. The red curve denotes the cumulated number of days during which a given longitude is visible for EWS observations, assuming transit observations, and converted to equivalent area (1 day corresponding to 10 deg2). The available time is not uniform, it is determined after the stage-1 schedule is defined (see Sect. 6.3), which creates a strong variation along the year (i.e. wiggles in the red curve). For example, the absolute minimum corresponds to the highly booked longitudes of the EDFF and EDFS, where less time is left for EWS observations.

thumbnail Fig. 48.

Area of the RoI as a function of (folded) longitude (blue). The equivalent area is visible as a function of (folded) longitude (red). The area is the available RoI area at transit per degree of longitude. The red curve assumes that the EWS observations are done at transit. Unallocated time arises when the red curve is above the blue curve. See text for details.

In longitudes where the red curve is above the blue curve, there is a deficit of area for the time available for EWS, leading to unallocated time. Conversely, in longitudes where the blue curve is above the red curve, there is an excess of area for the time available for EWS, leading to unobserved areas in the RoI. In Fig. 46 this corresponds to the areas with no patches, which clearly are on the areas of the RoI of worst quality.

We note that the presence of unallocated time in the EWS schedule does not mean that there will be any idle time, because some areas of the EWS may be re-observed or new areas that do not qualify for the EWS, but have scientific value nonetheless, may be observed instead. In doing so, we can either consider fields that are observable within the thermal and pointing constraints enforced for the EWS, or we can operate outside these constraints, thus with the risk of perturbing the continuation of the EWS afterwards. Therefore, one needs to have the real in-flight characteristics to get a solid picture of the possibilities and constraints.

Some examples that maintain the strict survey limitations are: (i) possible multiple exposures (4) on areas in the ecliptic plane that do not qualify for S/N with ROS single visit; (ii) possible decontamination procedures; (iii) repeat some suitable but lower-quality EWS regions to boost their S/N, or fill possible gaps due to unexpected events that might interrupt the basic scheduling; (iv) increase the depth of the self-cal field, which could become a reference field for a dedicated supernova programme in a possible extension of the mission; (v) build a Euclid medium deep field (EMDF) covering 200 − 300 deg2; (vi) that is one magnitude deeper than the EWS (this would require five additional visits that would preferentially be done once every year); observe suitable astronomical objects that are located outside the RoI that would benefit from localised re-visits, such as observations on the ecliptic or Galactic plane, Galactic bulge (e.g. microlensing), specific low surface brightness objects, nearby galaxies, or clusters of galaxies; (vii) use of the blue grism on targets during new or repeated visits on specific targets. How to best use the time that cannot be used for single pass EWS will be decided at a later stage upon a consolidated scenario by the Euclid Science Team and ESA.

8.3. Euclid footprint area

The area of the RSD_2021A Euclid footprint is 14 514 deg2. Table 9 lists the observed EWS area at the end of each year of the mission, while the growth with time of the area covered by the EWS is shown in Fig. 49.

thumbnail Fig. 49.

Growth of the area surveyed by the RSD_2021A EWS with time (upper curve), and growth of unallocated time with time (lower curve), here in units of EWS equivalent area (one month corresponds to about 300 deg2).

Table 9.

EWS area observed in RSD_2021A per year and cumulatively at the end of each year. The last column shows the cumulative unallocated time (see text).

During the first year many calibration observations are scheduled to support the first data releases. As a result, the EWS initially progresses slower than in the second and third year, but it still reaches an area in excess of 2500 deg2. At the end of the third year, more periods occur when standard EWS observations cannot be made due to the increasing paucity of available unobserved areas within the RoI; the slope of the growth in time flattens there, causing the staircase-like pattern seen in Fig. 49. These periods of ‘unallocated time’ (see Sect. 8.2) have a 6-month periodicity due to the intersection of the Galactic plane with the ecliptic plane, have an increasing duration and slow down the progression of the EWS.

The time available to implement the EWS, once the other mandatory observations are carried out and the geometry of the RoI is taken into account, is an important input for the construction of the EWS and the final covered areas. In RSD_2021A this time is 2190 − 427 − 39 − 292 = 1432 days, where 427 days are used for calibrations, EDF and EAF observations, 39 days are reserved for SOP and there are 292 unallocated days (Fig. 50 depicts the time allocation breakdown). Given the length of the ROS, close to 4400 s per field (Sect. 4.1), and the FoV of 0.53 deg2 (Sect. 3.3), the effective EWS area increases by 10.1 deg2 per day, (considering the effective average field overlap of 3%). This means that, with the current calibration, EDFs, EAFs, and SOP required times, the EWS can reach an area of 15 000 deg2 only if the unallocated time is shorter than 240 days.

thumbnail Fig. 50.

Breakdown of the total time of the Euclid mission (2190 days) in EWS and non-EWS usage, for RSD_2021A, corresponding to the schedule depicted in Fig. 29.

The time available for the EWS has generally decreased over the years, as the mission matured. Figure 51 depicts this evolution, indicating the driving changing factor for each transition (modifications in calibrations, RoI or ROS). In earlier surveys, such as the RSD_2015A prepared for PDR, the time needed for calibrations was small and the RoI was larger (driven by WL counts). In early 2018, RSD_2018A used for SPV2 introduced a smaller RoI based on the more restrictive conditions for GC, which led to a larger distance from the Galactic plane (because of effects from star density on spectra), with a smaller area to be observed and an increase in unallocated time (cf. Sect. 8.2). Shortly after, a major redefinition of the calibration plan, increased its allocated time by 100 days. At the same time the redefinition of the slew concept (see Sect. 7.4.2) made it no longer possible to schedule the EWS at high latitudes with the global tessellation, leading to a faster filling of the EWS and a faster buildup of unallocated time. These two factors led to a strong decrease in the EWS available time in RSD_2018B for CDR (Laureijs et al. 2020). In 2019, RSD_2019A introduced a new algorithm to schedule the high-latitude regions. Given the fact that the throughputs were measured to be larger than specified, the RoI was reverted to the limits defined in Laureijs et al. (2011). The larger area of this RoI led to an increase in the available time for the EWS, and to a solution that covered 15 000 deg2. The year 2020 introduced the new, better high-quality RoI, described in Sect. 5, based on the latest estimates of S/N. This again decreased the available EWS time, now causing the RSD_2020A to barely reach 15 000 deg2. In RSD_2020B the diffusion algorithm was introduced to increase thermal stability. This was achieved without decreasing the schedule efficiency and the available time for the EWS. Nevertheless the EWS time decreased due to a further increase in time needed for the calibrations. Finally, in RSD_2021A we implemented the skipping of bright stars, together with the latest revised times for both ROS (significant increase in the duration of the dither steps) and SOP (halved).

thumbnail Fig. 51.

Evolution of the time available for EWS observations in the RSDs produced from 2015 to 2021. The changes are due to the increase or decrease (denoted by the symbols + or −; larger changes have double symbol) in calibrations, RoI area (directly impacting the unallocated time), and ROS time. See text for details.

RSD_2021A falls short of covering the desired 15 000 deg2 area for the EWS by ∼3%. The missing 486 deg2 are the equivalent of one and a half months of EWS observing time. In part this can be recovered by making local tessellations, to allow shifting the position of the centre of the tiles affected by stars, so as to avoid the blinding stars falling on the detectors, instead of skipping the whole tile. Moreover, it is expected that the EWS available time will further decrease sightly with the insertion of the complete set of PSF calibrations. In fact, the expected cadence is not yet fully respected in RSD_2021A, because the stringent stability constraints could not always be fulfilled, especially in the final years of the survey. Study is underway to tackle this issue and, once this problem will be solved, the time needed amounts to ≈10 days (100 deg2) to be taken off from the EWS present coverage. It must be noticed, however, that a 3% reduction in survey area translates into a comparable reduction in the final DE FoM (actually even less since the still uncovered areas of the RoI are of lower quality than the average). An area a few percent lower than originally expected is therefore of no consequence for achieving the original goal of the mission, that is, a FoM larger than 400 from the two main probes for a w0 − wa CDM model with no other prior information (Laureijs et al. 2011). Indeed, for the present coverage, one preliminarily expects the FoM to scale down from 500 (Euclid Collaboration 2020a) to ≈480.

Nevertheless, the final area of the Euclid footprint is not yet fixed and its value can change in either direction. On the one hand, possible further synergies between targets and optimisations will likely increase the area covered by the EWS. On the other hand, modifications in the ROS, the calibrations plan, and in the amount of overlap between tiles might further reduce the final area.

9. Summary

We have presented the status of the reference survey for the Euclid mission at the beginning of the year 2021. The reference survey encompasses all six years of the mission baseline, fulfilling all the constraints while combining observations of the EWS, EDFs, EAFs and other calibration fields, as well as all the slews. We discussed and provided models of the main backgrounds that affect the space observations. We also presented the complementary ground-based observations needed for photometric redshifts.

This is a non-trivial achievement, because of the complex nature of the mission: the Euclid step-and-stare strategy is severely constrained by pointing limitations with respect to the Sun, whilst it also has to carry out a large number of dedicated observations for calibration purposes and sample characterisation. This leads to a complex interplay of timing observations, visibilities, and manoeuvres. Although the two main instruments are operated in concert, they require three separate calibration strategies, thus adding to the complexity. This makes Euclid not only different from a typical observatory mission (e.g. Herschel) but also from all-sky missions that have a single scanning strategy (e.g. Planck and Gaia).

Despite these challenges, we have found a highly optimised solution that takes into account the main spacecraft and instrument characteristics and limitations, the current models for the expected background, the various calibrations, the dithering pattern, and the methods to cover the wide area expected to be observed. The latest version of the Euclid RSD (RSD_2021A) fulfils practically all requirements, resulting in an EWS that covers ≈14 500 deg2 of the extragalactic sky. This is 3% short of the initial target, mainly because of the paucity of good sky for Euclid and of the severe constraints on the pointing (if feasible, 10 months of presently unallocated time would in principle allow 3000 square degrees to be added to the EWS). A companion paper will detail the rationale and results for the EDFs that cover about 40 deg2 but at much greater depth.

The definition of the reference survey is an ongoing process, and some of the results presented will continue to evolve because the Euclid Consortium has developed both the expertise and the specific tools that allow one to probe significant variations in any of the multiple boundary conditions that originate from the knowledge of the spacecraft, the instruments, or the astronomical environment. The current solution is, however, a close proxy for the final survey. Crucially, it demonstrates the feasibility of the core mission within the numerous constraints. Future plans include:

  • further optimisations to reach the nominal 15 000 deg2 (e.g. double exposures on the ecliptic plane to compensate for the local high background);

  • implementation of refined simulations;

  • implementation of refined background models;

  • incorporation of possible future changes in parameters (e.g. changes in ROS or in the calibration plan);

  • consideration of ‘what if’ scenarios for non-recurrent operations (decontamination, phase diversity calibrations) and possible failures (electronics or other systems);

  • revisions and updates based on real measurements and in-flight performance.

These further improvements will present challenges in their own right, but given the current level of maturity, there is little doubt that Euclid will dramatically advance our understanding of the nature of DM and DE whilst impacting many aspects of astronomy thanks to the tremendous legacy value of these unique data.


1

This is a revision of the Laureijs et al. (2011) value after the removal in 2014 of blue-grism exposures from the wide survey.

2

In three-dimensional space, according to Euler’s rotation theorem, any rotation or sequence of rotations of a rigid body or coordinate system about a fixed point is equivalent to a single rotation by a given angle θ about a fixed axis (called the Euler axis) that runs through the fixed point. Therefore, the quaternion fully describes the spacecraft attitude and a single rotation quaternion relates one pointing to another.

3

In particular Ebv_xgal_ns2048_REL5.fits, found at http://hyperstars.lmpa.eu/mamd/planck_dust_model.html.

8

A summary of the SPV2 exercise can be found in pages 4 and 5 of the EC newsletter at: https://www.euclid-ec.org/Documents/Newsletter/EC-Newsletter_issue08.pdf

9

The properties of the three EDFs can be found at https://www.cosmos.esa.int/web/euclid/euclid-survey

10

The only planned ‘holes’ are the tiles skipped because of the presence of blinding stars.

Acknowledgments

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 website (http://www.euclid-ec.org).

References

  1. Albrecht, A., Bernstein, G., Cahn, R., et al. 2006, ArXiv e-prints [arXiv:astro-ph/0609591] [Google Scholar]
  2. Aldering, G. 2001, SNAP Sky Background at the North Ecliptic Pole, www-supernova.lbl.gov/~aldering/zodi.pdf [Google Scholar]
  3. Aldering, G., Akerlof, C. W., Amanullah, R., et al. 2002, Proc. SPIE, 4835, 146 [NASA ADS] [CrossRef] [Google Scholar]
  4. Amendola, L., Appleby, S., Avgoustidis, A., et al. 2018, Liv. Rev. Rel., 21, 2 [Google Scholar]
  5. Amiaux, J., Scaramella, R., Mellier, Y., et al. 2012, Proc. SPIE, 8442, 84420Z [NASA ADS] [CrossRef] [Google Scholar]
  6. Bagley, M. B., Scarlata, C., Mehta, V., et al. 2020, ApJ, 897, 98 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bertin, E., & Arnouts, S. 1996, A&AS, 317, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. BOSS Collaboration (Alam, S., et al.) 2017, MNRAS, 470, 2617 [NASA ADS] [CrossRef] [Google Scholar]
  9. Buffington, A., Bisi, M. M., Clover, J. M., et al. 2016, Icarus, 272, 88 [NASA ADS] [CrossRef] [Google Scholar]
  10. Cenarro, A. J., Ederoclite, A., & Íñiguez, C. 2018, Proc. SPIE, 10700, 107000D [NASA ADS] [Google Scholar]
  11. Chambers, K., Magnier, E., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
  12. Costille, A., Beaumont, F., Prieto, E., et al. 2016, Proc. SPIE, 9912, 99124I [NASA ADS] [CrossRef] [Google Scholar]
  13. Cropper, M., Hoekstra, H., Kitching, T., et al. 2013, MNRAS, 431, 3103 [Google Scholar]
  14. Cropper, M., Pottinger, S., Niemi, S.-M., et al. 2014, Proc. SPIE, 9143, 91430J [NASA ADS] [CrossRef] [Google Scholar]
  15. Cropper, M., Pottinger, S., Niemi, S., et al. 2016, Proc. SPIE, 9904, 99040Q [Google Scholar]
  16. de Bernardis, P., Ade, P. A. R., Bock, J. J., et al. 2000, Nature, 404, 955 [Google Scholar]
  17. de Bernardis, P., Ade, P. A. R., Bock, J. J., et al. 2002, ApJ, 564, 559 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dark Energy Survey Collaboration (Abbott, T. M. C., et al.) 2019, Phys. Rev. D, 99, 123505 [Google Scholar]
  19. Dark Energy Survey Collaboration (Abbott, T. M. C., et al.) 2022, Phys. Rev. D, 105, 023520 [CrossRef] [Google Scholar]
  20. Er, X., Hoekstra, H., Schrabback, T., et al. 2018, MNRAS, 476, 5645 [NASA ADS] [CrossRef] [Google Scholar]
  21. Euclid Collaboration (Blanchard, A., et al.) 2020a, A&A, 642, A191 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Euclid Collaboration (Desprez, G., et al.) 2020b, A&A, 644, A31 [EDP Sciences] [Google Scholar]
  23. Gil-Marin, H., Bautista, J., Paviot, R., et al. 2020, MNRAS, 498, 2492 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gómez-Alvarez, P., Dupac, X., Buenadicha, G., et al. 2018, SPIE Conf. Ser., 10707, 1070712 [Google Scholar]
  25. Gordon, K. D. 2004, ASP Conf. Ser., 309, 77 [NASA ADS] [Google Scholar]
  26. Gordon, K. D., Clayton, G. C., Misselt, K. A., et al. 2003, AJ, 594, 279 [NASA ADS] [CrossRef] [Google Scholar]
  27. Hanany, S., Ade, P., Balbi, A., et al. 2000, ApJ, 545, L5 [Google Scholar]
  28. Hildebrandt, H., Köhlinger, F., van den Busch, J. L., et al. 2020, A&A, 633, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [Google Scholar]
  30. Hoekstra, H., Viola, M., & Herbonnet, R. 2017, MNRAS, 468, 3295 [NASA ADS] [CrossRef] [Google Scholar]
  31. Huchra, J., Macri, L., Masters, K., et al. 2012, ApJS, 199, 26 [Google Scholar]
  32. Ibata, R., McConnachie, A., Cuillandre, J.-C., et al. 2017, ApJ, 848, 128 [NASA ADS] [CrossRef] [Google Scholar]
  33. Ivezić, Ž., Kahn, S., Tyson, A., et al. 2019, ApJ, 873, 111 [NASA ADS] [CrossRef] [Google Scholar]
  34. Jamal, S., Le Brun, V., Le Fèvre, O., et al. 2018, A&A, 611, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Kelsall, T., Weiland, J. L., Franz, B. A., et al. 1998, ApJ, 508, 44 [Google Scholar]
  36. Kilbinger, M. 2015, Rep. Progr. Phys., 78, 086901 [Google Scholar]
  37. Korsch, D. 1972, Appl. Opt., 11, 2986 [NASA ADS] [CrossRef] [Google Scholar]
  38. Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv e-prints [arXiv:1110.3193] [Google Scholar]
  39. Laureijs, R., Racca, G. D., Mellier, Y., et al. 2020, Proc. SPIE, 11443, 114430F [NASA ADS] [Google Scholar]
  40. Leinert, C., Bowyer, S., Haikala, L. K., et al. 1998, ApJS, 127, 1 [Google Scholar]
  41. Lupton, R., Gunn, J., & Szalay, A. 1999, AJ, 118, 1406 [NASA ADS] [CrossRef] [Google Scholar]
  42. Maciaszek, T., Ealet, A., Jahnke, K., et al. 2014, Proc. SPIE, 9143, 91430K [NASA ADS] [CrossRef] [Google Scholar]
  43. Maciaszek, T., Ealet, A., Jahnke, K., et al. 2016, Proc. SPIE, 9904, 99040T [NASA ADS] [CrossRef] [Google Scholar]
  44. Maris, M., Burigana, C., & Fogliani, S. 2006, A&A, 452, 865 [Google Scholar]
  45. Maris, M., Scaramella, R., Burigana, C., et al. 2016, ASP Conf. Ser., 521, 531 [NASA ADS] [Google Scholar]
  46. Markovič, K., Percival, W. J., Scodeggio, M., et al. 2017, MNRAS, 467, 3677 [Google Scholar]
  47. Massey, R., Hoekstra, H., Kitching, T., et al. 2013, MNRAS, 429, 661 [Google Scholar]
  48. Mihos, J. C., Harding, P., Spengler, C., Rudick, C., & Feldmeier, J. J. 2013, ApJ, 762, 13 [NASA ADS] [CrossRef] [Google Scholar]
  49. Miville-Deschênes, M.-A., Duc, P.-A., Marleau, F., et al. 2016, A&A, 593, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Miyazaki, S., Komiyama, Y., Kawanomoto, S., et al. 2018, PASJ, 70, S1 [NASA ADS] [Google Scholar]
  51. Pamplona, T., Gimenez, J. L., Febvre, A., et al. 2016, SPIE Conf. Ser., 9912, 99120M [NASA ADS] [Google Scholar]
  52. Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJ, 517, 565 [Google Scholar]
  53. Pickles, A., & Depagne, É. 2010, PASP, 122, 1437 [NASA ADS] [CrossRef] [Google Scholar]
  54. Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Planck Collaboration I. 2020, A&A, 641, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Pozzetti, L., Hirata, C. M., Geach, J. E., et al. 2016, A&A, 590, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Pryke, C., Halverson, N. W., Leitch, E. M., et al. 2002, ApJ, 568, 46 [NASA ADS] [CrossRef] [Google Scholar]
  59. Pyo, J., Matsumoto, T., Jeong, W.-S., et al. 2012, ApJ, 760, 102 [NASA ADS] [CrossRef] [Google Scholar]
  60. Racca, G. D., Laureijs, R., Stagnaro, L., et al. 2016, Proc. SPIE, 9904, 99040O [Google Scholar]
  61. Rassat, A., Amara, A., Amendola, L., et al. 2008, ArXiv e-prints [arXiv:0810.0003] [Google Scholar]
  62. Refregier, A., & Douspis, M. 2008, Proc. SPIE, 7010, 701018 [NASA ADS] [CrossRef] [Google Scholar]
  63. Rhodes, J., Nichol, R. C., Aubourg, É., et al. 2017, ApJS, 233, 21 [Google Scholar]
  64. Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009 [Google Scholar]
  65. Robin, A., Marshall, D., Schultheis, M., & Reylé, C. 2012, A&A, 538, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. Robin, A., Reylé, C., Fliri, J., et al. 2014, A&A, 569, A13 [CrossRef] [EDP Sciences] [Google Scholar]
  67. Rowe, B., Hirata, C., & Rhodes, J. 2011, ApJ, 741, 46 [NASA ADS] [CrossRef] [Google Scholar]
  68. Scaramella, R., Mellier, Y., Amiaux, J., et al. 2015, IAU Symp., 306, 375 [NASA ADS] [Google Scholar]
  69. Semboloni, E., Hoekstra, H., & Huang, Z. 2013, MNRAS, 432, 2385 [NASA ADS] [CrossRef] [Google Scholar]
  70. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  71. Tereno, I., Carvalho, C. S., Dinis, J., et al. 2015, IAU Symp., 306, 379 [NASA ADS] [Google Scholar]
  72. Tonry, J. L., Denneau, L., Flewelling, H., et al. 2018, ApJ, 867, 105 [NASA ADS] [CrossRef] [Google Scholar]
  73. Venancio, L. M. G., Pachot, C., Carminati, L., et al. 2016, Proc. SPIE, 9904, 99040P [NASA ADS] [Google Scholar]
  74. Venancio, L. M. G., Carminati, L., Amiaux, J., et al. 2020, Proc. SPIE, 11443, 114430G [NASA ADS] [Google Scholar]
  75. Viola, M., Joachimi, B., & Kitching, T. 2014, MNRAS, 439, 1909 [NASA ADS] [CrossRef] [Google Scholar]
  76. Weinberg, S. 1989, Rev. Mod. Phys., 61, 1 [NASA ADS] [CrossRef] [Google Scholar]
  77. Wright, N. 1998, ApJ, 496, 1 [NASA ADS] [CrossRef] [Google Scholar]
  78. Zakharov, A., Prokhorov, M., Tuchin, M., & Zhukov, A. 2013, Astrophys. Bull., 68, 4 [Google Scholar]
  79. Zoubian, J., Kümmel, M., Kermiche, S., et al. 2014, ASP Conf. Ser., 485, 509 [NASA ADS] [Google Scholar]

Appendix A: Acronyms

Table A.1.

Acronyms used in the paper.

All Tables

Table 1.

Sizes of the VIS and NISP focal planes and their corresponding FoVs, as well as the angular sizes of detectors and detector gaps and the width of the VIS charge injection lines (see Fig. 2).

Table 2.

Statistics (mean and standard deviation) on the number, X, of exposures per pixel for the ‘S’ dither pattern: percentages of covering for individual (X = 0, 1, 2, 3, 4) and cumulative (X > 4, X ≥ 3, X ≥ 4) bins.

Table 3.

Ecliptic latitude dependence of the simple zodiacal background model of Leinert et al. (1998).

Table 4.

Normalisation of the basic zodiacal model at the average wavelength of Euclid bands (see Fig. 7).

Table 5.

Extinction scaling coefficients Aλ for the Euclid passband central wavelengths and total extinction A at the NEP.

Table 6.

Four separate quadrants of the RoI (see Fig. 17).

Table 7.

S/N statistics for the RoI, for each channel: IE band (boldface values refer to extended objects), NIR bands (JE, YE, and HE values refer to point-like objects), and red grism band, RGE (italic values in the last column refer to 0 . 5 $ 0{{\overset{\prime\prime}{.}}}5 $ diameter sources).

Table 8.

Time budget for calibration, EDF, and EAF observations.

Table 9.

EWS area observed in RSD_2021A per year and cumulatively at the end of each year. The last column shows the cumulative unallocated time (see text).

Table A.1.

Acronyms used in the paper.

All Figures

thumbnail Fig. 1.

Euclid timeline. The EWS has to be carried out within the six year mission baseline and will start 3 months after the launch, following a commissioning period of 1 month and a PV period of 2 months. An extension of operations beyond the six years is possible and will be decided in due time. During the first six months of the main mission (the early operations phase) a faster replanning of the survey is allowed. The three main data releases (DR#) are shown. The plan is to have 2500 deg2 made public in DR1, to grow to 7500 in DR2, and be complete at DR3 for 15 000 deg2. In addition, four quick data releases (Q#) are foreseen, each of ∼50 deg2.

In the text
thumbnail Fig. 2.

Euclid reference frames. Left panel: Euclid spacecraft reference frame. XSC points towards the Sun disk centre. The edge and centre of the joint FoV are offset by 0 . ° 47 $ 0{{\overset{\circ}{.}}}47 $ and 0 . ° 82 $ 0{{\overset{\circ}{.}}}82 $, respectively. The longer side of the FoV is typically aligned with ecliptic meridians during observations. Right panel: Euclid optical reference frame. We note that the XoptYopt plane is defined as looking onto the sky with z ̂ $ \hat z $ towards the spacecraft. There are four additional chips used as fine guidance sensors that are placed on each side of the VIS FoV.

In the text
thumbnail Fig. 3.

Illustration of the relevant pointing angles defined in the text. In the top right image the minimal SAA (pointing towards the Sun) is indicated, while the maximum SAA (pointing away from the Sun) adopted in the survey is 110° (see also Sect. 2.5). In the lower right image, AA is shown; the maximum range allowed for the survey is ±5°, which corresponds to a large margin with respect to the spacecraft capabilities.

In the text
thumbnail Fig. 4.

Euclid’s orbit and operation mode. Left panel: Euclid’s transfer orbit following launch at time T0, and the subsequent Lissajous orbit at the Sun-Earth Lagrangian point L2 (SEL2). Right panel: Euclid’s main step-and-stare operating mode, showing the steps as rotations around XSC.

In the text
thumbnail Fig. 5.

Strip of sky observable by Euclid at a given epoch (in light red). The dashed blue line is the great circle orthogonal to the Sun with the satellite observing orthogonal to the Sun for a transit at 0 deg. The asymmetry of the red region with respect to the transit blue lines is due to the asymmetry of the usable SAA range with respect to 90° (wider when pointing away from the Sun).

In the text
thumbnail Fig. 6.

VIS and NISP focal plane arrays. Left panel: VIS FPA, illustrating the detector placement. The dashed line shows the joint FoV of both instruments. Two narrow strips at the extremes of the Y axis are outside the joint FoV. Right panel: NISP FPA, with two narrow strips at the extremes of the X axis outside the FoV (the reference frame is XFoV − YFoV; see also Fig. 2 and the definition at the end of Sect. 2.1).

In the text
thumbnail Fig. 7.

Total sensitivity of Euclid’s photometric IE, YE, JE, and HE, and spectroscopic BGE, RGE (solid and dashed lines, respectively). The sensitivity (in electrons per photon) includes all optical surfaces as well as the detectors’ average quantum efficiency, all considered at their expected EOL performance after six years at L2 (degraded by radiation damage and contamination; see e.g. Venancio et al. 2020).

In the text
thumbnail Fig. 8.

(Simplified) ROS executed for each field, observing four dithered pointings. Each pointing results in a IE, YE, JE, and HE image and one spectral exposure, RGE. Filter and grism wheel movements are shown, together with stabilisation times (‘Stab’). Not shown are wheel movements during the dither slew, readout and processing done simultaneously with other operations and slews, overheads, and various inline calibrations.

In the text
thumbnail Fig. 9.

Seven dither patterns that we have analysed shown in the spacecraft reference frame (see Fig. 2). For clarity, a shift of 5″ (or of −5″) in the Xsc direction is applied to the ‘J’ (or to the ‘N’) pattern in the figure so as to distinguish it from the ‘S’ pattern. Only the ‘J’ and the ‘S’ patterns, the first two steps of which are identical, are fully compliant with the stringent survey requirements.

In the text
thumbnail Fig. 10.

Illustration of the four stacked integration time maps for ‘J’ (top) and ‘S’ (bottom) dithering patterns. The left half displays the maps for VIS and the right half those for NISP. The colour scaling gives the number of stacked exposures per map pixel. Details are given in the text.

In the text
thumbnail Fig. 11.

Adopted ‘S’ dither pattern projected on the sky for a leading (left) and a trailing (right) pointing. See Fig. 2 for reference frames. The red frame is the first to be observed, and the dots mark the centre of each frame. We note how the pattern on the sky is invariant if the telescope is flipped.

In the text
thumbnail Fig. 12.

Comparison of the zodiacal emission at 1.25 μm. For a random subset of the pointings in the reference survey, we show results from the conservative basic model (red triangles) and the estimates by the more complex model (blue diamonds). The latter (Maris et al. 2016) depends on the observation time and angular direction, and it no longer has cylindrical symmetry. This causes the scatter in the blue points. The EuclidYJH bands show a similar behaviour.

In the text
thumbnail Fig. 13.

1034 blinding stars present in Euclid’s RoI. There are 1003 stars over the RoI that have Y, or J, or H < 4 AB mag (in red), and 31 stars with i < 4 AB mag and all three bands Y, J, H > 4 AB mag (in blue). The size of the dots scales down with the magnitude from −2 to +4. Those nearby bright stars are evenly distributed over the RoI, uncorrelated with the Galactic latitude. The three EDFs (see Sect. 6.2), free of blinding stars, are shown in green.

In the text
thumbnail Fig. 14.

Impact from bright stars on the S/N achievable by a target galaxy IE = 24.5 AB mag (extended source) when placed over a full Euclid FoV ( 0 . ° 7 × 0 . ° 7 $ {\sim}0{{\overset{\circ}{.}}}7 \times 0{{\overset{\circ}{.}}}7 $, with a scheme of the 6 × 6 CCD mosaic layout superimposed) in a typical part of the Euclid RoI. The colour scale is the S/N of a target source placed in the field at that position. The degradation in S/N when the target galaxy is placed closer and closer to a bright star is evident. Left panel: IE = 4 AB mag star with effects on almost 50% of the FoV, but science impacted on less than 10% of the FoV. The solid yellow area is where the target galaxy is at S/N < 14 (see text). Right panel: common stars of AB magnitude IE = 8, 9, and 10, which have a very limited impact overall with respect to the FoV size: less than 0.1% of the FoV lost on average across the RoI.

In the text
thumbnail Fig. 15.

RoI outline (17 354 deg2) with the accepted ranges of the stellar density, dust extinction, and zodiacal light.

In the text
thumbnail Fig. 16.

Distribution of reddening values (left), stellar density (middle), and zodiacal background (right) across the RoI. Below each graph, the box-and-whiskers plot marks the mean (red line), median (black line), and interquartile range (empty box), plus the minimum and maximum values.

In the text
thumbnail Fig. 17.

Sky constraints as experienced by Euclid from L2, and the RoI avoiding the worst regions. The four quadrants of the RoI are in blue. The solid yellow lines trace the ecliptic latitude threshold of |β| = 10°. The solid red lines trace the galactic latitude threshold of |b| = 23°. For declination δ ≥ +30°, the threshold is |b| = 25° to match the northern ground surveys (Sect. 5.2.1).

In the text
thumbnail Fig. 18.

Available RoI area as a function of the ecliptic longitude. The troughs are produced by the intersection of the Galactic and ecliptic planes and lead to unallocated time.

In the text
thumbnail Fig. 19.

Euclid S/N maps at detection limit in equi-rectangular celestial projection (same referential as Fig. 17). Top: VIS IE band S/N; the entire RoI is within specifications (S/N ≥ 10, the first level of light green). Middle: NISP S/N for the JE band (the YE and HE bands are similar; see Table 7 for overall statistics). Bottom: NISP RGE band S/N. The more convoluted isocontours of the IE S/N are caused by the stronger dust extinction versus the NIR. Greyed regions outside the RoI are excluded due to extreme reddening, high stellar density, and/or high zodiacal background (cf. RoI definition and Fig. 17).

In the text
thumbnail Fig. 20.

Best S/N areas in each Galactic cap. The best 1300 deg2 in each Galactic cap is shown in yellow, and the best 2600 deg2 in red (including the yellow area). These are similar to the areas planned for the first and second data releases, 2500 deg2 and 5000 deg2, respectively (cf. Fig. 46).

In the text
thumbnail Fig. 21.

EuclidIE and JE band 1σ asinh limiting surface brightness at the 10″ × 10″ scale (low surface brightness performance) and the nearby Universe: Local Group galaxies are represented by large cyan dots and galaxies up to z = 0.03 by red dots, while the purple ellipses indicate the actual sky area covered by the nearest clusters of galaxies (Virgo, Coma, Fornax, and Hydra are within the RoI). The RoI exhibits a 0.4 mag range in sensitivity for the four Euclid bands from the best area down to the ecliptic plane.

In the text
thumbnail Fig. 22.

E(B − V) over the RoI at the native 5′ resolution of the Planck/IRAS map, and surface brightness of the DGL (cirrus) in the VIS and J bands. The two islands (smaller parts of the RoI) collect the worst areas of the Euclid sky due to the combined presence of cirrus and the proximity of the Galactic plane. The small green ellipses show the three EDFs.

In the text
thumbnail Fig. 23.

Expected average number density of IE galaxies that satisfy the WL requirements (we notice that the local impact of foreground galaxies or clusters is not considered here). Upper panel: counts as a function of the local background (in mag arcsec−2) for different levels of E(B − V) reddening. Lower panel: counts as a function of reddening for different levels of the background. Lines are linear fits to the simulations.

In the text
thumbnail Fig. 24.

Expected counts of WL galaxies in arcmin−2 (see text; we note that the local impact of foreground galaxies is not considered here) in ecliptic coordinates with the Galactic plane removed (black). The decrease in the number counts is evident going towards the ecliptic plane, where the zodiacal background is higher. See also Scaramella et al. (2015).

In the text
thumbnail Fig. 25.

Expected counts of GC galaxies (with reliable spectroscopic redshifts) in deg−2 (see text) in ecliptic coordinates. Here the number counts decrease faster with latitude than for the WL counts (cf. Fig. 24) since the impact of stars (spectra) from our Galaxy on the GC one-dimensional spectra is more severe than the impact of stars (point-like) on the two-dimensional WL data.

In the text
thumbnail Fig. 26.

Weight (defined as the GC expected number of good redshifts per deg2) of observed fields along a typical EWS progression (see Sect. 7). The moving median of the weights in the previous year is superimposed (red line, starting on year 1). The approximately regular cusps every 6 months are due to the paucity of observable fields close to the Galactic and ecliptic planes and their increased background and extinction. The decrease in the median and the upper envelope in the second half of the mission is due to a relatively low quality of the sky areas still available for observations in those periods. In fact, the best areas are observed in the first half of the mission (see Figs. 20 and 46).

In the text
thumbnail Fig. 27.

Simulated spectrum of a typical galaxy (in terms of mass and star formation rate) from the Euclid weak lensing population at z = 1 (D. Masters, priv. comm.), overlapping the IE, YE, JE and HE bands. The Vera C. Rubin Observatory photometric bands (u, g, r, i, z, y) are also shown, highlighting the importance of the ground-based data to the photometric redshift derivation as the 4000 Å break falls within the i band. The total sensitivity for the ground is given at 1.2 airmass.

In the text
thumbnail Fig. 28.

Current plan (2021) of the ground-based coverage of the Euclid sky based on the six most powerful ground-based wide-field imagers and the Euclid RoI (in solid yellow). Vera C. Rubin Observatory (LSST and suggested northern and southern extensions: purple, green, and pink rectangles) and Blanco (DES, red) in the south. Subaru (WISHES), JST (JEDIS-g), Pan-STARRS, and CFHT (CFIS) in the north (orange). The best S/N 2600 deg2 areas for Euclid in each Galactic cap (black dotted contours) are key to steering the ground efforts with respect to the mission’s first data release (DR1). The Galactic referential is shown in light green.

In the text
thumbnail Fig. 29.

Example of the full schedule, including the stage-1 schedule discussed in Sect. 6.3.3. Time goes from left to right for each row, which spans one year of the survey. Subsequent years go from top to bottom. See text for a description.

In the text
thumbnail Fig. 30.

Example EWS scheduling with a rigid tiling and requiring |AA|< 1°, covering parts of the northern ecliptic hemisphere. Different colours indicate different observing epochs. Left panel: empty dark regions caused by interruptions of the EWS. Right panel: empty regions only partially recovered at a later stage, and new gaps formed at lower latitudes.

In the text
thumbnail Fig. 31.

Steps and concepts required for the computation of the EWS, linking the RoI and stage-1 schedule to patches (defined later).

In the text
thumbnail Fig. 32.

Example of a few rows of a tile-driven tessellation drawn in a cylindrical projection. Individual tiles are indicated by the red rectangles, whereas the (distorted) footprints of the FoV for non-zero latitudes are shown in blue (sizes are exaggerated for clarity). The number of tiles per row varies with latitude.

In the text
thumbnail Fig. 33.

Subset of the tessellation observable at a given epoch, with the Sun pointing at 0 deg and the transit meridians at 90 deg (trailing side) and 270 deg (leading side). This shows that the need to observe fields aligned with meridians strongly restricts visibility (see Fig. 5).

In the text
thumbnail Fig. 34.

Example of adjacent FoV slew ability at high latitudes (left panel) and low latitudes (right panel). The plots show the slew from a central tile (blue) to its neighbours: up to 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $ in light red, and from 1 . ° 2 $ 1{{\overset{\circ}{.}}}2 $ to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $ in dark red. At high latitudes, slew movements between FoVs of a tessellation are severely restricted.

In the text
thumbnail Fig. 35.

Tessellated RoI (in green), covering the four quadrants. The dashed red line represents the RoI.

In the text
thumbnail Fig. 36.

Zoomed-in view of the tessellation for the EWS at low latitudes (left panel) and high latitudes (right panel). The increased overlap of neighbouring tiles at high latitudes is evident. Orange tiles contain a blinding star (black dots) and are excluded from the observation.

In the text
thumbnail Fig. 37.

Coverage of the northern and southern ecliptic polar caps (left and right panels, respectively), with the design of polar caps drawn in red and patches drawn in blue. Extra patches below 79° are also added to polar caps (two in the north and one in the south) to observe corners of the RoI that would be difficult to observe with the regular wide strategy. The presence of the Galactic plane affects both caps. In addition, the southern boundary avoids the LMC. Blinding stars are also shown; in particular, the small gap in the southern cap is due to the presence of the extremely bright star R-Doradus.

In the text
thumbnail Fig. 38.

Example of four patch sources (in red) generated by a survey window. We notice how the patch source in the northern island intersects that quadrant only partially.

In the text
thumbnail Fig. 39.

Example of a patch source reduced by previously scheduled patches (plot of region IV only, no star skipping for clarity). Previously scheduled patches are shown in blue, with the red region depicting the patch source of some windows. We notice how the bottom side of this quadrant acquired a lat-long shape after the stacking of some patches.

In the text
thumbnail Fig. 40.

Example of sibling patch segments from a fragmented patch source (of Fig. 39). New patch segments are displayed in shades of light blue, while previously observed patches are in dark blue. The discontinuities in longitude of the previously observed patches induce the fragmentation of this patch source. Tiles containing a blinding star are excluded from the observation, here displayed in orange. Tiles in the vicinity of an extremely bright star are also excluded (even if not containing it), as is the case of the star marked with a circle at ∼(323, −36).

In the text
thumbnail Fig. 41.

Example output of the look-ahead algorithm. The sequence starts at mid-height, going downwards first. It traverses all fields running up and down, from right to left. We notice how it sometimes shows an irregular behaviour, zigzagging to the sides instead of smoothly following the expected up and down directions.

In the text
thumbnail Fig. 42.

Example run of the diffusion algorithm. The green and blue squares represent the fields in a patch. The blue ones mark the base row. The red line segments depict the progress of the algorithm as it traverses all the fields. Panel a: diffusion after fours iterations, extending the threads upwards. The remaining panels illustrate several points of the diffusion process: in panel b two threads skip across the hole left by a bright star, while two other threads merge, in panel c two threads merge while the remaining ones reach the top, in panel d the threads are extended two rows downwards, and panel e shows the final sequence, after the linking of adjacent threads. The plots are in elliptic coordinates.

In the text
thumbnail Fig. 43.

Variation in AA when observing the same EWS patch with the two different algorithms: look-ahead (left panel) and diffusion (right panel). The latter has a comparatively more monotonic progression with much smaller field-to-field AA variations and spikes. Bottom panels: are the time derivative of the top panels.

In the text
thumbnail Fig. 44.

Example of three sibling patch segments (cf. Fig. 40). Patch segments are displayed in alternating colours. Previously observed tiles are displayed in dark blue. Red arrows represent large slews (moving from patch segment to patch segment).

In the text
thumbnail Fig. 45.

Example of the ECTile output, showing the buildup of the coverage of the EWS in RSD_2021A (cf. Sect. 8) over time in ecliptic coordinates. In the main panel the EWS patches are coloured according to different years, in correspondence to the time bars below. The part of the RoI that is not observed is in light green (we recall that the RoI covers an area larger than the one needed to be observed). Under the main panel, the progression over the six years is unfolded, with two rows for each year indicating the trailing and leading directions of the telescope. Different colours represent the progression in time, from deep blue to light blue, whereas the thin line segments represent time intervals reserved for calibrations and deep fields. The light grey boxes for the leading or trailing directions mark when the telescope is actually observing in the opposite direction (for each grey box there is a corresponding coloured box). Towards the end of the mission, some of the survey windows are completely or partially light red, showing that that particular time slot has run out of sky areas to observe.

In the text
thumbnail Fig. 46.

RSD 2021A (14 514 deg2) chronology shown in celestial coordinates. RoI boundaries are shown as solid red lines. Blinding stars cause 809 avoidance areas within the reference survey, with an average of 0.785 deg2 per avoidance area, totalling 635 deg2. Dashed lines (1300 deg2 in white per Galactic cap) delimit the highest S/N areas. The ecliptic referential is over-plotted in red. The three EDFs (bright green) and the six EAFs (red diamonds, not to scale) are shown.

In the text
thumbnail Fig. 47.

Distributions of angles and slews. Left panel: AA angle. Central panel: SAA angle. Right panel: field slews up to 3 . ° 6 $ 3{{\overset{\circ}{.}}}6 $. All quantities comply with the allowed ranges.

In the text
thumbnail Fig. 48.

Area of the RoI as a function of (folded) longitude (blue). The equivalent area is visible as a function of (folded) longitude (red). The area is the available RoI area at transit per degree of longitude. The red curve assumes that the EWS observations are done at transit. Unallocated time arises when the red curve is above the blue curve. See text for details.

In the text
thumbnail Fig. 49.

Growth of the area surveyed by the RSD_2021A EWS with time (upper curve), and growth of unallocated time with time (lower curve), here in units of EWS equivalent area (one month corresponds to about 300 deg2).

In the text
thumbnail Fig. 50.

Breakdown of the total time of the Euclid mission (2190 days) in EWS and non-EWS usage, for RSD_2021A, corresponding to the schedule depicted in Fig. 29.

In the text
thumbnail Fig. 51.

Evolution of the time available for EWS observations in the RSDs produced from 2015 to 2021. The changes are due to the increase or decrease (denoted by the symbols + or −; larger changes have double symbol) in calibrations, RoI area (directly impacting the unallocated time), and ROS time. See text for details.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.