The Bayesian inference approach to the inverse problem of spectral signal recovery has been extended to mixtures of Gaussian probability distributions of a training dataset in order to increase the efficiency of estimating the spectral signal from the response of a transformation system. Bayesian (BIC) and Akaike (AIC) information criteria were assessed in order to provide the Gaussian mixture model (GMM) with the optimum number of clusters within the spectral space. The spectra of 2600 solar illuminations measured in Granada (Spain) were recovered over the range of 360-830 nm from their corresponding tristimulus values using a linear model of basis functions, the Wiener inverse (WI) method, and the Bayesian inverse approach extended to the GMM (BGMM). A model of Gaussian mixtures for solar irradiance was deemed to be more appropriate than a single Gaussian distribution for representing the probability distribution of the solar spectral data. The results showed that the estimation performance of the BGMM method was better than either the linear model or the WI method for the spectral approximation of daylight from the three-dimensional tristimulus values.