Identification of z~>2 Herschel 500 micron sources using
  color-deconfusion
  We present a new method to search for candidate z~>2 Herschel 500{\mu}m
sources in the GOODS-North field, using a S500{\mu}m/S24{\mu}m "color
deconfusion" technique. Potential high-z sources are selected against
low-redshift ones from their large 500{\mu}m to 24{\mu}m flux density ratios.
By effectively reducing the contribution from low-redshift populations to the
observed 500{\mu}m emission, we are able to identify counterparts to high-z
500{\mu}m sources whose 24{\mu}m fluxes are relatively faint. The recovery of
known z~4 starbursts confirms the efficiency of this approach in selecting
high-z Herschel sources. The resulting sample consists of 34 dusty star-forming
galaxies at z~>2. The inferred infrared luminosities are in the range
1.5x10^12-1.8x10^13 Lsun, corresponding to dust-obscured star formation rates
(SFRs) of ~260-3100 Msun/yr for a Salpeter IMF. Comparison with previous SCUBA
850{\mu}m-selected galaxy samples shows that our method is more efficient at
selecting high-z dusty galaxies with a median redshift of z=3.07+/-0.83 and 10
of the sources at z~>4. We find that at a fixed luminosity, the dust
temperature is ~5K cooler than that expected from the Td-LIR relation at z<1,
though different temperature selection effects should be taken into account.
The radio-detected subsample (excluding three strong AGN) follows the
far-infrared/radio correlation at lower redshifts, and no evolution with
redshift is observed out to z~5, suggesting that the far-infrared emission is
star formation dominated. The contribution of the high-z Herschel 500{\mu}m
sources to the cosmic SFR density is comparable to that of SMG populations at
z~2.5 and at least 40% of the extinction-corrected UV samples at z~4
(abridged).
