Somatosensory Evoked BOLD ‐Signals With Ultra‐High Temporal Resolution

The blood oxygen level-dependent (BOLD) signal has been instrumental in characterizing brain activity. While the spatial resolution of fMRI continues to improve, relatively few methods have focused on enhancing and leveraging temporal resolution to investigate the spatiotemporal dynamics of the hemodynamic response. In this study, we applied a reordering method to achieve ultra-high temporal resolution (60 ms) in data acquired during a somatosensory stimulation paradigm. We then used a finite impulse response model (FIR) for each participant (N = 31) to preserve the temporal dynamics in the statistical analysis. At the group level, we employed an ANOVA combined with 4D nonparametric permutation testing to identify significant signal changes in time across the whole brain. Our results characterize the hemodynamic response in terms of both its temporal and spatial patterns and reveal distinct differences in response shapes within the somatosensory system. This method introduces a time-resolved approach to BOLD signal analysis, drawing inspiration from grand-average techniques commonly used in EEG research.