Separating slow BOLD from non-BOLD baseline drifts using multi-echo fMRI

The functional magnetic resonance (fMRI) baseline is known to drift over the course of an experiment and is often attributed to hardware instability. These ultraslow fMRI fluctuations are inseparable from blood oxygenation level dependent (BOLD) changes in standard single echo fMRI and they are therefore typically removed before further analysis in both resting-state and task paradigms. However, some part of these fluctuations may be of neuronal origin, as neural activity can indeed fluctuate at the scale of several minutes or even longer, such as after the administration of drugs or during the ultradian rhythms. Here, we show that it is possible to separate the slow BOLD and non-BOLD drifts automatically using multi-echo fMRI and multi-echo independent components analysis (ME-ICA) denoising by demonstrating the detection of a visual signal evoked from a flickering checkerboard with slowly changing contrast.