We introduce the largest to date survey of massive quiescent galaxies at redshift z~1.6. With these data, which cover 27.6 deg^2, we can find significant numbers of very rare objects such as ultra-massive quiescent galaxies that populate the extreme massive end of the galaxy mass function, or dense environments that are likely to become present-day massive galaxy clusters. In this paper, the first in a series, we apply our gzK adaptation of the BzK technique to select our z~1.6 galaxy catalog and then study the quiescent galaxy stellar mass function with good statistics over M_stars ~ 10^10.2 - 10^11.7 M_sun --- a factor of 30 in mass --- including 60 ultra-massive z~1.6 quiescent galaxies with M_stars > 10^11.5 M_sun. We find that the stellar mass function of quiescent galaxies at z~1.6 is well represented by the Schechter function over this large mass range. This suggests that the mass quenching mechanism observed at lower redshifts must have already been well established by this epoch, and that it is likely due to a single physical mechanism over a wide range of mass. This close adherence to the Schechter shape also suggests that neither merging nor gravitational lensing significantly affect the observed quenched population. Finally, comparing measurements of M* parameters for quiescent and star-forming populations (ours and from the literature), we find hints of an offset (M*_SF > M*_PE), that could suggest that the efficiency of the quenching process evolves with time.

We perform spatially-resolved, pixel-by-pixel SED fitting on galaxies up to $z\sim2.5$ in the Hubble Extreme Deep Field (XDF). Comparing stellar mass estimates from spatially resolved and spatially unresolved photometry we find that unresolved masses can be systematically underestimated by factors of up to 5. The ratio of the unresolved to resolved mass measurement depends on the galaxy's specific star formation rate (sSFR): at low sSFRs the bias is small, but above sSFR$\ \sim 10^{-9.5}$ yr$^{-1}$ the discrepancy increases rapidly such that galaxies with sSFRs$\ \sim 10^{-8}$ yr$^{-1}$ have unresolved mass estimates of only one half to one fifth of the resolved value. This result indicates that stellar masses estimated from spatially-unresolved datasets need to be systematically corrected, in some cases by large amounts, and we provide an analytic prescription for applying this correction. We show that correcting stellar mass measurements for this bias changes the normalization and slope of the star-forming main sequence and reduces its intrinsic width; most dramatically, correcting for the mass bias increases the stellar mass density of the Universe at high redshift and can resolve the long-standing discrepancy between the directly-measured cosmic star formation rate density at $z\gtrsim1$ and that inferred from stellar mass densities ("the missing mass problem").

We fit model spectral energy distributions to each pixel in 67 nearby (<z>=0.0057) galaxies using broadband photometry from the Sloan Digital Sky Survey and GALEX. For each galaxy, we compare the stellar mass derived by summing the mass of each pixel to that found from fitting the entire galaxy treated as an unresolved point source. We find that, while the pixel-by-pixel and unresolved masses of galaxies with low specific star formation rates (such as ellipticals and lenticulars) are in rough agreement, the unresolved mass estimate for star-forming galaxies is systematically lower then the measurement from spatially-resolved photometry. The discrepancy is strongly correlated with sSFR, with the highest sSFRs in our sample having masses underestimated by 25% (0.12 dex) when treated as point sources. We found a simple relation to statistically correct mass estimates derived from unresolved broad-band SED fitting to the resolved mass estimates: m_{resolved} = m_{unresolved}/(-0.057log(sSFR) + 0.34) where sSFR is in units of yr^{-1}. We study the effect of varying spatial resolution by degrading the image resolution of the largest images and find a sharp decrease in the pixel-by-pixel mass estimate at a physical scale of approximately 3 kpc, which is comparable to spiral arm widths. The effects we observe are consistent with the "outshining" idea which posits that the youngest stellar populations mask more massive, older -- and thus fainter -- stellar populations. Although the presence of strong dust lanes can also lead to a drastic difference between resolved and unresolved mass estimates (up to 45% or 0.3 dex) for any individual galaxy, we found that resolving dust does not affect mass estimates on average. The strong correlation between mass discrepancy and sSFR is thus most likely due to the outshining systematic bias.

We study how the addition of on-board optical photometric bands to future space-based weak lensing instruments could affect the photometric redshift estimation of galaxies, and hence improve estimations of the dark energy parameters through weak lensing. Basing our study on the current proposed Euclid configuration and using a mock catalog of galaxy observations, various on-board options are tested and compared with the use of ground-based observations from the Large Synoptic Survey Telescope (LSST) and Pan-STARRS. Comparisons are made through the use of the dark energy Figure of Merit, which provides a quantifiable measure of the change in the quality of the scientific results that can be obtained in each scenario. Effects of systematic offsets between LSST and Euclid photometric calibration are also studied. We find that adding two (U and G) or even one (U) on-board optical band-passes to the space-based infrared instrument greatly improves its photometric redshift performance, bringing it close to the level that would be achieved by combining observations from both space-based and ground-based surveys while freeing the space mission from reliance on external datasets.

We explore the feasibility and limitations of using the 1.6um bump as a photometric redshift indicator and selection technique and use it to study the rest-frame H-band galaxy luminosity and stellar mass functions at redshift z~2. We use publicly available Spitzer/IRAC images in the GOODS fields and find that color selection in the IRAC bandpasses alone is comparable in completeness and contamination to BzK selection. We find that the shape of the 1.6um bump is robust, and photometric redshifts are not greatly affected by choice of model parameters. Comparison with spectroscopic redshifts shows photometric redshifts to be reliable. We create a rest-frame NIR selected catalog of galaxies at z~2 and construct a galaxy stellar mass function (SMF). Comparisons with other SMFs at approximately the same redshift but determined using shorter wavelengths show good agreement. This agreement suggests that selection at bluer wavelengths does not miss a significant amount of stellar mass in passive galaxies. Comparison with SMFs at other redshifts shows evidence for the downsizing scenario of galaxy evolution. We conclude by pointing out the potential for using the 1.6um technique to select high-redshift galaxies with the JWST, whose lambda > 0.6 um coverage will not be well suited to selecting galaxies using techniques that require imaging at shorter wavelengths.