Until recently, only a handful of dusty, star-forming galaxies (DSFGs) were known at z>4, most of them significantly amplified by gravitational lensing. Here, we have increased the number of such DSFGs substantially, selecting galaxies from the uniquely wide 250, 350, and 500{mu}m Herschel-ATLAS imaging survey on the basis of their extremely red far-infrared colors and faint 350 and 500{mu}m flux densities, based on which, they are expected to be largely unlensed, luminous, rare, and very distant. The addition of ground-based continuum photometry at longer wavelengths from the James Clerk Maxwell Telescope and the Atacama Pathfinder Experiment allows us to identify the dust peak in their spectral energy distributions (SEDs), with which we can better constrain their redshifts. We select the SED templates that are best able to determine photometric redshifts using a sample of 69 high-redshift, lensed DSFGs, then perform checks to assess the impact of the CMB on our technique, and to quantify the systematic uncertainty associated with our photometric redshifts, {sigma}=0.14(1+z), using a sample of 25 galaxies with spectroscopic redshifts, each consistent with our color selection. For Herschel-selected ultrared galaxies with typical colors of S_500_/S_250_~2.2 and S_500_/S_350_~1.3 and flux densities, S500~50mJy, we determine a median redshift, z_phot_=3.66, an interquartile redshift range, 3.30-4.27, with a median rest-frame 8-1000{mu}m luminosity, L_IR_, of 1.3x10^13^L_{sun}. A third of the galaxies lie at z>4, suggesting a space density, {rho}_z>4, of ~6x10^-7^Mpc^-3^. Our sample contains the most luminous known star-forming galaxies, and the most overdense cluster of starbursting proto-ellipticals found to date.

Cone search capability for table J/ApJ/832/78/dsfg (Targets and their photometric properties (Tables 1 and 2))