The upcoming $\textit{James Webb Space Telescope}$ ($\textit{JWST}$) will allow observations of high-redshift galaxies at fainter detection levels than ever before, and $\textit{JWST}$ surveys targeting gravitationally lensed fields are expected to bring $z\gtrsim 6$ objects with very low star formation rate (SFR) within reach of spectroscopic studies. As galaxies at lower and lower star formation activity are brought into view, many of the standard methods used in the analysis of integrated galaxy spectra are at some point bound to break down, due to violation of the assumptions of a well-sampled stellar initial mass function (IMF) and a slowly varying SFR. We argue that galaxies with SFR$\sim 0.1\ M_\odot$ yr$^{-1}$ are likely to turn up at the spectroscopic detection limit of $\textit{JWST}$ in lensed fields, and investigate to what extent star formation sampling may affect the spectral analysis of such objects. We use the $\small{\text{SLUG}}$ spectral synthesis code to demonstrate that such effects are likely to have significant impacts on spectral diagnostics of, for example, the Balmer emission lines. These effects are found to stem primarily from SFRs varying rapidly on short ($\sim$ Myr) timescales due to star formation in finite units (star clusters), whereas the effects of an undersampled IMF is deemed insignificant in comparison. In contrast, the ratio between the HeII- and HI-ionizing flux is found to be sensitive to IMF-sampling as well as ICMF-sampling (sampling of the initial cluster mass function), which may affect interpretations of galaxies containing Population III stars or other sources of hard ionizing radiation.

Photometric observations of the spectroscopically confirmed $z\approx 9.1$ galaxy MACS1149-JD1 have indicated the presence of a prominent Balmer break in its spectral energy distribution, which may be interpreted as due to very large fluctuations in its past star formation activity. In this paper, we investigate to what extent contemporary simulations of high-redshift galaxies produce star formation rate variations sufficiently large to reproduce the observed Balmer break of MACS1149-JD1. We find that several independent galaxy simulations are unable to account for Balmer breaks of the inferred size, suggesting that MACS1149-JD1 either must be a very rare type of object or that our simulations are missing some key ingredient. We present predictions of spectroscopic Balmer break strength distributions for $z\approx 7-9$ galaxies that may be tested through observations with the upcoming James Webb Space Telescope and also discuss the impact that various assumptions on dust reddening, Lyman continuum leakage and deviations from a standard stellar initial mass function would have on the results.

New challenges in the modelling of galaxy spectra are bound to emerge as upcoming telescopes like the James Webb Space Telescope will allow us to detect galaxies at fainter flux levels and higher redshifts than ever before. Here, we highlight three modelling problems that may become relevant for upcoming observations in the rest-frame ultraviolet/optical at the high-redshift frontier: stellar initial mass function sampling effects, Population III signatures and the leakage of ionizing radiation into the intergalactic medium.