The present functional magnetic resonance imaging (fMRI) study investigated human brain regions subserving the discrimination of vibrotactile frequency. An event-related adaptation paradigm was used in which blood-oxygen-level-dependent (BOLD) responses are lower to same compared with different pairs of stimuli (BOLD adaptation). This adaptation effect serves as an indicator for feature-specific responding of neuronal subpopulations. Subjects had to discriminate two vibrotactile stimuli sequentially applied with a delay of 600 ms to their left middle fingertip. The stimulus frequency was in the flutter range of 18–26 Hz. In half of the trials, the two stimuli possessed identical frequency (same), whereas in the other half, a frequency difference of ±2 Hz was used (diff). As a result, BOLD adaptation was observed in the contralateral primary somatosensory cortex (S1), precentral gyrus, superior temporal gyrus (STG); ipsilateral insula as well as bilateral secondary somatosensory cortex and supplementary motor area. When statistically comparing the BOLD time courses between same and diff trials in these cortical areas, it was found that the vibrotactile BOLD adaptation is initiated in the contralateral S1 and STG simultaneously. These findings suggest that the cortical areas responsive to the frequency difference between two serially presented stimuli sequentially process the frequency of a vibrotactile stimulus and constitute a putative neuronal network underlying human vibrotactile frequency discrimination.