The characteristics of blood oxygenation level-dependent (BOLD) fMRI and magnetoencephalographic (MEG) responses to vibrotactile stimuli in humans were studied and compared. The stimuli, presented with interstimulus intervals (ISIs) ranging from 1 to 5 s, yielded highly reproducible MEG responses, with current dipoles in the primary somatosensory (SI) cortex in all subjects. BOLD fMRI responses to similar stimuli showed substantial intrasubject variation in the activation sites around the SI cortex. BOLD responses were detected in all subjects in the secondary somatosensory (SII) cortices as well, with comparable BOLD response amplitudes to those in the SI cortex. Current dipoles, used to model the MEG signals, were stronger at longer ISIs than shorter ISIs. The BOLD response amplitudes did not show a similar dependence on ISI, but the activated brain area was larger when longer ISIs or longer stimuli were applied. Our results support the view that combined use of brain mapping methods provides complementary information and should be considered in functional brain examinations.