We give an overview and describe the rationale, methods, and first results from NIRCam images of the JWST “Prime Extragalactic Areas for Reionization and Lensing Science” (PEARLS) project. PEARLS uses up to eight NIRCam filters to survey several prime extragalactic survey areas: two fields at the North Ecliptic Pole (NEP); seven gravitationally lensing clusters; two high redshift protoclusters; and the iconic backlit VV 191 galaxy system to map its dust attenuation. PEARLS also includes NIRISS spectra for one of the NEP fields and NIRSpec spectra of two high-redshift quasars. The main goal of PEARLS is to study the epoch of galaxy assembly, active galactic nucleus (AGN) growth, and First Light. Five fields—the JWST NEP Time-Domain Field (TDF), IRAC Dark Field, and three lensing clusters—will be observed in up to four epochs over a year. The cadence and sensitivity of the imaging data are ideally suited to find faint variable objects such as weak AGN, high-redshift supernovae, and cluster caustic transits. Both NEP fields have sightlines through our Galaxy, providing significant numbers of very faint brown dwarfs whose proper motions can be studied. Observations from the first spoke in the NEP TDF are public. This paper presents our first PEARLS observations, their NIRCam data reduction and analysis, our first object catalogs, the 0.9–4.5 μm galaxy counts and Integrated Galaxy Light. We assess the JWST sky brightness in 13 NIRCam filters, yielding our first constraints to diffuse light at 0.9–4.5 μm. PEARLS is designed to be of lasting benefit to the community.

We give an overview and describe the rationale, methods, and testing of the Hubble Space Telescope (HST) Archival Legacy project “SKYSURF.” SKYSURF uses HST’s unique capability as an absolute photometer to measure the ∼0.2–1.7 μm sky-surface brightness (sky-SB) from 249,861 WFPC2, ACS, and WFC3 exposures in ∼1400 independent HST fields. SKYSURF’s panchromatic data set is designed to constrain the discrete and diffuse UV to near-IR sky components: Zodiacal Light (ZL), Kuiper Belt Objects (KBOs), Diffuse Galactic Light (DGL), and the discrete plus diffuse Extragalactic Background Light (EBL). We outline SKYSURF’s methods to: (1) measure sky-SB levels between detected objects; (2) measure the discrete EBL, most of which comes from AB≃17–22 mag galaxies; and (3) estimate how much truly diffuse light may exist. Simulations of HST WFC3/IR images with known sky values and gradients, realistic cosmic ray (CR) distributions, and star plus galaxy counts were processed with nine different algorithms to measure the “Lowest Estimated Sky-SB” (LES) in each image between the discrete objects. The best algorithms recover the LES values within 0.2% when there are no image gradients, and within 0.2%–0.4% when there are 5%–10% gradients. We provide a proof of concept of our methods from the WFC3/IR F125W images, where any residual diffuse light that HST sees in excess of zodiacal model predictions does not depend on the total object flux that each image contains. This enables us to present our first SKYSURF results on diffuse light in Carleton et al.

The first James Webb Space Telescope (JWST) data on the massive colliding cluster El Gordo allow for 23 known families of multiply lensed images to be confirmed and for eight new members of these families to be identified. Based on these families, which have been confirmed spectroscopically by MUSE, we derived an initial lens model. This model guided the identification of 37 additional families of multiply lensed galaxies, among which 28 are entirely new systems, and nine were previously known. The initial lens model determined geometric redshifts for the 37 new systems. The geometric redshifts agree reasonably well with spectroscopic or photometric redshifts when those are available. The geometric redshifts enable two additional models that include all 60 families of multiply lensed galaxies spanning a redshift range 2 < z < 6. The derived dark-matter distribution confirms the double-peak configuration of mass found by earlier work with the southern and northern clumps having similar masses. We confirm that El Gordo is the most massive known cluster at z > 0.8 and has an estimated virial mass close the maximum mass allowed by standard cosmological models. The JWST images also reveal the presence of small-mass perturbers that produce small lensing distortions. The smallest of these is consistent with being a dwarf galaxy at z = 0.87 and has an estimated mass of 3.8 × 109 M⊙, making it the smallest substructure found at z > 0.5. The JWST images also show several candidate caustic-crossing events. One of them is detected at high significance at the expected position of the critical curve and is likely a red supergiant star at z = 2.1878. This would be the first red supergiant found at cosmological distances. The cluster lensing should magnify background objects at z > 6, making more of them visible than in blank fields of a similar size, but there appears to be a deficiency of such objects.

We report the results of James Webb Space Telescope/NIRCam observations of 19 (sub)millimeter sources detected by the Atacama Large Millimeter Array (ALMA). The accurate ALMA positions allowed unambiguous identifications of their NIRCam counterparts. Taking gravitational lensing into account, these represent 16 distinct galaxies in three fields and constitute the largest sample of its kind to date. The counterparts’ spectral energy distributions cover from rest-frame ultraviolet to near-IR and provide photometric redshifts (1 < z < 4.5) and stellar masses (M * > 1010.5 M ⊙), which are similar to submillimeter galaxies (SMGs) studied previously. However, our sample is fainter in (sub)millimeter than the classic SMG samples are, and our sources exhibit a wider range of properties. They have dust-embedded star formation rates as low as 10 M ⊙ yr−1, and the sources populate both the star-forming main sequence and the quiescent categories. The deep NIRCam data allow us to study the rest-frame near-IR morphologies. Excluding two multiply imaged systems and one quasar, the majority of the remaining sources are disk-like and show either little or no disturbance. This suggests that secular growth is a potential route for the assembly of high-mass disk galaxies. While a few objects have large disks, the majority have small disks (median half-mass radius of 1.6 kpc). At this time, it is unclear whether this is due to the prevalence of small disks at these redshifts or some unknown selection effects of deep ALMA observations. A larger sample of ALMA sources with NIRCam observations will be able to address this question.

The massive galaxy cluster El Gordo (z = 0.87) imprints multitudes of gravitationally lensed arcs onto James Webb Space Telescope Near-Infrared Camera (NIRCam) images. Eight bands of NIRCam imaging were obtained in the “Prime Extragalactic Areas for Reionization and Lensing Science” (“PEARLS”) program. Point-spread function–matched photometry across Hubble Space Telescope and NIRCam filters supplies new photometric redshifts. A new light-traces-mass lens model based on 56 image multiplicities identifies the two mass peaks and yields a mass estimate within 500 kpc of (7.0 ± 0.30) × 1014 M ⊙. A search for substructure in the 140 cluster members with spectroscopic redshifts confirms the two main mass components. The southeastern mass peak that contains the brightest cluster galaxy is more tightly bound than the northwestern one. The virial mass within 1.7 Mpc is (5.1 ± 0.60)×1014 M ⊙, lower than the lensing mass. A significant transverse velocity component could mean the virial mass is underestimated. We contribute one new member to the previously known z = 4.32 galaxy group. Intrinsic (delensed) positions of the five secure group members span a physical extent of ∼60 kpc. 13 additional candidates selected by spectroscopic/photometric constraints are small and faint, with a mean intrinsic luminosity ∼2.2 mag fainter than L *. NIRCam imaging admits a fairly wide range of brightnesses and morphologies for the group members, suggesting a more diverse galaxy population in this galaxy overdensity.