Sampling fluctuations in stellar populations give rise to dispersions in observables when a small number of sources contribute effectively to the observables. This is the case for a variety of linear functions of the spectral energy distribution (SED) in small stellar systems, such as galactic and extragalactic HII regions, dwarf galaxies or stellar clusters. In this paper we show that sampling fluctuations also introduce systematic biases and multi-modality in non-linear functions of the SED, such as luminosity ratios, magnitudes and colours. In some cases, the distribution functions of rational and logarithmic quantities are bimodal, hence complicating considerably the interpretation of these quantities in terms of age or evolutionary stages. These biases can be only assessed by Monte Carlo simulations. We find that biases are usually negligible when the effective number of stars, Neff, which contribute to a given observable is larger than 10. Bimodal distributions may appear when Neff is between 10 and 0.1. Predictions from any model of stellar population synthesis become extremely unreliable for small Neff values, providing an operational limit to the applicability of such models for the interpretation of integrated properties of stellar systems. In terms of stellar masses, assuming a Salpeter Initial Mass Function in the range 0.08 -- 120 Mo, Neff=10 corresponds to about 10**5 Mo (although the exact value depends on the age and the observable). This bias may account, at least in part, for claimed variations in the properties of the stellar initial mass function in small systems, and arises from the discrete nature of small stellar populations.

We use the deepest and most complete redshift catalog currently available (the Hubble Deep Field (HDF) North supplemented by new HDF South redshift data) to minimize residuals between photometric and spectroscopic redshift estimates. The good agreement at z_spec < 1.5 shows that model libraries provide a good description of the galaxy population. At z_spec > 2.0, the systematic shift between photometric and spectroscopic redshifts decreases when the modeling of the absorption by the interstellar and intergalactic media is refined. As a result, in the entire redshift range 0 < z < 6, residuals between photometric and spectroscopic redshifts are roughly halved. For objects fainter than the spectroscopic limit, the main source of uncertainty in photometric redshifts is related to photometric errors, and can be assessed with Monte Carlo simulations.

Ranging from track interpolation techniques through model atmospheres to the stochastic nature of the IMF, there are many uncertainties which need to be taken into account when modelling HR diagrams or performing population synthesis, particularly if comparison with actual data is sought. In this paper, we highlight (1) the problem of discontinuities along evolutionary tracks of massive stars (M> 8 Mo), showing that inconsistencies appear in the computation of the corresponding isochrones, and (2) the sampling fluctuations produced by the stochastic nature of the IMF, presenting a statistical formalism to estimate the dispersion in any given observable of a stellar population due to sampling effects which bypasses the need of performing Monte Carlo simulations.