We developed a data-driven method to spatiotemporally and spectrally characterize the dynamics of brain oscillations in resting-state magnetoencephalography (MEG) data. The method, called envelope spatial Fourier independent component analysis (eSFICA), maximizes the spatial and spectral sparseness of Fourier energies of a cortically constrained source current estimate. We compared this method using a simulated data set against 5 other variants of independent component analysis and found that eSFICA performed on par with its temporal variant, eTFICA, and better than other ICA variants, in characterizing dynamics at time scales of the order of minutes.

We then applied eSFICA to real MEG data obtained from 9 subjects during rest. The method identified several networks showing within- and cross-frequency inter-areal functional connectivity profiles which resemble previously reported resting-state networks, such as the bilateral sensorimotor network at ~ 20 Hz, the lateral and medial parieto-occipital sources at ~ 10 Hz, a subset of the default-mode network at ~ 8 and ~ 15 Hz, and lateralized temporal lobe sources at ~ 8 Hz.

Finally, we interpreted the estimated networks as spatiospectral filters and applied the filters to obtain the dynamics during a natural stimulus sequence presented to the same 9 subjects. We observed occipital alpha modulation to visual stimuli, bilateral rolandic mu modulation to tactile stimuli and video clips of hands, and the temporal lobe network modulation to speech stimuli, but no modulation of the sources in the default-mode network.

We conclude that (1) the proposed method robustly detects inter-areal cross-frequency networks at long time scales, (2) the functional relevance of the resting-state networks can be probed by applying the obtained spatiospectral filters to data from measurements with controlled external stimulation.