Line intensity mapping is a superb tool to study the collective radiation from early galaxies. However, the method is hampered by the presence of strong foregrounds, mostly produced by low-redshift interloping lines. We present here a general method to overcome this problem which is robust against foreground residual noise and based on the cross-correlation function $\psi_{\alpha L}(r)$ between diffuse line emission and Ly$\alpha$ emitters (LAE). We compute the diffuse line (Ly$\alpha$ is used as an example) emission from galaxies in a $(800{\rm Mpc})^3$ box at $z = 5.7$ and $6.6$. We divide the box in slices and populate them with $14000(5500)$ LAEs at $z = 5.7(6.6)$, considering duty cycles from $10^{-3}$ to $1$. Both the LAE number density and slice volume are consistent with the expected outcome of the Subaru HSC survey. We add gaussian random noise with variance $\sigma_{\rm N}$ up to 100 times the variance of the Ly$\alpha$ emission, $\sigma_\alpha$, to simulate foregrounds and compute $\psi_{\alpha L}(r)$. We find that the signal-to-noise of the observed $\psi_{\alpha L}(r)$ does not change significantly if $\sigma_{\rm N} \le 10 \sigma_\alpha$ and show that in these conditions the mean line intensity, $I_{Ly\alpha}$, can be precisely recovered independently of the LAE duty cycle. Even if $\sigma_{\rm N} = 100 \sigma_\alpha$, $I_\alpha$ can be constrained within a factor $2$. The method works equally well for any other line (e.g. HI 21 cm, [CII], HeII) used for the intensity mapping experiment.

Intensity mapping (IM) is sensitive to the cumulative line emission of galaxies. As such it represents a promising technique for statistical studies of galaxies fainter than the limiting magnitude of traditional galaxy surveys. The strong hydrogen Ly-alpha line is the primary target for such an experiment, as its intensity is linked to star formation activity and the physical state of the interstellar (ISM) and intergalactic (IGM) medium. However, to extract the meaningful information one has to solve the confusion problems caused by interloping lines from foreground galaxies. We discuss here the challenges for a Ly-alpha IM experiment targeting z > 4 sources. We find that the Ly-alpha power spectrum can be in principle easily (marginally) obtained with a 40 cm space telescope in a few days of observing time up to z < 8 (z = 10) assuming that the interloping lines (e.g. H-alpha, [O II], [O III] lines) can be efficiently removed. We show that interlopers can be removed by using an ancillary photometric galaxy survey with limiting AB mag 26 in the NIR bands (Y, J, H, or K). This would enable detection of the Ly-alpha signal from 5 < z < 9 faint sources. However, if a [C II] IM experiment is feasible, by cross-correlating the Ly-alpha with the [C II] signal the required depth of the galaxy survey can be decreased to AB mag 24. This would bring the detection at reach of future facilities working in close synergy.

We present a study of the cosmological Ly$\alpha$ emission signal at $z > 4$. Our goal is to predict the power spectrum of the spatial fluctuations that could be observed by an intensity mapping survey. The model uses the latest data from the HST legacy fields and the abundance matching technique to associate UV emission and dust properties with the halos, computing the emission from the interstellar medium (ISM) of galaxies and the intergalactic medium (IGM), including the effects of reionization, self-consistently. The Ly$\alpha$ intensity from the diffuse IGM emission is 1.3 (2.0) times more intense than the ISM emission at $z = 4(7)$; both components are fair tracers of the star-forming galaxy distribution. However the power spectrum is dominated by ISM emission on small scales ($k > 0.01 h{\rm Mpc}^{-1}$) with shot noise being significant only above $k = 1 h{\rm Mpc}^{-1}$. At very lange scales ($k < 0.01h{\rm Mpc}^{-1}$) diffuse IGM emission becomes important. The comoving Ly$\alpha$ luminosity density from IGM and galaxies, $\dot \rho_{{\rm Ly}\alpha}^{\rm IGM} = 8.73(6.51) \times 10^{40} {\rm erg}{\rm s}^{-1}{\rm Mpc}^{-3}$ and $\dot \rho_{{\rm Ly}\alpha}^{\rm ISM} = 6.62(3.21) \times 10^{40} {\rm erg}{\rm s}^{-1}{\rm Mpc}^{-3}$ at $z = 4(7)$, is consistent with recent SDSS determinations. We predict a power $k^3 P^{{\rm Ly}\alpha}(k, z)/2\pi^2 = 9.76\times 10^{-4}(2.09\times 10^{-5}){\rm nW}^2{\rm m}^{-4}{\rm sr}^{-2}$ at $z = 4(7)$ for $k = 0.1 h {\rm Mpc}^{-1}$.