A new methodology has been applied to the experimental data obtained about a white siliceous conglomerate from Zamora (Spain), which was subjected to 25 cycles of 2 types of aging [freezing/thawing with cooling/heating (T1) and freezing/thawing with cooling/heating + phosphate crystallization (T2)]. Our model (multivariate Gaussian subspatial regression) allows the behavior and prediction of the chromatic coordinates (L*,a*,b*), including more than 25 cycles, to be analyzed. This model is much more flexible than classical models as it allows multiple variable combinations to be predicted in a dynamic way. The final result showed that the conglomerate experiences darkening, yellowing, and reddening, as the number of cycles increase and that the darkening is much less pronounced in T2 due to phosphate crystallization.