The James Webb Space Telescope will allow to spectroscopically study an unprecedented number of galaxies deep into the reionization era, notably by detecting [OIII] and H$\beta$ nebular emission lines. To efficiently prepare such observations, we photometrically select a large sample of galaxies at $z\sim8$ and study their rest-frame optical emission lines. Combining data from the GOODS Re-ionization Era wide-Area Treasury from Spitzer (GREATS) survey and from HST, we perform spectral energy distribution (SED) fitting, using synthetic SEDs from a large grid of photoionization models. The deep Spitzer/IRAC data combined with our models exploring a large parameter space enables to constrain the [OIII]+H$\beta$ fluxes and equivalent widths for our sample, as well as the average physical properties of $z\sim8$ galaxies, such as the ionizing photon production efficiency with $\log(\xi_\mathrm{ion}/\mathrm{erg}^{-1}\hspace{1mm}\mathrm{Hz})\geq25.77$. We find a relatively tight correlation between the [OIII]+H$\beta$ and UV luminosity, which we use to derive for the first time the [OIII]+H$\beta$ luminosity function (LF) at $z\sim8$. The $z\sim8$ [OIII]+H$\beta$ LF is higher at all luminosities compared to lower redshift, as opposed to the UV LF, due to an increase of the [OIII]+H$\beta$ luminosity at a given UV luminosity from $z\sim3$ to $z\sim8$. Finally, using the [OIII]+H$\beta$ LF, we make predictions for JWST/NIRSpec number counts of $z\sim8$ galaxies. We find that the current wide-area extragalactic legacy fields are too shallow to use JWST at maximal efficiency for $z\sim8$ spectroscopy even at 1hr depth and JWST pre-imaging to $\gtrsim30$ mag will be required.

We present an exhaustive census of Lyman alpha (Ly$\alpha$) emission in the general galaxy population at $3<z<4.6$. We use the Michigan/Magellan Fiber System (M2FS) spectrograph to study a stellar mass (M$_*$) selected sample of 625 galaxies homogeneously distributed in the range $7.6<\log{\mbox{M$_*$/M$_{\odot}$}}<10.6$. Our sample is selected from the 3D-HST/CANDELS survey, which provides the complementary data to estimate Ly$\alpha$ equivalent widths ($W_{Ly\alpha}$) and escape fractions ($f_{esc}$) for our galaxies. We find both quantities to anti-correlate with M$_*$, star-formation rate (SFR), UV luminosity, and UV slope ($\beta$). We then model the $W_{Ly\alpha}$ distribution as a function of M$_{UV}$ and $\beta$ using a Bayesian approach. Based on our model and matching the properties of typical Lyman break galaxy (LBG) selections, we conclude that the $W_{Ly\alpha}$ distribution in such samples is heavily dependent on the limiting M$_{UV}$ of the survey. Regarding narrowband surveys, we find their $W_{Ly\alpha}$ selections to bias samples toward low M$_*$, while their line-flux limitations preferentially leave out low-SFR galaxies. We can also use our model to predict the fraction of Ly$\alpha$-emitting LBGs at $4\leqslant z\leqslant 7$. We show that reported drops in the Ly$\alpha$ fraction at $z\geqslant6$, usually attributed to the rapidly increasing neutral gas fraction of the universe, can also be explained by survey M$_{UV}$ incompleteness. This result does not dismiss reionization occurring at $z\sim7$, but highlights that current data is not inconsistent with this process taking place at $z>7$.

In this work, we show how the stellar mass (M) of galaxies affects the 3<z<4.6 Ly-alpha equivalent width (EW) distribution. To this end, we design a sample of 629 galaxies in the M range 7.6 < logM/Msun < 10.6 from the 3D-HST/CANDELS survey. We perform spectroscopic observations of this sample using the Michigan/Magellan Fiber System, allowing us to measure Ly-alpha fluxes and use 3D-HST/CANDELS ancillary data. In order to study the Ly-alpha EW distribution dependence on M, we split the whole sample in three stellar mass bins. We find that, in all bins, the distribution is best represented by an exponential profile of the form dN(M)/dEW= A(M)exp(-EW/W0(M))/W0(M). Through a Bayesian analysis, we confirm that lower M galaxies have higher Ly-alpha EWs. We also find that the fraction A of galaxies featuring emission and the e-folding scale W0 of the distribution anti- correlate with M, recovering expressions of the forms A(M)= -0.26(.13) logM/Msun+3.01(1.2) and W0(M)= -15.6(3.5) logM/Msun +166(34). These results are crucial for proper interpretation of Ly-alpha emission trends reported in the literature that may be affected by strong M selection biases.