Strength in Numbers: Red Galaxies Bolster the Cosmic Star Formation Rate Density at z > 3

A comprehensive account of the cosmic star-formation history demands an accurate census of dust-enshrouded star formation over cosmic time. We provide strong new constraints from a large sample of 777 red galaxies, selected based on their dust-reddened, rest-frame UV-optical emission. This sample of 777 galaxies spans $1 < z < 8$ and is selected from PRIMER JWST NIRCam and HST COSMOS optical data, ensuring robust colour criteria. The SEDs indicate that these dust-reddened galaxies are star-forming, with median $\mathrm{SFR \sim 40M_{\odot}yr^{-1}}$ and stellar mass $\log(M_{*}/M_{\odot}) = 10.3^{+0.6}_{-0.8}$; each exceeds the corresponding medians of the full JWST-detected population by over two dex. Our sample thus clearly shows that red galaxies dominate the high-mass end: they comprise 72 \% of galaxies with $\log(M/M_{\odot}) > 10$ at $z = 3.3$, rising to 91\% by $z \sim 7$ (albeit with large uncertainties at the highest redshifts). Crucially, we find that the number density of massive red star-forming galaxies at $z \sim 6$ is sufficient to explain the abundance of quiescent galaxies at $z > 3$, consistent with typical quenching timescales allowed in the $\mathrm{\sim 1Gyr}$ interval from $z \sim 6$ to $z \sim 3$. This large abundance yields a substantial contribution to the cosmic star-formation rate density: at $z \sim 4$, red galaxies provide $\mathrm {ρ_{SFR} = 3.9^{+0.6}_{-0.5} \times 10^{-2} M_{\odot} yr^{-1}Mpc^{-3}}$, and at $z \sim 5$ they supply nearly 40 \% of the total $ρ_{SFR}$. This exceeds the contribution of bright sub(mm)-selected dusty star-forming galaxies by more than an order of magnitude. Future deeper and wider ALMA surveys will provide further opportunities to strengthen and extend our results in our quest to fully quantify the contribution of dust-obscured activity to $ρ_{\mathrm{SFR}}$ at high redshifts.