There Is More to Outshining: 2D Dust Effects on Stellar Mass Estimates at $3 \leq z < 9$ with JWST in the JADES Field

TL;DR

This study investigates how dust attenuation affects stellar mass estimates at high redshift using JWST data, revealing that integrated SED fitting underestimates masses due to dust and age effects, with resolved analyses providing more accurate results.

Contribution

It introduces a detailed comparison between resolved and integrated SED fitting for high-redshift galaxies, highlighting the impact of dust and age on stellar mass estimates.

Findings

Resolved stellar masses are systematically higher than integrated estimates.

Dust attenuation and age differences drive discrepancies in mass estimates.

Integrated fits underestimate stellar ages by 23-31%.

Abstract

Dust attenuation modifies the observed spectral energy distribution (SED), leading to biases in the properties inferred from integrated SED fitting. As spatially resolved SED modeling becomes feasible for large high-redshift samples, it is increasingly important to assess how dust attenuation affects resolved mass estimates. We evaluate the impact of dust attenuation on stellar mass estimates derived from integrating spatially resolved SED fitting results. We perform spatially resolved and integrated SED fitting on a sample of 3408 galaxies at $3 \leq z < 9$ from the GOODS South field, combining deep NIRCam from the JWST Advanced Deep Extragalactic Survey (JADES) and HST/ACS imaging from GOODS and CANDELS. We compare galaxy-integrated properties derived from fitting the summed SED with those obtained from spatially resolved SED modeling. Using a two-component dust attenuation model with a variable slope, we investigate how the dust attenuation slope, A(V), and stellar population properties contribute to discrepancies in the resulting stellar mass estimates. Resolved stellar masses are systematically higher than integrated estimates, with a median offset of +0.24 dex. Resolved analyses recover higher dust attenuations ($\Delta A(V)\approx +0.08$ mag), lower birth cloud fractions ($\Delta\mu \approx -0.28$), and grayer attenuation curves ($\Delta\delta_{\mathrm{ISM}} = +0.08$), arising from preferential sampling of compact star-forming regions. Integrated fits underestimate stellar ages by $\sim23\%$ at $z < 5$ and 31$\%$ at $z \gtrsim 5$. The stellar mass offset correlates strongly with the age difference and the attenuation slope difference, indicating that age-dependent outshining and spatially varying dust geometry are primary drivers of the discrepancy between resolved and integrated stellar masses.