The JWST FRESCO Survey: Legacy NIRCam/Grism Spectroscopy and Imaging in the two GOODS Fields

We present the JWST Cycle 1 53.8hr medium program FRESCO, short for"First Reionization Epoch Spectroscopically Complete Observations". FRESCO covers 62 arcmin$^2$ in each of the two GOODS/CANDELS fields for a total area of 124 arcmin$^2$ exploiting JWST's powerful new grism spectroscopic capabilities at near-infrared wavelengths. By obtaining ~2 hr deep NIRCam/grism observations with the F444W filter, FRESCO yields unprecedented spectra at R~1600 covering 3.8 to 5.0 $\mu$m for most galaxies in the NIRCam field-of-view. This setup enables emission line measurements over most of cosmic history, from strong PAH lines at z~0.2-0.5, to Pa$\alpha$ and Pa$\beta$ at z~1-3, HeI and [SIII] at z~2.5-4.5, H$\alpha$ and [NII] at z~5-6.5, up to [OIII] and H$\beta$ for z~7-9 galaxies, and possibly even [OII] at z~10-12. FRESCO's grism observations provide total line fluxes for accurately estimating galaxy stellar masses and calibrating slit-loss corrections of NIRSpec/MSA spectra in the same field. Additionally, FRESCO results in a mosaic of F182M, F210M, and F444W imaging in the same fields to a depth of ~28.2 mag (5 $\sigma$ in 0.32"diameter apertures). Together with this publication, the v1 imaging mosaics are released as high-level science products via MAST. Here, we describe the overall survey design and the key science goals that can be addressed with FRESCO. We also highlight several, early science results, including: spectroscopic redshifts of Lyman break galaxies that were identified almost 20 years ago, the discovery of broad-line active galactic nuclei at z>4, and resolved Pa$\alpha$ maps of galaxies at z~1.4. These results demonstrate the enormous power for serendipitous discovery of NIRCam/grism observations.


INTRODUCTION
Revealing the dramatic build-up of galaxies from  > 6 to the peak of star formation at z∼2-3 is one of astronomy's great achievements with the Hubble and Spitzer Space Telescopes.Data from these observatories revealed that the first 1 Gyr of cosmic history ( > 6) was a time of rapid change: soon after the birth of the first stars from metal-free primordial gas, the galaxy population grew rapidly both in star-formation and stellar mass (e.g., Finkelstein 2016;Stark 2016;Oesch et al. 2018;Bouwens et al. 2023).After  ∼ 8, the star-formation and stellar mass density grew ∼10-30× up to the peak of cosmic star formation at  ∼ 2 − 4 (see Fig 7; e.g., Madau & Dickinson 2014;Bouwens et al. 2015a;Song et al. 2016;Davidzon et al. 2017;Furtak et al. 2021;Stefanon et al. 2021).
Even though this general picture of early galaxy build-up is well established, the foundation it stands on remains uncertain.In particular, HST-based analyses only probe the rest-frame ultra-violet ★ E-mail: pascal.oesch@unige.ch† Hubble Fellow (UV) of  > 3 galaxies, which can lead to a bias against dusty or old galaxies (e.g., Casey et al. 2014Casey et al. , 2018;;Wang et al. 2019;Xiao et al. 2023;Barrufet et al. 2023;Rodighiero et al. 2023).Additionally, broad-band rest-frame optical imaging can be contaminated by extremely strong emission lines that appear to be common at  ≥ 4, complicating accurate measurements of the underlying continuum and inferences of stellar masses and ages (e.g., Schaerer & de Barros 2009;Labbé et al. 2013;Stark et al. 2013;Bisigello et al. 2019;Endsley et al. 2023;Stefanon et al. 2022).
The most significant shortcoming before the advent of JWST, however, remained the lack of spectroscopically confirmed redshifts.In particular, studies exploring the earliest galaxies in the Epoch of Reionization (EoR) with spectroscopy were sparse before JWST due to disappearing Ly lines in the neutral IGM (e.g., Stark et al. 2010;Treu et al. 2013;Pentericci et al. 2014), and due to low Ly fractions in general.Consequently, at  > 7, only a handful of luminous galaxies had spectroscopic confirmations, based on Ly or weak rest-frame UV emission (e.g., Oesch et al. 2015;Stark et al. 2017), or from ALMA ISM lines (e.g., Inoue et al. 2016;Hashimoto et al. 2018;Bakx et al. 2020;Schouws et al. 2022;Bouwens et al. 2022).FRESCO was designed to cover the majority of the two CANDELS/Deep regions in the GOODS-S and -N fields, which are among the most valuable extragalactic legacy fields in the full sky.Among others, the FRESCO fields overlap with the medium-band survey JEMS (orange) that covers the original HUDF/XDF (black).The vast amount of ancillary data available in these fields made these fields the obvious choice for spectroscopic follow-up.At the highest redshifts, more than 40% of all known  ∼ 7 − 8 candidate galaxies lie in these fields, with GOODS-N featuring an overdensity of  ∼ 7 − 8 candidates.Additional overdensities are known to exist at lower redshifts, which will be mapped out with these data.The FRESCO mosaics are obtained in a 4×2 mosaic with partial overlap for optimal wavelength coverage.The total area amounts to 7.3×8.5 arcmin 2 both with grism spectroscopy and medium band imaging, which further enhances the legacy value of these fields.
This caused extra uncertainty given that the inferences on galaxy abundances at high redshift are statistical in nature.The agreement between candidate catalogs from different teams often remained limited, even when based on the same data: <60% overall, dropping to ∼ 30% at the faint end (e.g., compare catalogs from Bouwens et al. 2015a;Finkelstein et al. 2015).Similarly, blind spectroscopic surveys, e.g., with VIMOS or MUSE, up to  ∼ 6 consistently showed that a non-negligible fraction (up to 60%) of genuine high-redshift galaxies can be missed in traditional Lyman Break color-color selections due to scatter in their observed colors (e.g., Le Fèvre et al. 2015;Inami et al. 2017).
In this paper, we present the JWST FRESCO survey, short for "First Reionization Epoch Spectroscopically Complete Observations" (GO-1895; PI Oesch).FRESCO fully exploits the unprecedented grism spectroscopic ability of   at ∼ 4 − 5 m to obtain a complete sample of early star-forming emission line galaxies down to ∼ 2 × 10 −18 erg s −1 cm −2 across cosmic history in two ∼62 arcmin 2 F444W NIRCam/grism mosaics in the extragalactic legacy fields CANDELS/Deep in GOODS-South and North (Figure 1).FRESCO is designed to obtain a spectrum with  ∼ 1600 for every reionization-era galaxy in the field down to ∼0.2-0.5  *  , probing the strong rest-frame optical emission lines [OIII]+H at  = 6.7−9.0  and H+[NII] at  = 4.9 − 6.6, in addition to other emission lines across the full cosmic history.This includes dust-insensitive Paschen lines at  ∼ 1 − 3, as well as PAH 3.3 m lines at  < 0.5.Therefore, the FRESCO observations enable accurate measurements of early star formation and stellar mass build up; provide detailed insights into early, low-metallicity star formation; reveal small scale 3D clustering and measure the contribution of mergers to early galaxy assembly.These key science cases will be discussed in this paper, after first introducing the survey design.
This paper is organized as follows: in Section 2, we introduce the observational design of the FRESCO grism and imaging survey.In Section 3, we describe the main scientific goals of FRESCO, before we end with a summary in Section 4. Throughout this paper a standard ΛCDM cosmology is adopted with H 0 = 70 km s −1 Mpc −1 , Ω  = 0.3, and Ω Λ = 0.7.

OBSERVATIONAL DESIGN
The primary goals of FRESCO are to resolve two of the key limitations of our current understanding of the galaxy stellar mass build-up during the epoch of reionization at  > 6: (1) uncertain redshifts and outlier fractions, and (2) unknown emission line contribution to restframe optical broad-band photometry.In the following, we outline the survey design to overcome these limitations.

Slitless Grism Spectra
The core of FRESCO is NIRCam slitless grism observations, which result in a spectrum for every source in the field-of-view (see, e.g., Sun et al. 2022).FRESCO follows in the footsteps of numerous successful grism surveys conducted with HST such as 3D-HST (Brammer et al. 2012;Momcheva et al. 2016), GLASS (Treu et al. 2015), FIGS (Pirzkal et al. 2018), CLEAR (Estrada-Carpenter et al. 2019), WISP (Atek et al. 2010), and GRAPES (Pirzkal et al. 2004), but with greater sensitivity, at longer wavelength, and with spectral resolution a factor 10× higher.
Unlike traditional slit-based spectroscopy that can only target preselected sources, slitless spectra effectively provide a complete sample of emission line selected galaxies.With JWST, the NIRISS instrument provides slitless spectra up to ∼2m.The NIRCam camera (Rieke et al. 2005(Rieke et al. , 2023) ) includes a grism mode at longer wavelengths, providing spectra at relatively high resolution, R∼1600, up to 5m.In cycle 1, this mode is only used by a few programs, including the GTO/EIGER survey (Kashino et al. 2023;Matthee et al. 2023b), as well as the ERS/CEERS program (Finkelstein et al. 2022), who are all using the F356W filter.As we will show below, such grism observations at 3-5 m are extremely powerful probes of galaxy build-up across cosmic history.

Choice of Filter and Grism
Grism observations can use the same set of medium and wide filters as NIRCam imaging in its long-wavelength (LW) channel.The choice of filter determines the wavelength range of the spectra.The key redshift range that FRESCO targets is  > 6, i.e., the reionization epoch, where spectroscopic redshifts can be determined from strong restoptical emission lines such as H+[NII] and [OIII]+H.FRESCO therefore obtains 4-5m grism spectra with the F444W filter.This setup allows us to probe H+[NII] lines for galaxies at  = 4.9 − 6.6 and [OIII]+H lines from  = 6.7 − 9.0.If bright enough, the [OII] line could in principle even be detected in sources at  = 9.3 − 12.5.Hence, the F444W grism spectra result in near complete coverage of the epoch of reionization out to some of the most distant galaxies known in the Universe.
NIRCam is equipped with two different grisms in the LW channel that provide the same spectral resolution, but have a dispersion direction that is rotated by 90 degree relative to each other (GrismR and GrismC).FRESCO is using only one of them: GrismR.The advantages of taking data with both grisms are a better handle on contami-nation along the dispersion direction and unambiguous identification of sources in the direct image from which the spectra originate.Unfortunately, using GrismC results in a disproportionate increase of overheads, given that out-of-field images need to be taken at 35 ′′ offsets, which results in separate visits.As an example, if both GrismC and GrismR were used, FRESCO's survey time would have increased by 15hrs (∼30%) without any gain in sensitivity.Even with only one grism observation, most shortcomings can be mitigated, as demonstrated by the extremely successful HST grism survey 3D-HST, from which the community has built extensive experience to deal with overlapping spectra (Brammer et al. 2012).Due to these excessive overhead costs, FRESCO is obtained with GrismR observations only, similarly to the GTO/EIGER survey (GTO-1243; PI Lilly).

NIRCam/Grism Mosaic and Wavelength Coverage
The wavelength coverage of the NIRCam/grism depends on the location of the source on the detector.At a wavelengths of 3.95 m, the F444W grism and the direct image coincide.However, at other wavelengths, the positions are shifted by 9.85Å/pixel, varying only slightly across the field.This means that only a certain fraction of sources from the direct image will obtain complete spectral coverage.Additionally, the grism spectra cover sources from outside the field-of-view.The GrismR spectra of detector A and B are dispersed in the opposite direction.By obtaining a second pointing offset by ∼ 1.8 arcmin, FRESCO achieves coverage over 3.9-4.4m over almost the full field of view and partial coverage from 3.9-5 m over 4.9 arcmin in x-direction.In the y-direction, FRESCO obtained 4 pointings with minimal overlap (thus spanning 8.5 arcmin).Therefore, each FRESCO field covers 6.7×8.5 arcmin 2 (57 arcmin 2 ) at 3.9-4.4m and 4.9×8.5 arcmin 2 (42 arcmin 2 ) at 3.9-5 m (corresponding to 73% of the full mosaic shown in Figure 1).More details on the spectral coverage and the spectral extraction will be provided in Brammer et al., in prep.

Exposure Times and Sensitivity
FRESCO is designed to reach a 5 emission-line sensitivity for NIRCam/grism observations of ∼2×10 −18 erg s −1 cm −2 (for compact sources), which is the expected [OIII]5008 flux of sources at 1 mag below L ★ at  = 7 − 9 (based on the LF derived from De Barros et al. 2019; see also Matthee et al. 2023b).This is achieved with eight grism exposures of 880s taken with the MEDIUM2 readout mode with 9 groups for each pointing.This results in a total grism exposure time of 7043 s at each position (see also Table 2).Four large scale dithers are used, but no sub-pixel dithering.The grism exposures are taken in two sets of these four-point dithers, in each set using a different short-wavelength filter for imaging (see Sect. 2.2).
An example of the resulting FRESCO imaging and spectroscopy in the GOODS-S field is shown in Figure 2.

Grism Data Reduction
Details of the NIRCam/grism data reduction will be provided in a data release paper (Brammer et al, in prep).Briefly, we use the publicly available grizli code 1 , with slightly modified sensitivity curves and spectral traces based on the v4 grism configuration files provided by Nor Pirzkal 2 .grizli is used to reduce all images and spectra and to align them to Gaia-matched reference frames.Given that we are mostly interested in emission line detections, continuum subtracted spectra are created using a running median filter along each row.Following Kashino et al. (2023), the filter uses a 12 pixel central gap to minimize self-subtraction.For each source in the imaging catalog, optimally extracted spectra are then produced based on the individual input exposures.The spectral catalogs will be released together with individual science papers and with a future data release paper.

FRESCO's Imaging Legacy
Simultaneous with the grism spectra, FRESCO obtains NIRCam medium band images in two short-wavelength filters (F182M and F210M) to extend the space-based legacy data in these fields to 1.8 and 2.1 m down to 28.3 mag (at 5).These medium-band filters seamlessly extend the NIR wavelength coverage of the HST data in the CANDELS/Deep area, complementing the GTO/JADES NIR-Cam imaging.These data improve the UV spectral slope measurement of  > 7 galaxies and, at lower redshift, allow the community to push rest-frame optical imaging analyses beyond the peak of cosmic star formation, from  ∼ 2.5 to beyond  ∼ 4.After the grism exposures, three direct and out-of-field images are taken in the F444W filter.These images are needed to associate the spectra with individual galaxies.The exposure times of these direct images are set to produce a > 4 detection of every single galaxy for which a significant emission-line detection is expected.This ensures that spectra can be associated with the imaged galaxies in the field.FRESCO thus obtained 934s exposures in F444W, split over three offset positions, leading to a 5-depth of 28.2 mag (see Table 2).
The medium-band images are taken with the same readout mode as the grism exposures (MEDIUM2, 9 groups), resulting in an exposure 1 https://github.com/gbrammer/grizli 2 https://github.com/npirzkal/GRISMCONFtime of 3522s each.The F182M filter is used during the three F444W imaging exposures, which are taken with the SHALLOW4 readout pattern with 6 groups.It thus receives an additional integration of 934s.The total exposure times per pointing for each filter are listed in Table 2.

Target Fields
FRESCO is split into two separate fields with an equal setup covering two of the most well-studied fields in the sky: the central regions of the two GOODS fields (North and South; Giavalisco et al. 1996).These were also covered by the Deep tier of the CANDELS survey (Koekemoer et al. 2011;Grogin et al. 2011).In GOODS-South, the FRESCO field further covers the HUDF/XDF field (Beckwith et al. 2006;Illingworth et al. 2013;Ellis et al. 2013).Additional imaging over these fields was taken over the years by a very large number of programs.A complete listing of HST programs can be found on the Hubble Legacy Field (HLF) release page3 (see also Whitaker et al. 2019 andIllingworth et al. 2016).
These fields are thus among the most valuable extragalactic legacy fields, especially for distant galaxy science.More than 40% of all HST-selected  ∼ 7 − 8 candidate galaxies from blank-field surveys lie in these two areas alone (e.g.Bouwens et al. 2015a;Finkelstein et al. 2015).One of the first sources spectroscopically confirmed with Ly at  > 7 lies in GOODS-N ( = 7.5, Finkelstein et al. 2013) within a spectroscopically confirmed overdensity of  ∼ 7.5 galaxies (Jung et al. 2020).GOODS-N also contains the most distant confirmed galaxy with HST: GN-z11 (Oesch et al. 2016  The two galaxies shown here have already been detected with the NICMOS camera in Bouwens et al. (2004).The left stamps correspond to the F444W direct image, while the right panels show the continuum-subtracted grism data.The [OIII] doublet is significantly detected in both cases resulting in an unambiguous redshift confirmation.The H line is found to be weaker for both sources.This example demonstrates the power of NIRCam/grism observations to obtain an emission line selected census of galaxies in the epoch of reionization.et al. 2023), a dusty, red QSO at  ∼ 7 (Fujimoto et al. 2022), as well as an enigmatic dusty source, HDF850.1,embedded in another overdensity at  ∼ 5.2 (Walter et al. 2012).By targeting both GOODS fields, FRESCO thus enables a first estimate of the diversity of cosmic structures across cosmis history, well into the EoR.
Both of the FRESCO fields are (partially) observed with the proprietary NIRCam/GTO data from the JADES survey (Eisenstein et al. 2023; see also Helton et al. 2023), and with two of the MIRI GTO programs (Rieke et al, in prep.; Ostlin et al, in prep.; see also Rieke et al. 2015;Rinaldi et al. 2023).Additionally, the FRESCO data overlap with the public programs JEMS (Williams et al. 2023) and NGDEEP (Bagley et al. 2023).

Data Acquisition
The FRESCO survey design results in a total science time of 35.5 hrs, for a total of 53.8 hrs including overheads.The FRESCO-S data in the GOODS-South field were acquired between Nov 13 and Nov 18, 2022, at an orientation of ∼0 deg (V3PA in the range 353.9 to 0.6).The FRESCO-N data in the GOODS-North field were obtained between Feb 11 and 13, 2023, at an orientation of 230.5 (V3PA in the range 230.37-230.59).The central position of FRESCO-N was offset slightly from the center of the CANDELS/Deep region in order to capture a few high-profile objects and an expected overdensity of  ∼ 7.5 galaxies, as discussed in the previous section.

SCIENTIFIC GOALS
FRESCO exploits the unique capability of the NIRCam/grism to obtain deep spectra at ∼4-5m for the entire survey field.This enables a wealth of both targeted science, as well as serendipitous discovery.Below, we discuss some of the most important scientific questions the community can address with these data.(Labbé et al. 2010).This enables extremely efficient spectroscopic confirmation (inset).However, such strong lines contaminate broad-band photometry, as clearly seen in the IRAC CH2 photometry of this source.This can result in stellar mass uncertainties of up to 5 − 10×.FRESCO's spectra enable us to correct JWST's broad-band photometry to obtain true rest-frame optical continuum and hence stellar mass measurements.

Spectroscopic Redshifts in the Heart of Reionization
The first 1 Gyr of cosmic history that constitute the cosmic reionization epoch remain a key unknown in our understanding of the Universe's evolution (see recent reviews Dayal & Ferrara 2018;Robertson 2022).Despite enormous efforts, the galaxies in the EoR have eluded almost all attempts at spectroscopic characterization before JWST, apart from a handful of especially bright (or lensed) galaxies.Only two dozen "normal" galaxies had confirmed spectroscopic redshifts at  > 7, out of more than 2000 candidates that were known from prime extragalactic HST legacy fields alone (Fig. 2; e.g., Bouwens et al. 2015b;Finkelstein et al. 2015;Atek et al. 2018).The main reason for this difficulty was that the primary line for redshift confirmations, Ly, is severely attenuated at  > 6 due to absorption in the largely neutral intergalactic medium (e.g., Dĳkstra 2014;Mason et al. 2018).With JWST, this critical shortcoming can finally be overcome through rest-frame optical spectroscopy.FRESCO probes the [OIII]+H emission lines of galaxies at  = 6.7 − 9.0 as well as H for sources at  = 4.9−6.7.The spectroscopically confirmed redshifts allow us to probe the small scale clustering and 3D correlation function of galaxies during the reionization epoch (e.g., Endsley et al. 2020).From simulations, the progress of reionization is expected to be correlated with the underlying density field (see Fig. 3; Qin et al. 2022;Leonova et al. 2022).
Figure 4 demonstrates the power of FRESCO's NIRCam/grism spectra to measure the redshifts of EoR galaxies.The two highlighted sources are Lyman break galaxies (LBGs) that were identified almost 20 years ago by Bouwens et al. (2004) from the HST/NICMOS images over the Hubble Ultra Deep Field (Thompson et al. 2005).With FRESCO, their high-redshift nature is confirmed through the detection of the [OIII]4960,5008 doublet at  grism = 7.238 and  grism = 7.223, respectively.Importantly, these objects were not selected prior to the observations: all objects in the field are observed spectroscopically, without the need to pre-select targets.This is the main advantage of grism observations and enables a very wide array of science.

11
>2700 km s −1 Figure 6.An example of a serendipitously discovered broad-line, low-mass AGN from FRESCO spectra.The extremely broad H+[NII] lines identify these galaxies as hosting massive black holes (for more information, see Matthee et al. 2023a).NIRCam/grism observations are uniquely capable of discovering such sources and confirming their number densities in the future.Given the high spectral resolution of  ∼ 1600, such observations also have the power to reveal kinematic information for less extreme galaxies.

Accurate Stellar Masses in the Heart of Reionization
From HST+Spitzer observations, it has been known for several years that the emission line strengths of [OIII] or H are rapidly increasing toward higher redshift (e.g.Schaerer & de Barros 2010;Fumagalli et al. 2012;Labbé et al. 2013), resulting in observed-frame equivalent widths that are >5000 Å on average (e.g., Smit et al. 2014;Roberts-Borsani et al. 2016;De Barros et al. 2019;Endsley et al. 2021).This evolution has dramatic consequences for stellar mass estimates that can remain uncertain by factors up to 5-10× (Fig 5 ) a problem that still applies to NIRCam imaging.Bisigello et al. 2019 estimate median correction factors up to 0.87dex, if strong emission lines are not accounted for.FRESCO provides emission line measurements for all sources in the field, enabling precise corrections of the broad-band fluxes for emission line contamination on a sourceby-source basis.A prominent example is shown in Figure 5.The source in question still has one of the most distant Ly detections at  Ly = 7.5 (Finkelstein et al. 2013;Jung et al. 2020).Thanks to this, it was clear that the excess in the IRAC 4.5m band was due to strong rest-frame optical emission lines.Nevertheless, the strength of these lines was uncertain.The FRESCO data shows that the source is actually composed of two components with slightly offset velocities.For both components, the contribution of these emission lines to the broad-band photometry can now be assessed separately.With these corrections, FRESCO thus provides accurate stellar mass measurements over two key deep fields.Additionally, it will enable the derivation of statistical correction factors for other imaging surveys.
While these strong rest-optical emission lines can complicate stellar mass estimates, they enable efficient spectroscopic confirmations.Therefore, FRESCO's spectra continue to propel the spectroscopic frontier into the heart of cosmic reionization by obtaining redshifts down to 0.2-0.5  * at  ∼ 7 − 9 (see Fig 2).This critically enables measurements of the UV luminosity and stellar mass functions based on pure spectroscopic samples at  ≳ 5.

Metal-Free Star Formation and the Build-up of Metals at 𝑧 > 3
A major question for extragalactic surveys with JWST remains whether it is possible to detect primordial (PopIII/zero-metallicity) star formation.While candidate PopIII galaxies at  ∼ 7 have been claimed in the literature (e.g., Sobral et al. 2015), none have been confirmed (see, e.g., Bowler et al. 2017;Matthee et al. 2017).FRESCO has a unique discovery potential for extremely low metallicity candidates at  > 6.The medium-band imaging at 1.8-2 m will result in improved UV continuum slope measurements at  ∼ 7 − 9 enabling the identification of especially blue galaxies for which the grism [OIII] and H line ratios provide an initial gas-phase metallicity estimate ([OIII]/H decreases at < 0.2 ⊙ ; e.g., Maiolino et al. 2008;Inoue 2011;Curti et al. 2020).FRESCO thus has the capability to provide promising, extremely metal-poor candidates for future NIRSpec follow-up.At  ∼ 4.9 − 6.6, FRESCO continues to trace the build-up of metals through [NII]/H line ratios in order to constrain the mass-metallicity relations and models of early chemical enrichment.As such, FRESCO data allows us to trace the transition from extremely metal poor star formation at  > 6 to more enriched conditions that are found at later times.

The Prevalence of AGN among Early Galaxies
Another central question of extragalactic astronomy is the abundance of active galactic nuclei (AGN) in the early Universe.These sources could provide a non-negligible contribution to reionization (e.g., Madau & Haardt 2015;Giallongo et al. 2015), and are expected to explain the abundances of super-massive black holes at later times (e.g., Volonteri 2010).Interestingly, several of the most luminous  ∼ 7 − 9 galaxies show indications for AGN activity, based on their ground-based rest-UV spectra, suggesting that the AGN fraction could be significant (Laporte et al. 2017;Sobral et al. 2018).Owing to the R∼1600 spectra, FRESCO can test this scenario directly by enabling the identification of potential AGN through broad emission lines across the full redshift range  ∼ 5 − 9 based on H, H, or [OIII] as well as through indirect methods such as the mass excitation diagram (after a recalibration to  ∼ 7 − 9 based on   spectra; Juneau et al. 2011;Trump et al. 2013).
First JWST spectra have already revealed a significant population of low-mass AGN (e.g., Kocevski et al. 2023;Übler et al. 2023;Harikane et al. 2023), and more candidates are identified through NIRCam imaging (e.g., Labbe et al. 2023).NIRCam/grism observations have the potential to provide a complete sample of such sources.A prominent example from the FRESCO data is shown in Figure 6, with an extremely broad H emission line at  ∼ 5.This source is part of a larger sample of such galaxies identified in NIRCam/grism data (Matthee et al. 2023a).

Unbiased Star-Formation Rate Indicators:
Star formation rate (SFR) is a fundamental observable property of galaxies that is required to trace the growth and formation of galaxies throughout cosmic time.Gold standard SFR indicators are Hydrogen recombination lines, as they emerge from Hii regions around the most massive and recently born stars and, hence, trace the almostinstantaneous star formation rate in galaxies (Kennicutt & Evans 2012).Aside from Ly, which suffers from uncertainties associated to resonant scattering, Balmer optical lines, such as H and H are the bright H lines and good tracers of star formation activity.However, these lines are highly affected by dust attenuation, and the uncertainties associated with their attenuation correction factors (e.g., uncertainties in the nebular reddening measurements and attenuation curve assumptions) hamper their potential as accurate SFR diagnostics (e.g., Fanelli et al. 1988;Reddy et al. 2015;Shivaei et al. 2018, among many more).Pa, on the other hand, is an instantaneous and dust-insensitive SFR indicator, owing to its longer wavelength, making it an important SFR diagnostic, particularly in dusty starforming galaxies at the peak epoch of cosmic star formation history.FRESCO takes advantage of the unprecedented near-IR capabilities of JWST to observe Paschen lines (Pa and Pa) and trace opticallythick star formation in large samples of galaxies at  ∼ 1 − 3, which is typically missed in optical surveys (see also Finkelstein et al. 2011;Cleri et al. 2022;Reddy et al. 2023).Another tracer of obscured star formation activity that is accessible to FRESCO is the 3.3 m feature at  ≲ 0.5, which is the emission from PAH dust grains (see, e.g., Genzel & Cesarsky 2000;Kim et al. 2012).
The FRESCO spectra at ∼ 4 − 5m thus have the capability to trace star-formation across the entire cosmic history (see Fig. 7).To illustrate this, Figure 8 shows a collage of extracted emission line spectra in the FRESCO-S field.This includes 333 sources that previously had spectroscopic redshifts measured from the literature up to  ∼ 6, in addition to 27 sources for which FRESCO detected bright [OIII]+H lines at  ∼ 7. The spectra are sorted by redshift.Given the availability of previous spectroscopic redshifts the Pa lines at  ∼ 1−2 dominate the figure.Several overdensities are clearly visible.The figure shows the power of NIRCam/grism observations to probe star-formation and clustering across cosmic history.

Spatial Distribution of Star Formation at Cosmic Noon
A unique feature of grism observations is that they result in spatially resolved maps of emission lines and hence star formation.For Figure 10.FRESCO's imaging and spectroscopy has the power to reveal the nature of the enigmatic population of optically-faint galaxies.Shown is an example galaxy from Xiao et al. (in prep.) in the FRESCO-S field.The source was known from ALMA, but its counterpart could not be identified in the low-spatial resolution IRAC imaging.FRESCO's F444W images clearly reveal this galaxy, and with the NIRCam/grism spectra its redshift can be determined through the H emission line.
instance, H maps from HST grism observations of  ∼ 1 − 2 galaxies have been used to reveal star-forming clumps and centrally suppressed star formation in massive galaxies, consistent with inside-out growth (e.g.Wuyts et al. 2013;Nelson et al. 2016).This appears to hold, even when accounting for dust gradients inferred from resolved broadband colors or stacked Balmer decrements (e.g.Tacchella et al. 2018).However, ALMA observations of the 870m dust continuum emission show a different picture: the dust continuum in many massive galaxies is centrally concentrated (Tadaki et al. 2020).These observations support a different scenario, in which violent dissipative events fuel nuclear star formation.Given the uncertainties in H dust corrections and complications in the interpretation of the dust continuum emission (e.g., due to different heating mechanisms, dust temperature gradients, insensitivity to diffuse emission), an instantaneous and dust-insensitive tracer of star formation is critically needed.FRESCO provides exactly this, with spatially resolved measurements of the Paschen series lines at  ∼ 1 − 3.These measurements allow us to finally answer the question: do massive galaxies at cosmic noon form their last stars in their centers or in their outskirts before they quench?
A first example of a resolved Pa map in the FRESCO-S field is shown in Figure 9.The grizli-generated Pa map indeed reveals a very compact star-forming core for this galaxy that is completely absent in the shorter wavelength imaging.This is embedded in more extended star-formation throughout the galaxy.NIRCam/grism observations such as from FRESCO enable systematic studies of the spatial distribution of star-formation in cosmic noon galaxies.

Illuminating the Complex ISM in Assembling Galaxies
At  > 5, galaxies are rapidly assembling, resulting in highly complex interstellar media, as demonstrated by ALMA observations of EoR galaxies, where spatial offsets between the rest-UV continuum, and different ISM lines such as [CII] or [OIII]88m are common (e.g., Maiolino et al. 2015;Carniani et al. 2017).FRESCO provides detailed H (H) maps of a large sample of galaxies at  = 4.9 − 6.6 ( = 6.7 − 9.0) to constrain how and where within galaxies stars are forming during the epoch of reionization (for first examples see Good arguments for Medium Band Filters? -These two medium band filters would provide accurate sampling of the Balm identify quiescent galaxies at these redshifts. -Seamlessly extend HST studies of quiescent galaxies from z~2.5 to z~4 to th -Would enable z>10 galaxy selection based on CANDELS data.(don't expect get lucky..) -Not too much to gain for z~6-8 population, except improved UV slope measu Medium Band Filter Science?
Figure 11.FRESCO obtains images in two medium-band filters at ∼ 2m (and a direct image in F444W).These data seamlessly extend HST's legacy imaging over these fields and complement the GTO data resulting in vast improvements in UV-continuum slope measurements for  > 6 galaxies and enabling the selection of quiescent galaxies up to  ∼ 4.
the EIGER program Matthee et al. 2023b).Comparisons between H and rest-optical continuum sizes provide a first test of inside-out growth during an epoch when galaxies are still assembling -rather than fully settled in -(cold) disks, or may reveal homologous growth with a similar pace of stellar build-up across all radii.

Illuminating the Dark Side of Star-Formation
Over the past few years it has become more and more clear that HST-based surveys have been missing a potentially important population of massive, star-forming galaxies at  > 3 that are sufficiently obscured such that they remained undetected at rest-UV wavelengths (e.g., Caputi et al. 2015;Franco et al. 2018;Alcalde Pampliega et al. 2019;Wang et al. 2019;Fudamoto et al. 2021;Xiao et al. 2023).These galaxies have been identified either through Spitzer/IRAC or through ALMA observations.However, given the limited photometry available, their photometric redshifts and thus also their nature remained highly uncertain.With JWST we can now finally probe the rest-frame optical wavelengths up to  ∼ 9, which brings this optically-faint galaxy population into view (e.g., Barrufet et al. 2023;Rodighiero et al. 2023;Pérez-González et al. 2023).FRESCO's imaging and spectroscopy is especially powerful to characterize this enigmatic population of galaxies.The grism spectra enable -for the first time -to obtain spectroscopic redshifts for such galaxies thanks to emission line searches.An example of this is shown in Fig 10 of a galaxy that has previously been detected only at sub-millimeter wavelengths (Cowie et al. 2018;Yamaguchi et al. 2019;Gómez-Guĳarro et al. 2022;Xiao et al. 2023, Xiao et al. in prep).

The Emergence of Quiescent Galaxies
At  < 2, the massive galaxy population is dominated by dead, quiescent galaxies, with no significant ongoing star formation (e.g., Muzzin et al. 2013;Davidzon et al. 2017).However, it is still unknown when such quiescent galaxies first appeared in the Universe.Deep ground-based K-band surveys demonstrated that quiescent galaxies existed since  ∼ 4, however, their exact number densities remain debated (Straatman et al. 2015;Davidzon et al. 2017;Forrest et al. 2020;Valentino et al. 2020), mainly because HST-based surveys are limited to rest-frame optical observations at  < 2.5.First searches with JWST are already finding and confirming quiescent, low starformation rate galaxies up to  ∼ 5 − 7 (see, e.g., Valentino et al. 2023;Carnall et al. 2023;Looser et al. 2023).A large population of such sources in the very early Universe challenges simulations (e.g., Wellons et al. 2015;Merlin et al. 2019;Hartley et al. 2023), such that the timing of when massive galaxies grow, shut off their star formation and turn quiescent is a very powerful constraint on galaxy evolution models.The short-wavelength imaging from FRESCO is designed to complement the JADES/GTO broad-band data in these fields and seamlessly extends HST's legacy data to enable accurate selections of quiescent galaxies up to  ∼ 4 from Balmer break measurements (see Fig. 11).In the SW channel, the imaging filters are the same as obtained with the JEMS survey over the HUDF/XDF.

Legacy Science
The CANDELS/Deep fields in the centers of GOODS-S and -N are two of the most studied areas in the sky with the most comprehensive ancillary data.FRESCO builds on these data and further increases the legacy value of these fields beyond that of the GTO teams by enabling a more complete census of emission line galaxies through blind spectroscopy.Among others, FRESCO provides emission samples for future spectroscopic follow-up with NIRSpec.Thus, FRESCO's zero-proprietary data enable a wealth of legacy science to be performed by the community.First science papers by the community using FRESCO imaging and spectroscopy have already been published as pre-prints (Laporte et al. 2022;Helton et al. 2023).

SUMMARY
This paper provided an overview of the medium program FRESCO -a NIRCam grism spectroscopy and imaging survey of two fields in the GOODS-North and -South fields, respectively.The power of NIR-Cam/grism observations are clearly demonstrated with this survey: the grism provides  ∼ 1600 spectroscopy of all sources in the field of view in its long-wavelength channel, while imaging is obtained in a short-wavelength filter at the same time.Adopting this strategy with the F444W filter, FRESCO obtained deep spectra covering 3.8-5.0m, reaching average line sensitivities of 2×10 −18 erg s −1 cm −2 (5).At the same time, deep F182M and F210M images are obtained reaching ∼ 28.2 mag rms (5 in 0. ′′ 32 diameter apertures), together with F444W direct images reaching similar depth.The FRESCO observations thus significantly enhance the rich ancillary dataset in these key legacy fields, enabling a vast amount of science by the community.The spectral coverage allows one to probe emission line galaxies across almost the full cosmic history.We have highlighted a few science cases, that showcase the enormous power for serendipitous discovery of such NIRCam/grism observations.Future, widerarea grism programs would thus be able to obtain large, complete samples of rare classes of galaxies that will be very difficult to followup, e.g., with NIRSpec spectroscopy through targeted programs.We hope that the community will use the FRESCO dataset for a large number of scientific discoveries as well as as inspiration for many further JWST programs, including NIRCam/grism observations.
Figure1.FRESCO was designed to cover the majority of the two CANDELS/Deep regions in the GOODS-S and -N fields, which are among the most valuable extragalactic legacy fields in the full sky.Among others, the FRESCO fields overlap with the medium-band survey JEMS (orange) that covers the original HUDF/XDF (black).The vast amount of ancillary data available in these fields made these fields the obvious choice for spectroscopic follow-up.At the highest redshifts, more than 40% of all known  ∼ 7 − 8 candidate galaxies lie in these fields, with GOODS-N featuring an overdensity of  ∼ 7 − 8 candidates.Additional overdensities are known to exist at lower redshifts, which will be mapped out with these data.The FRESCO mosaics are obtained in a 4×2 mosaic with partial overlap for optimal wavelength coverage.The total area amounts to 7.3×8.5 arcmin 2 both with grism spectroscopy and medium band imaging, which further enhances the legacy value of these fields.

Figure 2 .
Figure 2. Example of FRESCO's spectroscopic strategy.The left panel shows a portion of the F444W direct image in the FRESCO-S field.The middle panel shows the associated GrismR data, from which FRESCO can identify EoR galaxies as emission line sources.Down to UV continuum magnitudes of 26.8,  = 6.7 − 9.0 sources are expected to show three emission lines from the [OIII] doublet as well as H, resulting in unambiguous redshift identifications.For example, the right panels show the 2D and the extracted 1D spectra of a source in the FRESCO-S field with an unambiguous grism redshift at  = 7.6.

Figure 3 .
Figure 3.The average spectroscopic sensitivity of the FRESCO survey to observe strong rest-frame optical emission lines during the epoch of reionization.Before the advent of JWST, less than 1% of known sources at  > 6 had spectroscopically confirmed redshifts.FRESCO is obtaining complete samples with spec- down to   = −19 (z∼6) to −20 (z∼8) thanks to strong H and [OIII]+H lines, respectively.This allows us to map out the large-scale structure during cosmic reionization.Simulations predict that ionized regions first appear around overdensities as shown by a slice through the DRAGONS simulation (top panel).The 3D large scale structure during the epoch of reionization will for the first time be revealed thanks to FRESCO redshifts.

Figure 4 .
Figure 4.An example of how the NIRCam/grism data from the FRESCO survey can finally obtain spectroscopic redshift confirmation for sources that have been identified as high-redshift galaxy candidates almost 20 years prior.The two galaxies shown here have already been detected with the NICMOS camera inBouwens et al. (2004).The left stamps correspond to the F444W direct image, while the right panels show the continuum-subtracted grism data.The [OIII] doublet is significantly detected in both cases resulting in an unambiguous redshift confirmation.The H line is found to be weaker for both sources.This example demonstrates the power of NIRCam/grism observations to obtain an emission line selected census of galaxies in the epoch of reionization.

Figure 5 .
Figure 5.At  > 5, the rest-frame optical [OIII]+H lines are found to have observed-frame EWs of > 5000 Å -on average(Labbé et al. 2010).This enables extremely efficient spectroscopic confirmation (inset).However, such strong lines contaminate broad-band photometry, as clearly seen in the IRAC CH2 photometry of this source.This can result in stellar mass uncertainties of up to 5 − 10×.FRESCO's spectra enable us to correct JWST's broad-band photometry to obtain true rest-frame optical continuum and hence stellar mass measurements.

Figure 7 .
Figure 7. FRESCO is designed to return a complete sample of star-forming galaxies based on emission lines across the majority of cosmic history.FRESCO 4-5m spectra (a) directly constrain the debated accelerated evolution of the SFRD at  > 8 based on [OIII]+H samples at  =6.9-9.0;(b) follow galaxy build-up through the reionization epoch with H+[NII] measurements at  =4.9-6.6; and (c) probe the peak of cosmic SF (z∼1-3) through dust-insensitive Pa and Pa maps.In the local Universe, FRESCO even covers PAH lines at 3.3 m for a small number of galaxies (see also Fig 8).

Figure 8 .
Figure 8. Line coverage at different redshifts in the NIRCam/grism observations of FRESCO.Shown is a stack of 333 2D spectra in the GOODS-South field for sources with pre-existing spectroscopic redshifts from the literature up to  ∼ 6.5.Beyond that, we show examples of new H and [OIII] lines measured for galaxies at  > 6.8.Given the large number of spectroscopic redshifts in this field at  ∼ 1 − 2, Pa lines currently dominate this figure.Galaxy overdensities are clearly visible.Full spectral extractions of the FRESCO NIRCam/grism data will be described in Brammer et al. (in prep).

Figure 9 .
Figure 9.The spatial distribution of SFR from emission line maps across cosmic history as probed by FRESCO.NIRCam/grism data uniquely provide spatially resolved line maps at 4-5m.Shown is the Pa map from a dusty  = 1.38 galaxy in the FRESCO-S field.The left panel shows an RGB composite of the FRESCO F182M, F210M, and F444W images.The central panel shows the NIRCam/grism data, from which we extract an emission line map at the Paschen- wavelength using the grizli tool.The FRESCO dataset thus reveals the sites of star-formation in galaxies.

Table 1 .
Summary of the FRESCO fields and areas.

Table 2 .
× 10 −18 erg s −1 cm −2 Summary of exposure times per FRESCO pointing.In total, 2 mosaics of 4×2 pointings were obtained.Line sensitivity estimated for compact sources, integrated over full extension of the line.Imaging depths measured as the 5 rms in 0. ′′ 32 diameter apertures, without corrections to total fluxes.