Abstract
While the analysis of temporal signal fluctuations and co-fluctuations has long been a fixture of blood oxygenation-level dependent (BOLD) functional magnetic resonance imaging (fMRI) research, the role and implications of spatial propagation within the 4D neurovascular BOLD signal has been almost entirely neglected. As part of a larger research program aimed at capturing and analyzing spatially propagative dynamics in BOLD fMRI, we report here a method that exposes large-scale functional attractors of spatial flows formulated as Markov processes defined at the voxel level. The brainwide stationary distributions of these voxel-level Markov processes represent patterns of signal accumulation toward which we find evidence that the brain exerts a probabilistic propagative undertow. These probabilistic propagative attractors are spatially structured and organized interpretably over functional regions. They also differ significantly between schizophrenia patients and controls.