Although blood oxygenation level-dependent (BOLD) fMRI has been widely used to map brain responses to external stimuli and to delineate functional circuits at rest, the extent to which BOLD signals correlate spatially with underlying neuronal activity, the spatial relationships between stimulus-evoked BOLD activations and local correlations of BOLD signals in a resting state, and whether these spatial relationships vary across functionally distinct cortical areas are not known. To address these critical questions, we directly compared the spatial extents of stimulated activations and the local profiles of intervoxel resting state correlations for both high-resolution BOLD at 9.4 T and local field potentials (LFPs), using 98-channel microelectrode arrays, in functionally distinct primary somatosensory areas 3b and 1 in nonhuman primates. Anatomic images of LFP and BOLD were coregistered within 0.10 mm accuracy. We found that the point spread functions (PSFs) of BOLD and LFP responses were comparable in the stimulus condition, and both estimates of activations were slightly more spatially constrained than local correlations at rest. The magnitudes of stimulus responses in area 3b were stronger than those in area 1 and extended in a medial to lateral direction. In addition, the reproducibility and stability of stimulus-evoked activation locations within and across both modalities were robust. Our work suggests that the intrinsic resolution of BOLD is not a limiting feature in practice and approaches the intrinsic precision achievable by multielectrode electrophysiology.BOLD fMRI | local field potential | point spread function | resting state correlations | primary somatosensory cortex F unctional MRI (fMRI) is well established as a neuroimaging technique for detecting and delineating regions in the brain that change their levels of activity in response to specific experimental conditions (1-3). In addition, the discovery and analysis of synchronized fluctuations of low-frequency MRI signals between different brain regions at rest have provided a powerful approach to probe functional connectivity between regions and to delineate functional circuits (4-7). However, stimulus-evoked fMRI responses usually rely on detecting blood oxygenation leveldependent (BOLD) signal changes, which reflect hemodynamic processes, and thus are indirect indicators of neuronal activity. The measured extents of BOLD activations depend on the integrated contributions from the intrinsic spatial distributions of the neural activity involved, the effects of converting neural electrical activity to spatial distributions of metabolic and hemodynamic changes that then affect MRI signals, and the effects of image acquisitions and reconstruction with limited resolution, but the relative contributions of these to detected signals remain obscure. (14,15). No study, to our knowledge, has directly compared the spatial distributions of BOLD and LFP signals in both information processing (to external stimuli) and their correlation profiles in a resting...