Despite recent progress in understanding Lyman-alpha (Lya) emitters (LAEs), relatively little is known regarding their typical black hole activity across cosmic time. Here, we study the X-ray and radio properties of ~4000 LAEs at 2.2<z<6 from the SC4K survey in the COSMOS field. We detect 254 (6.8+-0.4%) of the LAEs individually in the X-rays (S/N>3) and find an average luminosity of 10^{44.31+-0.01} erg/s and an average black hole accretion rate (BHAR) of 0.72+-0.01 Msun/yr, consistent with moderate to high accreting AGN. We detect 120 sources in deep radio data (radio AGN fraction of 3.2%+-0.3%). The global AGN fraction (8.6%+-0.4%) rises with Lya luminosity and declines with increasing redshift. For X-ray detected LAEs, Lya luminosities correlate with the BHARs, suggesting that Lya luminosity becomes an accretion rate indicator. Most LAEs (93.1+-0.6%) at 2<z<6 have no detectable X-ray emission (BHARs<0.017 Msun/yr). The median star formation rate (SFR) of star-forming LAEs from Lya and radio luminosities is 7.6^{+6.6}_{-2.8} Msun/yr. The black hole to galaxy growth ratio (BHAR/SFR) for LAEs is <0.0022, consistent with typical star forming galaxies and the local BHAR/SFR relation. We conclude that LAEs at 2<z<6 include two different populations: an AGN population, where Lya luminosity traces BHAR, making them bright in Lya, and another with low SFRs which remain undetected in even the deepest X-ray stacks but is detected in the radio stacks.

We explore deep rest-frame UV to FIR data in the COSMOS field to measure the individual spectral energy distributions (SED) of the ~4000 SC4K (Sobral et al. 2018) Lyman-alpha (Lya) emitters (LAEs) at z~2-6. We find typical stellar masses of 10$^{9.3\pm0.6}$ M$_{\odot}$ and star formation rates (SFR) of SFR$_{SED}=4.4^{+10.5}_{-2.4}$ M$_{\odot}$/yr and SFR$_{Lya}=5.9^{+6.3}_{-2.6}$ M$_{\odot}$/yr, combined with very blue UV slopes of beta=-2.1$^{+0.5}_{-0.4}$, but with significant variations within the population. M$_{UV}$ and beta are correlated in a similar way to UV-selected sources, but LAEs are consistently bluer. This suggests that LAEs are the youngest and/or most dust-poor subset of the UV-selected population. We also study the Lya rest-frame equivalent width (EW$_0$) and find 45 "extreme" LAEs with EW$_0>240$ A (3 $\sigma$), implying a low number density of $(7\pm1)\times10^{-7}$ Mpc$^{-3}$. Overall, we measure little to no evolution of the Lya EW$_0$ and scale length parameter ($w_0$) which are consistently high (EW$_0=140^{+280}_{-70}$ A, $w_0=129^{+11}_{-11}$ A) from z~6 to z~2 and below. However, $w_0$ is anti-correlated with M$_{UV}$ and stellar mass. Our results imply that sources selected as LAEs have a high Lya escape fraction (f$_{esc, Lya}$) irrespective of cosmic time, but f$_{esc, Lya}$ is still higher for UV-fainter and lower mass LAEs. The least massive LAEs ($<10^{9.5}$ M$_{\odot}$) are typically located above the star formation "Main Sequence" (MS), but the offset from the MS decreases towards z~6 and towards $10^{10}$ M$_{\odot}$. Our results imply a lack of evolution in the properties of LAEs across time and reveals the increasing overlap in properties of LAEs and UV-continuum selected galaxies as typical star-forming galaxies at high redshift effectively become LAEs.

We investigate the clustering and halo properties of $\sim 5000$ Ly$\alpha$-selected emission line galaxies (LAEs) from the Slicing COSMOS 4K (SC4K) and from archival NB497 imaging of SA22 split in 15 discrete redshift slices between $z \sim 2.5 - 6$. We measure clustering lengths of $r_0 \sim 3 - 6\ h^{-1}$ Mpc and typical halo masses of $\sim 10^{11}$ M$_\odot$ for our narrowband-selected LAEs with typical $L_{\rm{Ly}\alpha} \sim 10^{42 - 43}$ erg s$^{-1}$. The intermediate band-selected LAEs are observed to have $r_0 \sim 3.5 - 15\ h^{-1}$ Mpc with typical halo masses of $\sim 10^{11 - 12}$ M$_\odot$ and typical $L_{\rm{Ly}\alpha} \sim 10^{43 - 43.6}$ erg s$^{-1}$. We find a strong, redshift-independent correlation between halo mass and Ly$\alpha$ luminosity normalized by the characteristic Ly$\alpha$ luminosity, $L^\star(z)$. The faintest LAEs ($L \sim 0.1\ L^\star(z)$) typically identified by deep narrowband surveys are found in $10^{10}$ M$_\odot$ halos and the brightest LAEs ($L \sim 7\ L^\star(z)$) are found in $\sim 5 \times 10^{12}$ M$_\odot$ halos. A dependency on the rest-frame 1500 \AA~UV luminosity, M$_\rm{UV}$, is also observed where the halo masses increase from $10^{11}$ to $10^{13}$ M$_\odot$ for M$_\rm{UV} \sim -19$ to $-23.5$ mag. Halo mass is also observed to increase from $10^{9.8}$ to $10^{12.3}$ M$_\odot$ for dust-corrected UV star formation rates from $\sim 0.6$ to $10$ M$_\odot$ yr$^{-1}$ and continues to increase up to $10^{13.5}$ M$_\odot$ in halo mass, where the majority of those sources are AGN. All the trends we observe are found to be redshift-independent. Our results reveal that LAEs are the likely progenitors of a wide range of galaxies depending on their luminosity, from dwarf-like, to Milky Way-type, to bright cluster galaxies. LAEs therefore provide unique insight into the early formation and evolution of the galaxies we observe in the local Universe.

We present the second data release of the Large Early Galaxy Astrophysics Census (LEGA-C), an ESO 130-night public spectroscopic survey conducted with VIMOS on the Very Large Telescope. We release 1988 spectra with typical continuum S / N ~= 20 /Angstrom of galaxies at 0.6 ~< z ~< 1.0, each observed for ~20 hours and fully reduced with a custom-built pipeline. We also release a catalog with spectroscopic redshifts, emission line fluxes, Lick/IDS indices, and observed stellar and gas velocity dispersions that are spatially integrated quantities including both rotational motions and genuine dispersion. To illustrate the new parameter space in the intermediate redshift regime probed by LEGA-C we explore relationships between dynamical and stellar population properties. The star-forming galaxies typically have observed stellar velocity dispersions of ~150 km/s and strong Hdelta absorption (Hd_A ~ 5 Angstrom), while passive galaxies have higher observed stellar velocity dispersions ~200 km/s and weak Hdelta absortion (Hd_A ~ 0 Angstrom). Strong [O III]5007 / Hbeta ratios tend to occur mostly for galaxies with weak Hd_A or galaxies with higher observed velocity dispersion. Beyond these broad trends, we find a large diversity of possible combinations of rest-frame colors, absorption line strengths and emission line detections, illustrating the utility of spectroscopic measurements to more accurately understand galaxy evolution. By making the spectra and value-added catalogs publicly available we encourage the community to take advantage of this very substantial investment in telescope time provided by ESO.

Using high resolution spectra from the VLT LEGA-C program, we reconstruct the star formation histories (SFHs) of 607 galaxies at redshifts $z = 0.6-1.0$ and stellar masses $\gtrsim10^{10}$M$_{\odot}$ using a custom full spectrum fitting algorithm that incorporates the emcee and FSPS packages. We show that the mass-weighted age of a galaxy correlates strongly with stellar velocity dispersion ($\sigma_*$) and ongoing star-formation (SF) activity, with the stellar content in higher-$\sigma_*$ galaxies having formed earlier and faster. The SFHs of quiescent galaxies are generally consistent with passive evolution since their main SF epoch, but a minority show clear evidence of a rejuvenation event in their recent past. The mean age of stars in galaxies that are star-forming is generally significantly younger, with SF peaking after $z<1.5$ for almost all star-forming galaxies in the sample: many of these still have either constant or rising SFRs on timescales $>100$Myrs. This indicates that $z>2$ progenitors of $z\sim1$ star-forming galaxies are generally far less massive. Finally, despite considerable variance in the individual SFHs, we show that the current SF activity of massive galaxies ($>$L$_*$) at $z\sim1$ correlates with SF levels at least $3$Gyrs prior: SFHs retain `memory' on a large fraction of the Hubble time. Our results illustrate a novel approach to resolve the formation phase of galaxies, and in identifying their individual evolutionary paths, connects progenitors and descendants across cosmic time. This is uniquely enabled by the high-quality continuum spectroscopy provided by the LEGA-C survey.

We present stellar rotation curves and velocity dispersion profiles for 104 quiescent galaxies at $z=0.6-1$ from the Large Early Galaxy Astrophysics Census (LEGA-C) spectroscopic survey. Rotation is typically probed across 10-20kpc, or to an average of 2.7${\rm R_e}$. Combined with central stellar velocity dispersions ($\sigma_0$) this provides the first determination of the dynamical state of a sample selected by a lack of star formation activity at large lookback time. The most massive galaxies ($M_{\star}>2\times10^{11}\,M_{\odot}$) generally show no or little rotation measured at 5kpc ($|V_5|/\sigma_0<0.2$ in 8 of 10 cases), while ${\sim}64\%$ of less massive galaxies show significant rotation. This is reminiscent of local fast- and slow-rotating ellipticals and implies that low- and high-redshift quiescent galaxies have qualitatively similar dynamical structures. We compare $|V_5|/\sigma_0$ distributions at $z\sim0.8$ and the present day by re-binning and smoothing the kinematic maps of 91 low-redshift quiescent galaxies from the CALIFA survey and find evidence for a decrease in rotational support since $z\sim1$. This result is especially strong when galaxies are compared at fixed velocity dispersion; if velocity dispersion does not evolve for individual galaxies then the rotational velocity at 5kpc was an average of ${94\pm22\%}$ higher in $z\sim0.8$ quiescent galaxies than today. Considering that the number of quiescent galaxies grows with time and that new additions to the population descend from rotationally-supported star-forming galaxies, our results imply that quiescent galaxies must lose angular momentum between $z\sim1$ and the present, presumably through dissipationless merging, and/or that the mechanism that transforms star-forming galaxies also reduces their rotational support.

We investigate the rest-frame UV morphologies of a large sample of Lyman-a emitters (LAEs) from z~2 to z~6, selected in a uniform way with 16 different narrow- and medium-bands over the full COSMOS field. We use 3045 LAEs with HST coverage in a stacking analysis and find that they have M_UV~-20, below M*_UV at these redshifts. We also focus our analysis on a subsample of 780 individual galaxies with i_AB<25 for which GALFIT converges for 429 of them. The individual median size (re~1 kpc), ellipticities (slightly elongated with (b/a)~0.45), S\'ersic index (disk-like with n<2) and light concentration (comparable to that of disk or irregular galaxies, with C~2.7) of LAEs show mild evolution from z~2 to z~6. LAEs with the highest rest-frame equivalent widths (EW) are the smallest/most compact (re~0.8 kpc, compared to re~1.5 kpc for the lower EW LAEs). When stacking our samples in bins of fixed Lya luminosity and Lya EW we find evidence for redshift evolution in n and C, but not in galaxy sizes. The evolution seems to be stronger for LAEs with 25<EW<100 {\AA}. When compared to other SFGs, LAEs are found to be smaller at all redshifts. The difference between the two populations changes with redshift, from a factor of ~1 at z>5 to SFGs being a factor of ~2-4 larger than LAEs for z<2. This means that at the highest redshifts, where typical sizes approach those of LAEs, the fraction of galaxies showing Lya in emission (and with a high Lya escape fraction) should be much higher, consistent with observations.

Drawing from the LEGA-C dataset, we present the spectroscopic view of the stellar population across a large volume- and mass-selected sample of galaxies at large lookback time. We measure the 4000\AA\ break (D$_n$4000) and Balmer absorption line strengths (probed by H$\delta$) from 1019 high-quality spectra of $z=0.6 - 1.0$ galaxies with $M_\ast = 2 \times 10^{10} M_\odot - 3 \times 10^{11} M_\odot$. Our analysis serves as a first illustration of the power of high-resolution, high-S/N continuum spectroscopy at intermediate redshifts as a qualitatively new tool to constrain galaxy formation models. The observed D$_n$4000-EW(H$\delta$) distribution of our sample overlaps with the distribution traced by present-day galaxies, but $z\sim 0.8$ galaxies populate that locus in a fundamentally different manner. While old galaxies dominate the present-day population at all stellar masses $> 2\times10^{10} M_\odot$, we see a bimodal D$_n$4000-EW(H$\delta$) distribution at $z\sim0.8$, implying a bimodal light-weighted age distribution. The light-weighted age depends strongly on stellar mass, with the most massive galaxies $>1\times10^{11}M_\odot$ being almost all older than 2 Gyr. At the same time we estimate that galaxies in this high mass range are only $\sim3$ Gyr younger than their $z\sim0.1$ counterparts, at odd with pure passive evolution given a difference in lookback time of $>5$ Gyr; younger galaxies must grow to $>10^{11}M_\odot$ in the meantime, and/or small amounts of young stars must keep the light-weighted ages young. Star-forming galaxies at $z\sim0.8$ have stronger H$\delta$ absorption than present-day galaxies with the same D$_n$4000, implying larger short-term variations in star-formation activity.

We present and explore deep narrow- and medium-band data obtained with the Subaru and the Isaac Newton telescopes in the ~2 deg$^2$ COSMOS field. We use these data as an extremely wide, low-resolution (R~20-80) IFU survey to slice through the COSMOS field and obtain a large sample of ~4000 Lyman-$\alpha$ (Lya) emitters from z~2 to z~6 in 16 redshift slices (SC4K). We present new Lya luminosity functions (LFs) covering a co-moving volume of ~10$^8$Mpc$^3$. SC4K extensively complements ultra-deep surveys, jointly covering over 4 dex in Lya luminosity and revealing a global (2.5<z<6) synergy LF with $\alpha=-1.93\pm0.12$, $\log_{10}\Phi^*=-3.45^{+0.22}_{-0.29}$ Mpc$^{-3}$ and $\log_{10}L^*=42.93^{+0.15}_{-0.11}$ erg/s. The Schechter component of the Lya LF reveals a factor ~5 rise in $L^*$ and a ~7x decline in $\Phi^*$ from z~2 to z~6. The data reveal an extra power-law (or Schechter) component above L~10$^{43.3}$ erg/s at z~2.2-3.5 and we show that it is partially driven by X-ray and radio AGN, as their Lya LF resembles the excess. The power-law component vanishes and/or is below our detection limits above z>3.5, likely linked with the evolution of the AGN population. The Lya luminosity density rises by a factor ~2 from z~2 to z~3 but is then found to be roughly constant (~$1.1\times10^{40}$ erg s$^{-1}$ Mpc$^{-3}$) to z~6, despite the ~0.7 dex drop in UV luminosity density. The Lya/UV luminosity density ratio rises from $4\pm1$% to $30\pm6$% from z~2.2 to z~6. Our results imply a rise of a factor of ~2 in the global ionisation efficiency ($\xi_{\rm ion}$) and a factor ~$4\pm1$ in the Lya escape fraction from z~2 to z~6, hinting for evolution in both the typical burstiness/stellar populations and even more so in the typical ISM conditions allowing Ly$\alpha$ photons to escape.

We investigate the stellar kinematics and stellar populations of 58 radio-loud galaxies of intermediate luminosities (L$_{3 GHz}$ $>$ 10$^{23}$ W Hz$^{-1}$ ) at 0.6 < z < 1. This sample is constructed by cross-matching galaxies from the deep VLT/VIMOS LEGA-C spectroscopic survey with the VLA 3 GHz dataset. The LEGA-C continuum spectra reveal for the first time stellar velocity dispersions and age indicators of z $\sim$ 1 radio galaxies. We find that $z\sim 1$ radio-loud AGN occur exclusively in predominantly old galaxies with high velocity dispersions: $\sigma_*>$ 175 km s$^{-1}$, corresponding to black hole masses in excess of $10^8$ M$_{\odot}$. Furthermore, we confirm that at a fixed stellar mass the fraction of radio-loud AGN at z $\sim$ 1 is 5 - 10 times higher than in the local universe, suggesting that quiescent, massive galaxies at z $\sim$ 1 switch on as radio AGN on average once every Gyr. Our results strengthen the existing evidence for a link between high black-hole masses, radio loudness and quiescence at z $\sim$ 1.

While traditionally associated with active galactic nuclei (AGN), the properties of the CII] ($\lambda=2326$\,{\AA}), CIII] ($\lambda,\lambda=1907,1909$\,{\AA}) and CIV ($\lambda,\lambda=1549, 1551$\,{\AA}) emission lines are still uncertain as large, unbiased samples of sources are scarce. We present the first blind, statistical study of CII], CIII] and CIV emitters at $z\sim0.68,1.05,1.53$, respectively, uniformly selected down to a flux limit of $\sim4\times10^{-17}$ erg s$^{-1}$ cm$^{-1}$ through a narrow band survey covering an area of $\sim1.4$ deg$^2$ over COSMOS and UDS. We detect 16 CII], 35 CIII] and 17 CIV emitters, whose nature we investigate using optical colours as well as HST, X-ray, radio and far infra-red data. We find that $z\sim0.7$ CII] emitters are consistent with a mixture of blue (UV slope $\beta=-2.0\pm0.4$) star forming galaxies with disky HST structure and AGN with Seyfert-like morphologies. Bright CII] emitters have individual X-ray detections as well as high average black hole accretion rates (BHAR) of $\sim0.1$ $M_{\odot}$ yr$^{-1}$. CIII] emitters at $z\sim1.05$ trace a general population of SF galaxies, with $\beta=-0.8\pm1.1$, a variety of optical morphologies, including isolated and interacting galaxies and low BHAR ($<0.02$ $M_{\odot}$ yr$^{-1}$). Our CIV emitters at $z\sim1.5$ are consistent with young, blue quasars ($\beta\sim-1.9$) with point-like optical morphologies, bright X-ray counterparts and large BHAR ($0.8$ $M_{\odot}$ yr$^{-1}$). We also find some surprising CII], CIII] and CIV emitters with rest-frame equivalent widths which could be as large as $50-100$ {\AA}. AGN or spatial offsets between the UV continuum stellar disk and the line emitting regions may explain the large EW. These bright CII], CIII] and CIV emitters are ideal candidates for spectroscopic follow up to fully unveil their nature.

Recently, the CIII] and CIV emission lines have been observed in galaxies in the early Universe ($z>5$), providing new ways to measure their redshift and study their stellar populations and AGN. We explore the first blind CII], CIII] and CIV survey ($z\sim0.68, 1.05, 1.53$, respectively) presented in Stroe et al. (2017). We derive luminosity functions (LF) and study properties of CII], CIII] and CIV line emitters through comparisons to the LFs of H$\alpha$ and Ly$\alpha$ emitters, UV selected star forming (SF) galaxies and quasars at similar redshifts. The CII] LF at $z\sim0.68$ is equally well described by a Schechter or a power-law LF, characteristic of a mixture of SF and AGN activity. The CIII] LF ($z\sim1.05$) is consistent to a scaled down version of the Schechter H$\alpha$ and Ly$\alpha$ LF at their redshift, indicating a SF origin. In stark contrast, the CIV LF at $z\sim1.53$ is well fit by a power-law, quasar-like LF. We find that the brightest UV sources ($M_{UV}<-22$) will universally have CIII] and CIV emission. However, on average, CIII] and CIV are not as abundant as H$\alpha$ or Ly$\alpha$ emitters at the same redshift, with cosmic average ratios of $\sim0.02-0.06$ to H$\alpha$ and $\sim0.01-0.1$ to intrinsic Ly$\alpha$. We predict that the CIII] and CIV lines can only be truly competitive in confirming high redshift candidates when the hosts are intrinsically bright and the effective Ly$\alpha$ escape fraction is below 1 per cent. While CIII] and CIV were proposed as good tracers of young, relatively low-metallicity galaxies typical of the early Universe, we find that, at least at $z\sim1.5$, CIV is exclusively hosted by AGN/quasars, especially at large line equivalent widths.

We present the first results from the largest H$\alpha$ survey of star formation and AGN activity in galaxy clusters. Using 9 different narrow band filters, we select $>3000$ H$\alpha$ emitters within $19$ clusters and their larger scale environment over a total volume of $1.3\times10^5$ Mpc$^3$. The sample includes both relaxed and merging clusters, covering the $0.15-0.31$ redshift range and spanning from $5\times10^{14}$ $M_{\odot}$ to $30\times10^{14}$ $M_{\odot}$. We find that the H$\alpha$ luminosity function (LF) for merging clusters has a higher characteristic density $\phi^*$ compared to relaxed clusters. $\phi^*$ drops from cluster core to cluster outskirts for both merging and relaxed clusters, with the merging cluster values $\sim0.3$ dex higher at each projected radius. The characteristic luminosity $L^*$ drops over the $0.5-2.0$ Mpc distance from the cluster centre for merging clusters and increases for relaxed objects. Among disturbed objects, clusters hosting large-scale shock waves (traced by radio relics) are overdense in H$\alpha$ emitters compared to those with turbulence in their intra-cluster medium (traced by radio haloes). We speculate that the increase in star formation activity in disturbed, young, massive galaxy clusters can be triggered by interactions between gas-rich galaxies, shocks and/or the intra-cluster medium, as well as accretion of filaments and galaxy groups. Our results indicate that disturbed clusters represent vastly different environments for galaxy evolution compared to relaxed clusters or average field environments.

Understanding galaxy formation and evolution requires studying the interplay between the growth of galaxies and the growth of their black holes across cosmic time. Here we explore a sample of Ha-selected star-forming galaxies from the HiZELS survey and use the wealth of multi-wavelength data in the COSMOS field (X-rays, far-infrared and radio) to study the relative growth rates between typical galaxies and their central supermassive black holes, from z=2.23 to z=0. Typical star-forming galaxies at z~1-2 have black hole accretion rates (BHARs) of 0.001-0.01 Msun/yr and star formation rates (SFRs) of ~10-40 Msun/yr, and thus grow their stellar mass much quicker than their black hole mass (~3.3 orders of magnitude faster). However, ~3% of the sample (the sources detected directly in the X-rays) show a significantly quicker growth of the black hole mass (up to 1.5 orders of magnitude quicker growth than the typical sources). BHARs fall from z=2.23 to z=0, with the decline resembling that of star formation rate density or the typical SFR. We find that the average black hole to galaxy growth (BHAR/SFR) is approximately constant for star-forming galaxies in the last 11 Gyrs. The relatively constant BHAR/SFR suggests that these two quantities evolve equivalently through cosmic time and with practically no delay between the two.

Cluster star-forming galaxies are found to have an excess of Far-Infrared emission relative to H-alpha (Ha), when compared to those in the field, which could be caused by intense AGN activity, dust and/or declining star formation histories. Here we present spectroscopic observations of Ha emitters in the Cl 0939+4713 (Abell 851) super-cluster at z=0.41, using AF2+WYFFOS on the WHT. We measure [OII], Hbeta (Hb), [OIII], Ha and [NII] for a sample of 119 Ha emitters in and around the cluster. We find that 17+-5% of the Ha emitters are AGN, irrespective of environment. For star-forming galaxies, we obtain Balmer decrements, metallicities and ionisation parameters with different methods, individually and by stacking. We find a strong mass-metallicity relation at all environments, with no significant dependence on environment. The ionisation parameter declines with increasing stellar mass for low-mass galaxies. Ha emitters residing in intermediate environments show the highest ionisation parameters (along with high [OIII]/Ha and high [OIII]/[OII] line ratios, typically twice as large as in the highest and lowest densities), which decline with increasing environmental density. Dust extinction (A$_{H\alpha}$) correlates strongly with stellar mass, but also with environmental density. Star-forming galaxies in the densest environments are found to be significantly dustier (A$_{H\alpha}$~1.5-1.6) than those residing in the lowest density environments (A$_{H\alpha}$~0.6), deviating significantly from what would be predicted given their stellar masses.

We use new near-infrared spectroscopic observations to investigate the nature and evolution of the most luminous H\alpha (Ha) emitters at z~0.8-2.23, which evolve strongly in number density over this period, and compare them to more typical Ha emitters. We study 59 luminous Ha emitters with $L_{H\alpha}>L_{H\alpha}^*$, roughly equally split per redshift slice at z~0.8, 1.47 and 2.23 from the HiZELS and CF-HiZELS surveys. We find that, overall, 30$\pm$8% are AGN (80$\pm$30% of these AGN are broad-line AGN, BL-AGN), and we find little to no evolution in the AGN fraction with redshift, within the errors. However, the AGN fraction increases strongly with Ha luminosity and correlates best with $L_{H\alpha}/L_{H\alpha}^*(z)$. While $L_{H\alpha}<L_{\rm H\alpha}^*(z)$ Ha emitters are largely dominated by star-forming galaxies (>80%), the most luminous Ha emitters ($L_{H\alpha}>10L_{H\alpha}^*(z)$) at any cosmic time are essentially all BL-AGN. Using our AGN-decontaminated sample of luminous star-forming galaxies, and integrating down to a fixed Ha luminosity, we find a factor of ~1300x evolution in the star formation rate density from z=0 to z=2.23. This is much stronger than the evolution from typical Ha star-forming galaxies and in line with the evolution seen for constant luminosity cuts used to select "Ultra-Luminous" Infrared Galaxies and/or sub-millimetre galaxies. By taking into account the evolution in the typical Ha luminosity, we show that the most strongly star-forming Ha-selected galaxies at any epoch ($L_{H\alpha}>L^*_{H\alpha}(z)$) contribute the same fractional amount of ~15% to the total star-formation rate density, at least up to z=2.23.

We present results from the largest contiguous narrow-band survey in the near-infrared. We have used WIRCam/CFHT and the lowOH2 filter (1.187$\pm$0.005 $\mu$m) to survey ~10 deg$^2$ of contiguous extragalactic sky in the SA22 field. A total of ~6000 candidate emission-line galaxies are found. We use deep ugrizJK data to obtain robust photometric redshifts. We combine our data with the High-redshift Emission Line Survey (HiZELS), explore spectroscopic surveys (VVDS, VIPERS) and obtain our own spectroscopic follow-up with KMOS, FMOS and MOSFIRE to derive large samples of high-redshift emission-line selected galaxies: 3471 H$\alpha$ emitters at z=0.8, 1343 [OIII]+H$\beta$ emitters at z=1.4 and 572 [OII] emitters at z=2.2. We probe co-moving volumes of >10$^6$ Mpc$^3$ and find significant over-densities, including an 8.5$\sigma$ (spectroscopically confirmed) over-density of H$\alpha$ emitters at z=0.81. We derive H$\alpha$, [OIII]+H$\beta$ and [OII] luminosity functions at z=0.8,1.4,2.2, respectively, and present implications for future surveys such as Euclid. Our uniquely large volumes/areas allow us to sub-divide the samples in thousands of randomised combinations of areas and provide a robust empirical measurement of sample/cosmic variance. We show that surveys for star-forming/emission-line galaxies at a depth similar to ours can only overcome cosmic-variance (errors <10%) if they are based on volumes >5x10$^{5}$ Mpc$^{3}$; errors on $L^*$ and $\phi^*$ due to sample (cosmic) variance on surveys probing ~10$^4$ Mpc$^{3}$ and ~10$^5$ Mpc$^{3}$ are typically very high: ~300% and ~40-60%, respectively.