We employ the hydrodynamical simulation IllustrisTNG-300-1 to explore the
halo occupation distribution (HOD) and environmental dependence of luminous
star-forming emission-line galaxies (ELGs) at $z \sim 1$. Such galaxies are key
targets for current and upcoming cosmological surveys. We select model galaxies
through cuts in color-color space allowing for a direct comparison with the
Extended Baryon Oscillation Spectroscopic Survey and the Dark Energy
Spectroscopic Instrument (DESI) surveys and then compare them with galaxies
selected based on star-formation rate (SFR) and stellar mass. We demonstrate
that the ELG populations are twice more likely to reside in lower-density
regions (sheets) compared with the mass-selected populations and twice less
likely to occupy the densest regions of the cosmic web (knots). We also show
that the color-selected and SFR-selected ELGs exhibit very similar occupation
and clustering statistics, finding that the agreement is best for lower
redshifts. In contrast with the mass-selected sample, the occupation of halos
by a central ELG peaks at $\sim$20\%. We furthermore explore the dependence of
the HOD and the auto-correlation on environment, noticing that at fixed halo
mass, galaxies in high-density regions cluster about 10 times more strongly
than low-density ones. This result suggests that we should model carefully the
galaxy-halo relation and implement assembly bias effects into our models
(estimated at $\sim$4\% of the clustering of the DESI color-selected sample at
$z = 0.8$). Finally, we apply a simple mock recipe to recover the clustering on
large scales ($r \gtrsim 1 \ {\rm Mpc}/h$) to within 1\% by augmenting the HOD
model with an environment dependence, demonstrating the power of adopting
flexible population models.