The JWST/AURORA Survey: Multiple Balmer and Paschen Emission Lines for Individual Star-forming Galaxies at z=1.5-4.4. I. A Diversity of Nebular Attenuation Curves and Evidence for Non-Unity Dust Covering Fractions

We present the nebular attenuation curves and dust covering fractions for 24 redshift z=1.5-4.4 star-forming galaxies using multiple Balmer and Paschen lines from the JWST/AURORA survey. Nebular reddening derived from Paschen lines exceeds that from Balmer lines for at least half the galaxies in the sample when assuming the commonly-adopted Galactic extinction curve, implying the presence of optically-thick star formation. The nebular attenuation curves exhibit a broad range of normalizations (Rv ~ 3.2-16.4). Motivated by the offsets in reddening deduced from the Balmer and Paschen lines, and the high Rv values for the individual nebular attenuation curves, both of which suggest variations in the dust-stars geometry, we propose a model with a subunity dust covering fraction (fcov). Fitting such a model to the HI recombination line ratios indicates fcov ~ 0.6-1.0. The normalizations of the nebular attenuation curves, Rv, are driven primarily by fcov and the mix of optically-thick and thin OB associations. Thus, the diversity of nebular attenuation curves can be accommodated by assuming dust grain properties similar to that of Milky Way sightlines but with a subunity covering fraction of dust. Integrated measurements of multiple Balmer and Paschen lines can be used to place novel constraints on the dust covering fraction towards OB associations. These, in turn, provide new avenues for exploring the role of dust and gas covering fraction in a number of relevant aspects of high-redshift galaxies, including the impact of stellar feedback on ISM porosity and the escape of Ly-alpha and Lyman continuum radiation.