The quest to discover the most distant galaxies has developed rapidly in the last decade. We are now exploring redshifts of 6 and beyond, when the Universe was less than a billion years old, an epoch when the previously-neutral intergalactic medium was reionized. The continuing discovery of galaxies at progressively higher and higher redshifts has been driven by the availability of large telescopes on the ground and in space, improvements in detector technology, and new search strategies. Over the past 4 years, the Lyman break technique has been shown to be effective in isolating z~6 star-forming i'-drop galaxies through spectroscopic confirmation with large ground-based telescopes (Keck, Gemini and the ESO VLTs). Narrow-band imaging, notably with the wide field of the Subaru telescope, has also produced samples of Lyman-alpha emitters at these redshifts. A The discovery of this i'-drop galaxy population has been used to infer the global star formation rate density at this epoch (z~6), and we are now beginning to constrain the contribution to reionization of the UV flux from these galaxies. Infrared data from the Spitzer Space Telescope has been used to determine the spectral energy distributions (SEDs) from the rest-frame UV to the optical, and constrain the previous star formation histories, masses and ages. The indications are that much of the stellar mass of these galaxies might have formed in vigorous bursts at z>6. The next big advances would be to test the population synthesis modelling of these z~6 galaxies through spectroscopy of the rest-frame optical (rather than crude broad-band SEDs), and also to push the observational horizon for galaxies further to directly explore star formation during the reionization epoch. JWST is likely to have a profound impact on realising these goals.

We present a new measurement of the integrated stellar mass per comoving volume at redshift 5 determined via spectral energy fitting drawn from a sample of 214 photometrically-selected galaxies with z'<26.5 in the southern GOODS field. Following procedures introduced by Eyles et al. (2005), we estimate stellar masses for various sub-samples for which reliable and unconfused Spitzer IRAC detections are available. A spectroscopic sample of 14 of the most luminous sources with <z>=4.92 provides a firm lower limit to the stellar mass density of 1e6 Msun/Mpc^3. Several galaxies in this sub-sample have masses of order 10^11 Msun implying significant earlier activity occurred in massive systems. We then consider a larger sample whose photometric redshifts in the publicly-available GOODS-MUSIC catalog lie in the range 4.4 <z 5.6. Before adopting the GOODS-MUSIC photometric redshifts, we check the accuracy of their photometry and explore the possibility of contamination by low-z galaxies and low-mass stars. After excising probable stellar contaminants and using the z'-J color to exclude any remaining foreground red galaxies, we conclude that 196 sources are likely to be at z~5. The implied mass density from the unconfused IRAC fraction of this sample, scaled to the total available, is 6e6 Msun/Mpc^3. We discuss the uncertainties as well as the likelihood that we have underestimated the true mass density. Including fainter and quiescent sources the total integrated density could be as high as 1e7 Msun/Mpc^3. Using the currently available (but highly uncertain) rate of decline in the star formationhistory over 5 <z< 10, a better fit is obtained for the assembled mass at z~5 if we admit significant dust extinction at early times or extend the luminosity function to very faint limits. [abridged]

We measure the ages, stellar masses, and star formation histories of z~6 galaxies, observed within 1Gyr of the Big Bang. We use imaging from HST and Spitzer from the public "Great Observatories Origins Deep Survey", coupled with ground-based near-infrared imaging, to measure the spectral energy distributions from 0.8-5microns, spanning the rest-frame UV and optical. From our sample of ~50 i-drop Lyman-break star-forming galaxies in GOODS-South with z'(AB)<27mag, we focus on ~30 with reliable photometric or spectroscopic redshifts. Half of these are confused with foreground sources at Spitzer resolution, but from the 16 with clean photometry we find that a surprisingly large fraction (40%) have evidence for substantial Balmer-breaks. This indicates the presence of old underlying stellar populations that dominate the stellar masses. For these objects, we find ages of 200-700Myr, implying formation redshifts of 7<z<14, and large stellar masses in the range 1-3x10^10M_sun. Analysis of 7 i-drops that are undetected at 3.6microns indicates that these are younger, considerably less massive systems. We calculate that line contamination should not severely affect our photometry or derived results. Using data out to 8microns, we find little evidence for substantial intrinsic dust reddening. Correcting for incompleteness in our sample, we find a lower limit on the comoving stellar mass density at z~6 to be 2.5x10^6M_sun/Mpc^3. We are able to explore the star formation histories of our selected galaxies, and we infer that the past global star formation rate may have been much higher than that observed at z~6. The associated UV flux we infer at z>7 could have played a major role in reionizing the universe.

We determine the abundance of i'-band drop-outs in the HST/ACS GOODS surveys and the Hubble Ultra Deep Field (UDF). The majority of these sources are likely to be z~6 galaxies whose flux decrement arises from Lyman-alpha absorption. We have shown with Keck/DEIMOS spectroscopy that this technique does indeed select high redshift galaxies, and we discovered Lyman-alpha emission for about a third of the galaxies with z'_AB<25.6 The increased depth of UDF enables us to reach a ~10sigma limiting magnitude of z'_AB=28.5 (equivalent to 1.5M_sun/yr at z=6.1, or 0.1L*_UV for the z~3 U-drop population). The star formation rate at z~6 was approximately x6 LESS than at z~3. This declining comoving star formation rate poses an interesting challenge for models which suggest that L_UV>0.1L* star forming galaxies at z~6 reionized the universe. The short-fall in ionizing photons might be alleviated by galaxies fainter than our limit, or a radically different IMF. Alternatively, the bulk of reionization might have occurred at z>>6. We have recently discovered evidence of an early epoch of star formation in some of the i'-drops at z~6. Spitzer images with IRAC at 3.6-4.5microns show evidence of the age-sensitive Balmer/4000Ang break, dominated by stars older than 100Myr (and most probably 400Myr old). This pushes the formation epoch for these galaxies to z_form=7.5-13.5. There are at least some galaxies already assembled with stellar masses ~3x10^10M_sun (equivalent to 0.2M* today) within the first billion years. The early formation of such systems may have played a key role in reionizing the Universe at z~10.

We present new evidence for mature stellar populations with ages >100Myr in massive galaxies (M_stellar>10^10M_sun) seen at a time when the Universe was less than 1Gyr old. We analyse the prominent detections of two z~6 star-forming galaxies (SBM03#1 & #3) made at wavelengths corresponding to the rest-frame optical using the IRAC camera onboard the Spitzer Space Telescope. We had previously identified these galaxies in HST/ACS GOODS images of Chandra Deep Field South through the "i-drop" Lyman break technique, and subsequently confirmed spectroscopically with the Keck telescope. The new Spitzer photometry reveals significant Balmer/4000Ang discontinuities, indicative of dominant stellar populations with ages >100Myr. Fitting a range of population synthesis models (for normal initial mass functions) to the HST/Spitzer photometry yields ages of 250-650Myr and implied formation redshifts z~7.5-13.5 in presently-accepted world models. Remarkably, our sources have best-fit stellar masses of 1.3-3.8x10^10M_sun (95% confidence) assuming a Salpeter initial mass function. This indicates that at least some galaxies with stellar masses >20% of those of a present-day L* galaxy had already assembled within the first Gyr after the Big Bang. We also deduce that the past average star formation rate must be comparable to the current observed rate (SFR_UV~5-30M_sun/yr), suggesting that there may have been more vigorous episodes of star formation in such systems at higher redshifts. Although a small sample, limited primarily by Spitzer's detection efficiency, our result lends support to the hypothesis advocated in our earlier analyses of the Ultra Deep Field and GOODS HST/ACS data. The presence of established systems at z~6 suggests long-lived sources at earlier epochs (z>7) played a key role in reionizing the Universe.