We present a synopsis of the project to establish 32 new faint (16.5 <= V <= 19.8) DA white dwarf (DAWD) spectrophotometric standards distributed over the whole sky. Our results validate the use of fully radiative pure hydrogen model fluxes for hot DA white dwarfs to predict the observed broadband fluxes from near-ultraviolet through the near-infrared to accuracies of a few parts per thousand. After fitting the line of sight reddenings simultaneously with the model spectral energy distributions of these stars against spectroscopic and multiband photometric observations, we have shown that residuals have an rms of typically 0.4%. This indicates that the complications from interstellar dust extinction have been adequately mitigated. Our stars supplement the three brighter DAWDs that define the flux scale of CALSPEC. The consequent photometric accuracy, their all-sky coverage, and their brightness range that matches the dynamic range of large telescopes, constitute an unprecedented ensemble of standard stars for both ground as well as space based use. This paper targets readers who may wish to use these as standard stars and provides them with the essential content to understand their strengths and limitations, without traversing the technical details of analysis that are already captured in a series of papers since 2016. The narrative here describes the motivation, justification, and evolution of the analysis methods; the input data that constrain the modeling; as well as the stability of our results in the face of future improvements in models.

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (16.5 < V < 19.5) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zeropoints and WD parameters (surface gravity log g, effective temperature T_eff, extinction A_V, dust relation parameter R_V) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy and ground-based optical spectroscopy to sub-percent precision. Photometric residuals for the sample are the lowest yet yielding <0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory's Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and VISible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic survey fields: GOODS-N, GOODS-S, EGS, and COSMOS. We introduce this survey by presenting a comprehensive analysis of the absolute escape fraction ( fescabs ) of Lyman continuum radiation through stacking the UV images of a population of star-forming galaxies with secure redshifts at 2.4 <= z <= 3.0. Our stacking benefits from the catalogs of high-quality spectroscopic redshifts compiled from archival ground-based data and HST slitless spectroscopy, carefully vetted by dedicated visual inspection efforts. We develop a robust stacking method to apply to 10 samples of in total 56 galaxies, and perform detailed Monte Carlo simulations of the intergalactic medium (IGM) attenuation, to take into account the sample variance of the mean IGM transmission when measuring fescabs . The full stack at z approximate to 2.44 from 28 galaxies places a stringent 1 sigma upper limit of fescabs less than or similar to 5% , whereas the full stack at z approximate to 2.72 of an equal number of galaxies gives an upper limit of fescabs less than or similar to 26% at 1 sigma confidence level. These new F275W and F435W imaging mosaics from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes.

The SKYSURF project constrained extragalactic background light and diffuse light (DL) with the vast archive of Hubble Space Telescope (HST) images. Thermal emission from HST itself introduces an additional uncertain background and hinders accurate measurement of the DL level. Here, we use archival Wide Field Camera 3 (WFC3)/IR engineering data to investigate and model changes in the temperature of various components in HSTs optical path as a function of time (solar cycle) and time of the year (Earth–Sun distance). We also specifically investigate changes in temperature with HST's orbital phase and time since Earth occultation. We investigate possible correlations between HST component temperature and year, and temperature and month. The thermal background changes by less than one Kelvin in the WFC3 pickoff mirror, one of the most important contributors to the thermal background. We model these data to describe the impact that orbital phase, year, and time of year have on the HST and WFC3 component temperatures, and use this to derive the impact on the thermal dark signal and the resulting DL measurements. Based on this improved modeling, we provide new upper limits on the level of DL of 21, 32, and 25 nW m ^−2 sr ^−1 for F125W, F140W, and F160W. Additionally, by accounting for all known sources of measurement uncertainty, we report lower limits on the level of DL of 12, 20, and 2 nW m ^−2 sr ^−1 for F125W, F140W, and F160W.

UVCANDELS is a Hubble Space Telescope Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields-GOODS-N, GOODS-S, EGS, and COSMOS-covering a total area of similar to 426 arcmin(2). This is similar to 2.7 times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag (5 sigma in 0."2 apertures) for F275W and F435W, respectively. Along with new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method, yielding a factor of 1.5 increase over H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well-tested photometric redshift measurements, we identify 5810 galaxies at redshifts 0.6 < z < 1, down to an absolute magnitude of M-UV = -14.2. In order to minimize the effect of uncertainties in estimating the completeness function, especially at the faint end, we restrict our analysis to sources above 30% completeness, which provides a final sample of 4726 galaxies at -21.5 < M-UV < -15.5. We performed a maximum likelihood estimate to derive the best-fit parameters of the UV LF. We report a best-fit faint-end slope of alpha=-1.359(-0.041)(+0.041) at z similar to 0.8. Creating subsamples at z similar to 0.7 and z similar to 0.9, we observe a possible evolution of alpha with redshift. The unobscured UV luminosity density at M-UV < -10 is derived as rho(UV)=1.339(-0.030)(+0.027)(x10(26)ergs(-1)Hz(-1)Mpc(-3)) using our best-fit LF parameters. The new F275W and F435 photometric catalogs from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes.

We present the results of our search for Lyman continuum (LyC)-emitting (weak) active galactic nuclei (AGN) at redshifts 2.3 ≲ z ≲ 4.9 from Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) F275W observations in the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS) fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W imaging found in Early Release Science (ERS) and Hubble Deep UV Legacy Survey data. We performed exhaustive queries of the Vizier database to locate AGN with high-quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Out of these 51, we find 12 AGN that had ≥4 σ detected LyC flux in the WFC3/UVIS images. Using a wide variety of space-based plus ground-based data, ranging from X-ray to radio wavelengths, we fit the multiwavelength photometric data of each AGN to a CIGALE spectral energy distribution (SED) using AGN models and correlate various SED parameters to the LyC flux. Kolmogorov–Smirnov tests of the SED parameter distributions for the LyC-detected and nondetected AGN showed they are likely not distinct samples. However, we find that the X-ray luminosity, star formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlations of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission and 6 μ m luminosity, and anticorrelations with metallicity and A _FUV . We simulate the LyC escape fraction ( f _esc ) using the CIGALE and intergalactic medium transmission models for the LyC-detected AGN and find an average f _esc ≃ 18%, weighted by uncertainties. We stack the LyC fluxes of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our submillimeter-detected F336W sample (3.15 < z < 3.71) shows the brightest stacked LyC flux. These findings indicate that LyC production and escape in AGN are more complicated than the simple assumption of thermal emission and a 100% escape fraction. Further testing of AGN models with larger samples than presented here is needed.

Abstract The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey provided ultraviolet F275W imaging with coordinated parallel optical F435W imaging in four of the five CANDELS fields: GOODS-N, GOODS-S, EGS, and COSMOS, covering a total area of ∼426 arcmin2. UVCANDELS takes primary WFC3/UVIS F275W exposures at a uniform 3-orbit depth and ACS F435W exposures (in parallel) at slightly varying depth due to the roll angle constraints and the overlap from the increased field of view of the ACS camera, reaching a limiting magnitude of ∼27 and ∼28 ABmag (5σ in 0.″2 apertures) for F275W and F435W, respectively. We present the results of the UVCANDELS observations, custom calibrations, and the creation of F275W and F435W imaging mosaics, which have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes.

The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides deep Hubble Space Telescope (HST) F275W and F435W imaging over four CANDELS fields (GOODS-N, GOODS-S, COSMOS, and Extended Groth Strip). We combine this newly acquired UV imaging with existing HST imaging from CANDELS as well as existing ancillary data to obtain robust photometric redshifts and reliable estimates for galaxy physical properties for over 150,000 galaxies in the similar to 430 arcmin2 UVCANDELS area. Here, we leverage the power of the new UV photometry to not only improve the photometric redshift measurements in these fields, but also constrain the full redshift probability distribution combining multiple redshift-fitting tools. Furthermore, using the full UV-to-IR photometric data set, we measure the galaxy physical properties by fitting templates from population synthesis models with two different parameterizations (flexible and fixed form) of the star formation histories (SFHs). Compared to the flexible SFH parameterization, we find that the fixed-form SFHs systematically underestimate the galaxy stellar masses, both at the low-mass (less than or similar to 109 M circle dot) and high-mass (greater than or similar to 1010 M circle dot) end, by as much as similar to 0.5 dex. This underestimation is primarily due the limited ability of fixed-form SFH parameterization to simultaneously capture the chaotic nature of star formation in these galaxies.

We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3≲z≲4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Of these 51, we find 12 AGN had ≥4σ detected LyC flux in the WFC3/UVIS images. Using space- and ground-based data from X-ray to radio, we fit the multi-wavelength photometric data of each AGN to a CIGALE SED and correlate various SED parameters to the LyC flux. KS-tests of the SED parameter distributions for the LyC-detected and non-detected AGN showed they are likely not distinct samples. However, we find that X-ray luminosity, star-formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlation of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission, 6 μ m luminosity, and anti-correlation with metallicity and A_FUV. We simulate the LyC escape fraction (f_esc) using the CIGALE and IGM transmission models for the LyC-detected AGN and find an average f_esc≃18 weighted by uncertainties. We stack the LyC flux of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our sub-mm detected F336W sample shows the brightest stacked LyC flux. These findings indicate that LyC-production and -escape in AGN is more complicated than the simple assumption of thermal emission and a 100 AGN models with larger samples than presented here is needed.

Auxiliary data files for Axelrod et al. 2023 (in prep.). "HST_passbands" - Hubble Space Telescope photometric filter passthroughs. First column is vacuum wavelength in nm, second column is response. "SED" - calibrated spectral flux density distributions. First column is vacuum wavelength in Angstrom, second column is SED in erg/sec/cm**2/Angstrom. "Tables" - paper tables in LaTeX.

Hot DA white dwarfs have fully radiative pure hydrogen atmospheres that are the least complicated to model. Pulsationally stable, they are fully characterized by their effective temperature Teff, and surface gravity log g, which can be deduced from their optical spectra and used in model atmospheres to predict their spectral energy distribution (SED). Based on this, three bright DAWDs have defined the spectrophotometric flux scale of the CALSPEC system of HST. In this paper we add 32 new fainter (16.5 < V < 19.5) DAWDs spread over the whole sky and within the dynamic range of large telescopes. Using ground based spectra and panchromatic photometry with HST/WFC3, a new hierarchical analysis process demonstrates consistency between model and observed fluxes above the terrestrial atmosphere to < 0.004 mag rms from 2700 Å to 7750 Å and to 0.008 mag rms at 1.6μm for the total set of 35 DAWDs. These DAWDs are thus established as spectrophotometric standards with unprecedented accuracy from the near ultraviolet to the near-infrared, suitable for both ground and space based observatories. They are embedded in existing surveys like SDSS, PanSTARRS and GAIA, and will be naturally included in the LSST survey by Rubin Observatory. With additional data and analysis to extend the validity of their SEDs further into the IR, these spectrophotometric standard stars could be used for JWST, as well as for the Roman and Euclid observatories.

The diffuse, unresolved sky provides most of the photons that the Hubble Space Telescope (HST) receives, yet remains poorly understood. The HST Archival Legacy program SKYSURF aims to measure the 0.2-1.6 mu m sky surface brightness (sky-SB) from over 140,000 HST images. We describe a sky-SB measurement algorithm designed for SKYSURF that is able to recover the input sky-SB from simulated images to within 1% uncertainty. We present our sky-SB measurements estimated using this algorithm on the entire SKYSURF database. Comparing our sky-SB spectral energy distribution (SED) to measurements from the literature shows general agreements. Our SKYSURF SED also reveals a possible dependence on the Sun angle, indicating either nonisotropic scattering of solar photons off interplanetary dust or an additional component to zodiacal light. Finally, we update the diffuse light limits in the near-IR based on the methods from Carleton et al., with values of 0.009 MJy sr(-1) (22 nW m(-2) sr(-1)) at 1.25 mu m, 0.015 MJy sr(-1) (32 nW m(-2) sr(-1)) at 1.4 mu m, and 0.013 MJy sr(-1 )(25 nW m(-2) sr(-1)) at 1.6 mu m. These estimates provide the most stringent all-sky constraints to date in this wavelength range. SKYSURF sky-SB measurements are made public on the official SKYSURF website and will be used to constrain diffuse light in future papers.

We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in 0.2 < z < 0.8 groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20 k group catalog. Limiting our analysis to a complete sample of log ( M * / M ⊙ ) > 10.56 group galaxies, we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow (∼Gyr) timescales. We find a quenching time of approximately 4.91 − 1.47 + 0.91 Gyr, consistent with the slow processes of strangulation and delayed-then-rapid quenching although more data are needed to confirm this result.

The NASA/ESA Hubble Space Telescope provides us with many lessons regarding the interplay between science and technology that will be applicable to future telescopes. In this talk, we will consider several “case studies” where technological choices and advances had profound impacts on the science from HST and, conversely, where science needs drove technical and operational plans. HST is unique amongst astronomy space missions by virtue of its multiple upgrades during servicing missions. Furthermore, HST has been fortunate to have operated during a period of dramatic improvements in computational and communications capabilities which have changed how the scientific community interacts with both the observatory and its data. We will discuss how this has influenced both the evolution of astronomical research, the continuing development HST’s systems, and implications for future missions.

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.

We verified for photometric stability a set of DA white dwarfs with Hubble Space Telescope magnitudes from the near-ultraviolet to the near-infrared and ground-based spectroscopy by using time-spaced observations from the Las Cumbres Observatory network of telescopes. The initial list of 38 stars was whittled to 32 final ones which comprise a high quality set of spectrophotometric standards. These stars are homogeneously distributed around the sky and are all fainter than r 16.5 mag. Their distribution is such that at least two of them would be available to be observed from any observatory on the ground at any time at airmass less than two. Light curves and different variability indices from the Las Cumbres Observatory data were used to determine the stability of the candidate standards. When available, Pan-STARRS1, Zwicky Transient Facility and TESS data were also used to confirm the star classification. Our analysis showed that four DA white dwarfs may exhibit evidence of photometric variability, while a fifth is cooler than our established lower temperature limit, and a sixth star might be a binary. In some instances, due to the presence of faint nearby red sources, care should be used when observing a few of the spectrophotometric standards with ground-based telescopes. Light curves and finding charts for all the stars are provided.

We present the first results from the HST Archival Legacy project "SKYSURF." As described in Windhorst et al. 2022, SKYSURF utilizes the large HST archive to study the diffuse UV, optical, and near-IR backgrounds and foregrounds in detail. Here we utilize SKYSURF's first sky-surface brightness measurements to constrain the level of near-IR diffuse Extragalactic Background Light (EBL). Our sky-surface brightness measurements have been verified to an accuracy of better than 1%, which when combined with systematic errors associated with HST, results in sky brightness uncertainties of $\sim$2-4% $\simeq$ 0.005 MJy/sr in each image. We put limits on the amount of diffuse EBL in three near-IR filters (F125W, F140W, and F160W) by comparing our preliminary sky measurements of $> 30,000$ images to Zodiacal light models, carefully selecting the darkest images to avoid contamination from stray light. In addition, we investigate the impact that instrumental thermal emission has on our measurements, finding that it has a limited impact on F125W and F140W measurements, whereas uncertainties in the exact thermal state of HST results in significant uncertainties in the level of astrophysical diffuse light in F160W images. When compared to the Kelsall et al. (1998) Zodiacal model, an isotropic diffuse background of $30$ nW m$^{-2}$ sr$^{-1}$ remains, whereas using the Wright (1998) Zodiacal model results in no discernible diffuse background. Based primarily on uncertainties in the foreground model subtraction, we present limits on the amount of diffuse EBL of 29 nW m$^{-2}$ sr$^{-1}$, 40 nW m$^{-2}$ sr$^{-1}$, and 29 nW m$^{-2}$ sr$^{-1}$ for F125W, F140W, and F160W respectively. While this light is generally isotropic, our modeling at this point does not distinguish between a cosmological origin or a Solar System origin (such as a dim, diffuse, spherical cloud of cometary dust).

We present the H_160 morphological catalogs for the COSMOS-DASH survey, the largest area near-IR survey using HST-WFC3 to date. Utilizing the "Drift And SHift" observing technique for HST-WFC3 imaging, the COSMOS-DASH survey imaged approximately 0.5 deg^2 of the UltraVISTA deep stripes (0.7 deg^2 when combined with archival data). Global structural parameters are measured for 51,586 galaxies within COSMOS-DASH using GALFIT (excluding the CANDELS area) with detection using a deep multi-band HST image. We recover consistent results with those from the deeper 3D-HST morphological catalogs, finding that, in general, sizes and Sérsic indices of typical galaxies are accurate to limiting magnitudes of H_160<23 and H_160<22 ABmag, respectively. In size-mass parameter space, galaxies in COSMOS-DASH demonstrate robust morphological measurements out to z∼2 and down to log(M_⋆/M_⊙)∼9. With the advantage of the larger area of COSMOS-DASH, we measure a flattening of the quiescent size-mass relation below log(M_⋆/M_⊙)∼10.5 that persists out to z∼2. We show that environment is not the primary driver of this flattening, at least out to z=1.2, whereas internal physical processes may instead govern the structural evolution.

We report observations from HST of Cepheids in the hosts of 42 SNe Ia used to calibrate the Hubble constant (H0). These include all suitable SNe Ia in the last 40 years at z<0.01, measured with >1000 orbits, more than doubling the sample whose size limits the precision of H0. The Cepheids are calibrated geometrically from Gaia EDR3 parallaxes, masers in N4258 (here tripling that Cepheid sample), and DEBs in the LMC. The Cepheids were measured with the same WFC3 instrument and filters (F555W, F814W, F160W) to negate zeropoint errors. We present multiple verifications of Cepheid photometry and tests of background determinations that show measurements are accurate in the presence of crowding. The SNe calibrate the mag-z relation from the new Pantheon+ compilation, accounting here for covariance between all SN data, with host properties and SN surveys matched to negate differences. We decrease the uncertainty in H0 to 1 km/s/Mpc with systematics. We present a comprehensive set of 70 analysis variants to explore the sensitivity of H0 to selections of anchors, SN surveys, z range, variations in the analysis of dust, metallicity, form of the P-L relation, SN color, flows, sample bifurcations, and simultaneous measurement of H(z). Our baseline result from the Cepheid-SN sample is H0=73.04+-1.04 km/s/Mpc, which includes systematics and lies near the median of all analysis variants. We demonstrate consistency with measures from HST of the TRGB between SN hosts and NGC 4258 with Cepheids and together these yield 72.53+-0.99. Including high-z SN Ia we find H0=73.30+-1.04 with q0=-0.51+-0.024. We find a 5-sigma difference with H0 predicted by Planck+LCDM, with no indication this arises from measurement errors or analysis variations considered to date. The source of this now long-standing discrepancy between direct and cosmological routes to determining the Hubble constant remains unknown.

Extragalactic background light (EBL) studies have revealed a significant discrepancy between direct measurements—via instruments measuring “bare” sky from which Zodiacal and Galactic light models are subtracted—and measurements of the integrated galaxy light (IGL). This discrepancy could lie in either method, whether it be an incomplete Zodiacal model or missed faint galaxies in the IGL calculations. It has been proposed that the discrepancy is due to deep galaxy surveys, such as those with the Hubble Space Telescope, missing up to half of the faint galaxies with 24 ≲ m AB ≲ 29 mag. We address this possibility by simulating higher number densities of galaxies, and so assess incompleteness due to object overlap, with three replications of the Hubble UltraDeep Field (HUDF). SourceExtractor is used to compare the recovered counts and photometry to the original HUDF, allowing us to assess how many galaxies may have been missed due to confusion, i.e., due to blending with neighboring faint galaxies. This exercise reveals that, while up to 50% of faint galaxies with 28 ≲ m AB ≲ 29 mag were missed or blended with neighboring objects in certain filters, not enough were missed to account for the EBL discrepancy alone in any of the replications.