Sub-gigayear variability around the star formation main sequence and its contribution to the scatter

Aims. We aim to measure the evolution of individual galaxies around the star formation main sequence (SFMS) during the last gigayear as a function of their stellar mass to quantify how much of its scatter is due to short-term variability. Methods. We derived star formation histories using full spectral fitting for a sample of 8960 galaxies from the MaNGA survey to track the position of the galaxies in the SFMS over the past gigayear. Results. The variability correlates with both the stellar mass of the galaxies and their current position in both the SFMS and the mass-metallicity relation (MZR), with the position in the latter strongly affecting variability in the star formation rate. While most of the fluctuations are compatible with stochasticity, there is a very weak but statistically significant preference for ∼135 − 150 Myr timescales. Conclusions. These results indicate that the star formation rate is strongly self-regulated within galaxies, establishing characteristic intensities and timescales for bursts of star formation and quenching episodes. We also find that short-term variability cannot account for the entirety of the scatter in the SFMS. It appears to originate to a similar degree in short-term variability and long-term (halo-level) differentiation and fits predictions from models.