It is well-known that foreground subtraction in 21cm surveys removes large scale power. We investigate associated systematic biases. We show that removing line-of-sight fluctuations on large scales aliases into suppression of the 3D power spectrum across a broad range of scales. This bias can be eliminated by marginalizing over small k in the 1D power spectrum; however, the unbiased estimator will have unavoidably larger variance. We also show that Gaussian realizations of the power spectrum permit accurate and extremely rapid Monte-Carlo simulations for error analysis; repeated realizations of the fully non-Gaussian field are unnecessary. We perform Monte-Carlo maximum-likelihood simulations of foreground removal which yield unbiased, minimum variance estimates of the power spectrum in agreement with Fisher matrix estimates. Foreground removal also distorts the 21cm PDF, reducing the contrast between neutral and ionized regions. We show that it is the subtraction of large-scales modes which is responsible for this distortion, and that it is less severe in the earlier stages of reionization. It can be reduced by using larger bandwidths for foreground removal. In the late stages of reionization, the largest ionized regions (which consist of foreground emission only) provides calibration points which potentially allow recovery of large-scale modes. Finally, we also show that: (i) the broad frequency response of synchrotron and free-free emission will smear out any features in the electron momentum distribution and ensure spectrally smooth foregrounds; (ii) extragalactic radio recombination lines should be negligible foregrounds.

We assess the impact of Galactic synchrotron foreground removal on the observation of high-redshift quasar HII regions in redshifted 21-cm emission. We consider the case where a quasar is observed in an intergalactic medium (IGM) whose ionisation structure evolves slowly relative to the light crossing time of the HII region, as well as the case where the evolution is rapid. The latter case is expected towards the end of the reionisation era where the highest redshift luminous quasars will be observed. In the absence of foregrounds the fraction of neutral hydrogen in the IGM could be measured directly from the contrast between the HII region and surrounding IGM. However, we find that foreground removal lowers the observed contrast between the HII region and the IGM. This indicates that measurement of the neutral fraction would require modelling to correct for this systematic effect. On the other hand, foreground removal does not modify the most prominent features of the 21-cm maps. Using a simple algorithm we demonstrate that measurements of the size and shape of observed HII regions will not be affected by continuum foreground removal. Moreover, measurements of these quantities will not be adversely affected by the presence of a rapidly evolving IGM.

We analyze the photometric information contained in individual pixels of galaxies in the Hubble Deep Field North (HDFN) using a new technique, _pixel-z_, that combines predictions of evolutionary synthesis models with photometric redshift template fitting. Each spectral energy distribution template is a result of modeling of the detailed physical processes affecting gas properties and star formation efficiency. The criteria chosen to generate the SED templates is that of sampling a wide range of physical characteristics such as age, star formation rate, obscuration and metallicity. A key feature of our method is the sophisticated use of error analysis to generate error maps that define the reliability of the template fitting on pixel scales and allow for the separation of the interplay among dust, metallicity and star formation histories. This technique offers a number of advantages over traditional integrated color studies. As a first application, we derive the star formation and metallicity histories of galaxies in the HDFN. Our results show that the comoving density of star formation rate, determined from the UV luminosity density of sources in the HDFN, increases monotonically with redshift out to at least redshift of 5. This behavior can plausibly be explained by a smooth increase of the UV luminosity density with redshift coupled with an increase in the number of star forming regions as a function of redshift. We also find that the information contained in individual pixels in a galaxy can be linked to its morphological history. Finally, we derive the metal enrichment rate history of the universe and find it in good agreement with predictions based on the evolving HI content of Lyman-alpha QSO absorption line systems.